Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large eddy simulation”

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 1,603 records · Page 89

Impacts of solar-absorbing aerosol layers on the transition of stratocumulus to trade cumulus clouds

Abstract. The effects of an initially overlying layer of solar-absorbing aerosol on the transition of stratocumulus to trade cumulus clouds are examined using large-eddy simulations. For lightly drizzling cloud the transition is generally hastened, resulting mainly from increased cloud droplet number concentration (Nc) induced by entrained aerosol. The increased Nc slows sedimentation of cloud droplets and shortens their relaxation time for diffusional growth, both of which accelerate entrainment of overlying air and thereby stratocumulus breakup. However, the decrease in albedo from cloud breakup is more than offset by redistributing cloud water over a greater number of droplets, such that the diurnal-average shortwave forcing at the top of the atmosphere is negative. The negative radiative forcing is enhanced by sizable longwave contributions, which result from the greater cloud breakup and a reduced boundary layer height associated with aerosol heating. A perturbation of moisture instead of aerosol aloft leads to a greater liquid water path and a more gradual transition. Adding absorbing aerosol to that atmosphere results in substantial reductions in liquid water path (LWP) and cloud cover that lead to positive shortwave and negative longwave forcings on average canceling each other. Only for heavily drizzling clouds is the breakup delayed, as inhibition of precipitation overcomes cloud water loss from enhanced entrainment. Considering these simulations as an imperfect proxy for biomass burning plumes influencing Namibian stratocumulus, we expect regional indirect plus semi-direct forcings to be substantially negative to negligible at the top of the atmosphere, with its magnitude sensitive to background and perturbation properties.

54 ENVIRONMENTAL SCIENCES↗

The relative impact of cloud condensation nuclei and ice nucleating particle concentrations on phase partitioning in Arctic mixed-phase stratocumulus clouds

This study investigates the interactions between cloud dynamics and aerosols in idealized large-eddy simulations (LES) of Arctic mixed-phase stratocumulus clouds (AMPS) observed at Oliktok Point, Alaska, in April 2015. This case was chosen because it allows the cloud to form in response to radiative cooling starting from a cloud-free state, rather than requiring the cloud ice an dliquid to adjust to an initial cloudy state. Sensitivity studies are used to identify whether there are buffering feedbacks that limit the impact of aerosol perturbations. The results of this study indicate that perturbations in ice nucleating particles (INPs) dominate over cloud condensation nuclei(CCN) perturbations; i.e., an equivalent fractional decrease in CCN and INPs results in an increase in the cloud-top long wave cooling rate, even though the droplet effective radius increases and the cloud emissivity decreases.The dominant effect of ice in the simulated mixed-phase cloud is a thinning rather than a glaciation, causing the mixed-phase clouds to radiate as a grey body and the radiative properties of the cloud to be more sensitive to aerosol perturbations. It is demonstrated that allowing prognostic CCN and INPs causes a layering of the aerosols, with increased concentrations of CCN above cloud top and increased concentrations of INPs at the base of the cloud-driven mixed layer. This layering contributes to the maintenance of the cloud liquid, which drives the dynamics of the cloud system.

54 ENVIRONMENTAL SCIENCES↗

Environmental sensitivities of shallow-cumulus dilution – Part 2: Vertical wind profile

This second part of a numerical study on shallow-cumulus dilution focuses on the sensitivity of cloud dilution to changes in the vertical wind profile. Insights are obtained through large-eddy simulations of maritime and continental cloud fields. In these simulations, the speed of the initially uniform geostrophic wind and the strength of geostrophic vertical wind shear in the cloud and subcloud layer are varied. Increases in the cloud-layer vertical wind shear (up to 9 ms -1 km -1 ) lead to 40%–50% larger cloud-core dilution rates compared to their respective unsheared counterparts. When the background wind speed, on the other hand, is enhanced by up to 10 ms -1 and subcloud-layer vertical wind shear develops or is initially prescribed, the dilution rate decreases by up to 25%. The sensitivities of the dilution rate are linked to the updraft strength and the properties of the entrained air. Increases in the wind speed or vertical wind shear result in lower vertical velocities across all sets of experiments with stronger reductions in the cloud-layer wind shear simulation (27%–47%). Weaker updrafts are exposed to mixing with the drier surrounding air for a longer time period, allowing more entrainment to occur (i.e., the “core-exposure effect”). However, reduced vertical velocities, in concert with increased cloud-layer turbulence, also assist in widening the humid shell surrounding the cloud cores, leading to entrainment of more humid air (i.e., the “core–shell dilution effect”). In the experiments with cloud-layer vertical wind shear, the core-exposure effect dominates and the cloud-core dilution increases with increasing shear. Conversely, when the wind speed is increased and subcloud-layer vertical wind shear develops or is imposed, the core–shell dilution effect dominates to induce a buffering effect. The sensitivities are generally stronger in the maritime simulations, where weaker sensible heat fluxes lead to narrower, more tilted, and, therefore, more suppressed cumuli when cloud-layer shear is imposed. Moreover, in the experiments with subcloud wind shear, the weaker baseline turbulence in the maritime case allows for a larger turbulence enhancement, resulting in a widening of the transition zones between the cores and their environment, leading to the entrainment of more humid air.

54 ENVIRONMENTAL SCIENCES↗

Cloud water adjustments to aerosol perturbations are buffered by solar heating in non-precipitating marine stratocumuli

Abstract. Marine low-level clouds are key to the Earth's energy budget due to their expansive coverage over global oceans and their high reflectance of incoming solar radiation. Their responses to anthropogenic aerosol perturbations remain the largest source of uncertainty in estimating the anthropogenic radiative forcing of climate. A major challenge is the quantification of the cloud water response to aerosol perturbations. In particular, the presence of feedbacks through microphysical, dynamical, and thermodynamical pathways at various spatial and temporal scales could augment or weaken the response. Central to this problem is the temporal evolution in cloud adjustment, governed by entangled feedback mechanisms. We apply an innovative conditional Monte Carlo subsampling approach to a large ensemble of diurnal large-eddy simulation of non-precipitating marine stratocumulus to study the role of solar heating in governing the evolution in the relationship between droplet number and cloud water. We find a persistent negative trend in this relationship at night, confirming that the role of microphysically enhanced cloud-top entrainment. After sunrise, the evolution in this relationship appears buffered and converges to ∼-0.2 in the late afternoon. This buffering effect is attributed to a strong dependence of cloud-layer shortwave absorption on cloud liquid water path. These diurnal cycle characteristics further demonstrate a tight connection between cloud brightening potential and the relationship between cloud water and droplet number at sunrise, which has implications for the impact of the timing of advertent aerosol perturbations.

Zhang, Jianhao (ORCID:0000000169882935)↗

Above-cloud concentrations of cloud condensation nuclei help to sustain some Arctic low-level clouds

Abstract. Previous studies have found that low-level Arctic clouds often persist for long periods even in the face of very low surface cloud condensation nuclei (CCN) concentrations. Here, we investigate whether these conditions could occur due to continuous entrainment of aerosol particles from the free troposphere (FT). We use an idealized large eddy simulation (LES) modeling framework, where aerosol concentrations are low in the boundary layer (BL) but increased up to 50× in the free troposphere. We find that the tests with higher tropospheric aerosol concentrations simulated clouds, which persisted for longer and maintained higher liquid water paths (LWPs). This is due to direct entrainment of the tropospheric aerosol into the cloud layer, which results in a precipitation suppression from the increase in cloud droplet number and in stronger cloud-top radiative cooling, which causes stronger circulations maintaining the cloud in the absence of surface forcing. Together, these two responses result in a more well-mixed boundary layer with a top that remains in contact with the tropospheric aerosol reservoir and can maintain entrainment of those aerosol particles. The surface aerosol concentrations, however, remained low in all simulations. The free-tropospheric aerosol concentration necessary to maintain the clouds is consistent with concentrations that are frequently seen in observations.

Environmental Sciences & Ecology↗

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗

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↗

A robust aerosol impact on clouds along the subtropical to tropical transition

Marine clouds undergo a transition from subtropical stratocumulus (Sc) to shallow cumulus (Cu) and eventually to deep convective (DC) systems as air masses advect from the subtropics towards the deep tropics. How aerosols modulate this Lagrangian cloud evolution remains largely uncertain. Here we use both 5 years of satellite observations mapped along 8 d Lagrangian trajectories and complementary large-eddy simulations from 9 initiation locations across the Northeast Pacific, Southeast Pacific, and Southeast Atlantic. This Lagrangian framework allows us to quantify the aerosol effect and its co-variability with meteorological conditions on cloud microphysics, macrophysics, and top-of-atmosphere radiation through the full Sc-Cu-DC transition. We show that increasing aerosol concentrations leads to deeper and more reflective clouds throughout this cloud transition. Examining the thermodynamic evolution along the trajectory indicates a well-known trend: enhanced moistening near the boundary-layer top and lower free troposphere under polluted conditions, suggesting that part of the co-variability between aerosol and meteorological conditions may be internally driven. The agreement between model simulations and satellite data alongside the multi-basin coherence of the results indicates that aerosols systematically amplify cloud depth and reflectivity during the subtropical–to–tropical cloud transition.

54 ENVIRONMENTAL SCIENCES↗

Improving the estimate of higher-order moments from lidar observations near the top of the convective boundary layer

Abstract. Ground-based lidar data have proven extremely useful for profiling the convective boundary layer (CBL). Many groups have derived higher-order moments (e.g., variance, skewness, fluxes) from high-temporal-resolution lidar data using an autocovariance approach. However, these analyses are highly uncertain near the CBL top when the depth of the CBL (zi) is changing during the analysis period. This is because the autocovariance approach is usually applied to constant height levels and the character of the eddies is changing on either side of the changing CBL top. Here, a new approach is presented wherein the autocovariance analysis is performed on a normalized height grid, with a temporally smoothed zi. Output from a large eddy simulation model demonstrates that deriving higher-order moments from time series on a normalized height grid has better agreement with the slab-averaged quantities than the moments derived from the original height grid.

Rosenberger, Tessa E. (ORCID:0000000333205873)↗

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗

A new subgrid-scale representation of hydrometeor fields using a multivariate PDF

Abstract. The subgrid-scale representation of hydrometeor fields is important for calculating microphysical process rates. In order to represent subgrid-scale variability, the Cloud Layers Unified By Binormals (CLUBB) parameterization uses a multivariate probability density function (PDF). In addition to vertical velocity, temperature, and moisture fields, the PDF includes hydrometeor fields. Previously, hydrometeor fields were assumed to follow a multivariate single lognormal distribution. Now, in order to better represent the distribution of hydrometeors, two new multivariate PDFs are formulated and introduced.The new PDFs represent hydrometeors using either a delta-lognormal or a delta-double-lognormal shape. The two new PDF distributions, plus the previous single lognormal shape, are compared to histograms of data taken from large-eddy simulations (LESs) of a precipitating cumulus case, a drizzling stratocumulus case, and a deep convective case. Finally, the warm microphysical process rates produced by the different hydrometeor PDFs are compared to the same process rates produced by the LES.

97 MATHEMATICS AND COMPUTING↗

Revisiting Kelvin Helmholtz Instabilities and von Kármán Vortices in Canopy Turbulence

Studying turbulence in vegetation canopies is vital in the context of a number of micrometeorological and hydrological applications. While recent focus has shifted more towards exploring different kinds of canopy heterogeneities, there are still gaps in the existing knowledge on the multiple types of dynamics involved in the case of horizontally homogeneous canopies. For example, experimental studies have indicated that turbulence in the canopy sublayer (CSL) can be divided into three regimes. In the deep-zone, the flow-field is dominated by von Kármán vortex streets and interrupted by strong sweep events. The second zone near the canopy top is dominated by attached eddies and Kelvin-Helmholtz waves associated with the velocity inflection point in the mean longitudinal velocity profile. Above the canopy, the flow resembles classical boundary layer flow. In this study, these different kinds of dynamics are studied together by means of a large eddy simulation (LES). The main theme of this work is to address the question whether the parametrization of the canopy by a distributed drag force in numerical simulations instead of placing real solid obstacles is consistent with the three layer conceptual model. Unique techniques such as measures from information theory and coupled oscillator analysis are used to extract the coherent structures associated with the two motions. It can be said that a better understanding of the rich dynamics associated with the simplest case of canopy turbulence can lead to more efficient simulations and more importantly improve the interpretation of more complex scenarios.

canopy turbulence↗

Spectral proper orthogonal decomposition of active wake mixing dynamics in a stable atmospheric boundary layer

Recent advancements in the use of active wake mixing (AWM) to reduce wake effects on downstream turbines open new avenues for increasing power generation in wind farms. However, a better understanding of the fluid dynamics underlying AWM is still needed to make wake mixing a reliable strategy for wind farm flow control. In this work, a spectral proper orthogonal decomposition (SPOD) is used to analyze the dynamics of coherent flow structures that are induced in the wake through blade pitch actuation. The data are generated using the ExaWind software suite to perform large eddy simulations of an NREL 2.8 MW turbine operating in a stable atmospheric boundary layer. SPOD tracks the modal behavior of flow structures from their generation in the turbine induction field through their growth in the near-wake region and to their subsequent evolution and energy transfers in the far wake. SPOD is shown to be a useful tool in the context of AWM because it translates the wavenumber and frequency inputs to the turbine controller to structures in the wake. A decomposition of the radial shear stress flux in the wake is also developed using SPOD to measure the contribution of coherent flow structures to mean flow turbulent entrainment and wake recovery. The effectiveness of AWM is connected to its ability to excite inherent structures in the wake of the turbine that arise using baseline controls. The effects of AWM on blade loading are also analyzed by connecting the axial force along the blade to the SPOD analysis of the turbine induction field. Lastly, the performance of different AWM strategies is demonstrated in a two-turbine array.

17 WIND ENERGY↗

Observations of wind farm wake recovery at an operating wind farm

Abstract. The interplay of momentum surrounding wind farms significantly influences wake recovery, affecting the speed at which wakes return to their freestream velocities. Under stable atmospheric conditions, wind farm wakes can extend over considerable distances, leading to sustained vertical momentum flux downstream, with variations observed throughout the diurnal cycle. Particularly in regions such as the US Great Plains, stable conditions can induce low-level jets (LLJs), impacting wind farm performance and power output. This study examines the implications of wake recovery using long-term observations of vertical momentum flux profiles across diverse atmospheric conditions. In these observations, several key findings were observed, such as (a) LLJ heights being altered downstream of a wind farm, especially when the LLJs are below 250 m above ground level; (b) a notable impact of LLJ height on wake recovery being observed using momentum flux profiles at upwind and downwind locations, wherein LLJs between 250 and 500 m above ground level resulted in larger momentum transfer within the wake (i.e., smaller velocity deficit) compared to LLJs below 250 m above ground level; (c) the largest momentum flux variability being observed during stable atmospheric conditions, with non-negligible variability observed during neutral and unstable atmospheric conditions; (d) detection of wake effects almost always being observed throughout the atmospheric boundary layer height; and finally (e) enhancement of wake recovery being observed in the presence of propagating gravity waves. These insights deepen our understanding of the intricate dynamics governing wake recovery in wind farms, advancing efforts to model and predict their behavior across varying atmospheric contexts. In addition, the performance of large-eddy-simulation-based semi-empirical internal boundary layer height model estimates incorporating real-world atmospheric and turbine inputs was evaluated using observations during LLJ conditions.

17 WIND ENERGY↗

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗

Introduction to and comparison of deep learning and optimization approaches to analytical wake modeling of a tilted wind turbine

This paper introduces innovative optimization and deep learning techniques to enhance the prediction of complex wake dynamics in the downstream wind velocity of tilted wind turbines. Traditional methods for calibrating the Bastankhah wake model often lead to increased errors in wind velocity distribution due to overfitting of the local wake characteristics. To address this issue, we propose an additional global optimization step to reduce errors in wind velocity predictions with respect to various wake parameters. Despite this improvement, the Bastankhah model's axisymmetric Gaussian wake shape limits its accuracy for complex wake structures. Therefore, we also propose a deep learning approach, which demonstrates promising results by accurately modeling complex wake shapes across a broader range of tilt angles with minimal computational cost. The deep learning approach achieves near-identical predictions to high-fidelity large-eddy simulations, representing a promising advancement in wake modeling.

17 WIND ENERGY↗