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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Large eddy simulations of time-dependent and buoyancy-driven channel flows

The primary goal of this work has been to assess the performance of the dynamic SGS model in the large eddy simulation (LES) of channel flows in a variety of situations, viz., in temporal development of channel flow turned by a transverse pressure gradient and especially in buoyancy-driven turbulent flows such as Rayleigh-Benard and internally heated channel convection. For buoyancy-driven flows, there are additional buoyant terms that are possible in the base models, and one objective has been to determine if the dynamic SGS model results are sensitive to such terms. The ultimate goal is to determine the minimal base model needed in the dynamic SGS model to provide accurate results in flows with more complicated physical features. In addition, a program of direct numerical simulation (DNS) of fully compressible channel convection has been undertaken to determine stratification and compressibility effects. These simulations are intended to provide a comparative base for performing the LES of compressible (or highly stratified, pseudo-compressible) convection at high Reynolds number in the future.

Cabot, William H.↗
NASA NTRS1993

Large-eddy simulation of flow in a plane, asymmetric diffuser

Recent improvements in subgrid-scale modeling as well as increases in computer power make it feasible to investigate flows using large-eddy simulation (LES) which have been traditionally studied with techniques based on Reynolds averaging. However, LES has not yet been applied to many flows of immediate technical interest. Preliminary results from LES of a plane diffuser flow are described. The long term goal of this work is to investigate flow separation as well as separation control in ducts and ramp-like geometries.

Kaltenbach, Hans-Jakob↗
NASA NTRS1982

An experimental search for near-wall boundary conditions for large eddy simulation

Instantaneous wall shear stress and streamwise velocities have been measured simultaneously in a flat plate, turbulent boundary layer at moderate Reynolds number in an effort to provide experimental support for large eddy simulations. Data were obtained by using a buried-wire wall shear gage and a hot-wire rake positioned in the log region of the flow. All data processing was accomplished with digital data analysis techniques on a minicomputer. Fluctuations of the instantaneous U plus versus Y plus profiles about a mean law of the wall are shown to be significant and complex. Peak cross-correlation values between wall shear stress and the velocities are high and reflect the passage of a large structure inclined at a small angle to the wall. Estimates of this angle are consistent with those made by other investigators. Conditional sampling techniques were used to detect the passage of various sizes and types of flow disturbances (events) and to estimate their mean frequency of occurrence. Events characterized by large and sudden streamwise accelerations were found to be highly coherent throughout the log region and were strongly correlated with large fluctuations in wall shear-stress. Phase randomness between the near-wall quantities and the outer velocities was small. The results suggest that the flow events detected by conditional sampling applied to velocities in the log region may be related to the bursting process.

Robinson, S. K.↗
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 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 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↗
NASA NTRS1993

A large eddy simulation scheme for turbulent reacting flows

The recent development of the dynamic subgrid-scale (SGS) model has provided a consistent method for generating localized turbulent mixing models and has opened up great possibilities for applying the large eddy simulation (LES) technique to real world problems. Given the fact that the direct numerical simulation (DNS) can not solve for engineering flow problems in the foreseeable future (Reynolds 1989), the LES is certainly an attractive alternative. It seems only natural to bring this new development in SGS modeling to bear on the reacting flows. The major stumbling block for introducing LES to reacting flow problems has been the proper modeling of the reaction source terms. Various models have been proposed, but none of them has a wide range of applicability. For example, some of the models in combustion have been based on the flamelet assumption which is only valid for relatively fast reactions. Some other models have neglected the effects of chemical reactions on the turbulent mixing time scale, which is certainly not valid for fast and non-isothermal reactions. The probability density function (PDF) method can be usefully employed to deal with the modeling of the reaction source terms. In order to fit into the framework of LES, a new PDF, the large eddy PDF (LEPDF), is introduced. This PDF provides an accurate representation for the filtered chemical source terms and can be readily calculated in the simulations. The details of this scheme are described.

Gao, Feng↗
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↗
NASA NTRS1995

A family of dynamic models for large-eddy simulation

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Carati, D.↗
NASA NTRS2001

Large Eddy Simulation of Gravitational Effects on Transitional and Turbulent Gas-Jet Diffusion Flames

The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.

Givi, Peyman↗
NASA NTRS1999

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗
NASA NTRS1993

Large-eddy simulation of turbulent flow in a channel with a surface-mounted two-dimensional obstacle using a dynamic subgrid-scale model

We apply the dynamic subgrid-scale model (DSGSM) to large-eddy simulation (LES) of turbulent channel flow with a square rib mounted on one wall. The Reynolds number (Re) is 3210 based on the mean velocity above the obstacle and the obstacle height. Near-wall structures are resolved with the no-slip boundary condition. The results show better agreement with direct numerical simulation (DNS) than LES with a fixed model constant, verifying the value of the DSGSM for simulating complex turbulent flows.

Yang, Kyung-Soo↗
NASA NTRS2003

Large Eddy Simulation Of Gravitational Effects In Transitional And Turbulent Gas-Jet Diffusion Flames

The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.

Jaberi, Farhad A.↗
NASA NTRS1993

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.↗
DOE OSTI2024

Multicycle large-eddy simulations of a direct-injection hydrogen-fueled optical engine

Hydrogen (H 2 ) is a carbon-free chemical energy carrier and one promising solution for achieving effective decarbonization of the transportation sector, particularly for internal combustion engines (ICEs). With a focus on ICEs, and compared to port-fuel injection, direct injection (DI) of gaseous H 2 during the compression stroke offers potential advantages, which include backfire avoidance and reduction of preignition occurrence. In these last two decades, much research, experimental and numerical, has been devoted to understanding H 2 's mixing and combustion processes in ICEs. Computational fluid dynamics modeling efforts commonly rely on unsteady Reynolds-averaged Navier Stokes (URANS) turbulence frameworks, mostly due to their computational affordability. However, many authors have pointed out the opportunity to perform large-eddy simulations (LESs) to investigate the cyclic variability of H 2 engines and assess potential advantages of using LES in place of URANS, especially for lean operation. This study addresses this knowledge gap and presents a computational fluid dynamics (CFD) study of the H 2 DI process in an optical engine operating at relatively low tumble conditions, using multicycle LESs. In conclusion, the manuscript presents a thorough validation of the results against experimental data available from the literature as well as direct comparison with URANS, demonstrating the feasibility of multicycle LESs for CFD modeling of DI H 2 -fueled ICEs.

Direct injection↗
DOE OSTI2024

Diurnal evolution of non-precipitating marine stratocumuli in a large-eddy simulation ensemble

Abstract. We explore the cloud system evolution of non-precipitating marine stratocumuli with a focus on the impacts of the diurnal cycle and free-tropospheric (FT) humidity based on an ensemble of 244 large-eddy simulations generated by perturbing initial thermodynamic profiles and aerosol conditions. Cases are categorized based on their degree of decoupling and the cloud liquid water path (LWPc, based on model columns with cloud optical depths greater than one). A budget analysis method is proposed to analyze the evolution of cloud water in both coupled and decoupled boundary layers. More coupled clouds start with a relatively low LWPc and cloud fraction (fc) but experience the least decrease in LWPc and fc during the daytime. More decoupled clouds undergo greater daytime reduction in LWPc and fc, especially those with higher LWPc at sunrise because they suffer from faster weakening of net radiative cooling. During the nighttime, a positive correlation between FT humidity and the LWPc emerges, consistent with higher FT humidity reducing both radiative cooling and the humidity jump, both of which reduce entrainment and increase LWPc. The LWPc is more likely to decrease during the nighttime for a larger LWPc and greater inversion base height (zi), conditions under which entrainment dominates as turbulence develops. In the morning, the rate of the LWPc reduction depends on the LWPc at sunrise, zi, and the degree of decoupling, with distinct contributions from subsidence and radiation.

Chen, Yao-Sheng (ORCID:0000000208354132)↗
NASA NTRSDate not supplied

Wall-Modeled Large-Eddy Simulations of a Swept Wing with Leading-Edge Ice

The aerodynamic performance of a swept wing derived from the Common Research Model wing with high-fidelity and smooth leading-edge ice shapes is evaluated using wall-modeled large-eddy simulations. Two mesh paradigms are considered within the Launch, Ascent, and Vehicle Aerodynamics framework: structured hybrid curvilinear body-fitted/immersed-boundary and unstructured Voronoi. The curvilinear approach models the effect of the ice on the flow through penalty terms and the Voronoi strategy includes the ice representation directly in the mesh.Aerodynamic loads, surface pressure profiles, and skin-friction streamlines are compared against experimental results collected at the Wichita State University wind tunnel at a Mach number of 0.18 and a Reynolds number per mean aerodynamic chord of 1.6 million. For the curvilinear simulations with high-fidelity ice, results are sensitive to spanwise mesh resolution near the leading edge, where insufficient spanwise resolution leads to a spanwise running vortexaft of the high-fidelity ice, which is partially removed with increasing spanwise resolution. For the Voronoi simulations with high-fidelity ice, refining the mesh led to good convergence towards the experimental results with the best case demonstrating a maximum discrepancy of three lift counts relative to the experiment across a range of angles of attack from about 5 to 24 degrees, validating the shrink-wrapping procedure used for the ice. For the smooth ice, the curvilinear approach consistently under-predicted the upper surface suction leading to an under-prediction in lift, but did not demonstrate the pronounced over-prediction in upper surface suction near the leading edge observed in the Voronoi simulations.

AATT↗