Significance: Near a solid wall, the mean velocity of an incompressible fluid is a nearly universal function of the distance from the wall when appropriately transformed (nondimensionalized). For high-speed compressible flows, due to the wall-normal variations of density and viscosity, it has not been established how to transform a dimensional velocity profile into a universal profile. We propose a transformation that is more accurate and more broadly applicable than existing approaches. This improvement results from deploying different physical arguments for the two critical subdomains of a boundary layer. The transformation can be used to extend incompressible turbulence models to compressible turbulent flows, such as those encountered by vehicles for high-speed transportation and planetary reentry.
1 November 2017 in Journal of Computational Physics
While the computation of the boundary-layer thickness is straightforward for canonical equilibrium flows, there are no established definitions for general nonequilibrium flows. In this paper, a method is developed based on a local reconstruction of the inviscid velocity profile resulting from the application of the Bernoulli equation in the wall-normal direction. The boundary-layer thickness is then defined as the location where , which is consistent with its classical definition for the zero-pressure-gradient boundary layers. The proposed local-reconstruction method is parameter-free and can be deployed for both internal and external flows without resorting to an iterative procedure, numerical integration, or numerical differentiation. The superior performance of the local-reconstruction method over various existing methods is demonstrated by applying the methods to laminar and turbulent boundary layers and two flows over airfoils. Numerical experiments reveal that the local-reconstruction method is more accurate and more robust than existing methods, and it is applicable for flows over a wide range of Reynolds numbers.
The interaction between an incident shock wave and a Mach-6 undisturbed hypersonic laminar boundary layer over a cold wall is addressed using direct numerical simulations (DNS) and wall-modelled large-eddy simulations (WMLES) at different angles of incidence. At sufficiently high shock-incidence angles, the boundary layer transitions to turbulence via breakdown of near-wall streaks shortly downstream of the shock impingement, without the need of any inflow free-stream disturbances. The transition causes a localized significant increase in the Stanton number and skin-friction coefficient, with high incidence angles augmenting the peak thermomechanical loads in an approximately linear way. Statistical analyses of the boundary layer downstream of the interaction for each case are provided that quantify streamwise spatial variations of the Reynolds analogy factors and indicate a breakdown of the Morkovin's hypothesis near the wall, where velocity and temperature become correlated. A modified strong Reynolds analogy with a fixed turbulent Prandtl number is observed to perform best. Conventional transformations fail at collapsing the mean velocity profiles on the incompressible log law. The WMLES prompts transition and peak heating, delays separation and advances reattachment, thereby shortening the separation bubble. When the shock leads to transition, WMLES provides predictions of DNS peak thermomechanical loads within at a computational cost lower than DNS by two orders of magnitude. Downstream of the interaction, in the turbulent boundary layer, the WMLES agrees well with DNS results for the Reynolds analogy factor, the mean profiles of velocity and temperature, including the temperature peak, and the temperature/velocity correlation.