The Reynolds-averaged Navier–Stokes equations with one-equation turbulence models are used to simulate the flow field past a cone–flare geometry in the Mach number range from 5 to 8 with emphasis on the interaction region of the flare shock with the upstream boundary layer. A model based on the physics of shock unsteadiness is used to correct the standard Spalart–Allmaras turbulence model to improve the prediction of the extent of the separation bubble arising from the shock/turbulent boundary-layer interaction and its accompanying peak pressure and aerothermal loads on the surface. The computed results are validated against the experimental data. The limitations of the shock-unsteadiness model and the extent of the improvement in predicting the heat flux are discussed.