The spatial structure of the logarithmic region in very-high-Reynolds-number rough wall turbulent boundary layers

2018 ◽  
Vol 857 ◽  
pp. 704-747 ◽  
Author(s):  
Michael Heisel ◽  
Teja Dasari ◽  
Yun Liu ◽  
Jiarong Hong ◽  
Filippo Coletti ◽  
...  

Using super-large-scale particle image velocimetry (SLPIV), we investigate the spatial structure of the near-wall region in the fully rough atmospheric surface layer with Reynolds number$Re_{\unicode[STIX]{x1D70F}}\sim O(10^{6})$. The field site consists of relatively flat, snow-covered farmland, allowing for the development of a fully rough turbulent boundary layer under near-neutral thermal stability conditions. The imaging field of view extends from 3 m to 19 m above the ground and captures the top of the roughness sublayer and the bottom of an extensive logarithmic region. The SLPIV technique uses natural snowfall as seeding particles for the flow imaging. We demonstrate that SLPIV provides reliable measurements of first- and second-order velocity statistics in the streamwise and wall-normal directions. Our results in the logarithmic region show that the structural features identified in laboratory studies are similarly present in the atmosphere. Using instantaneous vector fields and two-point correlation analysis, we identify vortex structures sharing the signature of hairpin vortex packets. We also evaluate the zonal structure of the boundary layer by tracking uniform momentum zones (UMZs) and the shear interfaces between UMZs in space and time. Statistics of the UMZs and shear interfaces reveal the role of the zonal structure in determining the mean and variance profiles. The velocity difference across the shear interfaces scales with the friction velocity, in agreement with previous studies, and the size of the UMZs scales with wall-normal distance, in agreement with the attached eddy framework.

2012 ◽  
Vol 225 ◽  
pp. 109-117 ◽  
Author(s):  
Zambri Harun ◽  
Mohamad Dali Isa ◽  
Mohammad Rasidi Rasani ◽  
Shahrir Abdullah

Single normal hot-wire measurements of the streamwise component of velocity were taken in boundary layer flows subjected to pressure gradients at matched friction Reynolds numbers Reτ ≈ 3000. To evaluate spatial resolution effects, the sensor lengths are varied in both adverse pressure gradient (APG) and favorable pressure gradient (FPG). A control boundary layer flow in zero pressure gradient ZPG is also presented. It is shown here that, when the sensor length is maintained a constant value, in a contant Reynolds number, the near-wall peak increases with (adverse) pressure gradient. Both increased contributions of the small- and especially large-scale features are attributed to the increased broadband turbulence intensities. The two-mode increase, one centreing in the near-wall region and the other one in the outer region, makes spatial resolution studies in boundary layer flow more complicated. The increased large-scale features in the near-wall region of an APG flow is similar to large-scales increase due to Reynolds number in ZPG flow. Additionally, there is also an increase of the small-scales in the near-wall region when the boundary layer is exposed to adverse pressure gradient (while the Reynolds number is constant). In order to collapse the near-wall peaks for APG, ZPG and FPG flows, the APG flow has to use the longest sensor and conversely, the FPG has to use the shortest sensor. This study recommends that the empirical prediction by Huthins et. al. (2009) to be reevaluated if pressure gradient flows were to be considered such that the magnitude of the near-wall peak is also a function of the adverse pressure gradient parameter.


2011 ◽  
Vol 666 ◽  
pp. 573-604 ◽  
Author(s):  
M. GUALA ◽  
M. METZGER ◽  
B. J. McKEON

Simultaneous streamwise velocity measurements across the vertical direction obtained in the atmospheric surface layer (Reτ ≃ 5 × 105) under near thermally neutral conditions are used to outline and quantify interactions between the scales of turbulence, from the very-large-scale motions to the dissipative scales. Results from conditioned spectra, joint probability density functions and conditional averages show that the signature of very-large-scale oscillations can be found across the whole wall region and that these scales interact with the near-wall turbulence from the energy-containing eddies to the dissipative scales, most strongly in a layer close to the wall, z+ ≲ 103. The scale separation achievable in the atmospheric surface layer appears to be a key difference from the low-Reynolds-number picture, in which structures attached to the wall are known to extend through the full wall-normal extent of the boundary layer. A phenomenological picture of very-large-scale motions coexisting and interacting with structures from the hairpin paradigm is provided here for the high-Reynolds-number case. In particular, it is inferred that the hairpin-packet conceptual model may not be exhaustively representative of the whole wall region, but only of a near-wall layer of z+ = O(103), where scale interactions are mostly confined.


2016 ◽  
Vol 795 ◽  
pp. 356-379 ◽  
Author(s):  
T. Khapko ◽  
P. Schlatter ◽  
Y. Duguet ◽  
D. S. Henningson

Turbulence in the asymptotic suction boundary layer is investigated numerically at the verge of laminarisation using direct numerical simulation. Following an adiabatic protocol, the Reynolds number $Re$ is decreased in small steps starting from a fully turbulent state until laminarisation is observed. Computations in a large numerical domain allow in principle for the possible coexistence of laminar and turbulent regions. However, contrary to other subcritical shear flows, no laminar–turbulent coexistence is observed, even near the onset of sustained turbulence. High-resolution computations suggest a critical Reynolds number $Re_{g}\approx 270$, below which turbulence collapses, based on observation times of $O(10^{5})$ inertial time units. During the laminarisation process, the turbulent flow fragments into a series of transient streamwise-elongated structures, whose interfaces do not display the characteristic obliqueness of classical laminar–turbulent patterns. The law of the wall, i.e. logarithmic scaling of the velocity profile, is retained down to $Re_{g}$, suggesting a large-scale wall-normal transport absent in internal shear flows close to the onset. In order to test the effect of these large-scale structures on the near-wall region, an artificial volume force is added to damp spanwise and wall-normal fluctuations above $y^{+}=100$, in viscous units. Once the largest eddies have been suppressed by the forcing, and thus turbulence is confined to the near-wall region, oblique laminar–turbulent interfaces do emerge as in other wall-bounded flows, however only transiently. These results suggest that oblique stripes at the onset are a prevalent feature of internal shear flows, but will not occur in canonical boundary layers, including the spatially growing ones.


2009 ◽  
Vol 619 ◽  
pp. 79-94 ◽  
Author(s):  
CARLO COSSU ◽  
GREGORY PUJALS ◽  
SEBASTIEN DEPARDON

The optimal energy growth of perturbations sustained by a zero pressure gradient turbulent boundary is computed using the eddy viscosity associated with the turbulent mean flow. It is found that even if all the considered turbulent mean profiles are linearly stable, they support transient energy growths. The most amplified perturbations are streamwise uniform and correspond to streamwise streaks originated by streamwise vortices. For sufficiently large Reynolds numbers two distinct peaks of the optimal growth exist, respectively scaling in inner and outer units. The optimal structures associated with the peak scaling in inner units correspond well with the most probable streaks and vortices observed in the buffer layer, and their moderate energy growth is independent of the Reynolds number. The energy growth associated with the peak scaling in outer units is larger than that of the inner peak and scales linearly with an effective turbulent Reynolds number formed with the maximum eddy viscosity and a modified Rotta–Clauser length based on the momentum thickness. The corresponding optimal perturbations consist of very large–scale structures with a spanwise wavelength of the order of 8δ. The associated optimal streaks scale in outer variables in the outer region and in wall units in the inner region of the boundary layer, in which they are proportional to the mean flow velocity. These outer streaks protrude far into the near wall region, having still 50% of their maximum amplitude at y+ = 20. The amplification of very large–scale structures appears to be a robust feature of the turbulent boundary layer: optimal perturbations with spanwise wavelengths ranging from 4δ to 15δ can all reach 80% of the overall optimal peak growth.


2010 ◽  
Vol 662 ◽  
pp. 409-446 ◽  
Author(s):  
G. SILANO ◽  
K. R. SREENIVASAN ◽  
R. VERZICCO

We summarize the results of an extensive campaign of direct numerical simulations of Rayleigh–Bénard convection at moderate and high Prandtl numbers (10−1 ≤ Pr ≤ 104) and moderate Rayleigh numbers (105 ≤ Ra ≤ 109). The computational domain is a cylindrical cell of aspect ratio Γ = 1/2, with the no-slip condition imposed on all boundaries. By scaling the numerical results, we find that the free-fall velocity should be multiplied by $1/\sqrt{{\it Pr}}$ in order to obtain a more appropriate representation of the large-scale velocity at high Pr. We investigate the Nusselt and the Reynolds number dependences on Ra and Pr, comparing the outcome with previous numerical and experimental results. Depending on Pr, we obtain different power laws of the Nusselt number with respect to Ra, ranging from Ra2/7 for Pr = 1 up to Ra0.31 for Pr = 103. The Nusselt number is independent of Pr. The Reynolds number scales as ${\it Re}\,{\sim}\,\sqrt{{\it Ra}}/{\it Pr}$, neglecting logarithmic corrections. We analyse the global and local features of viscous and thermal boundary layers and their scaling behaviours with respect to Ra and Pr, and with respect to the Reynolds and Péclet numbers. We find that the flow approaches a saturation state when Reynolds number decreases below the critical value, Res ≃ 40. The thermal-boundary-layer thickness increases slightly (instead of decreasing) when the Péclet number increases, because of the moderating influence of the viscous boundary layer. The simulated ranges of Ra and Pr contain steady, periodic and turbulent solutions. A rough estimate of the transition from the steady to the unsteady state is obtained by monitoring the time evolution of the system until it reaches stationary solutions. We find multiple solutions as long-term phenomena at Ra = 108 and Pr = 103, which, however, do not result in significantly different Nusselt numbers. One of these multiple solutions, even if stable over a long time interval, shows a break in the mid-plane symmetry of the temperature profile. We analyse the flow structures through the transitional phases by direct visualizations of the temperature and velocity fields. A wide variety of large-scale circulation and plume structures has been found. The single-roll circulation is characteristic only of the steady and periodic solutions. For other regimes at lower Pr, the mean flow generally consists of two opposite toroidal structures; at higher Pr, the flow is organized in the form of multi-jet structures, extending mostly in the vertical direction. At high Pr, plumes mainly detach from sheet-like structures. The signatures of different large-scale structures are generally well reflected in the data trends with respect to Ra, less in those with respect to Pr.


2013 ◽  
Vol 715 ◽  
pp. 477-498 ◽  
Author(s):  
Zambri Harun ◽  
Jason P. Monty ◽  
Romain Mathis ◽  
Ivan Marusic

AbstractResearch into high-Reynolds-number turbulent boundary layers in recent years has brought about a renewed interest in the larger-scale structures. It is now known that these structures emerge more prominently in the outer region not only due to increased Reynolds number (Metzger & Klewicki, Phys. Fluids, vol. 13(3), 2001, pp. 692–701; Hutchins & Marusic, J. Fluid Mech., vol. 579, 2007, pp. 1–28), but also when a boundary layer is exposed to an adverse pressure gradient (Bradshaw, J. Fluid Mech., vol. 29, 1967, pp. 625–645; Lee & Sung, J. Fluid Mech., vol. 639, 2009, pp. 101–131). The latter case has not received as much attention in the literature. As such, this work investigates the modification of the large-scale features of boundary layers subjected to zero, adverse and favourable pressure gradients. It is first shown that the mean velocities, turbulence intensities and turbulence production are significantly different in the outer region across the three cases. Spectral and scale decomposition analyses confirm that the large scales are more energized throughout the entire adverse pressure gradient boundary layer, especially in the outer region. Although more energetic, there is a similar spectral distribution of energy in the wake region, implying the geometrical structure of the outer layer remains universal in all cases. Comparisons are also made of the amplitude modulation of small scales by the large-scale motions for the three pressure gradient cases. The wall-normal location of the zero-crossing of small-scale amplitude modulation is found to increase with increasing pressure gradient, yet this location continues to coincide with the large-scale energetic peak wall-normal location (as has been observed in zero pressure gradient boundary layers). The amplitude modulation effect is found to increase as pressure gradient is increased from favourable to adverse.


2018 ◽  
Vol 856 ◽  
pp. 958-983 ◽  
Author(s):  
Jinyul Hwang ◽  
Hyung Jin Sung

Wall turbulence is a ubiquitous phenomenon in nature and engineering applications, yet predicting such turbulence is difficult due to its complexity. High-Reynolds-number turbulence arises in most practical flows, and is particularly complicated because of its wide range of scales. Although the attached-eddy hypothesis postulated by Townsend can be used to predict turbulence intensities and serves as a unified theory for the asymptotic behaviours of turbulence, the presence of coherent structures that contribute to the logarithmic behaviours has not been observed in instantaneous flow fields. Here, we demonstrate the logarithmic region of the turbulence intensity by identifying wall-attached structures of the velocity fluctuations ($u_{i}$) through the direct numerical simulation of a moderate-Reynolds-number boundary layer ($Re_{\unicode[STIX]{x1D70F}}\approx 1000$). The wall-attached structures are self-similar with respect to their heights ($l_{y}$), and in particular the population density of the streamwise component ($u$) scales inversely with $l_{y}$, reminiscent of the hierarchy of attached eddies. The turbulence intensities contained within the wall-parallel components ($u$ and $w$) exhibit the logarithmic behaviour. The tall attached structures ($l_{y}^{+}>100$) of $u$ are composed of multiple uniform momentum zones (UMZs) with long streamwise extents, whereas those of the cross-stream components ($v$ and $w$) are relatively short with a comparable width, suggesting the presence of tall vortical structures associated with multiple UMZs. The magnitude of the near-wall peak observed in the streamwise turbulent intensity increases with increasing $l_{y}$, reflecting the nested hierarchies of the attached $u$ structures. These findings suggest that the identified structures are prime candidates for Townsend’s attached-eddy hypothesis and that they can serve as cornerstones for understanding the multiscale phenomena of high-Reynolds-number boundary layers.


Author(s):  
Khaled J. Hammad

Particle Image Velocimetry (PIV) was used to study the flow structure and turbulence, upstream, over, and downstream a shallow open cavity. Three sets of PIV measurements, corresponding to a turbulent incoming boundary layer and a cavity length-to-depth ratio of four, are reported. The cavity depth based Reynolds numbers were 21,000; 42,000; and 54,000. The selected flow configuration and well characterized inflow conditions allow for straightforward assessment of turbulence models and numerical schemes. All mean flow field measurements display a large flow recirculation region, spanning most of the cavity and a smaller, counter-rotating, secondary vortex, immediately downstream of the cavity leading edge. The Galilean decomposed instantaneous velocity vector fields, clearly demonstrate two distinct modes of interaction between the free shear and the cavity trailing edge. The first corresponds to a cascade of vortical structures emanating from the tip of the leading edge of the cavity that grow in size as they travel downstream and directly interact with the trailing edge, i.e., impinging vortices. The second represents vortices that travel above the trailing edge of the cavity, i.e., non-impinging vortices. In the case of impinging vortices, a strong, large scale region of recirculation forms inside the cavity and carries the flow disturbances, arising from the impingement of vortices on the trailing edge of the cavity, upstream in a manner that interacts with and influences the flow as it separates from the cavity leading edge.


2019 ◽  
Vol 875 ◽  
pp. 44-70 ◽  
Author(s):  
Karin Blackman ◽  
Laurent Perret ◽  
Romain Mathis

Urban-type rough-wall boundary layers developing over staggered cube arrays with plan area packing density, $\unicode[STIX]{x1D706}_{p}$, of 6.25 %, 25 % or 44.4 % have been studied at two Reynolds numbers within a wind tunnel using hot-wire anemometry (HWA). A fixed HWA probe is used to capture the outer-layer flow while a second moving probe is used to capture the inner-layer flow at 13 wall-normal positions between $1.25h$ and $4h$ where $h$ is the height of the roughness elements. The synchronized two-point HWA measurements are used to extract the near-canopy large-scale signal using spectral linear stochastic estimation and a predictive model is calibrated in each of the six measurement configurations. Analysis of the predictive model coefficients demonstrates that the canopy geometry has a significant influence on both the superposition and amplitude modulation. The universal signal, the signal that exists in the absence of any large-scale influence, is also modified as a result of local canopy geometry suggesting that although the nonlinear interactions within urban-type rough-wall boundary layers can be modelled using the predictive model as proposed by Mathis et al. (J. Fluid Mech., vol. 681, 2011, pp. 537–566), the model must be however calibrated for each type of canopy flow regime. The Reynolds number does not significantly affect any of the model coefficients, at least over the limited range of Reynolds numbers studied here. Finally, the predictive model is validated using a prediction of the near-canopy signal at a higher Reynolds number and a prediction using reference signals measured in different canopy geometries to run the model. Statistics up to the fourth order and spectra are accurately reproduced demonstrating the capability of the predictive model in an urban-type rough-wall boundary layer.


2011 ◽  
Vol 1 (4) ◽  
Author(s):  
Ziemowit Malecha

AbstractIn this paper, a new computational benchmark test for fluid dynamics is presented. The new benchmark is based on the interaction of a single vortex structure (vortex patch) with a wall. It will be shown that it is possible to distinguish two critical or threshold values of the Reynolds number in the considered flow. The increase of the Reynolds number causes the appearance of the vortex bubble in the near-wall region first, and then next, the eruption of the boundary layer phenomenon. Further increase of the Reynolds number causes the flow to be more complex. The eruption phenomenon becomes more intense and also shows its regenerative nature.


Sign in / Sign up

Export Citation Format

Share Document