Search NASA⌕ Search

SEARCH · Search NASA

Results for “spherical maps”

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.

66 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
DOE OSTI2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Various Low-Symmetry Phases in High- χ and Conformationally Asymmetric PDMS- b -PTFEA Copolymers

Since the discovery of the Frank–Kasper (FK) phase a decade ago, single-component block copolymer systems have offered unique opportunities to explore the structural and phase-transition properties of topologically close-packed FK phases composed of point-particle motifs. Here, we present the formation of complex low-symmetry phases and their thermodynamic stability using high-χ and high-ε linear diblock copolymers, where χ and ε are the Flory–Huggins interaction parameter and conformational asymmetry between the two blocks, respectively. In this work, a series of polydimethylsiloxane-b-poly(2,2,2-trifluoroethyl acrylate) (PDMS-b-PTFEA) copolymers are prepared by tuning the volume fraction ($f_{\text{PTFEA}}$) of a minor PTFEA block. The packing structures of PDMS-b-PTFEAs are mapped in the space of temperature versus $f_{\text{PTFEA}}$ using X-ray scattering measurements. In addition to the packing structures reported earlier, we identify another A15 phase at higher $f_{\text{PTFEA}}$ close to a hexagonally packed cylinder (HEX) phase and establish the phase sequence of HEX–A15–σ–C14–BCC with decreasing $f_{\text{PTFEA}}$. Theoretical calculation of relative free-energy densities confirms the experimental phase sequence of the composition-dependent sphere-packing structures. Interestingly, the domain sizes of three-dimensional (3D) sphere-packing structures are nearly temperature-invariant compared to that of a two-dimensional HEX structure, and we attribute this dimensional stability to the strong enthalpic effects in the formation of 3D packing structures and the geometrical characteristics of spherical domains. Our results demonstrate that the interaction parameter χ, together with the known role of conformational asymmetry ε, serves as the primary thermodynamic parameter for selecting and stabilizing various low-symmetry packing structures in sphere-forming linear diblock copolymer systems.

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

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

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

Theoretical Assessment of the Transition Between Electron Emission Mechanisms for Nonplanar Diodes

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

Conformal mapping↗
DOE OSTI2021

Direct numerical simulation of turbulent dispersion of evaporative aerosol clouds produced by an intense expiratory event

Airborne particles are a major route for transmission of COVID-19 and many other infectious diseases. When a person talks, sings, coughs, or sneezes, nasal and throat secretions are spewed into the air. After a short initial fragmentation stage, the expelled material is mostly composed of spherical particles of different sizes. While the dynamics of the largest droplets are dominated by gravitational effects, the smaller aerosol particles, mostly transported by means of hydrodynamic drag, form clouds that can remain afloat for long times. In subsaturated air environments, the dependence of pathogen-laden particle dispersion on their size is complicated due to evaporation of the aqueous fraction. Particle dynamics can significantly change when ambient conditions favor rapid evaporation rates that result in a transition from buoyancyto-drag dominated dispersion regimes. To investigate the effect of particle size and evaporation on pathogen-laden cloud evolution, a direct numerical simulation of a mild cough was coupled with an evaporative Lagrangian particle advection model. The results suggest that while the dispersion of cough particles in the tails of the size distribution are unlikely to be disrupted by evaporative effects, preferential aerosol diameters (30-40 μm) may exhibit significant increases in the residence time and horizontal range under typical ambient conditions. Using estimations of the viral concentration in the spewed fluid and the number of ejected particles in a typical respiratory event, we obtained a map of viral load per volume of air at the end of the cough and the number of virus copies per inhalation in the emitter vicinity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Experimental and Computational Analysis of NEAMS Pebble Bed Reactors

The Advanced High-Temperature Pebble Bed Reactor concept leverages a particle-based fuel format consisting of discrete spherical graphite pebbles arrayed in a packed bed architecture. Thermal regulation achieved via flow of gas (e.g., helium) or liquid (e.g., molten salt) coolants through the void spaces between pebbles in the bed (characteristic pebble diameters: ~ 6 cm, gas cooled; ~ 3 cm, liquid cooled). This design is of interest owing to inherent advantages including passive safety features and highly efficient heat transfer characteristics that enable high power densities to be achieved. Fully leveraging these advantages requires a fundamental understanding of the complex coolant flow structure within the pebble-shaped fuel bed. Unfortunately, complexities associated with non-invasively probing the micro-scale flow phenomena within the complex network of void spaces between pebbles have hindered efforts to characterize the underlying transport phenomena. Without this fundamental foundation, it becomes challenging to rationally construct macroscopic transport models based upon volume averages of micro-scale parameters that accurately capture the observed flow behavior. Consequently, this lack of fundamental knowledge has limited previous work to coarse comparisons between theory and extremely limited experimental data. Moreover, these comparisons have been made only in terms of integral or macroscopic results. Therefore, multipoint measurements with a high level of spatial and temporal resolution are critically needed to fully map the complex flow patterns and to provide data at high spatial density to permit accurate volume averaging in the pebble bed. For this project, we performed a coordinated experimental and computational effort to quantitatively map the full-field velocity and temperature fields in the interstitial spaces within the pebble bed. The major accomplishments of this project included the important knowledge acquired for the flow mixing and heat transfer in the interstitial pore scales, the effects of flow mixing in the near-wall and far-wall regions, and the by-pass flows of pebble bed systems. The produced high-resolution experimental data (not broadly available in existing literature) were then used for the numerical studies enabling the refinement and validation of standard simulation tools (Nek5000).

42 ENGINEERING↗
DOE OSTI2020

Second Annual Progress Report on Correlation Between Microstructure and Mechanical Properties of Neutron-Irradiated Ferritic-Martensitic and Austenitic Steels

Ferritic-martensitic steels G92-2b (an optimized Grade 92 heat), NF616 and T91, and austenitic stainless steel 800H and its Grain Boundary Engineering (GBE)-treated version 800H-TMP (ThermoMechanical Processing) were irradiated in the High Flux Isotope Reactor (HFIR) of Oak Ridge National Laboratory (ORNL) and the Advanced Test Reactor (ATR) of Idaho National Laboratory (INL). Selected G92-2b samples were irradiated up to 14.66 dpa in the HFIR at two temperature ranges: 400–496.7°C and 683.3– 720°C. NF616 and T91 were irradiated in the ATR up to 8.16 dpa with the irradiation temperatures ranged from 241°C to 447.5°C. Alloy 800H and 800H-TMP samples were irradiated in both the HFIR and the ATR. Selected 800H and 800H-TMP samples had HFIR irradiation to 1.28 dpa at 580°C and ATR irradiation up to 9.12 dpa at 359°C to 431°C. Vickers hardness measurements, fractography, and microstructural characterization were performed on the selected samples in the Low Activation Materials Design and Analysis (LAMDA) laboratory. Radiation-hardening of G92-2b was observed at the lower doses and lower irradiation temperatures (400- 496.7°C), with GB03 (0.52 dpa at 400°C) and GB04 (7.44 dpa at ~490°C) showing ~12% and ~8% hardening, respectively. Softening by ~14% was observed for GB05 (14.66 dpa at 496.7°C). Radiationsoftening of G92-2b was more prevalent at the higher irradiation temperatures (683.3-~720°C), with GB10 (0.46 dpa at 683.3°C), GB11 (7.44 dpa at ~720°C), and GB12 (14.63 dpa at ~720°C) showing ~8%, ~8%, and ~40% softening, respectively. Radiation-hardening of NF616 and T91 was observed with the hardness increased by ~37% to ~65% depending on the irradiation doses and irradiation temperatures. Within the studied irradiation conditions of NF616 and T91, samples with a higher dose had a larger hardness after irradiation. All the tested alloy 800H and 800H-TMP samples in this work showed radiation-hardening by ~96±7% to ~152±10%. Alloy 800H-TMP tended to have slightly smaller radiation-hardening than alloy 800H. The fractography results of G92-2b sample GB03, GB10, and GB11, together with the previously characterized fractography of GB04, GB05, and GB12, indicated that the ductility of G92-2b was maintained up to 14.66 dpa at the lower irradiation temperatures of 400-496.7°C, while some loss of ductility (less necking) was observed for higher doses at the higher irradiation temperatures of 683.3-720°C. This agrees with the previously reported tensile test results of G92-2b, where the elongation of G92-2b was reduced at higher doses at the higher irradiation temperatures. Dimple sizes increased at higher doses, which are more evident at the higher irradiation temperatures of 683.3-~720°C. The fractography of NF616 sample D2 (2.96 dpa at 291.5°C), D4 (5.91 dpa at 359°C), and D6 (8.16 dpa at 431°C) indicated loss of ductility with negligible necking for sample D2, while ductile failure for samples D4 and D6. Fractography of alloy 800H and 800H-TMP samples in various irradiation conditions showed ductile failure with obvious necking. Dimples were observed, with some of them containing large Ti-rich particles, in all the 800H/800H-TMP samples. Electron backscatter diffraction characterization of GB12 indicated the recovery of the lath structure, which was generally replaced by an equiaxed grain structure. Transmission electron microscopy (TEM) characterization showed the presence of frequent M 23 C 6 (M = primarily Cr), MX (M = primarily V), spherical Nb(C,N) precipitates, and occasional Laves phase precipitates in the G92-2b samples. MX precipitates with sizes of 20-30 nm were observed at boundaries of smaller grains, indicating the pinning effect of the V-rich precipitates. The lath structure recovery was more evident at the higher irradiation temperatures (683.3-~720°C), with decreased densities of line dislocations and M23C6 precipitates. The irradiated T91 (TA04) showed the growth of M23C6 precipitates to 101 ± 40 nm from the initial 68 ± 22 nm in the unirradiated condition. Dislocation loops of both {100} and {111} types were present in TA04. TEM characterization was also performed on the irradiated 800H (N4, N5, N6, and AR2) and 800H-TMP (P4, P5, P6, and HG1). Accumulation of M 23 C 6 precipitates at grain boundaries was observed in all the xii 800H/800H-TMP samples, and the presence of Ti(C,N) precipitates at grain boundaries and in the matrix was observed in the irradiated 800H-TMP. Some Ti(C,N) precipitates are embedded in the M 23 C 6 precipitates, maintaining specific orientation relationships between the precipitates and between the precipitate and matrix. In addition, nanoscale Si-rich clusters were observed in the matrix of all the 800H/800H-TMP samples, with EDS Si maps tending to have a lower contrast in 800H-TMP samples. Atom probe tomography was conducted on the same samples, supported by a Rapid Turnaround Examination project under Nuclear Science User Facilities. The results are being analyzed to be integrated with the TEM results for a confident description of the γ’ precipitates. Dislocation loops also formed in all the 800H/800H-TMP samples. The density and the average size of dislocation loops were quantified to be in the order of 10 22 – 10 23 m -3 and 11.7 – 15.9 nm, respectively, in the ATR-irradiated 800H/800H-TMP samples. The loop densities in irradiated 800H were higher than that in irradiated 800H-TMP under the same irradiation conditions. Further systematic data analyses, together with some complementary experiments, will be pursued for these samples to foster peer-reviewed journal article publications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2025

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nanoscale trace-element zoning in pyrite framboids and implications for paleoproxy applications

Pyrite framboids (spherical masses of nanoscale pyrite) are among the earliest textures of pyrite to form in sediments. It has been proposed that their trace-element (TE) contents can be used to track the TE composition of the water column in which they formed. However, it is not clear how these TEs are associated with the framboidal pyrite grains. For instance, it is important to know whether they are incorporated uniformly or are enriched in different regions of the framboid. Here, we used high-resolution scanning transmission electron microscopy to identify chemical zoning within pyrite framboids. We found that initial, nanoscale pyrite euhedral crystals, which make up the volumetric majority of the framboids, are covered/infilled by later pyrite that templates on the earlier pyrite. Further, this later pyrite is enriched in TEs, suggesting that many TEs are incorporated in pyrite relatively late (during early diagenesis; not in the water column). This observation suggests that although chemical analyses of pyrite framboids may provide ocean-water chemistry trends through time, the details are complex. Specifically, the TEs found in pyrite may be linked to adsorption onto organic matter, detrital material, and authigenic minerals such as Fe- and Mn-oxide phases followed by desorption in the sediments or release via dissolution and incorporation into pyrite as overgrowths on the initial nanoscale euhedral crystals that make up framboids. While the use of pyrite chemistry to understand past ocean conditions remains promising, and even diagenetic additions may not preclude the utility of pyrite for reconstructing ancient ocean conditions, care must be taken in interpretations because the end concentration may be influenced by diagenesis.

58 GEOSCIENCES↗