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 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

Here, we use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Application of LES-CFD for predicting pulverized-coal working conditions after installation of NOx control system

The upcoming environmental requirements outlined at the 2017 European Commission Integrated Pollution Prevention and Control (IPPC) are becoming more and more strict in comparison to the existing standards – IED 2010/75/EU. Questions arise about changes in boiler operating conditions after the adoption of the regulations to fulfill specified emission limits. Of initial concern is the question of which technology to pursue in order to reach the specified level of NOx emissions at 150 mg/Nm 3 @6%O 2 . By introducing an air staging technique the NOx limit can be partially achieved; however these changes affect the boiler operation in a detrimental way. In order to investigate impact of de-NOx installation on the temperature and heat flux distribution, numerical techniques can be used. In order to determine the impact of the changes in the boiler operation, here this paper presents the application of an advanced, multiphase open-source code Arches, developed at the University of Utah for modeling pulverized-coal combustion using Large-Eddy Simulation (LES). Numerical simulations showed the effect on the boiler's working parameters e.g., heat-transfer rate, temperature distribution, species concentration, and NOx emission before and after proposed modifications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Gaseous jet through an outward opening injector: Details of mixing characteristic and turbulence scales

Direct injection (DI) strategy of natural gas (NG) into internal combustion engines (ICE) has led to higher thermal efficiency and lower exhaust emissions. In order to thoroughly understand the most relevant phenomena affecting the performances of such engines, computational fluid dynamics (CFD) plays a key role as an accurate description of the jet evolution and interaction within the combustion chamber is required to that aim. Accurate description of high-pressure gaseous jets is rather challenging at high Mach numbers, as the injected gas is strongly under-expanded once in the ambient, giving room to shocks due to compressibility effects. Also the interaction between shock waves and mixing layers needs to be carefully represented with a multi-dimensional model, calling for substantial computational resources requirements. In this paper a numerical investigation of the behavior of a gaseous jet (Argon) through an outward opening injector has been carried out. A Large Eddy Simulation (LES) approach has been used in order to track the structures derived by the interaction of the injected fuel with the surrounding ambient. Although already good results were obtained using a Reynolds Averaged Navier-Stokes (RANS) approach, the adoption of LES is required to characterize more accurately the jet properties in terms of vortex structures and mixing effectiveness. Finally, the effect of the Nozzle Pressure Ratio (NPR) on the jet evolution has been highlighted in the paper, showing how a higher NPR would give a faster injection process, compromising however the homogeneity of the mixture.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Spray process of multi-component gasoline surrogate fuel under ECN Spray G conditions

As modern gasoline direct injection (GDI) engines utilize sophisticated injection strategies, a detailed understanding of the air-fuel mixing process is crucial to further improvements in engine emission and fuel economy. In this study, a comprehensive evaluation of the spray process of single-component iso-octane (IC8) and multi-component gasoline surrogate E00 (36 % n-pentane, 46 % iso-octane, and 18 % n-undecane, by volume) fuels was conducted using an Engine Combustion Network (ECN) Spray G injector. High-speed extinction, schlieren, and microscopy imaging campaigns were carried out under engine-like ambient conditions in a spray vessel. Experimental results including liquid/vapor penetration, local liquid volume fraction, droplet size, and projected liquid film on the nozzle tip were compared under ECN G1 (573 K, 3.5 kg/m 3 ), G2 (333 K, 0.5 kg/m 3 ), and G3 (333 K, 1.01 kg/m 3 ) conditions. In addition to the experiments, preferential evaporation process of the E00 fuel was elucidated by Large–Eddy Simulations (LES). The three-dimensional liquid volume fraction measurement enabled by the computed tomographic reconstruction showed substantial plume collapse for E00 under the G2 and G3 conditions having wider plume growth and plume-to-plume interaction due to the fuel high vapor pressure. Here, the CFD simulation of E00 showed an inhomogeneity in the way fuel components vaporized, with more volatile components carried downstream in the spray after the end of injection. The high vapor pressure of E00 also results in ~4 μm smaller average droplet diameter than IC8, reflecting a higher rate of initial vaporization even though the final boiling point temperature is higher. Consistent with high vapor pressure, E00 had a wider plume cone angle and enhanced interaction with the wall to cover the entire surface of the nozzle tip in a film. However, the liquid fuel underwent faster evaporation, so the final projected tip wetting area was smaller than the IC8 under the flash-boiling condition.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2019

Compressive sensing adaptation for polynomial chaos expansions

Basis adaptation in Homogeneous Chaos spaces rely on a suitable rotation of the underlying Gaussian germ. Several rotations have been proposed in the literature resulting in adaptations with different convergence properties. Here, we present a new adaptation mechanism that builds on compressive sensing algorithms, resulting in a reduced polynomial chaos approximation with optimal sparsity. The developed adaptation algorithm consists of a two-step optimization procedure that computes the optimal coefficients and the input projection matrix of a low dimensional chaos expansion with respect to an optimally rotated basis. Finally, we demonstrate the attractive features of our algorithm through several numerical examples including the application on Large-Eddy Simulation (LES) calculations of turbulent combustion in a HIFiRE scramjet engine.

42 ENGINEERING↗

DOE OSTI2019

Design optimization of a scramjet under uncertainty using probabilistic learning on manifolds

Here, we demonstrate, on a scramjet combustion problem, a constrained probabilistic learning approach that augments physics-based datasets with realizations that adhere to underlying constraints and scatter. The constraints are captured and delineated through diffusion maps, while the scatter is captured and sampled through a projected stochastic differential equation. The objective function and constraints of the optimization problem are then efficiently framed as non-parametric conditional expectations. Different spatial resolutions of a large-eddy simulation filter are used to explore the robustness of the model to the training dataset and to gain insight into the significance of spatial resolution on optimal design.

42 ENGINEERING↗

DOE OSTI2023

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗

DOE OSTI2019

A comparison of two- and three-dimensional single-mode reshocked Richtmyer–Meshkov instability growth

The growth dynamics of two- and three-dimensional single-mode reshocked Richtmyer–Meshkov instability are compared systematically using data from high-resolution implicit large-eddy simulations of a model of the Mach 1.3 air(acetone) and sulfur hexafluoride (Jacobs and Krivets, 2005) shock tube experiment. The vorticity deposition by the incident shock and the dynamics of interface evolution are examined quantitatively and qualitatively. The perturbation amplitudes from the two- and three-dimensional simulations are compared to the experimental data and to the predictions of several nonlinear instability growth models. It is shown that the perturbation amplitudes from the two- and three-dimensional simulations with matching initial Richtmyer velocity are in excellent agreement with the experimental data. In addition, the dynamics of reshock (not considered in the experiment) are described in detail, and the post-reshock mixing layer amplitude growth rate is compared to the predictions of several reshock models. It is shown that using two-dimensional simulations to understand three-dimensional dynamics is valid only at early-to-intermediate times before reshock; at intermediate-to-late times after reshock the three-dimensional growth is generally larger than the corresponding two-dimensional growth. Furthermore, the reshock dynamics are also different between two and three dimensions. The quantitative results, together with visualizations of the flow field, were also used to contrast the difference between two- and three-dimensional vorticity and enstrophy dynamics.

42 ENGINEERING↗

DOE OSTI2018

Reproducing curvature effects due to differential diffusion in tabulated chemistry for premixed flames

Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Modeling stratified flames with and without shear using multiple mapping conditioning

A stochastic sparse particle approach is coupled with an artificial thickening flame (ATF) model for large eddy simulations (LES) to predict a series of turbulent premixed-stratified flames with and without shear and stratification. The thickened reaction progress variable serves as reference variable for the multiple mapping conditioning (MMC) mixing model which emulates turbulent mixing of the stochastic particles. The key feature of MMC is to enforce localness in this reference space when particle pairs are mixed and prevents unphysical mixing of burnt and unburnt fluid across the flame front. We apply MMC-ATF to three flames of a series of turbulent stratified flames and validate the method by comparison with experimental data. The new measurements feature increased accuracy in comparison to previously published data of the same flames due to a better signal-to-noise ratio and a setup which is less prone to beam steering. All flame locations are well predicted by the LES-ATF approach and an analysis of the MMC particle statistics demonstrates that MMC preserves the flamelet-like behaviour in regions where the experiments show low scatter around the flamelet solution. Predicted (local) deviations from the flamelet-solution are comparable to deviations observed in the measurements and variations in the flame structure due to differences in stratification and shear are reasonably well captured by the method.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Effects of detailed geometry and real fluid thermodynamics on Spray G atomization

We present recent results toward the quantification of spray characteristics at engine conditions for an eight-hole counter-bored (stepped) GDI injector-Spray G in the ECN denomination. This computational study is characterized by two novel features: the detailed description of a real injector's internal surfaces via tomographic reconstruction; and a general equation of state that represents the thermodynamic properties of homogeneous liquid-vapor mixtures. The combined level-set moment-of-fluid approach, coupled to an embedded boundary formulation for moving solid walls, makes it possible to seamlessly connect the injector's internal flow to the spray. The Large Eddy Simulation (LES) discussed here presents evidence of partial hydraulic flipping and, during the closing transient, string cavitation. Results are validated by measurements of spray density profiles and droplet size distribution.

42 ENGINEERING↗

DOE OSTI2020

A priori analysis of a power-law sub-filter scale mixing model for transported PDF modeling of high Karlovitz turbulent premixed flames

Accurate modeling of mixing in large-eddy simulation (LES) / transported probability density function (PDF) modeling of turbulent combustion remains an outstanding issue. The issue is particularly salient in turbulent premixed combustion under extreme conditions such as high-Karlovitz number Ka. Here, the present study addresses this issue by conducting an a priori analysis of a power-law scaling based mixing timescale model for the transported PDF model. A recently produced DNS dataset of a high-Ka turbulent jet flame is used for the analysis. A power-law scaling is observed for a scaling factor used to model the sub-filter scale mixing timescale in this high-Ka turbulent premixed DNS flame when the LES filter size is much greater than the characteristic thermal thickness of a laminar premixed flame. The sensitivity of the observed power-law scaling to the different viewpoints (local or global) and to the different scalars for the data analysis is examined and the dependence of the model parameters on the dimensionless numbers Ka and Re (the Reynolds number) is investigated. Different model formulations for the mixing timescale are then constructed and assessed in the DNS flame. The proposed model is found to be able to reproduce the mixing timescale informed by the high-Ka DNS flame significantly better than a previous model.

42 ENGINEERING↗

DOE OSTI2023

A Three-Dimensional, Analytical Wind Turbine Wake Model: Flow Acceleration, Empirical Correlations, and Continuity

A new, three-dimensional, analytical, steady-state wake model is presented that includes local flow acceleration near the rotor, improving the wake description compared to existing models. Wake structures such as the momentum deficit and regions of accelerated flow are concisely described with compound and normal Gaussian functions. Large-eddy simulations (LES) are used as training data to develop the model using two in-line turbines under various inflow conditions parameterized by hub-height wind speed and turbulence intensity. Mass conservation is considered by fixing two components of the wake velocity model and optimizing the third to best satisfy continuity; after which, the model performs comparably if not better than existing work with regards to both relative error and mass consistency. The final model demonstrates a high degree of flexibility, making use of empirical correlations to scale across different inflow conditions. The inclusion of these effects is capable of revealing unused opportunities for enhanced power generation by aligning wake trajectories with these regions of accelerated flow.

modeling↗

DOE OSTI2016

Coherent dynamics in the rotor tip shear layer of utility-scale wind turbines

Here, recent field experiments conducted in the near wake (up to 0.5 rotor diameters downwind of the rotor) of a Clipper Liberty C96 2.5 MW wind turbine using snow-based super-large-scale particle image velocimetry (SLPIV) were successful in visualizing tip vortex cores as areas devoid of snowflakes. The so-visualized snow voids, however, suggested tip vortex cores of complex shape consisting of circular cores with distinct elongated comet-like tails. We employ large-eddy simulation (LES) to elucidate the structure and dynamics of the complex tip vortices identified experimentally. We show that the LES, with inflow conditions representing as closely as possible the state of the flow approaching the turbine when the SLPIV experiments were carried out, reproduce vortex cores in good qualitative agreement with the SLPIV results, essentially capturing all vortex core patterns observed in the field in the tip shear layer. The computed results show that the visualized vortex patterns are formed by the tip vortices and a second set of counter-rotating spiral vortices intertwined with the tip vortices. To probe the dependence of these newly uncovered coherent flow structures on turbine design, size and approach flow conditions, we carry out LES for three additional turbines: (i) the Scaled Wind Farm Technology (SWiFT) turbine developed by Sandia National Laboratories in Lubbock, TX, USA; (ii) the wind turbine developed for the European collaborative MEXICO (Model Experiments in Controlled Conditions) project; and (iii) the model turbine, and the Clipper turbine under varying inflow turbulence conditions. We show that similar counter-rotating vortex structures as those observed for the Clipper turbine are also observed for the SWiFT, MEXICO and model wind turbines. However, the strength of the counter-rotating vortices relative to that of the tip vortices from the model turbine is significantly weaker. We also show that incoming flows with low level turbulence attenuate the elongation of the tip and counter-rotating vortices. Sufficiently high turbulence levels in the incoming flow, on the other hand, tend to break up the coherence of spiral vortices in the near wake. To elucidate the physical mechanism that gives rise to such rich coherent dynamics we examine the stability of the turbine tip shear layer using the theory. We show that for all simulated cases the theory consistently indicates the flow to be unstable exactly in the region where counter-rotating spirals emerge. We thus postulate that centrifugal instability of the rotating turbine tip shear layer is a possible mechanism for explaining the phenomena we have uncovered herein.

17 WIND ENERGY↗

DOE OSTI2018

Subgrid-scale effects in compressible variable-density decaying turbulence

We present that many turbulent flows are characterized by complex scale interactions and vorticity generation caused by compressibility and variable-density effects. In the large-eddy simulation of variable-density flows, these processes manifest themselves as subgrid-scale (SGS) terms that interact with the resolved-scale flow. This paper studies the effect of the variable-density SGS terms and quantifies their relative importance. We consider the SGS terms appearing in the density-weighted Favre-filtered equations and in the unweighted Reynolds-filtered equations. The conventional form of the Reynolds-filtered momentum equation is complicated by a temporal SGS term; therefore, we derive a new form of the Reynolds-filtered governing equations that does not contain this term and has only double-correlation SGS terms. The new form of the filtered equations has terms that represent the SGS mass flux, pressure-gradient acceleration and velocity-dilatation correlation. To evaluate the dynamical significance of the variable-density SGS effects, we carry out direct numerical simulations of compressible decaying turbulence at a turbulent Mach number of 0.3. Two different initial thermodynamic conditions are investigated: homentropic and a thermally inhomogeneous gas with regions of differing densities. The simulated flow fields are explicitly filtered to evaluate the SGS terms. The importance of the variable-density SGS terms is quantified relative to the SGS specific stress, which is the only SGS term active in incompressible constant-density turbulence. It is found that while the variable-density SGS terms in the homentropic case are negligible, they are dynamically significant in the thermally inhomogeneous flows. Investigation of the variable-density SGS terms is therefore important, not only to develop variable-density closures but also to improve the understanding of scale interactions in variable-density flows.

42 ENGINEERING↗

DOE OSTI2018

Similarity of wake meandering for different wind turbine designs for different scales

The wake meandering characteristics of four different wind turbine designs with diameters ranging from a few centimetres (wind tunnel scale) to a hundred metres (utility scale) are investigated using large-eddy simulation with the turbine blades and nacelle parametrised using a new actuator surface model. Different velocity fields and meandering behaviours are observed at near-wake locations. At far-wake locations, on the other hand, the mean velocity deficit profiles begin to collapse when scaled by the centreline velocity deficit based on the incoming wind speed at turbine hub height, suggesting far-wake similarity across scales. The turbine-added turbulence kinetic energy profiles are shown to also nearly collapse with each other in the far wake when normalised using a velocity scale defined by the thrust on the turbine rotor. Moreover, we show that at far-wake locations, the simulated flow fields for all four turbine designs exhibit similar wake meandering characteristics in terms of (1) a Strouhal number independent of rotor designs of different sizes and (2) the distributions of wake meandering wavelengths and amplitudes when normalised by the rotor diameter and a length scale defined by the turbine thrust respectively.

Mechanics↗

DOE OSTI2019

Effect of wind turbine nacelle on turbine wake dynamics in large wind farms

Wake meandering, a phenomenon of large-scale lateral oscillation of the wake, has significant effects on the velocity deficit and turbulence intensities in wind turbine wakes. Previous studies of a single turbine (Kang et al. , J. Fluid. Mech. , vol. 774, 2014, pp. 374–403; Foti et al. , Phys. Rev. Fluids , vol. 1 (4), 2016, 044407) have shown that the turbine nacelle induces large-scale coherent structures in the near field that can have a significant effect on wake meandering. However, whether nacelle-induced coherent structures at the turbine scale impact the emergent turbine wake dynamics at the wind farm scale is still an open question of both fundamental and practical significance. We take on this question by carrying out large-eddy simulation of atmospheric turbulent flow over the Horns Rev wind farm using actuator surface parameterisations of the turbines without and with the turbine nacelle taken into account. While the computed mean turbine power output and the mean velocity field away from the nacelle wake are similar for both cases, considerable differences are found in the turbine power fluctuations and turbulence intensities. Furthermore, wake meandering amplitude and area defined by wake meanders, which indicates the turbine wake unsteadiness, are larger for the simulations with the turbine nacelle. The wake influenced area computed from the velocity deficit profiles, which describes the spanwise extent of the turbine wakes, and the spanwise growth rate, on the other hand, are smaller for some rows in the simulation with the nacelle model. Our work shows that incorporating the nacelle model in wind farm scale simulations is critical for accurate predictions of quantities that affect the wind farm levelised cost of energy, such as the dynamics of wake meandering and the dynamic loads on downwind turbines.

Mechanics↗