scholarly journals Self-similar compressible turbulent boundary layers with pressure gradients. Part 1. Direct numerical simulation and assessment of Morkovin’s hypothesis

2019 ◽  
Vol 880 ◽  
pp. 239-283 ◽  
Author(s):  
Christoph Wenzel ◽  
Tobias Gibis ◽  
Markus Kloker ◽  
Ulrich Rist

A direct numerical simulation study of self-similar compressible flat-plate turbulent boundary layers (TBLs) with pressure gradients (PGs) has been performed for inflow Mach numbers of 0.5 and 2.0. All cases are computed with smooth PGs for both favourable and adverse PG distributions (FPG, APG) and thus are akin to experiments using a reflected-wave set-up. The equilibrium character allows for a systematic comparison between sub- and supersonic cases, enabling the isolation of pure PG effects from Mach-number effects and thus an investigation of the validity of common compressibility transformations for compressible PG TBLs. It turned out that the kinematic Rotta–Clauser parameter $\unicode[STIX]{x1D6FD}_{K}$ calculated using the incompressible form of the boundary-layer displacement thickness as length scale is the appropriate similarity parameter to compare both sub- and supersonic cases. Whereas the subsonic APG cases show trends known from incompressible flow, the interpretation of the supersonic PG cases is intricate. Both sub- and supersonic regions exist in the boundary layer, which counteract in their spatial evolution. The boundary-layer thickness $\unicode[STIX]{x1D6FF}_{99}$ and the skin-friction coefficient $c_{f}$, for instance, are therefore in a comparable range for all compressible APG cases. The evaluation of local non-dimensionalized total and turbulent shear stresses shows an almost identical behaviour for both sub- and supersonic cases characterized by similar $\unicode[STIX]{x1D6FD}_{K}$, which indicates the (approximate) validity of Morkovin’s scaling/hypothesis also for compressible PG TBLs. Likewise, the local non-dimensionalized distributions of the mean-flow pressure and the pressure fluctuations are virtually invariant to the local Mach number for same $\unicode[STIX]{x1D6FD}_{K}$-cases. In the inner layer, the van Driest transformation collapses compressible mean-flow data of the streamwise velocity component well into their nearly incompressible counterparts with the same $\unicode[STIX]{x1D6FD}_{K}$. However, noticeable differences can be observed in the wake region of the velocity profiles, depending on the strength of the PG. For both sub- and supersonic cases the recovery factor was found to be significantly decreased by APGs and increased by FPGs, but also to remain virtually constant in regions of approximated equilibrium.

2019 ◽  
Vol 877 ◽  
pp. 167-195 ◽  
Author(s):  
Feng-Yuan Zuo ◽  
Antonio Memmolo ◽  
Guo-ping Huang ◽  
Sergio Pirozzoli

Direct numerical simulation of the Navier–Stokes equations is carried out to investigate the interaction of a conical shock wave with a turbulent boundary layer developing over a flat plate at free-stream Mach number $M_{\infty }=2.05$ and Reynolds number $Re_{\unicode[STIX]{x1D703}}\approx 630$, based on the upstream boundary layer momentum thickness. The shock is generated by a circular cone with half opening angle $\unicode[STIX]{x1D703}_{c}=25^{\circ }$. As found in experiments, the wall pressure exhibits a distinctive N-wave signature, with a sharp peak right past the precursor shock generated at the cone apex, followed by an extended zone with favourable pressure gradient, and terminated by the trailing shock associated with recompression in the wake of the cone. The boundary layer behaviour is strongly affected by the imposed pressure gradient. Streaks are suppressed in adverse pressure gradient (APG) zones, but re-form rapidly in downstream favourable pressure gradient (FPG) zones. Three-dimensional mean flow separation is only observed in the first APG region associated with the formation of a horseshoe vortex, whereas the second APG region features an incipient detachment state, with scattered spots of instantaneous reversed flow. As found in canonical geometrically two-dimensional wedge-generated shock–boundary layer interactions, different amplification of the turbulent stress components is observed through the interacting shock system, with approach to an isotropic state in APG regions, and to a two-component anisotropic state in FPG. The general adequacy of the Boussinesq hypothesis is found to predict the spatial organization of the turbulent shear stresses, although different eddy viscosities should be used for each component, as in tensor eddy-viscosity models, or in full Reynolds stress closures.


2001 ◽  
Vol 429 ◽  
pp. 187-216 ◽  
Author(s):  
THIERRY MAEDER ◽  
NIKOLAUS A. ADAMS ◽  
LEONHARD KLEISER

The present paper addresses the direct numerical simulation of turbulent zero-pressure-gradient boundary layers on a flat plate at Mach numbers 3, 4.5 and 6 with momentum-thickness Reynolds numbers of about 3000. Simulations are performed with an extended temporal direct numerical simulation (ETDNS) method. Assuming that the slow streamwise variation of the mean boundary layer is governed by parabolized Navier–Stokes equations, the equations solved locally in time with a temporal DNS are modified by a distributed forcing term so that the parabolized Navier–Stokes equations are recovered for the spatial average. The correct mean flow is obtained without a priori knowledge, the streamwise mean-flow evolution being approximated from its upstream history. ETDNS reduces the computational effort by up to two orders of magnitude compared to a fully spatial simulation.We present results for a constant wall temperature Tw chosen to be equal to its laminar adiabatic value, which is about 2.5 T∞, 4.4 T∞ and 7 T∞, respectively, where T∞ is the free-stream temperature for the three Mach numbers considered. The simulations are initialized with transition-simulation data or with re-scaled turbulent data at different parameters. We find that the ETDNS results closely match experimental mean-flow data. The van Driest transformed velocity profiles follow the incompressible law of the wall with small logarithmic regions.Of particular interest is the significance of compressibility effects in a Mach number range around the limit of M∞ ≃ 5, up to which Morkovin's hypothesis is believed to be valid. The results show that pressure dilatation and dilatational dissipation correlations are small throughout the considered Mach number range. On the other hand, correlations derived from Morkovin's hypothesis are not necessarily valid, as is shown for the strong Reynolds analogy.


1992 ◽  
Vol 238 ◽  
pp. 699-722 ◽  
Author(s):  
P. A. Durbin ◽  
S. E. Belcher

An asymptotic analysis is developed for turbulent boundary layers in strong adverse pressure gradients. It is found that the boundary layer divides into three distinguishable regions: these are the wall layer, the wake layer and a transition layer. This structure has two key differences from the zero-pressure-gradient boundary layer: the wall layer is not exponentially thinner than the wake; and the wake has a large velocity deficit, and cannot be linearized. The mean velocity profile has a y½ behaviour in the overlap layer between the wall and transition regions.The analysis is done in the context of eddy viscosity closure modelling. It is found that k-ε-type models are suitable to the wall region, and have a power-law solution in the y½ layer. The outer-region scaling precludes the usual ε-equation. The Clauser, constant-viscosity model is used in that region. An asymptotic expansion of the mean flow and matching between the three regions is carried out in order to determine the relation between skin friction and pressure gradient. Numerical calculations are done for self-similar flow. It is found that the surface shear stress is a double-valued function of the pressure gradient in a small range of pressure gradients.


2007 ◽  
Vol 570 ◽  
pp. 347-364 ◽  
Author(s):  
M. PINO MARTIN

A systematic procedure for initializing supersonic and hypersonic turbulent boundary layers at controlled Mach number and Reynolds number conditions is described. The initialization is done by locally transforming a true direct numerical simulation flow field, and results in a nearly realistic initial magnitude of turbulent fluctuations, turbulence structure and energy distribution. The time scales necessary to forget the initial condition are studied. The experimental conditions of previous studies are simulated. The magnitude of velocity and temperature fluctuations, as well as the turbulent shear stresses given by the direct numerical simulations are in agreement with the experimental data.


1983 ◽  
Vol 27 (03) ◽  
pp. 147-157 ◽  
Author(s):  
A. J. Smits ◽  
N. Matheson ◽  
P. N. Joubert

This paper reports the results of an extensive experimental investigation into the mean flow properties of turbulent boundary layers with momentum-thickness Reynolds numbers less than 3000. Zero pressure gradient and favorable pressure gradients were studied. The velocity profiles displayed a logarithmic region even at very low Reynolds numbers (as low as Rθ = 261). The results were independent of the leading-edge shape, and the pin-type turbulent stimulators performed well. It was found that the shape and Clauser parameters were a little higher than the correlation proposed by Coles [10], and the skin friction coefficient was a little lower. The skin friction coefficient behavior could be fitted well by a simple power-law relationship in both zero and favorable pressure gradients.


1966 ◽  
Vol 33 (2) ◽  
pp. 429-437 ◽  
Author(s):  
J. C. Rotta

A review is given of the recent development in turbulent boundary layers. At first, for the case of incompressible flow, the variation of the shape of velocity profile with the pressure gradient is discussed; also the temperature distribution and heat transfer in incompressible boundary layers are treated. Finally, problems of the turbulent boundary layer in compressible flow are considered.


2007 ◽  
Vol 594 ◽  
pp. 59-69 ◽  
Author(s):  
MATTHEW J. RINGUETTE ◽  
MINWEI WU ◽  
M. PINO MARTÍN

We demonstrate that data from direct numerical simulation of turbulent boundary layers at Mach 3 exhibit the same large-scale coherent structures that are found in supersonic and subsonic experiments, namely elongated, low-speed features in the logarithmic region and hairpin vortex packets. Contour plots of the streamwise mass flux show very long low-momentum structures in the logarithmic layer. These low-momentum features carry about one-third of the turbulent kinetic energy. Using Taylor's hypothesis, we find that these structures prevail and meander for very long streamwise distances. Structure lengths on the order of 100 boundary layer thicknesses are observed. Length scales obtained from correlations of the streamwise mass flux severely underpredict the extent of these structures, most likely because of their significant meandering in the spanwise direction. A hairpin-packet-finding algorithm is employed to determine the average packet properties, and we find that the Mach 3 packets are similar to those observed at subsonic conditions. A connection between the wall shear stress and hairpin packets is observed. Visualization of the instantaneous turbulence structure shows that groups of hairpin packets are frequently located above the long low-momentum structures. This finding is consistent with the very large-scale motion model of Kim & Adrian (1999).


Sign in / Sign up

Export Citation Format

Share Document