Search NASA⌕ Search

SEARCH · Search NASA

Results for “large-eddy simulation (LES)”

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 OSTI2026

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗
DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗
DOE OSTI2025

Comparison of wind-farm control strategies under realistic offshore wind conditions: wake quantities of interest

Wind-farm control strategies aim to increase the efficiency, and therefore lower the levelized cost of energy, of wind farms. This is done by using turbine settings such as the yaw angle, blade pitch angles, or generator torque to manipulate the wake that negatively affects downstream turbines in the farm. Two inherently different wind-farm control methods have been identified in the literature: wake steering (WS) and active wake mixing (AWM). As one of two companion papers focused on understanding practical aspects of these two wind-farm control strategies using large-eddy simulation (LES), we below analyze the wake quantities of interest for a single wind turbine performing WS and AWM, while the companion article (Frederik et al., 2025) focuses on turbine quantities of interest including power and structural loads for the same computational setup and also includes two-turbine arrays with full and partial wake overlap. The simulations, which are based in the LES solver AMR-Wind, are tailored to have inflow conditions representative of measurements from a site off the East Coast of the US, including with strong veer and low turbulence. The turbine, which is modeled in OpenFAST and coupled to the LES, is the IEA 15 MW, an open-source offshore design. After presenting an overview of the wake recovery for the different wake-control cases, the analysis probes the fluid-dynamic causes for the different performance of the arrays reported in the companion article by examining control volumes around the wakes and the budget of the mean-flow kinetic energy (MKE) within these volumes. In the high veer environment considered, the MKE recovery is dominated by mean convection, and this is shown to especially benefit the WS strategy when a neighboring turbine is directly downstream: there is ≈65 % more available power for a downstream turbine than in the baseline case, and this power is gained primarily through mean convection on the left-tip and top-tip faces of the control volume. However, the case with imperfect knowledge of the exact wind direction favors the pulse-type AWM strategy, largely because of ≈9 % increased turbulent entrainment from aloft versus the baseline that could be related to an apparent resistance to skewing in the pulsed wake. The general reduced effectiveness of helix-type and other individual-pitch-based AWM strategies for inflow with high veer and low turbulence as reported in the companion paper is due, in part, to low magnitudes of phase-averaged turbulent entrainment. Two main findings of this study are thus that veer has a significant impact on the effectiveness of different wake-control strategies and that pulse-type AWM may be a useful strategy when the objective is power maximization in realistic, offshore flow environments with imperfect knowledge of the exact wake overlap position on the downstream turbine.

17 WIND ENERGY↗
DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗
DOE OSTI2026

Dependence of wind-farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration

This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .

17 WIND ENERGY↗
DOE OSTI2026

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗
DOE OSTI2026

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗
DOE OSTI2026

Power output of turbines mounted on tension-leg platforms subjected to fully developed ocean gravity waves

A concern in the deployment of large wind turbines on ocean floating platforms is the effect of floating-platform motions on their electrical power generation. Further, it is not clear how floating motions influence waking, which might affect the combined power generation of collections of turbines. We examine the average power output of a single and a collection of NREL 5 MW wind turbines mounted on a tension-leg platform (TLP) under the action of fully developed ocean wave motions, coupling floating motions with the large-eddy simulation (LES) of atmospheric and rotor dynamics. The ocean dynamics enter as fully developed waves derived from the Pierson–Moskowitz spectrum. To assess the influence of ocean motions, we performed simulations over the full range of wind speeds in the operational range of the turbine, reporting comparisons of average power output when the platforms are allowed to move to when they are held rigidly in place. In all simulations, we find that the effects of the TLP floating-platform's induced motions have a minor effect on single and multiple turbine power production and wake deficits. Even when using coherent and large amplitude harmonic-floating-induced perturbations, any significant wake modifications from floating motions are confined to the near-wake region, where downstream turbines are unlikely to be located. The relatively small amplitude of TLP motions relative to pre-existing turbulent fluctuations are the primary reason for low wake and power modifications downstream.

Restrepo, Juan [ORNL] (ORCID:0000000326092882)↗
DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

LER25 (Large Eddy Reconstruction 2025) [SWR-25-65]

Large-eddy simulation codes struggle to reconstruct the state of a real-world atmosphere. To date, LES codes have been able to reconstruct time-averaged statistics from field campaigns in the real world. LER25 presents an approach to zoom in and reconstruct the second-by-second dynamics of a real-world atmosphere.

Rybchuk, Alex [National Renewable Energy Laborator↗
DOE OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗
DOE OSTI2026

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗
NASA NTRSDate not supplied

The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗
NASA NTRSDate not supplied

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗
DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗
DOE OSTI2025

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗
DOE OSTI2025

Resolving Low Cloud Feedbacks Globally With E3SM High‐Res MMF: Agreement With LES but Stronger Shortwave Effects

This study investigates low cloud feedback in a warmer climate using global simulations from the High-Resolution Multi-scale Modeling Framework (HR-MMF), which explicitly simulates small-scale eddies globally. Two 5-year simulations—one with present-day sea surface temperatures (SSTs) and a second with SSTs warmed uniformly by 4 K—reveal a positive global shortwave cloud radiative effect (SWCRE = 0.3 W/m 2 /K), comparable to estimates from CMIP models. As the climate warms, significant reductions in low cloud cover occur over stratocumulus regions. This study is the first attempt to compare HR-MMF results with predictions from idealized large-eddy simulations from the CGILS intercomparison. Despite different underlying assumptions, we find qualitative agreement in SWCRE and inversion height changes between HR-MMF and CGILS predictions. This suggests reasonable credibility for the CGILS framework in predicting cloud responses under the out-of-sample conditions found in HR-MMF. However, the HR-MMF exhibits stronger SWCRE changes than predicted by CGILS. We explore potential causes for this discrepancy, examining variations in cloud-controlling factors (CCFs) and cloud conditions. Our results show a fairly homogeneous SWCRE response, with little systematic variation tied to the variations in CCFs. This reveals a dominant role for SST forcing in modulating SWCRE.

boundary-layer clouds↗