Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY LAYER STABILITY”

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
DOE OSTI2023

Interplay of molten salt and dopants in tuning the performance of single-crystalline LiNiO 2

Single-crystalline layered oxide cathodes have seen great research interest in recent years. Adoption of single-crystal morphology improves mechanical robustness and eliminates intergranular boundaries, improving the chemical stability of the cathode. Furthermore, single-crystal cathodes also suffer from worsened lithium-diffusion kinetics stemming from the significantly increased diffusion path lengths. Sluggish lithiation kinetics induce concentration gradients within the cathode particle during cell operation, resulting in the accumulation of lattice strain, which deteriorates cell performance. Foreign dopants have been shown to greatly improve ion diffusion in many high-Ni layered oxide cathode systems. In particular, we identify Ti and Te as promising dopants to ameliorate the slow diffusion kinetics. We find that doping 1 mol.% Te into single crystal LiNiO 2 greatly improves Li + diffusivity, as determined through galvanostatic intermittent titration technique. Further incorporating 2 mol.% Al and pairing with lithium-metal anode in a localized high concentration electrolyte (LHCE) delivers 206 mA h g -1 at C/3 rate and retains 94% capacity after 100 cycles in Li metal half cells. By addressing kinetically mediated degradation with Te doping, and surface-based degradation with Al and LHCE, the cathode experiences very little polarization and overpotential growth, enabling stable cycling with high accessible capacity.

25 ENERGY STORAGE↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
DOE OSTI2024

The Role of Non-Local Effects on Surface Sensible Heat Flux Under Different Types of Thermal Structures Over the Arctic Sea-Ice Surface

The effects of atmospheric thermal structure on the surface energy flux are poorly understood over the Arctic sea-ice surface. Here, we explore the mechanism of sensible heat exchange under different types of thermal structures over the Arctic sea-ice surface by using data collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. The quadrant analysis indicates that strong surface temperature inversions below 100 m enhance non-local effects on the positive (upward) sensible heat flux ($\overline{w'θ'}$) through entrainment of large eddies from the convective boundary layer aloft. However, strong surface inversions restrict the contributions of large eddies to the negative (downward) $\overline{w'θ'}$ due to intensified surface stability. By inspecting the existing parameterization schemes, we found that the European Center for Medium-Range Weather Forecasts Integrated Forecasting System scheme fails to predict the impacts of non-local processes on the positive $\overline{w'θ'}$, and an adjustment term to correct the bias of parameterized $\overline{w'θ'}$ is proposed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Editorial: Resolving atmospheric flow in complex environments: recent experiments in terrain and forest canopies

The characterization of atmospheric flows in complex environments, which may include steep terrain slopes and heterogeneous vegetation and/or forest cover, is a long-standing challenge in boundary-layer meteorology. Atmospheric observations are complicated by the presence of transient, terrain-induced flow features, forest-canopy-atmosphere interactions, and atmospheric stability effects, not to mention the logistical hurdles involved with instrument deployment, data analysis, and quality control. Furthermore, challenges in atmospheric modeling arise due to numerical errors associated with complex terrain flows, as well as reliance on simplified parameterizations for unresolved processes such as turbulent mixing and land-surface or forest-canopy-atmosphere interactions. These modeling challenges are exacerbated in the so-called “gray zone,” wherein features of interest have length scales that are similar to the model grid spacing, or when the principal flow layer is smaller than the grid spacing (e.g., slope flows).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterizing Seasonal Variation of the Atmospheric Mixing Layer Height Using Machine Learning Approaches

As machine learning becomes more integrated into atmospheric science, XGBoost has gained popularity for its ability to assess the relative contributions of influencing factors in the atmospheric boundary layer height. To examine how these factors vary across seasons, a seasonal analysis is necessary. However, dividing data by season reduces the sample size, which can affect result reliability and complicate factor comparisons. To address these challenges, this study replaces default parameters with grid search optimization and incorporates cross-validation to mitigate dataset limitations. Using XGBoost with four years of data from the atmospheric radiation measurement (ARM) (Southern Great Plains (SGP) C1 site, cross-validation stabilizes correlation coefficient fluctuations from 0.3 to within 0.1. With optimized parameters, the R value can reach 0.81. Analysis of the C1 site reveals that the relative importance of different factors changes across seasons. Lower tropospheric stability (LTS, ~0.53) is the dominant factor at C1 throughout the year. However, during DJF, latent heat flux (LHF, 0.44) surpasses LTS (0.22). In SON, LTS (0.58) becomes more influential than LHF (0.18). Further comparisons among the four long-term SGP sites (C1, E32, E37, and E39) show seasonal variations in relative importance. Notably, during JJA, the differences in the relative importance of the three factors across all sites are lower than in other seasons. This suggests that boundary layer development in the summer is not dominated by a single factor, reflecting a more intricate process likely influenced by seasonal conditions such as enhanced convective activity, higher temperatures, and humidity, which collectively contribute to a balanced distribution of parameter impacts. Furthermore, the relative importance of LTS gradually increases from morning to noon, indicating that LTS becomes more significant as the boundary layer approaches its maximum height. Consequently, the LTS in the early morning in autumn exhibits greater relative importance compared to other seasons. This reflects a faster development of the mixing layer height (MLH) in autumn, suggesting that it is easier to retrieve the MLH from the previous day during this period. The findings enhance understanding of boundary layer evolution and contribute to improved boundary layer parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Experimental and Theoretical Studies of the Surface Oxidation Process of Rare–Earth Tritellurides

Recent studies have established Van der Waals (vdW) layered and 2D rare-earth tritellurides (RTe 3 ) as superconductors and near room-temperature charge density wave (CDW) materials. Their environmental stability raises natural concern owing to aging/stability effects observed in other tellurium-based layered crystals. Here, the results establish the stability and environmental aging characteristics of these RTe 3 systems involving a variety of metals such as La, Nd, Sm, Gd, Dy, and Ho. The atomic force microscopy (AFM) and scanning electron microscopy (SEM) results show that all the RTe 3 sheets oxidize to form thin TeO x layers that are primarily confined to the surface, edges, and grain boundaries. Time-resolved in situ Raman spectroscopy measurements are used to understand the kinetics of the oxidization process for different lanthanide metal cations and establish their relative stability/resilience to oxidization. Overall results indicate that the vdW layers show higher air stability as the 4f electron number decreases going from Ho to La, resulting in the most stable LaTe 3 compared to the least stable HoTe 3 . Comprehensive quantum mechanical simulations reveal that environmental degradation originates from a strong oxidizing reaction with O 2 molecules, while humidity (H 2 O) plays a negligible role unless Te vacancies are present. Moreover, the simulations explain the effects of 4f electrons on the work function and Te vacancies formation, which directly impact the aging characteristics of RTe 3 layers. Interestingly, optical and electrical measurements show that the CDW response is still observed in aged RTe 3 layers owing to the presence of underlying pristine/nonoxidized RTe 3 layers, except CDW transition temperatures increase due to the thickness effect. Overall results offer the first in-depth environmental aging studies on these materials, which can be applied to engineer and design their chemical stability, surface properties, and overall CDW characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhancing Crystallization in Hybrid Perovskite Solar Cells Using Thermally Conductive 2D Boron Nitride Nanosheet Additive

Controlling crystallization and grain growth is crucial for realizing highly efficient hybrid perovskite solar cells (PSCs). In this work, enhanced PSC photovoltaic performance and stability by accelerating perovskite crystallization and grain growth via 2D hexagonal boron nitride (hBN) nanosheet additives incorporated into the active perovskite layer are demonstrated. Specifically, in situ X-ray scattering and infrared thermal imaging during the perovskite annealing process revealed the highly thermally conductive hBN nanosheets promoted the phase conversion and grain growth in the perovskite layer by facilitating a more rapid and spatially uniform temperature rise within the perovskite film. Complementary structural, physicochemical, and electrical characterizations further showed that the hBN nanosheets formed a physical barrier at the perovskite grain boundaries and the interfaces with charge transport layers, passivating defects, and retarding ion migration. As a result, the power conversion efficiency of the PSC is improved from 17.4% to 19.8%, along with enhanced device stability, retaining ≈90% of the initial efficiency even after 500 h ambient air storage. Here, the results not only highlight 2D hBN as an effective additive for PSCs but also suggest enhanced thermal transport as one of the pathways for improved PSC performance by 2D material additives in general.

2D materials↗
DOE OSTI2021

Reactive boride infusion stabilizes Ni-rich cathodes for lithium-ion batteries

Engineered polycrystalline electrodes are critical to the cycling stability and safety of lithium-ion batteries, yet it is challenging to construct high-quality coatings at both the primary- and secondary-particle levels. Here, we present a room-temperature synthesis route to achieve full surface coverage of secondary particles and facile infusion into grain boundaries, thus offering a complete “coating-plus-infusion” strategy. Cobalt boride metallic glass is successfully applied to Ni-rich layered cathode LiNi 0.8 Co 0.1 Mn 0.1 O 2 . Here, it dramatically improves the rate capability and cycling stability, including under high-discharge-rate and elevated-temperature conditions and in pouch full cells. The superior performance originates from simultaneous suppression of microstructural degradation of intergranular cracking and side reactions with electrolyte. Atomistic simulations identified the critical role of strong selective interfacial bonding, which offers not only a large chemical driving force to ensure uniform reactive wetting and facile infusion but also lowered the surface/interface oxygen activity, contributing to the exceptional mechanical and electrochemical stabilities of the infused electrode.

25 ENERGY STORAGE↗
DOE OSTI2024

Engineering 3D Floquet Codes by Rewinding

Floquet codes are a novel class of quantum error-correcting codes with dynamically generated logical qubits arising from a periodic schedule of noncommuting measurements. We utilize the interpretation of measurements in terms of condensation of topological excitations and the rewinding of measurement sequences to engineer new examples of Floquet codes. In particular, rewinding is advantageous for obtaining a desired set of instantaneous stabilizer groups on both toric and planar layouts. Our first example is a Floquet code with instantaneous stabilizer codes that have the same topological order as the three-dimensional (3D) toric code(s). This Floquet code also exhibits a splitting of the topological order of the 3D toric code under the associated sequence of measurements, i.e., an instantaneous stabilizer group of a single copy of the 3D toric code in one round transforms into an instantaneous stabilizer group of two copies of the 3D toric code up to nonlocal stabilizers in the following round. We further construct boundaries for this 3D code and argue that stacking it with two copies of the 3D subsystem toric code allows for a transversal implementation of the logical non-Clifford controlled-controlled- Z gate. We also show that the coupled-layer construction of the X-cube Floquet code can be modified by a rewinding schedule such that each of the instantaneous stabilizer codes is finite depth equivalent to the X-cube model up to toric codes; the X-cube Floquet code exhibits a splitting of the X-cube model into a copy of the X-cube model and toric codes under the measurement sequence. Our final 3D example is a generalization of the 2D Floquet toric code on the honeycomb lattice to three dimensions, which has instantaneous stabilizer codes with the same topological order as the 3D fermionic toric code. Published by the American Physical Society2024

Physics↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Addressing Chemical and Thermomechanical Instabilities in Metal Halide Perovskites Across Length Scales

The chemical and thermomechanical instabilities of metal halide perovskites represent barriers to long operational lifetimes of perovskite-based devices. There is no single magic bullet to address these instabilities, and so multiple complementary strategies must be undertaken. I will discuss the work by our research group to address these instabilities at multiple length scales, including the effects of the molecular structure and functionality of additives on their interactions with perovskite surfaces, the enhancement of mechanical toughness of perovskite layers by compositing with polymeric nanofibers, and the use of confined volume printing to rapidly produce perovskite films with large grains and reduced surface area of grain boundaries detrimental to stability. Together, these approaches provide one possibly pathway toward significantly increasing the operational lifetimes of perovskite-based devices.

Printz, Adam↗
DOE OSTI2023

Enhancing the performance of the perovskite solar cells by modifying the SnO 2 electron transport layer

The tin oxide (SnO 2 ) electron transport layer (ETL) plays a vital role in the photo-conversion efficiency (PCE) and stability of organic-inorganic perovskite solar cells (PSCs). However, SnO 2 ETL-induced defects such as hydroxyl groups, oxygen vacancies, exposed Sn atoms, and dangling bonds hinder device performance. In this study, rubidium chloride (RbCl) has been used to modify the SnO 2 ETL. Perovskite film formed on the RbCl-modified SnO 2 ETL exhibits improved crystallinity with enlarged grain size and reduced grain boundaries and enhanced optical absorption. Further, the Hall-effect measurements indicate the improved carrier mobility, and the dark J-V curve shows the increment of electrical conductivity for the RbCl-modified SnO 2 ETL. X-ray photoelectron spectroscopy (XPS) results demonstrate the surface defects passivation of the perovskite layer by modifying the SnO 2 ETL. A champion PCE of 19.35% has been achieved for the RbCl-modified SnO 2 ETL-based devices with improved stability, while the control devices with unmodified SnO 2 ETL show a PCE of 17.18%.

14 SOLAR ENERGY↗
DOE OSTI2025

Nucleation and Antiphase Twin Control in Bi 2 Se 3 via Step‐Terminated Al 2 O 3 Substrates

The epitaxial synthesis of high-quality 2D layered materials is an essential driver of both fundamental physics studies and technological applications. Bi 2 Se 3 , a prototypical 2D layered topological insulator, is sensitive to defects imparted during the growth, either thermodynamically or due to the film-substrate interaction. Here, in this study, it is shown that step-terminated Al 2 O 3 substrates with a high miscut angle (3°) can effectively suppress a particular hard-to-mitigate defect, the antiphase twin. Systematic investigations across a range of growth temperatures and substrate miscut angles confirm that atomic step edges act as preferential nucleation sites, stabilizing a single twin domain. First-principles calculations suggest that there is a significant energy barrier for twin boundary formation at step edges, supporting the experimental observations. Detailed structural characterization indicates that this twin-selectivity is lost through the mechanism of the 2D layers overgrowing the step edges, leading to higher twin density as the thickness increases. These findings highlight the complex energy landscape unique to 2D materials that is driven by the interplay between substrate properties, nucleation dynamics, and defect formation, and overcoming and controlling these are critical to improve material quality for quantum and electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Characterization of Flashback and Flame-Holding in a Jet-In-Crossflow Mixing Configuration with Methane-Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence (CL) diagnostics techniques. The JICF burner has an optically accessible pre-mixing tube allowing the optical diagnostics. The testing was conducted for varied pre-mixer velocities (V) and equivalence ratio (ϕ) for 90%-100% (H2, by mole) H2/CH4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashbacks were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower flux ratio which bends the jet sharply, closer toward the injection plane. The flashback was found to occur along the boundary layer, assisted by the inner wall of the pre-mixer tube. The mean OH-PLIF images characterized the flame-holding behavior where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port.

flame holding↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗