Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magma”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Lost Circulation Materials in Geothermal Drilling: Thermal Degradation, Viscosity, Compression, and Flow Loop Tests

This dataset provides a set of experimental data on the performance of lost circulation materials (LCMs) used in geothermal drilling operations. It includes results from four distinct tests: thermal degradation, viscosity measurements, compression tests, and high-temperature flow loop experiments. The materials tested, such as cedar fiber, cotton seed hulls, magma fiber, microcellulose, and sawdust, were analyzed under a range of conditions to assess their effectiveness in mitigating fluid losses in geothermal wells. Data on thermal degradation explores mass loss, water loss, and gas release for each material when subjected to temperatures between 90C and 250C. Viscosity measurements capture the behavior of LCM mixtures under different temperatures, while compression tests investigate the structural integrity of both undegraded and thermally degraded materials. The high-temperature flow loop experiments simulate geothermal well conditions, analyzing how different materials and mixtures affect fluid flow and sealing performance in a fractured reservoir environment. In addition to the experimental data, the dataset includes photos documenting the materials, schematics of test setups, and several published manuscripts that provide detailed insights into the experimental methods and findings.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗
DOE OSTI2024

Sulfur speciation in dacitic melts using X-ray absorption near-edge structure spectroscopy of the S K -edge (S-XANES): Consideration of radiation-induced changes and the implications for sulfur in natural arc systems

Abstract The synchrotron technique of micro X-ray absorption near-edge structure spectroscopy at the sulfur K-edge (S-XANES) provides a unique opportunity to measure the proportion of different oxidation states of sulfur (S) in silicate glasses. Although applied extensively in the analysis of basaltic silicate glasses, few S-XANES studies have investigated variations in S oxidation states with fO2 in felsic silicate glasses. In addition, no study has systematically compared the S-XANES results obtained from the same samples at different photon flux densities to quantify the relationship between exposure time and changes in S speciation in silicate glass, as has been done for Fe and V. This study evaluates observed differences in S speciation measured in experimentally produced H2O-saturated dacitic glasses over a range of reducing to oxidizing conditions (from log fO2 = ΔFMQ-0.7 to ΔFMQ+3.3; FMQ is the fayalite-magnetite-quartz mineral redox buffer) and equilibrated at 1000 °C and 300 MPa. S-XANES spectra were collected at three different photon flux densities using three microspectroscopy beamlines. As is observed in S-XANES analyses of basaltic silicate glasses, beam-induced changes to the S6+/ΣS are observed as a function of photon flux density and beam exposure time. Our results demonstrate that silicate glasses of dacitic composition undergo beam-induced photo-reduction in samples equilibrated at ΔFMQ > +1.75 and photo-oxidation if equilibrated at ΔFMQ < +1. The time required to observe beam-induced changes in the spectra varies as a function of flux density, and our study establishes an upper photon density limit at ~1.0 × 1012 photons/µm2. The S6+/ΣS calculated from spectra collected below this absorbed photon limit at intermediate flux densities (~1–4 × 109 photons/s per µm2) are affected by beam damage, as no conditions were found to be completely free of beam-induced changes. However, the S6+/ΣS ratios calculated below the limit at intermediate flux densities are consistent with thermodynamic constraints, demonstrating that S6+/ΣS ratios calculated from S-XANES spectra can be considered reliable for estimating the oxygen fugacity. Our results carry important implications for the S budget of felsic magmas and dissolution mechanisms in evolved melts. While our results from all three flux densities show the presence of S4+ dissolved in relatively oxidized (ΔFMQ > +1.75) dacitic glass, even in the spectra exposed to the lowest photon densities, we are unable to rule out the possibility that the S4+ signal is the result of instantaneous X-ray irradiation induced beam damage using S-XANES alone. When our spectra are compared to S-XANES spectra from basaltic silicate glasses, important differences exist in the solubility of S2– and S6+ between dacitic silicate melts, pointing to differences in solubility mechanisms as melt composition changes. This study highlights the need for further investigation into beam damage systematics, presence of S4+, and the solubility mechanisms of different oxidation states of S as silicate melt composition changes.

Geochemistry & Geophysics↗
DOE OSTI2020

Towards Use of Mixed Precision in ECP Math Libraries [Exascale Computing Project]

The use of multiple types of precision in mathematical software has the potential to increase its performance on new heterogeneous architectures. The xSDK project focuses both on the investigation and development of multiprecision algorithms as well as their inclusion into xSDK member libraries. This report summarizes current efforts on including and/or using mixed precision capabilities in the math libraries Ginkgo, heFFTe, hypre, MAGMA, PETSc/TAO, SLATE, SuperLU, and Trilinos, including KokkosKernels. It contains both numerical results from libraries that already provide mixed precision capabilities, as well as descriptions of the strategies to incorporate multiprecision into established libraries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Towards Use of Mixed Precision in ECP Math Libraries

The use of multiple types of precision in mathematical software has the potential to increase its performance on new heterogeneous architectures. The xSDK project focuses both on the investigation and development of multiprecision algorithms as well as their inclusion into xSDK member libraries. This report summarizes current efforts on including and/or using mixed precision capabilities in the math libraries Ginkgo, heFFTe, hypre, MAGMA, PETSc/TAO, SLATE, SuperLU, and Trilinos, including KokkosKernels. It contains both numerical results from libraries that already provide mixed precision capabilities, as well as descriptions of the strategies to incorporate multiprecision into established libraries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Play Fairway Analysis of the Snake River Plain, Idaho (Final Report)

The Snake River Plain (SRP) volcanic province overlies a thermal anomaly that extends deep into the mantle; it represents one of the highest heat flow provinces in North America. The Yellowstone hotspot continues to feed a magma system that underlies much of southern Idaho and has produced basaltic volcanism as young as 2000 years old. It has been estimated to host up to 855 MW of potential geothermal power production, most of which is associated with the Snake River Plain volcanic province in Idaho, which lies outside the area of Yellowstone National Park (Neely, K.W. and Galinato, G., 2007, Geothermal power generation in Idaho: an overview of current developments and future potential, Open File Report, Idaho Office of Energy Resources).

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Batched Sparse Linear Algebra (Final Report for Subcontract B648960)

This report finalizes design specifications for developing batched kernels for small tensor operations for unassembled matrix-free iterative solvers, batched solvers for partially assembled operators, and batched solvers with support for various sparse formats. The outcome of the project milestones is a set of interfaces to Batched Sparse LA solvers running on hardware accelerators for use in ECP Libraries and Applications. It is part of the development of sparse batched kernels, solvers/preconditioners as well as creating interoperability in xSDK libraries with sparse and dense batched functions to benefit ECP applications. The participants included representatives from ECP libraries (not limited to the xSDK project), applications, and vendors (AMD, Intel, and NVIDIA). Batched sparse linear algebra solvers form the new frontier for algorithmic development and performance engineering. Many applications (ECP and non-ECP alike) require simultaneous solutions of small linear systems of equations that are structurally sparse. To move towards high hardware utilization, it is important to provide these applications with appropriate interfaces to efficient batched sparse solvers running on modern hardware accelerators. We present interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the software portable between the major hardware accelerators from AMD, Intel, and NVIDIA. The presented interface specifications includes batched band, sparse iterative, and sparse direct solvers. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, SUNDIALS, and SuperLU_dist.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Milestone 49 Report: Batched Sparse LA Phase 5 Implementation

Batched sparse linear algebra operations in general, and solvers in particular, have become the major algorithmic development activity and foremost performance engineering effort in the numerical software libraries work on modern hardware with accelerators such as GPUs. Many applications, ECP and non-ECP alike, require simultaneous solutions of many small linear systems of equations that are structurally sparse in one form or another. In order to move towards high hardware utilization levels, it is important to provide these applications with appropriate interface designs to be both functionally efficient and performance portable and give full access to the appropriate batched sparse solvers running on modern hardware accelerators prevalent across DOE supercomputing sites since the inception of ECP. To this end, we present here a summary of recent advances on the interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the corresponding software portable between the major hardware accelerators from AMD, Intel, and NVIDIA, while maintaining the appropriate disclosure levels conforming to the active NDA agreements. The presented interface specifications include a mix of batched band, sparse iterative, and sparse direct solvers with their accompanying functionality that is already required by the application codes or we anticipated to be needed in the near future. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, PETSc, and SuperLU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
DOE OSTI2020

Structure and Stability of Iron Fluoride at High Pressure–Temperature and Implication for a New Reservoir of Fluorine in the Deep Earth

Fluorine (F) is the most abundant halogen in the bulk silicate Earth. F plays an important role in geochemical and biological systems, but its abundance and distribution in the terrestrial mantle are still unclear. Recent studies suggested that F reservoirs in the deep mantle are potentially hosted in terrestrial oxide minerals, especially in aluminous bridgmanite. However, the knowledge about the formation and stability field of fluoride in the Earth’s interior is rare. In this study, we combine in situ laser-heated diamond anvil cell, synchrotron X-ray diffraction, and first-principles structure search to show that a new tetragonal structure of FeF 3 is stable at pressures of 78–130 GPa and temperatures up to ~1900 K. Simulation predicted the tetragonal phase takes a much denser structure due to the rotation of FeF 6 octahedral units. The equations of states of tetragonal FeF3 are determined by experiment and verified by simulation. Our results indicate that FeF 3 can be a potential key phase for storing F in the Earth’s lower mantle and may explain some mantle-derived magma with high F concentration.

58 GEOSCIENCES↗
DOE OSTI2021

Composition and Pressure Effects on Partitioning of Ferrous Iron in Iron-Rich Lower Mantle Heterogeneities

Both seismic observations of dense low shear velocity regions and models of magma ocean crystallization and mantle dynamics support enrichment of iron in Earth’s lowermost mantle. Physical properties of iron-rich lower mantle heterogeneities in the modern Earth depend on distribution of iron between coexisting lower mantle phases (Mg,Fe)O magnesiowüstite, (Mg,Fe)SiO3 bridgmanite, and (Mg,Fe)SiO3 post-perovskite. The partitioning of iron between these phases was investigated in synthetic ferrous-iron-rich olivine compositions (Mg0.55Fe0.45)2SiO4 and (Mg0.28Fe0.72)2SiO4 at lower mantle conditions ranging from 33–128 GPa and 1900–3000 K in the laser-heated diamond anvil cell. The resulting phase assemblages were characterized by a combination of in situ X-ray diffraction and ex situ transmission electron microscopy. The exchange coefficient between bridgmanite and magnesiowüstite decreases with pressure and bulk Fe# and increases with temperature. Thermodynamic modeling determines that incorporation and partitioning of iron in bridgmanite are explained well by excess volume associated with Mg-Fe exchange. Partitioning results are used to model compositions and densities of mantle phase assemblages as a function of pressure, FeO-content and SiO2-content. Unlike average mantle compositions, iron-rich compositions in the mantle exhibit negative dependence of density on SiO2-content at all mantle depths, an important finding for interpretation of deep lower mantle structures.

Dorfman, Susannah M. (ORCID:0000000239689592)↗
DOE OSTI2022

Reduced Late Bombardment on Rocky Exoplanets around M Dwarfs

Ocean-vaporizing impacts of chemically reduced planetesimals onto the early Earth have been suggested to catalyze atmospheric production of reduced nitrogen compounds and trigger prebiotic synthesis despite an oxidized lithosphere. While geochemical evidence supports a dry, highly reduced late veneer on Earth, the composition of late-impacting debris around lower-mass stars is subject to variable volatile loss as a result of their hosts' extended pre-main-sequence phase. We perform simulations of late-stage planet formation across the M-dwarf mass spectrum to derive upper limits on reducing bombardment epochs in Hadean-analog environments. We contrast the solar system scenario with varying initial volatile distributions due to extended primordial runaway greenhouse phases on protoplanets and the desiccation of smaller planetesimals by internal radiogenic heating. We find a decreasing rate of late-accreting reducing impacts with decreasing stellar mass. Young planets around stars ≤0.4 M ⊙ experience no impacts of sufficient mass to generate prebiotically relevant concentrations of reduced atmospheric compounds once their stars have reached the main sequence. For M-dwarf planets to not exceed Earth-like concentrations of volatiles, both planetesimals, and larger protoplanets must undergo extensive devolatilization processes and can typically emerge from long-lived magma ocean phases with sufficient atmophile content to outgas secondary atmospheres. Our results suggest that transiently reducing surface conditions on young rocky exoplanets are favored around FGK stellar types relative to M dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Iodine speciation in basaltic melts at depth

The speciation of iodine in basalts has been investigated by combining in situ X-ray diffraction at high pressures and temperatures up to 4.9 GPa and 1600 °C, and Raman spectroscopy on recovered high pressure glasses at ambient conditions. Both methods point to iodine being oxidized in basalts, whether molten or quenched as glasses. Observed interatomic distances and Raman vibrational modes are consistent with iodine being dissolved as complex iodate groups alike polyiodates or periodates, not as IO$^-_3$ groups. Iodine speciation in basalts therefore seems to reflect a trend amongst halogens, with lighter chlorine bonding to network modifying cations, and bromine changing affinity from network modifying cations to oxygen anions under pressure. In the absence of a fluid aqueous phase, iodine could thus reach the Earth’s surface in basaltic magmas as an oxide, not as a reduced species.

58 GEOSCIENCES↗
DOE OSTI2021

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗
DOE OSTI2021

Porting the Nonlinear Optimization Library HiOp to Accelerator-Based Hardware Architectures

While interior point method has been the centerpiece of nonlinear programming tools used in science and engineering, its reliance on linear solvers that can tackle sparse symmetric indefinite and highly ill-conditioned problems made it difficult to implement it effectively on hardware accelerators. HiOp optimization package attempts to provide an implementation of the interior point method suitable for hardware accelerators by compressing the original sparse problem to produce an underlying linear problem that is dense and of manageable size. Implementations of dense linear solvers are more mature and utilize hardware accelerators better than their sparse counterparts. There is a number of important domain problems, such as optimal power flow analysis for power grids, where the sparse problem can be effectively compressed and deploying dense linear solver within the interior point method can improve performance. Here we describe a portable implementation of HiOp optimization engine, which uses a linear solver from Magma library and runs entirely on hardware accelerators. To compress the problem, HiOp uses customized mixed dense-sparse linear algebra. All HiOp kernels are implemented using Umpire and RAJA portability libraries. We describe details of the implementation and discuss trade-offs between performance, portability and development cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The timing of lunar solidification and mantle overturn recorded in ferroan anorthosite 62237

Ferroan anorthosite suite (FAS) rocks are widely interpreted to represent primordial lunar crust. Despite their importance in pinpointing the timing of lunar crust formation, robust chronological investigations for this rock type are scarce. Here, we report the Ar-Ar, Rb-Sr, and Sm-Nd isotopic systematics for the FAS troctolitic anorthosite 62237. The Ar-Ar isotopic system has been reset by a thermal event at 3710 ± 48 Ma, and the Rb-Sr isotopic systematics has been disturbed such that a Rb-Sr isochron age cannot be determined. However, an internal isochron for the Sm-Nd isotopic system has yielded an age of 4350 ± 73 Ma (MSWD = 2.0) with an initial ε 143 Nd CHUR of -0.53 ± 0.26. The mineral and whole-rock fractions of 62237 plot on the same internal isochron as FAS sample 60025. The combined datasets define an age of 4372 ± 35 Ma (MSWD = 4.0) with an initial ε 143 Nd CHUR of -0.17 ± 0.22. Literature Sm-Nd data for FAS and Mg-suite whole-rocks also plot on the 60025-62237 isochron. The coherence of data from both FAS and Mg-suite rocks examined thus far suggests that both rock suites formed contemporaneously from identical, or nearly identical, sources. In addition, the concordance of FAS and Mg-suite ages suggests that primordial crust solidification either involved both magmatic suites, or that Mg-suite magmatism was contemporaneous with FAS magmatism within resolution of the Sm-Nd chronometer. The ages for FAS and Mg-suite also coincide with the formation ages of the mare basalt source regions and urKREEP. Ferroan anorthosite suite rocks and urKREEP are thought to represent primordial LMO solidification products, whereas Mg-suite and the mare basalt source regions are argued to represent mixtures of various LMO crystallization products that were formed during density-driven overturn of the LMO. The concordance of ages implies that the 4372 ± 35 Ma Sm-Nd isochron records the age of mantle overturn, and that overturn occurred during, or shortly after, solidification of the LMO.

58 GEOSCIENCES↗
DOE OSTI2026

A collective trigger for widespread planetesimal formation revealed by accretion ages

The formation of planetesimals was an integral part of the cascading series of processes that built the terrestrial planets. To illuminate planetesimal formation, here we develop a refined thermal evolution model to calculate the formation ages of meteorite parent planetesimals. This model includes chemical reactions and phase changes during heating, as well as natural variations in the proportions of the constituent phases of these planetesimals. We find that the parent bodies of non-carbonaceous (NC) and carbonaceous (CC) iron meteorites start forming at very similar times ( ∼ 0.95 Myr after calcium-aluminium-rich inclusion [CAI] formation) and occupy overlapping time windows. NC and CC chondrite parent bodies formed later during non-overlapping periods. We combine these ages with proportions of isotopic end-members we recover from mixing models to construct records of motion throughout the protoplanetary disk. These records argue that NC and CC material traversed the barrier in the disk after ∼ 0.95 Myr after CAI formation. The onset of this motion coincided with planetesimal formation, indicating that the phenomenon that drove motion also triggered planetesimal formation. We argue that this feature also served as the semi-permeable barrier in the disk. Although its identity is uncertain, the effects this phenomenon had on the timing of planetesimal formation and motion through the disk can now serve as constraints on models of disk evolution. Models that reproduce these effects would elucidate the nature and implications of this phenomenon, which is key to unlocking a holistic model of terrestrial planet building.

Astronomy and AstroPhysics↗