Abstract

The need for better understanding of the low-frequency unsteadiness observed in shock wave/turbulent boundary layer interactions has been driving research in this area for several decades. We present here a large-eddy simulation investigation of the interaction between an impinging oblique shock and a Mach 2.3 turbulent boundary layer. Contrary to past large-eddy simulation investigations on shock/turbulent boundary layer interactions, we have used an inflow technique which does not introduce any energetically-significant low frequencies into the domain, hence avoiding possible interference with the shock/boundary layer interaction system. The large-eddy simulation has been run for much longer times than previous computational studies making a Fourier analysis of the low frequency possible. The broadband and energetic low-frequency component found in the interaction is in excellent agreement with the experimental findings. Furthermore, a linear stability analysis of the mean flow was performed and a stationary unstable global mode was found. The long-run large-eddy simulation data were analyzed and a phase change in the wall pressure fluctuations was found to coincide with the global-mode structure, leading to a possible driving mechanism for the observed low-frequency motions.Keywords shock boundary layer interaction · global mode · compressible turbulence · LES · low-frequency unsteadiness · separation bubble · digital filter · inflow turbulence PACS 47.40.Nm · 47.40.Ki · 47.27.ep

Direct numerical simulations (DNS) of laminar separation bubbles on a NACA-0012 airfoil at Re c = 5 × 10 4 and incidence 5• are presented. Initially volume forcing is introduced in order to promote transition to turbulence. After obtaining sufficient data from this forced case, the explicitly added disturbances are removed and the simulation run further. With no forcing the turbulence is observed to self-sustain, with increased turbulence intensity in the reattachment region. A comparison of the forced and unforced cases shows that the forcing improves the aerodynamic performance whilst requiring little energy input. Classical linear stability analysis is performed upon the time-averaged flow field; however no absolute instability is observed that could explain the presence of self-sustaining turbulence. Finally, a series of simplified DNS are presented that illustrate a three-dimensional absolute instability of the twodimensional vortex shedding that occurs naturally. Three-dimensional perturbations are amplified in the braid region of developing vortices, and subsequently convected upstream by local regions of reverse flow, within which the upstream velocity magnitude greatly exceeds that of the time-average. The perturbations are convected into the braid region of the next developing vortex, where they are amplified further, hence the cycle repeats with increasing amplitude. The fact that this transition process is independent of upstream disturbances has implications for modelling separation bubbles.

Direct numerical simulation databases have been used to study the effect of compressibility on mixing layers. The simulations cover convective Mach numbers from 0.2 to 1.2 and all contain a fully resolved turbulent energy cascade to small spatial scales. Statistical information is extracted from the databases to determine reasons for the reduced growth rate that is observed as the convective Mach number is increased. It is found that the dilatational contribution to dissipation is negligible even when eddy shocklets are observed in the flow. Also pressure-dilatation is not found to be significant. Using an accurate relation between the momentum thickness growth rate and the production of turbulence kinetic energy together with integrated equations for the Reynolds stress tensor it is shown that reduced pressure fluctuations are responsible for the changes in growth rate via the pressure-strain term. A deterministic model for the required pressure fluctuations is given based on the structure of variable-density vortices and the assumption that the limiting eddies are sonic. Simple anisotropy considerations are used to close the averaged equations. Good agreement with turbulence statistics obtained from the simulations is found.

Direct numerical simulation of the incompressible Navier-Stokes equations is used to study flows where laminar boundary-layer separation is followed by turbulent reattachment forming a closed region known as a laminar separation bubble. In the simulations a laminar boundary layer is forced to separate by the action of a suction profile applied as the upper boundary condition. The separated shear layer undergoes transition via oblique modes and Λ-vortex-induced breakdown and reattaches as turbulent flow, slowly recovering to an equilibrium turbulent boundary layer. Compared with classical experiments the computed bubbles may be classified as 'short', as the external potential flow is only affected in the immediate vicinity of the bubble. Near reattachment budgets of turbulence kinetic energy are dominated by turbulence events away from the wall. Characteristics of near-wall turbulence only develop several bubble lengths downstream of reattachment. Comparisons are made with two-dimensional simulations which fail to capture many of the detailed features of the full three-dimensional simulations. Stability characteristics of mean flow profiles are computed in the separated flow region for a family of velocity profiles generated using simulation data. Absolute instability is shown to require reverse flows of the order of 15-20%. The three-dimensional bubbles with turbulent reattachment have maximum reverse flows of less than 8% and it is concluded that for these bubbles the basic instability is convective in nature.