Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate boundary”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗
DOE OSTI2025

Turbulent Processes that Influence Boundary-Layer Cloud Structure (Final Technical Report)

Many models of the atmosphere simulate vertical profiles of observable fields, such as temperature, moisture, and vertical velocity. To assess the realism of such model output, it is important to evaluate the modeled profiles against observations. However, many extant evaluations consist of comparing time averaged profiles from the model and the observations. Such time averages can smooth out unrealistic features that appear in instantaneous profiles. In this way, the time averages may obscure certain model errors from detection. For instance, many models misrepresent turbulent processes, leading to profiles that are too smoothed out or too noisy. E.g., if a model has coarse vertical resolution and overly strong diffusion, then the model may be unable to produce observed sharp peaks. On the other hand, a model may have artificial approximations that produce spurious jumps in the vertical. Much can be learned by browsing through instantaneous snapshots of model output. But the volume of output produced by a global model is far too great to rely purely on visual inspection. In addition, we need a way to statistically characterize the profiles, so that anomalies can be detected more conveniently.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Evaluation of thermal energy storage for central receiver concentrating solar power plants under charging cycles through structural assessment

Electricity production by concentrated solar power (CSP) systems has stood out among energy transition alternatives due to their large generation capacity (between 30 MW and 400 MW) and high capacity factor (80%) when incorporating thermal storage systems (TES). Although TES tanks in operation worldwide are made of austenitic steels with remarkable mechanical and anticorrosive properties, operational failures have been reported due to low-cycle fatigue, creep, and the abrupt release of residual stresses generated during tank manufacturing. A critical operation for salt tanks is the initial daily charging operation, given the low fluid level and the thermal gradients between the inventory and the salt entering the tank. This study presents a structural simulation of a 39.6 m-diameter molten salt tank for central receiver CSP plants, accounting for temperature and pressure variations during charging and different sparger ring inlet configurations. The initial pre-stressed condition of the floor, resulting from the manufacture of welded plates, is used as a boundary condition. The fatigue and creep life analyses were performed in accordance with the ASME BPVC and API standards. It was found that the tank's floor is in a less favourable mechanical condition than the wall, due to its initial manufacturing condition and the maximum temperature differences during charging. A critical zone was identified in an area affected by residual stress on the floor and under the sparger ring. The baseline sparger ring configuration, with 52 2-inch-diameter orifices and a flow direction of 90 degrees , results in an average floor stress of approximately 25 MPa and creep damage after 1.2 years of tank operation.

14 SOLAR ENERGY↗
DOE OSTI2024

Three-Dimensional Aerodynamics and Vortex-Shedding Characteristics of Wind Turbine Airfoils over 360-Degree Angles of Attack

In this work, we present the first three-dimensional (3D) computational investigation of wind turbine airfoils over 360° angles of attack to predict unsteady aerodynamic loads and vortex-shedding characteristics. To this end, static–airfoil simulations are performed for the FFA-W3 airfoil family at a Reynolds number of 107 with the Improved Delayed Detached Eddy Simulation turbulence model. Aerodynamic forces reveal that the onset of boundary-layer instabilities and flow separation does not necessarily coincide with the onset of stall. In addition, a comparison with two-dimensional simulation data and flat plate theory extension of airfoil polars, suggest that, in the deep stall regime, 3D effects remain critical for predicting both the unsteady loads and the vortex-shedding dynamics. For all airfoils, the vortex-shedding frequencies are found to be inversely proportional to the wake width. In the case of slender airfoils, the frequencies are nearly independent of the airfoil thickness, and their corresponding Strouhal number St is approximately 0.15. Based on the calculated St, the potential for shedding frequencies to coincide with the natural frequencies of the International Energy Agency 15 MW reference wind turbine blades is investigated. The analysis shows that vortex-induced vibrations occur primarily at angles of attack of around ±90° for all airfoils.

17 WIND ENERGY↗
DOE OSTI2021

Modeling Coupled Chemo-Mechanical Behavior of Randomly Oriented NMC811 Polycrystalline Li-Ion Battery Cathodes

This paper develops a three-dimensional, transient, chemo-mechanical model that predicts the performance of single secondary particle Li-ion battery cathodes. The secondary particles are composed of numerous (approximately 60) randomly oriented single-crystal primary particles. The model incorporates concentration-dependent and anisotropic material properties. As much as possible, electrochemical, transport, and structural properties for crystalline NMC811 (Li x Ni 0.8 Mn 0.1 Co 0.1 O 2 ) are taken from prior publications. Weak Van der Waals bonding between primary particles is modeled empirically using a spring analogy, which enables local primary-particle separations (disintegration) and subsequent reattachments. The model fully couples Li diffusion and the mechanical response. Results include predictions of local Li-concentrations and stresses. High stresses are found near grain boundaries, especially when the lattice orientations are greatly misaligned. Particle separations are characterized in terms of a damage parameter. The model is used to predict the effects of design and operating conditions, including charge/discharge rates, cycling scenarios, and particle sizes.

25 ENERGY STORAGE↗
DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

CROCUS Sodar Measurements of Lower Atmospheric Wind Profiles at Argonne Testbed for Multiscale Observational Science (ATMOS) Site

The Scintec MFAS Sodar (Multiple-Frequency Acoustic Sounder) is an autonomous, ground-based acoustic remote sensing system designed to measure vertical profiles of horizontal wind speed, wind direction, and vertical velocity in the lower atmosphere. The instrument transmits sequences of acoustic pulses and detects the Doppler-shifted sound waves backscattered by small-scale temperature and velocity fluctuations caused by atmospheric turbulence. From these Doppler shifts, the system derives three-dimensional wind vectors by combining radial velocities from multiple beam orientations.The MFAS Sodar operates with a first usable range gate beginning at approximately 30 m above ground level and a configurable vertical resolution of 10 m. Under favorable acoustic conditions, the system provides wind profiles extending up to 600 m above ground level. Measurements are processed into 15-minute averaged profiles containing wind speed, direction, vertical velocity, and diagnostic quantities such as signal-to-noise ratio and echo strength.This dataset was collected at the Argonne Testbed for Multiscale Observational Science (ATMOS) facility in Lemont, Illinois, as part of DOE's CROCUS Urban Integrated Field Laboratory (UIFL) initiative. The purpose of these observations is to characterize the vertical wind structure and boundary-layer evolution across the urban–suburban gradient of the greater Chicago region. In particular, these data are intended to improve understanding of how local meteorology, such as lake-breeze penetration, nocturnal jets, and daytime mixing, varies between the densely built urban core and the suburban periphery. The MFAS observations provide critical context for evaluating high-resolution model simulations and for integrating with complementary lidar, radar, and in-situ meteorological measurements within the broader CROCUS UIFL network.All data are archived in NetCDF (Network Common Data Form) format and include wind and diagnostic parameters. The files can be accessed and analyzed using standard software that supports NetCDF, such as Python (e.g., xarray, netCDF4), MATLAB, R (e.g., ncdf4, raster), or Panoply (NASA’s NetCDF visualization application).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A reduced model for edge localized mode control by supersonic molecular beam injection and pellet injection

We develop a diffusive, bistable, tri-unstable cellular automata (CA) model to study the dynamics of H-mode pedestal with edge localized modes (ELMs) and their control by supersonic molecular beam injection (SMBI) and pellet injection (PI). It is shown that the new CA model can reproduce the key features of H-mode pedestals with various types of ELM, including Type-I ELM. SMBI and PI are modeled as additional grain injections into pedestal with varying degrees of injected materials and profiles. It is found that H-mode pedestal responds to SMBI differently depending on the baseline fueling. If the baseline fueling is large enough to allow Type-I ELM, SMBI enhances large transport avalanches caused by ballooning instabilities. These avalanches prevent the total pedestal current from reaching the boundary for peeling instability. On the other hand, if the baseline fueling is low to avoid Type-I ELM, SMBI enhances small scale avalanches, which prevent the pedestal from growing to profiles globally vulnerable to ballooning instabilities. These imply that SMBI can mitigate different types of ELM by converting them to more benign types. From CA modeling of pellet injection, it is shown that Type-I ELM can be triggered by pellet injection with sufficient strength and depth. Scanning the frequency of pellet injection, it is found that a maximum efficiency of pellet pacing is achieved when the injection frequency is approximately ten times the natural frequency of Type-I ELM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An Infinite Domain 3D Poisson Solver Based on the Barnes-Hut Algorithm

We present a domain decomposition method for the solution of the 3D Poisson equation with infinite domain boundary conditions. The method is based on an application of the Barnes-Hut tree particle scheme adapted to gridded data. Long range interactions are computed using the first two terms in the Cartesian multipole expansion of Green’s function convoluted with the charge while short range computations are performed using Hockney’s domain doubling algorithm. A standard domain decomposition strategy requires O(N 2 ) applications of Hockney’s algorithm, where N is the number of subdomains that intersect that charge support, while in the present approach only O(Nlog 2 N) such computations suffice. The discretization scheme employed is a sixth order Mehrstellen approximation of the 3D Laplace opera tor. The method exhibits satisfactory accuracy at a substantially reduced computational cost compared to the full domain decomposition Hockney’s algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Exploration of Sn 70 Li 30 alloy as possible material for flowing liquid metal plasma facing components

As an advanced alternative to solid materials, Liquid Metals (LM) may offer more resilient and feasible Plasma Facing Components (PFCs). Particularly, regarding the unavoidable material erosion/degradation produced by particle/heat fluxes in future fusion devices where much longer duty cycles are expected. Furthermore, configurations that propose a flowing LM surface can add the advantage of a continuously fresh and clean layer facing the plasma. Although lithium is the most widely tested option, tin-lithium (SnLi) alloys have been proposed to attempt to combine the positive characteristics of both pure elements and ameliorate the specific issues of lithium. In this work, the potential use of Sn70Li30 alloy in such flowing concepts has been explored by addressing several preliminary and mandatory aspects for its utilization. Key issues such as wettability and compatibility of the alloy with relevant substrates have been studied in a multidisciplinary approach. The data obtained from deposited liquid tin-lithium droplets indicates approximate wetting temperatures of 360 °C, 390 °C and 405 °C for the fresh alloy on 316 stainless steel, molybdenum, and tungsten, respectively. However, the alloy contamination appeared to strongly affect the wetting characteristics of materials, increasing their wetting temperature by ~130 °C in the worst observed cases. Interestingly, in some instances, the instability of the liquid alloy surface was observed in the form of sudden gaseous ejection. The deposited droplets were posteriorly characterized in terms of absolute composition and depth profile by Inductively Coupled PlasmaOptical Emission Spectroscopy (ICP-OES) and Secondary Ion Mass Spectrometry (SIMS-ToF). Additionally, the nature and composition of the boundaries between the substrates and alloy microparticles was investigated by Scanning Electron Microscopy (SEM), Energy Dispersive X-Ray Spectroscopy (EDS), and 3D Laser microscopy. The overall results of this post-mortem characterization revealed that first signs of corrosion induced by both alloy elements (lithium-chromium association and iron-tin intermetallic mixing) were present on 316 stainless steel after short exposures (≤3 h) at temperatures lower than 550 °C. Conversely, molybdenum and tungsten showed good compatibility with the alloy in equivalent conditions. The global implications of these results are finally addressed, focusing on the future perspectives and the more viable scenarios for the eventual utilization of these alloys in flowing liquid metal configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The effect of microstructure on the dynamic shock response of 1045 steel

Herein, the effect of microstructure on the shock response of 1045 steel is investigated via plate impact experiments and postmortem characterization. Three unique microstructures are explored: ferrite-pearlite, martensite, and ferrite with spheroidal cementite (i.e. spheroidized). Two spall recovery experiments, at approximate peak pressures of 3.2 and 3.5 GPa, are conducted to assess the Hugoniot elastic limit (HEL), spall strength, and damage morphology of the various microstructures. The ferrite-pearlite and martensite microstructures exhibit continuous yielding at both quasi-static and dynamic rates, while the spheroidized condition displays discontinuous yielding. Discontinuous yielding of the spheroidized microstructure is attributed to a combined low initial dislocation density coupled with a low dislocation nucleation rate. The spall strength of ferrite-pearlite is consistently lower than the spheroidized microstructure, attributed to elongated cementite that is more susceptible to cracking than more spherical cementite precipitates. Despite a high density of boundaries, martensite exhibits the highest spall strength. A large percentage of the boundaries within the martensite microstructure are found to be low energy (i.e. Σ3 or low angle), and are thus less susceptible to spall damage. Overall, the high spall strength of martensite is likely linked to traditional strengthening mechanisms that limit dislocation motion.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Acoustic Codes in 2D Spherical Coordinate

Finite-difference methods are widely used to simulate infrasound propagation in the atmosphere. Flexibility of finite-difference scheme allows implementation of highly heterogeneous media for sound propagation as well as complex source models for sound generation. While full 3-D finite-difference methods have been utilized for local infrasound propagation with pronounced topography, 2-D modeling approach has been preferred for regional and global propagation as full 3-D methods generally require enormous computational resources. Infrasound propagation is often simulated with a second-order finite difference scheme. This lowest-order finite-difference scheme is robust and straightforward to implement complex boundary conditions, but the solution includes large error with numerical dispersion and dissipation. This large numerical error may make the second-order finite-difference not optimal for long range infrasound propagation modeling as the numerical dispersion degrades the accuracy of the solution unacceptably. Here, we developed a high-order finite-difference solver for long-range infrasound simulation. The high-order scheme is particularly popular for linear wave modeling in aeroacoustics owing to its low-dispersive and low-dissipative behavior. We develop and evaluate a high-order finite difference scheme in 2-D axisymmetric coordinates. The axisymmetry allows to approximate 3-D spherical sound propagation and amplitude attenuation by a 2-D method. AC2Dr is developed to simulate infrasound propagation in realistic atmosphere, but can be used for linear acoustic waves in general materials with background flow. AC2Dr in an axisymmetric coordinates allows for spherical radiation of acoustic waves from compact sources.

Sjogreen, Bjorn↗
DOE OSTI2025

An apologia for islands

Entanglement islands have played a key role in the recent derivation of the Page curve and other progress on the black hole information problem. Arising from the inclusion of connected wormhole saddles in a gravitational replica trick, islands signal that degrees of freedom in the black hole interior are not microscopically independent of the exterior Hawking radiation. Islands were originally discovered in the context of AdS/CFT coupled to an external, nongravitating reservoir, where the coupling gives graviton excitations an anomalous boundary scaling dimension (or “mass”). It has been claimed in the literature that this mass is crucial for the existence of islands and even the Page curve itself. In this paper, however, we explain how entanglement islands can also appear in setups with massless gravitons and no external reservoir, giving a number of examples including the entanglement wedges of boundary CFT regions, of radiation at null infinity in asymptotically flat spacetimes, and of radiation inside a semiclassical but gravitating spacetime. In each case, the Page curve is physically observable and can be determined with sufficiently careful experiments on many copies of the black hole. We give general arguments for the existence of gauge-invariant operators in gravity which are compactly supported to all orders in perturbation theory (whenever no isometries of the background spacetime exist) and refine a recently-proposed explicit construction of such operators. When applied to islands, these results — together with entanglement wedge reconstruction — guarantee that semiclassical operators in the island can be approximated by nonperturbative operators on the Hawking radiation.

AdS-CFT Correspondence↗
DOE OSTI2020

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Acoustic Tomography of the Atmosphere: A Large-Eddy Simulation Sensitivity Study

Accurate measurement of atmospheric turbulent fluctuations is critical for understanding environmental dynamics and improving models in applications such as wind energy. Advanced remote sensing technologies are essential for capturing instantaneous velocity and temperature fluctuations. Acoustic tomography (AT) offers a promising approach that utilizes sound travel times between an array of transducers to reconstruct turbulence fields. This study presents a systematic evaluation of the time-dependent stochastic inversion (TDSI) algorithm for AT using synthetic travel-time measurements derived from large-eddy simulation (LES) fields under both neutral and convective atmospheric boundary-layer conditions. Unlike prior work that relied on field observations or idealized fields, the LES framework provides a ground-truth atmospheric state, enabling quantitative assessment of TDSI retrieval reliability, sensitivity to travel-time measurement noise, and dependence on covariance model parameters and temporal data integration. A detailed sensitivity analysis was conducted to determine the best-fit model parameters, identify the tolerance thresholds for parameter mismatch, and establish a maximum spatial resolution. The TDSI algorithm successfully reconstructed large-scale velocity and temperature fluctuations with root mean square errors ( RMSE s) below 0.35 m/s and 0.12 K, respectively. Spectral analysis established a maximum spatial resolution of approximately 1.4 m, and reconstructions remained robust for travel-time measurement uncertainties up to 0.002 s. These findings provide critical insights into the operational limits of TDSI and inform future applications of AT for atmospheric turbulence characterization and system design.

17 WIND ENERGY↗
DOE OSTI2025

Comparative Evaluation of the Ability of the MYNN‐EDMF PBL Scheme in WRF Model to Reproduce Near Surface Wind Speed Over Different Topographical Types

Abstract This study systematically evaluates the performance of the Mellor‐Yamada‐Nakanishi‐Niino‐Eddy‐Diffusion‐Mass‐Flux planetary boundary layer (PBL) scheme within the Weather Research and Forecasting (WRF) model in simulating near‐surface wind speeds across various topographies in New York State (NYS). Simulated wind speeds are compared with in‐situ measurements from 22 surface sites, grouped into six topographic categories: continental plain (CT), lakeside (LS), river valley (RV), Long Island (LI), Block Island (BI), and offshore ocean (OO). A quantitative evaluation based on Relative Euclidean Distance shows that wind speeds at the OO site are the most accurately reproduced, followed by those at LI sites, while the model performs less accurately for the remaining topographic groups. Wind speeds over CT sites tend to be overestimated by approximately 1 m/s, although their diurnal variability (DV) is well captured. In contrast, the model underestimates wind DV at LS, RV, LI, and BI sites, with the largest biases occurring at LI and BI, resulting in underestimated daytime wind speed and/or overestimated nighttime wind speed. The OO winds exhibit minimal diurnal variation, accurately captured by our WRF model. The surface wind diurnal variation is closely linked to PBL development. Among the indicators of PBL development, surface potential temperature biases most strongly correlate with wind speed biases. Our WRF model faces challenges in capturing the distinctions between winds influenced by local circulations and those over continental plains, and the significantly stronger winds at OO compared to BI. Potential causes for these biases are discussed, offering pathways for improving surface wind simulations in future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗