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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗

DOE OSTI2025

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗

NASA NTRS1993

Parallel spatial direct numerical simulations on the Intel iPSC/860 hypercube

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube is documented. The direct numerical simulation approach is used to compute spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows. The feasibility of using the PSDNS on the hypercube to perform transition studies is examined. The results indicate that the direct numerical simulation approach can effectively be parallelized on a distributed-memory parallel machine. By increasing the number of processors nearly ideal linear speedups are achieved with nonoptimized routines; slower than linear speedups are achieved with optimized (machine dependent library) routines. This slower than linear speedup results because the Fast Fourier Transform (FFT) routine dominates the computational cost and because the routine indicates less than ideal speedups. However with the machine-dependent routines the total computational cost decreases by a factor of 4 to 5 compared with standard FORTRAN routines. The computational cost increases linearly with spanwise wall-normal and streamwise grid refinements. The hypercube with 32 processors was estimated to require approximately twice the amount of Cray supercomputer single processor time to complete a comparable simulation; however it is estimated that a subgrid-scale model which reduces the required number of grid points and becomes a large-eddy simulation (PSLES) would reduce the computational cost and memory requirements by a factor of 10 over the PSDNS. This PSLES implementation would enable transition simulations on the hypercube at a reasonable computational cost.

Joslin, Ronald D.↗

NASA NTRS1991

Simulations of curved turbulent boundary layers

The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Orlandi, Paolo↗

NASA NTRS1991

Simulation of the inertial-conductive subrange

The inertial-conductive subrange spectrum of a passive temperature field of a fluid of small Prandtl number is determined by large-eddy simulation. Results for simulations of both freely decaying and forced turbulence are presented. In the simulations of freely decaying turbulence, a subgrid model is used to simulate a decaying turbulent velocity field with a k exp -5/3 inertial subrange energy spectrum convecting eight different decaying temperature fields with well-resolved conductive subranges. In the simulations of forced turbulence, a subgrid model is again used to simulate an inertial subrange velocity field; however, the velocity field is now forced at the lowest wave numbers of the simulation, and an external uniform mean temperature gradient is imposed. Statistically stationary velocity and temperature fluctuations are generated. The results of the decaying and forced simulations are in excellent agreement with the Batchelor et al. (1959) k exp -17/3 spectrum in the far inertial-conductive subrange, whereas significant departures from the Batchelor spectrum are observed in the near inertial-conductive subrange.

Chasnov, Jeffrey R.↗

DOE OSTI2024

The Role of the Toroidal Vortex in Cumulus Clouds' Entrainment and Mixing

Abstract Shallow convective clouds play a crucial role in Earth's energy budget, as they modulate the radiative transfer in the atmosphere and participate in the vertical transport of aerosols, energy, and humidity. The parameterizations representing these complex, vital players in weather and climate models are mostly based on a description of steady‐state plumes and are a source of major uncertainty. Recently, several studies have shown that buoyant thermals are inherent in atmospheric convection and contain a toroidal (ring) vortex. This work studies those vortices in growing shallow cumulus (Cu) clouds using high‐resolution (10 m) Large Eddy Simulations that resolve these vortices in much detail. Recent analysis of such data showed that small‐scale turbulent diffusion is unable to explain the large diluted portion of the cloud. Here we advocate for the important role of the Cu toroidal vortex (TV) in cloud dilution and present the complex dynamics and structure of a Cu TV. Nevertheless, since the vortex dominates the cloud's dilution, simplicity emerges when considering the cloud's lateral mass flux profile. The cloud mixing is quantified using direct flux calculations and Eulerian tracers. In addition, Lagrangian tracers are used to identify the origin of the entrained air and its thermodynamic properties. It shows that most of the air entrained by the vortex is not recycled by the vortex, yet is significantly more humid than the environment. We suggest that the development of new models describing thermals, together with their toroidal vortices, might improve cloud parameterizations in weather and climate models.

Meteorology & Atmospheric Sciences↗

NASA NTRS1984

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗

DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗

DOE OSTI2026

SANDWake3D

SAND2026-23958O SANDWake3D models atmospheric surface layers and turbine wakes using a 3D parabolic RANS solver. This Python-based tool uses an efficient method to simulate wakes, capturing complex wake behaviors such as wind shear, veer, and atmospheric stratification, naturally. It has been calibrated and validated against higher-fidelity large eddy simulations (LES), providing reliable results for wind energy and atmospheric research applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

NASA NTRSDate not supplied

Aeroacoustic Computations of a Generic Low Boom Concept in Landing Configuration: Part 2 - Airframe Noise Simulations

Aeroacoustic simulations of a sub-scale, generic, low-boom supersonic aircraft model at low speeds were carried out using the Lattice-Boltzmann Very Large Eddy Simulation software PowerFLOW® and a Ffowcs Williams and Hawkings approach. The objective was to characterize the airframe noise signature of the configuration and determine the similarities and differences to conventional subsonic aircraft. We quantified the noise produced by the flap and landing gear by simulating the aircraft with and without deployment of these components. Mach number scaling of the far-field noise levels was examined by performing simulations at different flow speeds. We demonstrate that a solid surface formulation does not capture the aircraft acoustic field properly, while a carefully constructed permeable surface formulation yields far-field spectra that compare very favorably to results from direct probing of the simulated flow field.

low boom↗

NASA NTRS2019

Exit Boundary Layer Data for a Round Convergent Nozzle in Support of Numerical Simulation Efforts

An experiment is conducted with hot-wire anemometry to document the exit boundary layer characteristics of two nozzle configurations at jet Mach numbers up to 0.82. Far-field noise and jet plume experimental data from these two configurations have been used in Large Eddy Simulations (LES) of jets by colleagues at other Institutions. The current experiment provides the boundary layer data which have been identified as being critical for validation of the simulations since the initial conditions can significantly affect subsequent jet evolution and its radiated noise. The data exhibit fully turbulent boundary layers for the case with a pipe attached upstream of the nozzle. The case without the pipe involves ‘Blasius-like’ mean velocity profiles but a ‘highly disturbed laminar state’ with large turbulence intensities in a range of subsonic Mach numbers.

Transition↗

NASA NTRS2019

Exit Boundary Layer Data for a Round Convergent Nozzle in Support of Numerical Simulation Efforts

Transition↗

NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗

NASA NTRSDate not supplied

CFD Predictions of Fluctuating Pressure Environments on NASA’s SLS Ascent Unsteady Aerodynamics Wind Tunnel Test

This paper details work that is being conducted on predicting fluctuating pressure environments on NASA launch vehicles using computational fluid dynamics (CFD). The accurate characterization of these aeroacoustics environments is necessary in the analysis of vehicle structural health and aerodynamic performance and in the vehicle design process. High-fidelity simulations of the unsteady flow over NASA’s SLS Ascent Unsteady Aerodynamics wind tunnel Test (AUAT) using the hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulations (LES) methodologies in Mississippi States’ Loci/CHEM solver are presented. The two available methods, Nelson-Nichols and dynamic hybrid RANS-LES (DHRL), are compared head-to-head on identical grids at two Mach numbers. At the high-subsonic Mach number, locally supersonic flow expanding over a shoulder induces a separation-reattachment system that is predicted poorly by the Nelson-Nichols method due to an under-prediction in unsteady content. The DHRL solution predicts the spatial extent and frequency domain response of this system well but under-predicts the peak noise levels due to a delay in the RANS-to-LES transition. At the low-supersonic Mach number, compression corner dynamics dominate the unsteady flow. While the Nelson Nichols solution captures this phenomenon well, another solution on a refined grid indicates that some model stress depletion is occurring. The DHRL solution again performs well at this Mach number and matches the wind tunnel data fairly closely, but perhaps most impressively, is fairly agnostic to increasing grid resolution. The findings in this paper are expected to be generally applicable and will be used to guide the prediction of aeroacoustics environments of other NASA launch vehicles.

CFD↗

NASA NTRS2008

Hybrid LES/RANS Simulation of Transverse Sonic Injection into a Mach 2 Flow

A computational study of transverse sonic injection of air and helium into a Mach 1.98 cross-flow is presented. A hybrid large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) turbulence model is used, with the two-equation Menter baseline (Menter-BSL) closure for the RANS part of the flow and a Smagorinsky-type model for the LES part of the flow. A time-dependent blending function, dependent on modeled turbulence variables, is used to shift the closure from RANS to LES. Turbulent structures are initiated and sustained through the use of a recycling / rescaling technique. Two higher-order discretizations, the Piecewise Parabolic Method (PPM) of Colella and Woodward, and the SONIC-A ENO scheme of Suresh and Huyhn are used in the study. The results using the hybrid model show reasonably good agreement with time-averaged Mie scattering data and with experimental surface pressure distributions, even though the penetration of the jet into the cross-flow is slightly over-predicted. The LES/RANS results are used to examine the validity of commonly-used assumptions of constant Schmidt and Prandtl numbers in the intense mixing zone downstream of the injection location.

Boles, John A.↗

NASA NTRS2014

DNS of Flows over Periodic Hills using a Discontinuous-Galerkin Spectral-Element Method

Direct numerical simulation (DNS) of turbulent compressible flows is performed using a higher-order space-time discontinuous-Galerkin finite-element method. The numerical scheme is validated by performing DNS of the evolution of the Taylor-Green vortex and turbulent flow in a channel. The higher-order method is shown to provide increased accuracy relative to low-order methods at a given number of degrees of freedom. The turbulent flow over a periodic array of hills in a channel is simulated at Reynolds number 10,595 using an 8th-order scheme in space and a 4th-order scheme in time. These results are validated against previous large eddy simulation (LES) results. A preliminary analysis provides insight into how these detailed simulations can be used to improve Reynoldsaveraged Navier-Stokes (RANS) modeling

DNS↗

NASA NTRS2014

Pre-Test CFD for the Design and Execution of the Enhanced Injection and Mixing Project at NASA Langley Research Center

With the increasing costs of physics experiments and simultaneous increase in availability and maturity of computational tools it is not surprising that computational fluid dynamics (CFD) is playing an increasingly important role, not only in post-test investigations, but also in the early stages of experimental planning. This paper describes a CFD-based effort executed in close collaboration between computational fluid dynamicists and experimentalists to develop a virtual experiment during the early planning stages of the Enhanced Injection and Mixing project at NASA Langley Research Center. This projects aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships relevant to flight Mach numbers greater than 8. The purpose of the virtual experiment was to provide flow field data to aid in the design of the experimental apparatus and the in-stream rake probes, to verify the nonintrusive measurements based on NO-PLIF, and to perform pre-test analysis of quantities obtainable from the experiment and CFD. The approach also allowed for the joint team to develop common data processing and analysis tools, and to test research ideas. The virtual experiment consisted of a series of Reynolds-averaged simulations (RAS). These simulations included the facility nozzle, the experimental apparatus with a baseline strut injector, and the test cabin. Pure helium and helium-air mixtures were used to determine the efficacy of different inert gases to model hydrogen injection. The results of the simulations were analyzed by computing mixing efficiency, total pressure recovery, and stream thrust potential. As the experimental effort progresses, the simulation results will be compared with the experimental data to calibrate the modeling constants present in the CFD and validate simulation fidelity. CFD will also be used to investigate different injector concepts, improve understanding of the flow structure and flow physics, and develop functional relationships. Both RAS and large eddy simulations (LES) are planned for post-test analysis of the experimental data.

Drozda, Tomasz G.↗