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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2003

Limits to the Indirect Aerosol Forcing in Stratocumulus

The indirect radiative forcing of aerosols is poorly constrained by the observational data underlying the simple cloud parameterizations in GCMs. signal of cloud response to increased aerosol concentrations from meteorological noise. Recent satellite observations indicate a significant decrease of cloud water in ship tracks, in contrast to an ensemble of in situ measurements showing no average change in cloud water relative to the surrounding clouds. Both results contradict the expectation of cloud water increasing in polluted clouds. We find through large-eddy simulations of stratocumulus that the trend in the satellite data is likely an artifact of sampling only overcast clouds. The simulations instead show cloud cover increasing with droplet concentrations. The simulations also show that increases in cloud water from suppressing drizzle by increased droplet concentrations are favored at night or at extremely low droplet concentrations. At typical droplet concentrations we find that the Twomey effect on cloud albedo is amplified very little by the secondary indirect effect of drizzle suppression, largely because the absorption of solar radiation by cloud water reduces boundary-layer mixing in the daytime and thereby restricts any possible increase in cloud water from drizzle suppression. The cloud and boundary layer respond to radiative heating variations on a time scale of hours, and on longer time scales respond to imbalances between large-scale horizontal advection and the entrainment of inversion air. We analyze the co-varying response of cloud water, cloud thickness, width of droplet size distributions, and dispersion of the optical depth, as well as the overall response of cloud albedo, to changes in droplet concentrations. We also dissect the underlying physical mechanisms through sensitivity studies. Ship tracks represent an ideal natural laboratory to extricate the

Ackerman, Andrew↗
NASA NTRS2002

Studying Turbulence Using Numerical Simulation Databases - IX: Proceedings of the 2002 Summer Program

The ninth Summer Program of the Center for Turbulence Research was held during the period July 29th - August 23rd, 2002. The increase in number of participants, noted in the Preface to the Proceedings of the 2000 Program, continues: this year there were 50 participants from ten countries, and 30 hosts from Stanford and NASA-Ames. This Proceedings volume contains 32 papers that span a wide range of topics and an enormous range of physical scales. The papers have been divided into seven groups: Acoustics, RANS modeling, Combustion, Large-eddy simulation (LES), LES Numerics, Stratified Flows, and Fundamentals, In several cases, a paper could have fitted in more than one group so the classification is somewhat arbitrary.

Bradshaw, Peter↗
DOE OSTI2026

The remarkable inefficiency of stratocumulus

Marine stratocumulus clouds play a central role in Earth's climate system by reflecting incoming solar radiation and exerting a strong cooling effect. Their organization into open and closed mesoscale cellular morphologies can be thought of as an example of bistable dynamics driven by aerosol–cloud interactions and mesoscale processes. From the perspective of non-equilibrium thermodynamics, these structures are an example of a far-from-equilibrium open system that continuously produces and exports entropy. While entropy production has been studied in idealized deep convective systems, it has not yet been quantified for shallow clouds. Here, we compute and decompose the internal entropy production of open- and closed-cell stratocumulus using an ensemble of large-eddy simulations. We show that the overall entropy production of stratocumulus is low, reflecting the limited vertical extent and corresponding reduced ability to utilize the energy fluxes at the system's boundaries. Moist processes dominate the overall irreversibility, which, combined with their low entropy production, leads to a mechanical efficiency about an order of magnitude smaller than in deep convective systems. Although the dominant irreversible processes differ between open- and closed-cell regimes, the distributions of total entropy production largely overlap across the ensemble, limiting the ability to distinguish the dynamics of individual cases based solely on total entropy production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating the Conservation of Energy Variables in Simulations of Deep Moist Convection

Abstract It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy ( B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraft B diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorer B predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Dynamic Smagorinsky model on anisotropic grids

Large Eddy Simulation (LES) of complex-geometry flows often involves highly anisotropic meshes. To examine the performance of the dynamic Smagorinsky model in a controlled fashion on such grids, simulations of forced isotropic turbulence are performed using highly anisotropic discretizations. The resulting model coefficients are compared with a theoretical prediction (Scotti et al., 1993). Two extreme cases are considered: pancake-like grids, for which two directions are poorly resolved compared to the third, and pencil-like grids, where one direction is poorly resolved when compared to the other two. For pancake-like grids the dynamic model yields the results expected from the theory (increasing coefficient with increasing aspect ratio), whereas for pencil-like grids the dynamic model does not agree with the theoretical prediction (with detrimental effects only on smallest resolved scales). A possible explanation of the departure is attempted, and it is shown that the problem may be circumvented by using an isotropic test-filter at larger scales. Overall, all models considered give good large-scale results, confirming the general robustness of the dynamic and eddy-viscosity models. But in all cases, the predictions were poor for scales smaller than that of the worst resolved direction.

Scotti, A.↗
DOE OSTI2020

Improving the Physical Basis for Updraft Dynamics in Deep Convection Parameterizations

This article presents a new deep convective parameterization that determines cloud characteristics based on a specified cloud size distribution. The vertical profiles of cloud properties are determined by analytic equations, which formulate entrainment with an inverse relationship to cloud width. In line with recent studies of large eddy simulations (LES), cloud widths are assumed to be constant with height and vertical mass flux (M) characteristics of the clouds are therefore regulated by the vertical velocity profile. The parameterization is configured to work with existing cloud base M closure formulations, with the closure predicting the total cloud area rather than the cloud base M directly. Analytic formula are also used to connect the vertical wind shear magnitude to the cloud size distribution, wherein larger shear magnitudes result in more numerous large updrafts than weaker shear magnitudes, which is in line with recent research results. The parameterization is compared against 10 deep convective LES with varying thermodynamic and vertical wind shear profiles. Results show dramatic improvements in the prediction of normalized M, detrainment, and the properties of detrained air over the existing Zhang and McFarlane (1995) scheme. In particular, the new model is able to correctly portray the transition from a bottom-heavy M profile in weakly sheared environments, to a top-heavy M profile in strongly sheared environments.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Filtering analysis of a direct numerical simulation of the turbulent Rayleigh-Benard problem

A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.

Eidson, T. M.↗
DOE OSTI2017

Application of Gauss's theorem to quantify localized surface emissions from airborne measurements of wind and trace gases

Airborne estimates of greenhouse gas emissions are becoming more prevalent with the advent of rapid commercial development of trace gas instrumentation featuring increased measurement accuracy, precision, and frequency, and the swelling interest in the verification of current emission inventories. Multiple airborne studies have indicated that emission inventories may underestimate some hydrocarbon emission sources in US oil- and gas-producing basins. Consequently, a proper assessment of the accuracy of these airborne methods is crucial to interpreting the meaning of such discrepancies. We present a new method of sampling surface sources of any trace gas for which fast and precise measurements can be made and apply it to methane, ethane, and carbon dioxide on spatial scales of ~1000 m, where consecutive loops are flown around a targeted source region at multiple altitudes. Using Reynolds decomposition for the scalar concentrations, along with Gauss's theorem, we show that the method accurately accounts for the smaller-scale turbulent dispersion of the local plume, which is often ignored in other average mass balance methods. With the help of large eddy simulations (LES) we further show how the circling radius can be optimized for the micrometeorological conditions encountered during any flight. Furthermore, by sampling controlled releases of methane and ethane on the ground we can ascertain that the accuracy of the method, in appropriate meteorological conditions, is often better than 10 %, with limits of detection below 5 kg h -1 for both methane and ethane. Because of the FAA-mandated minimum flight safe altitude of 150 m, placement of the aircraft is critical to preventing a large portion of the emission plume from flowing underneath the lowest aircraft sampling altitude, which is generally the leading source of uncertainty in these measurements. Finally, we show how the accuracy of the method is strongly dependent on the number of sampling loops and/or time spent sampling the source plume.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
DOE OSTI2018

Response of the Atmospheric Boundary Layer to Submesoscale Sea Surface Temperature Fronts

Submesoscale sea surface temperature fronts are ubiquitous throughout much of the global ocean; however, the response of the marine atmospheric boundary layer (MABL) to the ocean submesoscale is not well understood. Here large-eddy simulation is used to explore the time-dependent response of the MABL to idealized submesoscale sea surface temperature fronts, with an emphasis on how the dynamics of the MABL determine the strength and position of gradients in wind speed and air temperature. Results suggest that horizontal mixing only becomes important in response to frontogenesis by horizontally convergent ageostrophic flows, contrary to the common assumption that the MABL response will be strongly dependent on horizontal turbulent mixing. The fronts that develop in the MABL are also associated with large vertical relative vorticity, suggesting the possibility that submesoscale fronts may induce inertial instability in the MABL. Finally, these results provide guidance for high-resolution ocean and atmosphere modeling and for interpreting observations.

58 GEOSCIENCES↗
NASA NTRS1993

Computation of large-scale statistics in decaying isotropic turbulence

We have performed large-eddy simulations of decaying isotropic turbulence to test the prediction of self-similar decay of the energy spectrum and to compute the decay exponents of the kinetic energy. In general, good agreement between the simulation results and the assumption of self-similarity were obtained. However, the statistics of the simulations were insufficient to compute the value of gamma which corrects the decay exponent when the spectrum follows a k(exp 4) wave number behavior near k = 0. To obtain good statistics, it was found necessary to average over a large ensemble of turbulent flows.

Chasnov, Jeffrey R.↗
NASA NTRS2011

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

Hadjadj, A↗
NASA NTRS2012

Electrical Characteristics of Simulated Tornadoes and Dust Devils

It is well known that tornadoes and dust devils have the ability to accumulate significant, visible clouds of debris. Collisions between sand-like debris species produce different electric charges on different types of grains, which convect along different trajectories around the vortex. Thus, significant charge separations and electric currents are possible, which as the vortex fluctuates over time are thought to produce ULF radiation signatures that have been measured in the field. These electric and magnetic fields may contain valuable information about tornado structure and genesis, and may be critical in driving electrochemical processes within dust devils on Mars. In the present work, existing large eddy simulations of debris-laden tornadoes performed at West Virginia University are coupled with a new debris-charging and advection code developed at Goddard Space Flight Center to investigate the detailed (meter-resolution) fluid-dynamic origins of electromagnetic fields within terrestrial vortices. First results are presented, including simulations of the electric and magnetic fields that would be observed by a near-surface, instrument-laden probe during a direct encounter with a tornado.

Zimmerman, Michael I.↗
DOE OSTI2023

Optimization and Evaluation of Stochastic Unified Convection Using Single-Column Model Simulations at Multiple Observation Sites

We extend the previously developed stochastic unified convection scheme (UNICON) for shallow convection to deep convection by parameterizing the impact of mesoscale organized flow on updraft properties. The extended stochastic UNICON parameterizes thermodynamic properties of updrafts at the near-surface as a multivariate Gaussian distribution, where the variances of the distribution are the summation of variances from non-organized turbulence and mesoscale organized flow. The distribution of updraft radius is parameterized as a power-law distribution with a scale break which is parameterized as a linear function of the strength of mesoscale organized flow. The proposed parameterization is validated using a series of large-eddy simulations of deep convection. The free parameters introduced in the formulation of stochastic UNICON are optimized using 10 cases of single-column model simulations over the ocean. Stochastic UNICON with the optimized parameters significantly reduces the biases of thermodynamic profiles and surface precipitation rates simulated in the original UNICON for tropical convection cases. The simulation of the variation in anomalies of temperature and moisture associated with the Madden-Julian oscillation is also improved. The overall improvements in simulated thermodynamic profiles are found to be due to the increased heating and drying tendencies by convective processes in stochastic UNICON. An additional simulation of an idealized deep convection case shows that stochastic UNICON produces enhanced cloud variabilities with dependency on updraft radius, indicating its ability to represent the coexistence of shallow and deep convection.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

The NASA Low-Pressure Turbine Flow Physics Program

An overview of the NASA Lewis Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The program was established in response to the aero-engine industry's need for improved LPT efficiency and designs. Modern jet engines have four to seven LPT stages, significantly contributing to engine weight. In addition, there is a significant efficiency degradation between takeoff and cruise conditions, of up to 2 points. Reducing the weight and part count of the LPT and minimizing the efficiency degradation will translate into fuel savings. Accurate prediction methods of LPT flows and losses are needed to accomplish those improvements. The flow in LPT passages is at low Reynolds number, and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The affecting parameters traditionally considered are Reynolds number, freestream turbulence intensity, wake frequency parameter, and the pressure distribution (loading). Three-dimensional effects and additional parameters, particularly turbulence characteristics like length scales, spectra and other statistics, as well as wake turbulence intensity and properties also play a role. The flow of most interest is on the suction surface, where large losses are generated as the flow tends to separate at the low Reynolds numbers. Ignoring wakes, a common flow scenario, there is laminar separation, followed by transition on the separation bubble and turbulent reattachment. If transition starts earlier the separation will be eliminated and the boundary layer will be attached leading to the well known bypass transition issues. In contrast, transition over a separation bubble is closer to free shear layer transition and was not investigated as well, particularly in the turbine environment. Unsteadiness created by wakes complicates the picture. Wakes induce earlier transition, and the calmed regions trailing the induced turbulent spots can delay or eliminate separation via shear stress modification. Three-dimensional flow physics and geometry will have strong effects. Altogether a very complex and challenging problem emerges. The objective of the program is to provide improved models and physical understanding of the complex flow, which are essential for accurate prediction of flow and losses in the LPT. Experimental, computational and analytical work as complementing and augmenting approaches are used. The program involves industry, universities and research institutes, and other government laboratories. It is characterized by strong interaction among participants, quick dissemination of results, and responsiveness to industry's needs. The presentation will describe the work elements. Highlighting some activities in progress are experiments on simulated blade suction surface in low-speed wind tunnels, on curved wall, and on a flat-plate, both with pressure gradient. In the area of computation, assessment of existing models is performed using RANS (Reynolds Averaged Navier Stokes) simulations. Laminar flow DNS was completed. Analytical studies of instability and receptivity in attached and separated flows were started. In the near future the program is moving to include wake effects and development of improved modeling. Experimental work in preparation stages are: (1) Addition of wakes to the curved tunnel experiment; (2) Low-speed rotating rig experiment on GE90 engine LPT; and (3) Transonic cascade. In the area of computation, it is expected to move from model assessment towards development of improved models. In addition, a new project of Large Eddy Simulation (LES) of LPT is to begin and will provide numerical data bases. It is planned to implement the emerging improved models in a multistage turbomachinery code and to validate against the GE90 engine LPT.

Ashpis, David E.↗
NASA NTRS2015

RACORO Continental Boundary Layer Cloud Investigations: 3. Separation of Parameterization Biases in Single-Column Model CAM5 Simulations of Shallow Cumulus

Climatically important low-level clouds are commonly misrepresented in climate models. The FAst-physics System TEstbed and Research (FASTER) Project has constructed case studies from the Atmospheric Radiation Measurement Climate Research Facility's Southern Great Plain site during the RACORO aircraft campaign to facilitate research on model representation of boundary-layer clouds. This paper focuses on using the single-column Community Atmosphere Model version 5 (SCAM5) simulations of a multi-day continental shallow cumulus case to identify specific parameterization causes of low-cloud biases. Consistent model biases among the simulations driven by a set of alternative forcings suggest that uncertainty in the forcing plays only a relatively minor role. In-depth analysis reveals that the model's shallow cumulus convection scheme tends to significantly under-produce clouds during the times when shallow cumuli exist in the observations, while the deep convective and stratiform cloud schemes significantly over-produce low-level clouds throughout the day. The links between model biases and the underlying assumptions of the shallow cumulus scheme are further diagnosed with the aid of large-eddy simulations and aircraft measurements, and by suppressing the triggering of the deep convection scheme. It is found that the weak boundary layer turbulence simulated is directly responsible for the weak cumulus activity and the simulated boundary layer stratiform clouds. Increased vertical and temporal resolutions are shown to lead to stronger boundary layer turbulence and reduction of low-cloud biases.

parameterization↗
DOE OSTI2025

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗