Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuity Factors”

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.

72 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Unveiling Swift Heavy Ion Track Morphology in Sr-Based High-Entropy Perovskites

The incorporation of multiple cations on a single lattice site in the high-entropy oxides is considered the key driving factor for modifying the known atomic-level response to energetic ion irradiation due to the presence of structural disorder; however, these effects are not well-understood yet. In this work, we present atomic-level insight into irradiation-induced nanoscale phase transformations in a perovskite-structured high-entropy oxide, Sr(Zr 0.2 Sn 0.2 Ti 0.2 Hf 0.2 Nb 0.2 )O 3 (Sr(HE)O 3 ), subjected to 774 MeV swift Xe heavy ions, where damage is dominated by inelastic ion−lattice interactions. While these ions generally are known to create nanoscale disordered channels, “ion tracks”, along the penetration direction in the material, this study shows the formation of discontinuous and partially recrystallized ion tracks in Sr(HE)O 3 . Compared to SrTiO 3 irradiated under identical energy loss conditions, the ion tracks in Sr(HE)O 3 exhibit significantly reduced diameters and a markedly different interfacial structure. Notably, the crystalline−amorphous interface in Sr(HE)O 3 shows minimal lattice distortion, confined to approximately 2−3 monolayers, in contrast to the extended disordered shell commonly observed in SrTiO 3 . Using in situ atomic-resolution electron microscopy, we further demonstrate that the amorphous/disordered regions within Sr(HE)O 3 ion tracks remain highly stable under electron irradiation, whereas tracks in SrTiO 3 readily recrystallize. This enhanced stability is attributed to the dominance of structural and chemical complexity arising from multiple B-site cations, which suppress defect migration and templated recrystallization driven by electronic excitations and local heating. Overall, this study highlights how high-entropy oxide chemistry fundamentally reshapes irradiation damage evolution, offering insights into defect formation and phase stability under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Euler–Lagrange stochastic modeling of droplet breakup and impact in supersonic flight

Blunt bodied aircraft traveling supersonically in weather environments may be damaged by impacts with water droplets and other airborne particles, such as snow and ice. Prior to an impact, these particles will encounter a bow shock that causes a discontinuity in their relative velocity with the gas phase, which can lead droplets to breakup into smaller droplets. These smaller droplets are more easily diverted from colliding with the blunt body due to their significantly reduced inertia relative to the initial rain droplets. One-way coupled Euler–Lagrange simulations are used to study the dynamics of droplets approaching a blunt body in steady two dimensional and axi-symmetric flow fields using a stochastic version of the Taylor analogy breakup model for the breakup dynamics. Ultimately, the dominant mechanism determining engineering quantities of interest was observed to be a competition between breakup time and the time available for a droplet to reach the body after encountering the bow shock. At Mach numbers 2, 3, and 6, the competition between these mechanisms was the dominant factor determining the momentum transfer to the blunt body via droplet collisions, which can be well characterized by a scaling relation.

Briney, Sam (ORCID:0000000179061081)↗
DOE OSTI2022

Mechanisms of Retention and Transport of Fission Products in Virgin and Irradiated Nuclear Graphite

To fulfill the design and licensing requirements for advanced reactors, the U.S. Department of Energy (DOE) is specifically interested in targeted research to quantify the mechanisms of fission product transport and retention in graphitic grades that will likely be used in next-generation very/high temperature reactors, fluoride salt-cooled high temperature reactors and molten salt reactors. We have examined an integrated experimental-computational approach to determine the diffusivities of selected fission products – ruthenium (Ru) and silver (Ag) in five nuclear graphite grades – POCO AXF-5Q, POCO ZXF-5Q, PCEA, IG110, and NBG-18. Experimental investigation for diffusion of other fission products (iodine, cesium, and krypton) in graphite could not be carried out due to extreme evaporation due to high vapor pressure, despite several attempts incorporating various experimental schemes. Experimental investigation with irradiated graphite could not be carried out because necessary approval to transport irradiated graphite could not be obtained due to procedural reasons at ORNL and national policy change in UK. An effort was made to “produce” by ion sputtering, but it was discontinued because the outcome did not adequately represent irradiated graphite. Thin film and “cup-and-cap” methods were employed to deposit fission products, diffusion annealed were carried out in an argon atmosphere, and concentration profiles were determined by dynamic secondary ion mass spectroscopy depth profiling technique. Regardless of graphite type, anomalous diffusion behavior of Ru (i.e., lower than expected based on Arrhenius temperature dependence) at an intermediate temperature, around 600 to 700 °C, was observed. Diffusion behavior of silver on the other hand agreed well to the Arrhenius temperature dependence although the excessive evaporation of silver hampered the accurate determination of diffusivities at temperature higher than 900 °C. In general, silver was observed to diffuse faster than Ru in the temperature range compared from 500 to 900 °C, had consistently had higher pre-exponential factor corresponding to higher jump frequency. To understand the abnormal diffusional behavior of Ru in graphite in intermediate temperature range, a detailed microstructural analysis was carried out. Clustering of Ru into particles on the surface of graphite was observed around 700 °C corresponding to de-wetting and spheroidization. This change would influence the boundary condition of diffusion, although features of clustering was much smaller (<1 mm) than the SIMS raster size (120 x 120 mm). Raman spectroscopy also demonstrated that Raman peak intensity for Ru decreased significantly for samples annealed higher than 700 °C and presence of Ru in graphite altered the crystallinity of graphite. To further elucidate the abnormal diffusion behavior observed in intermediate temperature investigated, the vibrational properties of graphene, bilayer graphene and Kr atoms in bilayer graphene were investigated with molecular dynamics. Krypton, which is a fission product, was introduced between a bilayer graphene system to evaluate (i) whether the phonon dispersion curves change with high temperature in bilayer graphene (ii) whether the addition of Kr atoms between the bilayer introduces any significant change in the vibrational properties, and (iii) whether any change in the vibrational properties can be correlated to the diffusive behavior of Kr atoms. The results show that the Kr atoms indeed exhibit an anomalous diffusive behavior at 550°C. Phonon analysis shows that a few phonon frequencies are enhanced as much as 20 to 30 % with the addition of just 0.2 at.% Kr. So, this may explain the unexpected self-diffusivity of fission products.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Impacts of Different Operation Conditions and Geological Formation Characteristics on CO2 Sequestration in Citronelle Dome, Alabama

Major concerns of carbon dioxide (CO2) sequestration in subsurface formations are knowledge of the well injectivity and gas storage capacity of the formation, the CO2 pressure and saturation plume extensions during and after injection, and the risks associated with CO2 leakage and fault reactivation. Saline reservoirs are considered as one of the target formations for CO2 sequestration through structural, residual, dissolution, and mineral trapping mechanisms. The boundary condition of the saline reservoir dictates the pressure and saturation plume extension of the injected supercritical CO2 that could expand over large distances. This can lead to sources of risk, e.g., leakage and/or fault reactivation due to presence of wells, thief zones, and geological discontinuities. Therefore, there is a critical need to develop a model that describes how risk-related performance metrics (i.e., the CO2 saturation plume size, the pressure differential plume area, and the pressure differential at specific locations) vary as a function of the size of injection, time following injection, injection operations, and geologic environment. In this study, a systematic reservoir modeling studies of anthropogenic CO2 sequestration in Citronelle dome, Alabama, was performed where all relevant scenarios and conditions to address the questions of the saturation and pressure plume size in the area of review (AoR) and post-injection site care (PISC) are considered. The objective for this study was firstly to systematically simulate CO2 sequestration, i.e., saturation dynamics, and pressure behavior over a range of operational and geological conditions and to derive conclusions about the factors influencing saturation and pressure plume size, post-injection behavior, and the risk associated with them, by developing third-generation reduced order models (ROMs) for reservoir behavior. Finally, to assess the uncertainty associated with our studies, Latin Hypercube Sampling (LHS) together with an experimental design technique, i.e., Plackett–Burman design, was used. Application of Pareto charts and respond surfaces enabled us to determine the most important parameters impacting saturation and pressure plume sizes and to quantify the auto- and cross-correlation among different parameters in both history-matched and upscaled models.

58 GEOSCIENCES↗
DOE OSTI2020

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

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

A Hydro-based MCMC Analysis of SNR 0509–67.5: Revealing the Explosion Properties from Fluid Discontinuities Alone

Using Hubble Space Telescope/ACS Hα images of SNR 0509–67.5 taken ~10 yr apart, we measure the forward shock (FS) proper motions (PMs) at 231 rim locations. The average shock radius and velocity are 3.66 ± 0.036 pc and 6315 ± 310 km s –1 . Hydrodynamic simulations, recast as similarity solutions, provide models for the supernova remnant’s expansion into a uniform ambient medium. These are coupled to a Markov chain Monte Carlo (MCMC) analysis to determine explosion parameters, constrained by the FS measurements. For our baseline model, the MCMC posteriors yield an age of 315.5 ± 1.8 yr, a dynamical explosion center at 5 h 09 m 31. s 16, $-67^\circ 31^{\prime} 17\buildrel{\prime\prime}\over{.} 1$ and ambient medium densities at each azimuth ranging over 3.7–8.0 × 10 –25 g cm –3 . The age uncertainty can be an order of magnitude larger when considering other models, or subsets of the data. We detect stellar PMs corresponding to speeds in the Large Magellanic Cloud ≥ 770 km s –1 . Five stars in the remnant show measurable PMs but none are moving radially from the dynamical center. There are four stars 1."4 from the center, including three faint, previously unidentified ones. Using coronal [Fe XIV ] λ5303 emission as a proxy for the reverse shock location, we constrain the explosion energy (for a compression factor of 4) to a value of E = (1.30 ± 0.41) × 10 51 erg for the first time from shock kinematics alone. Higher compression factors (7 or more) are strongly disfavored based on multiple criteria, arguing for inefficient particle acceleration in the Balmer shocks of SNR 0509–67.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reserve and energy scarcity pricing in United States power markets: A comparative review of principles and practices

Here, errors in forecasting load and renewable-based generation in restructured power systems mean that independent system operators (ISOs) must procure sufficient operating reserves to keep the real-time operation of the system reliable and secure. But when procured reserves turn out to be insufficient in real-time due to the lack of resource capacity or ramp capability, operators often set higher prices for reserves and energy to encourage more supply, and to motivate consumers to decrease usage or shift it to other times. This procedure, which is called scarcity or shortage pricing, is a core feature of U.S. electricity markets. It is receiving increased attention from market designers and stakeholders because scarcity will become more important for spot price formation in the future with the increased penetration of zero-marginal cost renewables, and the shrinking role of fuel costs in setting prices. Scarcity pricing is implemented in various ways by different ISOs. These differences have practical implications for the level of prices and incentives for investment, operations, and demand modification. In this paper, general approaches and specific calculation procedures for reserve and energy scarcity pricing practices and calculations across the seven ISO-based U.S. power markets are reviewed and compared. A consistent terminology is used to facilitate the comparison. Current scarcity pricing practices are grouped into three approaches: (1) imposing an adder after the spot market is run; (2) including stepwise demand curves within market clearing procedures for non-contingency reserve products (e.g., the novel flexiramp product), which tends to yield longer right tails for energy scarcity premium curves; and (3) having stepwise demand curves for traditional contingency reserve products only, which results in shorter right tails in energy scarcity curves. A generic numerical example is presented to highlight the large practical differences among the reserve scarcity pricing approaches and specific implementations. To further investigate factors that contribute the most to demand curves differences among ISOs, a sensitivity analysis is performed. This analysis shows that the largest source of differences among the curves is the scarcity prices assumed in the case of severe scarcity, while the number of steps used and whether flexiramp is considered also yields important differences in scarcity prices. As renewable penetration increases, it will become increasingly crucial to employ administrative demand curves so that spot prices more effectively motivate supply and demand adjustments exactly when and where they are needed. This study shows that the different assumptions yield very different scarcity premiums for reserves and energy, and are likely to provide divergent incentives for resources to respond to shortages. It is concluded that to promote market efficiency, a reserve shortage demand curve should have at least three features: inclusion of the marginal value of reserve products at each shortage level, consideration of the magnitude and probability of supply contingencies, and avoidance of abrupt price discontinuities that can cause excessively volatile market outcomes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Microstructural, Corrosion and Mechanical Characterization of Friction-Stir Welded Joints Between Aluminum and Magnesium Alloys

Friction Stir Welding (FSW) of Aluminum-Magnesium (Al-Mg) alloys has garnered increasing attention due to its wide-ranging applications in the automotive and aerospace industries. The microstructural, corrosion, and mechanical properties of FSW Al-Mg welds under selected conditions are critical factors that must be well-investigated to ensure their durability and performance. In this study, electrochemical tests, including Open Circuit Potential (OCP), Linear Polarization Resistance (LPR), Potentiodynamic Polarization Scan (PD), and Electrochemical Impedance Spectroscopy (EIS), have been used to analyze the corrosion resistance of these welds, along with weight loss measurements through cyclic corrosion testing (CCT) and surface evolution via SEM over an 8-week period. Mechanical properties are evaluated through microindentation, nanoindentation, fatigue, and lap-shear tests. The results of the electrochemical tests reveal that the corrosion properties of the center weld zone differ from those of the two base alloys. Indentation tests indicate that the weld zone behaves differently from the base alloys. Surprisingly, the lap shear strength did not decrease significantly even after 8 weeks of cyclic corrosion testing, indicating that the joint prevented any ingress of the electrolyte. Fatigue testing showed that applying an adhesive between the Al and Mg joint improved the fatigue life to some extent, although the corrosion resistance reduced due to the discontinuous bonding of the adhesive layer. Overall, these findings highlight the importance of investigating the corrosion and mechanical properties of FSW Al-Mg welds and provide valuable insights for future modeling works.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The Effect of Homogenization Heat Treatment on the Texture Evolution in U-10Mo Alloy

Uranium alloyed with 10 wt% Mo (U-10Mo) undergoes discontinuous precipitation (DP) when subjected to sub-eutectoid heat-treatment. It is hypothesized that the crystallographic texture of the U-10Mo can be modified to minimize the degree of DP, however the evolution of texture throughout the manufacturing process is not well understood. Crystallographic textures of two U-10Mo samples homogenized at 900°C and 1000°C were studied after uniaxial rolling and annealing. A total of 6 samples were analyzed upon rolling using a combination of hot and cold rolling conditions to strains of 83%, 90% and 97% and subsequently annealed at 700°C for one hour. Microstructures and textures were analyzed via electron backscatter diffraction (EBSD) for each homogenization, strain and post-annealed states. Results show that the U-10Mo exhibits typical rolling textures seen in body-centered cubic materials with increasing intensities of both a- and ?-fibers with increasing strain. After recrystallization, the texture appeared to be effectively randomized in all samples. Homogenization temperature showed no effect on the rolling and recrystallization behavior. Potential factors leading to the randomization of texture are discussed along with the implications of textures on the minimization of DP.

nuclear fuel, texture, rolling↗
DOE OSTI2022

Development and performance of INCONEL® alloy 740H® seam-welded piping

INCONEL® alloy 740H® is an age-hardenable nickel-based alloy approved for pressure vessels and piping within ASME Section I, VIII, and B31.1. Currently, the code applies a weld strength reduction factor (WSRF) of 0.7 to the allowable stresses for longitudinal seam welds in the time dependent creep regime. In this work, a full-scale seam weld was successfully produced using typical industrial practices. The component was solution heat-treated and aged after fabrication to improve the performance of the weldment. Tensile, bend, impact, and low-cycle fatigue tests showed the component met all the specification requirements and elevated temperature properties were within base metal expectations with failures predominately in the weld metal. Long-term creep tests, including large samples more representative of the entire weldment, were fabricated and tested to times in excess of 10,000 h at multiple temperatures. Analysis of the creep data supports a WSRF of 0.9 for the solution annealed + aged weldments in contrast to the WSRF of 0.7 currently applied to welded + aged weldments. Metallurgical analysis shows that the solution annealing causes recrystallization of the weld metal and reduces the chemical compositional and microstructural gradients in the weldments resulting in acceptable short-term performance. In long-term creep, evidence for coarsened zones in the weldments due to discontinuous coarsening reactions were identified as the mechanism leading to accelerated creep damage formation in the weld metal and sample failure.

14 SOLAR ENERGY↗
DOE OSTI2024

Properties of a fading AGN from SDSS-IV MaNGA

ABSTRACT We identify a fading AGN SDSS J220141.64+115124.3 from the internal Product Launch-11 (MPL-11) in Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The central region with a projected radius of $\sim$2.4 kpc is characterized as LINER-like line ratios while the outskirts extended to $\sim$15 kpc show Seyfert-like line ratios. The ${\rm{[O {\small III}]}}$$\lambda$5007 luminosity of the Seyfert regions is a factor of 37 (2) higher than the LINER regions without (with) dust attenuation correction, suggesting that the AGN activity decreases at least $\sim$8 $\times$ 10$^3$ yr ($\sim$2.4 kpc/light-speed) ago. We model the emission line spectra in the central region with double Gaussian components (a narrow core and a broad wing) and analyse the properties of each component. The narrow core component mostly co-rotates with the stellar disc, whereas the broad wing component with a median of the velocity dispersion $\sim$300 km s$^{-1}$ is related to a wind outflow. The kinematic position angle (PA) of the ionized gas shows an $\sim 20^{\circ }$ twist from the galaxy centre to 1.5 effective radius. The median of the PA difference between the gas and stellar components is as large as $\sim 50^{\circ }$ within 0.4 effective radius. The tidal feature in DESI image and star–gas misalignment suggest this galaxy is a merger remnant. Combining all these observational results as well as public available X-ray and MIR luminosities, we confirm this is a fading AGN, the merger process kick-started the central engine to quasar phase which ionized gas composed of tidal debris, and now the activity of the central black hole decreases. The discontinuity in ${\rm{[O {\small III}]}}$$\lambda$5007 flux and EQW maps is due to multiple AGN outbursts triggered by merger remnant gas inflows.

Mo, Hao (ORCID:0000000234430768)↗
DOE OSTI2020

Hysteretic temperature sensitivity of wetland CH4 fluxes explained by substrate availability and microbial activity: Model Archive

This Modeling Archive is in support of an NGEE Arctic publication "Hysteretic temperature sensitivity of wetland CH4 fluxes explained by substrate availability and microbial activity" in the Journal Biogeosciences (https://doi.org/10.5194/bg-17-5849-2020), which includes the model data used in the publication. CH4 emissions from terrestrial systems are posited to increase, which can offset mitigation efforts and accelerate climate change. Yet, the accuracy of modeled CH4 emissions is sensitive to the prescribed CH4 production (or emission) temperature dependencies that are currently uncertain. Here, we use a comprehensive biogeochemistry model (ecosys) to investigate factors modulating CH4 production and emission rates across a permafrost thaw gradient encompassing a partly thawed bog and a fully thawed fen. We find that seasonally varying substrate availability drives lower and higher modeled methanogen biomass and activity, and thereby CH4 production, during the earlier and later periods of the thawed season, respectively. Package follows the Model data archiving guidelines with data in a *.zip file with three subfolders containing *.csv files and raw model output files; a data dictionary table (data_dictionary.csv) and two model output description tables (ecosys_plantspecies_output_notes.csv and ecosys_soil_ouput_notes.csv) to explain the format and meaning of individual output variables; and a user guide as a *.pdf. A detailed model description can be found in the supplement of (Grant, 2013). The ecosys source code is available at https://doi:10.5281/zenodo.3906642. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Performance of reanalysis and mesoscale models off the coast of Hawai'i

The eastern Hawai'i coast in the United States is characterized by considerable wind resource fuelled by persistent trade winds, making it an important area for energy research. The need is strong for reanalyses and higher-resolution regional simulations where observations have been historically limited, such as Hawai'i's offshore environments. However, studies using offshore observations in other parts of the world have shown that significant errors can occur in reanalyses and wind datasets, which can lead to inaccurate estimates of wind energy generation, payback periods, and extreme weather risks at project locations. The degree of such errors is influenced by a number of factors, including spatial resolution and the handling of processes within the planetary boundary layer (PBL). In this work, we provide a wind resource characterization from year-long lidar buoy measurements off the eastern coast of O'ahu, Hawai'i, an environment previously unobserved at the rotor level, and use the characterization to evaluate the performance of two simulation datasets. The O'ahu deployment location is meteorologically unique and less complex than land-based wind resource characterizations, being strongly characterized by trade winds with minimal land–atmosphere interaction influences. Despite the unique and fairly consistent meteorological conditions, we hypothesize that distinct simulation datasets will exhibit diverse ranges of errors similar to those that have been seen for other offshore locations. We find the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) to strongly underestimate observed wind speeds at the O'ahu location (bias = −1.54 m s −1 at a height of 140 m above sea level), while a regional Weather Research and Forecasting Model (WRF) simulation produced by the University of Hawai'i (UH-WRF) provides a significantly smaller wind speed bias (−0.25 m s −1 ), highlighting the value of running regional, higher-resolution simulations. The large bias noted for ERA5 is driven by significant underestimation of fast wind speeds (>9 m s −1 ), which the study site is largely characterized by, along with discontinuities in the ERA5 diurnal cycle. We also speculate that the relative sparsity of observations for data assimilation in this remote part of the world could influence the performance of ERA5 and that challenges with characterizing island effects could impact the performance of both datasets.

17 WIND ENERGY↗
DOE OSTI2021

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Recent Thermal and Hydraulic Stimulation Results at Raft River, ID EGS Site

The Raft River geothermal system, located 160 km northwest of Salt Lake City, Utah is the site of a Department of Energy Enhanced Geothermal System (EGS) stimulation demonstration project. The target well RRG-9 ST-1, is roughly 1.6 km (1 mile) south of the power plant. This well has undergone a series of thermal and hydraulic stimulation treatments. Furthermore, the well passes through approximately 1,524 m MD (5,000 ft. MD) of laterally discontinuous Quaternary and Tertiary volcanoclastic and volcanic rocks that unconformably overlie the Precambrian basement. The Elba Quartzite formation located in the Precambrian basement is the target zone for the stimulation program. Fluid injected into the well exits through a fracture zone in the Elba Quartzite at 1,567 m TVD (5,140 ft. TVD). In June 2013 low-rate stimulation began at RRG-9 ST-1 using injection water from the geothermal plant. The average injection rate, temperature, and well head pressure of the low-rate stimulation was 163 Lpm (43 gpm), 34°C, and 1,931kPa (280 psig), respectively. This injection continued through late August of 2013 when wellhead pressures were intentionally increased to approximately 5,516 kPa (800 psig) at approximately 984 Lpm (260 gpm). From September 12 to September 24, 2013 cold well water was injected. The average surface injection temperature was approximately 12°C. On September 25, 2013 injection of the warmer plant water was resumed. This was continued through the rest of 2013 and into 2014. Between April 1, 2014 and April 4, 2014, a higher rate hydraulic stimulation was conducted. Colder water (~10°C) was injected at rates of 20 bpm to 30 bpm (3,178 Lpm to 4,770 Lpm) at wellhead pressures of 5,861 kPa (850 psig) over approximately 24 hours. Microseismic activity was recorded near the Tertiary-Precambrian contact. Subsequently, injection of water from the plant was resumed. As of mid-January, 2015, over 666 million liters (176 million gallons) have been injected into the well. The fate of this water is the subject of ongoing numerical analyses and continuous site monitoring. Various analyses have been implemented to monitor and assess the ongoing stimulation program. Modified Hall plot analysis shows increasing conductivity and/or a decreasing skin factor around the wellbore. Since April, 2014, the injectivity index has increased significantly from 0.28 Lpm/kPa (0.5 gpm/psig) to 0.93 Lpm/kPa (1.7 gpm/psig).

enhanced geothermal system↗
DOE OSTI2020

Carbon Fiber Prepreg Recycling- Automated Perform Manufacturing Equipment (Final Technical Report)

Discontinuous fiber composites have a compelling business case as related to light metals such as aluminum and magnesium machined or die cast components. DFC’s have been in commercial production using carbon fiber/epoxy prepreg for both commercial aircraft and for military applications for a number of years and their usage is expanding into automotive and consumer goods components. Reasons for this include significantly lower cost as compared to continuous fiber composites easier design and understanding of the material by metal traditionalists, higher degree of part complexity with relatively simple and low cost molds when compared to high-pressure die casting, and potential for near zero waste. IACMI’s 10-year goals of reduction in carbon fiber reinforced plastics manufacturing cost by 25%, reduction in their embodied energy by 50%, and recyclability of 95% can all be advanced through robust recycling. With recycling We have shown that a 15% contribution in manufacturing cost reduction is feasible, that a significant reduction in embodied energy (over 90% for the recycled fraction) is achievable, and that when combined with primary pyrolysis recycling, the 95% target can be met. Large-scale application of advanced composites began with the aerospace industry, and CRTC’s location in Washington State provided access to a waste stream of nearly 900 tonnes/year that was going to landfill. As the aerospace industry has automated parts production, their material forms have become significantly easier to recycle, and some evidence of this exists for the future automotive and wind turbine components as well. CRTC decided to focus on re-purposing the carbon fiber/epoxy aerospace scrap streams through materials reformatting and reuse, and not on primary recycling (pyrolysis mainly) as this technology is commercially available. The main reasons for this were that it offers a route to effective use of the most widely available and consistent feedstock called edge-trim, as well as allows creation of complex and low-cost parts using rapid and cost-effective compression molding technologies. The reasons also include leading the way for potential automotive pre-preg recycling as the combination of roughly 85% virgin pre-preg and 15% recycled pre-preg directly supports IACMI’s goals and enables significant part cost reduction and a path to zero-waste composite manufacturing at the OEMs or Tier 1’s. The purpose of the project was to develop equipment and technology to transform continuous fiber pre-preg in roll forms into smaller chip-format DFC feedstock, and to determine impacts on moldability and performance based on lessons learned along the path. The availability of scrap from pre-preg manufacturers was a significant value to the program even though the aerospace pre-preg is designed to be autoclave molded thus has handling, tackiness, and processing factors that are not ideal for DFCs. The scope of the project was to develop a Phase 1 system that could identify and resolve some of the major hurdles, determine where to focus effort on building a commercial-scale Phase 2 system, and develop technology for handling a wide variety of incoming scrap materials to produce the most usable product. A chopping system was developed and reduced to practice via an experimental approach that enabled all forms of scrap coming from the pre-preg manufacturers and some forms of the scrap coming from OEM production lines to be very effectively re-purposed. The systems that were developed were not considered to be production robust but rather identified and resolved some of the major technical hurdles, and demonstrated the path towards effective full-scale commercialization.

36 MATERIALS SCIENCE↗