An official website of the United States government
Idealized large-eddy simulations of shallow convection often utilize horizontally periodic computational domains. The development of precipitation in shallow cumulus convection changes the spatial structure of convection and creates large-scale organization. However, the limited periodic domain constrains the horizontal variability of the atmospheric boundary layer. Small computational domains cannot capture the mesoscale boundary layer organization and artificially constrain the horizontal convection structure. The effects of the horizontal domain size on large-eddy simulations of shallow precipitating cumulus convection are investigated using four computational domains, ranging from 40×40km2 to 320×320km2 and fine grid resolution (40 m). The horizontal variability of the boundary layer is captured in computational domains of 160×160km2. Small LES domains (≤40 km) cannot reproduce the mesoscale flow features, which are about 100km long, but the boundary layer mean profiles are similar to those of the larger domains. Turbulent fluxes, temperature and moisture variances, and horizontal length scales are converged with respect to domain size for domains equal to or larger than 160×160km2. Vertical velocity flow statistics, such as variance and spectra, are essentially identical in all domains and show minor dependence on domain size. Characteristic horizontal length scales (i.e., those relating to the mesoscale organization) of horizontal wind components, temperature and moisture reach an equilibrium after about hour 30.
more »
« less
On the Comparison of Flow Physics between Minimal and Extended Flow Units in Turbulent Channels
Direct numerical simulations were performed to study the effects of the domain size of a minimal flow unit (MFU) and its inherent periodic boundary conditions on flow physics of a turbulent channel flow in a range of 200≤Reτ≤1000. This was accomplished by comparing turbulent statistics with those computed in sub-domains (SD) of extended domain simulations. The dimensions of the MFU and SD were matched, and SD dynamics were set to minimize artificial periodicities. Streamwise and spanwise dimensions of healthy MFUs were found to increase linearly with Reynolds number. It was also found that both MFU and SD statistics and dynamics were healthy and in good agreement. This suggests that healthy MFU dynamics represent extended-domain dynamics well up to Reτ=1000, indicating a nearly negligible effect of periodic conditions on MFUs. However, there was a small deviation within the buffer layer for the MFU at Reτ=200, which manifested in an increased mean velocity and a tail in the Q2 quadrant of the u′-v′ plane. Thus, it should be noted that when considering an MFU domain size, stricter criteria may need to be put in place to ensure healthy turbulent dynamics.
more »
« less
- Award ID(s):
- 1832976
- PAR ID:
- 10296804
- Date Published:
- Journal Name:
- Fluids
- Volume:
- 6
- Issue:
- 5
- ISSN:
- 2311-5521
- Page Range / eLocation ID:
- 192
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Understanding spatial development of a turbulent mixing layer is essential for many engineering applications. However, the flow development is difficult to replicate in physical or numerical experiments. For this reason, the most attractive method for the mixing layer analysis is the direct numerical simulation (DNS), with the most control over the simulation inputs and free from modeling assumptions. On the other hand, the DNS cost often prevents conducting the sensitivity analysis of the simulation results to variations in the numerical procedure and thus, separating numerical and physical effects. In this paper, effects of the computational domain dimensions on statistics collected from DNS of a spatially developing incompressible turbulent mixing layer are analyzed with the focus on determining the domain dimensions suitable for studying the flow asymptotic state. In the simulations, the mixing layer develops between two coflowing laminar boundary layers formed on two sides of a sharp-ended splitter plate of a finite thickness with characteristics close to those of the untripped boundary layers in the experiments by Bell and Mehta, AIAA J., 28(12), 2034 (1990). The simulations were conducted using the spectral-element code Nek5000.more » « less
-
The relation between the helicity and the rate of dissipation of turbulent kinetic energy in turbulent flows has been a matter of debate. Herein, direct numerical simulations of turbulent Poiseuille and Couette flow were used in combination with the tracking of helicity, helicity density, and dissipation along the trajectories of passive scalar markers to probe the correlation between helicity and dissipation in anisotropic turbulence. The Schmidt number of the scalar markers varied between 0.7, 6, and infinite (i.e., fluid particles), while the friction Reynolds number for both simulations was 300. The probing tools were the autocorrelation coefficients, the cross correlation coefficients between helicity and dissipation, and the joint probability density function calculated in the Lagrangian framework along the positions of the scalar markers. These markers were released at different locations within the flow field, including the viscous wall sublayer, the transition layer, the logarithmic region, and the outer flow. In addition, conditional statistics for scalar markers that dispersed most or least in the flow field were also calculated. It was found that helicity and dissipation changed along the trajectories of scalar markers; however, helicity and dissipation were not correlated in the Lagrangian framework. There was anticorrelation between helicity and dissipation in the near wall region, which was less obvious in the logarithmic region. More importantly, helicity could be used to characterize the alignment of the fluctuating velocity and vorticity vectors along the trajectories of scalar markers that disperse the farthest in the direction normal to the channel wall.more » « less
-
In this study, we conduct a parametric analysis to evaluate the sensitivities of wall-modeled large-eddy simulation (LES) with respect to subgrid-scale (SGS) models, mesh resolution, wall boundary conditions and mesh anisotropy. While such investigations have been conducted for attached/flat-plate flow configurations, systematic studies specifically targeting turbulent flows with separation are notably sparse. To bridge this gap, our study focuses on the flow over a two-dimensional Gaussian-shaped bump at a moderately high Reynolds number, which involves smooth-body separation of a turbulent boundary layer under pressure-gradient and surface- curvature effects. In the simulations, the no-slip condition at the wall is replaced by three different forms of boundary condition based on the thin boundary layer equations and the mean wall-shear stress from high-fidelity numerical simulation to avoid the additional complexity of modeling the wall-shear stress. Various statistics, including the mean separation bubble size, mean velocity profile, and dissipation from SGS model, are compared and analyzed. The results reveal that capturing the separation bubble strongly depends on the choice of SGS model. While simulations approach grid convergence with resolutions nearing those of wall-resolved LES meshes, above this limit, the LES predictions exhibit intricate sensitivities to mesh resolution. Furthermore, both wall boundary conditions and the anisotropy of mesh cells exert discernible impacts on the turbulent flow predictions, yet the magnitudes of these impacts vary based on the specific SGS model chosen for the simulation.more » « less
-
High-resolution large eddy simulations and complementary laboratory experiments using particle image velocimetry were performed to provide a detailed quantitative assessment of flow response to gaps in cylinder arrays. The base canopy consists of a dense array of emergent rigid cylinders placed in a regular staggered pattern. The gaps varied in length from [Formula: see text] to 24, in intervals of 4 d, where d is the diameter of the cylinders. The analysis was performed under subcritical conditions with Froude numbers [Formula: see text] and bulk Reynolds numbers [Formula: see text]. Results show that the gaps affect the flow statistics at the upstream and downstream proximity of the canopy. The affected zone was [Formula: see text] for the mean flow and [Formula: see text] for the second-order statistics. Dimensionless time-averaged streamwise velocity within the gap exhibited minor variability with gap spacing; however, in-plane turbulent kinetic energy, k, showed a consistent decay rate when normalized with that at [Formula: see text] from the beginning of the gap. The emergent canopy acts as a passive turbulence generator for the gap flow for practical purposes. The streamwise dependence of k follows an exponential trend within [Formula: see text] and transitions to a power-law at [Formula: see text]. The substantially lower maximum values of k within the gap compared to k within the canopy evidence a limitation of gap measurements representative of canopy flow statistics. We present a base framework for estimating representative in-canopy statistics from measurements in the gap.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/fluids6050192
