scholarly journals Multiscale Shear Forcing of Turbulence in the Nocturnal Boundary Layer: A Statistical Analysis

Author(s):  
Vyacheslav Boyko ◽  
Nikki Vercauteren

AbstractThe lower nocturnal boundary layer is governed by intermittent turbulence which is thought to be triggered by sporadic activity of so-called sub-mesoscale motions in a complex way. We analyze intermittent turbulence based on an assumed relation between the vertical gradients of the sub-mean scales and turbulence kinetic energy. We analyze high-resolution nocturnal eddy-correlation data from 30-m tower collected during the Fluxes over Snow Surfaces II field program. The non-turbulent velocity signal is decomposed using a discrete wavelet transform into three ranges of scales interpreted as the mean, jet and sub-mesoscales. The vertical gradients of the sub-mean scales are estimated using finite differences. The turbulence kinetic energy is modelled as a discrete-time autoregressive process with exogenous variables, where the latter ones are the vertical gradients of the sub-mean scales. The parameters of the discrete model evolve in time depending on the locally-dominant turbulence-production scales. The three regimes with averaged model parameters are estimated using a subspace-clustering algorithm which illustrates a weak bimodal distribution in the energy phase space of turbulence and sub-mesoscale motions for the very stable boundary layer. One mode indicates turbulence modulated by sub-mesoscale motions. Furthermore, intermittent turbulence appears if the sub-mesoscale intensity exceeds $$10 \%$$ 10 % of the mean kinetic energy in strong stratification.

2020 ◽  
Vol 77 (5) ◽  
pp. 1661-1681
Author(s):  
Qingfang Jiang ◽  
Qing Wang ◽  
Shouping Wang ◽  
Saša Gaberšek

Abstract The characteristics of a convective internal boundary layer (CIBL) documented offshore during the East Coast phase of the Coupled Air–Sea Processes and Electromagnetic Ducting Research (CASPER-EAST) field campaign has been examined using field observations, a coupled mesoscale model (i.e., Navy’s COAMPS) simulation, and a couple of surface-layer-resolving large-eddy simulations (LESs). The Lagrangian modeling approach has been adopted with the LES domain being advected from a cool and rough land surface to a warmer and smoother sea surface by the mean offshore winds in the CIBL. The surface fluxes from the LES control run are in reasonable agreement with field observations, and the general CIBL characteristics are consistent with previous studies. According to the LESs, in the nearshore adjustment zone (i.e., fetch < 8 km), the low-level winds and surface friction velocity increase rapidly, and the mean wind profile and vertical velocity skewness in the surface layer deviate substantially from the Monin–Obukhov similarity theory (MOST) scaling. Farther offshore, the nondimensional vertical wind shear and scalar gradients and higher-order moments are consistent with the MOST scaling. An elevated turbulent layer is present immediately below the CIBL top, associated with the vertical wind shear across the CIBL top inversion. Episodic shear instability events occur with a time scale between 10 and 30 min, leading to the formation of elevated maxima in turbulence kinetic energy and momentum fluxes. During these events, the turbulence kinetic energy production exceeds the dissipation, suggesting that the CIBL remains in nonequilibrium.


1980 ◽  
Vol 101 (1) ◽  
pp. 159-178 ◽  
Author(s):  
S. Yavuzkurt ◽  
R. J. Moffat ◽  
W. M. Kays

Hydrodynamic data are reported in the companion paper (Yavuzkurt, Moffat & Kays 1980) for a full-coverage film-cooling situation, both for the blown and the recovery regions. Values of the mean velocity, the turbulent shear stress, and the turbulence kinetic energy were measured at various locations, both within the blown region and in the recovery region. The present paper is concerned with an analysis of the recovery region only. Examination of the data suggested that the recovery-region hydrodynamics could be modelled by considering that a new boundary layer began to grow immediately after the cessation of blowing. Distributions of the Prandtl mixing length were calculated from the data using the measured values of mean velocity and turbulent shear stresses. The mixing-length distributions were consistent with the notion of a dual boundary-layer structure in the recovery region. The measured distributions of mixing length were described by using a piecewise continuous but heuristic fit, consistent with the concept of two quasi-independent layers suggested by the general appearance of the data. This distribution of mixing length, together with a set of otherwise normal constants for a two-dimensional boundary layer, successfully predicted all of the observed features of the flow. The program used in these predictions contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length. The program is a two-dimensional, finite-difference program capable of predicting the mean velocity and turbulence kinetic energy profiles based upon initial values, boundary conditions, and a closure condition.


2019 ◽  
Vol 142 (1) ◽  
Author(s):  
Minghan Chu ◽  
Donald J. Bergstrom

Abstract The present paper reports a numerical study of fully developed turbulent flow over a flat plate with a step change from a smooth to a rough surface. The Reynolds number based on momentum thickness for the smooth flow was Reθ=5950. The focus of the study was to investigate the capability of the Reynolds-averaged Navier–Stokes (RANS) equations to predict the internal boundary layer (IBL) created by the flow configuration. The numerical solution used a two-layer k−ε model to implement the effects of surface roughness on the turbulence and mean flow fields via the use of a hydrodynamic roughness length y0. The prediction for the mean velocity field revealed a development zone immediately downstream of the step in which the mean velocity profile included a lower region affected by the surface roughness below and an upper region with the characteristics of the smooth-wall boundary layer above. In this zone, both the turbulence kinetic energy and Reynolds shear stress profiles were characterized by a significant reduction in magnitude in the outer region of the flow that is unaffected by the rough surface. The turbulence kinetic energy profile was used to estimate the thickness of the IBL, and the resulting growth rate closely matched the experimental results. As such, the IBL is a promising test case for assessing the ability of RANS models to predict the discrete roughness configurations often encountered in industrial and environmental applications.


2007 ◽  
Vol 46 (12) ◽  
pp. 2098-2109 ◽  
Author(s):  
Yansen Wang ◽  
Cheryl L. Klipp ◽  
Dennis M. Garvey ◽  
David A. Ligon ◽  
Chatt C. Williamson ◽  
...  

Abstract Boundary layer wind data observed by a Doppler lidar and sonic anemometers during the mornings of three intensive observational periods (IOP2, IOP3, and IOP7) of the Joint Urban 2003 (JU2003) field experiment are analyzed to extract the mean and turbulent characteristics of airflow over Oklahoma City, Oklahoma. A strong nocturnal low-level jet (LLJ) dominated the flow in the boundary layer over the measurement domain from midnight to the morning hours. Lidar scans through the LLJ taken after sunrise indicate that the LLJ elevation shows a gradual increase of 25–100 m over the urban area relative to that over the upstream suburban area. The mean wind speed beneath the jet over the urban area is about 10%–15% slower than that over the suburban area. Sonic anemometer observations combined with Doppler lidar observations in the urban and suburban areas are also analyzed to investigate the boundary layer turbulence production in the LLJ-dominated atmospheric boundary layer. The turbulence kinetic energy was higher over the urban domain mainly because of the shear production of building surfaces and building wakes. Direct transport of turbulent momentum flux from the LLJ to the urban street level was very small because of the relatively high elevation of the jet. However, since the LLJ dominated the mean wind in the boundary layer, the turbulence kinetic energy in the urban domain is correlated directly with the LLJ maximum speed and inversely with its height. The results indicate that the jet Richardson number is a reasonably good indicator for turbulent kinetic energy over the urban domain in the LLJ-dominated atmospheric boundary layer.


2015 ◽  
Vol 782 ◽  
pp. 333-355 ◽  
Author(s):  
Maarten van Reeuwijk ◽  
John Craske

We discuss energetic restrictions on the entrainment coefficient${\it\alpha}$for axisymmetric jets and plumes. The resulting entrainment relation includes contributions from the mean flow, turbulence and pressure, fundamentally linking${\it\alpha}$to the production of turbulence kinetic energy, the plume Richardson number$\mathit{Ri}$and the profile coefficients associated with the shape of the buoyancy and velocity profiles. This entrainment relation generalises the work by Kaminskiet al. (J. Fluid Mech., vol. 526, 2005, pp. 361–376) and Fox (J. Geophys. Res., vol. 75, 1970, pp. 6818–6835). The energetic viewpoint provides a unified framework with which to analyse the classical entrainment models implied by the plume theories of Mortonet al.(Proc. R. Soc. Lond.A, vol. 234, 1955, pp. 1–23) and Priestley & Ball (Q. J. R. Meteorol. Soc., vol. 81, 1954, pp. 144–157). Data for pure jets and plumes in unstratified environments indicate that to first order the physics is captured by the Priestley and Ball entrainment model, implying that (1) the profile coefficient associated with the production of turbulence kinetic energy has approximately the same value for pure plumes and jets, (2) the value of${\it\alpha}$for a pure plume is roughly a factor of$5/3$larger than for a jet and (3) the enhanced entrainment coefficient in plumes is primarily associated with the behaviour of the mean flow and not with buoyancy-enhanced turbulence. Theoretical suggestions are made on how entrainment can be systematically studied by creating constant-$\mathit{Ri}$flows in a numerical simulation or laboratory experiment.


1998 ◽  
Vol 367 ◽  
pp. 163-203 ◽  
Author(s):  
FARZAD MASHAYEK

Several important issues pertaining to dispersion and polydispersity of droplets in turbulent flows are investigated via direct numerical simulation (DNS). The carrier phase is considered in the Eulerian context, the dispersed phase is tracked in the Lagrangian frame and the interactions between the phases are taken into account in a realistic two-way (coupled) formulation. The resulting scheme is applied for extensive DNS of low-Mach-number, homogeneous shear turbulent flows laden with droplets. Several cases with one- and two-way couplings are considered for both non-evaporating and evaporating droplets. The effects of the mass loading ratio, the droplet time constant, and thermodynamic parameters, such as the droplet specific heat, the droplet latent heat of evaporation, and the boiling temperature, on the turbulence and the droplets are investigated. The effects of the initial droplet temperature and the initial vapour mass fraction in the carrier phase are also studied. The gravity effects are not considered as the numerical methodology is only applicable in the absence of gravity. The evolution of the turbulence kinetic energy and the mean internal energy of both phases is studied by analysing various terms in their transport equations. The results for the non-evaporating droplets show that the presence of the droplets decreases the turbulence kinetic energy of the carrier phase while increasing the level of anisotropy of the flow. The droplet streamwise velocity variance is larger than that of the fluid, and the ratio of the two increases with the increase of the droplet time constant. Evaporation increases both the turbulence kinetic energy and the mean internal energy of the carrier phase by mass transfer. In general, evaporation is controlled by the vapour mass fraction gradient around the droplet when the initial temperature difference between the phases is negligible. In cases with small initial droplet temperature, on the other hand, the convective heat transfer is more important in the evaporation process. At long times, the evaporation rate approaches asymptotic values depending on the values of various parameters. It is shown that the evaporation rate is larger for droplets residing in high-strain-rate regions of the flow, mainly due to larger droplet Reynolds numbers in these regions. For both the evaporating and the non-evaporating droplets, the root mean square (r.m.s.) of the temperature fluctuations of both phases becomes independent of the initial droplet temperature at long times. Some issues relevant to modelling of turbulent flows laden with droplets are also discussed.


2020 ◽  
Vol 77 (7) ◽  
pp. 2605-2626 ◽  
Author(s):  
Bowen Zhou ◽  
Yuhuan Li ◽  
Kefeng Zhu

AbstractBased on a priori analysis of large-eddy simulations (LESs) of the convective atmospheric boundary layer, improved turbulent mixing and dissipation length scales are proposed for a turbulence kinetic energy (TKE)-based planetary boundary layer (PBL) scheme. The turbulent mixing length incorporates surface similarity and TKE constraints in the surface layer, and makes adjustments for lateral entrainment effects in the mixed layer. The dissipation length is constructed based on balanced TKE budgets accounting for shear, buoyancy, and turbulent mixing. A nongradient term is added to the TKE flux to correct for nonlocal turbulent mixing of TKE. The improved length scales are implemented into a PBL scheme, and are tested with idealized single-column convective boundary layer (CBL) cases. Results exhibit robust applicability across a broad CBL stability range, and are in good agreement with LES benchmark simulations. It is then implemented into a community atmospheric model and further evaluated with 3D real-case simulations. Results of the new scheme are of comparable quality to three other well-established PBL schemes. Comparisons between simulated and radiosonde-observed profiles show favorable performance of the new scheme on a clear day.


2017 ◽  
Vol 10 (8) ◽  
pp. 3021-3039 ◽  
Author(s):  
Timothy A. Bonin ◽  
Aditya Choukulkar ◽  
W. Alan Brewer ◽  
Scott P. Sandberg ◽  
Ann M. Weickmann ◽  
...  

Abstract. Measurements of turbulence are essential to understand and quantify the transport and dispersal of heat, moisture, momentum, and trace gases within the planetary boundary layer (PBL). Through the years, various techniques to measure turbulence using Doppler lidar observations have been proposed. However, the accuracy of these measurements has rarely been validated against trusted in situ instrumentation. Herein, data from the eXperimental Planetary boundary layer Instrumentation Assessment (XPIA) are used to verify Doppler lidar turbulence profiles through comparison with sonic anemometer measurements. For 17 days at the end of the experiment, a single scanning Doppler lidar continuously cycled through different turbulence measurement strategies: velocity–azimuth display (VAD), six-beam scans, and range–height indicators (RHIs) with a vertical stare.Measurements of turbulence kinetic energy (TKE), turbulence intensity, and stress velocity from these techniques are compared with sonic anemometer measurements at six heights on a 300 m tower. The six-beam technique is found to generally measure turbulence kinetic energy and turbulence intensity the most accurately at all heights (r2  ≈  0.78), showing little bias in its observations (slope of  ≈  0. 95). Turbulence measurements from the velocity–azimuth display method tended to be biased low near the surface, as large eddies were not captured by the scan. None of the methods evaluated were able to consistently accurately measure the shear velocity (r2 =  0.15–0.17). Each of the scanning strategies assessed had its own strengths and limitations that need to be considered when selecting the method used in future experiments.


2021 ◽  
Author(s):  
Francesco Barbano ◽  
Luigi Brogno ◽  
Francesco Tampieri ◽  
Silvana Di Sabatino

<p>The presence of waves in the nocturnal boundary layer has proven to generate complex interaction with turbulence. On complex terrain environments, where turbulence is observed in a weak but continuous state of activity, waves can be a vehicle of additional production/loss of turbulence energy. The common approach based on the Reynolds decomposition is unable to disaggregate turbulence and wave motion, thus revealing impaired to explicit the terms of this additional interaction. In the current investigation, we adopt a triple-decomposition approach to separate mean, wave, and turbulence motions within near-surface boundary-layer flows, with the aim of unveiling the role of wave motion as source and/or sink of turbulence kinetic and potential energies in the respective explicit budgets. This investigation reveals that the waves contribute to the kinetic energy budget where the production is not shear-dominated and the budget equation does not reduce to a shear-dissipation balance (e.g., as it occurs close to a surface). Away from the surface, the buoyancy effects associated with the wave motion become a significant factor in generating a three-terms balance (shear-buoyancy-dissipation). Similar effects can be found in the potential energy budget, as the waves affect for instance the production associated with the vertical heat flux. On this basis, we develop a simple interpretation paradigm to distinguish two layers, namely near-ground and far-ground sublayer, estimating where the turbulence kinetic energy can significantly feed or be fed by the wave. To prove this paradigm and evaluate the explicit contributions of the wave motion on the turbulence kinetic and potential energies, we investigate a nocturnal valley flow observed under weak synoptic forcing in the Dugway Valley (Utah) during the MATERHORN Program. From this dataset, the explicit kinetic and potential energy budgets are calculated, relying on a variance-covariance analysis to further comprehend the balance of energy production/loss in each sublayer. With this investigation, we propose a simple interpretation scheme to capture and interpret the extent of the complex interaction between waves and turbulence in nocturnal stable boundary layers.</p>


Sign in / Sign up

Export Citation Format

Share Document