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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Cloud Adjustments From Large-Scale Smoke–Circulation Interactions Strongly Modulate the Southeastern Atlantic Stratocumulus-to-Cumulus Transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (∼1–2 K) and greatly enhanced moisture (>2 g kg−1) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

southeast Atlantic Ocean↗
NASA NTRS1993

Dynamic subgrid-scale modeling for high-speed transitional boundary layers

The subgrid scales are modeled dynamically in a large-eddy simulation of transitional boundary-layer flow along a hollow cylinder at a Mach number of 4.5. The behavior of the dynamic-model coefficients, which is determined from local information in the resolved field, is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both contractions proposed by Germano et al. (1991) and Lilly (1992) are used for the unique determination of the coefficients of the dynamic model, and their results are compared and assessed. The behavior, as well as the energy cascade of the subgridscale field structure, is investigated at various stages of the transition process.

El-Hady, Nabil M.↗
NASA NTRS1979

Evaluation of subgrid-scale models using an accurately simulated turbulent flow

A calculation of periodic homogeneous isotropic turbulence is used to simulate the experimental decay of grid turbulence. The calculation is found to match the experiment in a number of important aspects and the computed flow field is then treated as a realization of a physical turbulent flow. From this flow, a calculation is conducted of the large eddy field and the various averages of the subgrid-scale turbulence that occur in the large eddy simulation equations. These quantities are compared with the predictions of the models that are usually applied in large eddy simulation. The results show that the terms which involve the large-scale field are accurately modeled but the subgrid-scale Reynolds stresses are only moderately well modeled. It is also possible to use the method to predict the constants of the models without reference to experiment. Attempts to find improved models have not met with success.

Clark, R. A.↗
NASA NTRS1996

A New Class of Finite Difference Schemes

Fluid flows in the transitional and turbulent regimes possess a wide range of length and time scales. The numerical computation of these flows therefore requires numerical methods that can accurately represent the entire, or at least a significant portion, of this range of scales. The inaccurate representation of small scales is inherent to non-spectral schemes. This can be detrimental to computations where the energy in the small scales is comparable to that in the larger scales, e.g. large-eddy simulations of high Reynolds number turbulence. The inaccurate numerical representation of the small scales in these large-eddy simulations can result in the numerical error overwhelming the contribution of the subgrid-scale model.

Mahesh, K.↗
NASA NTRS2001

Effects of Aircraft Wake Dynamics on Measured and Simulated NO(x) and HO(x) Wake Chemistry

High-resolution numerical large-eddy simulations of the near wake of a B757 including simplified NOx and HOx chemistry were performed to explore the effects of dynamics on chemistry in wakes of ages from a few seconds to several minutes. Dilution plays an important basic role in the NOx-O3 chemistry in the wake, while a more interesting interaction between the chemistry and dynamics occurs for the HOx species. These simulation results are compared with published measurements of OH and HO2 within a B757 wake under cruise conditions in the upper troposphere taken during the Subsonic Aircraft Contrail and Cloud Effects Special Study (SUCCESS) mission in May 1996. The simulation provides a much finer grained representation of the chemistry and dynamics of the early wake than is possible from the 1 s data samples taken in situ. The comparison suggests that the previously reported discrepancy of up to a factor of 20 - 50 between the SUCCESS measurements of the [HO2]/[OH] ratio and that predicted by simplified theoretical computations is due to the combined effects of large mixing rates around the wake plume edges and averaging over volumes containing large species fluctuations. The results demonstrate the feasibility of using three-dimensional unsteady large-eddy simulations with coupled chemistry to study such phenomena.

Lewellen, D. C.↗
NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗
NASA NTRS1992

LES versus DNS: A comparative study

We have performed Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of forced isotropic turbulence at moderate Reynolds numbers. The subgrid scale model used in the LES is based on an eddy viscosity which adjusts instantaneously the energy spectrum of the LES to that of the DNS. The statistics of the large scales of the DNS (filtered DNS field or fDNS) are compared to that of the LES. We present results for the transfer spectra, the skewness and flatness factors of the velocity components, the PDF's of the angle between the vorticity and the eigenvectors of the rate of strain, and that between the vorticity and the vorticity stretching tensor. The above LES statistics are found to be in good agreement with those measured in the fDNS field. We further observe that in all the numerical measurements, the trend was for the LES field to be more gaussian than the fDNS field. Future research on this point is planned.

Shtilman, L.↗
NASA NTRS2015

Effects of Unsteady Flow Interactions on the Performance of a Highly-Loaded Transonic Compressor Stage

The primary focus of this paper is to investigate the loss sources in an advanced GE transonic compressor design with high reaction and high stage loading. This advanced compressor has been investigated both experimentally and analytically in the past. The measured compressor efficiency is significantly lower than the efficiency calculated with various existing tools based on RANS and URANS. The general understanding is that some important flow physics in this modern compressor design are not represented in the current tools. To pinpoint the source of the efficiency miss, an advanced test with detailed flow traverse was performed for the front one and a half stage at the NASA Glenn Research Center. In the present paper, a Large Eddy Simulation (LES) is employed to determine whether a higher-fidelity simulation can pick up any additional flow physics that can explain past efficiency miss with RANS and URANS. The results from the Large Eddy Simulation were compared with the NASA test results and the GE interpretation of the test data. LES calculates lower total pressure and higher total temperature on the pressure side of the stator, resulting in large loss generation on the pressure side of the stator. On the other hand, existing tools based on the RANS and URANS do not calculate this high total temperature and low total pressure on the pressure side of the stator. The calculated loss through the stator from LES seems to match the measured data and the GE data interpretation. Detailed examination of the unsteady flow field from LES indicates that the accumulation of high loss near the pressure side of the stator is due to the interaction of the rotor wake with the stator blade. The strong rotor wake interacts quite differently with the pressure side of the stator than with the suction side of the stator blade. The concave curvature on the pressure side of the stator blade increases the mixing of the rotor wake with the pressure side boundary layer significantly. On the other hand, the convex curvature on the suction side of the stator blade decreases the mixing and the suction side blade boundary layer remains thin. The jet velocity in the rotor wake in the stator frame seems to magnify the curvature effect in addition to inviscid redistribution of wake fluid toward the pressure side of the blade.

LES↗
NASA NTRS2015

RACORO Continental Boundary Layer Cloud Investigations: 1. Case Study Development and Ensemble Large-Scale Forcings

Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.

Boundary layers↗
NASA NTRS2011

Density Weighted FDF Equations for Simulations of Turbulent Reacting Flows

In this report, we briefly revisit the formulation of density weighted filtered density function (DW-FDF) for large eddy simulation (LES) of turbulent reacting flows, which was proposed by Jaberi et al. (Jaberi, F.A., Colucci, P.J., James, S., Givi, P. and Pope, S.B., Filtered mass density function for Large-eddy simulation of turbulent reacting flows, J. Fluid Mech., vol. 401, pp. 85-121, 1999). At first, we proceed the traditional derivation of the DW-FDF equations by using the fine grained probability density function (FG-PDF), then we explore another way of constructing the DW-FDF equations by starting directly from the compressible Navier-Stokes equations. We observe that the terms which are unclosed in the traditional DW-FDF equations are now closed in the newly constructed DW-FDF equations. This significant difference and its practical impact on the computational simulations may deserve further studies.

Shih, Tsan-Hsing↗
NASA NTRS1994

Structure of supersonic jet flow and its radiated sound

The present paper explores the use of large-eddy simulations as a tool for predicting noise from first principles. A high-order numerical scheme is used to perform large-eddy simulations of a supersonic jet flow with emphasis on capturing the time-dependent flow structure representating the sound source. The wavelike nature of this structure under random inflow disturbances is demonstrated. This wavelike structure is then enhanced by taking the inflow disturbances to be purely harmonic. Application of Lighthill's theory to calculate the far-field noise, with the sound source obtained from the calculated time-dependent near field, is demonstrated. Alternative approaches to coupling the near-field sound source to the far-field sound are discussed.

Mankbadi, Reda R.↗
NASA NTRSDate not supplied

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.

Bryson Frank↗
NASA NTRS2013

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

DeBonis, James R.↗
NASA NTRS1998

The NASA-Langley Wake Vortex Modelling Effort in Support of an Operational Aircraft Spacing System

Two numerical modelling efforts, one using a large eddy simulation model and the other a numerical weather prediction model, are underway in support of NASA's Terminal Area Productivity program. The large-eddy simulation model (LES) has a meteorological framework and permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, humidity, and atmospheric turbulence. Results from the numerical simulations are being used to assist in the development of algorithms for an operational wake-vortex aircraft spacing system. A mesoscale weather forecast model is being adapted for providing operational forecast of winds, temperature, and turbulence parameters to be used in the terminal area. This paper describes the goals and modelling approach, as well as achievements obtained to date. Simulation results will be presented from the LES model for both two and three dimensions. The 2-D model is found to be generally valid for studying wake vortex transport, while the 3-D approach is necessary for realistic treatment of decay via interaction of wake vortices and atmospheric boundary layer turbulence. Meteorology is shown to have an important affect on vortex transport and decay. Presented are results showing that wake vortex transport is unaffected by uniform fog or rain, but wake vortex transport can be strongly affected by nonlinear vertical change in the ambient crosswind. Both simulation and observations show that atmospheric vortices decay from the outside with minimal expansion of the core. Vortex decay and the onset three-dimensional instabilities are found to be enhanced by the presence of ambient turbulence.

Proctor, Fred H.↗
DOE OSTI2026

Haze–cloud correlations mediated by supersaturation fluctuations

Atmospheric aerosol particles that contain water-soluble components can absorb water vapor in humid environments and form either haze particles or cloud droplets, depending on supersaturation conditions. Laboratory and in situ measurements have shown that haze particles and cloud droplets often coexist and compete for available water vapor in shallow clouds and fogs, especially under polluted conditions. It is expected that more aerosol particles can form more haze particles and cloud droplets, so that the haze and cloud number concentrations are positively correlated. However, recent large-eddy simulations show that haze and cloud number concentration can be negatively correlated under extremely polluted conditions. Haze–cloud interactions across different environmental settings remain poorly understood. In this study, experiments in a convection cloud chamber with the same aerosol injection rate show that, as supersaturation forcing increases (i.e., changing from polluted to clean conditions), the covariance between haze and cloud droplet number concentration changes from negative to positive and finally to zero. Large-eddy simulations (LES) of the cloud chamber with a fixed supersaturation forcing but varying aerosol injection rates show a similar result for the haze–cloud correlation: near zero in clean, mean-supersaturation-dominated activation conditions; positive in moderate, supersaturation-fluctuation-influenced activation conditions; and negative in polluted, supersaturation-fluctuation-dominated activation conditions. A theoretical covariance framework was developed to interpret this behavior based on the relative magnitude and signs of the correlations between supersaturation and the populations of haze and cloud droplets. Significantly, experiments, LES, and theory all yield the same three-regime behavior for the sign of the haze–cloud covariance. Furthermore, our results show that the haze–cloud covariance remains robust and easily measurable, thereby providing a useful metric for regime identification in the atmosphere, improving regime-aware parameterizations, and informing aerosol interventions such as fog dispersion, rainfall enhancement, and albedo modification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2025

High-fidelity forced convection simulations of the University of Wisconsin–Madison air-cooled reactor cavity cooling system

Among the various types of advanced nuclear technologies that exist, high-temperature gas-cooled reactors (HTGRs) stand out for their co-generation capabilities and exceptional passive safety systems. HTGRs rely on a reactor cavity cooling system (RCCS) to remove decay heat from the reactor pressure vessel (RPV) during the course of long-term transients. Several studies have focused on experimental and numerical modeling of the RCCS; however, the literature reveals a notable gap when it comes to high-fidelity simulations such as direct numerical simulations (DNS) or large eddy simulation (LES) models of the RCCS. An LES model was developed for forced convection inside the RCCS so as to obtain a high-fidelity representation of the RCCS at a scale never before attempted. Furthermore, the LES model was compared against three different Reynolds-averaged Navier–Stokes (RANS) models (Realizable k − ϵ, k − ω SST, and Reynolds stress model [RSM]). This comparison verified the robustness and accuracy of the turbulence modeling approach of the RANS models for the RCCS. Finally, based on the comparison between the LES and the RANS models, the k − ω SST RANS model of the experimental facility was developed and compared against the available experimental data.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1984

Some physical and numerical aspects of computing the effects of viscosity on fluid flow

Some physical, analytical, and computational aspects of viscous flow are examined with reference to examples of computed flows. The discussion of the physical aspects covers the development of important scales used to reference flow phenomena in laminar and turbulent shear layers; the usefulness of the concepts of circulation and vorticity; and some relatively large-scale organized structures that have recently been identified in transitional and fully developed turbulent flows. Among the analytical aspects discussed are a compact presentation of the compressible Navier-Stokes equations, the Reynolds-averaged form of these equations, and a simplified description of some forms of turbulent models. Finally, results of a simulation of the onset of transition, direct turbulence simulations, and large-eddy simulations are given.

Lomax, H.↗