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 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the acceleration and separation regions of the flow. An unstructured finite-volume solver is used along with an equilibrium wall model. The Mach number of the oncoming flow selected for the simulation is 0.176, and two Reynolds numbers (Re) are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Spanwise periodic results for the lower Re case are compared with the available DNS data, while the full three-dimensional simulation results for the higher Re case are compared to available experimental data. Sensitivity of WMLES results are assessed for a number of factors including grid resolution, wall model exchange location, type of wall model, and unstructured grid topology. The WMLES results for the lower Re case agree well with available DNS data in terms of the wall pressure variation, velocity and turbulent stress profile comparisons. The skin-friction predictions show a reduced tendency to separate. This appears to be accentuated for the higher Re case, for the medium resolution grid used in this study, where WMLES does not predict any flow separation in the mean, which is in contrast to the large separation bubble observed in the experiments.

turbulence↗

DOE OSTI2018

Shallow Cumulus (SHALLOWCUMULUS) Value-Added Product Report

Periods of shallow cumulus clouds are identified automatically at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory in order to help determine times of interest for Large-Eddy Simulation (LES) and Observation (LASSO) modeling efforts. This value-added data product (VAP) uses input from the Cloud Type (CLDTYPE) VAP (with cloud types determined from Active Remotely Sensed Cloud Locations [ARSCL] VAP cloud boundaries), the ceilometer, and the total sky imager (TSI). The data set uses an hourly temporal resolution and only classifies events as shallow cumulus when low clouds are detected for at least 1.5 hours and meet specified criteria for cloud occurrence and cloud fraction. In addition to time periods that contain just shallow cumulus, cases with overlying cirrus and a number of transitional cases (transitions to/from cirrus, stratus, and altocumulus/altostratus) are also identified.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Numerical investigation of wind turbine wakes under high thrust coefficient

Abstract We study wind turbine wakes of rotors operating at high thrust coefficients ( C T > 24/25) using large‐eddy simulations with a rotating actuator disk model. Wind turbine wakes at high thrust coefficients are different from wakes at low thrust coefficients. Wakes behave differently at high thrust, with increased turbulence and faster recovery. Lower induction in the wake is achieved because wakes in high‐thrust conditions recover much faster than in normal operating conditions. This enhanced recovery is possible thanks to the turbulence generated in the near wake. We explore the mechanism behind this behavior and propose a simple model to reproduce it. We also propose a Gaussian fit for the wakes under high‐thrust conditions and use it use it to initialize an Ainslie type model within the FAST.Farm framework.

FAST.Farm↗

DOE OSTI2023

Surface‐Atmosphere Decoupling Prolongs Cloud Lifetime Under Warm Advection Due To Reduced Entrainment Drying

Abstract An initially well‐mixed stratocumulus deck can remain overcast for several tens of hours under warm‐advection conditions, although moisture supply is cut off from the ocean due to surface‐atmosphere decoupling (stabilization of the surface‐atmosphere interface). In this study, a set of idealized large‐eddy simulations were performed to investigate the physical mechanism of how warm‐air advection impacts the evolution of a pre‐existing stratocumulus deck. To mimic warm‐air advection, we decrease the sea surface temperature linearly over time in a doubly periodic domain. Given the same initial conditions, the stratocumulus deck is more persistent when experiencing warm‐air advection than cold‐air advection. This persistence is caused by reduced cloud‐top entrainment drying due to decoupling, a process more influential than the decoupling‐induced cutoff of moisture supply. This mechanism is more notable when the free troposphere becomes more humid. The relevance of the mechanism to previous observations of less low‐level cloudiness under warm‐advection conditions is discussed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Toward control co-design of utility-scale wind turbines: Collective vs. individual blade pitch control

A large-eddy simulation framework has been coupled with controller modules to systematically investigate the impacts of collective (CPC) and individual (IPC) pitch control strategies on utility-scale wind turbine energy production and fatigue loads. Wind turbine components were parameterized using an actuator surface model to simulate the rotor blades and the turbine nacelle. The baseline CPC and IPC algorithms, consisting of single-input single-output proportional–integral controllers and two integral controllers, respectively, were incorporated into the numerical framework. A series of simulations were carried out to investigate the relative performance of the two controllers under various turbulent inflow conditions, spanning hub-height velocities of 7 to 14 m/s. The numerical simulation results of this study showed that, in comparison to the CPC, the IPC controller could successfully reduce the damage equivalent loads of utility-scale turbines at regions 2 and 3 of turbine operation by about 3% and 40%, respectively, without any penalty on the power production of the turbine. It was also shown that, despite its minor impact on the turbulence kinetic energy of the wake, the IPC controller did not influence the recovery of the turbine wake.

17 WIND ENERGY↗

DOE OSTI2019

Parametric investigation of the transition to turbulence in Rayleigh-Taylor mixing

In this study, a database is developed consisting of highly resolved large-eddy simulations of the transition process in a Rayleigh-Taylor mixing layer. Over 300 simulations are performed with parametric variations in Grashof number, Schmidt number, Atwood number, and initial perturbation spectrum. The database is then analyzed to explore how the factors that influence the transition-to-turbulence process impact mixing evolution and growth rate prior to the development of late-time, steady-state growth. The competition between entrainment-driven mixing and diffusion-driven mixing is observed, and two different mechanisms are suggested by which the homogeneity of mixing might increase, rather than decrease, at early time. The impact of parametric variation on growth rate is additionally examined, and all variations in initial conditions are found to converge to approximately the same growth parameter around α ≈ 0.025.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

An analysis of the potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) indicates that satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with WMLES computations of subharmonic transition in a flat-plate boundary layer that yield reasonable agreement with DNS results. However, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation↗

NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

Wall-modeled large-eddy simulation↗

DOE OSTI2024

Three-step LES-C models for flows at high Reynolds numbers

Here, we investigate the need for the second correction step in the recently proposed LES-C models for fluid flows at high Reynolds numbers. These models use a predictor-corrector idea to enhance the efficiency of the existing Large Eddy Simulation models. Different three-step(one defect step, two corrections) LES-C models, based on the Leray-α, ADM and NS-ω LES models, are tested in three different situations. The new Leray-α-C2 model (C2 stands for two correction steps) is applied to the Navier–Stokes equations; the ADC2 is applied to the MagnetoHydroDynamic flow; and the NS-ω-C2 is used in the fluid-fluid interaction problem. We evaluate the effectiveness of the second correction step in all these settings, using qualitative and quantitative numerical tests.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Mass and momentum transport in the Tilted Rocket Rig experiment

Here, high-fidelity large eddy simulations of the inclined Rayleigh-Taylor instability in the Tilted Rocket Rig experimental configuration [Smeeton and Youngs, AWE Report No. O 35/87 (1987)] are performed using a tenth-order compact finite difference code. These simulations are analyzed for spatial distributions of turbulent kinetic energy, turbulent mass flux velocity, species mass fraction flux, species mass fraction variance, and Favre-averaged Reynolds stresses at two time instances, t = 37 ms and t = 52 ms. Additionally, the vertical distribution of the components of the unclosed budget equations over the center of the domain for these quantities are examined. The dominant terms of these budget equations are further decomposed to examine the principal contributions to these terms in each axis. Notably, the principal contribution to the horizontal turbulent mass flux velocity budget is found to be from a term which is commonly neglected in many Reynolds-averaged Navier-Stokes models.

42 ENGINEERING↗

DOE OSTI2022

Computation of conventional and alternative jet fuel sensitivity to lean blowout

Large Eddy Simulations (LES) are performed to compute the sensitivity of a conventional (A-2) and an alternate bio-jet (C-1) fuel to Lean Blowout (LBO). A realistic aviation gas turbine engine combustor configuration is considered. Reliable experimental LBO data and OH* chemiluminescence data for the conventional and alternate jet fuel in the combustor configuration have recently become available. The present work utilizes a highly automated, on-the-fly meshing strategy, along with adaptive mesh refinement, to demonstrate the feasibility of capturing the realistic combustion processes. A Lagrangian framework, with initial conditions specified using measurements of spray statistics, is used to model the fuel spray. Newly developed compact reaction mechanisms based on fuel surrogates are validated for the A-2 and the C-1 fuels. The compact reaction mechanisms are implemented using a detailed finite rate chemistry solver. Spray statistics computed by the present LES simulations compare well with available measurements at stable flame conditions near the lean blowout limit. The computed shape of the stable flame as represented by line integrated OH concentrations compares well with the experimental OH* chemiluminescence data. Lean blowout is reached by gradually decreasing the fuel flow rate in the computations, similar to that in the experiments. The results of the LES simulations effectively capture the fuel composition effects and estimate the sensitivity of the LBO limits to the fuel type. The computed trends in LBO limits agree within engineering accuracy with the experimental results for conventional and alternative aviation fuels. The methodology for predicting the fuel composition effects on the lean blowout limits in a fully resolved realistic, complex combustor is established for the first time.

Adaptive mesh refinement↗

DOE OSTI2020

Deep learning-based model for progress variable dissipation rate in turbulent premixed flames

A deep neural network (DNN) based large eddy simulation (LES) model for progress variable dissipation rate in turbulent premixed flames is presented. The DNN model is trained using filtered data from direct numerical simulations (DNS) of statistically planar turbulent premixed flames with n-heptane as fuel. Training data was comprised of flames with varying turbulence levels leading to a range of Karlovitz numbers. Through a-priori tests the DNN model is shown to predict the subfilter contribution to progress variable dissipation rate accurately over a range of filter widths and for all Karlovitz numbers examined in this study. Superior performance of the DNN model relative to an established physics-based model is also demonstrated. Additionally, transferability of the DNN model is highlighted by a-priori evaluation of the model using filtered DNS data from multiple cases with different Karlovitz numbers and fuel species than those that were used for training the model.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2019

Chemical explosive mode analysis of lean blowout in a gas turbine combustor

In this study, stabilization and lean blowout (LBO) mechanisms in a gas turbine combustor are numerically studied based on large eddy simulation (LES) data and the chemical explosive mode analysis (CEMA). Three CEMA-based local combustion modes, that is, auto-ignition, diffusion-assisted ignition, and extinction, are employed to quantify interactions between local diffusion and chemical processes in complex turbulent fields. The capability of CEMA-based local modes in extinction detection is first demonstrated in a one-dimensional non-premixed flamelet undergoing transient extinction. Flame structure and dynamics in the gas turbine combustor are then examined in detail using CEMA and the local combustion modes for both the stably burning and lean blowout conditions. Local structure comprised of all the three modes are identified, and each mode is found to play an important role in determining the overall burning rate. The heat release rate attributed to each local mode remains nearly constant at the stably burning condition. In contrast, the flame during LBO shows a two-stage extinction behavior: at the initial stage, the extinction mode is largely enhanced and becomes dominant over the assisted-ignition model; at the second stage, all three modes are suppressed, signifying global extinction. It is therefore suggested that the excessive local extinction mode, due to an increased amount of heat and radical losses through diffusion, finally leads to lean flame blowout.

chemical explosive mode↗

DOE OSTI2019

Parameterization of Vertical Profiles of Governing Microphysical Parameters of Shallow Cumulus Cloud Ensembles Using LES with Bin Microphysics

Shallow convection is a subgrid process in cloud-resolving models for which their grid box is larger than the size of small cumulus clouds (Cu). At the same time such Cu substantially affect radiation properties and thermodynamic parameters of the low atmosphere. The main microphysical parameters used for calculation of radiative properties of Cu in cloud-resolving models are liquid water content (LWC), effective droplet radius, and cloud fraction (CF). In this study, these parameters of fields of small, warm Cu are calculated using large-eddy simulations (LESs) performed using the System for Atmospheric Modeling (SAM) with spectral bin microphysics. Despite the complexity of microphysical processes, several fundamental properties of Cu were found. First, despite the high variability of LWC and droplet concentration within clouds and between different clouds, the volume mean and effective radii per specific level vary only slightly. Second, the values of effective radius are close to those forming during adiabatic ascent of air parcels from cloud base. These findings allow for characterization of a cloud field by specific vertical profiles of effective radius and of mean liquid water content, which can be calculated using the theoretical profile of adiabatic liquid water content and the droplet concentration at cloud base. Using the results of these LESs, a simple parameterization of cloud-field-averaged vertical profiles of effective radius and of liquid water content is proposed for different aerosol and thermodynamic conditions. These profiles can be used for calculation of radiation properties of Cu fields in large-scale models. The role of adiabatic processes in the formation of microstructure of Cu is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Reconciling Chord Length Distributions and Area Distributions for Fields of Fractal Cumulus Clouds

While the total cover of broken cloud fields can in principle be obtained from one-dimensional measurements, the cloud size distribution normally differs between two-dimensional (area) and one-dimensional retrieval (chord length) methods. In this study, we use output from high-resolution Large Eddy Simulations to generate a transfer function between the two. We retrieve chord lengths and areas for many clouds, and plot the one as a function of the other, and vice versa. We find that the cloud area distribution conditional on the chord length behaves like a gamma distribution with well-behaved parameters, with a mean μ=1.1L and a shape parameter β=L−0.645. Using this information, we are able to generate a transfer function that can adjust the chord length distribution so that it comes much closer to the cloud area distribution. Our transfer function improves the error in predicting the mean cloud size, and is performs without strong biases for smaller sample sizes. However, we find that the method is still has difficulties in accurately predicting the frequency of occurrence of the largest cloud sizes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Process-oriented stochastic perturbations applied to the parametrization of turbulence and shallow convection for ensemble prediction

Weather forecasting nowadays often requires some estimation of uncertainties associated with the output of meteorological models, in order to better inform decision making, especially in the context of intense weather events. Ensemble prediction systems provide such information through sets of possible scenarios which are designed to represent various uncertainty sources, including model uncertainties. A wide variety of methods have been proposed to estimate model uncertainties, among which perturbation methods targeting uncertain processes are a promising research field. Here, we focus on the representation of small-scale variability by process-oriented perturbation schemes applied to two key physical processes, namely turbulence and shallow convection. The perturbations are applied to a single-column version of the convection-permitting AROME model, in three idealized boundary-layer cases. Large-eddy simulations (LESs) of the same cases serve as a reference for the subgrid variability that has to be represented, and the results are also compared to those given by the Stochastically Perturbed Parametrization Tendencies (SPPT) method, which is a method commonly used by weather forecast centres to represent model uncertainty. The spread produced by our process-oriented perturbations of turbulence and shallow convection does not represent all the small-scale variability implied by the LESs for temperature and humidity. However, it is of a similar order of magnitude for the wind, thanks to perturbations generated by the stochastic turbulence scheme. The dispersion is structurally different from what is obtained with SPPT. It is non-negligible in the lower levels, where SPPT perturbations are usually suppressed because of numerical instabilities, indicating a possible complementarity between the schemes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Description of the ARM Large-Scale Forcing Data from the Constrained Variational Analysis (VARANAL) Version 2

This technical report represents an update of the previous technical report written by Zhang et al., which described the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility constrained variational analysis (VARANAL) that is used to develop the large-scale forcing data for driving single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LES). The VARANAL algorithm was originally developed by Zhang and Lin and Zhang et al. at Stony Brook University and was migrated to the Lawrence Livermore National Laboratory (LLNL) as the ARM operational objective analysis system in May 1999. Since then, the algorithm has been evolved with time along with the availability of new observations and techniques to meet various modeling needs. Major updates include: 1. The method used to develop multi-year continuous forcing data, 2. The incorporation of eddy correlation flux measurement system (ECOR) turbulent fluxes into the analysis, and 3. Improvements to the workflow (e.g., implementing part of the code into the ARM Data Integrator [ADI]) to increase efficiency. The ARM large-scale forcing data have been widely used for SCM/CRM/LES to understand and improve physical processes in models. Zhang et al. has provided a comprehensive review of the SCM concept, early efforts to derive forcing data for SCM studies, efforts of the ARM constrained variational analysis, and previous SCM studies using ARM cases. This technical report focuses on the constrained variational analysis algorithm and the introduction of the ARM VARANAL products. We also extended the VARANAL algorithm into a three-dimensional constrained variational analysis (3DCVA) and designed an ensemble framework to address the forcing uncertainty. The 3D large-scale forcing data are released as another datastream named “varanal3d”. Please refer to the ARM varanal3d technical report for more information.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗