Search NASA⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRS2014

Intercomparison of Large-Eddy Simulations of Arctic Mixed-Phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds.

clouds↗
NASA NTRS1996

Evaluation of Subgrid-Scale Models for Large Eddy Simulation of Compressible Flows

The objective of this project was to evaluate and develop subgrid-scale (SGS) turbulence models for large eddy simulations (LES) of compressible flows. During the first phase of the project results from LES using the dynamic SGS model were compared to those of direct numerical simulations (DNS) of compressible homogeneous turbulence. The second phase of the project involved implementing the dynamic SGS model in a NASA code for simulating supersonic flow over a flat-plate. The model has been successfully coded and a series of simulations has been completed. One of the major findings of the work is that numerical errors associated with the finite differencing scheme used in the code can overwhelm the SGS model and adversely affect the LES results. Attached to this overview are three submitted papers: 'Evaluation of the Dynamic Model for Simulations of Compressible Decaying Isotropic Turbulence'; 'The effect of the formulation of nonlinear terms on aliasing errors in spectral methods'; and 'Large-Eddy Simulation of a Spatially Evolving Compressible Boundary Layer Flow'.

Blaisdell, Gregory A.↗
DOE OSTI2025

Toward Understanding the Differences between Mesoscale and Large-Eddy Simulations of Tropical Cyclones

In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.

17 WIND ENERGY↗
NASA NTRS2018

Wall-Resolved Large-Eddy Simulations of Transonic Shock-Induced Flow Separation

Wall-resolved large-eddy simulations of shock-induced flow separation over an axisymmet- ric bump for a Mach number of 0.875 and a chord-based Reynolds number of 2.763 million are performed. The incoming boundary layer has a momentum-thickness Reynolds number of 6600 at 1.5 chords upstream of the bump. The calculations, which employ up to 24 billion grid points, simulate the experimental model of Bachalo and Johnson (AIAA Journal, Vol. 24, No. 3, 1986), except that the tunnel walls are ignored and free air is assumed. The effects of domain span and grid resolution are examined along with the main flowfield features. The predicted shock position as well as separation and reattachment locations agree well with the experiment. Grid convergence is observed in the attached region well upstream of separation. Two-point azimuthal correlations suggest that a span of at least 20 degrees is needed in the attached region, while a 120-degree span might be barely large enough in the separated region. Computed root-mean-square surface pressure fluctuations at the shock foot reach about 5 percent of ambient pressure and the level at reattachment location drops to about half of the primary peak. Simulations reveal evidence of the low-frequency shock unsteadiness; however, a much longer statistical sample is needed to investigate this phenomenon.

Uzun, Ali↗
NASA NTRS1995

Preliminary large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids

Large-eddy simulation (LES) has matured to the point where application to complex flows is desirable. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential. In the past two years, important building blocks have been put into place making possible a careful study of LES on unstructured grids. The first building block was an efficient mesh generator which allowed the placement of points according to smooth variation of physical length scales. This variation of length scales is in all three directions independently, which allows a large reduction in points when compared to structured-grid methods, which can only vary length scales in one direction at a time. The second building block was the development of a dynamic model appropriate for unstructured grids. The principle obstacle was the development of an unstructured-grid filtering operator. In the past year, some of the new filters developed by Jansen have been implemented into a highly parallelized finite element code based on the Galerkin/least-squares finite element method. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation for a variety of reasons. First, it is a problem of significant interest since it would be the first LES of an aircraft component. Second, this flow has been the subject of three experimental studies. The third reason for considering this flow is the variety of flow features which provide an important test of the dynamic model. Only the dynamic model can be expected to perform satisfactorily in this variety of situations: from the laminar regions where it must not modify the flow at all to the turbulent boundary layers and wake where it must represent a wide variety of subgrid-scale structures. The flow configuration we have chosen is that of Wadcock (1987) at Reynolds number based on chord Re(sub c) = u(sub infinity)c/v = 1.64 x 10(exp 6), Mach number M = 0.2, and 12 deg angle of attack.

Jansen, Kenneth↗
NASA NTRS1996

Large-Eddy Simulation of Flow Around a NACA 4412 Airfoil Using Unstructured Grids

Large-eddy simulation (LES) has matured to the point where application to complex flows is described. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential.

Jansen, Kenneth↗
NASA NTRS2003

Dynamic Turbulence Modelling in Large-eddy Simulations of the Cloud-topped Atmospheric Boundary Layer

The use of large eddy simulation, or LES, to study the atmospheric boundary layer dates back to the early 1970s when Deardor (1972) used a three-dimensional simulation to determine velocity and temperature scales in the convective boundary layer. In 1974 he applied LES to the problem of mixing layer entrainment (Deardor 1974) and in 1980 to the cloud-topped boundary layer (Deardor 1980b). Since that time the LES approach has been applied to atmospheric boundary layer problems by numerous authors. While LES has been shown to be relatively robust for simple cases such as a clear, convective boundary layer (Mason 1989), simulation of the cloud-topped boundary layer has proved more of a challenge. The combination of small length scales and anisotropic turbulence coupled with cloud microphysics and radiation effects places a heavy burden on the turbulence model, especially in the cloud-top region. Consequently, over the past few decades considerable effort has been devoted to developing turbulence models that are better able to parameterize these processes. Much of this work has involved taking parameterizations developed for neutral boundary layers and deriving corrections to account for buoyancy effects associated with the background stratification and local buoyancy sources due to radiative and latent heat transfer within the cloud (see Lilly 1962; Deardor 1980a; Mason 1989; MacVean & Mason 1990, for example). In this paper we hope to contribute to this effort by presenting a number of turbulence models in which the model coefficients are calculated dynamically during the simulation rather than being prescribed a priori.

Kirkpatrick, M. P.↗
NASA NTRS1993

Large-eddy simulation of laminar-turbulent breakdown at high speeds with dynamic subgrid-scale modeling

The laminar-turbulent breakdown of a boundary-layer flow along a hollow cylinder at Mach 4.5 is investigated with large-eddy simulation. The subgrid scales are modeled dynamically, where the model coefficients are determined from the local resolved field. The behavior of the dynamic-model coefficients is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both formulations proposed by Germano et al. and Lilly are used for the determination of unique coefficients for the dynamic model and their results are compared and assessed. The behavior and the energy cascade of the subgrid-scale field structure are investigated at various stages of the transition process. The investigations are able to duplicate a high-speed transition phenomenon observed in experiments and explained only recently by the direct numerical simulations of Pruett and Zang, which is the appearance of 'rope-like' waves. The nonlinear evolution and breakdown of the laminar boundary layer and the structure of the flow field during the transition process were also investigated.

El-Hady, Nabil M.↗
NASA NTRS2021

Large-Eddy Simulations of Marine Boundary-Layer Clouds Associated With Cold Air Outbreaks During the Activate Campaign – Part 1: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary-layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LES simulations. We derive these quantities from measurements and reanalysis obtained for two cold air outbreak (CAO) events during Phase I of the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) in February-March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from the reanalysis data ERA5 agree reasonablywell with those derived fromACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time evolving heat fluxes, wind and advective tendencies profiles from ERA5 reanalysis data to drive the LES.We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 reanalysis data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LES simulations. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path imulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol-cloud-meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

Large-eddy simulation↗
NASA NTRS2011

Large Eddy Simulation Study for Fluid Disintegration and Mixing

A new modeling approach is based on the concept of large eddy simulation (LES) within which the large scales are computed and the small scales are modeled. The new approach is expected to retain the fidelity of the physics while also being computationally efficient. Typically, only models for the small-scale fluxes of momentum, species, and enthalpy are used to reintroduce in the simulation the physics lost because the computation only resolves the large scales. These models are called subgrid (SGS) models because they operate at a scale smaller than the LES grid. In a previous study of thermodynamically supercritical fluid disintegration and mixing, additional small-scale terms, one in the momentum and one in the energy conservation equations, were identified as requiring modeling. These additional terms were due to the tight coupling between dynamics and real-gas thermodynamics. It was inferred that if these terms would not be modeled, the high density-gradient magnitude regions, experimentally identified as a characteristic feature of these flows, would not be accurately predicted without the additional term in the momentum equation; these high density-gradient magnitude regions were experimentally shown to redistribute turbulence in the flow. And it was also inferred that without the additional term in the energy equation, the heat flux magnitude could not be accurately predicted; the heat flux to the wall of combustion devices is a crucial quantity that determined necessary wall material properties. The present work involves situations where only the term in the momentum equation is important. Without this additional term in the momentum equation, neither the SGS-flux constant-coefficient Smagorinsky model nor the SGS-flux constant-coefficient Gradient model could reproduce in LES the pressure field or the high density-gradient magnitude regions; the SGS-flux constant- coefficient Scale-Similarity model was the most successful in this endeavor although not totally satisfactory. With a model for the additional term in the momentum equation, the predictions of the constant-coefficient Smagorinsky and constant-coefficient Scale-Similarity models were improved to a certain extent; however, most of the improvement was obtained for the Gradient model. The previously derived model and a newly developed model for the additional term in the momentum equation were both tested, with the new model proving even more successful than the previous model at reproducing the high density-gradient magnitude regions. Several dynamic SGS-flux models, in which the SGS-flux model coefficient is computed as part of the simulation, were tested in conjunction with the new model for this additional term in the momentum equation. The most successful dynamic model was a "mixed" model combining the Smagorinsky and Gradient models. This work is directly applicable to simulations of gas turbine engines (aeronautics) and rocket engines (astronautics).

Bellan, Josette↗
NASA NTRSDate not supplied

Improved Prediction of Losses with Large Eddy Simulation in a Low-Pressure Turbine

There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of freestream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Reex=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.

CFD↗
NASA NTRS1988

Model consistency in large eddy simulation of turbulent channel flows

Combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are examined by a priori tests and numerical simulations. The structure of the subgrid scales is found to depend strongly on the type of filter used, and consistency between model and filter is essential to ensure accurate results. The implementation of consistent combinations of filter and model gives more accurate turbulence statistics than those obtained in previous investigations in which the models were chosen independently from the filter. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.

Piomelli, Ugo↗
NASA NTRSDate not supplied

Wall-Modeled Large-Eddy Simulations of Jet Noise in Flight Conditions

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics (LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction (PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology (CST) project. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several operating points both for static (no ambient co-flow), and in-flight (M ͚ =0.3 co-flow) conditions. The operating conditions range in jet exit Mach number from 0.38 to 1.1 with nozzle temperature ratios (NTR) from 0.84 to 2.7. The effect of the flight-stream on far-field noise sound spectra is assessed. Comparison of LES predictions to microphone array measurements demonstrate good agreement within the resolved frequency range. A dip in the predicted low frequency noise spectra for observers between 120° and 145° is observed. This dip is smaller for lower Mach numbers and seems to be correlated to the Mach wave radiation angle. The Mach 1.1 jet shows broadband-shock associated noise. While the onset of BBSN appears to be captured correctly in WMLES, some differences in its magnitude and the prominent frequency at which it occurs persist between the experiment and the simulations. The effect of the outer nozzle boundary layer state created by the co-flow is assessed and shows to be important for accurate comparisons with experiments. A change of 2.5dB between a slip-wall condition and a artificially thickened turbulent boundary layer was observed. In addition, simulations were performed with an alternative nozzle geometry that includes a internal plug and has twice the nozzle exit diameter. These two configurations resulted in very comparable spectra which is consistent with experimental observations. A generally stronger deviation from experimental results is observed for in-flight cases compared to static conditions. The applicability and correct usage of acoustic analogies used for far-field propagation with strong turbulent co-flows needs to be investigated more systematically using canonical problems to improve comparisons between experiments and WMLES. This is especially true for coherent noise sources seen in BBSN.

CST↗
NASA NTRSDate not supplied

Large Eddy Simulations of Upper-Ocean Turbulence Across the Equatorial Pacific Cold Tongue

Turbulence in the equatorial Pacific cold tongue plays multiple roles in the sea-surface temperature budget. However, observations of turbulence are limited to a handful of field experiments lasting from a few days to one month at 0 N , 140 W, and the key turbulent fluxes can only be inferred. Although explicit simulations of turbulence resolve these fluxes, the required large eddy simulations (LES) are costly and thus also limited to a few weeks. A two-fold hypothesis motivates this contribution to the PUMP pre-field modelling experiments in support of TPOS: 1) that there exist useful relationships between the surface wind, wave, and buoyancy forcing, oceanic shear and stratification, and turbulent mixing, but 2) these relationships differ between regions of equatorial dynamics (f~0) and higher latitudes, such that the state space is insufficiently sampled to understand these relationships and represent them in ocean general circulation models (OGCM). To address the lack of observations, we simulated 35 continuous days of upper ocean turbulence using LES during October-November 1985 at sites in two regions (70 days in total) of the cold tongue: 0 N, 140 W and 3 N, 140 W. The model reproduces many features of the observed turbulence at 0,140W, including the “deep-cycle” observed by Gregg et al. (1985) Nature and the turbulent pulses every few hours during the night and in the early morning. Furthermore, the results at 3 N support the hypothesis advanced by Cherian et al. (2021) JPO that strong, but intermittent deep-cycle turbulence is observed to the north of the equator in the presence of strong shear associated with Tropical Instability Waves. Comparisons with the K-profile parameterization (KPP) scheme show that additional tuning of the KPP parameters can reduce differences between the turbulent fluxes produced by KPP and LES at both sites. The results will suggest modifications to OGCMs and already emphasize the need for a more extensive measurements of vertical profiles of shear, stratification and turbulent mixing to reveal the spatio-temporal structure of turbulence and how it relates to forcing, shear and stratification across the equatorial Pacific Cold tongue.

Large Eddy Simulations↗
DOE OSTI2025

Numerical Investigation of Ammonia/n-Heptane Dual-Fuel Spray Flames Using Large Eddy Simulations

Diesel engines are extensively used in heavy-duty transportation, power generation, and marine vehicles due to their superior thermal efficiency and extended high-load operability compared to spark ignition (SI) engines. However, combustion in diesel engines is generally characterized by locally rich fuel–air mixtures and high combustion temperatures, causing significant amounts of soot and NO x emissions from these engines. Utilizing carbon-free alternative fuels and enhancing fuel efficiency represent promising strategies to mitigate greenhouse gas (GHG) and other emissions in the heavy-duty transportation sector. In this context, ammonia (NH 3 ), as a hydrogen carrier, has received significant attention as a viable substitute for hydrocarbon fuels due to its carbon-free composition, relatively high energy density, and well-established infrastructure. Many previous studies have considered combustion and emission characteristics of ammonia-hydrocarbon fuel blends in engines and simplified flames. But, detailed investigations on the effects of ammonia on the performance of hydrocarbon fuels under engine conditions are lacking. In the present study, we perform large eddy simulations (LES) of the ignition and flame processes in a constant-volume combustion reactor, where n-heptane is injected in an ammonia/air ambient mixture in a diesel-like environment. A detailed and validated reaction mechanism containing 302 species and 1981 reactions is employed. The Engine Combustion Network Spray H experimental data is used to validate the spray model under both non-reacting and reacting conditions. Dual-fuel combustion is simulated using the well-stirred reactor (WSR) approach. Results are presented for two spray cases: (1) single fuel (SF) with n-heptane injected into a mixture of air and combustion products and (2) dual-fuel (DF) with the injection of n-heptane in a mixture of air, ammonia, and combustion products. It is observed that the presence of ammonia has a significant effect on the ignition and flame development processes. With ammonia addition, both the first- and second-stage ignition delay times increase, but the effect of ammonia on the second-stage ignition is significantly more prominent. In addition, the ignition kernel size and growth rate decrease noticeably. For SF spray, the main ignition is characterized by multiple ignition kernels near the spray tip, whereas for DF spray, a single relatively small ignition kernel forms and grows slowly in the downstream direction. The flame development and the final quasi-steady flame structure are also modified due to ammonia. Here, the outcome of this research would enable a better understanding of ammonia–diesel dual-fuel spray flame behavior and guide the development of associated engine combustion strategies.

ammonia↗
DOE OSTI2026

Differences in cluster and internal wake effects from mesoscale and large-eddy simulations off the US East Coast

Mesoscale simulations are increasingly used to estimate wake effects within and between large wind farms, despite limited validation for large-scale wake effects. This study evaluates the capabilities and limitations of mesoscale simulations in capturing wake-induced impacts on wind turbine power production through a direct comparison with large-domain large-eddy simulations (LESs) for three planned offshore wind farms under realistic atmospheric conditions and a range of atmospheric stabilities. We assess mesoscale performance in replicating wake characteristics behind single and multiple turbine clusters and quantify the resulting variability in mean turbine power. Results show that mesoscale Weather Research and Forecasting simulations with the Fitch wind farm parameterization capture key features of the velocity deficit downstream of both single and multiple wind farms, with mean root-mean-square errors near 5 % and good agreement with stability-driven wake behavior. However, in these simulations, the mesoscale Fitch parameterization underestimates power losses from internal wake effects, particularly when turbines align with the prevailing wind direction or under stable stratification. In these conditions, individual wakes persist and dominate downstream power deficits. The coarse resolution of the mesoscale simulations limits their ability to resolve individual wind turbine wakes that drive power fluctuations within wind farms. Nonetheless, mesoscale simulations can yield accurate estimates of combined wake losses from internal and cluster effects across some wind direction sectors, where errors in wake representation may cancel each other out. These findings underscore the strengths of mesoscale simulations for capturing broader wake patterns while highlighting their limitations for modeling turbine-level power losses. Future work should explore hybrid modeling approaches to capture both long-range cluster wake propagation and localized internal wake dynamics.

17 WIND ENERGY↗
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 ↗