Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grain refinement”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
NASA NTRS2007

NavP: Structured and Multithreaded Distributed Parallel Programming

We present Navigational Programming (NavP) -- a distributed parallel programming methodology based on the principles of migrating computations and multithreading. The four major steps of NavP are: (1) Distribute the data using the data communication pattern in a given algorithm; (2) Insert navigational commands for the computation to migrate and follow large-sized distributed data; (3) Cut the sequential migrating thread and construct a mobile pipeline; and (4) Loop back for refinement. NavP is significantly different from the current prevailing Message Passing (MP) approach. The advantages of NavP include: (1) NavP is structured distributed programming and it does not change the code structure of an original algorithm. This is in sharp contrast to MP as MP implementations in general do not resemble the original sequential code; (2) NavP implementations are always competitive with the best MPI implementations in terms of performance. Approaches such as DSM or HPF have failed to deliver satisfying performance as of today in contrast, even if they are relatively easy to use compared to MP; (3) NavP provides incremental parallelization, which is beyond the reach of MP; and (4) NavP is a unifying approach that allows us to exploit both fine- (multithreading on shared memory) and coarse- (pipelined tasks on distributed memory) grained parallelism. This is in contrast to the currently popular hybrid use of MP+OpenMP, which is known to be complex to use. We present experimental results that demonstrate the effectiveness of NavP.

navigational programming (NavP)↗
DOE OSTI2022

Evolution of carbide precipitates in Haynes® 282 superalloy processed by wire arc additive manufacturing

Precipitation hardening regulated with heat treatments is a critical technique for tuning the mechanical properties of Ni-based superalloys fabricated with additive manufacturing, where mechanical deformation is no longer viable for microstructural refinement. As one of the key precipitates, carbides prevail in the additively manufactured alloys, whose precipitation kinetics remain largely vague. This work studies the influence of heat treatments on the evolution of MC- and M 23 C 6 -type carbides within a Haynes® 282 superalloy prepared by wire arc additive manufacturing (WAAM). SEM-EDS analysis was performed to identify the evolution pathways of the carbides, and TEM and XRD characterizations were carried out to track the crystallographic features. The MC-type carbides are originally present in the as-printed alloy, which distribute uniformly across the whole material and remain largely unchanged in size during the heat treatments. By contrast, the incipiently non-existent M 23 C 6 -type carbides precipitate during an aging treatment at 1010 °C, which majorly extend along the grain boundaries as kinetically feasible sites. Further experimental and computational analysis confirms the enhancing effect of the M 23 C 6 -type carbides on the Vickers hardness. For the first time, this work revealed the evolution pathways of carbides within a Haynes® 282 superalloy prepared by WAAM, providing critical information for tunning the microstructure and mechanical properties. Finally, the carbide evolution kinetics established in this work can be potentially extended to other alloys prepared with WAAM, providing critical information to prevent mechanical failures associated with carbide coarsening.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Origin of SNC kaersutitic amphibole: Experimental data

The SNC meteorites, a group of cumulus textured, fine grained diabases, pyroxenites, and dunites, appear to have crystallized at relatively shallow depths on the same SNC parent body. Hydrous minerals generally are not present among the cumulus and intercumulus minerals in these meteorites except for some iddingsite alteration of olivine. The presence of hydrous magmatic amphibole in the SNC melt inclusions indicates that crystallization of the melt inclusions had to take place at significant pressure, probably greater than 1 kb based on previous amphibole stability data. If experimental data for kaersutite amphibole were to be obtained, it should be possible to estimate this pressure more precisely then previously, and to estimate the volatile (H2O) content of the parent magma. At this point, the factors controlling the chemistry and stability of high TiO2, kaersutitic amphiboles are not known. In an attempt to determine the factors which control the stability and chemistry of TiO2-rich amphibole, data was refined and extrapolated from four experimental studies of amphibole-melt equilibria recently completed. At the same time, hypothermal experiments were performed on a composition considered to be an early melt in the Shergotty magma liquid line of descent. The latter experiments were an attempt to reproduce crystallization of the amphibole-bearing melt inclusions.

Rutherford, M. J.↗
NASA NTRS2012

In Situ Carbon and Sulfur Isotope Analysis of Archean Organic Matter and Pyrite

Stable isotopic compositions of biologically important elements (e.g., C and S) in sedimentary rocks are valuable biosignatures to the extent that they indicate the presence and variable expression of microbial metabolisms in space and time. Strong interactions between the carbon and sulfur cycles (e.g., via organic matter remineralization during microbial sulfate reduction) make coordinated, in situ C and S isotope analysis by secondary ion mass spectrometry (SIMS) a particularly powerful tool. In rocks ranging in age from 2.7-2.5 Ga, expansions in the ranges of delta C-13 of organic matter and delta S-34 of pyrite likely reflect the increasing influence of oxygenic photosynthesis in the surface ocean (as well as methane and sulfur metabolisms in deeper waters), whereas the large range of mass independent sulfur isotope fractionation (Delta S-33) suggests that the atmosphere remained anoxic until approx 2.4 Gyr ago. We report in situ delta C-13 measurements of organic matter in the approx 2.7-2.6 Ga Carawine Dolomite, Marra Mamba Iron Formation, and Jeerinah, Wittenoom, and Tumbiana Formations, as well as the approx 2.5 Ga Mount McRae Shale. We also report in situ delta S-24 and Delta S-33 measurements of pyrite associated with organic matter in a subset of these samples. In a single square cm sample of the Tumbiana Formation with bulk delta C-13(sub org) of -49.7% (VPDB), two distinct kerogen types have delta C-13 values, measured in situ, consistent with oxygenic photosynthesis (-33%) and methane metabolism (-52%). In a sample from the ABDP-9 core, radiobitumen associated with a uraniferous mineral grain is C-13-enriched by 8% (-26.8%0) relative to average in situ kerogen (-34.9%0) and similar in delta C-13 to solvent extractable hydrocarbons from the Mount McRae Shale (avg delta C-13 = -27.1 %). Average reproducibility for delta C-13 was 0.4% (2 SD) using a 6 micron spot and 0.8% using a 3 micron spot. In situ sulfur isotope analyses of 33 authigenic pyrite grains in 3 samples of the ABDP-9 core using a 10 micron spot (2 SD reproducibility = 0.4% for delta S-34 and -0.1% for Delta S-33) show a range of 28.1 % in delta S-34 (-10.3 to 17.8%) and 13.3%0 in Delta S-33 (-3.8 to 9.5%), whereas the range from 132 bulk analyses across 84 m of core is 19.4% for delta S-34 and 11.5% for Delta S-33. Coordinated, in situ carbon and sulfur isotope analyses in one ABDP-9 sample are shown. In situ values from this kerogen-pyrite association are within 0.1% of the bulk value in the case of delta C-13 and higher by several permil in the case of delta S-34 and Delta S-33. Coordinated in situ carbon and sulfur isotope analyses in rocks deposited during key intervals of environmental change (e.g., the Great Oxidation Event) can refine our understanding of the mode and tempo of change. In Earth's oldest sedimentary rocks, and in extraterrestrial samples, these coordinated in situ analyses may reveal biosignatures in the form of isotopic correlations at the scale of individual microorganisms and their microhabitats.

Williford, K. H.↗
DOE OSTI2022

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗
NASA NTRS2015

NASAs Evolvable Mars Campaign: Mars Moons Robotic Precursor

Human exploration missions to the moons of Mars are being considered within NASA's Evolvable Mars Campaign (EMC) as an intermediate step for eventual human exploration and pioneering of the surface of Mars. A range of mission architectures is being evaluated in which human crews would explore one or both moons for as little as 14 days or for as long as 500 days with a variety of orbital and surface habitation and mobility options being considered. Relatively little is known about the orbital, surface, or subsurface characteristics of either moon. This makes them interesting but challenging destinations for human exploration missions during which crewmembers must be able to effectively conduct scientific exploration without being exposed to undue risks due to radiation, dust, micrometeoroids, or other hazards. A robotic precursor mission to one or both moons will be required to provide data necessary for the design and operation of subsequent human systems and for the identification and prioritization of scientific exploration objectives. This paper identifies and discusses considerations for the design of such a precursor mission based on current human mission architectures. Objectives of a Mars' moon precursor in support of human missions are expected to include: 1) identifying hazards on the surface and the orbital environment at up to 50-km distant retrograde orbits; 2) collecting data on physical characteristics for planning of detailed human proximity and surface operations; 3) performing remote sensing and in situ science investigations to refine and focus future human scientific activities; and 4) prospecting for in situ resource utilization. These precursor objectives can be met through a combination or remote sensing (orbital) and in-situ (surface) measurements. Analysis of spacecraft downlink signals using radio science techniques would measure the moon's mass, mass distribution, and gravity field, which will be necessary to enable trajectory planning. Laser altimetry would precisely measure the moon's shape and improve the accuracy of radio science measurements. A telescopic imaging camera would map the moon at submeter resolution and photograph selected areas of interest at subcentimeter resolution and a visible and near-infrared (0.4-3.0 mm) imaging spectrograph would produce a global map of mineral composition variations at a resolution of tens of meters and maps of selected areas of interest at meter resolution. Additional remote sensing capabilities could include a thermal infrared imager (heat flow, thermal inertia, and grain size distributions), a gamma-ray and neutron detector (atomic composition), a ground-penetrating radar (internal structure), and a magnetometer and Langmuir probe (magnetic properties and plasma field). Once on the surface of Phobos or Deimos, necessary instrumentation would include a penetrometer (regolith compressive strength), a motion-imagery camera (to observe the penetrometer tests before, during, and after contact), a dust-adhesion witness plate and camera (dust levitation), a microimager (dust particle sizes and shapes), and an alpha-proton-X-ray, X-ray fluorescence, Mossbauer, or Raman spectrometer (atomic and mineral composition of surface materials) and an optional temperature probe (regolith thermal properties). A variety of robotic mission design options to enable both orbital and surface measurements are being considered that include fully integrated and modular approaches. In-situ measurements from at least one surface location would be required, with additional measurement locations possible through use of multiple landers, through propulsive relocation of a single lander, or through electromechanical surface translation by a walking or hopping lander vehicle, which could also serve to evaluate such mobility capabilities for subsequent human missions. Preliminary orbital analysis suggests that remote sensing would likely be performed while in a distant retrograde orbit around the target moon. Mission design options to enable characterization of both Mars’ moons in a single mission are also being studied.

Gernhardt, Michael L.↗
NASA NTRS2022

Local Scale (3-M) Soil Moisture Mapping Using SMAP and Planet Superdove

A capability for mapping meter-level resolution soil moisture with frequent temporal sampling over large regions is essential for quantifying local-scale environmental heterogeneity and eco-hydrologic behavior. However, available surface soil moisture (SSM) products generally involve much coarser grain sizes ranging from 30 m to several 10s of kilometers. Hence a new method is proposed to estimate 3-m resolution SSM using a combination of multi-sensor fusion, machine- learning (ML) and Cumulative Distribution Function (CDF) matching approaches. This method established favorable SSM correspondence between 3-m pixels and overlying 9-km grid cells from overlapping Planet SuperDove (PSD) observations and NASA Soil Moisture Active-Passive (SMAP) mission products. The resulting 3-m SSM predictions showed improved accuracy by reducing ab- solute bias and RMSE by ~0.01 cm3/cm3 over the original SMAP data in relation to in-situ soil moisture measurements for the Australian Yanco region, while preserving the high sampling frequency (1-3 day global revisit) and sensitivity to surface wetness (R 0.865) from SMAP. Heterogeneous soil moisture distributions varying with vegetation biomass gradients and irrigation regimes were generally captured within a selected study area. Further algorithm refinement and implementation for regional applications will allow for improvement in water resources management, precision agriculture, and disaster forecasts and responses.

soil moisture↗
DOE OSTI2022

Unified many-worlds browsing of arbitrary physics-based animations

Manually tuning physics-based animation parameters to explore a simulation outcome space or achieve desired motion outcomes can be notoriously tedious. This problem has motivated many sophisticated and specialized optimization-based methods for fine-grained (keyframe) control, each of which are typically limited to specific animation phenomena, usually complicated, and, unfortunately, not widely used. In this paper, we propose Unified Many-Worlds Browsing (UMWB), a practical method for sample-level control and exploration of physics-based animations. Our approach supports browsing of large simulation ensembles of arbitrary animation phenomena by using a unified volumetric WORLDPACK representation based on spatiotemporally compressed voxel data associated with geometric occupancy and other low-fidelity animation state. Beyond memory reduction, the WORLDPACK representation also enables unified query support for interactive browsing: it provides fast evaluation of approximate spatiotemporal queries, such as occupancy tests that find ensemble samples ("worlds") where material is either IN or NOT IN a user-specified spacetime region. WORLDPACKS also support real-time hardware-accelerated voxel rendering by exploiting the spatially hierarchical and temporal RLE raster data structure. Our UMWB implementation supports interactive browsing (and offline refinement) of ensembles containing thousands of simulation samples, and fast spatiotemporal queries and ranking. We show UMWB results using a wide variety of physics-based animation phenomena---not just JELL-O ® .

Computer Science↗
DOE OSTI2019

Microstructure and Properties of Inconel 718 Fabricated by Directed Energy Deposition with In-Situ Ultrasonic Impact Peening

Many inherent issues, such as the detrimental residual stress, columnar grains with anisotropy, and weak mechanical properties, have severely impeded the adoption of metal additive manufacturing (AM) techniques including powder bed fusion and directed energy deposition (DED) processes. In this study, a hybrid AM process that consists of layer-wise laser metal deposition (i.e., a DED process) and in-situ ultrasonic impact peening (UIP) was applied to obtain Inconel 718 superalloy workpieces. Also, for further property enhancement, a post-heat treatment was applied to the deposited material obtained by the hybrid AM process. Scanning electron microscopy and transmission electron microscope were used to investigate the microstructure morphology and reveal the underlying strengthening mechanism. Electron backscatter diffraction was employed to quantitatively study the microstructure resulted from the hybrid AM process and the post-heat treatment. The profile of residual stress along the depth direction was obtained through X-ray diffraction. The results demonstrate that this hybrid AM process is capable of producing high-quality metal parts with significantly refined microstructure, and beneficial compressive residual stress along the depth into surface. Severe plastic strains are introduced by UIP, and the resulted mechanical twinning and dynamic recrystallization play an important role in refining microstructure. The material microstructure is further refined down to 100 µm, and the texture anisotropy is significantly diminished after solution treatment at 980 °C for 1 hour. Under the as-built condition, in-situ ultrasonic peening alters the residual stress component from a tensile state to an overall compressive state with a maximum value of − 190 MPa within the range of measurement depth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructural evolution in Cu-Nb processed via friction consolidation

Immiscible alloys, whether in well-mixed or layered forms, are of increasing interest based on their novel structural and functional properties, such as enhanced thermal stability against grain growth or radiation-induced defect trapping at the interfaces. To address the need for new approaches to tailor microstructures, the microstructural development of an immiscible Cu-Nb alloy processed via friction consolidation of elemental powders is investigated. Friction consolidation is a solid-phase processing technique that imparts severe plastic strain into a deforming volume resulting in elevated temperatures below the melting temperature of the alloy. Two distinct processing pathways were chosen to understand the effect of thermomechanical conditions on the final microstructure. The microstructure was characterized using scanning electron microscopy, scanning transmission electron microscopy, and X-ray diffraction techniques. Path 1 exhibited larger strain, strain rate, and temperature as compared with path 2. In path 1, agglomerated Nb particles were present in the recrystallized ultrafine-grained Cu matrix, while in path 2 extremely fine and dispersed Nb particles were present in a highly deformed Cu matrix. In both pathways, supersaturation of Cu in Nb lattices was noted, but not vice-versa. The asymmetry in mixing is explained based on deformation-based, thermodynamic and kinetic factors. These findings provide a pathway for creation of novel tailored microstructures and improved properties in any number of binary immiscible alloy systems.

Cu-Nb alloys, Microstructural refinement, Friction↗
DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fleetite, Cu 2 RhIrSb 2 , a New Species of Platinum-Group Mineral from the Miass Placer Zone, Southern Urals, Russia

Fleetite, Cu 2 RhIrSb 2 , a new species of platinum-group mineral (PGM), was discovered intergrown with an Os–Ir–Ru alloy in the Miass Placer Zone (Au–PGE), southern Urals, Russia. A single grain 50 μm across was found. Osmium, ruthenium, and iridium are the main associated minerals; also present are Pt–Fe alloys, laurite, Sb-rich irarsite, Rh-rich tolovkite, kashinite, anduoite, ferronickelplatinum, heazlewoodite, PGE-bearing pentlandite and digenite, as well as micrometric inclusions of forsterite (Fo 93.7 ), chromite–magnesiochromite, and Mg-rich edenite. In reflected light, fleetite is light gray; it is opaque, isotropic, non-pleochroic, and non-bireflectant. We report reflectance values measured in air. A mean of seven point-analyses (wavelength-dispersive spectrometry) gave Cu 13.93, Ni 8.60, Fe 0.10, Ir 28.07, Rh 7.91, Ru 1.96, Sb 39.28, total 99.85 wt.%, corresponding to (Cu 1.41 Ni 0.58 Fe 0.01 ) Σ2.00 (Rh 0.49 Ni 0.36 Ru 0.12 ) Σ0.97 Ir 0.95 Sb 2.08 on the basis of six atoms per formula unit, taking into account the structural results. The calculated density is 10.83 g/cm 3 . Single-crystal X-ray studies show that fleetite is cubic, space group Fdm (#227), a = 11.6682(8) Å, V = 1588.59(19) Å 3 , and Z = 16. A least-squares refinement of X-ray powder-diffraction data gave a = 11.6575(5) Å and V = 1584.22(19) Å 3 . The strongest five reflections in the powder pattern [d in Å(I)(hkl)] are: 6.70(75)(111), 4.13(100)(220), 3.52(30)(311), 2.380(50)(422), 2.064(40)(440). Results of synchrotron micro-Laue diffraction experiments are consistent [a = 11.66(2) Å]. The crystal structure of fleetite was solved and refined to R = 0.0340 based upon 153 reflections with F o > 4σ(F o ). It is isotypic with Pd 11 Bi 2 Se 2 and best described as intermetallic, with all metal atoms in 12-fold coordination. Fleetite and other late exotic phases were formed by reaction of the associated alloy phases with a fluid phase enriched in Sb, As, and S in circulation in the cooling ophiolite source-rock. The mineral is named after Michael E. Fleet (1938–2017) in recognition of his significant contributions to the Earth Sciences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Mineralogical, Elemental and Tomographic Reconnaissance Investigation for CLPS (“METRIC”)

METRIC comprises a suite of two instruments in a Lunar lander that perform X-ray diffraction (XRD) for mineral structure, X-ray fluorescence (XRF) for elemental composition, and X-ray Micro Computed Tomography (XCT) for 3D internal micromorphology. The instruments are accompanied by optical-near IR cameras to provide local geologic context. The Honeybee Robotics PlanetVac pneumatic sampling and transfer system [1],positioned on a lander footpad, will deliver sieved regolith to the X-ray instruments for analysis. The instrument suite is intended for delivery to the lunar surface on a Commercial Lunar Payload Services (CLPS) lander. The METRICXRD/F instrument draws on heritage from the Mars Science Laboratory CheMin instrument [2] and improves upon the design in multiple ways [3]. Like CheMin, Rietveld refinement and full-pattern fitting of METRIC XRD data can identify minerals at a detection limit of ~1 wt.%, quantify their abundances when present at >3 wt.%, and determine mineral composition (e.g., Fo#in olivine) from lattice parameters for minerals present at >5 wt.%[4,5]. An optimized XRF geometry provides for improved detection and quantification of major, minor and trace elements. The METRIC XCT instrument is a miniaturized X-ray Computed Tomography scanner[6]. Lunar regolith is delivered to a 3 mm diameter, 10 mm long graphite tube inside the instrument. The tube is rotated through 360°in 0.9deg.incrementsand a divergent point source X-ray beam is directed through the material. A Charge Coupled Device (CCD)records attenuation images whose brightness and contrast area function of average atomic number and density. Quantitative data, including particle and void sizes, 3D particle shape parameters, modal volumes and pore geometry can be derived from the resulting 3D reconstructions(voxel resolution: 30 μm). Crystal morphologies derived from METRIC XCT data complement the bulk mineralogy determined by the METRIC XRD/F and provide a measure of grain size distribution for the different phases. Taken together, the METRIC instrument suite determines crystal structure, elemental composition and morphology, three principal characteristics of geological materials that are highly useful in determining the origin and subsequent processing of lunar regolith.

Moon↗
DOE OSTI2023

Refined views of ancient ocean chemistry: Tracking trace element incorporation in pyrite framboids using atom probe tomography

The trace element chemistry of pyrite can be used to determine the origin, timing, and conditions of formation of ore deposits; as a vectoring tool for mineral exploration; and to determine the evolution of the Earth’s oceanic and atmospheric chemistry. However, little is known about whether trace elements are held with the pyrite structure or within nano-inclusions of other phases. This distinction is important for two primary reasons. First, trace element incorporation into the pyrite structure can affect the partitioning of other trace elements. Second, if trace elements are held within nano-inclusions, the partition coefficient of the mineral phases that make up the nano-inclusion, rather than pyrite, are the critical consideration in related interpretations. Previous studies addressing this topic have used laser ablation inductively coupled plasma mass spectrometry; however, the resulting data do not provide sufficient spatial resolution to delineate trace element distributions unless the inclusions are large. Further, they use these flat element profiles in time resolved laser ablation output graphs to argue that pyrite trace element content provides a direct relationship between trace element content of seawater and pyrite trace element content. To improve resolution, we have used atom probe tomography to characterize trace element distributions in pyrite framboids from the Cariaco Basin and Demerara Rise in three-dimensions at sub-nanometer resolution. Manganese was found to be concentrated in the pyrite nanocrystal part of the original framboid structure. In contrast, Ni was mostly found along the grain boundary, though it still appeared to be contained within the pyrite structure. Copper was concentrated in later pyrite overgrowths, and As varied in its location. These observations suggest that some important trace elements are incorporated into pyrite during early diagenesis, even in euxinic settings dominated by water-column pyrite formation. Statistical analysis was used to determine whether trace elements were incorporated in the lattice or within nanoscale inclusions (referred to here are nano-inclusions). We found that As, Ni, Cu, and Mn were commonly held within the pyrite structure, but As, Ni, and Cu can also be held as nano-inclusions or within grain interfaces. Incorporation of As is known to enhance the incorporation of other trace elements and in this case appears to correlated to elevated Ni and Cu concentrations in the Cariaco Basin samples. Furthermore, understanding these relationships strongly impacts our ability to utilize pyrite trace metal concentrations to analyze and quantify early ocean chemistry and its evolution through time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

ROADRUNNER uranium nitride MiniFuel: Experimental design, fabrication and pre-irradiation baseline characterization for accelerated burnup testing

Uranium nitride (UN) is a promising fuel candidate for advanced reactor systems owing to its high uranium density and thermal conductivity; however, its qualification remains constrained by the scarcity of well-controlled irradiation performance data. Here, to address this limitation, the ROADRUNNER (Research On ADvancing the peRformance of UraNium Nitrides in Extreme enviRonments) campaign employs the MiniFuel platform in the High Flux Isotope Reactor (HFIR) to enable accelerated burnup irradiation testing under tightly controlled and largely isothermal conditions. This paper presents the experimental design, fuel fabrication, and pre-irradiation baseline characterization of the ROADRUNNER UN MiniFuel campaign. Thirty-six UN minidisc specimens were fabricated with systematically varied as-fabricated density (86–96% of theoretical density), carbon impurity content (961–5240 ppm), oxygen content (≤ ∼2000 ppm), and grain size (2.5–24 μm). The irradiation matrix spans nominal fuel temperatures of 873 K, 1173 K, and 1473 K and target burnups of 3.75%, 6.0%, and 7.5% fissions per initial metal atom (FIMA). Neutronic and thermal analyses were performed to define specimen-specific burnup accumulation and temperature histories, establishing the boundary conditions for subsequent in-pile behavior. Comprehensive pre-irradiation characterization—including dimensional metrology, density verification, impurity analysis, X-ray diffraction, Raman spectroscopy, scanning electron microscopy, X-ray computed tomography, and confocal profilometry—provides a detailed baseline for post-irradiation examination. Pre-irradiation data were further used to generate predictive estimates of fission gas release and swelling using existing empirical correlations. This quantitative comparison reveals substantial inter-model divergence at intermediate and elevated temperatures that exceeds propagated input uncertainties, highlighting structural gaps in the historical irradiation database. The ROADRUNNER irradiation campaign is currently underway in HFIR, with initial firs cycle completed in late 2025 and remaining targets scheduled through 2027. The experimental design and baseline dataset presented here establish the framework needed to interpret forthcoming post-irradiation measurements and to provide discriminating data for the validation and refinement of physics-based UN fuel performance models.

Lopes, Denise Adorno [Oak Ridge National Laborator↗
DOE OSTI2021

Efficient exascale discretizations: High-order finite element methods

Efficient exploitation of exascale architectures requires rethinking of the numerical algorithms used in many large-scale applications. These architectures favor algorithms that expose ultra fine-grain parallelism and maximize the ratio of floating point operations to energy intensive data movement. One of the few viable approaches to achieve high efficiency in the area of PDE discretizations on unstructured grids is to use matrix-free/partially assembled high-order finite element methods, since these methods can increase the accuracy and/or lower the computational time due to reduced data motion. In this paper we provide an overview of the research and development activities in the Center for Efficient Exascale Discretizations (CEED), a co-design center in the Exascale Computing Project that is focused on the development of next-generation discretization software and algorithms to enable a wide range of finite element applications to run efficiently on future hardware. CEED is a research partnership involving more than 30 computational scientists from two US national labs and five universities, including members of the Nek5000, MFEM, MAGMA and PETSc projects. We discuss the CEED co-design activities based on targeted benchmarks, miniapps and discretization libraries and our work on performance optimizations for large-scale GPU architectures. We also provide a broad overview of research and development activities in areas such as unstructured adaptive mesh refinement algorithms, matrix-free linear solvers, high-order data visualization, and list examples of collaborations with several ECP and external applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗