SEARCH · Search NASA

Results for “Boundary layer structure”

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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Nanometer-Scale Imaging of Inhomogeneous Active Charge Carriers in Arsenic-Doped CdTe Thin Films

We report nanometer-scale imaging of active carrier distribution of As-doped CdTe films by scanning capacitance microscopy (SCM). We developed SCM sample preparation for CdTe by ion-milling followed by thermal processing. The nanometer-resolution carrier delineation for CdTe was validated by imaging on a CdTe cross-section sample made by a molecular beam epitaxy layer stack with As-doping concentrations of 10^15 ~10^18 /cm 3 . We found that the carrier distribution in As-doped films was significantly nonuniform, with inhomogeneity ranging from sub-um to a few um and concentration variation of one order of magnitude (low 10^16 to low 10^17 /cm 3 ). This nonuniformity is distributed randomly, independent of grain structure and grain boundary (GB). We used Kelvin probe force microscopy (KPFM) and cathodoluminescence (CL) to further map the surface potential and radiative illumination on the same area as the SCM image. Higher potential and lower CL intensity were found on GBs but not on SCM contrast, illustrating positive GB charging and GB recombination but not GB-distinguished doping. The overall KPFM potential image is in rough agreement with the SCM carrier distribution, in terms of Fermi-level position relative to the bandgap edge - thus resulting in the band-edge potential fluctuation. Nonuniform carrier concentration, potential fluctuation, and defect recombination can all together cause the Voc deficit of the As-doped CdTe device.

14 SOLAR ENERGY↗

DOE OSTI2024

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗

DOE OSTI2021

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis

In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.

58 GEOSCIENCES↗

DOE OSTI2022

The complex structural and chemical nature of monolithic U-10Mo fuel and Zr barrier layer

Nanoscale microstructural characterization by advanced transmission electron microscopy techniques on a U-10Mo/Zr barrier layer monolithic fuel plate was performed to evaluate the microstructural evaluation after high burn-up. Gas bubble superlattice evolution, grain restructuring, and evolution of the Zr interaction layer is investigated through detailed electron microscopy characterization. The use of automated crystallographic orientation mapping to irradiated U-10Mo fuel highlights that the restructured ultra-fine grains are separated by high angle grain boundaries at a burn up of 4.42 × 10 21 fissions/cm 3 . Additionally, advanced chemical analysis and multi-variable statistical analysis shows spatial clustering of solid fission product precipitates. Finally, characterization of a newly observed porous nanocrystalline Zr region in the barrier layer is studied. Finally, this work provides insights into the grain subdivision and restructuring process while using advanced microscopy techniques to analyze fission products in neutron-irradiated U-10Mo fuel.

36 MATERIALS SCIENCE↗

DOE OSTI2020

He ion irradiation response of a gradient T91 steel

Metallic materials with a gradient microstructure usually exhibit excellent mechanical properties. How- ever, the radiation response of gradient structural materials is less well understood. In this work, room- temperature He ion irradiation up to ~4.5 dpa with ~10 at% He injection was performed on a gradient T91 steel processed by surface severe plastic deformation. In comparison to the coarse-grained ferritic T91, the gradient T91 with the nanocrystalline layers shows improved radiation tolerance in terms of less bubble swelling and radiation hardening. Additionally, bubble distribution along grain boundaries depends on misorientation angle suggesting a strong capacity of non-equilibrium grain boundaries in storing He atoms. The present study provides insight into the design of radiation tolerant gradient steels for nuclear industry applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effects of Extended Aqueous Processing on Structure, Chemistry, and Performance of Polycrystalline LiNi x Mn y Co z O 2 Cathode Powders

The prospect of aqueous processing of LiNi x Mn y Co z O 2 (NMC) cathodes has significant appeal to battery manufacturers for the reduction in materials cost, toxicological risk, and environmental impact compared to conventional N -methyl-2-pyrrolidone (NMP)-based processing. However, the effects of aqueous processing of NMC powders at industrial timescales are not well studied, with prior studies mostly focusing on relatively brief water washing processes. Here, we investigate the bulk and surface impacts of extended aqueous processing of polycrystalline NMC powders with different compositions. We demonstrate that at timescales of several hours, polycrystalline NMC is susceptible to intergranular fracture, with the severity of fracture scaling with the NMC nickel content. While bulk crystallinity and composition are unchanged, surface sensitive techniques such as X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) indicate that the exposure of water leads to a level of delithiation, nickel reduction, and reconstruction from the layered to rock-salt structure at the surface of individual grains. Dynamic single NMC microparticle compression testing suggests that the resulting mechanical stresses weaken the integrity of the polycrystalline particle and increases susceptibility of intergranular fracture. The initially degraded surfaces along with the increased surface area lead to faster capacity fade and impedance growth during electrochemical cycling. From this work, it is demonstrated that NMC powders require surface or grain boundary modifications to make industrial-scale aqueous cathode processing viable, especially for next-generation nickel-rich NMC chemistries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Variations in GARS powder microstructure as a function of powder chemistry and particle size

The properties of metal produced through powder metallurgy depends on the feedstock used. Powders produced via gas atomization reaction synthesis (GARS) are used to produce oxide dispersion strengthened alloys. The desired powder size range can vary for each consolidation technique. However, powder microstructure also can vary with powder particle size, which in turn can impact the microstructure and properties of the consolidated parts. In this study, GARS powders are characterized via inductively coupled plasma mass spectroscopy, inert gas fusion, and high-resolution x-ray diffraction to determine variations in elemental and phase compositions. Transmission electron microscopy was used to understand microstructure variations as a function of chemistry and size. Across the three batches tested intermetallic content was 0.73–1.35 wt% in the 0-20 μm powder batch and increased to 2.46–3.80 wt% in the coarse 45-106 μm batch. Across all batches, volume percent of surface oxidation decreased with powder diameter, with volume percents within the range of 0.75–1.2 % across 10 μm powder particles, and below 0.4 % across coarse powder particles approximately 100 μm in diameter. These observations were supported by inert gas fusion measurements. However, the oxide layer was thicker in coarse powder particles due to a slower cooling rate. Increasing oxygen content in atomization gas to 2000 ppm and adding yttrium increased both the surface oxidation content and yttrium intermetallic content. Lastly, intermetallic phases within the powder coarsened with powder size. Intermetallic morphology changed from fine spherical intermetallic and columnar dendritic growth to a cellular structure with finer spherical intermetallic, to coarse irregular intermetallic and intermetallic along grain boundaries as a result of slower cooling rate and solidification rate in coarse powder particles. Furthermore, the addition of zirconium does not appear to significantly change intermetallic morphology, but the composition changed from a Y-Fe rich intermetallic to a Y-Zr-Fe intermetallic.

42 ENGINEERING↗

DOE OSTI2024

Enhancing the Reaction Kinetics and Stability of Co‐Free Li‐Rich Cathode Materials via a Multifunctional Strategy

Co-free Li-rich layered oxides (CFLLOs) with anionic redox activity are among the most promising cathode materials for high-energy-density and low-cost lithium-ion batteries (LIBs). However, irreversible oxygen release often causes severe structural deterioration, electrolyte decomposition, and the formation of unstable cathode-electrolyte interface (CEI) film with high impedance. Additionally, the elimination of cobalt elements further deteriorates the reaction kinetics, leading to reduced capacity and poor rate performance. Here, in this study, a multifunctional strategy is proposed, incorporating Li 2 MnO 3 phase content regulation, micro-nano structure design, and heteroatom substitution. The increased content of Li 2 MnO 3 phase enhances the capacity through oxygen redox. The smaller nanoscale primary particles induce greater tensile strain and introduce more grain boundaries, thereby improving the reaction kinetics and reactivity, while the larger micron-sized secondary particles help to reduce interfacial side reactions. Furthermore, Na⁺ doping modulates the local coordination environment of oxygen, stabilizing both the anion framework and the crystal structure. As a result, the designed cathode exhibits enhanced rate performance, delivering a capacity of 158 mAh g⁻¹ at 5.0 C and improved cyclic stability, with a high capacity retention of 99% after 400 cycles at 1.0 C. This multifunctional strategy holds great promise for advancing the practical application of CFLLOs in next-generation LIBs.

Co-free Li-rich layered oxide↗

DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗

DOE OSTI2022

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗

DOE OSTI2023

Distribution of Copper States, Phases, and Defects across the Depth of a Cu-Doped CdTe Solar Cell

Copper has been used as a p-type dopant in cadmium telluride (CdTe) for decades. However, the density of Cu atoms in the finished device is much higher than that of holes, which means that most Cu atoms are not activated as acceptors during incorporation. Furthermore, studies have demonstrated that the distribution of copper (Cu) atoms across the device is highly inhomogeneous, with reports citing Cu substitution on Cd sites and segregation to grain boundaries. Fast diffusion along these boundaries and Cu accumulation at the CdTe/CdS interface have also been observed and validated computationally. These levels of inhomogeneity make it difficult to accurately characterize and correlate the performance with the nature of the Cu atomic species present. To address this challenge, we utilize X-ray microscopy and, specifically, nanoscale fluorescence-mode X-ray absorption near-edge structure to resolve the atomic Cu environment throughout the depth of the CdTe layer. Our results suggest that the majority of Cu atoms are in the form of Cu x Te phases (or similar local environments) near the ZnTe|CdTe interface, Cu x O phases in the CdTe absorber, and present in various oxidation states, including Cu 1+ and Cu 2+ , near the CdS/CdTe junction. Here this work also provides experimental evidence for the first time of the presence of CuS around the ZnTe|CdTe interface and the hypothesized Cu Cd -Cl i complex in the CdTe absorber.

14 SOLAR ENERGY↗

DOE OSTI2026

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗

DOE OSTI2020

Evidence of Inverse Hall-Petch Behavior and Low Friction and Wear in High Entropy Alloys

We present evidence of inverse Hall-Petch behavior for a single-phase high entropy alloy (CoCrFeMnNi) in ultra-high vacuum and show that it is associated with low friction coefficients (~0.3). Grain size measurements by STEM validate a recently proposed dynamic amorphization model that accurately predicts grain size-dependent shear strength in the inverse Hall-Petch regime. Wear rates in the initially soft (coarse grained) material were shown to be remarkably low (~10 –6 mm 3 /N-m), the lowest for any HEA tested in an inert environment where oxidation and the formation of mixed metal-oxide films is mitigated. The combined high wear resistance and low friction are linked to the formation of an ultra-nanocrystalline near-surface layer. The dynamic amorphization model was also used to predict an average high angle grain boundary energy (0.87 J/m 2 ). This value was used to explain cavitation-induced nanoporosity found in the highly deformed surface layer, a phenomenon that has been linked to superplasticity.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

DOE OSTI2025

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

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

17 WIND ENERGY↗