AbstractAn upwind finite element technique that uses cell‐centred quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive triangular grids. The approach is an extension to unstructured grids of the LAURA algorithm due to Gnoffo. A structured grid of quadrilaterals is laid out near a solid surface. For inviscid flows the method is stable at Courant numbers of over 100000. A first‐order basic scheme and a higher‐order flux‐corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view to simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are distinctly captured with a sequence of adaptively generate