Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.
The subgrid-scale (SGS) scalar variance represents the “unmixedness” of the unresolved small scales in large-eddy simulations (LES) of turbulent flows. Supersaturation variance can play an important role in the activation, growth, and evaporation of cloud droplets in a turbulent environment, and therefore efforts are being made to include SGS supersaturation fluctuations in microphysics models. We present results from a priori tests of SGS scalar variance models using data collected in turbulent Rayleigh–Bénard convection in the Michigan Tech Pi chamber for Rayleigh numbers Ra ∼ 108–109. Data from an array of 10 thermistors were spatially filtered and used to calculate the true SGS scalar variance, a scale-similarity model, and a gradient model for dimensionless filter widths of
- Award ID(s):
- 2113060
- NSF-PAR ID:
- 10504008
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- Volume:
- 81
- Issue:
- 5
- ISSN:
- 0022-4928
- Format(s):
- Medium: X Size: p. 839-853
- Size(s):
- p. 839-853
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract The Pi Cloud Chamber offers a unique opportunity to study aerosol‐cloud microphysics interactions in a steady‐state, turbulent environment. In this work, an atmospheric large‐eddy simulation (LES) model with spectral bin microphysics is scaled down to simulate these interactions, allowing comparison with experimental results. A simple scalar flux budget model is developed and used to explore the effect of sidewalls on the bulk mixing temperature, water vapor mixing ratio, and supersaturation. The scaled simulation and the simple scalar flux budget model produce comparable bulk mixing scalar values. The LES dynamics results are compared with particle image velocimetry measurements of turbulent kinetic energy, energy dissipation rates, and large‐scale oscillation frequencies from the cloud chamber. These simulated results match quantitatively to experimental results. Finally, with the bin microphysics included the LES is able to simulate steady‐state cloud conditions and broadening of the cloud droplet size distributions with decreasing droplet number concentration, as observed in the experiments. The results further suggest that collision‐coalescence does not contribute significantly to this broadening. This opens a path for further detailed intercomparison of laboratory and simulation results for model validation and exploration of specific physical processes.
-
more » « less
Abstract Advancing our understanding of astrophysical turbulence is bottlenecked by the limited resolution of numerical simulations that may not fully sample scales in the inertial range. Machine-learning (ML) techniques have demonstrated promise in upscaling resolution in both image analysis and numerical simulations (i.e., superresolution). Here we employ and further develop a physics-constrained convolutional neural network ML model called “MeshFreeFlowNet” (MFFN) for superresolution studies of turbulent systems. The model is trained on both the simulation images and the evaluated partial differential equations (PDEs), making it sensitive to the underlying physics of a particular fluid system. We develop a framework for 2D turbulent Rayleigh–Bénard convection generated with the
Dedalus code by modifying the MFFN architecture to include the full set of simulation PDEs and the boundary conditions. Our training set includes fully developed turbulence sampling Rayleigh numbers (Ra ) ofRa = 106–1010. We evaluate the success of the learned simulations by comparing the power spectra of the directDedalus simulation to the predicted model output and compare both ground-truth and predicted power spectral inertial range scalings to theoretical predictions. We find that the updated network performs well at allRa studied here in recovering large-scale information, including the inertial range slopes. The superresolution prediction is overly dissipative at smaller scales than that of the inertial range in all cases, but the smaller scales are better recovered in more turbulent than laminar regimes. This is likely because more turbulent systems have a rich variety of structures at many length scales compared to laminar flows. -
more » « less
We developed a novel autonomously dynamic nonlocal turbulence model for the large and very large eddy simulation (LES, VLES) of homogeneous isotropic turbulent flows. The model is based on a generalized (integer-to-noninteger)-order Laplacian of the filtered velocity field, and a novel dynamic model has been formulated to avoid the need for tuning the model constant. Three data-driven approaches were introduced for the determination of the fractional-order to have a model that is totally free of any tuning parameter. Our analysis includes both the a priori and the a posteriori tests. In the former test, using a high-fidelity and well-resolved dataset from direct numerical simulations (DNSs), we computed the correlation coefficients for the stress components of the subgrid-scale (SGS) stress tensor and the one we get directly from the DNS results. Moreover, we compared the probability density function of the ensemble-averaged SGS forces for different filter sizes. In the latter, we employed our new model along with other conventional models including the static and dynamic Smagorinsky models into our pseudo-spectral solver and tested the final predicted quantities. The results of the newly developed model exhibit an expressive agreement with the ground-truth DNS results in all components of the SGS stress and forces. Also, the model exhibits promising results in the VLES region as well as the LES region, which could be remarkably important for cost-efficient nonlocal turbulence modeling, e.g., in meteorological and environmental applications.
-
more » « less
In this study, we revisit the spectral transfer model for the turbulent intensity in passive scalar transport (under large-scale anisotropic forcing), and a subsequent modification to the scaling of scalar variance cascade is presented. From the modified spectral transfer model, we obtain a revised scalar transport model using a fractional-order Laplacian operator that facilitates the robust inclusion of the non-local effects originating from large-scale anisotropy transferred across the multitude of scales in the turbulent cascade. We provide an
a priori estimate for the non-local model based on the scaling analysis of the scalar spectrum, and later examine our developed model through direct numerical simulation. We present a detailed analysis on the evolution of the scalar variance, high-order statistics of the scalar gradient and important two-point statistical metrics of the turbulent transport to make a comprehensive comparison between the non-local model and its standard version. Finally, we present an analysis that seamlessly reconciles the similarities between the developed model with the fractional-order subgrid-scale scalar flux model for large-eddy simulation (Akhavan-Safaeiet al. ,J. Comput. Phys. , vol. 446, 2021, 110571) when the filter scale approaches the dissipative scales of turbulent transport. In order to perform this task, we employ a Gaussian process regression model to predict the model coefficient for the fractional-order subgrid model.
