Engineering Papers⌕ Search

Engineering topics

Ovchinnikov, Mikhail

Publications and source records attributed to Ovchinnikov, Mikhail.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Microphysics regimes due to haze–cloud interactions: cloud oscillation and cloud collapse

Abstract. It is known that aqueous haze particles can be activated into cloud droplets in a supersaturated environment. However, haze–cloud interactions have not been fully explored, partly because haze particles are not represented in most cloud-resolving models. Here, we conduct a series of large-eddy simulations (LESs) of a cloud in a convection chamber using a haze-capable Eulerian-based bin microphysics scheme to explore haze–cloud interactions over a wide range of aerosol injection rates. Results show that the cloud is in a slow microphysics regime at low aerosol injection rates, where the cloud responds slowly to an environmental change and droplet deactivation is negligible. The cloud is in a fast microphysics regime at moderate aerosol injection rates, where the cloud responds quickly to an environmental change and haze–cloud interactions are important. More interestingly, two more microphysics regimes are observed at high aerosol injection rates due to haze–cloud interactions. Cloud oscillation is driven by the oscillation of the mean supersaturation around the critical supersaturation of aerosol due to haze–cloud interactions. Cloud collapse happens under weaker forcing of supersaturation where the chamber transfers cloud droplets to haze particles efficiently, leading to a significant decrease (collapse) in cloud droplet number concentration. One special case of cloud collapse is the haze-only regime. It occurs at extremely high aerosol injection rates, where droplet activation is inhibited, and the sedimentation of haze particles is balanced by the aerosol injection rate. Our results suggest that haze particles and their interactions with cloud droplets should be considered, especially in polluted conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of Autoconversion Representation in E3SMv2 Using an Ensemble of Large-Eddy Simulations of Low-Level Warm Clouds

In numerical atmospheric models that treat cloud and rain droplet populations as separate condensate categories, precipitation initiation in warm clouds is often represented by an autoconversion rate (Au), which is the rate of formation of new rain droplets through the collisions of cloud droplets. Being a function of the cloud droplet size distribution (DSD), the local Au is commonly parameterized as a function of DSD moments: cloud droplet number (n c ) and mass (q c ) concentrations. When applied in a large-scale model, the grid-mean Au must also include a correction, or enhancement factor, to account for the horizontal variability of the cloud properties across the model grid. In this study, we evaluate the Au representation in the Energy Exascale Earth System Model version 2 (E3SMv2) climate model using large-eddy simulations (LES), which explicitly resolve cloud droplet spectra, and therefore the local Au, as well as its spatial variability. The analysis of an ensemble of warm low-level cloud cases shows that the E3SMv2 formulation represents the Au reasonably well compared to the horizontally averaged explicitly computed rate from LES. The agreement, however, comes from a combination of an underestimated E3SM-tuned local Au rate and an overestimated subgrid cloud variability enhancement factor. The latter bias is traced to neglecting the horizontal variability of n c and its co-variability with q c in parameterizing the grid-mean Au.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

EAGLES Liquid Cloud Testbed Large Eddy Simulation Library (v2)

This library consists of large eddy simulation (LES) model output using the PINACLES codebase coupled to the Hebrew University Fast Spectral Bin Microphysics scheme representing shallow, liquid phase clouds from a range of global liquid cloud testbed regions as well as some well known LES model intercomparison cases. This data is particularly suitable for examining microphysical modeling assumptions in coarser-scale models, and was used for this purpose in the work "“Evaluation of Autoconversion Representation in E3SM.v2 using an Ensemble of Large-Eddy Simulations of Low-Level Warm Clouds” by M. Ovchinnikov, P.-L. Ma, C. M. Kaul, K. G. Pressel, M. Huang, J. Shpund, and S. Tang

Kaul, Colleen M↗
DOE OSTI2024

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Designing a Convection‐Cloud Chamber for Collision‐Coalescence Using Large‐Eddy Simulation With Bin Microphysics

Abstract Collisional growth of cloud droplets is an essential yet uncertain process for drizzle and precipitation formation. To improve the quantitative understanding of this key component of cloud‐aerosol‐turbulence interactions, observational studies of collision‐coalescence in a controlled laboratory environment are needed. In an existing convection‐cloud chamber (the Pi Chamber), collisional growth is limited by low liquid water content and short droplet residence times. In this work, we use numerical simulations to explore various configurations of a convection‐cloud chamber that may intensify collision‐coalescence. We employ a large‐eddy simulation (LES) model with a size‐resolved (bin) cloud microphysics scheme to explore how cloud properties and the intensity of collision‐coalescence are affected by the chamber size and aspect ratio, surface roughness, side‐wall wetness, side‐wall temperature arrangement, and aerosol injection rate. Simulations without condensation and evaporation within the domain are first performed to explore the turbulence dynamics and wall fluxes. The LES wall fluxes are used to modify the Scalar Flux‐budget Model, which is then applied to demonstrate the need for non‐uniform side‐wall temperature (two side walls as warm as the bottom and the two others as cold as the top) to maintain high supersaturation in a tall chamber. The results of LES with full cloud microphysics reveal that collision‐coalescence is greatly enhanced by employing a taller chamber with saturated side walls, non‐uniform side‐wall temperature, and rough surfaces. For the conditions explored, although lowering the aerosol injection rate broadens the droplet size distribution, favoring collision‐coalescence, the reduced droplet number concentration decreases the frequency of collisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Enhancements in Cloud Condensation Nuclei Activity From Turbulent Fluctuations in Supersaturation

Abstract The effect of aerosols on the properties of clouds is a large source of uncertainty in predictions of weather and climate. These aerosol‐cloud interactions depend critically on the ability of aerosol particles to form cloud droplets. A challenge in modeling aerosol‐cloud interactions is the representation of interactions between turbulence and cloud microphysics. Turbulent mixing leads to small‐scale fluctuations in water vapor and temperature that are unresolved in large‐scale atmospheric models. To quantify the impact of turbulent fluctuations on cloud condensation nuclei (CCN) activation, we used a high‐resolution Large Eddy Simulation of a convective cloud chamber to drive particle‐based cloud microphysics simulations. We show small‐scale fluctuations strongly impact CCN activity. Once activated, the relatively long timescales of evaporation compared to fluctuations causes droplets to persist in subsaturated regions, which further increases droplet concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating Shallow Convection Parameterization Assumptions With a q t – w Quadrant Analysis

Abstract Uncertainties associated with the response of shallow clouds to global warming remain a great challenge for climate projection. Due to the small size of these clouds, parameterizations are required to represent them in both current‐ and next‐generation climate models. We present a quantitative evaluation of several important assumptions used in both mass‐flux and Assumed‐PDF Higher‐Order Closure parameterizations of shallow convection. We use large‐eddy simulations (LESs) of four different shallow convection regimes as benchmarks and apply a q t (total moisture)– w (vertical velocity) quadrant analysis technique to identify the “coherent structures” of moist and dry up/down‐drafts. The statistics of these coherent structures and the environment are then used to evaluate commonly used assumptions in mass‐flux parameterizations related to (a) the convective downdrafts and (b) the subplume variability of convective updrafts. For quantitative evaluation of the PDF closure used in AP‐HOCs, we perform offline calculations with the PDF closure used in the Cloud Layers Unified By Bi‐normals and the Simplified Higher‐Order Closure schemes by supplying the closure with statistical quantities directly calculated from the LESs. Then, the impact of parameterization assumptions embedded in the PDF closure in the same two categories is analyzed and understood through comparisons of the output from the closure with the statistics of the moist and dry up/down‐drafts from the LESs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An Intercomparison of Large‐Eddy Simulations of a Convection Cloud Chamber Using Haze‐Capable Bin and Lagrangian Cloud Microphysics Schemes

Abstract Recent in situ observations show that haze particles exist in a convection cloud chamber. The microphysics schemes previously used for large‐eddy simulations of the cloud chamber could not fully resolve haze particles and the associated processes, including their activation and deactivation. Specifically, cloud droplet activation was modeled based on Twomey‐type parameterizations, wherein cloud droplets were formed when a critical supersaturation for the available cloud condensation nuclei (CCN) was exceeded and haze particles were not explicitly resolved. Here, we develop and adapt haze‐capable bin and Lagrangian microphysics schemes to properly resolve the activation and deactivation processes. Results are compared with the Twomey‐type CCN‐based bin microphysics scheme in which haze particles are not fully resolved. We find that results from the haze‐capable bin microphysics scheme agree well with those from the Lagrangian microphysics scheme. However, both schemes significantly differ from those from a CCN‐based bin microphysics scheme unless CCN recycling is considered. Haze particles from the recycling of deactivated cloud droplets can strongly enhance cloud droplet number concentration due to a positive feedback in haze‐cloud interactions in the cloud chamber. Haze particle size distributions are more realistic when considering solute and curvature effects that enable representing the complete physics of the activation process. Our study suggests that haze particles and their interactions with cloud droplets may have a strong impact on cloud properties when supersaturation fluctuations are comparable to mean supersaturation, as is the case in the cloud chamber and likely is the case in the atmosphere, especially in polluted conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Scaling of Turbulence and Microphysics in a Convection–Cloud Chamber of Varying Height

Abstract The convection–cloud chamber enables measurement of aerosol and cloud microphysics, as well as their interactions, within a turbulent environment under steady‐state conditions. Increasing the size of a convection–cloud chamber, while holding the imposed temperature difference constant, leads to increased Rayleigh, Reynolds and Nusselt numbers. Large–eddy simulation coupled with a bin microphysics model allows the influence of increased velocity, time, and spatial scales on cloud microphysical properties to be explored. Simulations of a convection–cloud chamber, with fixed aspect ratio and increasing heights of H = 1, 2, 4, and (for dry conditions only) 8 m are performed. The key findings are: Velocity fluctuations scale as H 1/3 , consistent with the Deardorff expression for convective velocity, and implying that the turbulence correlation time scales as H 2/3 . Temperature and other scalar fluctuations scale as H −3/7 . Droplet size distributions from chambers of different sizes can be matched by adjusting the total aerosol injection rate as the horizontal cross‐sectional area (i.e., as H 2 for constant aspect ratio). Injection of aerosols at a point versus distributed throughout the volume makes no difference for polluted conditions, but can lead to cloud droplet size distribution broadening in clean conditions. Cloud droplet growth by collision and coalescence leads to a broader right tail of the distribution compared to condensation growth alone, and this tail increases in magnitude and extent monotonically as the increase of chamber height. These results also have implications for scaling within turbulent, cloudy mixed‐layers in the atmosphere, such as fog layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling Land-Atmosphere Coupling at Cloud-Resolving Scale Within the Multiple Atmosphere Multiple Land (MAML) Framework in SP-E3SM

Representing subgrid variabilities of land surface processes and their upscaled effects is crucial for global climate modeling. Here, we implement a Multiple Atmosphere Multiple Land (MAML) framework in the superparamaterized version of E3SM (SP32 E3SM) to explicitly simulate the subgrid variabilities of land states and fluxes at cloud-resolving scale and their interactions with atmosphere. Comparing to the standard SP-E3SM in which all the atmospheric columns of the cloud resolving model embedded within the global atmospheric model grid interact with the same land surface (i.e., Multiple Atmosphere Single Land (MASL)), the impact of MAML on the strength of land-atmosphere coupling is limited, partly because the current implementation mainly facilitates one-way coupling between the cloud-resolving model and the land surface model. Despite such limitation, MAML increases the surface latent heat flux at the expense of sensible heat flux, and increases precipitation in India, Amazon, and Central Africa, reducing the model dry bias compared to the standard SP-E3SM. By employing a normalized gross moist stability (NGMS) diagnostic framework, we find that the increase in precipitation minus evaporation (PE) is primarily driven by the change in large-scale moisture convergence, particularly by the increase of water vapor in the lower atmosphere, while the local effect of total surface energy flux plays a minor role in the P-E change. More specifically, MAML changes the surface energy partitioning (evaporative fraction), increases the atmosphere water vapor, and further increases P-E by decreasing the NGMS. Finally, future development in the MAML framework is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Checkerboard patterns in E3SMv2 and E3SM-MMFv2

Abstract. An unphysical checkerboard pattern is identified in E3SMv2 and E3SM-MMF that is detectable across a wide range of timescales, from instantaneous snapshots to multi-year averages. A detection method is developed to quantify characteristics of the checkerboard signal by cataloguing all possible configurations of the eight adjacent neighbors for each cell on the model's cubed sphere grid using daily mean data. The checkerboard pattern is only found in cloud-related quantities, such as precipitation and liquid water path. Instances of pure and partial checkerboard are found to occur more often in E3SMv2 and E3SM-MMF when compared to satellite data regridded to the model grid. Continuous periods of partial checkerboard state are found to be more persistent in both models compared to satellite data, with E3SM-MMF exhibiting more persistence than E3SMv2. The checkerboard signal in E3SMv2 is found to be a direct consequence of the recently added deep convective trigger condition based on dynamically generated CAPE (DCAPE). In E3SM-MMF the checkerboard signal is found to be associated with the “trapping” of cloud-scale fluctuations within the embedded cloud-resolving model. Solutions to remedy this issue are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The COMBLE campaign: a study of marine boundary-layer clouds in Arctic cold-air outbreaks

One of the most intense air mass transformations on Earth happens when cold air flows from frozen surfaces to much warmer open water in cold-air outbreaks (CAOs), a process captured beautifully in satellite imagery. Despite the ubiquity of the CAO cloud regime over high-latitude oceans, we have a rather poor understanding of its properties, its role in energy and water cycles, and its treatment in weather and climate models. The Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) was conducted to better understand this regime and its representation in models. COMBLE aimed to examine the relations between surface fluxes, boundary-layer structure, aerosol, cloud and precipitation properties, and mesoscale circulations in marine CAOs. Processes affecting these properties largely fall in a range of scales where boundary-layer processes, convection, and precipitation are tightly coupled, which makes accurate representation of the CAO cloud regime in numerical weather prediction and global climate models most challenging. COMBLE deployed an Atmospheric Radiation Measurement Mobile Facility at a coastal site in northern Scandinavia (69°N), with additional instruments on Bear Island (75°N), from December 2019 to May 2020. CAO conditions were experienced 19% (21%) of the time at the main site (on Bear Island). A comprehensive suite of continuous in situ and remote sensing observations of atmospheric conditions, clouds, precipitation, and aerosol were collected. Because of the clouds’ well-defined origin, their shallow depth, and the broad range of observed temperature and aerosol concentrations, the COMBLE dataset provides a powerful modeling test bed for improving the representation of mixed-phase cloud processes in large-eddy simulations and large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Large‐Eddy Simulations of a Convection Cloud Chamber: Sensitivity to Bin Microphysics and Advection

Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.

54 ENVIRONMENTAL SCIENCES↗