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.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2025

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↗

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 OSTI2026

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↗

DOE OSTI2024

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)↗

DOE OSTI2024

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↗

DOE OSTI2026

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↗

DOE OSTI2025

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↗

DOE OSTI2025

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↗

DOE OSTI2026

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↗

DOE OSTI2026

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↗

NASA NTRS1975

Computation of turbulent flows

Substantial advances made over the past decade in the prediction of turbulent flows are discussed. There was extensive work in the development of turbulence models, particularly for use in boundary layer calculations. Basic aspects of several important methods based on partial differential equations for the mean velocity field and turbulence quantities, including the relationship between the methods and suggestions for future development were reviewed. Work on three-dimensional time-dependent large eddy simulations is discussed. The emphasis is on the hydrodynamics of incompressible flows, but sources for consideration of heat transfer and compressibility are mentioned.

Reynolds, W. C.↗

NASA NTRS1978

Levels of turbulence prediction

A classification according to the level of detail of description the computational method provides is probably most useful. Flow calculations can be classified into five categories: (1) correlations; (2) zonal method; (3) time-averaged equations; (4) large-eddy simulation; and (5) Navier-Stokes solution. There are methods that fall into more than one category, and there are sub-divisions of each category. A discussion of the advantages and disadvantages of each of these five categories is given.

Ferziger, J. H.↗

NASA NTRS1980

Developments in the computation of turbulent boundary layers

Two methods of turbulence computation are discussed in terms of their basic simularities. It is shown that the two methods are interrelated and that each can gain from advances in the other. The degree of success of a pair of increasingly complex Reynolds stress models to broaden their range of applicability is examined through comparison with experimental data for a variety of flow conditions. An example of a large eddy simulation is presented, compared with experimental results, and used to evaluate the models for pressure rate of strain correlation and dissipation in the Reynolds averaged equations.

Rubesin, M. W.↗

NASA NTRS1981

State of the art in subgrid scale modeling

The present state of the art in subgrid scale modeling for large eddy simulation is reviewed. It is found that the eddy viscosity models are sufficient for simulating homogeneous flows but are probably insufficient for inhomogeneous shear flows. One- and two-equation subgrid scale models probably will not be of great value in subgrid scale modeling, and full Reynolds stress models may be required. In the latter, the convection and diffusion terms are expected to be of greater importance than they are in time-average modeling, but the pressure redistribution terms may be of less importance. A scale similarity model may provide an alternative.

Ferziger, J. H.↗

NASA NTRS1981

Trends and pacing items in computational aerodynamics

A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.

Chapman, D. R.↗

NASA NTRS1983

On the structure of wall-bounded turbulent flows

The variable-interval time-averaging (VITA) technique developed by Blackwelder and Kaplan is applied to data obtained from large-eddy simulation of turbulent channel flow in an investigation of the organized structures associated with the bursting phenomenon in the near-wall region. Conditionally averaged velocities, shear stress, pressure, and vorticity are discussed in conjunction with the bursting phenomenon detected by the VITA technique. The conditionally averaged pressure reveals that the ejection process is associated with a localized adverse pressure gradient. In the plane perpendicular to the flow direction, the conditionally averaged vorticity field indicates that a pair of counterrotating streamwise vorticity is being lifted through the ejection process.

Kim, J.↗

NASA NTRS1983

On the structure of wall-bounded turbulent flows

Kim, J.↗

NASA NTRS1984

Turbulence structures associated with the bursting event

Turbulence structures in a wall bounded shear layer during the bursting event detected by a conditional sampling technique are investigated using data obtained from large eddy simulation of turbulent channel flow. Streamlines are constructed from the ensemble-averaged velocity field to illustrate the flow patterns associated with the bursting event. The spatting motions during the sweep event and the existence of a pair of counterrotating streamwise vortices during the ejection process are exhibited.

Kim, J.↗