SEARCH · Search NASA

Results for “Explainability”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Comets and carbonaceous chondrites delivered noble gases to the Moon

Trapped xenon isotopes in two Apollo 17 mare basalt fragments are similar to those in primitive meteorites. Xenon would have been effectively excluded from the Moon in the aftermath of its formation in a giant impact. As such, lunar mantle xenon trapped in the mare basalts is best explained as being derived from late accretion occurring before or concurrent with the formation of the lunar crust. Here, the isotopic composition indicates that a combination of comets and carbonaceous chondrites delivered this xenon to the lunar mantle. The inferred mass of accreted cometary ice would have delivered significantly less than a part per million of water to the lunar mantle. The inferred mass of accreted carbonaceous chondrites would have supplied at least a half of a part per million of water. The data further indicate that enstatite chondrites are unlikely to have supplied the majority of late accreted mass.

Chemistry↗

DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗

DOE OSTI2024

What to expect when you're expecting engagement: Delivering procedural justice in large-scale solar energy deployment

Community engagement in the planning process to build large-scale solar (LSS) projects can win local support and advance procedural justice. However, an understanding of community engagement in current LSS development is lacking. Using responses from a U.S. nationwide survey (n = 979) of residential neighbors living within 3 miles (4.8 km) of completed LSS projects (i.e. “solar neighbors”) and project details from the U.S. Large-Scale Solar Photovoltaic Database (USPVDB), this study seeks to answer the following questions: How are solar neighbors' perceptions of community engagement associated with their attitudes toward their LSS projects? How do solar neighbors' perceptions of community engagement compare to their expectations? And, how do neighbors explain what they perceived about the planning process? We answer these questions using mixed methods, including regression modeling, a new gap analysis technique, and qualitative coding. We find that higher perceived engagement is associated with more positive attitudes toward the project, even when controlling for respondents who acted in opposition. Supporters and opponents alike expect more engagement than they perceived and information about projects both before construction and after operation is lacking. Awareness and engagement expectations increase at certain project size and proximity thresholds. However, most neighbors expect the public to offer input during engagement, but not make decisions. We contextualize these findings with explanatory comments from respondents.

14 SOLAR ENERGY↗

DOE OSTI2025

Go slow to go fast? A review of the impacts of permitting on large-scale solar project development

State and local permitting challenges could impede the ability of large-scale solar (LSS) to meet growing electricity demand in the United States. Here, we review research that explores LSS permitting and its impacts on the pace and scale of LSS project development. Research on LSS permitting is relatively scarce, such that we support our review with research in the context of wind permitting, where appropriate. Further, few studies attempt to rigorously quantify the effects of permitting on the pace and scale of LSS project development. The available evidence allows us to identify various hypotheses and identify gaps for further research. Our review suggests that differences in permitting policies, regulations, and ordinances explain relatively little variation in LSS permitting and development outcomes across jurisdictions, except where jurisdictions implement rules designed to impede LSS. The evidence suggests LSS permitting challenges largely accrue during the implementation of permitting processes. Recent research suggests that community opposition to project development is a key driver of LSS permitting challenges, given that project opponents often use permitting processes to translate opposition into legal action. We call on future researchers to more concretely describe the LSS permitting challenge and to identify the specific actors responsible for implementing solutions.

Community opposition↗

DOE OSTI2026

Governing in Time: Temporal Capacity and the Feasibility of Energy Transitions

Energy systems function as both technological systems and temporal institutions that shape how societies coordinate, justify, and support collective choices over time. This paper introduces the concept of governance horizons to explain why energy transitions can remain morally supported yet become institutionally weak under increasing pressure. We argue that governability depends on institutions' capacity to synchronize across multiple timeframes - aligning short-term decisions with intermediate coordination and long-term commitments. When this synchronization fails, transitions struggle not because their goals are dismissed, but because governance lacks sufficient time to justify, coordinate, and uphold decisions. Comparative analysis of San Antonio, Texas, and Interior Alaska reveals how energy system pressures generate distinct temporal configurations: San Antonio exhibits governance horizon stretching, where institutions must simultaneously meet near-term reliability demands and long-term transformation goals, while Interior Alaska exhibits horizon compression, where extreme environmental constraints force decision-making into short stabilization cycles. In both contexts, public support for sustainability goals coexists with institutional strain because evaluative judgments are unevenly distributed over time. A temporal configuration analysis is introduced as a diagnostic analytic stance for identifying these patterns. By treating temporal alignment as an explanatory variable rather than a background condition, this approach clarifies how feasibility, sequencing, and legitimacy are shaped by constraints on institutional time. The analysis demonstrates that successful energy transitions depend not only on technological innovation or institutional support, but on governance systems’ ability to sustain credible coordination across multiple time horizons.

Comparative case study↗

DOE OSTI2024

Understanding particle size effect on fast-charging behavior of graphite anode using ultra-thin-layer electrodes

Extreme fast charging (≤ 15 min) of lithium-ion batteries is highly desirable to accelerate mass-market adoption of electric vehicles. However, significant capacity fading, as well as safety issues due to the lithium plating caused by the fast charging rate, limit its implementation. In this study, we investigated the fast-charging capability of graphite materials with various particle sizes. To eliminate the Li + ion concentration gradient effect across the thickness of the electrode, ultra-thin-layer graphite electrodes were developed to investigate the "real" fast-charging capability of graphite at the particle level. Electrochemical assessments as well as microscopic characterizations revealed that smaller particles exhibited superior fast-charging performance, featuring enhanced capacity reversibility, faster charging rate, and less lithium plating under the same fast-charging conditions. It is shown that small-particle graphite (mean radius of 3.3 μm) could withstand a 4C charge (to 80 % state-of-charge) without plating, with minimal plating occurring at 6C. Thicker particles exhibited plating at lower C-rates. Since the experimental data could not directly explain whether intra-particle diffusion limitations or interfacial reaction limitations dominated the plating mechanism, the pseudo-2-dimensional model was used to evaluate the most likely plating mechanism. The model suggested that particle-level diffusion is the dominant mechanism contributing to plating at high rates. Finally, this work provides comprehensive insights into the particle size effects on fast-charging capability, offering a better understanding of fast-charging behavior and valuable guidance for designing optimal electrode architecture for high-rate lithium-ion batteries.

25 ENERGY STORAGE↗

DOE OSTI2025

Mechanisms and stability of Li dynamics in amorphous Li-Ti-P-S-based mixed ionic–electronic conductors: A machine learning molecular dynamics study

Mixed ionic–electronic conductors (MIECs) exhibit both high ionic and electronic conductivity to improve the battery performance. In this work, we investigate the mechanism and stability of transport channels in our recently developed MIEC material, amorphous Ti-doped lithium phosphorus sulfide (LPS), using molecular dynamics (MD) simulations with a 99% accurate machine-learning force field (MLFF) trained on ab initio MD data. The achieved MLFF helps efficient large-scale MD simulations on LPS with three Ti concentrations (10%, 20%, and 30%) and six temperatures (25°C to 225°C) to calculate ionic conductivity, activation energy, Li-ion transport mechanism, and configurational entropy. Results show that ionic conductivities and activation energies are consistent with our recent experimental values. Moreover, Li-ion transport occurs via free-volume diffusion facilitated by the formation of disordered Li-S polyhedra. The enhanced stability of transport channels at 10% and 20% Ti doping, compared to 0% and 30%, is observed by analyzing the vibrational and configurational entropy of these disordered Li-S polyhedra. Overall, this study highlights the utility of MLFF-based large-scale MD simulations in explaining the transport mechanism and the stability of Li-ion in Ti-doped LPS electrolyte with significant computational efficiency.

And configuration entropy↗

DOE OSTI2024

Computational analysis of flame initiation, quenching, and re-ignition in a prechamber natural gas engine under varying EGR-dilution levels

The on-road natural-gas (NG) fueled transportation relies on stoichiometric spark-ignition engines for the advantages of simple after-treatment system despite the efficiency penalty relative to lean-burn combustion strategies. Exhaust gas recirculation (EGR) has the potential to reduce this efficiency gap at low to moderate loads without the need for complex lean-exhaust aftertreatment systems. However, EGR dilution leads to reduced combustion stability and increased cycle-to-cycle variability. A promising technology that has the potential to achieve reliable operation under diluted conditions is the prechamber ignition (or turbulent jet ignition) which uses chemically active turbulent jets generated from combustion inside a prechamber to initiate, stabilize and accelerate combustion of the mixture inside the main chamber. The present work focusses on developing a RANS-based CFD approach to accurately reproduce in-cylinder phenomena in a stoichiometric NG prechamber-assisted heavy-duty engine without relying on complex combustion models that account for turbulence-chemistry interactions. This is necessary because reactive prechamber jets at high EGR dilution tend to extinguish while emerging into the main chamber, which is followed by a phase of re-ignition — a phenomenon that conventional G-equation or well-stirred reactor combustion models cannot reproduce. With addition of a damping multiplier to the well-stirred reactor model, the predictions are seen to show good agreement with experimental pressure evolution and combustion images acquired from a single cylinder Cummins N-14 optical diesel engine retrofitted with a prechamber ignition system. Model predictions of local heat release in the flame and temperature evolution inside the flame are used to investigate combustion dynamics in the prechamber and the main chamber. It is seen that the well-stirred reactor model with the inclusion of damping is able to reproduce the temporary reduction in heat release within the flame, which can be considered equivalent to quenching of jets, and the subsequent re-ignition of the flame inside the main chamber. The delay between quenching and re-ignition depends on the amount of dilution, as explained by an illustration of flame evolution in a Borghi diagram.

Prechamber ignition↗

DOE OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A survey on checkpointing strategies: Should we always checkpoint à la Young/Daly?

The Young/Daly formula provides an approximation of the optimal checkpointing period for a parallel application executing on a supercomputing platform. It was originally designed to handle fail-stop errors for preemptible tightly-coupled applications, but has been extended to other application and resilience frameworks. Here, we provide some background and survey various scenarios to assess the usefulness and limitations of the formula, both for preemptible applications and workflow applications represented as a graph of tasks. We also discuss scenarios with uncertainties, and extend the study to silent errors. We exhibit cases where the optimal period is of a different order than that dictated by the Young/Daly formula, and finally we explain how checkpointing can be further combined with replication.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Reassessing the proposed “CY chondrites”: Evidence for multiple meteorite types and parent bodies from Cr-Ti-H-C-N isotopes and bulk elemental compositions

Here, we report a coordinated bulk Cr-Ti-H-C-N isotopic and compositional study of six carbonaceous chondrites from Antarctica that are often considered to be related and termed Yamato-like carbonaceous (CY) chondrites. These meteorites are known to have undergone extensive aqueous alteration followed by different degrees of thermal alteration, to be similar to one another in regard to mineralogy, and share affinities with both the Ivuna-like carbonaceous (CI) and Mighei-like carbonaceous (CM) chondrites. While mineralogically similar, a key difference among these samples is that chondrules have been found in some of these samples, but not in others. The aim of this study is to evaluate the relationship of these meteorites to one another, and investigate how they relate to the CI and CM chondrite groups. We find that with the addition of the isotopic compositions of these ‘CY’ chondrites, there is now a continuum of isotopic compositions among the carbonaceous chondrites. The CI chondrites are no longer separate in O isotopic compositions from the other carbonaceous chondrite groups in plots of Cr-O and Ti-O. We also find that the ‘CY’ chondrites represent two distinct populations, which correlate with their heating stage. However, the peak temperatures experienced by each population can only explain the differences in H and C isotopes and abundances and N abundances between samples, and cannot have caused the differences in N, Cr, and Ti isotopes, or all the volatile element depletions of the sample’s bulk compositions. Instead, we conclude that the compositional and isotopic data of these two populations correlate with their known chondrule abundances, indicating distinct precursors for each population. We find it most likely that these samples originate from two distinct asteroids, implying that among the six samples studied here, there are not five related samples to constitute a new meteorite group. The chondrule-free heating Stage III samples are most likely heated CI chondrites, while the chondrule-bearing heating Stage IV samples could be heated CM chondrites.

58 GEOSCIENCES↗

DOE OSTI2025

Mn(II)-induced phase transformation of Mn(IV) oxide in seawater

Manganese (Mn) oxides are key components of oceanic and lacustrine Mn nodules and influence metal cycling through oxidation and adsorption processes. Layered Mn oxides (LMOs) are the most common minerals in these nodules and the immediate products of microbially mediated Mn(II) oxidation by O 2 . LMOs can transform into tunneled Mn oxides (TMOs), Mn oxyhydroxides (MnOOH), or Mn(II,III) phase (Mn 3 O 4 ). LMOs often concur with Mn(II) in the environment and the adsorption and oxidation of Mn(II) by LMOs can greatly promote the transformation of LMOs to those phases. However, the Mn(II)-promoted transformation of LMOs in seawater—rich in various cations (300 mM Na + , 10 mM K + , 50 mM Ca 2+ , and 10 mM Mg 2+ ) remains poorly understood. We examined the transformation of δ-MnO 2 in artificial seawater (pH 8.2) under anoxic conditions with the Mn(II)/MnO 2 ratio (r) ranging from 0.08 to 3.83. To assess the effect of ionic strength (IS), parallel experiments were conducted in a mixed 530 mM NaCl and 10 mM KCl solution (having seawater ionic strength but without Ca 2+ and Mg 2+ ) and in 100 mM NaCl solution as a control. At low r (0.08), δ-MnO 2 transformed into triclinic birnessite and a 4 × 4 TMO in 100 mM NaCl solution, which, however, was suppressed in seawater due to strong interactions of Ca 2+ /Mg 2+ with δ-MnO 2 . In the mixed 530 mM NaCl and 10 mM KCl solution (the same ionic strength as of seawater), the transformation occurred extensively but the products had lower crystallinity compared to in 100 mM NaCl solution. At the high Mn(II)/MnO 2 ratios (0.5 ≤ r ≤ 3.83), δ-MnO 2 transformed extensively into MnOOH phases and hausmannite (Mn 3 O 4 ) in 100 mM NaCl solution. The seawater suppressed the transformation, but the suppression became weaker with increasing Mn(II)/MnO 2 ratio. For example, the transformation was completely suppressed at r = 0.5 but essentially negligible at r = 3.83. The suppression at these high Mn(II)/MnO 2 ratios was mainly ascribed to the influence of Ca 2+ and Mg 2+ rather than of the high IS, and the weaker suppression at the higher Mn(II)/MnO 2 ratio suggests stronger competition of Mn(II) with Ca 2+ /Mg 2+ for interacting with δ-MnO 2 . Moreover, the composition and crystallinity of the transformation products (i.e., the relative abundance of MnOOH (α, β, and γ) and Mn 3 O 4 ) were influenced by both the high ionic strength and the presence of Ca 2+ and Mg 2+ . Therefore, even though Ca 2+ and Mg 2+ concentrations are much lower than Na + in seawater, their impacts on the transformation are dominant. Our study explains why MnOOH phases, hausmannite, and TMOs are less common than LMOs in oceanic environments, partially because seawater chemistry suppresses their formation. LMOs are the most reactive for metal adsorption and oxidation among all Mn oxides. Thus, the high stability of LMOs in an oceanic environment confers the high impacts of Mn oxides on metal cycling in the ocean.

Divalent manganese↗

DOE OSTI2025

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗

DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗

DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗

DOE OSTI2024

Delivering laser performance conditions to enable fusion ignition, and beyond at the National Ignition Facility

On December 5th, 2022, controlled fusion ignition was demonstrated for the first time at the National Ignition Facility (NIF), a major achievement in the field of Inertial Confinement Fusion (ICF) requiring a multi-decadal effort involving broad national and international collaborations. To drive the fusion ignition reaction with the compressed fuel capsule, that yielded 3.15 MJ of nuclear energy [1], the NIF laser delivered a high-precision pulse shape with 2.05 MJ of ultra-violet (UV) laser energy and a peak power of 440 TW. This laser energy was an increase of ~8 % compared to that delivered on the previous “threshold of ignition” record yield experiment (1.37 MJ of yield for 1.89 MJ of laser energy) on August 8th, 2021. Here we explain how the results of our extensive research in laser technology and UV optics damage mitigation led to major improvements in the NIF laser, enabling this energy increase along with additional accuracy, precision, and power balance enhancements. Furthermore, we will discuss on-going efforts that have enabled operations at 2.2 MJ of UV energy as well as potential new initiatives to push the laser performance –accuracy and delivered energy– to even higher levels in the future as previously demonstrated on a small subset of NIF beams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Cooper minima of atoms in a plasma environment

Recent efforts have extended the calculations of photon ionization and associated Cooper minima to atoms in a plasma environment simulated by screening the atomic potential. Furthermore, in this paper previously unrecognized features of bound electron photon excitations in screened atomic potentials are examined. First, the Cooper minima location as a function of photon energy is shown to be weakly dependent on the screening, which is explained by stationary perturbation theory. Second, the screening transfers transition probability from high to low photon energies relative to the isolated atom result.

Cooper minima↗