SEARCH · Engineering Papers
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.
Infographic on the Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Initiative, an Energy & Environmental Research Center–U.S. Department of Energy project, featuring flowcharts depicting potential pathways from unmined lignite to final products.
Uranium is a very important resource when using nuclear power. Only a small fraction of the uranium we use in the US is domestically sourced. Our goal is to effectively find and mine uranium in a way that is generally accurate and not difficult. Using computer science and machine learning, we want to automate a reasoning system that geologists use to analyze where the uranium ore bodies are. Presenting a general explanation of the goals and impacts of this project for the High School Intern showcase.
This project explores the use of Artificial Intelligence (AI) and Machine Learning (ML) techniques to automate well log analysis for uranium mining. Geophysical log data—spontaneous potential, resistivity, and gamma ray—were used to classify lithology, correlate well logs and identify roll front zonation patterns, which are critical for locating uranium ore bodies. Supervised ML algorithms such as eXtreme Gradient Boosting (XGBoost), Categorical Boosting (CatBoost), and Random Forest were trained to classify lithology with high accuracy. Gradient Boosting Machines (GBM), XGBoost, Random Forest, and Neural Networks were also used for role front zone identification. Moreover, a Fast Dynamic Time Warping (FastDTW) algorithm was employed for well log correlation. Additionally, sample lag was addressed using dynamic programming. Results demonstrate the potential of AI and ML to streamline well log analysis and enhance uranium exploration workflows.
Understanding how point defects in the bulk govern redox transformations is essential for advancing hydrogen-based metal production and designing high-performance oxide materials. This study reveals the atomic-scale mechanisms driving hydrogen-induced reduction of γ-Fe 2 O 3 to Fe 3 O 4 , focusing on how bulk vacancy dynamics dictate structural evolution and reaction kinetics. A key finding is the pronounced contrast in defect behavior between the two oxides: in γ-Fe 2 O 3 , intrinsic Fe vacancies promote oxygen vacancy clustering, destabilizing the local lattice and driving nanopore formation. In contrast, Fe 3 O 4 exhibits a higher oxygen vacancy formation energy and lacks intrinsic Fe vacancies, suppressing vacancy aggregation and maintaining a dense, pore-free structure. This divergence governs distinct reduction pathways—γ-Fe 2 O 3 undergoes an interface-reaction-limited transformation confined to the γ-Fe 2 O 3 /Fe 3 O 4 boundary, while Fe 3 O 4 supports a uniform increase in oxygen vacancy concentration, enabling bulk-phase reduction to lower-oxide FeO. Integrated in situ electron microscopy and density functional theory modeling uncover a vacancy-mediated mechanism, where synergistic cation-anion vacancy dynamics steer microstructure evolution and phase progression. These insights highlight the critical role of vacancy dynamics in controlling oxide reactivity and offer a pathway toward vacancy engineering to enhance reduction kinetics in hydrogen metallurgy and to tailor porosity, reactivity, and structural resilience in oxide-based catalysts and energy materials.
Deposits enriched in rare earth elements (REEs) are abundant and diverse in mineralogy, but those of current economic value have extremely limited geographic distribution. Due to rapidly increasing demand and concerns about supply chain security, new and improved methods of prospecting, beneficiation, separation, purification, and recycling are needed. By overviewing the types and geographic distribution of REE reserves globally, their petrogenesis and major mineralogic assemblages, this review aims to provide a foundation for posing new directions for research and development in these areas. In conclusion, our review also seeks to provide a concise summary of important background information that is accessible to researchers from a variety of disciplines who are seeking to engage in this rapidly expanding field.
Ammonium diuranate (ADU) is commonly encountered in the nuclear fuel cycle; however, previous investigations have shown that ADU is a complex mixture of distinct compounds. Moreover, production parameters are known to heavily influence the composition of the resulting ADU. Here, we examine four samples of ADU prepared at Oak Ridge National Laboratory (ORNL), and one sample of ADU made at Pacific Northwest National Laboratory (PNNL), with the goal of further characterizing and elucidating the effect of processing parameters such as stir rate, strike direction, and temperature on material composition. Process parameters during ADU precipitation at ORNL and PNNL were well documented, and we relate process variables to optical vibrational spectroscopic signatures observed using Raman and infrared (IR) spectroscopy. In addition, powder X-ray diffraction (PXRD) reveals differences in the solid-phase composition of ADU precipitates, but we find that the primary phase is similar to the uranyl oxyhydroxyhydrate mineral metaschoepite. Despite the significant phase contributions of a metaschoepite-like phase, spectroscopic evidence of both nitrate and ammonium are observed for all samples. To gain a more holistic understanding of spectroscopic features of process parameters in ADU, principal component analysis (PCA) is employed and results in observable signatures that relate to the stir rate used during synthesis. These results provide further information about the process-dependence of ADU precipitate composition.
Environmentally unsustainable and toxic chemical flocculants and dispersants from nickel (Ni) processing contribute to industrial effluents that greatly impact biodiversity and aquatic life. Despite the industry’s efforts to reduce its ecological footprint—primarily due to the lack of commercially available biodegradable, environmentally benign, and non-toxic reagents—these reagents continue to harm natural ecosystems. Chemicals used during the processing operations often target a specific unit operation that negatively affects downstream operations. In this study, we discover the concentration-dependent behavior of cellulose nanocrystals (CNCs) as a dispersant and flocculant, avoiding the use of various harmful chemicals in the nickel processing stages. Electrophoretic, optical brightfield microscopy, and quartz crystal microbalance-dissipation studies detected charge neutralization behavior by renewable and biodegradable CNCs, which can greatly benefit Ni processing operations. Microflotation experiments demonstrated that CNCs enhanced Ni recovery from 62 to 77 wt.% and concentrate grade from 15 to 20 wt.%. Settling and turbidity studies demonstrated the dual flocculant-dispersant behavior of CNCs resulting from the alignment of CNC fibers along the octahedral brucite basal and amphoteric edge plane of serpentine. CNCs can be used as dispersants during froth flotation to improve Ni beneficiation and as post-processing flocculants for tailings management and dewatering, which is one of the major environmental and social concerns facing the mining and mineral processing industry. This study paves the way for CNCs to be used as flocculants and dispersants in a range of industries from healthcare to pharmaceuticals to semiconductor devices.
Here, we present the results of the ‘Fast Calorimeter Simulation Challenge 2022’—the CaloChallenge. We study state-of-the-art generative models on four calorimeter shower datasets of increasing dimensionality, ranging from a few hundred voxels to a few tens of thousand voxels. The 31 individual submissions span a wide range of current popular generative architectures, including variational autoencoders (VAEs), generative adversarial networks (GANs), normalizing flows, diffusion models, and models based on conditional flow matching. We compare all submissions in terms of quality of generated calorimeter showers, as well as shower generation time and model size. To assess the quality we use a broad range of different metrics including differences in one-dimensional histograms of observables, KPD/FPD scores, AUCs of binary classifiers, and the log-posterior of a multiclass classifier. The results of the CaloChallenge provide the most complete and comprehensive survey of cutting-edge approaches to calorimeter fast simulation to date. In addition, our work provides a uniquely detailed perspective on the important problem of how to evaluate generative models. As such, the results presented here should be applicable for other domains that use generative AI and require fast and faithful generation of samples in a large phase space.
We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.
The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.
We present iHOMER, an iterative version of the HOMER method to extract Lund fragmentation functions from experimental data. Through iterations, we address the information gap between latent and observable phase spaces and systematically remove bias. To quantify uncertainties on the inferred weights, we use a combination of Bayesian neural networks and uncertainty-aware regression. We find that the combination of iterations and uncertainty quantification produces well-calibrated weights that accurately reproduce the data distribution. A parametric closure test shows that the iteratively learned fragmentation function is compatible with the true fragmentation function.
Direct reduction of hematite pellets with hydrogen (H 2 ) was used to produce directly reduced iron (DRI) in a pilot scale reactor at a pellet feed rate of 21.4 kg/h. At a steady state, operational parameters of the pilot plant (gas recycling rate and inlet temperature) along with key reactor output parameters, the pellet metallization, and the internal temperature profile of the reactor were reported for two scenarios with high recycle and low recycle rate of H 2 . Scenario 1, with a high recycle rate of 400 L/min H 2 along with external heating of 870 °C, gave an average metallization of 91.8%, while Scenario 2, with low recycle rate of 100 L/min H 2 and external heating of 850 °C gave a metallization of 67.8% due to the higher moles of H 2 available for reduction and the external energy required for the endothermic reduction reaction in Scenario 1 as compared with Scenario 2. Finite element analysis was used to build a model of the shaft reactor, which was validated against the metallization and internal temperature profile data. The average metallization values predicted by the model were very close to the metallization values obtained from the pilot plant samples, with 90.9% average metallization for Scenario 1 and 65.6% average metallization for Scenario 2. The internal temperature profiles in the lower region of the reactor obtained from the model were very close to these pilot plant data, with a maximum difference of 52.7 °C and 67.6 °C for Scenarios 1 and 2, respectively. The pilot plant reactor model was used extensively in the commissioning of the pilot plant and to predict the startup outcomes for a given set of operating parameters.
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.
Abstract The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,L jet ∼ 10 50 erg s −1 . We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,L jet ∼ 10 52 erg s −1 , but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2.
Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius, R B = 10 3 R g . The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outside R B and suppress BH accretion to 1.5% of the Bondi rate, M ̇ B . Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of M ̇ B reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% M ̇ B . Because the BAD and RAD states tangle up the jets and destroy them well inside R B , they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.
We expand the theoretical framework by O. Gottlieb et al., which connects binary merger populations with long and short binary gamma-ray bursts (lbGRBs and sbGRBs, respectively), incorporating kilonovae (KNe) as a key diagnostic tool. We show that lbGRBs, powered by massive accretion disks around black holes (BHs), should be accompanied by bright, red KNe. In contrast, sbGRBs—if also powered by BHs—would produce fainter, red KNe, potentially biasing against their detection. However, magnetized hypermassive neutron star (HMNS) remnants that precede BH formation can produce jets with power (P NS ≈ 10 51 erg s –1 ) and Lorentz factor (Γ > 10) likely compatible with sbGRB observations, and would result in distinctly bluer KNe, offering a pathway to identifying the sbGRB central engine. Recent modeling by J. C. Rastinejad et al. found luminous red KNe consistently accompany lbGRBs, supporting their origin in BH-massive disk systems, likely following a short-lived HMNS phase. The preferential association of sbGRBs with comparably luminous KNe argues against the BH engine hypothesis for sbGRBs, while the bluer hue of these KNe provides additional support for an HMNS-driven mechanism. Within this framework, BH–NS mergers likely contribute exclusively to the lbGRB population with red KNe. Our findings suggest that GW170817 may, in fact, have been an lbGRB to on-axis observers. Finally, we discuss major challenges faced by alternative lbGRB progenitor models, such as white dwarf–NS or white dwarf–BH mergers and accretion-induced collapse forming magnetars, which fail to align with observed GRB timescales, energies, and KN properties.
Long-duration γ-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth via black hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over six orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important, and launches bipolar jets. The jets reach the photosphere at ~10 12 cm with an opening angle θ j ~ 6° and a Lorentz factor Γ j ≲ 30, unbinding ≳90% of the star. We find that (i) the disk–jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle θ t ~ 12°, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection θ j + θ t suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs are rarer than currently thought and emerge in only ~0.1% of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where ~10% of the gas maintains magnetization σ ≳ 0.1. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.
Short γ-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is ≲10 -4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳1015 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta's isotropic equivalent mass along the pole must be small (≲10 -2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ~1.5).