Reynolds Number Dependence Research Articles

Practical computational fluid dynamic predictions of a cyclist in a time trial position

On a flat road, at race speeds, aerodynamic drag is the main resistive force a cyclist must overcome. Computational fluid dynamics (CFD) can be a useful tool to predict and understand the complex flow and, therefore, drive developments to reduce drag. However, cycling aerodynamics is complex. The effects of Reynolds number, surface roughness, boundary layer transition, flow separation, and turbulent wakes are challenging to accurately predict. High fidelity time-resolved computations, such as Large eddy simulations (LES), require high-performance computing and lengthy simulation times. This paper examines whether lower fidelity CFD, such as Reynolds averaged approaches, can predict the drag of a cyclist with sufficient accuracy and within practical timescales on a desktop PC. Wind tunnel tests of a rider model (without bicycle) were conducted at Reynolds numbers equivalent to speeds of ~ 20–70 km/h. Measured drag showed a notable Reynolds number dependency with the drag coefficient reducing almost linearly by ~ 20% from 0.88 to 0.71. The computational accurately replicated this relationship but only when employing a boundary layer transition model. The steady computations underpredicted the magnitude of the measured drag coefficient by ~ 3% but the unsteady computations were within ~ 2%. Examination of the predicted flow field revealed variations in boundary layer transition, separation, and wake formation from each body part which combine in a complex wake system. Overall, the data confirm validity and suitable accuracy of the CFD, and therefore this provides a practical time and cost-effective tool for further examination of drag reduction within cycling.

Growth characteristics of the mean shear layer in pipe bends with and without a guide vane

The growth characteristics of two identical pipe bends with and without a guide vane are investigated by means of large eddy simulation. The two pipe bends, with a radius of curvature slightly above the separation threshold, are subjected to two fully developed upstream flow conditions, with a corresponding Reynolds number of 11 700 and 24 000. A precursor computation method is employed to provide the fully developed turbulence inflow conditions for all cases. The growth of the mean shear layer in this work is characterized by the local momentum thickness, which measures the extent of momentum deficit confined under the mean shear layer. For both pipe bends, the initial growth of momentum thickness is observed in the first quarter of the bend. The onset location is almost independent of the Reynolds numbers. However, a clear Reynolds number dependence is observed in the onset magnitude, which strongly defines the growth rate thereafter. By examining the mean momentum balance in the bend section, the results show that rather than the adverse pressure gradient, the overall growth characteristics of the mean shear layer, which include the onset location and the growth rate, are better described by the balance between the centrifugal force and the radial pressure gradient. This balance manifests itself as a change in the swirling intensity of the secondary flow. The presence of guide vane significantly suppresses the swirling intensity in the bend section, leading to a noticeable reduction in the overall momentum thickness growth and the production of turbulence in the flow downstream of the bend. Further inspection also indicates that the initial mechanism leading to the suppression of separation at the inner bend is linked to the increasing dominance of the small vortices at the near-wall vicinity relative to the local adverse pressure gradient. Certain aspects pertaining to the turbulent statistics are also discussed.

Wall Shear Stress Characteristics Of A Turbulent Boundary Layer Obtained With Multi-Aperture Defocusing Micro-PTV And High-Speed Stereo Profile-PIV

A multi-aperture 3d 3c micro PTV system is presented which relies on single window optical access for both high-speed tracer illumination and image recording. Similar to the “defocusing” concept described by C. Willert & Gharib (1992), the wall distance of individual particles is obtained from the size of projected particle image triplets formed by a triplet of apertures on the entrance pupil of the microscope lens. The present article extends upon previously published material Klinner & Willert (2022) by applying multi-aperture micro particle tracking velocimetry (MA-micro PTV) on a canonical turbulent boundary layer (TBL) to track the near-wall motion of tracer particles with the aim of estimating the unsteady wall-shear stress from particle tracks. The technique is validated with measurements of a TBL inside the closed test section of the 1 m wind tunnel of DLR in Göttingen (1mWK) at free- stream velocities of 5.2 to 20 m/s with corresponding shear Reynolds numbers between 560 and 1630. In the viscous sublayer, the probability density distributions of stream- and span-wise wall shear stress (WSS) components could be reliably captured down to probability densities of 10E−3. As far as we know, this has not been achieved in the past, particularly not for the span-wise component. For the stream-wise component the measured Reynolds-number dependency of the root mean squared fluctuations agrees to correlations by Örlü & Schlatter (2011). Furthermore, the skewness of the stream-wise WSS agrees well to values from DNS. On the other hand, the span-wise WSS fluctuations consistently underestimate the predicted DNS values, for which it is assumed that the deviation can be explained by the rapid decrease of the span-wise velocity fluctuations with increasing wall distance. Wall-normal profiles of 3c velocity statistics were obtained by bin averaging of particle velocities. Up into the buffer layer, these profiles agree well with profiles from stereo PIV as well as from DNS and LES, indicating an accurate determination of both near-wall velocities and inner scaling.

Attempt to compensate for effects of energy conservation error on static pressure fluctuation using sub-grid scale components with the Reynolds number dependence of isotropic/anisotropic steady turbulence

Attempt to compensate for effects of energy conservation error on static pressure fluctuation using sub-grid scale components with the Reynolds number dependence of isotropic/anisotropic steady turbulence

Numerical analysis of static and dynamic wind turbine airfoil characteristics in transonic flow

This study performed an aerodynamic characterization of the FFA-W3-211 wind turbine tip airfoil in transonic flow using Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations, for both steady and dynamic operational conditions. First, the boundary between subsonic and supersonic flow in static conditions was identified, depending on the angle of attack, the approach flow Mach number, and the Reynolds number. The analysis points out that higher Reynolds numbers promote the occurrence of local supersonic flow. Thereafter, to investigate the dynamic behavior in the transonic flow regime, a sinusoidal pitching motion with representative values was imposed. A hysteresis, similar to but distinct from dynamic stall, was observed for entering and leaving the supersonic and subsonic regions. Elevated reduced frequencies widened the hysteresis loop, resulting in increased normal forces on the airfoil. The study indicated that an increase in reduced frequency leads to an earlier onset of transonic flow. In conclusion, the risk of transonic flow occurring during normal operation of the next generation wind turbines predicted in earlier studies could be corroborated. Moreover, dynamic effects and Reynolds number dependencies can be significant.

Drag on circular cylinders with porous outer layers in turbulent cross-flow

The effects of free-stream turbulence intensity and porosity on the drag on cylinders with porous outer layers in cross-flow was investigated experimentally. This work is motivated by the need to better model spotting—a forest fire propagation mechanism in which burning branches and other debris (termed firebrands) are transported away from the main fire by the prevailing wind and ignite new fires. Multiple levels of background turbulence were studied by using no grid, passive grids, and an active grid to generate turbulence intensities of 0.4%, 1.7%, 2.7% and 12.4%. The porous char-layer on the outer surface of firebrands was mimicked by wrapping wire meshes (of 10, 20 and 40 pores per cylinder diameter or PPD) around the cylinders, each to three different layer-thickness fractions ( 1/16 , 1/8 and 1/4 ) of the cylinder’s outer diameter. The drag on one smooth cylinder and nine cylinders with porous outer layers was measured for Reynolds numbers in the range 7000–17 000, at the aforementioned four turbulence intensities. The results showed that (i) the free-stream turbulence intensity and PPD of the wire meshes affect the Reynolds number dependence of the drag coefficient; (ii) the drag coefficient increases with free-stream turbulence intensity when it is relatively low (0.4%–2.7%), then decreases at high intensities (12.4%), and this decrease is more pronounced as the PPD increases; (iii) the drag coefficient increases with the thickness of the porous layer, and asymptotes to an effectively constant value after a critical thickness of about 1/8 of the cylinder’s diameter; and iv) the drag coefficient exhibits a non-monotonic dependence on the PPD of the wire mesh. These results demonstrate the importance of accounting for free-stream turbulence intensity and the parameters that characterize the porous outer layers of cylinders when modeling the drag on firebrands, or in other applications in which there is cross flow over cylinders with porous outer surfaces.

Reynolds number dependence of turbulent flows over a highly permeable wall

Direct numerical simulations of turbulent flows over highly permeable porous walls were performed at various Reynolds numbers to examine the effects of the Reynolds number on permeable wall turbulence. The porous medium consisted of Kelvin cell arrays with porosity $0.95$ , and the permeability Reynolds number $Re_K$ ranged from approximately 7 to 50. Simulations with thin and thick porous walls were performed to investigate the effects of spanwise roller vortices associated with the Kelvin–Helmholtz instability. The results show that the effect of the Kelvin–Helmholtz instability becomes more significant with increasing the permeability Reynolds number, and spanwise rollers, for which length scale is an order of channel height, dominate turbulence when $Re_K \gtrsim 30$ . Spanwise rollers reinforce the negative correlation between the wall-normal and streamwise velocity fluctuations close to the porous/fluid interface, and intensify the turbulent velocity fluctuations away from the porous walls, leading to increased frictional resistance. An investigation of the Reynolds number dependence of the modified logarithmic law indicates that the zero-plane displacement and equivalent roughness height are proportional to the square root of permeability, whereas the von Kármán constant increases with the permeability Reynolds number because of the increased mixing length resulting from the relatively large-scale velocity fluctuations induced by spanwise rollers. We developed a model for the modified log law for permeable wall turbulence based on permeability, and confirmed that the skin friction coefficient obtained from the model reasonably predicts the skin friction coefficient for several types of high-porosity porous media. Hence, permeability is a key parameter that characterizes the logarithmic mean velocity profiles over a variety of porous media with high porosity.

Reynolds number asymptotics of wall-turbulence fluctuations

In continuation of our earlier work (Chen & Sreenivasan, J. Fluid Mech., vol. 908, 2021, R3; Chen & Sreenivasan, J. Fluid Mech., vol. 933, 2022a, A20 – together referred to as CS hereafter), we present a self-consistent Reynolds number asymptotics for wall-normal profiles of variances of streamwise and spanwise velocity fluctuations as well as root-mean-square pressure, across the entire flow region of channel and pipe flows and flat-plate boundary layers. It is first shown that, when normalized by peak values, the Reynolds number dependence and wall-normal variation of all three profiles can be decoupled, in excellent agreement with available data, sharing the common inner expansion of the type $\phi (y^+)=f_0(y^+)+f_1(y^+)/Re^{1/4}_\tau$ , where $\phi$ is one of the quantities just mentioned, the functions $f_0$ and $f_1$ depend only on $y^+$ , and $Re_\tau$ is the friction Reynolds number. Here, the superscript $+$ indicates normalization by wall variables. We show that this result is completely consistent with CS. Secondly, by matching the above inner expansion and the outer flow similarity form, a bounded variation $\phi (y^\ast )=\alpha _\phi -\beta _{\phi }y^{{\ast {1}/{4}}}$ is derived for the outer region where, for each $\phi$ , the constants $\alpha _\phi$ and $\beta _{\phi }$ are independent of $Re_\tau$ and $y^\ast$ $\equiv y^+/Re_\tau$ – also in excellent agreement with simulations and experimental data. One of the predictions of the analysis is that, for asymptotically high Reynolds numbers, a finite plateau $\phi \approx \alpha _\phi$ appears in the outer region. This result sheds light on the intriguing issue of the outer shoulder of the variance of the streamwise velocity fluctuation, which should be bounded by the asymptotic plateau of approximately 10.

Reynolds number dependence of turbulent kinetic energy and energy balance of 3-component turbulence intensity in a pipe flow

Measurement data sets are presented for the turbulence intensity profile of three velocity components ( $u$ , $v$ and $w$ ) and turbulent kinetic energy (TKE, $k$ ) over a wide range of Reynolds numbers from $Re_\tau =990$ to 20 750 in a pipe flow. The turbulence intensity profiles of the $u$ - and $w$ -component show logarithmic behaviour, and that of the $v$ -component shows a constant region at high Reynolds numbers, $Re_\tau >10\,000$ . Furthermore, a logarithmic region is also observed in the TKE profile at $y/R=0.055\unicode{x2013}0.25$ . The Reynolds number dependences of peak values of $u$ -, $w$ -component and TKE fit to both a logarithmic law (Marusic et al., Phys. Rev. Fluids, vol. 2, 2017, 100502) and an asymptotic law (Chen and Sreenivasan, J. Fluid Mech., vol. 908, 2020, R3), within the uncertainty of measurement. The Reynolds number dependence of the bulk TKE $k^+_{bulk}$ , which is the total amount of TKE in the cross-sectional area of the pipe also fits to both laws. When the asymptotic law is applied to the $k^+_{bulk}$ , it asymptotically increases to the finite value $k^+_{bulk}=11$ as the Reynolds number increases. The contribution ratio $\langle u'^2\rangle /k$ reaches a plateau, and the value tends to be constant within $100< y^+<1000$ at $Re_\tau >10\ 000$ . Therefore, the local balance of each velocity component also indicates asymptotic behaviour. The contribution ratios are balanced in this region at high Reynolds numbers as $\langle u'^2\rangle /k\simeq 1.25$ , $\langle w'^2\rangle /k\simeq 0.5$ and $\langle v'^2\rangle /k \simeq 0.25$ .

Intermittency across Reynolds numbers – the influence of large-scale shear layers on the scaling of the enstrophy and dissipation in homogenous isotropic turbulence

Direct numerical simulations up to Reλ = 1445 show that the scaling exponents for the enstrophy and the dissipation rate extrema are different and depend on the Reynolds number. A similar Reynolds number dependence of the scaling exponents is observed for the moments of the dissipation rate, but not for the moments of the enstrophy. Significant changes in the exponents occur at approximately Reλ ≈ 250, where Reλ is the Taylor based Reynolds number, which coincides with structural changes in the flow, in particular the development of large-scale shear layers. A model for the probability density functions (PDFs) of the enstrophy and dissipate rate is presented, which is an extension of our existing model (Proc. R. Soc. A, vol. 476, 2020, p. 20200591) and is based on the mentioned development of large-scale layer regions within the flow. This model is able to capture the observed Reynolds number dependencies of the scaling exponents, in contrast to the existing theories which yield constant exponents. Moreover, the model reconciles the scaling at finite Reynolds number with the theoretical limit, where the enstrophy and dissipation rate scale identically at infinite Reynolds number. It suggests that the large-scale shear layers are vital for understanding the scaling of the extrema. Furthermore, to reach the theoretical limit, the scaling exponents must remain Reynolds number dependent beyond the present Reλ range.

Reynolds Number Dependency of Wall-Bounded Turbulence Over a Surface Partially Covered by Barnacle Clusters

The settlement of barnacles on a ship hull is a common form of marine biofouling. In this study, the Reynolds number dependency of turbulent flow over a surface partially covered by barnacle clusters is investigated using direct numerical simulations of turbulent channel flow at friction Reynolds numbers ranging from 180 to 720. Mean flow, Reynolds and dispersive stress statistics are evaluated and compared to the corresponding results for a generic irregular rough surface with a Gaussian height distribution. For the barnacle surface, distinctive features emerge in the velocity statistics due to the interplay between the barnacle clusters and the large, connected smooth-wall sections surrounding them. This aspect is further investigated by applying a rough-smooth decomposition to the local time-averaged flow statistics for the barnacle surface. Using this decomposition, the partial recovery of smooth-wall behaviour over the smooth sections of the barnacle surface can be observed in the Reynolds stress statistics with the streamwise Reynolds stresses exhibiting a similar behaviour as previously found for boundary layers over surfaces with a rough to smooth transition.

Reynolds number dependence of azimuthal and streamwise pipe flow structures

This paper revisits the study by Baileyet al.(J. Fluid Mech., vol. 615, 2008, pp. 121–138), adopting a higher-fidelity calibration approach to reveal subtle flow variations with Reynolds numbers that were not discernible previously. The paper aims therefore to provide insights into the characteristics of azimuthal and streamwise pipe flow structures adopting two-point joint statistics and spectral analysis for shear Reynolds numbers in the range$2\times {10^3}\le {Re_\tau }\le {16\times {10^3}}$, where${Re_\tau }$is based on the wall friction velocity$u_{\tau }$, the pipe radius$R$, and the fluid kinematic viscosity$\nu$. The streamwise velocity fluctuations were measured at four wall-normal locations,$0.1\le {x_{{2}}/R}\le {0.7}$, covering the logarithmic and core regions of fully developed turbulent pipe flow based on 35–41 azimuthal probe separations using, simultaneously, two single hot-wire probes. A uniquein situcalibration approach for both probes was adopted where a potential flow was insured, resulting in consistent and precise pipe flow data. The azimuthal velocity correlation, the cross-power spectral density and the coherence function of the streamwise velocity fluctuations are discussed, revealing a clear dependence of the azimuthal scales of the large and very large flow motions on the wall-normal location, the azimuthal separation, the streamwise wavenumber and the Reynolds number. Along the logarithmic region, a linear growth of the azimuthal scales of the large- and very-large-scale structures was observed; however, they do scale nonlinearly and reach their maximum sizes in the core region, i.e. near the centreline of the pipe. Additionally, the streamwise very-large- and large-scale motions were evaluated using the premultiplied energy spectra, showing wavelengths${\approx }{18R}$and${\approx }{3R}$for${Re_{\tau }}\approx {16\times {10^3}}$at half of the pipe radius, respectively.

Structured input–output analysis of oblique laminar–turbulent patterns in plane Couette–Poiseuille flow

In this work, we employ structured input–output analysis to study the flow patterns in transitional plane Couette–Poiseuille flow (CPF). First, we focus on the well-studied intermediate laminar profile, which balances the shear and pressure effects. We show that the highest structured gain corresponds to perturbations with wavelengths associated with the oblique turbulent bands observed in experiments. In addition, the inclination angles of these structures show a Reynolds number dependence consistent with experimentally observed trends. We then examine the Reynolds number scaling of the maximal structured frequency response as the velocity profile is varied from plane Couette to Poiseuille flow. Our results demonstrate that, as expected, the scaling exponent increases over this range, but this increase is not monotonic. We attribute the variation in the shape of the scaling exponent curve as a function of the velocity profile shape to changes in flow patterns that dominate each associated flow regime. We then focus on the particular case of plane Couette flow and compute the structured response modes. Their behavior and structural features are consistent with results obtained through direct numerical simulation (DNS) studies. Finally, we employ the structured analysis framework to examine the temporal evolution of the dominant structures. For the well-studied cases of plane Couette and plane Poiseuille flows, the computed advection speeds of the oblique structures are consistent with those predicted through DNS.

Low-speed, low induction multi-blade rotor for energy efficient small wind turbines

To promote the deployment of small wind turbines (SWTs), thorough understanding of design parameter implications is essential. In-depth research is required to comprehend the influence of design TSR on off-design wind speed performance of multi-blade SWTs. In-house built blade element momentum algorithm was employed, which considered Reynolds number dependence of aerodynamic coefficients and correlated well with experimental results. Peak power coefficients were produced for 0.9, 1.5, 2 and 3 m diameter S826 rotors with blade numbers from 2 to 12 and design TSR of 2–10. Remarkably, regardless of diameter, greatest CP,max values were achieved around design TSR of 4. For peak efficiency, smaller the diameter, narrower the blade number and design TSR range. High-speed rotors have wider TSR range for high power coefficient. Yet, it was shown that operating TSR of low-speed rotors deviates less from design TSR as wind speed varies. It was revealed that low-speed (with a threshold design TSR of 3), low-induction multi-blade rotors provide high CP,max, better efficiency at off-design, shorter starting time and lower wind speed than three-bladed high-speed rotor. A small boost in operational TSR was found to effectively mitigate loss in off-design performance. These are key features to maximize energy harvesting.

Quantitative Prediction of Sling Events in Turbulence at High Reynolds Numbers.

Collisional growth of droplets, such as occurring in warm clouds, is known to be significantly enhanced by turbulence. Whether particles collide depends on their flow history, in particular on their encounters with highly intermittent small-scale turbulent structures, which despite their rarity can dominate the overall collision rate. Here, we develop a quantitative criterion for sling events based on the velocity gradient history along particle paths. We show by a combination of theory and simulations that the problem reduces to a one-dimensional localization problem as encountered in condensed matter physics. The reduction demonstrates that the creation of slings is controlled by the minimal real eigenvalue of the velocity gradient tensor. We use fully resolved turbulence simulations to confirm our predictions and study their Stokes and Reynolds number dependence. We also discuss extrapolations to the parameter range relevant for typical cloud droplets, showing that sling events at high Reynolds numbers are enhanced by an order of magnitude for small Stokes numbers. Thus, intermittency could be a significant ingredient in the collisional growth of rain droplets.

Investigation of strong isothermal stratification effects on multi-mode compressible Rayleigh–Taylor instability

Rayleigh–Taylor instability, RTI, occurs at the interface separating two fluids subjected to acceleration when the density gradient and the acceleration are in opposite directions. Previous scientific research primarily considered RTI under the incompressible assumption, which may not be valid in many high-energy-density engineering applications and astrophysical phenomena. In this study, the compressibility effects of the background isothermal stratification strength on multi-mode two-dimensional RTI are explored using fully compressible multi-species direct numerical simulations. Cases under three different isothermal Mach numbers – Ma=0.15, 0.3, and 0.45 – are investigated to explore weakly, moderately, and strongly stratified compressible RTI, respectively, at an Atwood number of 0.04. Unlike incompressible RTI, an increase in the flow compressibility through the strength of the background stratification can suppress the RTI growth and can lead to a termination of the RTI mixing layer growth with a highly molecularly mixed state. Our findings suggest that even at the chosen relatively low Atwood number, the variable-density effects can be significantly enhanced due to an increase in the background stratification for the compressible RTI as different spatial profiles become noticeably asymmetric across the mixing layer for the strongly stratified case. In addition, this study compares the chaotic behavior of the cases by studying the transport of the turbulent kinetic energy as well as the vortex dynamics. The Reynolds number dependence of the results is also examined with three different Reynolds numbers, and the findings for the large-scale mixing and flow quantities of interest are shown to be universal in the range of the Reynolds numbers studied.

Forecasting small-scale dynamics of fluid turbulence using deep neural networks

Turbulence in fluid flows is characterized by a wide range of interacting scales. Since the scale range increases as some power of the flow Reynolds number, a faithful simulation of the entire scale range is prohibitively expensive at high Reynolds numbers. The most expensive aspect concerns the small-scale motions; thus, major emphasis is placed on understanding and modeling them, taking advantage of their putative universality. In this work, using physics-informed deep learning methods, we present a modeling framework to capture and predict the small-scale dynamics of turbulence, via the velocity gradient tensor. The model is based on obtaining functional closures for the pressure Hessian and viscous Laplacian contributions as functions of velocity gradient tensor. This task is accomplished using deep neural networks that are consistent with physical constraints and explicitly incorporate Reynolds number dependence to account for small-scale intermittency. We then utilize a massive direct numerical simulation database, spanning two orders of magnitude in the large-scale Reynolds number, for training and validation. The model learns from low to moderate Reynolds numbers and successfully predicts velocity gradient statistics at both seen and higher (unseen) Reynolds numbers. The success of our present approach demonstrates the viability of deep learning over traditional modeling approaches in capturing and predicting small-scale features of turbulence.

The hunt for the Kármán ‘constant’ revisited

The log law of the wall, joining the inner, near-wall mean velocity profile (MVP) in wall-bounded turbulent flows to the outer region, has been a permanent fixture of turbulence research for over hundred years, but there is still no general agreement on the value of the prefactor, the inverse of the Kármán ‘constant’ $\kappa$ , or on its universality. The choice diagnostic tool to locate logarithmic parts of the MVP is to look for regions where the indicator function $\varXi$ (equal to the wall-normal coordinate $y^+$ times the mean velocity derivative $\mathrm {d} U^+/\mathrm {d} y^+$ ) is constant. In pressure-driven flows, however, such as channel and pipe flows, $\varXi$ is significantly affected by a term proportional to the wall-normal coordinate, of order $O({Re}_{\tau }^{-1})$ in the inner expansion, but moving up across the overlap to the leading $O(1)$ in the outer expansion. Here we show that, due to this linear overlap term, ${Re}_{\tau }$ values well beyond $10^5$ are required to produce one decade of near constant $\varXi$ in channels and pipes. The problem is resolved by considering the common part of the inner asymptotic expansion carried to $O({Re}_{\tau }^{-1})$ , and the leading order of the outer expansion. This common part contains a superposition of the log law and a linear term $S_0 \,y^+{Re}_{\tau }^{-1}$ , and corresponds to the linear part of $\varXi$ , which, in channel and pipe, is concealed up to $y^+ \approx 500\unicode{x2013}1000$ by terms of the inner expansion. A new and robust method is devised to simultaneously determine $\kappa$ and $S_0$ in pressure-driven flows at currently accessible ${Re}_{\tau }$ values, yielding $\kappa$ values which are consistent with the $\kappa$ values deduced from the Reynolds number dependence of centreline velocities. A comparison with the zero-pressure-gradient turbulent boundary layer, further clarifies the issues and improves our understanding.

Patterns in transitional shear turbulence. Part 2. Emergence and optimal wavelength

Low Reynolds number turbulence in wall-bounded shear flows en route to laminar flow takes the form of oblique, spatially intermittent turbulent structures. In plane Couette flow, these emerge from uniform turbulence via a spatio-temporal intermittent process in which localised quasi-laminar gaps randomly nucleate and disappear. For slightly lower Reynolds numbers, spatially periodic and approximately stationary turbulent–laminar patterns predominate. The statistics of quasi-laminar regions, including the distributions of space and time scales and their Reynolds-number dependence, are analysed. A smooth, but marked transition is observed between uniform turbulence and flow with intermittent quasi-laminar gaps, whereas the transition from gaps to regular patterns is more gradual. Wavelength selection in these patterns is analysed via numerical simulations in oblique domains of various sizes. Via lifetime measurements in minimal domains, and a wavelet-based analysis of wavelength predominance in a large domain, we quantify the existence and nonlinear stability of a pattern as a function of wavelength and Reynolds number. We report that the preferred wavelength maximises the energy and dissipation of the large-scale flow along laminar–turbulent interfaces. This optimal behaviour is due primarily to the advective nature of the large-scale flow, with turbulent fluctuations playing only a secondary role.

Reynolds number dependence of the turbulent/non-turbulent interface in temporally developing turbulent boundary layers

Direct numerical simulations (DNS) of temporally developing turbulent boundary layers are performed with a wide range of Reynolds numbers based on the momentum thickness $Re_{\theta } = 2000$ – $13\,000$ for investigating the Reynolds number dependence of the turbulent/non-turbulent interface (TNTI) layer. The grid spacing in the DNS is determined carefully such that small-scale turbulent motions near the TNTI are well resolved. The outer edge of the TNTI layer, called the irrotational boundary, is detected with vorticity magnitude. The mean thicknesses of the TNTI layer, $\delta _{TNTI}$ , turbulent sublayer, $\delta _{TSL}$ , and viscous superlayer, $\delta _{VSL}$ , are found to be approximately $15\eta _{TI}$ , $10\eta _{TI}$ and $5\eta _{TI}$ , respectively, where $\eta _{TI}$ is the Kolmogorov scale taken in the turbulent region near the TNTI layer. The mean curvature of the irrotational boundary is also characterized by $\eta _{TI}$ . The shear parameter and the shear-to-vorticity ratio show that the mean shear effects near the TNTI layer are not significant for both large and small scales. The anisotropy tensors of Reynolds stress and vorticity suggest that the turbulence under the TNTI layer tends to be isotropic at high $Re_{\theta }$ , for which $\eta _{TI}/\delta \sim Re_{\theta }^{-3/4}$ is valid with the boundary layer thickness $\delta$ . The surface area of the irrotational boundary is consistent with the fractal analysis of the interface, where the fractal dimension $D_f$ is found to be $2.14$ – $2.20$ . The present results suggest that the mean entrainment rate per unit horizontal area normalized by the friction velocity varies slowly as $Re_{\theta }^{(3/4)(D_f-2)}$ for $Re_{\theta } \geq 4000$ .