Results for “large Eddy simulations”

DOE OSTI2022

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗

NASA NTRS1982

Turbulent flow simulation - A large eddy simulator's viewpoint

The state of the art of the data base on turbulent structures and numerical turbulence modelling is assessed. The failure of statistical methods to produce accurate results revealed that turbulence is more than simply random noise. Both randomness and coherence may be characteristics of turbulence, which may also feature noise that is a remnant of previous coherence. Five methods are presently used for numerically predicting turbulent flow behavior: correlations, integral methods, the k-epsilon model, and phenomenological and exact solutions to the Navier-Stokes equations. The higher level models include large eddy simulation and full simulation with numerical solutions to the Navier-Stokes equations; the two methods are applicable to experimental, rather than design work.

Ferziger, J. H.↗

NASA NTRS2010

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Shih, Tsan-Hsing↗

DOE OSTI2024

A Stereo Camera Simulator for Large-Eddy Simulations of Continental Shallow Cumulus Clouds Based on Three-Dimensional Path-Tracing

The complex spatial and temporal structure of cumulus clouds complicates their representation in weather and climate models. Classic meteorological instrumentation struggles to fully capture these features. Networks of multiple high-resolution hemispheric cameras are increasingly used to fill this data gap, and provide information on this missing multi-dimensional spatial information. In this study, a path-tracing algorithm is used to generate virtual camera images of resolved clouds in large-eddy simulations (LES). These images are then used as a camera network simulator, allowing reconstructions of three-dimensional cloud edges from the model output. Because the actual LES cloud field is fully known, the combined path-tracing and reconstruction method can be statistically analyzed. The method is applied to LES realizations of summertime shallow cumulus at the Jülich Observatory for Cloud Evolution (JOYCE), Germany, which also routinely operates a camera network. We find that the path-tracing method allows accurate reconstruction of up to 70% of the visible cloud edges. Additional sensitivity tests show that the method is robust for changes in its hyperparameters. The sensitivity to cloud optical thickness is also investigated, finding a cloud boundary placement error of approximately 182 m. This error can be considered typical for cloud boundary reconstruction using real stereo camera imagery. The results provide proof of principle for future use of the method for evaluating LES clouds against camera network imagery, and for further optimizing the configuration of such camera networks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Capturing differential diffusion effects in large eddy simulation of turbulent premixed flames

Not provided.

Energy & Fuels↗

DOE OSTI2024

Large-eddy simulations of the NACA23012 airfoil with laser-scanned ice shapes

Not provided.

Engineering↗

DOE OSTI2025

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗

DOE OSTI2025

Collaborative Proposal: Improving understanding of the internal structure and dynamics of deep convection using ARM observations and large eddy simulations

Recent observational and large eddy simulation (LES) modeling studies have nearly unanimously supported the view of deep cumulus convection being composed of a series of quasi-spherical bubbles of buoyant air, known as moist thermals. Despite the prevalence of moist thermals in deep convection, a comprehensive theory for the dynamics of these structures is lacking. Most current conceptual models for cumulus convection are based on canonical scaling theories for dry thermals or plumes; however, there is considerable evidence that the behavior of moist thermals differs markedly from these theories. Furthermore, the theoretical basis for most cumulus parameterizations originates from the plume conceptual model, and therefore these parameterizations are inconsistent with the real structure of moist convection. Motivated by the aforementioned knowledge gaps, this “end-to-end” research effort use theory, observations, numerical simulations, and direct improvements to the Zhang-McFarlane (ZM) convection scheme in the global climate Community Atmosphere Model (CAM) to address the following research questions: What key environmental parameters determine whether or not shallow convection will transition into deep convection, in the context of thermal-like updrafts? What factors regulate the size of thermals within cumulus updrafts? How does vertical wind shear influence thermal behavior, and as a consequence, vertical velocity and mass flux profiles and the shallow-to-deep convective transition? What are the critical processes that determine updraft vertical velocities and their connection to the vertical mass flux profile for thermal-like updrafts? Idealized LES modeling will be used in conjunction with theoretical models for the core properties of thermal-like updrafts to better understand key processes that regulate thermal ascent rates and entrainment properties. Thermal-tracking procedures will be used to characterize the behavior of thermals within the LES, and recently developed direct measures of entrainment and detrainment will be used to quantify entrainment/detrainment rates. Building from these results, we will analyze the structure of moist thermals from hemispheric range-height indicator scans taken during the Atmospheric Radiation Measurement Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, and from “real case” LES of CACTI events. This combined modeling and observational analysis will provide essential validation for the existing body of research on moist thermal dynamics, which is based primarily on modeling studies. With the insight gained from the aforementioned activities, we will modify the Zhang-McFarlane convection scheme to improve its representation of updraft vertical velocity and entrainment rate profiles. These process-level changes will be tested in the Community Atmosphere Model to assess the impact on global climate simulations.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Transonic Lift and Drag Predictions using Wall Modelled Large Eddy Simulations

Wall-modelled Large Eddy Simulations (WMLES) of the NASA Common Research Model (CRM) at transonic conditions and various angles of attacks leading up to and including shock-induced flow separation are performed using the LAVA computational framework. The simulations are shown to accurately predict the lift curve slope and the onset of separation characterized by the break in the pitching moment. Furthermore, careful assessment of skin-friction drag at cruise condition and its subsequent decrease with increasing angles of attack are shown to be in agreement with viscous sublayer resolving Reynolds Averaged Navier Stokes(RANS) simulations. The small differences between WMLES and RANS appear to be of the same order as differences seen between two RANS models at the cruise-point. Some sensitivity is observed to the coefficient used in the subgrid scale model, although this can be reconciled by noting the low chord- and shock-incidence Reynolds numbers (Rex≈106) seen in the outboard regions of the wing, along with the uncertainties associated with tripping and the numerical transition that occurs near the leading edge. Among the primary drawbacks seen in the WMLES predictions is the subdued increase in wave-drag with increasing angle-of-attack and the shock intensity when compared with experimental data. Preliminary assessment suggests that additional span- and stream-wise grid resolution is likely needed in mid- and out-board portions of the wing to better resolve the shock-induced separated flow and to further investigate prediction accuracy of unsteady temporal characteristics of the problem.

ARMD↗

DOE OSTI2021

Large-Scale Forcing Impact on the Development of Shallow Convective Clouds Revealed From LASSO Large-Eddy Simulations

Real-world large-eddy simulations (LES) are driven by time-varying large-scale forcings (LSF) - e.g., temperature advection, moisture advection, and subsidence - derived from large-scale weather models. This study investigates the impact of the uncertainty in LSF on real-world LES in terms of the development of shallow convection at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory for the 11 June 2016 case using LES provided by the U.S. Department of Energy's LES ARM Symbiotic Simulation and Observation (LASSO) activity. The LASSO dataset provides an ensemble of LES for the selected case, which consists of LES runs that were driven by different LSF. The two contrasting LES runs investigated here generate different types of convective clouds, i.e., nonprecipitating shallow clouds and precipitating cumulus congestus, mainly due to the difference of LSF in temperature advection in the free troposphere. The temperature advection modulates the strength of the capping inversion and therefore the buoyancy of the air parcels rising from the atmospheric boundary layer (ABL). The inversion, together with large-scale updrafts, controls the penetration of the ABL thermals into the free troposphere, leading to cumulus congestus in the case of a weaker inversion. In contrast, clouds remain shallow in the case of a strong inversion. Differences between the two simulations are amplified over time, as mixed-phase clouds are formed near the top of the congestus in the weaker inversion case. Furthermore, this high dependency of LES results to LSF stresses the importance of accurate LSF by large-scale models to real-world LES simulations.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Analysis of subgrid models using direct and large-eddy simulations of isotropic turbulence

Direct and large eddy simulations of forced and decaying isotropic turbulence have been performed using a pseudospectral and a finite-difference code. Subgrid models that include a one-equation subgrid kinetic energy model with and without a stochastic backscatter forcing term and a new scale similarity model have been analyzed in both Fourier space and physical space. The Fourier space analysis showed that the energy transfer across the cutoff wavenumber k(sub c) is dominated by local interaction. The correlation between the exact and the modeled (by a spectral eddy viscosity) nonlinear terms and the subgrid energy transfer in physical space was found to be quite low. In physical space, a similar correlation analysis was carried out using top hat filtering. Results show that the subgrid stress and the energy flux predicted by the subgrid models correlates very well with the exact data. The scale similarity model showed very high correlation for reasonable grid resolution. However, with decrease in grid resolution, the scale similarity model became more uncorrelated, when compared to the kinetic energy subgrid model. The subgrid models were then used for large-eddy simulations for a range of Reynolds number. It was determined that the dissipation was modeled poorly and that the correlation with the exact results was quite low for all the models. In general, for coarse grid resolution, the scale similarity model consistently showed very low correlation while the kinetic energy model showed a relatively higher correlation. These results suggest that to use the scale similarity model relatively fine grid resolution may be required, whereas, the kinetic energy model could be used even in coarse grid.

Menon, S.↗

NASA NTRS2001

Statistical Ensemble of Large Eddy Simulations

A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Carati, Daniele↗

NASA NTRS2011

Investigating the Impact of Surface Heterogeneity on the Convective Boundary Layer Over Urban Areas Through Coupled Large-Eddy Simulation and Remote Sensing

Large-eddy Simulation (LES) was used to study convective boundary layer (CBL) flow through suburban regions with both large and small scale heterogeneities in surface temperature. Constant remotely sensed surface temperatures were applied at the surface boundary at resolutions of 10 m, 90 m, 200 m, and 1 km. Increasing the surface resolution from 1 km to 200 m had the most significant impact on the mean and turbulent flow characteristics as the larger scale heterogeneities became resolved. While previous studies concluded that scales of heterogeneity much smaller than the CBL inversion height have little impact on the CBL characteristics, we found that further increasing the surface resolution (resolving smaller scale heterogeneities) results in an increase in mean surface heat flux, thermal blending height, and potential temperature profile. The results of this study will help to better inform sub-grid parameterization for meso-scale meteorological models. The simulation tool developed through this study (combining LES and high resolution remotely sensed surface conditions) is a significant step towards future studies on the micro-scale meteorology in urban areas.

Dominguez, Anthony↗

DOE OSTI2023

Department of Energy’s Atmospheric System Research (ASR) Program’s Workshop on the Future of Atmospheric Large Eddy Simulation (LES): Workshop Report

Large-eddy simulation (LES) is used as a tool to understand physical processes such as turbulence, aerosols, clouds, precipitation, radiation, the interactions among all these, and their interactions with the underlying surface. Over the next 10 years, LES will drive fundamental progress in open scientific questions in these areas as LES is increasingly used to gain understanding of complex interacting physical processes involving atmospheric turbulence. This growth will be driven both by scientific demand and the expansion of computational resources needed to conduct LES, and the form that the growth takes will largely be determined by how computational resources are leveraged for scientific gain. In particular, we suggest that computational resources are likely to be leveraged in two separate but not necessarily distinct ways. On one hand, growth in computational resources will allow LES to be made more routine, that is, performed more frequently, while on the other hand, the computational expense (measured in total floating point operations) afforded to individual LES will expand dramatically, allowing simulations to increase in both domain size and resolution as well as physical detail. Current U.S. Department of Energy (DOE) projects such as LES ARM Symbiotic Simulation and Observation Activity (LASSO) are leading the way in conducting routine LES, building large, public databases that are accessible for data science, sensitivity studies, and training for machine learning. LES will also become more routine as it becomes more accessible for individual researchers to address their scientific questions of interest. Scientific questions addressed by LES over the next 10 years are likely to include cloud organization and aggregation; aerosol cloud interactions and atmospheric chemistry (including geo-engineering); urban-scale LES; atmospheric extreme events, ranging from small-scale severe weather to wildfires; and ocean-wave-atmosphere interactions. Further LES-related research will likely grow significantly in areas related to societal impact studies of air quality and extreme weather events, applications to renewable energy forecasting and resource assessment, and aid in decision-making processes. The growth in the use of LES in atmospheric science research will drive the need for better physical process representations (e.g., cloud aerosol microphysics, radiation, and atmospheric chemistry) at the scales resolved by LES. To date, many of the process representations used by LES have been taken directly from coarser-resolution models. Promising methods for LES process representations include superdroplet and quadrature methods for microphysics, 3D approaches for radiation, and better representation of chemistry and aerosol processes. At LES resolution, land-atmosphere interactions for complex terrains, land cover/types, biogeochemistry, and plant canopy models are needed as an improvement beyond traditional and widely used Monin-Obuhkov similarity theory.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Department of Energy’s Atmospheric System Research (ASR) Program’s Workshop on the Future of Atmospheric Large Eddy Simulation (LES) (Workshop Report)

Large-eddy simulation (LES) is used as a tool to understand physical processes such as turbulence, aerosols, clouds, precipitation, radiation, the interactions among all these, and their interactions with the underlying surface. Over the next 10 years, LES will drive fundamental progress in open scientific questions in these areas as LES is increasingly used to gain understanding of complex interacting physical processes involving atmospheric turbulence. This growth will be driven both by scientific demand and the expansion of computational resources needed to conduct LES, and the form that the growth takes will largely be determined by how computational resources are leveraged for scientific gain. In particular, we suggest that computational resources are likely to be leveraged in two separate but not necessarily distinct ways. On one hand, growth in computational resources will allow LES to be made more routine, that is, performed more frequently, while on the other hand, the computational expense (measured in total floating point operations) afforded to individual LES will expand dramatically, allowing simulations to increase in both domain size and resolution as well as physical detail.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Department of Energy’s Atmospheric System Research (ASR) Program’s Workshop on the Future of Atmospheric Large Eddy Simulation (LES) (Workshop Report)

Large-eddy simulation (LES) is used as a tool to understand physical processes such as turbulence, aerosols, clouds, precipitation, radiation, the interactions among all these, and their interactions with the underlying surface. Over the next 10 years, LES will drive fundamental progress in open scientific questions in these areas as LES is increasingly used to gain understanding of complex interacting physical processes involving atmospheric turbulence. This growth will be driven both by scientific demand and the expansion of computational resources needed to conduct LES, and the form that the growth takes will largely be determined by how computational resources are leveraged for scientific gain. In particular, we suggest that computational resources are likely to be leveraged in two separate but not necessarily distinct ways. On one hand, growth in computational resources will allow LES to be made more routine, that is, performed more frequently, while on the other hand, the computational expense (measured in total floating point operations) afforded to individual LES will expand dramatically, allowing simulations to increase in both domain size and resolution as well as physical detail. Current U.S. Department of Energy (DOE) projects such as LES ARM Symbiotic Simulation and Observation Activity (LASSO) are leading the way in conducting routine LES, building large, public databases that are accessible for data science, sensitivity studies, and training for machine learning. LES will also become more routine as it becomes more accessible for individual researchers to address their scientific questions of interest. Scientific questions addressed by LES over the next 10 years are likely to include cloud organization and aggregation; aerosol cloud interactions and atmospheric chemistry (including geo-engineering); urban-scale LES; atmospheric extreme events, ranging from small-scale severe weather to wildfires; and ocean-wave-atmosphere interactions. Further LES-related research will likely grow significantly in areas related to societal impact studies of air quality and extreme weather events, applications to renewable energy forecasting and resource assessment, and aid in decision-making processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Adaptive Determination of the Optimal Exchange Location in Wall-Modeled Large-Eddy Simulation

Wall-modeled large-eddy simulation introduces a modeling interface (or exchange location) separating the wall-modeled layer from the rest of the domain. The current state-of-the-art is to rely on user expertise when choosing where to place this modeling interface, whether this choice is tied to the grid or not. This paper presents a postprocessing algorithm that determines the exchange location systematically. The algorithm is based on a model for the error in the predicted wall shear stress and a model for the computational cost, and then finds the exchange location that minimizes a combination of the two. Here, the algorithm is tested both a priori and a posteriori using an equilibrium wall model for the flow over a wall-mounted hump, a boundary layer in an adverse pressure gradient, and a shock/boundary-layer interaction. The algorithm produces exchange locations that mostly agree with what an experienced user would suggest, with thinner wall-modeled layers in nonequilibrium flow regions and thicker wall-modeled layers where the boundary layer is closer to equilibrium. This suggests that the algorithm should be useful in simulations of realistic and highly complex geometries.

42 ENGINEERING↗