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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Unstructured LES of Reacting Multiphase Flows in Realistic Gas Turbine Combustors

As part of the Accelerated Strategic Computing Initiative (ASCI) program, an accurate and robust simulation tool is being developed to perform high-fidelity LES studies of multiphase, multiscale turbulent reacting flows in aircraft gas turbine combustor configurations using hybrid unstructured grids. In the combustor, pressurized gas from the upstream compressor is reacted with atomized liquid fuel to produce the combustion products that drive the downstream turbine. The Large Eddy Simulation (LES) approach is used to simulate the combustor because of its demonstrated superiority over RANS in predicting turbulent mixing, which is central to combustion. This paper summarizes the accomplishments of the combustor group over the past year, concentrating mainly on the two major milestones achieved this year: 1) Large scale simulation: A major rewrite and redesign of the flagship unstructured LES code has allowed the group to perform large eddy simulations of the complete combustor geometry (all 18 injectors) with over 100 million control volumes; 2) Multi-physics simulation in complex geometry: The first multi-physics simulations including fuel spray breakup, coalescence, evaporation, and combustion are now being performed in a single periodic sector (1/18th) of an actual Pratt & Whitney combustor geometry.

Ham, Frank↗
NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗
NASA NTRS2002

The Role of Turbulence in Chemical and Dynamical Processes in the Near-Field Wake of Subsonic Aircraft

During this grant, covering the period from September 1998 to December 2001, we continued the investigation of the role of turbulent mixing in the wake of subsonic aircraft initiated in 1994 for NASA's Atmospheric Effects of Aviation Project. The goal of the research has been to provide sufficient understanding and quantitative analytical capability to assess the dynamical, chemical, and microphysical interactions in the near-field wake that have the greatest potential to influence the global atmospheric impact of the projected fleet of subsonic aircraft. Through large-eddy simulations we have shown that turbulence in the early wake dynamics can have a strong effect on both the ice microphysics of contrail evolution and on wake chemistry. The wake vortex dynamics are the primary determinant of the vertical extent of the contrail; this together with the local wind shear largely determines the horizontal extent. The fraction of the initial ice crystals surviving the wake vortex dynamics, their spatial distribution, and the ice mass distribution are all sensitive to the aircraft type, assumed initial ice crystal number, and ambient humidity and turbulence conditions. Our model indicates that there is a significant range of conditions for which a smaller aircraft such as a B737 produces as significant a persistent contrail as a larger aircraft such as a B747, even though the latter consumes almost five times as much fuel. Large-eddy simulations of the near wake of a B757 provided a fine-grained chemical-dynamical representation of simplified NOx - HOx chemistry in wakes of ages from a few seconds to several minutes. By sampling the simulated data in a manner similar to that of in situ aircraft measurements it was possible to provide a likely explanation for a puzzle uncovered in the 1996 SUCCESS flight measurements of OH and HO2 The results illustrate the importance of considering fluid dynamics effects in interpreting chemistry results when mixing rates and species fluctuations are large, and demonstrate the feasibility of using 3D unsteady LES with coupled chemistry to study such phenomena.

Lewellen, D. C.↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre-and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre-and post-buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the span wise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗
NASA NTRS2023

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES),and wall-modeled large-eddy simulation (WMLES)are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985).1The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model.The unsteady RANS simulations,only with the compressibility correction,predict the pre- and post-buffet characteristics,which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment.The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment

Transonic Buffet↗
NASA NTRS2010

Near Stall Flow Analysis in the Transonic Fan of the RTA Propulsion System

Turbine-based propulsion systems for access to space have been investigated at NASA Glenn Research center. A ground demonstrator engine for validation testing has been developed as a part of the program. The demonstrator, the Revolutionary Turbine Accelerator (RTA-1), is a variable cycle turbofan ramjet designed to transition from an augmented turbofan to a ramjet that produces the thrust required to accelerate the vehicle to Mach 4. The RTA-1 is designed to accommodate a large variation in bypass ratio from sea level static to Mach 4 flight condition. A key component of this engine is a new fan stage that accommodates these large variations in bypass ratio and flow ranges. In the present study, unsteady flow behavior in the fan of the RTA-1 is studied in detail with large eddy simulation (LES) and the numerical results are compared with measured data. During the experimental study of the fan stage, humming sound was detected at 100 % speed near stall operation. The main purpose of the study is to investigate details of the unsteady flow behavior at near stall operation and to identify a possible cause of the hum. The large eddy simulation of the current flow field reproduces main features of the measured flow very well. The LES simulation indicates that non-synchronous flow instability develops as the fan operates toward the stall limit. The FFT analysis of the calculated wall pressure shows that the rotating flow instability has the characteristic frequency that is about 50% of the blade passing frequency.

Hah, Chunill↗
NASA NTRS2014

Review of Idealized Aircraft Wake Vortex Models

Properties of three aircraft wake vortex models, Lamb-Oseen, Burnham-Hallock, and Proctor are reviewed. These idealized models are often used to initialize the aircraft wake vortex pair in large eddy simulations and in wake encounter hazard models, as well as to define matched filters for processing lidar observations of aircraft wake vortices. Basic parameters for each vortex model, such as peak tangential velocity and circulation strength as a function of vortex core radius size, are examined. The models are also compared using different vortex characterizations, such as the vorticity magnitude. Results of Euler and large eddy simulations are presented. The application of vortex models in the postprocessing of lidar observations is discussed.

Ahmad, Nashat N.↗
NASA NTRS2008

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗
NASA NTRS2018

Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Comparisons with experimental data showed a sensitivity of the current prediction to the modeling of the sidewalls. This was found to be common among various computational studies in the literature where periodic boundary conditions were used in the spanwise direction, as was the case in the present work. Thus, although the experiment was quasi-two-dimensional, the present simulation was determined to be two-dimensional. Quantities present in the exact equation of the Reynolds-stress transport, i.e., production, molecular diffusion, turbulent transport, pressure diffusion, pressure strain, dissipation, and turbulent mass flux were calculated. Reynolds-stress budgets were compared with past large-eddy simulation and direct numerical simulation datasets in the undisturbed portion of the turbulent boundary layer to validate the current approach. The budgets in SBLI showed the growth in the production term for the primary normal stress and energy transfer mechanism was led by the pressure strain term in the secondary normal stresses. The pressure diffusion term, commonly assumed as negligible by turbulence model developers, was shown to be small but non-zero in the normal stress budgets, however it played a key role in the primary shear stress budget.

shock wave/boundary-layer interactio↗
NASA NTRS2014

Extending 3D Near-Cloud Corrections from Shorter to Longer Wavelengths

Satellite observations have shown a positive correlation between cloud amount and aerosol optical thickness (AOT) that can be explained by the humidification of aerosols near clouds, and/or by cloud contamination by sub-pixel size clouds and the cloud adjacency effect. The last effect may substantially increase reflected radiation in cloud-free columns, leading to overestimates in the retrieved AOT. For clear-sky areas near boundary layer clouds the main contribution to the enhancement of clear sky reflectance at shorter wavelengths comes from the radiation scattered into clear areas by clouds and then scattered to the sensor by air molecules. Because of the wavelength dependence of air molecule scattering, this process leads to a larger reflectance increase at shorter wavelengths, and can be corrected using a simple two-layer model. However, correcting only for molecular scattering skews spectral properties of the retrieved AOT. Kassianov and Ovtchinnikov proposed a technique that uses spectral reflectance ratios to retrieve AOT in the vicinity of clouds; they assumed that the cloud adjacency effect influences the spectral ratio between reflectances at two wavelengths less than it influences the reflectances themselves. This paper combines the two approaches: It assumes that the 3D correction for the shortest wavelength is known with some uncertainties, and then it estimates the 3D correction for longer wavelengths using a modified ratio method. The new approach is tested with 3D radiances simulated for 26 cumulus fields from Large-Eddy Simulations, supplemented with 40 aerosol profiles. The results showed that (i) for a variety of cumulus cloud scenes and aerosol profiles over ocean the 3D correction due to cloud adjacency effect can be extended from shorter to longer wavelengths and (ii) the 3D corrections for longer wavelengths are not very sensitive to unbiased random uncertainties in the 3D corrections at shorter wavelengths.

aerosol optical depth↗
NASA NTRS1979

Developments in the computation of turbulent boundary layers

Computational techniques applicable to turbulent boundary layers are classified into solutions of Reynolds-averaged equations, in which all the effects of the turbulence are modelled, and solutions of three-dimensional, time dependent Navier-Stokes equations, in which the large eddies are calculated and only the turbulence at scales smaller than the computational mesh spacings has to be modelled. Current computation costs place engineering computations in the first of these categories; large eddy simulations are appropriate currently for special studies of the dynamical processes of turbulence in idealized flow fields. It is shown that the two methods are interrelated and that each can gain from advances in the other. The degree of success of a pair of increasingly complex Reynolds stress models to broaden their range of applicability is examined through comparisons with experimental data for a variety of flow conditions. An example of a large-eddy simulation is presented, compared with experimental results, and used to evaluate the models for pressure rate-of-strain correlations and dissipation in the Reynolds-averaged equations.

Rubesin, M. W.↗
NASA NTRS1990

On the Subgrid-Scale Modeling of Compressible Turbulence

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Squires, Kyle↗
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2024

Simplified Approximations of Direct Cumulus Entrainment and Detrainment

Abstract In recent years, direct calculations of simulated cumulus entrainment and detrainment have facilitated new physical insights into these highly elusive but critically important processes. However, these calculations require substantial computational resources that may limit their widespread usage. To facilitate such calculations, two simplified approximations of direct cumulus entrainment and detrainment are examined herein. The first approximation, termed the “semidirect” method, follows a standard bulk approach but makes more realistic assumptions about the sources of entrained and detrained air near the cloud edges. In contrast, the second approximation (the “projection” method) uses the governing equations of motion to project whether grid points near the cloud edge will entrain or detrain as the mean cloud ascends by one grid point. Verification exercises using large-eddy simulations reveal that both methods generally agree better with corresponding direct entrainment/detrainment estimates than the traditional bulk formulation, with the projection method outperforming the semidirect method. The two methods can be used in a synergistic fashion, with the semidirect method helping to optimize the projection method, to suit a wide range of applications. Because the latter incorporates the essential dynamics of entrainment and detrainment at the local scale, it can be used to gain physical insight into the causal mechanisms regulating these complex processes.

Meteorology & Atmospheric Sciences↗
NASA NTRS2016

Overview of the NASA Glenn Flux Reconstruction Based High-Order Unstructured Grid Code

A computational fluid dynamics code based on the flux reconstruction (FR) method is currently being developed at NASA Glenn Research Center to ultimately provide a large- eddy simulation capability that is both accurate and efficient for complex aeropropulsion flows. The FR approach offers a simple and efficient method that is easy to implement and accurate to an arbitrary order on common grid cell geometries. The governing compressible Navier-Stokes equations are discretized in time using various explicit Runge-Kutta schemes, with the default being the 3-stage/3rd-order strong stability preserving scheme. The code is written in modern Fortran (i.e., Fortran 2008) and parallelization is attained through MPI for execution on distributed-memory high-performance computing systems. An h- refinement study of the isentropic Euler vortex problem is able to empirically demonstrate the capability of the FR method to achieve super-accuracy for inviscid flows. Additionally, the code is applied to the Taylor-Green vortex problem, performing numerous implicit large-eddy simulations across a range of grid resolutions and solution orders. The solution found by a pseudo-spectral code is commonly used as a reference solution to this problem, and the FR code is able to reproduce this solution using approximately the same grid resolution. Finally, an examination of the code's performance demonstrates good parallel scaling, as well as an implementation of the FR method with a computational cost/degree- of-freedom/time-step that is essentially independent of the solution order of accuracy for structured geometries.

High-Order Methods↗
NASA NTRS1992

Physical processes within the nocturnal stratus-topped boundary layer

There are many physical processes involved in the stratus-topped boundary layer: longwave radiation cooling, entrainment, latent heating, surface heating, solar heating, and drizzling, for example. The manner in which the processes combine to maintain the turbulence within the stratus-topped boundary layer remains an unsolved problem. The large eddy simulation technique is used to examine the first four physical processes mentioned above. First, the contribution of each physical process to the thermodynamic differences between the updraft and downdraft branches of turbulent circulations is examined through a conditional sampling. Second, these mean thermodynamic differences are shown to express well the vertical distributions of heat and moisture fluxes within stratus-topped boundary layers. These provide a method to validate the process partitioning technique. (This technique assumes that the net flux profile can be partitioned into different component-flux profiles according to physical processes, and that each partitioned component flux is linear in height.) In this paper, the heat and moisture fluxes are process partitioned, and each component-flux is found to contribute to the net flux in a way that is consistent with its corresponding process contribution to the mean thermodynamic differences between updrafts and downdrafts. Also, the net flux obtained by summing all component-fluxes agrees very well with that obtained directly from the large-eddy simulations.

Moeng, Chin-Hoh↗