Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large Eddy simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical analysis of particle dispersion and deposition in coal combustion using large-eddy simulation

High-fidelity modeling provides a useful approach to investigate the particle dispersion and deposition mechanism in pulverized coal combustion. To be able to analyze these detailed mechanisms, this work couples detailed multiphysics models for poly-dispersed turbulent, particle-laden flow, particle and gas phase reaction chemistry, convective, conductive and radiative heat transfer, ash formation and deposition with a Large-Eddy Simulation (LES) approach to numerically simulate coal combustion in a downfired self-sustained oxy-fuel combustor (OFC). It is necessary to explicitly capture all but the highest frequency dynamics of the turbulence and its coupling with each of the other physical phenomena using LES. Effects of subgrid-scale unresolved turbulence on the particle motions are also analyzed by a novel Stokes number analysis. Due to the inherent relation ship between coal combustion and ash deposition, this study integrates the improved ash deposition model into the numerical simulation. Overall simulation results are compared with experimentally measured data from the OFC. The simulation and measured data for the averaged gas temperature and deposition rates agree with 5% and 28%. This study shows that high-fidelity LES coupled with other detailed multiphysics models running on a exascale computing facility can provide a good representation of complex coal combustion and deposition in a laboratory-scale furnace.

multiphysics model↗
DOE OSTI2023

A baseline for ensemble-based, time-resolved inflow reconstruction for a single turbine using large-eddy simulations and latent diffusion models

We are interested in reconstructing winds flowing through a turbine on a second-by-second basis over a 10 min window. Previously, we developed a machine learning algorithm that takes in a snapshot of wind speed measurements and generates ensembles of three-dimensional wind field estimates. Here, we use these estimates as initial conditions in large-eddy simulations and reconstruct atmospheric and turbine response dynamic quantities in a synthetic field campaign. In doing so, we establish a baseline for model validation that future time-aware data assimilation techniques will be compared to. In turbine-free case studies, ground truth wind speeds consistently fall within our estimated wind speed distribution for the first 100 s after the simulation start. In simulations with turbines, the wind estimates show a small bias of 0.10 m s 1 and good correlation of 0.80 during the first 100 s. During this window, our estimates of the Blade 1 bending moment and generator power typically span the ground truth, with the estimate of the former performing better overall. In summary, this approach shows promise as a stand-alone technique for reconstructing real-world inflow and turbine dynamics in 1-2 min windows and as a foundation for future time-aware data assimilation techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

DPM: A deep learning PDE augmentation method with application to large-eddy simulation

A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Investigating the Origins of Cyclic Variability in Internal Combustion Engines Using Wall-Resolved Large Eddy Simulations

Modern internal combustion engines (ICE) operate at the ragged edge of stable operation characterized by high cycle-to-cycle variations (CCV). A key scientific challenge for ICE is the understanding, modeling, and control of CCV in engine performance, which can contribute to partial burns, misfire, and knock. The objective of this study is to use high-fidelity numerical simulations to improve the understanding of the causes of CCV. Nek5000, a leading high-order spectral element, open source code, is used to simulate the turbulent flow in the engine combustion chamber. Multicycle, wall-resolved large-eddy simulations (LESs) are performed for the General Motors (GM), Transparent Combustion Chamber (TCC-III) optical engine under motored operating conditions. The mean and root-mean-square (rms) of the in-cylinder flow fields at various piston positions are validated using particle image velocimetry (PIV) measurements during the intake and compression strokes. The large-scale flow structures, including the swirl and tumble flow patterns, are analyzed in detail and the causes for cyclic variabilities in these flow features are explained. The energy distribution across the different scales of the flow are quantified using one-dimensional (1D) energy spectra, and the effect of the tumble breakdown process on the energy distribution is examined. Finally, the insights from this study can help us develop improved engine designs with reduced cyclic variabilities in the in-cylinder flow leading to enhanced engine performance.

33 ADVANCED PROPULSION SYSTEMS↗
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 OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Flower‐Type Organized Trade‐Wind Cumulus: A Multi‐Day Lagrangian Large Eddy Simulation Intercomparison Study

Shallow cumulus cloud fields in subtropical marine trade wind environments, particularly over the tropical Atlantic Ocean, show distinct organizational patterns. Among these, Flower‐type clouds are characterized by expansive stratiform cloud patches surrounded by regions of scattered convection. The objectives of this study were (a) to construct a case study of a time period during the EUREC 4 A/ATOMIC field campaign when Flower‐type organization was observed, (b) to evaluate the fidelity of a multi‐model ensemble of large eddy simulations of that case, and (c) to analyze the interaction between cloud and precipitation processes and mesoscale organization in the simulations. The simulations follow a quasi‐Lagrangian trajectory, allowing mesoscale features to develop over time in a domain that follows the boundary‐layer airmass. The results show a broad agreement in simulated thermodynamic properties across different LES codes, with Flower‐type cloud patches appearing within hours of each other. The consensus among models is consistent with observations made during the EUREC 4 A/ATOMIC field campaign on the specific day of interest. The cloud structure reveals three distinct peaks in the joint probability densities of cloud base and cloud top height, with the dominant peak at any given time influenced by the stage of cloud organization. The simulated cloud system evolution reveals consistent occurrence of maxima in liquid water path and rain rate before Flower reaches its maximum length scale. Targeted sensitivity tests reveal a weak relationship between Cloud Droplet Number concentration and the extent/degree/type of organization.

EUREC4A↗
DOE OSTI2023

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗
DOE OSTI2025

Acoustic Tomography of the Atmosphere: A Large-Eddy Simulation Sensitivity Study

Accurate measurement of atmospheric turbulent fluctuations is critical for understanding environmental dynamics and improving models in applications such as wind energy. Advanced remote sensing technologies are essential for capturing instantaneous velocity and temperature fluctuations. Acoustic tomography (AT) offers a promising approach that utilizes sound travel times between an array of transducers to reconstruct turbulence fields. This study presents a systematic evaluation of the time-dependent stochastic inversion (TDSI) algorithm for AT using synthetic travel-time measurements derived from large-eddy simulation (LES) fields under both neutral and convective atmospheric boundary-layer conditions. Unlike prior work that relied on field observations or idealized fields, the LES framework provides a ground-truth atmospheric state, enabling quantitative assessment of TDSI retrieval reliability, sensitivity to travel-time measurement noise, and dependence on covariance model parameters and temporal data integration. A detailed sensitivity analysis was conducted to determine the best-fit model parameters, identify the tolerance thresholds for parameter mismatch, and establish a maximum spatial resolution. The TDSI algorithm successfully reconstructed large-scale velocity and temperature fluctuations with root mean square errors ( RMSE s) below 0.35 m/s and 0.12 K, respectively. Spectral analysis established a maximum spatial resolution of approximately 1.4 m, and reconstructions remained robust for travel-time measurement uncertainties up to 0.002 s. These findings provide critical insights into the operational limits of TDSI and inform future applications of AT for atmospheric turbulence characterization and system design.

17 WIND ENERGY↗
DOE OSTI2025

Factors Controlling Precipitation Onset and Maintenance: Inferences From Large Eddy Simulations of Two Afternoon Deep Convective Regimes Over Amazon

To investigate the environmental factors controlling the onset and maintenance of afternoon precipitation over tropical rainforests, two contrasting cases are created for large eddy simulations (LES) using GoAmazon observations: one features a single-pulse rain dissipating quickly in earlier afternoon, while the other shows double pulses lasting until evening. For these two specific cases under consideration, LES confirms that early-morning relative humidity dominates afternoon rain patterns. The attributed impacts are distinct: moisture determines the rain onset timing while temperature affects the peak rain intensity. Single-pulse day observes one round of strong precipitation and cold pools, which further suppresses convection. On double-pulse day, the first precipitation peak results from an intermediate development of congestus clouds, whose detrainment leads to a gradual moistening of the lower-to-middle troposphere. This behavior of convection favors a second pulse of stronger precipitation with more convective organization, whose development decouples from surface fluxes and sustains until evening.

Amazon↗
DOE OSTI2020

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗
DOE OSTI2024

ECAR-7932 Rev 0 Large Eddy Simulation of MARVEL Reactor Core Subchannel to Evaluate Model Uncertainty of Reynolds-Averaged Navier-Stokes Equation Based Computational Fluid Dynamics Analysis

In the previous work (ECAR-7210), the peak cladding temperature of the MARVEL microreactor has been evaluated by steady-state Reynolds-Averaged Navier-Stokes (RANS) based computational fluid dynamics (CFD) simulations. Although numerical uncertainties of RANS-based CFD simulations has been assessed in ECAR-7210, the model uncertainty of RANS turbulence models must be investigated to resolve the issues related to inaccurate prediction of turbulent heat flux and flow pulsation in a tight lattice rod bundle using the steady-state RANS simulations. Consequently, this ECAR conducted a high-fidelity CFD analysis utilizing Large Eddy Simulation (LES) to generate reference data and investigated the model uncertainty of RANS-based CFD simulations.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

The FastEddy® Resident-GPU Accelerated Large-Eddy Simulation Framework: Moist Dynamics Extension, Validation and Sensitivities of Modeling Non-Precipitating Shallow Cumulus Clouds

Herein we describe the moist dynamics formulation implemented within the graphics processing unit-resident large-eddy simulation FastEddy® model, which includes a simple saturation adjustment scheme for condensation and evaporation processes. Two LES model intercomparison exercises for non-precipitating shallow cumulus clouds are simulated in order to validate this model extension, including a static forcing and a time-dependent forcing case. Overall, we find our dynamical, thermodynamical and microphysical quantities, along with turbulence variability and fluxes, to be commensurate with the corresponding model intercomparison results. In addition, sensitivities to specific model settings are investigated. Among these settings, it is shown that boundary layer and cloud layer structure and characteristics are sensitive to use of higher-order advection schemes impacting the vertical distribution of cloud content and associated turbulence statistics. Increasing the timescale of the saturation scheme leads to enhanced liquid water presence and decreases vertical velocity variance within the cloud deck. In some cases, these sensitivities agree with the model-to-model variability reported in the intercomparison exercises, highlighting the important role of specific model implementation choices in the context of shallow cumulus convection simulations. These analyses and findings also provide the basis for future extensions and applications of FastEddy® for modeling moist convection and precipitation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Measurement-driven large-eddy simulations of a diurnal cycle during a wake-steering field campaign

Abstract. High-fidelity flow modeling with data assimilation enables accurate representation of the wind farm operating environment under realistic, nonstationary atmospheric conditions. Capturing the temporal evolution of the turbulent atmospheric boundary layer is critical to understanding the behavior of wind turbines under operating conditions with simultaneously varying inflow and control inputs. This paper has three parts: the identification of a case study during a field evaluation of wake steering; the development of a tailored mesoscale-to-microscale coupling strategy that resolved local flow conditions within a large-eddy simulation (LES), using observations that did not completely capture the wind and temperature fields throughout the simulation domain; and the application of this coupling strategy to validate high-fidelity aeroelastic predictions of turbine performance and wake interactions with and without wake steering. The case study spans 4.5 h after midnight local time, during which wake steering was toggled on and off five times, achieving yaw offset angles ranging from 0 to 17°. To resolve nonstationary nighttime conditions that exhibited shear instabilities, the turbulence field was evolved starting from the diurnal cycle of the previous day. These background conditions were then used to drive wind farm simulations with two different models: an LES with actuator disk turbines and a steady-state engineering wake model. Subsequent analysis identified two representative periods during which the up- and downstream turbines were most nearly aligned with the mean wind direction and had observed yaw offsets of 0 and 15°. Both periods corresponded to partial waking on the downstream turbine, which had errors in the LES-predicted power of 4 % and 6 %, with and without wake steering. The LES was also able to capture conditions during which an upstream turbine wake induced a speedup at a downstream turbine and increased power production by up to 13 %.

17 WIND ENERGY↗
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↗