Perspectives on the Phenomenology and Modeling of Boundary Layer Transition

The characteristics of the coherent structures in a strongly decelerated large-velocity-defect boundary layer are analysed by direct numerical simulation. The simulated boundary layer starts as a zero-pressure-gradient boundary layer, decelerates under a strong adverse pressure gradient, and separates near the end of the domain, in the form of a very thin separation bubble. The Reynolds number at separation is R e _?? =3912 and the shape factor H=3.43. The three-dimensional spatial correlations of (u, u) and (u, v) are investigated and compared to those of a zero-pressure-gradient boundary layer and another strongly decelerated boundary layer. These velocity pairs lose coherence in the streamwise and spanwise directions as the velocity defect increases. In the outer region, the shape of the correlations suggest that large-scale u structures are less streamwise elongated and more inclined with respect to the wall in large-defect boundary layers. The three-dimensional properties of sweeps and ejections are characterized for the first time in both the zero-pressure-gradient and adverse-pressure-gradient boundary layers, following the method of Lozano-Durán et al. (J. Fluid Mech. 694, 100–130, [2012]). Although longer sweeps and ejections are found in the zero-pressure-gradient boundary layer, with ejections reaching streamwise lengths of 5 boundary layer thicknesses, the sweeps and ejections tend to be bigger in the adverse-pressure-gradient boundary layer. Moreover, small near-wall sweeps and ejections are much less numerous in the large-defect boundary layer. Large sweeps and ejections that reach the wall region (wall-attached) are also less numerous, less streamwise elongated and they occupy less space than in the zero-pressure-gradient boundary layer.     

Revisiting History Effects in Adverse-Pressure-Gradient Turbulent Boundary Layers

The goal of this study is to present a first step towards establishing criteria aimed at assessing whether a particular adverse-pressure-gradient (APG) turbulent boundary layer (TBL) can be considered well-behaved, i.e., whether it is independent of the inflow conditions and is exempt of numerical or experimental artifacts. To this end, we analyzed several high-quality datasets, including in-house numerical databases of APG TBLs developing over flat-plates and the suction side of a wing section, and five studies available in the literature. Due to the impact of the flow history on the particular state of the boundary layer, we developed three criteria of convergence to well-behaved conditions, to be used depending on the particular case under study. (i) In the first criterion, we develop empirical correlations defining the R e _?? -evolution of the skin-friction coefficient and the shape factor in APG TBLs with constant values of the Clauser pressure-gradient parameter β = 1 and 2 (note that β = δ ^?/τ _w dP _e /dx, where δ ^? is the displacement thickness, τ _w the wall-shear stress and dP _e /dx the streamwise pressure gradient). (ii) In the second one, we propose a predictive method to obtain the skin-friction curve corresponding to an APG TBL subjected to any streamwise evolution of β, based only on data from zero-pressure-gradient TBLs. (iii) The third method relies on the diagnostic-plot concept modified with the shape factor, which scales APG TBLs subjected to a wide range of pressure-gradient conditions. These three criteria allow to ensure the correct flow development of a particular TBL, and thus to separate history and pressure-gradient effects in the analysis.     

Streamwise Vortices and Velocity Streaks in a Locally Drag-Reduced Turbulent Boundary Layer

This work aims to understand the changes associated with the near-wall streaky structures in a turbulent boundary layer (TBL) where the local skin-friction drag is substantially reduced. The Reynolds number is R e _?? = 1000 based on the momentum thickness or R e _τ = 440 based on the friction velocity of the uncontrolled flow. The TBL is perturbed via a local surface oscillation produced by an array of spanwise-aligned piezo-ceramic (PZT) actuators and measurements are made in two orthogonal planes using particle image velocimetry (PIV). Data analyses are conducted using the vortex detection, streaky structure identification, spatial correlation and proper orthogonal decomposition (POD) techniques. It is found that the streaky structures are greatly modified in the near-wall region. Firstly, the near-wall streamwise vortices are increased in number and swirling strength but decreased in size, and are associated with greatly altered velocity correlations. Secondly, the velocity streaks grow in number and strength but contract in width and spacing, exhibiting a regular spatial arrangement. Other aspects of the streaky structures are also characterized; they include the spanwise gradient of the longitudinal fluctuating velocity and both streamwise and spanwise integral length scales. The POD analysis indicates that the turbulent kinetic energy of the streaky structures is reduced. When possible, our results are compared with those obtained by other control techniques such as a spanwise-wall oscillation, a spanwise oscillatory Lorentz force and a transverse traveling wave.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Evaluation of the Spectral Element Dynamic Model for Large-Eddy Simulation on Unstructured,Deformed Meshes

In this paper, the Spectral-Element Dynamic Model (SEDM), suited for Large-Eddy Simulation (LES) using Discontinuous Finite Element Methods (DFEM), is assessed using unstructured meshes. Five test cases of increasing complexity are considered, namely, the Taylor-Green vortex at Re =?5000, the turbulent channel flow at Re_τ =?587, the circular cylinder in cross-flow at Re_D =?3900, the square cylinder in cross-flow at Re_D =?22400 and the channel with periodic constrictions at Re_h =?10595. Various discretization parameters such as the grid spacing, polynomial degree and numerical flux are assessed and very accurate results are reported in all cases. This consistency in the results demonstrates the versatility of the SEDM approach and its ability to gage the actual resolution and quality of the mesh and, accordingly, to introduce an amount of sub-grid dissipation which is adapted to the spatial discretization considered.     

Effects of Lewis Number on Head on Quenching of Turbulent Premixed Flames: A Direct Numerical Simulation Analysis

The head on quenching of statistically planar turbulent premixed flames by an isothermal inert wall has been analysed using three-dimensional Direct Numerical Simulation (DNS) data for different values of global Lewis number Le(0.8, 1.0 and 1.2) and turbulent Reynolds number Re_t. The statistics of head on quenching have been analysed in terms of the wall Peclet number Pe (i.e. distance of the flame from the wall normalised by the Zel’dovich flame thickness) and the normalised wall heat flux Φ. It has been found that the maximum (minimum) value of Φ(Pe) for the turbulent Le=0.8 cases are greater (smaller) than the corresponding laminar value, whereas both Pe and Φ in turbulent cases remain comparable to the corresponding laminar values for Le=1.0 and 1.2. Detailed physical explanations are provided for the observed Le dependences of Pe and Φ. The existing closure of mean reaction rate \(\overline {\dot {\omega }}\) using the scalar dissipation rate (SDR) in the near wall region has been assessed based on a-priori analysis of DNS data and modifications to the existing closures of mean reaction rate and SDR have been suggested to account for the wall effects in such a manner that the modified closures perform well both near to and away from the wall.     

Effects of Fuel Lewis Number on Localised Forced Ignition of Globally Stoichiometric Stratified Mixtures: a Numerical Investigation

The influences of fuel Lewis number Le_F on localised forced ignition of globally stoichiometric stratified mixtures have been analysed using three-dimensional compressible Direct Numerical Simulations (DNS) for cases with Le_F ranging from 0.8 to 1.2. The globally stoichiometric stratified mixtures with different values of root-mean-square (rms) equivalence ratio fluctuation (i.e. ?^′= 0.2, 0.4 and 0.6) and the Taylor micro-scale l_? of equivalence ratio ? variation (i.e. l_?/l_f= 2.1, 5.5 and 8.3 with l_f being the Zel’dovich flame thickness of the stoichiometric laminar premixed flame) have been considered for different initial rms values of turbulent velocity u^′. A pseudo-spectral method is used to initialise the equivalence ratio variation following a presumed bi-modal distribution for prescribed values of ?^′ and l_?/l_f for global mean equivalence ratio 〈?〉=1.0. The localised ignition is accounted for by a source term in the energy transport equation that deposits energy for a stipulated time interval. It has been observed that the maximum values of temperature and the fuel reaction rate magnitude increase with decreasing Le_F during the period of external energy deposition. The initial values of Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f have been found to have significant effects on the extent of burning of the stratified mixtures following localised ignition. For a given value of u^′/S_b(?=1), the extent of burning decreases with increasing Le_F. An increase in u^′ leads to a monotonic reduction in the burned gas mass for all values of Le_F in all stratified mixture cases but an opposite trend is observed for the Le_F=0.8 homogeneous mixture. It has been found that an increase in ?^′ has adverse effects on the burned gas mass, whereas the effects of l_?/l_f on the extent of burning are non-monotonic and dependent on ?^′ and Le_F. Detailed physical explanations have been provided for the observed Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f dependences.     

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU _w ≡ 2U _w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u _τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U _w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u _τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.     

Experimental Control of Turbulent Boundary Layers with In-plane Travelling Waves

The experimental control of turbulent boundary layers using streamwise travelling waves of spanwise wall velocity, produced using a novel active surface, is outlined in this paper. The innovative surface comprises a pneumatically actuated compliant structure based on the kagome lattice geometry, supporting a pre-tensioned membrane skin. Careful design of the structure enables waves of variable length and speed to be produced in the flat surface in a robust and repeatable way, at frequencies and amplitudes known to have a favourable influence on the boundary layer. Two surfaces were developed, a preliminary module extending 152 mm in the streamwise direction, and a longer one with a fetch of 2.9 m so that the boundary layer can adjust to the new surface condition imposed by the forcing. With a shorter, 1.5 m portion of the surface actuated, generating an upstream-travelling wave, a drag reduction of 21.5% was recorded in the boundary layer with Re _τ =?1125. At the same flow conditions, a downstream-travelling produced a much smaller drag reduction of 2.6%, agreeing with the observed trends in current simulations. The drag reduction was determined with constant temperature hot-wire measurements of the mean velocity gradient in the viscous sublayer, while simultaneous laser Doppler vibrometer measurements of the surface recorded the wall motion. Despite the mechanics of the dynamic surface resulting in some out-of-plane motion (which is small in comparison to the in-plane streamwise movement), the positive drag reduction results are encouraging for future investigations at higher Reynolds numbers.     

Effect of Flow Structures on Turbulence Statistics of Taylor-Couette Flow in the Torque Transition State

Direct numerical simulations of Taylor-Couette flow from Re= 8000 to 25000 have been conducted to investigate changes of turbulence statistics in the transition of the Reynolds number dependency of the mean torque near Re= 10000. The velocity fluctuations are decomposed into the contributions of the Taylor vortex and remaining turbulent fluctuations. Significant Reynolds number dependencies of these components are observed in the radial profiles of the Reynolds stress and the transmission of the mean torque. The contributions of Taylor vortex and turbulent components in the net amount of mean torque are evaluated. The Taylor vortex component is overtaken by the turbulent counterpart around Re= 15000 when they are defined as the azimuthally averaged component and the remnants. The results show that the torque transition can be explained by the competition between the contributions of azimuthally averaged Taylor vortex and the remaining turbulent components.     

Assessment of Five SGS Models for Low-Rem MHD Turbulence

Assessment of three regularization-based and two eddy-viscosity-based subgrid-scale (SGS) turbulence models for large eddy simulations (LES) are carried out in the context of magnetohydrodynamic (MHD) decaying homogeneous turbulence (DHT) with a Taylor scale Reynolds number (Re_λ) of 120 and a MHD transition-to-turbulence Taylor-Green vortex (TGV) problems with a Reynolds number of 3000, through direct comparisons to direct numerical simulations (DNS). Simulations are conducted using the low-magnetic Reynolds number approximation (Re_m<<1). LES predictions using the regularization-based Leray- α,LANS- α, and Clark- α SGS models, along with the eddy viscosity-based non-dynamic Smagorinsky and the dynamic Smagorinsky models are compared to in-house DNS for DHT and previous results for TGV. With regard to the regularization models, this work represents their first application to MHD turbulence. Analyses of turbulent kinetic energy decay rates, energy spectra, and vorticity fields made between the varying magnetic field cases demonstrated that the regularization models performed poorly compared to the eddy-viscosity models for all MHD cases, but the comparisons improved with increase in magnitude of magnetic field, due to a decrease in the population of SGS eddies within the flow field.     

Improvement for expansion of parabolized stability equation method in boundary layer stability analysis

An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber α and its streamwise gradient dα/dx are unknown variables. This eigenvalue problem is solved for the eigenvalue dα/dx with an initial α, and the correction of α is performed with the conservation relation used in the PSE. The i EPSE is validated in several compressible and incompressible boundary layers. The computational results show that the prediction accuracy of the i EPSE is significantly higher than that of the ESPE, and it is in excellent agreement with the PSE which is regarded as the baseline for comparison. In addition, the unphysical multiple eigenmode problem in the EPSE is solved by using the i EPSE. As a local non-parallel stability analysis tool, the i EPSE has great potential application in the eNtransition prediction in general three-dimensional boundary layers.     

Consequences of the Translation Invariance on the Darcy Free Convection Flow Past a Vertical Surface

It is shown that the governing equation for the stream function of the Darcy free convection boundary layer flows past a vertical surface is invariant under arbitrary translations of the transverse coordinate y. The consequences of this basic symmetry property on the solutions corresponding to a prescribed surface temperature distribution T _w (x) are investigated. It is found that starting with a “primary solution” which describes the temperature boundary layer on an impermeable surface, infinitely many “translated solutions” can be generated which form a continuous group, the “translation group” of the given primary solution. The elements of this group describe free convection boundary layer flows from permeable counterparts of the original surface with a transformed temperature distribution \({\tilde {T}_w \left( x \right)}\), when simultaneously a suitable lateral suction/injection of the fluid is applied. It turns out in this way that several exact solutions discovered during the latter few decades are in fact not basically new solutions, but translated counterparts of some formerly reported primary solutions. A few specific examples are discussed in detail.     

Simulation of a Large-Eddy-Break-up Device (LEBU) in a Moderate Reynolds Number Turbulent Boundary Layer

A well-resolved large eddy simulation (LES) of a large-eddy break-up (LEBU) device in a spatially evolving turbulent boundary layer is performed with, Reynolds number, based on free-stream velocity and momentum-loss thickness, of R e _θ ≈ 4300. The implementation of the LEBU is via an immersed boundary method. The LEBU is positioned at a wall-normal distance of 0.8 δ (δ denoting the local boundary layer thickness at the location of the LEBU) from the wall. The LEBU acts to delay the growth of the turbulent boundary layer and produces global skin friction reduction beyond 180δ downstream of the LEBU, with a peak local skin friction reduction of approximately 12 %. However, no net drag reduction is found when accounting for the device drag of the LEBU in accordance with the towing tank experiments by Sahlin et al. (Phys. Fluids 31, 2814, 1988). Further investigation is performed on the interactions of high and low momentum bulges with the LEBU and the corresponding output is analysed, showing a ‘break-up’ of these large momentum bulges downstream of the LEBU. In addition, results from the spanwise energy spectra show consistent reduction in energy at spanwise length scales for \(\lambda _{z}^{+} > 1000\) independent of streamwise and wall-normal location when compared to the corresponding turbulent boundary layer without LEBU.     

Study of the Spectral Difference Numerical Dissipation for Turbulent Flows Using Unstructured Grids

In this paper, the numerical dissipation properties of the Spectral Difference (SD) method are studied in the context of vortex dominated flows and wall-bounded turbulence, using uniform and distorted grids. First, the validity of using the SD numerical dissipation as the only source of subgrid dissipation (the so-called Implicit-LES approach) is assessed on regular grids using various polynomial degrees (namely, p = 3, p = 4, p = 5) for the Taylor-Green vortex flow configuration at R e = 5 000. It is shown that the levels of numerical dissipation greatly depend on the order of accuracy chosen and, in turn, lead to an incorrect estimation of the viscous dissipation levels. The influence of grid distortion on the numerical dissipation is then assessed in the context of finite Reynolds number freely-decaying and wall-bounded turbulence. Tests involving different amplitudes of distortion show that highly skewed grids lead to the presence of small-scale, noisy structures, emphasizing the need of explicit subgrid modeling or regularization procedures when considering coarse, high-order SD computations on unstructured grids. Under-resolved, high-order computations of the turbulent channel flow at R e _τ = 1000 using highly-skewed grids are considered as well and present a qualitatively similar agreement to results obtained on a regular grid.     

A Simple Approach for Specifying Velocity Inflow Boundary Conditions in Simulations of Turbulent Opposed-Jet Flows

A new methodology is developed to specify inflow boundary conditions for the velocity field at the nozzle exit planes in turbulent counterflow simulations. The turbulent counterflow configuration consists of two coaxial opposed nozzles which emit highly-turbulent streams of varying species compositions depending on the mode considered. The specification of velocity inflow boundary conditions at the nozzle exits in the counterflow configuration is non-trivial because of the unique turbulence field generated by the turbulence generating plates (TGPs) upstream of the nozzle exits. In the method presented here, a single large-eddy simulation (LES) is performed in a large domain that spans the region between the TGPs of the nozzles, and the time series of the velocity fields at the nozzle exit planes are recorded. To provide inflow boundary conditions at the nozzle exit planes for simulations under other conditions (e.g., different stream compositions, bulk velocity, TGP location), transformations are performed on the recorded time series: the mean and r.m.s. (root-mean-square) quantities of velocity, as well as the longitudinal integral length scale on the centerline, at the nozzle exits in simulations are matched to those observed in experiments, thereby matching the turbulent Reynolds number R e _t . The method is assessed by implementing it in coupled large-eddy simulation/probability density function (LES/PDF) simulations on a small cylindrical domain between the nozzle exit planes for three different modes of the counterflow configuration: N ₂ vs. N ₂; N ₂ vs. hot combustion products; and C H ₄/N ₂ vs. O ₂. The inflow method is found to be successful as the first and second moments of velocity from the LES/PDF simulations agree well with the experimental data on the centerline for all three modes. This simple yet effective inflow strategy can be applied to eliminate the computational cost required to simulate the flow field upstream of the nozzle exits. It is also emphasized that, in addition to the predicted time series data, the availability of experimental data close to the nozzle exit planes plays a key role in the success of this method.     

Effect of Free-Stream Turbulence on the Reduction of Surface Friction in a Turbulent Boundary Layer Behind LEBU-Devices

The effect of increased free-stream turbulence on the reduction of the surface friction coefficient c _f in a turbulent boundary layer behind large-eddy break-up (LEBU) devices is investigated using a gravimetric method. The turbulence level was ε ≈ 1.9–4.9 % and the turbulence scale L _e ≈ 40–110 mm. The boundary layer Reynolds number Re** was varied from 2300 to 7500, with the boundary layer thickness being varied on the range δ = 33–44 mm. It is shown that an increase in the turbulence level ε has almost no impact on the relative reduction of friction behind LEBU-devices, whereas, under similar conditions of elevated free-stream turbulence, for another method, namely, the use of surface riblets, the friction reduction may be more strongly expressed.     

Kovasznay Mode Decomposition of Velocity-Temperature Correlation in Canonical Shock-Turbulence Interaction

The correlation coefficient R_uT between the streamwise velocity and temperature is investigated for the case of canonical shock-turbulence interaction, motivated by the fact that this correlation is an important component in compressible turbulence models. The variation of R_uT with the Mach number, the turbulent Mach number, and the Reynolds number is predicted using linear inviscid theory and compared to data from DNS. The contributions from the individual Kovasznay modes are quantified. At low Mach numbers, the peak post-shock R_uT is determined by the acoustic mode, which is correctly predicted by the linear theory. At high Mach numbers, it is determined primarily by the vorticity and entropy modes, which are strongly affected by nonlinear and viscous effects, and thus less well predicted by the linear theory.