Fine-scale turbulent bursts in stable atmospheric boundary layer in complex terrain

2017 ◽  
Vol 833 ◽  
pp. 745-772 ◽  
Author(s):  
E. Kit ◽  
C. M. Hocut ◽  
D. Liberzon ◽  
H. J. S. Fernando

Turbulence in the atmospheric boundary layer (ABL) is usually measured using sonic anemometers (sonics), but coarse spatial (${\sim}10$  cm) and temporal (${\sim}32$  Hz) resolutions of sonics preclude direct measurement of fine-scale parameters such as the turbulent kinetic energy (TKE) dissipation rate $\unicode[STIX]{x1D700}$. Instead, $\unicode[STIX]{x1D700}$ is estimated using techniques based on Kolmogorov theory. Fine-scale measurements of ABL turbulence down to Kolmogorov scale were made with a sonic and hot-film anemometer dyad (a ‘combo’ probe) during the field campaigns of the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) programme. The hot-film probe was located on a gimbal within the sonic probe volume, and was automated to rotate in the horizontal plane to align with the mean flow measured by sonic. This procedure not only helped satisfy the requirement of hot-film alignment with the mean flow, but also allowed in situ calibration of hot-film probes. This paper analyses a period of nocturnal flow that was similar to a stratified parallel shear flow. The combo-probe measurements showed an interesting phenomenon – the occurrence of strong bursts, characterized by short-term increase of velocity fluctuations and simultaneous increase of TKE dissipation rate by orders of magnitude. These bursts were indicative of unusual turbulence activity at finer (${\sim}0.1$–0.4 m) scales that are not captured by sonics since the smallest scales resolved by the latter are greater than 0.6 m. With bursting present, the spectra exhibited bumps at scales intermediate to inertial and dissipation subranges, resembling a bottleneck phenomenon. Its manifestation, although unequivocally related to bursts, may not convincingly fit into the framework of previous bottleneck-effect theories that allude to either viscous effects or buoyancy effects modifying the local energy cascade via non-local effects. The origins of burst are yet to be identified. Stratified ABL with bursts exhibits non-Kolmogorov behaviour, and hence should be modelled with caution.

2014 ◽  
Vol 755 ◽  
pp. 654-671 ◽  
Author(s):  
Enrico Deusebio ◽  
Erik Lindborg

AbstractHelicity, which is defined as the scalar product of velocity and vorticity, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}{\mathcal{H}} = {\boldsymbol {u}} \boldsymbol {\cdot }\boldsymbol{\omega}$, is an inviscidly conserved quantity in a barotropic fluid. Mean helicity is zero in flows that are parity invariant. System rotation breaks parity invariance and has therefore the potential of giving rise to non-zero mean helicity. In this paper we study the helicity dynamics in the incompressible Ekman boundary layer. Evolution equations for the mean field helicity and the mean turbulent helicity are derived and it is shown that pressure flux injects helicity at a rate $ 2 \varOmega G^2 $ over the total depth of the Ekman layer, where $ G $ is the geostrophic wind far from the wall and $ {\boldsymbol{\Omega}} = \varOmega {\boldsymbol {e}}_y $ is the rotation vector and $ {\boldsymbol {e}}_y $ is the wall-normal unit vector. Thus right-handed/left-handed helicity will be injected if $ \varOmega $ is positive/negative. We also show that in the uppermost part of the boundary layer there is a net helicity injection with opposite sign as compared with the totally integrated injection. Isotropic relations for the helicity dissipation and the helicity spectrum are derived and it is shown that it is sufficient to measure two transverse velocity components and use Taylor’s hypothesis in the mean flow direction in order to measure the isotropic helicity spectrum. We compare the theoretical predictions with a direct numerical simulation of an Ekman boundary layer and confirm that there is a preference for right-handed helicity in the lower part of the Ekman layer and left-handed helicity in the uppermost part when $ \varOmega > 0 $. In the logarithmic range, the helicity dissipation conforms to isotropic relations. On the other hand, spectra show significant departures from isotropic conditions, suggesting that the Reynolds number considered in the study is not sufficiently large for isotropy to be valid in a wide range of scales. Our analytical and numerical results strongly suggest that there is a turbulent helicity cascade of right-handed helicity in the logarithmic range of the atmospheric boundary layer when $\varOmega >0$, consistent with recent measurements by Koprov, Koprov, Ponomarev & Chkhetiani (Dokl. Phys., vol. 50, 2005, pp. 419–422). The isotropic relations which are derived may facilitate future measurements of the helicity spectrum in the atmospheric boundary layer as well as in controlled wind tunnel experiments.


Author(s):  
D. Furey ◽  
P. Atsavapranee ◽  
K. Cipolla

Stereo Particle Image velocimetry data was collected over high aspect ratio flexible cylinders (L/a = 1.5 to 3 × 105) to evaluate the axial development of the turbulent boundary layer where the boundary layer thickness becomes significantly larger than the cylinder diameter (δ/a>>1). The flexible cylinders are approximately neutrally buoyant and have an initial length of 152 m and radii of 0.45 mm and 1.25 mm. The cylinders were towed at speeds ranging from 3.8 to 15.4 m/sec in the David Taylor Model Basin. The analysis of the SPIV data required a several step procedure to evaluate the cylinder boundary flow. First, the characterization of the flow field from the towing strut is required. This evaluation provides the residual mean velocities and turbulence levels caused by the towing hardware at each speed and axial location. These values, called tare values, are necessary for comparing to the cylinder flow results. Second, the cylinder flow fields are averaged together and the averaged tare fields are subtracted out to remove strut-induced ambient flow effects. Prior to averaging, the cylinder flow fields are shifted to collocate the cylinder within the field. Since the boundary layer develops slowly, all planes of data occurring within each 10 meter increment of the cylinder length are averaged together to produce the mean boundary layer flow. Corresponding fields from multiple runs executed using the same experimental parameters are also averaged. This flow is analyzed to evaluate the level of axisymmetry in the data and determine if small changes in cylinder angle affect the mean flow development. With axisymmetry verified, the boundary flow is further averaged azimuthally around the cylinder to produce mean boundary layer profiles. Finally, the fluctuating velocity levels are evaluated for the flow with the cylinder and compared to the fluctuating velocity levels in the tare data. This paper will first discuss the data analysis techniques for the tare data and the averaging methods implemented. Second, the data analysis considerations will be presented for the cylinder data and the averaging and cylinder tracking techniques. These results are used to extract relevant boundary layer parameters including δ, δ* and θ. Combining these results with wall shear and momentum thickness values extracted from averaged cylinder drag data, the boundary layer can be well characterized.


2021 ◽  
Vol 929 ◽  
Author(s):  
N. Agastya Balantrapu ◽  
Christopher Hickling ◽  
W. Nathan Alexander ◽  
William Devenport

Experiments were performed over a body of revolution at a length-based Reynolds number of 1.9 million. While the lateral curvature parameters are moderate ( $\delta /r_s < 2, r_s^+>500$ , where $\delta$ is the boundary layer thickness and r s is the radius of curvature), the pressure gradient is increasingly adverse ( $\beta _{C} \in [5 \text {--} 18]$ where $\beta_{C}$ is Clauser’s pressure gradient parameter), representative of vehicle-relevant conditions. The mean flow in the outer regions of this fully attached boundary layer displays some properties of a free-shear layer, with the mean-velocity and turbulence intensity profiles attaining self-similarity with the ‘embedded shear layer’ scaling (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642). Spectral analysis of the streamwise turbulence revealed that, as the mean flow decelerates, the large-scale motions energize across the boundary layer, growing proportionally with the boundary layer thickness. When scaled with the shear layer parameters, the distribution of the energy in the low-frequency region is approximately self-similar, emphasizing the role of the embedded shear layer in the large-scale motions. The correlation structure of the boundary layer is discussed at length to supply information towards the development of turbulence and aeroacoustic models. One major finding is that the estimation of integral turbulence length scales from single-point measurements, via Taylor's hypothesis, requires significant corrections to the convection velocity in the inner 50 % of the boundary layer. The apparent convection velocity (estimated from the ratio of integral length scale to the time scale), is approximately 40 % greater than the local mean velocity, suggesting the turbulence is convected much faster than previously thought. Closer to the wall even higher corrections are required.


2017 ◽  
Vol 17 (1) ◽  
pp. 531-549 ◽  
Author(s):  
Sanjay Kumar Mehta ◽  
Madineni Venkat Ratnam ◽  
Sukumarapillai V. Sunilkumar ◽  
Daggumati Narayana Rao ◽  
Boddapaty V. Krishna Murthy

Abstract. The diurnal variation of atmospheric boundary layer (ABL) height is studied using high-resolution radiosonde observations available at 3 h intervals for 3 days continuously from 34 intensive campaigns conducted during the period December 2010–March 2014 over a tropical station Gadanki (13.5° N, 79.2° E; 375 m), in the Indian monsoon region. The heights of the ABL during the different stages of its diurnal evolution, namely, the convective boundary layer (CBL), the stable boundary layer (SBL), and the residual layer (RL) are obtained to study the diurnal variabilities. A clear diurnal variation is observed in 9 campaigns out of the 34 campaigns. In 7 campaigns the SBL did not form in the entire day and in the remaining 18 campaigns the SBL formed intermittently. The SBL forms for 33–55 % of the time during nighttime and 9 and 25 % during the evening and morning hours, respectively. The mean SBL height is within 0.3 km above the surface which increases slightly just after midnight (02:00 IST) and remains almost constant until the morning. The mean CBL height is within 3.0 km above the surface, which generally increases from morning to evening. The mean RL height is within 2 km above the surface which generally decreases slowly as the night progresses. The diurnal variation of the ABL height over the Indian region is stronger during the pre-monsoon and weaker during winter season. The CBL is higher during the summer monsoon and lower during the winter season while the RL is higher during the winter season and lower during the summer season. During all the seasons, the ABL height peaks during the afternoon (∼ 14:00 IST) and remains elevated until evening (∼ 17:00 IST). The ABL suddenly collapses at 20:00 IST and increases slightly in the night. Interestingly, it is found that the low level clouds have an effect on the ABL height variability, but the deep convective clouds do not. The lifting condensation level (LCL) is generally found to occur below the ABL for the majority of the database and they are randomly related.


Author(s):  
Sahib Singh Chawla

The laminar boundary layer on a magnetized plate, when the magnetic field oscillates in magnitude about a constant non-zero mean, is analysed. For low-frequency fluctuations the solution is obtained by a series expansion in terms of a frequency parameter, while for high frequencies the flow pattern is of the ‘skin-wave’ type unaffected by the mean flow. In the low-frequency range, the phase lead and the amplitude of the skin-friction oscillations increase at first and then decrease to their respective ‘skin-wave’ values. On the other hand the phase angle of the surface current decreases from 90° to 45° and its amplitude increases with frequency.


1997 ◽  
Vol 119 (2) ◽  
pp. 277-280 ◽  
Author(s):  
B. A. Singer

Models for the distribution of the wall-pressure under a turbulent boundary layer often estimate the coherence of the cross-spectral density in terms of a product of two coherence functions. One such function describes the coherence as a function of separation distance in the mean-flow direction, the other function describes the coherence in the cross-stream direction. Analysis of data from a large-eddy simulation of a turbulent boundary layer reveals that this approximation dramatically underpredicts the coherence for separation directions that are neither aligned with nor perpendicular to the mean-flow direction. These models fail even when the coherence functions in the directions parallel and perpendicular to the mean flow are known exactly. A new approach for combining the parallel and perpendicular coherence functions is presented. The new approach results in vastly improved approximations for the coherence.


1994 ◽  
Vol 116 (4) ◽  
pp. 586-596 ◽  
Author(s):  
P. L. Andrew ◽  
Wing-fai Ng

The turbulent character of the supersonic wake of a linear cascade of fan airfoils has been studied using a two-component laser-doppler anemometer. The cascade was tested in the Virginia Polytechnic Institute and State University intermittent wind tunnel facility, where the Mach and Reynolds numbers were 2.36 and 4.8 × 106, respectively. In addition to mean flow measurements, Reynolds normal and shear stresses were measured as functions of cascade incidence angle and streamwise locations spanning the near-wake and the far-wake. The extremities of profiles of both the mean and turbulent wake properties´ were found to be strongly influenced by upstream shock-boundary -layer interactions, the strength of which varied with cascade incidence. In contrast, the peak levels of turbulence properties within the shear layer were found to be largely independent of incidence, and could be characterized in terms of the streamwise position only. The velocity defect turbulence level was found to be 23 percent, and the generally accepted value of the turbulence structural coefficient of 0.30 was found to be valid for this flow. The degree of similarity of the mean flow wake profiles was established, and those profiles demonstrating the most similarity were found to approach a state of equilibrium between the mean and turbulent properties. In general, this wake flow may be described as a classical free shear flow, upon which the influence of upstream shock-boundary-layer interactions has been superimposed.


2007 ◽  
Vol 64 (9) ◽  
pp. 3363-3371 ◽  
Author(s):  
François Lott

Abstract The backward reflection of a stationary gravity wave (GW) propagating toward the ground is examined in the linear viscous case and for large Reynolds numbers (Re). In this case, the stationary GW presents a critical level at the ground because the mean wind is null there. When the mean flow Richardson number at the surface (J) is below 0.25, the GW reflection by the viscous boundary layer is total in the inviscid limit Re → ∞. The GW is a little absorbed when Re is finite, and the reflection decreases when both the dissipation and J increase. When J &gt; 0.25, the GW is absorbed for all values of the Reynolds number, with a general tendency for the GW reflection to decrease when J increases. As a large ground reflection favors the downstream development of a trapped lee wave, the fact that it decreases when J increases explains why the more unstable boundary layers favor the onset of mountain lee waves. It is also shown that the GW reflection when J &gt; 0.25 is substantially larger than that predicted by the conventional inviscid critical level theory and larger than that predicted when the dissipations are represented by Rayleigh friction and Newtonian cooling. The fact that the GW reflection depends strongly on the Richardson number indicates that there is some correspondence between the dynamics of trapped lee waves and the dynamics of Kelvin–Helmholtz instabilities. Accordingly, and in one classical example, it is shown that some among the neutral modes for Kelvin–Helmholtz instabilities that exist in an unbounded flow when J &lt; 0.25 can also be stationary trapped-wave solutions when there is a ground and in the inviscid limit Re → ∞. When Re is finite, these solutions are affected by the dissipation in the boundary layer and decay in the downstream direction. Interestingly, their decay rate increases when both the dissipation and J increase, as does the GW absorption by the viscous boundary layer.


2019 ◽  
Vol 865 ◽  
pp. 1085-1109 ◽  
Author(s):  
Yutaro Motoori ◽  
Susumu Goto

To understand the generation mechanism of a hierarchy of multiscale vortices in a high-Reynolds-number turbulent boundary layer, we conduct direct numerical simulations and educe the hierarchy of vortices by applying a coarse-graining method to the simulated turbulent velocity field. When the Reynolds number is high enough for the premultiplied energy spectrum of the streamwise velocity component to show the second peak and for the energy spectrum to obey the$-5/3$power law, small-scale vortices, that is, vortices sufficiently smaller than the height from the wall, in the log layer are generated predominantly by the stretching in strain-rate fields at larger scales rather than by the mean-flow stretching. In such a case, the twice-larger scale contributes most to the stretching of smaller-scale vortices. This generation mechanism of small-scale vortices is similar to the one observed in fully developed turbulence in a periodic cube and consistent with the picture of the energy cascade. On the other hand, large-scale vortices, that is, vortices as large as the height, are stretched and amplified directly by the mean flow. We show quantitative evidence of these scale-dependent generation mechanisms of vortices on the basis of numerical analyses of the scale-dependent enstrophy production rate. We also demonstrate concrete examples of the generation process of the hierarchy of multiscale vortices.


Sign in / Sign up

Export Citation Format

Share Document