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 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗

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 OSTI2026

Generating a stratocumulus-like cloud top in a convection-cloud chamber

Stratocumulus-topped boundary layers play a crucial role in influencing daily weather and earth energy balance. Entrainment at the stratocumulus cloud top affects the cloud’s lifetime, precipitation, and radiative properties, but our understanding remains limited due to the lack of resolution in both field observations and numerical simulations. A recently proposed convection-cloud chamber with detailed control of sidewall temperatures can provide a unique opportunity to explore this mechanism in a laboratory setting. In this work, we use numerical simulations to demonstrate that this design can produce a cloud top that mimics the entrainment interfacial layer in a stratocumulus cloud. Our results show that a steady-state cloud can be formed by cooling the lower portions of the sidewalls and warming the bottom surface, while a temperature inversion at the cloud top can be generated by keeping the upper sidewalls and top surface warmer than the bottom. The turbulent kinetic energy profile and budget are similar to those found in a convective boundary layer, and inhomogeneous mixing near the cloud top can be observed. These findings significantly enhance the scientific value of constructing the tall convection-cloud chamber.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Inverse Mapping of the Collision Kernel and Wall Flux Scaling in a Tall Convection‐Cloud Chamber Using Local Sensors and Knowledge‐Informed Deep Learning

Droplet collision–coalescence is a crucial process in cloud physics, but accurately representing this process under different dynamical conditions remains challenging. A proposed future convective‐cloud chamber aims to investigate this key process, but the method for observing it remains unclear, even though it is theoretically established that collision‐coalescence will occur. This study serves as a proof‐of‐concept demonstration of how knowledge‐informed deep learning, combined with measurement data from local sensors in the chamber, can be used to estimate the collision kernels, which determine how the droplet size distribution evolves during collision‐coalescence. In addition to estimating the collision kernel, we also address wall fluxes, another uncertain but important process that acts as a source of heat and moisture in the chamber. Ensemble runs of large‐eddy simulations are conducted by scaling the wall fluxes and the collision kernel, while the measured flow and cloud properties are used as inputs for a neural network. Results indicate that this approach successfully maps the scaling of wall fluxes and the collision kernel with biases of approximately 1% or less relative to the range of the target data. This proof‐of‐concept lays the groundwork for future applications; when the real measurements are available, real sensor data combined with the trained model presented in this work will enable estimation of the actual wall fluxes and collision kernel.

cloud chamber↗

DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗

DOE OSTI2025

Topography-Induced TKE Budget Behavior Over an Amazon Forest

Data from three different heights (35, 50 and 81 m) of one of the ATTO Project towers in the Amazon forest were used to calculate the TKE (turbulence kinetic energy) budget and some other statistics within the RSL (roughness sublayer). The statistical analyses were carried out for unstable and stable cases. The vertical transport and the vertical advection terms do not explain the imbalances found in the TKE budget, highlighting the impor- tance of horizontal terms in complex terrain. Indeed, patterns independent of the time of day were found for vertical transport, horizontal and vertical velocity skewness, the mean vertical component velocity, and momentum and vertical TKE flux. Here, the highest values of mean vertical velocity were observed from the direction where there is a valley near the tower, which may indicate a topographical effect. Monin-Obukhov Similarity Theory (MOST) predictions for the mean velocity gradient, strictly applicable only in the inertial sublayer (ISL), gave reasonable predictions at 35 and 50 m, but dimensionless dissipation rates displayed large scatter. At 81 m, there were clear deviations from MOST for the two functions, disclosing the effect of topography. In addition, a neutral LES (large-eddy simu- lation) study showed good agreement with tower data for the same general wind direction.

ATTO project↗

DOE OSTI2025

ERF: Energy Research and Forecasting Model

High performance computing (HPC) architectures have undergone rapid development in recent years. As a result, established software suites face an ever increasing challenge to remain performant on and portable across modern systems. Many of the widely adopted atmospheric modeling codes cannot fully (or in some cases, at all) leverage the acceleration provided by General-Purpose Graphics Processing Units, leaving users of those codes constrained to increasingly limited HPC resources. Energy Research and Forecasting (ERF) is a regional atmospheric modeling code that leverages the latest HPC architectures, whether composed of only Central Processing Units (CPUs) or incorporating GPUs. ERF contains many of the standard discretizations and basic features needed to model general atmospheric dynamics. The modular design of ERF provides a flexible platform for exploring different physics parameterizations and numerical strategies. ERF is built on a state-of-the-art, well-supported, software framework (AMReX) that provides a performance portable interface and ensures ERF's long-term sustainability on next generation computing systems. This paper details the numerical methodology of ERF, presents results for a series of verification/validation cases, and documents ERF's performance on current HPC systems. The roughly 5× speed up of ERF (using GPUs) over Weather Research and Forecasting (CPUs only) for a 3D squall line test case highlights the significance of leveraging GPU acceleration.

17 WIND ENERGY↗

DOE OSTI2025

Enhancing Turbulent Mixing and Microphysical Uniformity in a Tall Convection‐Cloud Chamber Through Idealized Heterogeneity of Boundaries

A large convection cloud chamber has been proposed for exploring aerosol–cloud–drizzle interactions under well‐controlled turbulent conditions. Recent theoretical and numerical studies suggest that a convection cloud chamber with two heated and two cooled sidewalls can significantly enhance the liquid water content and thus benefit drizzle initiation. However, a chamber with such a sidewall configuration develops stable stratification and extremely weak turbulence therein. In this study, we conduct large‐eddy simulations of a tall convection chamber with five different sidewall configurations consisting of alternating warm and cold patches. For each configuration, the total surface area of warm patches equals that of cold patches, resulting in the same expected cloud‐free supersaturation based on a flux budget model. Results show that changing the sidewall configuration, while keeping all other factors constant, can substantially enhance turbulent mixing and improve the uniformity of thermodynamic and cloud microphysical properties in the bulk region of the chamber. In addition, turbulence strength is positively correlated with liquid water content and negatively correlated with cloud droplet number concentration, consistent with theoretical predictions. Our results highlight the advantage of building a large cloud chamber using modular patches with individually controllable temperature and humidity to achieve well‐mixed conditions.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2002

Limits to Cloud Susceptibility

1-kilometer AVHRR observations of ship tracks in low-level clouds off the west coast of the U S. were used to determine limits for the degree to which clouds might be altered by increases in anthropogenic aerosols. Hundreds of tracks were analyzed to determine whether the changes in droplet radii, visible optical depths, and cloud top altitudes that result from the influx of particles from underlying ships were consistent with expectations based on simple models for the indirect effect of aerosols. The models predict substantial increases in sunlight reflected by polluted clouds due to the increases in droplet numbers and cloud liquid water that result from the elevated particle concentrations. Contrary to the model predictions, the analysis of ship tracks revealed a 15-20% reduction in liquid water for the polluted clouds. Studies performed with a large-eddy cloud simulation model suggested that the shortfall in cloud liquid water found in the satellite observations might be attributed to the restriction that the 1-kilometer pixels be completely covered by either polluted or unpolluted cloud. The simulation model revealed that a substantial fraction of the indirect effect is caused by a horizontal redistribution of cloud water in the polluted clouds. Cloud-free gaps in polluted clouds fill in with cloud water while the cloud-free gaps in the surrounding unpolluted clouds remain cloud-free. By limiting the analysis to only overcast pixels, the current study failed to account for the gap-filling predicted by the simulation model. This finding and an analysis of the spatial variability of marine stratus suggest new ways to analyze ship tracks to determine the limit to which particle pollution will alter the amount of sunlight reflected by clouds.

Coakley, James A., Jr.↗

NASA NTRS1985

Effect of rotation on isotropic turbulence - Computation and modelling

This paper uses numerical simulation to analyse the effects of uniform rotation on homogeneous turbulence. Both large-eddy and full simulations were made. The results indicate that the predominant effect of rotation is to decrease the rate of dissipation of the turbulence and increase the lengthscales, especially those along the axis of rotation. These effects are a consequence of the reduction, due to the generation of inertial waves, of the net energy transfer from large eddies to small ones. Experiments are also influenced by a more complicated interaction between the rotation and the wakes of the turbulence-generating grid which modifies the nominal initial conditions in the experiment. The latter effect is accounted for in simulations by modifying the initial conditions. Finally, a two-equation model is proposed that accounts for the effects of rotation and is able to reproduce the experimental decay of the turbulent kinetic energy.

Bardina, J.↗

DOE OSTI2025

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗

NASA NTRS1991

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗

NASA NTRS1990

Boundary layer transition in hypersonic flows

This paper summarizes some of the recent progress made at NASA Langley Research Center in the understanding, prediction and modeling of high speed boundary-layer transition. Linear and nonlinear theories together with large-eddy and direct numerical simulations have been used to understand various aspects of the transition problem while low disturbance 'quiet' tunnels provide means for validating the theoretical results.

Malik, M. R.↗

NASA NTRS2010

LES/RANS Simulation of a Supersonic Reacting Wall Jet

This work presents results from large-eddy / Reynolds-averaged Navier-Stokes (LES/RANS) simulations of the well-known Burrows-Kurkov supersonic reacting wall-jet experiment. Generally good agreement with experimental mole fraction, stagnation temperature, and Pitot pressure profiles is obtained for non-reactive mixing of the hydrogen jet with a non-vitiated air stream. A lifted flame, stabilized between 10 and 22 cm downstream of the hydrogen jet, is formed for hydrogen injected into a vitiated air stream. Flame stabilization occurs closer to the hydrogen injection location when a three-dimensional combustor geometry (with boundary layer development resolved on all walls) is considered. Volumetric expansion of the reactive shear layer is accompanied by the formation of large eddies which interact strongly with the reaction zone. Time averaged predictions of the reaction zone structure show an under-prediction of the peak water concentration and stagnation temperature, relative to experimental data and to results from a Reynolds-averaged Navier-Stokes calculation. If the experimental data can be considered as being accurate, this result indicates that the present LES/RANS method does not correctly capture the cascade of turbulence scales that should be resolvable on the present mesh. Instead, energy is concentrated in the very largest scales, which provide an over-mixing effect that excessively cools and strains the flame. Predictions improve with the use of a low-dissipation version of the baseline piecewise parabolic advection scheme, which captures the formation of smaller-scale structures superimposed on larger structures of the order of the shear-layer width.

Edwards, Jack R.↗

NASA NTRS2003

Aerosol-Cloud Interactions during Tropical Deep Convection: Evidence for the Importance of Free Tropospheric Aerosols

NASA's 2002 CRYSTAL-FACE field experiment focused on the formation and evolution of tropical cirrus cloud systems in southern Florida. Multiple aircraft extensively sampled cumulonimbus dynamical and microphysical properties, as well as characterizing ambient aerosol populations both inside and outside the full depth of the convective column. On July 18, unique measurements were taken when a powerful updraft was traversed directly by aircraft, providing a window into the primary source region of cumulonimbus anvil crystals. Observations of the updraft, entered at approximately l0 km altitude and -34 C, indicated more than 200 cloud particles per mL at vertical velocities exceeding 20 m/s and the presence of significant condensation nuclei and liquid water within the core. In this work, aerosol and cloud phase observations are integrated by simulating the updraft conditions using a large-eddy resolving model with 3 explicit multiphase microphysics, including treatment of size-resolved aerosol fields, aerosol activation and freezing, and evaporation of cloud particles back to the aerosol phase. Simulations were initialized with observed thermodynamic and aerosol size distributions profiles and convection was driven by surface fluxes assimilated from the ARPS forecast model. Model results are consistent with the conclusions that most crystals are homogeneously frozen droplets and that entrained free tropospheric aerosols may contribute a significant fraction of the crystals. Thus most anvil crystals appear to be formed aloft in updraft cores, well above cloud base. These conclusions are supported by observations of hydrometeor size distribution made while traversing the dore, as well as aerosol and cloud particle size distributions generally observed by aircraft below 4km and crystal properties generally observed by aircraft above 12km.

Ackerman, A.↗

NASA NTRS1999

Streamwise Vorticity Generation in Laminar and Turbulent Jets

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Demuren, Aodeji O.↗

NASA NTRS1992

Lumley's PODT definition of large eddies and a trio of numerical procedures

Lumley's 1967 Moscow paper provided, for the first time, a completely rational definition of the physically-useful term 'large eddy', popular for a half-century. The numerical procedures based upon his results are: (1) PODT (Proper Orthogonal Decomposition Theorem), which extracts the Large Eddy structure of stochastic processes from physical or computer simulation two-point covariances, and 2) LEIM (Large-Eddy Interaction Model), a predictive scheme for the dynamical large eddies based upon higher order turbulence modeling. Earlier Lumley's work (1964) forms the basis for the final member of the triad of numerical procedures: this predicts the global neutral modes of turbulence which have surprising agreement with both structural eigenmodes and those obtained from the dynamical equations. The ultimate goal of improved engineering design tools for turbulence may be near at hand, partly due to the power and storage of 'supermicrocomputer' workstations finally becoming adequate for the demanding numerics of these procedures.

Payne, Fred R.↗

NASA NTRS2008

Compressible Boundary Layer Predictions at High Reynolds Number using Hybrid LES/RANS Methods

Simulations of compressible boundary layer flow at three different Reynolds numbers (Re(sub delta) = 5.59x10(exp 4), 1.78x10(exp 5), and 1.58x10(exp 6) are performed using a hybrid large-eddy/Reynolds-averaged Navier-Stokes method. Variations in the recycling/rescaling method, the higher-order extension, the choice of primitive variables, the RANS/LES transition parameters, and the mesh resolution are considered in order to assess the model. The results indicate that the present model can provide good predictions of the mean flow properties and second-moment statistics of the boundary layers considered. Normalized Reynolds stresses in the outer layer are found to be independent of Reynolds number, similar to incompressible turbulent boundary layers.

Choi, Jung-Il↗