Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak turbulence”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coherent structures in stably stratified wall-bounded turbulent flows

To date, a growing body of literature has documented the existence and impacts of coherent structures known as large- and very-large-scale motions within wall-bounded turbulent flows under neutral and unstable thermal stratification. These coherent structures can account for a considerable fraction of the overall turbulent transport and have been found to modulate small-scale turbulent fluctuations near the wall. In the context of stably stratified flows, however, the examination of such coherent structures has garnered relatively little attention. Stable stratification limits vertical transport and turbulent mixing within flows, which makes it unclear the extent to which previous findings on coherent structures under unstable and neutral stratification are applicable to stably stratified flows. In this study, we investigate the existence and characteristics of coherent structures under stable stratification with a wide range of statistical and spectral analyses. Outer peaks in premultiplied spectrograms under weak stability indicate the presence of large-scale motions, but these peaks become weaker and eventually vanish with increasing stability. Quadrant analysis of turbulent transport efficiencies (the ratio of net fluxes to their respective downgradient components) demonstrates dependencies on both stability and height above ground, which is evidence of morphological differences in the coherent structures under increasing stability. Amplitude modulation by large-scale streamwise velocity was found to decrease with increasing gradient Richardson number, whereas modulation by large-scale vertical velocity was approximately zero across all stability ranges. For sufficiently stable stratification, large eddies are suppressed enough to limit any inner–outer scale interactions.

Greene, Brian R. (ORCID:0000000343766818)↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2022

Pressure–strain interaction. II. Decomposition in magnetic field-aligned coordinates

In weakly collisional and collisionless magnetized plasmas, the pressure–strain interaction describes the rate of conversion between bulk flow and thermal energy density. In this study, we derive an analytical expression for the pressure–strain interaction in a coordinate system with an axis aligned with the local magnetic field. The result is eight groups of terms corresponding to different physical mechanisms that can contribute to the pressure–strain interaction. We provide a physical description of each term. Here, the results are immediately of interest to weakly collisional and collisionless magnetized plasmas and the fundamental processes that happen therein, including magnetic reconnection, magnetized plasma turbulence, and collisionless shocks. The terms in the field-aligned coordinate decomposition are likely accessible to measurement with satellite observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic shear effect on plasma transport at T e / T i ~ 1 through electron cyclotron heating in DIII-D plasmas

The effect of magnetic shear on plasma transport for an electron to ion temperature ratio ( T e / T i ) near unity has been explored in DIII-D utilizing electron cyclotron heating (ECH). Previous reports showed that significant confinement degradation occurred at T e / T i ~ 1 in positive shear (PS) plasmas in DIII-D, whereas reduced confinement degradation was observed in negative central shear (NCS) plasmas. In this study, plasma transport in weak magnetic shear (WS) plasmas with ECH is investigated and compared with that in NCS and PS plasmas. Here the magnetic shears ($\hat s$) are $\hat s$ > 0.5, ~0 and <-0.1 in the core region (ρ~ 0.3–0.4) of PS, WS and NCS plasmas, respectively, and flat or negative inside ρ~ 0.4 in the WS and NCS plasmas. Weak magnetic shear is found to be effective in minimizing degradation of ion thermal confinement as T e / T i increases through ECH application, and an improved confinement factor of H 98y2 ~ 1.2 is maintained, similar to NCS plasmas. At T e / T i ~ 1, the ion thermal diffusivity around an internal transport barrier decreases when changing the magnetic shear from positive to weak or negative shear. Also, reduced local particle and momentum transport was indicated by steeper density and toroidal rotation profiles in the weak and negative shear regimes. Linear gyrokinetic simulations predict little change in growth rates of low-k turbulence with ECH application in the WS and NCS plasmas, which is consistent with the transport and profile analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Tethered balloon-borne profile measurements of atmospheric properties in the cloudy atmospheric boundary layer over the Arctic sea ice during MOSAiC: Overview and first results

The tethered balloon-borne measurement system BELUGA (Balloon-bornE moduLar Utility for profilinG the lower Atmosphere) was deployed over the Arctic sea ice for 4 weeks in summer 2020 as part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. Using BELUGA, vertical profiles of dynamic, thermodynamic, aerosol particle, cloud, radiation, and turbulence properties were measured from the ground up to a height of 1,500 m. BELUGA was operated during an anomalously warm period with frequent liquid water clouds and variable sea ice conditions. Three case studies of liquid water phase, single-layer clouds observed on 3 days (July 13, 23, and 24, 2020) are discussed to show the potential of the collected data set to comprehensively investigate cloud properties determining cloud evolution in the inner Arctic over sea ice. Simulated back-trajectories show that the observed clouds have evolved within 3 different air masses (“aged Arctic,” “advected over sea ice,” and “advected over open ocean”), which left distinct fingerprints in the cloud properties. Strong cloud top radiative cooling rates agree with simulated results of previous studies. The weak warming at cloud base is mostly driven by the vertical temperature profile between the surface and cloud base. In-cloud turbulence induced by the cloud top cooling was similar in strength compared to former studies. From the extent of the mixing layer, it is speculated that the overall cloud cooling is stronger and thus faster in the warm oceanic air mass. Larger aerosol particle number concentrations and larger sizes were observed in the air mass advected over the sea ice and in the air mass advected over the open ocean.

Arctic↗
DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: CCN Distribution

Abstract Following our previous investigation of the turbulence impact on cloud-base single-size CCN activation, this study considers a similar problem assuming CCN size distribution obtained from field measurements. The total CNN concentration is taken as either 200 cm −3 to represent clean conditions, or as 2000 cm −3 to represent polluted conditions. CCN is assumed to be sodium chloride. The CCN activation in the rising nonturbulent adiabatic parcel is contrasted with the activation within a rising adiabatic parcel filled with inertial-range homogeneous isotropic turbulence. The turbulent parcel of 64 3 m 3 and the turbulent kinetic energy dissipation rate of 10 −3 m −2 s −3 are used in most of the simulations. Results for a range of mean parcel ascent rates, between 0.125 and 8 m s −1 , are discussed. Overall, the adiabatic turbulent parcel simulations show results consistent with the adiabatic nonturbulent parcel, with higher activated CCN concentrations for stronger parcel ascent rates. The key difference is a blurriness of the separation between dry CCN size bins featuring activated and nonactivated (haze) CCN, especially for weak mean ascent rates. The blurriness comes from CCN getting activated and subsequently deactivated in the fluctuating supersaturation field, instead of all becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the nonturbulent parcel when activation is completed. Modeling results are discussed in the context of the impact of turbulent fluctuations on CCN activation documented in laboratory experiments using the Pi chamber.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Stability of superthermal strahl electrons in the solar wind

We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗
DOE OSTI2023

Edge radiated temperature fluctuations across confinement regime transitions in favorable and unfavorable drift configurations at ASDEX Upgrade

Abstract The physics behind the difference in heating power threshold for the low to high confinement (L–H) transition between favorable ion B × ∇ B drift (toward active X -point) and unfavorable ion B × ∇ B drift (away from active X -point) magnetic configurations is a longstanding open question. In this work we show differences in edge ( ρ p o l = 0.95 –1.0) turbulence in matched heating power favorable and unfavorable magnetic configuration discharges using a correlation electron cyclotron emission diagnostic to measure radiated temperature ( T r a d ) fluctuations. During power matched L-mode phases, the edge T r a d fluctuations have higher amplitude in the unfavorable configuration plasma as compared to the favorable configuration plasma. The lower fluctuation amplitude in favorable configuration is observed concomitant with increased E × B shear, as compared to the unfavorable case. Although this turbulence is different in amplitude, it is similar in radial and frequency structure between favorable and unfavorable configurations and grows to the weakly coherent mode during the I-mode phase of the unfavorable configuration discharge. A broadband pedestal fluctuation is seen during an inter-ELM H-mode period of the favorable configuration discharge.

Physics↗
DOE OSTI2023

Insights of Boundary Layer Turbulence Over the Complex Terrain of Central Himalaya from GVAX Field Campaign

Limited observations hinder understanding of turbulent characteristics in mountainous terrain resulting from heating or cooling of slopes, wind, vertical motions, and heat or moisture advection, which disperse aerosols and other pollutants over the region. In this study, the 1290 MHz radar wind profiler data are utilized to compute the boundary layer height (BLH), the refractive index structure constant (C n 2 ), and the energy dissipation rate (ε) over the central Himalayan site for the period of November 2011 to March 2012, from the intense Ganges Valley Aerosol Experiment (GVAX) field measurements. The radar wind profiler (RWP) based estimation of BLH and ε is validated against the radiosonde, representing the effectiveness of the datasets for further investigation. The strong seasonal variation of log C n 2 and log ε, with average values of ≈ -12 m -2/3 and -2 m 2 s -3 , respectively, is associated with the mountain-induced local circulations and stability in the atmospheric boundary layer. The weak stratification during weak flow is found to be responsible for deep mixing, particularly in the nocturnal boundary layer in spring. Furthermore, the level of cloud cover significantly impacts the strength of turbulence, with the highest cloud cover resulting in a substantial increase in log C n 2 (approximately -11 m -2/3 ) due to intense updraft and downdraft motions compared to clear skies. Additionally, the distribution of aerosol loading across the site, coupled with the behavior of BLH, atmospheric stability, and orographic-induced circulations, implies distinctive seasonal mechanisms for transporting aerosols toward the mountains. This study offers valuable insights into the diurnal and seasonal patterns of turbulent mixing and the mechanisms behind the transport of pollutants through boundary layer processes over the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2026

On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf

Type Ia supernovae are bright thermonuclear explosions that are important to numerous areas of astronomy. However, the origins of these events are poorly understood. One proposed setting is that of a near Chandrasekhar mass white dwarf that undergoes runaway carbon burning in the core. During the thousand years leading up to the explosion, the white dwarf undergoes a simmering phase where slow carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, which links convection with weak nuclear reactions. We use the low Mach number code MAESTROeX to perform full 3D simulations as is required to accurately capture the turbulent convection. We present simulations with and without the A=23 convective Urca process, which have relaxed to a steady state. We characterize the effects of the convective Urca process on the neutrino losses, the nuclear energy generation, and the convective boundary. We find that the size of the convection zone is substantially reduced by the convective Urca process, though convection still extends past the Urca shell. Our findings on the structure of the convective zone and the compositional changes can be used to inform 1D stellar models that track the longer-timescale evolution.

FOS: Physical sciences↗
DOE OSTI2023

Electrostatic turbulence in EAST plasmas with internal transport barrier

Based on first-principles nonlinear gyrokinetic simulations, the electrostatic turbulence properties in the internal transport barrier (ITB) region of an Experimental Advanced Superconducting Tokamak discharge (#93890) are investigated. Specifically, ITBs with steep density and temperature gradients are located in the weakly negative magnetic shear region at the plasma center. In the linear stage, the growth rate and frequency of the ion temperature gradient (ITG) mode increase significantly due to resonant excitation by trapped electrons. That is, the resonance between trapped electrons and the ITG becomes strong due to the precession drift reversal of trapped electrons by the negative magnetic shear and Shafranov shift. Meanwhile, the trapped electron mode is stable in the ITB region due to only a very small fraction of electrons precessing in the direction of the electron diamagnetic drift. Nonlinear simulations show that, after considering the non-adiabatic effect of trapped electrons, the heat conductivity of ions and the turbulence intensity increase by at least a factor of 7 compared with the results only considering the adiabatic effect of electrons. The zonal charge density of trapped electrons can partially cancel that of ions, which weakens the intensity of the zonal flow, and consequently reduces the zonal flow regulation and enhances the turbulent transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Are Supercells Resistant to Entrainment because of Their Rotation?

This research investigates a hypothesis posed by previous authors, which argues that the helical nature of the flow in supercell updrafts makes them more resistant to entrainment than nonsupercellular updrafts because of the suppressed turbulence in purely helical flows. It was further supposed that this entrainment resistance contributes to the steadiness and longevity of supercell updrafts. A series of idealized large-eddy simulations were run to address this idea, wherein the deep-layer shear and hodograph shape were varied, resulting in supercells in the strongly sheared runs, nonsupercells in the weakly sheared runs, and variations in the percentage of streamwise vorticity in updrafts among runs. Fourier energy spectrum analyses show well-developed inertial subranges in all simulations, which suggests that the percentages of streamwise and crosswise vorticity have little effect on turbulence in convective environments. Additional analyses find little evidence of updraft-scale centrifugally stable flow within updrafts, which has also been hypothesized to limit horizontal mass flux across supercell updrafts. Results suggest that supercells do have smaller fractional entrainment rates than nonsupercells, but these differences are consistent with theoretical dependencies of entrainment on updraft width, and with supercells being wider than nonsupercells. Thus, while supercells do experience reduced fractional entrainment rates and entrainment-driven dilution, this advantage is primarily attributable to increased supercell updraft width relative to ordinary convection, and has little to do with updraft helicity and rotation.

Peters, John M.↗
DOE OSTI2020

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Benchmarking of three DWM-based wake models at below-rated wind speeds

Wind turbine wake models are essential tools for predicting power losses and structural loads in wind farms. Among these, the dynamic wake meandering (DWM) model, included as a recommended approach in the International Electrotechnical Commission design standard, is a widely used engineering-fidelity method that balances accuracy and computational cost. This study compares the performance of three DWM-based wake model implementations (from the Technical University of Denmark, the National Renewable Energy Laboratory, and the Institute for Energy Technology) under below-rated wind speed conditions. Model predictions of wake flow, power output, and structural loads for a four-turbine row are evaluated across different ambient turbulence levels and wind-direction misalignments and compared against high-fidelity large-eddy simulation results. All three models captured the overall wake evolution and mean turbine performance with reasonable accuracy; their predicted time-averaged thrust and power were typically within 5 %–10 % of the large-eddy simulation benchmark. However, notable differences emerged in wake structure and unsteady load predictions, with discrepancies increasing for turbines further downstream. These differences highlight the importance of modelling choices such as wake summation and turbulence treatment, which strongly influence power-deficit and fatigue-load predictions. Comparison with large-eddy simulations reveals each approach's strengths and weaknesses, indicating where improvements are needed. Overall, the findings point to specific refinements for DWM models to improve their fidelity, ultimately enabling more robust wake predictions for wind farm design and operation.

17 WIND ENERGY↗
DOE OSTI2023

Representation of atmosphere-induced heterogeneity in land–atmosphere interactions in E3SM–MMFv2

In the Energy Exascale Earth System Model (E3SM) Multi-scale Modeling Framework (MMF), where parameterizations of convection and turbulence are replaced by a 2-D cloud-resolving model (CRM), there are multiple options to represent land–atmosphere interactions. Here, we propose three different coupling strategies, namely the (1) coupling of a single land surface model to the global grid (MMF), (2) coupling a single land copy directly to the embedded CRM (SFLX2CRM), and (3) coupling a single copy of land model to each column of the CRM grid (MAML). In the MAML (Multi-Atmosphere Multi-Land) framework, a land model is coupled to CRM at the CRM-grid scale by coupling an individual copy of a land model to each CRM grid. Therefore, we can represent intra-CRM heterogeneity in the land–atmosphere interaction processes. There are 5-year global simulations run using these three coupling strategies, and we find some regional differences but overall small changes with respect to whether a land model is coupled to CRM or a global atmosphere. In MAML, the spatial heterogeneity within CRM induces stronger turbulence, which leads to the changes in soil moisture, surface heat fluxes, and precipitation. However, the differences in the MAML from the other two cases are rather weak, suggesting that the impact of using MAML does not justify the increase in cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗