Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heavy elements”

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 451 records · Page 25

Atomic Resolution Imaging of Light Elements in a Crystalline Environment using Dynamic Hollow-Cone Illumination Transmission Electron Microscopy

Multiple electron scattering and the nonintuitive nature of image formation with coherent radiation complicate the interpretation of conventional transmission electron microscopy images. Precession of the illuminating beam in transmission electron microscopy (TEM) can lead to more robust and interpretable images with some penalty to image contrast, a technique known as dynamic hollow-cone illumination TEM. Here, we demonstrate direct and robust imaging of light and heavy atoms in a crystalline environment with this technique. This method is similar to the annular bright-field technique in scanning transmission electron microscopy, via the principle of reciprocity. Dynamic hollow-cone illumination TEM is challenging in practice due to sensitivity to the misalignment of the precession axis, microscope objective aperture, and crystal zone axis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hydrogen Storage System Modeling: Public Access, Maintenance, and Enhancements

The overall vision of this project is to provide ownership and support for maintaining existing material-based hydrogen storage systems models. This includes making models accessible to the research community through a public web page and updating and enhancing storage systems models to support material developers in assessing their materials relative to DOE vehicle-level targets. Key elements for FY24: Continue to update and enhance existing models for broader application and user-friendliness. Develop tools to evaluate the performance of hydrogen storage materials developed under HyMARC activities or other fundamental hydrogen storage materials discovery research. Expand the application of current hydrogen storage models beyond light-duty vehicles to include medium-, heavy-duty, and mining vehicles and stationary application(s). Develop models for alternatives, such as LOHCs, to material-based systems (liquefied & gaseous H2) and compare for various mobile and stationary use cases.

ENERGY STORAGE,HYDROGEN↗

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Organic compounds alter the preference and rates of heavy metal adsorption on ferrihydrite

The availability of heavy metals in terrestrial environments is largely controlled by their interactions with minerals and organic matter, with iron minerals having a particularly strong role in heavy metal fate. Because soil organic matter contains a variety of compounds that differ in their chemical properties, the underlying impact organic matter-soil mineral associations bestow on heavy metal binding is still unresolved. Here, we systematically examine the binding of Cd, Zn and Ni by a suite of organic-ferrihydrite assemblages, chosen to account for various compound chemistries within soil organic matter. We posited that organic compound functionality would dictate the extent of association with the organic-ferrihydrite assemblages. Increased heavy metal binding to the assemblages was observed and attributed to the introduction of additional binding sites by the organic functional groups with differing metal affinities. Furthermore, the relative increase depended on the metal's Lewis acidity and followed the order Cd > Zn > Ni, whereas the reverse order was obtained for metal binding by pristine ferrihydrite (Ni > Zn > Cd). Citric acid-, aspartic acid- and cysteine-ferrihydrite assemblages also enhanced the metal binding rate. X-ray absorption spectroscopy revealed that the organic coating contributed significantly to Zn binding by the assemblages, despite relatively low organic surface coverage. Our findings provide valuable information on the nature of heavy metal-organic-mineral interactions and metal adsorption processes regulating their bioavailability and transport.

54 ENVIRONMENTAL SCIENCES↗

Methods for Estimating Hydrogen Fuel Tank Characteristics

The pressure vessels needed to store hydrogen for next-generation hydrogen fuel cell vehicles are expected to be a substantial portion of the total system mass, volume, and cost. Gravimetric capacity, volumetric capacity, and cost per kilogram of usable hydrogen are key performance metrics that the U.S. Department of Energy (DOE) uses to determine the viability of hydrogen fuel cell systems. Research and development related to hydrogen storage systems covers a wide range of potential operating conditions, from cryogenic temperatures to high temperatures (above ambient) and low pressure to high pressure. Researchers at PNNL have developed methods for estimating these key pressure vessel characteristics to support on-board hydrogen storage system design and performance evaluation and to support decision-making about DOE hydrogen storage system research investments. This article describes the pressure tank estimation methodology that has been used as a stand-alone calculation and has been incorporated into larger system evaluation tools. The methodology estimates the geometry, mass, and material cost of type I, type III, and type IV pressure vessels based on operating pressure and material strength at the system's operating temperature, using classical thin-wall and thick-wall pressure vessel stress calculations. The geometry, mass, and material cost requirements of the pressure vessel have significant impacts on the total system performance. For example, hydrogen storage materials that can separately achieve a very high hydrogen density can be deemed impractical for use in fuel cell vehicle hydrogen storage systems because the pressure tank containing them is too large, heavy, or expensive. This article describes the design philosophy and analytical process of the tank characteristic estimation methodology, which has been implemented in spreadsheet calculation tools and system-level analysis tools used by DOE researchers. Each of the three tank types (type I, type III, and type IV) uses a different analysis methodology with some common elements. This article also provides examples of implementing the methodology to perform parametric studies of all three pressure vessel types. The goal of this article is to present the methodology in sufficient detail so it can be implemented in other hydrogen fuel cell vehicle design and analysis tools.

42 ENGINEERING↗

Data‐Driven Insights into Rare Earth Mineralization: Machine Learning Applications Using Functional Material Synthesis Data

Understanding rare‐earth element (REE) mineralization mechanisms is essential for developing efficient separation strategies. Although the geochemical pathways that generate REE deposits are qualitatively known, quantitative links between specific conditions and mineralization outcomes remain limited. Herein, the repurpose laboratory REE hydrothermal synthesis data—originally collected for functional‐materials fabrication—as a surrogate for studying mineralization with data‐driven methods. The compiled 1,200+ hydrothermal reaction records and trained three machine‐learning models—K‐nearest neighbors (KNN), random forest (RF), and extreme gradient boosting (XGB)—to predict product elements and phases from precursors, additives, reaction conditions, and engineered features. Validation shows XGB achieves the highest accuracy. Feature importance indicates thermodynamic properties of cations and anions dominate model decisions. Correlations reveal positive relationships among precursor concentration, reaction time, pH, and temperature, consistent with classical crystallization behavior. XGB‐based regressors are built to predict crystallization temperature and pH from precursor/product attributes. Performance is strongest when similar training examples exist, while accuracy declines for underrepresented reactions, notably REE carbonates and heavy‐REE systems. Overall, the study shows that functional‐materials datasets can illuminate REE mineralization and provide priors for exploration and processing. Expanding datasets with less‐studied chemistries and conditions will improve generality and support deposit discovery and more efficient REE recovery.

feature importance analysis↗

Low-Temperature Plasma Treatment for Enhanced Recovery of Highly Valued Critical REEs from Coal (Final Technical Report)

The principle objective of this project was to develop a novel process flowsheet using low-temperature plasma treatment (LTP) integrated with conventional hydrometallurgical processes to recover rare earth elements (REEs), especially highly valued REEs (i.e. scandium and critical REEs), from coal and coal byproducts. LTP treatment was used to pretreat the solid feed material prior to leaching in an effort to improve the leaching characteristics of the REEs and potentially reduce acid consumption. Initially, a detailed experimental program was performed on a laboratory LTP furnace to identify a set of process parameter values to be used for assessing the enhanced leaching performance achieved on four specific gravity fractions of West Kentucky No. 13 and Fire Clay coarse coal refuse materials. High-temperature oxidation (HTO)) was also evaluated as a means for improving REE leachability on the same fractions of both coal sources. The lixiviants used to assess leachability improvements included i) untreated water, ii) 0.1 M ammonium sulfate solution and iii) 1.2 M sulfuric acid solution at a temperature of 75°C. The laboratory data was used to develop multiple process flow sheets using HTO only to produce a rare earth oxide mix concentrate or utilizing both LTP and HTO treatment processes to generate a high purity heavy REO mix concentrate and a light REO mix concentrate. A techno-economic analysis was performed on each circuit.

01 COAL, LIGNITE, AND PEAT↗

Extending JuTrack’s capabilities to the FRIB accelerator to enhance online modeling

JuTrack is a Julia-based accelerator modeling and tracking package that utilizes compiler-level automatic differentiation (AD) to enable fast and accurate derivative calculations. While JuTrack provides a solid foundation for beam dynamics simulations, its capabilities must be extended to support the Facility for Rare Isotopes (FRIB) linac. This includes modeling heavy-ion linac accelerator components such as the liquid-lithium charge stripper, which facilitates efficient acceleration by remove electrons from heavy isotopes, and incorporating multi-charge state acceleration tracking, which allows for charge-dependent beam dynamics. These extensions address challenges such as the beam matching and optimization of multi charge state through various accelerating structures and beam-material interaction modeling while maintaining the auto differentiation capability. This work focuses on adapting JuTrack to incorporate these elements, enhancing its online modeling abilities. We present modifications to JuTrack’s framework and demonstrate their performance in FRIB simulations.

Accelerator Physics↗

Process design for recovering rare-earth elements from mine tailings with low rare-earth concentrations via sequential leaching and solvent extraction

Rare earth elements (REEs) are essential for advanced technologies and yet face significant supply chain risks due to their concentrated global production and limited domestic availability. Addressing this challenge requires efficient processes capable of upgrading low-grade secondary resources such as mine tailings. In this study, we developed a novel separation flowsheet that integrates sequential leaching and 2-stage solvent extraction (SX) processes to recover high-purity heavy REEs (HREEs) and light REEs (LREEs) from a simulated mine-tailing concentrate containing 2.4 wt% total REEs (TREEs; 0.6 wt% LREEs and 1.8 wt% HREEs). Sequential leaching with controlled pH adjustment selectively precipitated REEs while retaining the large amount of impurities in the solution, producing an REE-enriched leachate by following leaching processes with roughly twice the REE concentration and half the impurity concentration compared to that of single-step leaching. The optimized SX flowsheet employed Cyanex 572 to extract HREEs and Fe over LREEs, followed by Fe removal using tributyl phosphate (TBP), while the raffinate stream was processed by SX with di(2-ethylhexyl)phosphoric acid (D2EHPA) to recover LREEs under optimized conditions balancing both extraction efficiency and purity. Although increased extractant availability in the organic phase improved LREE recovery, it also increased co-extraction of Ca, underscoring trade-offs in process optimization. Both HREE- and LREE-rich solutions were subsequently precipitated into solid products via oxalate precipitation, resulting in high-purity REE solids containing ∼92.0 wt% HREEs (∼95.7 wt% TREEs) and ∼92.8 wt% LREEs (∼94.0 wt% TREEs). In conclusion, this proof-of-concept study using simulated mine tailings demonstrates a promising approach for upgrading low-grade REE resources, while highlighting the need for future validation with real materials.

Mine tailings↗

Atomic-scale mechanisms of void strengthening in tungsten

Neutron and heavy-ion irradiation of tungsten produces nanometer-size vacancy voids, gas-filled bubbles and dislocation loops. These defect features can affect mechanical properties and the impact can be quite significant because of their high density. Understanding the basic mechanisms of mechanical properties degradation is necessary for predicting radiation effects. Predictions can be made using discrete dislocation dynamics or/and finite element approaches which, however, need local interaction mechanisms as inputs. Such knowledge can be provided only by atomic-scale modeling. This paper reports the results of an extensive atomic-scale modeling study of the interactions between moving edge dislocations and voids in tungsten. Furthermore, the main focus is on the effects of the void size and ambient temperature. Critical resolved shear stress was calculated for voids up to 9 nm in diameter. Atomistic results are compared with the theoretical approach and with those obtained earlier for voids in body centered cubic (bcc) iron. An important role of the void surface has been revealed.

36 MATERIALS SCIENCE↗

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗

Investigation of post-breakup Coulomb acceleration using a trajectory model

Intermediate mass fragments ejected during the deexcitation of excited projectilelike fragments may promptly decay following ejection; the daughter particles that are subsequently produced are subject to interactions with the residual nucleus that affect final-state observables, a process herein referred to as post-breakup Coulomb acceleration. A simple classical Coulomb interaction model was used to study modification of 8 Be (2 + ), 5 Li (3/2 – ), 7 Li (7/2 – ), 7 Be (7/2 – ), and states in 12 B emitted following heavy-ion collisions of 28 Si + 12 C at 35 MeV/nucleon. Here, in contrast to previous work studying 8 Be (2 + ), excellent agreement between simulation and experiment was obtained using only Coulomb forces when either a Lorentzian or R-matrix line shape was used to describe the initial relative energy rather than a Gaussian. In consideration of the obtained results, improvements to the model and evaluation of experimental data are discussed as future directions, but it was concluded that the effects observed in the present data can be accurately described using only elements of classical mechanics and that the process is largely understood for a wide range of state lifetimes and mass and charge (a)symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Pore-resolved investigation of turbulent open channel flow over a randomly packed permeable sediment bed

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, Re K = 2.56 , 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). Here, the mean flow and shear penetration depths increase with Re K and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

turbulence simulation↗

Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. Here, we perform a thorough survey of the electron lepton number (ELN) crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $ν_e$ and $\bar{ν}_e$ angular distributions, or the ELN crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $\textit{r}$ ≈ 16 – 21 km . We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $ν_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy, and nucleosynthesis of heavy nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗

The MASSIVE Survey. XVI. The Stellar Initial Mass Function in the Center of MASSIVE Early-type Galaxies

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universe $\mathrm{log}({M}_{\star }/\,{M}_{\odot })\gt 11.2$. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01 μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture of R e /8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, α IMF = (M/L)/(M/L) MW , ranges from 1.1 to 3.1, with an average of $\langle$α IMF $\rangle$ = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σ uncertainty. We complement our sample with lower-mass galaxies from the literature and confirm that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with $\mathrm{log}(\sigma )$, $\mathrm{log}({M}_{\star })$, and $\mathrm{log}({M}_{\mathrm{dyn}})$. From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density Σ Kroupa might be responsible for the variation of α IMF . We conclude that σ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutrinoless double beta decay rates in the presence of light sterile neutrinos

We investigate neutrinoless double-beta decay (0νββ) in minimal extensions of the Standard Model of particle physics where gauge-singlet right-handed neutrinos give rise to Dirac and Majorana neutrino mass terms. We argue that the standard treatment of these scenarios, based on mass-dependent nuclear matrix elements, is missing important contributions to the 0νββ amplitude. First, new effects arise from the exchange of neutrinos with very small (ultrasoft) momenta, for which we compute the associated nuclear matrix elements for the decays of 76 Ge and 136 Xe. These contributions can dominate the 0νββ rate in cases with light sterile neutrinos. The ultrasoft terms are also relevant in the more standard scenario of just three light Majorana neutrinos where they lead to a 10% reduction of the total 0νββ amplitude. Secondly, we highlight the importance of short-range terms associated with medium-heavy sterile neutrinos and provide explicit formulae that can be used in phenomenological analyses. As examples we discuss impact of these new effects in several explicit scenarios, including a realistic 3 + 2 model with two right-handed gauge-singlet neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗