Structure Functions and Structure Parameters of Velocity Fluctuations in Numerically Simulated Atmospheric Convective Boundary Layer Flows

2020 ◽  
Vol 77 (10) ◽  
pp. 3619-3630
Author(s):  
Jeremy A. Gibbs ◽  
Evgeni Fedorovich

AbstractWe extend our previous study, which dealt with structure functions of potential temperature fluctuations, and focus on the characteristics of second-order velocity structure functions and corresponding structure parameters in the atmospheric convective boundary layer. We consider the three previously reported methods to compute the structure parameters of turbulent velocity fields: the direct method, the true spectral method, and the approximate spectral method. The methods are evaluated using high-resolution gridded numerical data from large-eddy simulations of shear-free and shear-driven convective boundary layers. Results indicate that the direct and true spectral methods are more suitable than the approximate spectral method, which overestimates the structure parameters of velocity as a result of assuming the inertial-subrange shape of the velocity spectrum for all turbulence scales. Results also suggest that structure parameters of vertical velocity fluctuations are of limited utility because of violations of local isotropy, especially in shear-free convective boundary layers.

2007 ◽  
Vol 64 (4) ◽  
pp. 1230-1248 ◽  
Author(s):  
A. Pier Siebesma ◽  
Pedro M. M. Soares ◽  
João Teixeira

Abstract A better conceptual understanding and more realistic parameterizations of convective boundary layers in climate and weather prediction models have been major challenges in meteorological research. In particular, parameterizations of the dry convective boundary layer, in spite of the absence of water phase-changes and its consequent simplicity as compared to moist convection, typically suffer from problems in attempting to represent realistically the boundary layer growth and what is often referred to as countergradient fluxes. The eddy-diffusivity (ED) approach has been relatively successful in representing some characteristics of neutral boundary layers and surface layers in general. The mass-flux (MF) approach, on the other hand, has been used for the parameterization of shallow and deep moist convection. In this paper, a new approach that relies on a combination of the ED and MF parameterizations (EDMF) is proposed for the dry convective boundary layer. It is shown that the EDMF approach follows naturally from a decomposition of the turbulent fluxes into 1) a part that includes strong organized updrafts, and 2) a remaining turbulent field. At the basis of the EDMF approach is the concept that nonlocal subgrid transport due to the strong updrafts is taken into account by the MF approach, while the remaining transport is taken into account by an ED closure. Large-eddy simulation (LES) results of the dry convective boundary layer are used to support the theoretical framework of this new approach and to determine the parameters of the EDMF model. The performance of the new formulation is evaluated against LES results, and it is shown that the EDMF closure is able to reproduce the main properties of dry convective boundary layers in a realistic manner. Furthermore, it will be shown that this approach has strong advantages over the more traditional countergradient approach, especially in the entrainment layer. As a result, this EDMF approach opens the way to parameterize the clear and cumulus-topped boundary layer in a simple and unified way.


2010 ◽  
Vol 27 (10) ◽  
pp. 1652-1664 ◽  
Author(s):  
Ewan J. O’Connor ◽  
Anthony J. Illingworth ◽  
Ian M. Brooks ◽  
Christopher D. Westbrook ◽  
Robin J. Hogan ◽  
...  

Abstract A method of estimating dissipation rates from a vertically pointing Doppler lidar with high temporal and spatial resolution has been evaluated by comparison with independent measurements derived from a balloon-borne sonic anemometer. This method utilizes the variance of the mean Doppler velocity from a number of sequential samples and requires an estimate of the horizontal wind speed. The noise contribution to the variance can be estimated from the observed signal-to-noise ratio and removed where appropriate. The relative size of the noise variance to the observed variance provides a measure of the confidence in the retrieval. Comparison with in situ dissipation rates derived from the balloon-borne sonic anemometer reveal that this particular Doppler lidar is capable of retrieving dissipation rates over a range of at least three orders of magnitude. This method is most suitable for retrieval of dissipation rates within the convective well-mixed boundary layer where the scales of motion that the Doppler lidar probes remain well within the inertial subrange. Caution must be applied when estimating dissipation rates in more quiescent conditions. For the particular Doppler lidar described here, the selection of suitably short integration times will permit this method to be applicable in such situations but at the expense of accuracy in the Doppler velocity estimates. The two case studies presented here suggest that, with profiles every 4 s, reliable estimates of ε can be derived to within at least an order of magnitude throughout almost all of the lowest 2 km and, in the convective boundary layer, to within 50%. Increasing the integration time for individual profiles to 30 s can improve the accuracy substantially but potentially confines retrievals to within the convective boundary layer. Therefore, optimization of certain instrument parameters may be required for specific implementations.


2011 ◽  
Vol 68 (12) ◽  
pp. 3112-3131 ◽  
Author(s):  
Rachel Honnert ◽  
Valéry Masson ◽  
Fleur Couvreux

Abstract Turbulence is well represented by atmospheric models at very fine grid sizes, from 10 to 100 m, for which turbulent movements are mainly resolved, and by atmospheric models with grid sizes greater than 2 km, for which those movements are entirely parameterized. But what happens at intermediate scales, Wyngaard’s so-called terra incognita? Here an original method is presented that provides a new diagnostic by calculating the subgrid and resolved parts of five variables at different scales: turbulent kinetic energy (TKE), heat and moisture fluxes, and potential temperature and mixing ratio variances. They are established at intermediate scales for dry and cumulus-topped convective boundary layers. The similarity theorem allows the determination of the dimensionless variables of the problem. When the subgrid and resolved parts are studied, a new dimensionless variable, the dimensionless mesh size , needs to be added to the Deardorff free convective scaling variables, where h is the boundary layer height and hc is the height of the cloud layer. Similarity functions for the subgrid and resolved parts are assumed to be the product of the similarity function of the total (subgrid plus resolved) variables and a “partial” similarity function that depends only on . In order to determine the partial similarity function form, large-eddy simulations (LES) of five dry and cloudy convective boundary layers are used. The resolved and subgrid parts of the variables at coarser grid sizes are then deduced from the LES fields. The evolution of the subgrid and resolved parts in the boundary layer with is as follows: fine grids mainly resolve variables. As the mesh becomes coarser, more eddies are subgrid. Finally, for very large meshes, turbulence is entirely subgrid. A scale therefore exists for which the subgrid and resolved parts are equal. This is obtained for in the case of TKE, 0.4 for the potential temperature variance, and 0.8 for the mixing ratio variance, indicating that the velocity structures are smaller than those for the potential temperature, which are smaller than those for the mixing ratio. Furthermore, boundary layers capped by convective clouds have structures larger than dry boundary layer ones as displayed by the scaling in the partial similarity functions. This new diagnostic gives a reference for evaluating current and future parameterizations at kilometric scales. As an illustration, the parameterizations of a mesoscale model are eventually evaluated at intermediate scales. In its standard version, the model produces too many resolved movements, as the turbulence scheme does not sufficiently represent the impact of the subgrid thermal. This is not true when a mass-flux scheme is introduced. However in this case, a completely subgrid thermal is modeled leading to an overestimation of the subgrid part.


2008 ◽  
Vol 8 (6) ◽  
pp. 19637-19677
Author(s):  
R. A. Verzijlbergh ◽  
H. J. J. Jonker ◽  
T. Heus ◽  
J. Vilà-Guerau de Arellano

Abstract. Compared to dry boundary layers, dispersion in cloud-topped boundary layers has received less attention. In this LES based numerical study we investigate the dispersion of a passive tracer in the form of Lagrangian particles for four kinds of atmospheric boundary layers: 1) a dry convective boundary layer (for reference), 2) a "smoke" cloud boundary layer in which the turbulence is driven by radiative cooling, 3) a stratocumulus topped boundary layer and 4) a shallow cumulus topped boundary layer. We show that the dispersion characteristics of the smoke cloud boundary layer as well as the stratocumulus situation can be well understood by borrowing concepts from previous studies of dispersion in the dry convective boundary layer. A general result is that the presence of clouds enhances mixing and dispersion – a notion that is not always reflected well in traditional parameterization models, in which clouds usually suppress dispersion by diminishing solar irradiance. The dispersion characteristics of a cumulus cloud layer turn out to be markedly different from the other three cases and the results can not be explained by only considering the well-known top-hat velocity distribution. To understand the surprising characteristics in the shallow cumulus layer, this case has been examined in more detail by 1) determining the velocity distribution conditioned on the distance to the nearest cloud and 2) accounting for the wavelike behaviour associated with the stratified dry environment.


2016 ◽  
Vol 73 (7) ◽  
pp. 2715-2727 ◽  
Author(s):  
Siddhartha Mukherjee ◽  
Jerôme Schalkwijk ◽  
Harmen J. J. Jonker

Abstract The predictability horizon of convective boundary layers is investigated in this study. Large-eddy simulation (LES) and direct numerical simulation (DNS) techniques are employed to probe the evolution of perturbations in identical twin simulations of a growing dry convective boundary layer. Error growth typical of chaotic systems is observed, marked by two phases. The first comprises an exponential error growth as , with δ0 as the initial error, δ(t) as the error at time t, and Λ as the Lyapunov exponent. This phase is independent of the perturbation wavenumber, and the perturbation energy grows following a self-similar spectral shape dominated by higher wavenumbers. The nondimensional error growth rate in this phase shows a strong dependence on the Reynolds number (Re). The second phase involves saturation of the error. Here, the error growth follows Lorenz dynamics with a slower saturation of successively larger scales. An analysis of the spectral decorrelation times reveals two regimes: an Re-independent regime for scales larger than the boundary layer height and an Re-dependent regime for scales smaller than , which are found to decorrelate substantially faster for increasing Reynolds numbers.


Sign in / Sign up

Export Citation Format

Share Document