SWBLI for High-Speed Intakes

The interaction of a shock wave with a boundary layer (SWBLI) occurs in many applications such as supersonic wind tunnel diffusers and supersonic (high-speed) intakes. Under specific operating conditions multiple shock wave boundary layer interactions (MSWBLIs) can form. MSWBLIs are often referred to as shock trains or pseudo-shocks. The formation of MSWBLIs affects the performance of high-speed intakes. An example is that the pressure recovery across MSWBLIs is lower than an equivalent single SWBLI. The project investigates the accuracy of the numerical computations of MSWBLIs, the sensitivity of MSWBLIs to modelling assumptions, and possible ways of MSWBLI control. At present Reynolds Averaged Navier Stokes (RANS), detached Eddy (DES), and scale-adaptive (SAS) simulations are used with the HMB3 CFD solver of Glasgow University to investigate the flow physics and the sensitivity to modelling assumptions of MSWBLIs in rectangular ducts.  Results are compared to the MSWBLI experiments by Carroll et al. [1] which offers extensive measurements. The use of linear turbulence models resulted in strong pressure oscillations (shown in the figure on the right of the page). The oscillations were attributed to the overprediction of the corner separations. The Explicit Algebraic Reynolds Stress Model (EARSM) resulted in better agreement with the experiment by reducing the size of the corner separations.  The addition of QCR extensions to the SST model also resulted in an improved agreement with the experiment. 

The SST model with the QCR extension and the EARSM model reduced the corner separations by accounting for the secondary flows near the duct corner. As the secondary flows are generated by the anisotropy of the normal Reynolds stresses only the models capable of accounting for the anisotropy will resolve the secondary flows. The wall shear stress visualisation, using friction lines just above the surface, below shows the difference in the predicted centreline and corner separations.

The linear models overpredict the corner separations and as a result the supersonic core-flow is closer to the duct wall. This causes pressure oscillations not observed in the experiment. The non-linear models reduce the pressure oscillations by reducing the corner separations and increasing the centreline separation. Once confidence was established in the simulations, simulations of shock trains in a more realistic geometry were considered. A fore-body intake geometry was created to study the effect of incidence angle, roll angle, cross-sectional area shape, and turbulence models on the shock train within the intake. The domain and boundary conditions are shown below.

The domain was discretised using multiblock structured grids. The resulting blocking and grid are showed below. 

Simulations were performed at different incidence and roll angles with the EARSM turbulence model. The SST model predicted just a single shock in some of the fore-body intake configurations which shows the importance of using non-linear models for complex shock wave boundary layer interactions like the presented one. The wall pressure and Mach number contours for one of the fore-body intake configuration are shown below.

The predicted shock train in the intake shares some similarity with the isolated shock trains in ducts. The main difference is the absence of symmetry due to the upstream interactions of the weak shocks with the boundary layer. The weak shocks are caused by the slant angle of the intake. For high total pressure recovery and low flow distortion, the shock train should be kept as far away as possible from the intake engine face. Ideally if a shock train forms in the intake it's length should be as short as possible in order to maximize the efficiency metrics. A informative presentation on the shock train phenomenon and numerical simulations is presented below. 


Boychev, K., Barakos, G.N., and Steijl, R. (2019). Flow physics and sensitivity to RANS modelling assumptions of a multiple shock wave/turbulent boundary layer interaction. Aerospace Science and Technology. 105640. 10.1016/j.ast.2019.105640. 

Boychev, K., Barakos, G.N., Steijl, R., and Shaw, S. (2020). High fidelity simulations of supersonic intakes. AIAA Scitech 2020 Forum and Exposition. 6.2020-2092. 10.2514/6.2020-2092.

K. Boychev, G.N. Barakos, R. Steijl, "Numerical Simulations of Multiple ShockWave Boundary Layer Interactions", Proceedings of the AIAA Scitech 2021 Forum and Exposition, 11-15 & 19-21 January 2021, Virtual Event.

K.Boychev, G.N. Barakos, R. Steijl, "Parametric Study of Multiple Shock Wave/Turbulent Boundary Layer Interactions with a Reynolds Stress Model", 8th European Congress on Computational Methods in Applied Science and Engineering (ECCOMAS 2020), 14th World Congress on Computational Mechanics (WCCM XIV), 11-15 January 2021, Virtual Congress.

N. Tylor et al., "The Prediction of Vortex Interactions on a Generic Missile Configuration Using CFD: Current Status of Activity in NATO AVT-316", Proceedings of the AVT-307 Research Symposium on Separated Flow: Prediction, Measurement and Assessment for Air and Sea Vehicles, 7-9 October 2019, Trondheim, Norway.


G. Barakos (Professor), George.Barakos@.gla.ac.uk
Rene Steijl (Senior lecturer), Rene.Steijl@gla.ac.uk
K. Boychev (PhD student), k.boychev.1@research.gla.ac.uk


[1] B.F. Carroll. “Numerical and Experimental Investigation of Multiple Shock Wave/Turbulent Boundary Layer Interactions in a Rectangular Duct”. PhD thesis. Urbana, IL: University of Illinois at Urbana-Champaign, 1988.