Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large Eddy simulations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An intercomparison of wall fluxes in a turbulent thermal convection chamber: Direct numerical simulations and wall-modeled large-eddy simulations enhanced by machine learning

Thermal convection in a closed chamber is driven by a warm bottom, a cold top, and side walls at various temperatures. Although wall fluxes are the source of convection energy, accurately modeling these fluxes (i.e., the wall model) is challenging. In large-eddy simulations (LESs), many wall models are traditionally derived from the canonical boundary layer, which may be unsuitable for thermal convection bounded by both horizontal and vertical walls. This study conducts a model intercomparison of dry convection in a cubic-meter chamber using three direct numerical simulations (DNSs) and four LESs with different wall models. The LESs employ traditional wall models, a new wall model employing physics-aware neural networks, and a refined grid near the walls. The experiment involves four cases with varying sidewall temperatures. Our results show that LESs capture the main flow features and the trends of mean fluxes. The physics-aware neural networks and refined wall grids can improve the temporally averaged local fluxes when the large-scale circulation has a preferred direction. Even without the local improvement of wall fluxes, the LES flow quantities (temperature and velocities) can still largely match those in DNSs, provided the mean flux largely matches the DNSs. Additionally, DNSs reveal that a variation in corner treatments has minimal impacts on the flow quantities away from corners. Finally, LESs underestimate the mean fluxes of the entire wall due to their inability to resolve corner regions, but their mean flux away from the corner can better match DNS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Using observational mean-flow data to drive large-eddy simulations of a diurnal cycle at the SWiFT site

Reproducing realistic date- and site-specific unsteady wind conditions in large-eddy simulations is becoming increasingly useful in wind energy. How to run a large-eddy simulation to match observed conditions, however, remains an open research question. One approach that has received considerable attention is mesoscale-to-microscale coupling, in which information about the mesoscale weather, most commonly acquired from a mesoscale numerical weather model, is passed on to a microscale model. In this paper, we demonstrate how the recently developed profile-assimilation technique, a form of mesoscale-to-microscale coupling, can be used to drive large-eddy simulations solely based on observed mean-flow profiles at a single location, bypassing the need for auxiliary mesoscale simulations. The new approach is evaluated for a diurnal cycle at the Scaled Wind Farm Technology site. Observed mean-flow profiles from the ground up to a height of 2 km are reconstructed by aggregating measurements from multiple instruments, and gaps in the data are infilled with natural neighbor interpolation. We perform nine simulations using various forcing approaches to deal with data limitations. The results show that it is indeed possible to drive microscale large-eddy simulation with observations using the profile-assimilation technique, notwithstanding large gaps in virtual potential temperature measurements. However, profile assimilation with vertical smoothing of the error between the desired and actual profiles is required. Without that smoothing, the microscale simulations develop unrealistically high turbulence levels under many situations. Finally, we show that simulated mesoscale data can account for missing observations, although care is needed as both data sources are not necessarily compatible.

17 WIND ENERGY↗
DOE OSTI2025

A Data Library of Liquid Clouds Modelled With a Large Eddy Simulation Framework

We describe a library of atmospheric large eddy simulations (LES) of liquid-phase boundary layer clouds constructed to enable aerosol–cloud–turbulence interaction studies, support parameterization evaluation and development, and provide training data for machine learning applications. The simulations use a modern LES framework designed for high numerical accuracy, coupled to a detailed spectral bin microphysical scheme. Case studies are configured to represent observed conditions in four key global cloud regions—the Northeastern Atlantic, Northeastern Pacific, Continental United States and Southern Ocean—following a semi-idealised approach. The library also includes aerosol concentration halving and doubling experiments to expose the sensitivities of the case studies to aerosol perturbations. Simulation results are compared to observations on a case-by-case basis, then the library's coverage is evaluated in terms of spreads in meteorological factors and atmospheric boundary layer attributes.

aerosol↗
DOE OSTI2024

EAGLES Liquid Cloud Testbed Large Eddy Simulation Library (v3)

This library consists of large eddy simulation (LES) model output using the PINACLES codebase coupled to the Hebrew University Fast Spectral Bin Microphysics scheme representing shallow, liquid phase clouds from four global liquid cloud testbed regions. This data is particularly suitable for examining interactions among cloud, aerosol, and turbulence in warm, boundary-layer clouds and for developing, assessing and refining parameterizations for coarser-scale models.

Kaul, Colleen M↗
DOE OSTI2021

A Graphics Processing Unit (GPU) Approach to Large Eddy Simulation (LES) for Transport and Contaminant Dispersion

Recent advances in the development of large eddy simulation (LES) atmospheric models with corresponding atmospheric transport and dispersion (AT&D) modeling capabilities have made it possible to simulate short, time-averaged, single realizations of pollutant dispersion at the spatial and temporal resolution necessary for common atmospheric dispersion needs, such as designing air sampling networks, assessing pollutant sensor system performance, and characterizing the impact of airborne materials on human health. The high computational burden required to form an ensemble of single-realization dispersion solutions using an LES and coupled AT&D model has, until recently, limited its use to a few proof-of-concept studies. An example of an LES model that can meet the temporal and spatial resolution and computational requirements of these applications is the joint outdoor-indoor urban large eddy simulation (JOULES). A key enabling element within JOULES is the computationally efficient graphics processing unit (GPU)-based LES, which is on the order of 150 times faster than if the LES contaminant dispersion simulations were executed on a central processing unit (CPU) computing platform. JOULES is capable of resolving the turbulence components at a suitable scale for both open terrain and urban landscapes, e.g., owing to varying environmental conditions and a diverse building topology. In this paper, we describe the JOULES modeling system, prior efforts to validate the accuracy of its meteorological simulations, and current results from an evaluation that uses ensembles of dispersion solutions for unstable, neutral, and stable static stability conditions in an open terrain environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

New insights on wind turbine wakes from large‐eddy simulation: Wake contraction, dual nature, and temperature effects

Abstract Large‐eddy simulation (LES) has been adopted to study wind turbine wakes because it can capture fine‐scale details of turbulent wind flows and interactions with wind turbines. Here, we use the LES version of the Weather Research and Forecasting (WRF) model with an actuator disk model to gain insights on several wake effects that have been traditionally difficult to measure. The first finding is that the wake has a “dual nature,” meaning that the wind speed deficit behaves differently from the added turbulent kinetic energy (TKE) and the two are not co‐located in space. For example, the wind speed deficit peaks at hub height and reaches the ground within 8D (D is the rotor diameter), but added TKE peaks near the rotor tip and generally remains aloft. Second, temperature changes near the ground are driven by the added TKE in the rotor area and by atmospheric stability. The combination of these two factors determines the sign and intensity of the vertical heat flux divergence below the rotor, with convergence and warming associated with stable conditions and weak divergence and modest cooling with unstable conditions. Third, wakes do not expand indefinitely, as suggested by similarity theory applied to the wind speed deficit, but eventually stop expanding and actually contract, at different rates depending on atmospheric stability. The implication of these findings is that, in order to study wakes, it is not sufficient to focus on wind speed deficit alone, because TKE is also important and yet behaves differently from the wind speed deficit.

17 WIND ENERGY↗
DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗
DOE OSTI2022

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Large-Eddy Simulations of Convection Initiation over Heterogeneous, Low Terrain

Abstract Large-eddy simulations are conducted to investigate and physically interpret the impacts of heterogeneous, low terrain on deep-convection initiation (CI). The simulations are based on a case of shallow-to-deep convective transition over the Amazon River basin, and use idealized terrains with varying levels of ruggedness. The terrain is designed by specifying its power-spectral shape in wavenumber space, inverting to physical space assuming random phases for all wave modes, and scaling the terrain to have a peak height of 200 m. For the case in question, these modest terrain fields expedite CI by up to 2–3 h, largely due to the impacts of the terrain on the size of, and subcloud support for, incipient cumuli. Terrain-induced circulations enhance subcloud kinetic energy on the mesoscale, which is realized as wider and longer-lived subcloud circulations. When the updraft branches of these circulations breach the level of free convection, they initiate wider and more persistent cumuli that subsequently undergo less entrainment-induced cloud dilution and detrainment-induced mass loss. As a result, the clouds become more vigorous and penetrate deeper into the troposphere. Larger-scale terrains are more effective than smaller-scale terrains in promoting CI because they induce larger enhancements in both the width and the persistence of subcloud updrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗
DOE OSTI2023

Multifidelity Approach to Sensitivity Estimation in Large-Eddy Simulation

A novel approach to compute affordable approximate sensitivities in a large-eddy simulation (LES) is proposed and assessed. The approach is based on solving a Reynolds-averaged Navier–Stokes (RANS) problem that has been linearized around the mean LES solution, with closure modeling required for the linearized changes in turbulent Reynolds stresses. In the present study, the closure modeling is based on a linearized algebraic turbulence model. The method is assessed for the flow over a NACA0012 airfoil at a fixed angle of attack, with the Reynolds number as the variable parameter. The results show that, provided an accurate linearized closure model, the method predicts the correct sensitivity of the skin friction coefficient and the mean velocity field at the cost of a linearized RANS, which provides an important proof-of-concept for this approach. Here, the linearized algebraic turbulence model with standard model coefficients produces a reasonably accurate sensitivity, but the results also suggest that accuracy could be gained from recalibrating the model coefficients for this new use of the model.

42 ENGINEERING↗
DOE OSTI2021

Large Eddy Simulation of the Southern Ocean

The data set contains seven Large Eddy Simulations (LES) at the Southern Ocean Flux Site for studies of deep turbulent ocean boundary layers, with and without surface wave effects, and with both idealized and observed forcing by wind, surface buoyancy flux and Stokes drift profiles. There are 20 days of hourly statistics computed every half-hour of turbulence quantities; namely, the vertical fluxes of buoyancy (temperature) and momentum, buoyancy and velocity variances, and the turbulent kinetic energy, its production terms and its dissipation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sensitivity Analysis of Wind and Turbulence Predictions With Mesoscale‐Coupled Large Eddy Simulations Using Ensemble Machine Learning

Abstract Coupling between mesoscale models and large‐eddy simulation (LES) models is increasingly used to more realistically represent the wide range of scales of atmospheric motions affecting boundary layer winds and turbulence that need to be simulated accurately for applications such as wind energy. However, such mesoscale‐to‐microscale coupled modeling frameworks are potentially affected by a large number of uncertain closure parameters. Here, we investigate the sensitivity associated with six closure parameters related to a 1.5‐order subgrid‐scale turbulence closure for an ensemble of mesoscale‐coupled LES. The simulations are performed using the Weather Research and Forecasting model nested from horizontal resolutions of greater than a kilometer down to tens of meters. Closure parameters are varied to generate perturbed parameter ensembles for two case studies of highly sheared, convective boundary layers observed in the Columbia Basin of Oregon and Washington during the Second Wind Forecast Improvement Project. Machine learning algorithms are used to explore the sensitivity of LES predictions, considering the effects of the perturbed physical parameters alongside categorical factors such as the case study identity, measurement location, and LES resolution. For the conditions we examine, a single parameter, the eddy viscosity coefficient, is the dominant source of parametric sensitivity and its importance is comparable to the categorical factors for several of the simulation response variables we examine.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coupling turbulent flow with blade aeroelastics and control modules in large-eddy simulation of utility-scale wind turbines

We present a large-eddy simulation framework capable of control co-design of large wind turbines, coupling the turbulent flow environment with blade aeroelastics and turbine controllers. The geometry and aerodynamics of the rotor blades and the turbine nacelle are parameterized using an actuator surface model. The baseline collective pitch control and individual pitch control (IPC) algorithms, consisting of a single-input, single-output proportional–integral controller and two integral controllers, respectively, are incorporated into the simulation framework. Furthermore, a second-order model based on the Euler–Bernoulli beam theory is implemented to describe the blade deformation. Simulations are carried out to investigate the impact of collective and individual pitch control strategies on the deflection of turbine blades. Our results show that the IPC reduces the blade tip deflection fluctuations in the out-of-plane direction, while the fluctuations of the blade tip deflection along the in-plane direction are barely affected by the IPC. Furthermore, the blade out-of-plane deformation fluctuation is underestimated by the one-way coupling approach compared to the two-way coupling approach. Furthermore, the findings of this study reveal the importance of advanced control systems in reducing the dynamic loads on wind turbine blades and underscore the potential of control co-design to reduce the levelized cost of wind energy.

42 ENGINEERING↗
DOE OSTI2026

High-Resolution Large-Eddy Simulations of Historical U.S. Gulf Coast Hurricanes

This dataset provides coupled mesoscale and microscale (large-eddy simulation; LES) atmospheric simulations of five historical U.S. Gulf Coast hurricanes. Meso-microscale coupled simulations of Hurricanes Harvey (2017), Nate (2017), Michael (2018), Laura (2020), and Delta (2020) were performed with the Weather Research and Forecasting (WRF) model v4.1.5. The mesoscale component captures the multi-day evolution of each storm's track, intensity, and large-scale structure across the Gulf of Mexico, while the LES component explicitly resolves the dominant turbulent eddies in the hurricane boundary layer at a horizontal grid spacing of 33.33 m. Together, these simulations characterize tropical cyclone wind fields across atmospheric scales ranging from ~100 km down to ~100 m, spanning storm intensities from Category 2 through Category 4 on the Saffir–Simpson scale. This dataset contains the Large-Eddy simulations. The data will be divided according to the storm intensity level at the simulation time. Category 2: Harvey, Nate Category 3: Michael, Delta Category 4: Michael, Laura

17 WIND ENERGY↗