SEARCH · Engineering Papers

Results for “anl”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Energy, greenhouse gas, and water life cycle analysis of synthetic graphite anode production in the United States

This study presents a comprehensive life cycle analysis of potential synthetic graphite battery anode material (BAM) production in the U.S. based on industrial-scale data. The analysis focuses on three impacts: greenhouse gas (GHG) emissions, total energy use, and water consumption. We also conducted sensitivity analyses to evaluate the effect of variation in process parameters and energy sources used for synthetic graphite BAM production on its life cycle GHG emissions. A detailed supply chain analysis of graphite BAM in the U.S. was also undertaken, along with a study of its associated GHG emissions. The results show GHG emissions of 29.7 kg CO 2 -eq. per kg BAM, total energy use of 580 MJ kg −1 BAM, and water consumption of 121 L kg −1 BAM for the baseline condition. The graphitization step is a major process hotspot, contributing to over 74% of all impacts. This is attributed to the energy and material input requirements for this step, particularly through the use of crucibles. Across the entire synthetic graphite production process, electricity is the primary contributor, followed by crucibles used in graphite block production, and then calcined petroleum coke. Sensitivity analyses indicate that improvement in micronization yield, reuse of crucibles, and use of low-carbon nuclear energy can significantly reduce GHG emissions of potential domestic graphite production (by ∼70%). Supply chain analysis identified major graphite BAM sources in the U.S. and showed that the U.S. has a competitive advantage in domestic production of synthetic graphite BAM in terms of reduced life cycle GHG emissions compared to present-day imported sources (by ∼20%).

Battery anode↗

DOE OSTI2025

Cryogenic hybrid magnonic circuits based on spalled YIG thin films

Yttrium iron garnet (YIG) magnonics has garnered significant research interest because of the unique properties of magnons (quasiparticles of collective spin excitation) for signal processing. In particular, hybrid systems based on YIG magnonics show great promise for quantum information science due to their broad frequency tunability and strong compatibility with other platforms. However, their broad applications have been severely constrained by substantial microwave loss in the gadolinium gallium garnet (GGG) substrate at cryogenic temperatures. In this study, we demonstrate that YIG thin films can be spalled from YIG/GGG samples. Our approach is validated by measuring hybrid devices comprising superconducting resonators and spalled YIG films, which exhibit anti-crossing features that indicate strong coupling between magnons and microwave photons. Such new capability of separating YIG thin films from GGG substrates via spalling and the integrated superconductor-YIG devices represent a significant advancement for integrated magnonic devices, paving the way for advanced magnon-based coherent information processing.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Ohmic contacts to nitrogen-incorporated n-type ultrananocrystalline diamond film grown on intrinsic single crystal diamond substrate

Nitrogen-incorporated ultrananocrystalline diamond (n-UNCD) films offer tremendous potential for diverse electronic applications. However, the absence of a reliable Ohmic contact at room temperature limits their practical integration and broader applicability. Here, in this study, we investigate the room temperature specific contact resistivity (ρ c ) of Ti/Pt/Au metal stack deposited on n-UNCD films grown on an intrinsic single crystal diamond substrate using a microwave plasma chemical vapor deposition system. We employ a circular transfer length model (c-TLM) and find the room temperature ρ c to be ∼4.67 × 10 −5 Ω cm 2 , which is among the lowest reported value for n-UNCD films. High temperature vacuum annealing conducted at 700 and 800 °C results in an initial improvement, followed by a minor degradation in ρ c values, respectively. The electrical contacts remain highly Ohmic for all measurements. Furthermore, cross-sectional transmission electron microscopy analysis suggests formation of conductive titanium carbide layer with no significant inter metallic diffusion. Overall, the electrical contacts demonstrate robust thermal stability, both of which are critical for attaining high-performance nanocrystalline diamond-based electronic devices.

Low resistance contacts↗

DOE OSTI2026

Orientation-tunable local crystallization of Si films enabled by atomic imprint crystallization

Here, in this paper, we demonstrate area-selective crystallization of amorphous Si into tunable crystal orientations enabled by atomic imprint crystallization (AIC), where an amorphous Si layer is crystallized by solid phase epitaxy (SPE) from an externally impressed single-crystalline Si template. Using micro-patterned single-crystalline Si templates, a limited area of an amorphous Si film, where the film surface and patterned template surface are in contact, is crystallized via SPE to create an array of crystallographically aligned dots embedded in amorphous matrix. Combining AIC from the top surface and conventional SPE from the substrate, we demonstrate the fabrication of an array of crystalline dots embedded in single-crystalline matrix with tunable in-plane rotation angle. The results indicate the high tunability of the crystallization process enabled by AIC, allowing precise control of crystallographic properties of thin films with area-selectivity; such capability opens opportunities for the design of new materials for a wide range of applications in materials science.

Amorphous↗

DOE OSTI2026

Emergent anisotropic three-phase order in critically doped superconducting diamond films

Two decades since its discovery, superconducting heavily boron-doped diamond (HBDD) still poses fundamental questions that need to be answered to unlock its full potential for quantum applications. We use electrical magnetotransport measurements of critically doped homoepitaxial single crystal HBDD films to reveal signatures of intrinsically granular superconductivity. By studying the dependence of electrical resistivity on temperature and magnetic field vector, we infer that this granularity arises from doping induced disorder. We observe an unexpected three-phase anisotropy in the magnetoresistance, accompanied by a spontaneous transverse voltage (Hall anomaly). Our findings indicate the emergence of an anisotropic order in an otherwise isotropic single crystal HBDD film, offering insights into the mechanism of superconductivity in this quantum material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

DOE OSTI2025

Species Transport Framework Development in SAM for System-Level Tritium Source Term Analysis

The SAM code is under development as a modern system-level modeling and simulation tool for advanced non–light water reactor safety analyses, with recent efforts to add capabilities to evaluate radiological source term risks in these novel reactor concepts. By leveraging the established system-level multiphysics thermal-hydraulic models in SAM, a framework for tightly coupled species transport modeling has been integrated into the code for engineering-scale source term evaluation. This species transport framework was first applied to the simulation of tritium, which is a well-known source term in conventional light water reactors. Tritium poses a unique risk in salt-cooled reactors, especially those with lithium-bearing salts such as the fluoride salt–cooled high-temperature reactor (FHR) concept, as tritium is generated in the salt coolant in significant quantities due to neutron interactions. A compounding factor is the increased mobility of tritium at high temperatures, which is able to permeate through metals while also potentially being retained in graphite pebbles and structures. Engineering-scale models for the tritium transport pathways in a FHR have been developed using the new species transport framework in SAM. The capabilities are assessed through analytical verification problems and validated with data from a graphite retention experiment. In conclusion, the system-level model is demonstrated by performing an initial estimate of baseline tritium generation and flows in a generic reference SAM FHR model, setting a foundation for future studies of source term transient analysis with the potential for further multiscale and multiphysics integration.

SAM↗

DOE OSTI2026

Technoeconomic Design Optimization for Fast Reactors. Part II: Impact of Technoeconomic Constraints on Optimal Design

There is a current drive toward optimizing reactors, particularly small/micro reactors to minimize cost and maximize performance. Previous work has investigated the development of technoeconomic workflows for the design optimization of pool-type fast reactors that aim to deploy into district energy grids. Initial scoping studies verified that the workflow was capable of capturing design trends throughout a variety of design configurations and problem formulations while remaining sufficiently flexible. In this paper, this methodology is applied to understand how cost functions and technoeconomic constraints can drive optimal reactor design. Specifically, the UPu10Zr-fueled fast reactor model from Part I is adapted to include changes in the fissile content limits, control rod worth limits, control rod drive cost, and assumed fuel form. In the case of constraint relaxation at fixed power (fissile content and control rod worth limits), cost sensitivities of 5% to 10% were uncovered. Multi-objective optimization at varying reactor power levels with individualized control rod drives for each assembly (as opposed to one operational and one safety drive) increased cost by $\$10$ to $\$25$ million and substantially altered the optimal core geometry, favoring geometries with substantially fewer control rod placements relative to baseline optimization. Finally, a multi-objective optimization was performed at varying power levels with the fuel form overhauled to metallic, high-assay low-enriched uranium–based U10Zr with more refined fuel cost models. In the case of uranium fueling, the costs increased by at least $50 million relative to the baseline case. Furthermore, economic fuel zoning and lower reactivity swing cores were recovered. Each case serves to demonstrate the value of applying technoeconomic workflows to initial reactor design scoping studies to better understand the trade-off for a proposed concept between different design options.

Argonne Reactor Computation (ARC) codes↗

DOE OSTI2026

Technoeconomic Design Optimization for Fast Reactors. Part I: Workflow Development and Case Study for Small LFR District Energy Application

The nuclear industry is developing small reactor designs that can target a variety of deployment locations and energy products. Smaller nuclear designs have traditionally struggled to handle the steep trade-offs between size and cost that have historically incentivized large reactors. This motivates computational optimization of small reactors to minimize costs and quantify the trade-off between size and cost. In this paper, the cost/size trade-off for a small fast reactor is derived using a multi-objective genetic algorithm optimization, with steady-state, transient, and cost analysis of the fast reactor being performed. Specifically, the method is demonstrated on a small 10- to 120-MW(thermal) U-Pu-Zr–fueled lead-cooled fast reactor with a 10-year core life for district energy applications, which can have a thermal load compatible with this range. The results reinforced that fast reactor cores at the lower end of this power range suffer cost penalties due to critical mass considerations. It was found that high power density cores with strong reactivity swings and many control rods were favored over designing to minimize reactivity swing. Furthermore, this contrasts with some traditional configurations designed using engineering judgment and demonstrates that optimizers can find nontraditional but realistic solutions, along with demonstrating the value of incorporating cost functions into whole-reactor design optimization.

Fast reactor↗

DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗

DOE OSTI2025

High-Fidelity CFD Assessments of Flow Resistance in a 61-Pin Wire-Wrapped Assembly with Partially Blocked Channels

The examination of thermal-hydraulic behaviors in wire-wrapped rod bundles continues to be an active area of research. The sodium fast reactor, a prominent candidate in next-generation nuclear designs, utilizes a hexagonal configuration of wire-wrapped fuel pins. Here, the potential for channel blockage within this compact arrangement poses a significant safety challenge, spurring a number of recent experimental and computational investigations to evaluate its impact on flow and heat transfer. The present work aims to benchmark the high-fidelity NekRS computational fluid dynamics (CFD) solver in predicting the pressure drops associated with substantial blockages, using available experimental data as a reference. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at a range of low to moderate Reynolds numbers (487 ≤ Re ≤ 14 600). The NekRS solver demonstrates an exponential reduction of spatial discretization error with increasing polynomial order. The high level of agreement between the numerical results and measured data confirms the accuracy and consistency of the present numerical approach. This benchmark study establishes the capability of NekRS to perform reliable hydrodynamic simulations for sodium fast reactor applications and supports its use in design, licensing, and safety analyses.

CFD Benchmarking↗

DOE OSTI2024

Systematic Uncertainties from Gribov Copies in Lattice Calculation of Parton Distributions in the Coulomb Gauge

Recently, a new method has been proposed to compute parton distributions using boosted correlators fixed in the Coulomb gauge (CG) within the framework of large-momentum effective theory. This approach, which does not involve Wilson lines, could greatly improve the efficiency and precision of lattice quantum chromodynamics calculations. However, concerns remain regarding whether systematic uncertainties from Gribov copies, which correspond to ambiguities in lattice gauge-fixing, are adequately controlled. This work assesses the effects of Gribov copies on Coulomb-gauge-fixed quark correlators. We utilize different strategies for Coulomb-gauge fixing, selecting two different groups of Gribov copies based on lattice gauge configurations. We examine the differences in the resulting spatial quark correlators in both vacuum and pion states. Our findings indicate that the statistical errors of the matrix elements from both Gribov copies, regardless of the correlation range, decrease proportionally to the square root of the number of gauge configurations. The difference between the strategies does not show statistical significance compared to the gauge noise, demonstrating that the effect of the Gribov copies can be neglected in practical lattice calculations of quark parton distributions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

First on-line commissioning experiments at the St. Benedict facility

Nuclear beta decays provide an excellent probe of fundamental symmetries due to their mediation by the weak interaction. In particular, precise measurements of these decays provide constraints on the unitarity of the Cabbibo-Kobayashi-Maskawa (CKM) quark-mixing matrix. While superallowed pure Fermi decays currently set the most precise limits, the alternative suite of superallowed mixed mirror decays has been ill-studied. These nuclei can provide an important consistency check of calculation and measurement methods employed for the pure Fermi decays, more critically needed now in the wake of a 2.4σ deviation from unitarity of the CKM matrix. In order to remedy the gap in data for mirror decays, the Superallowed Transition Beta-Neutrino Decay Ion Coincidence Trap (St. Benedict) facility is being commissioned at the University of Notre Dame’s Nuclear Science Laboratory (NSL). In this paper, we present first results of the commissioning of the St. Benedict facility on-line at the TwinSol radioactive beam facility. The results of initial commissioning experiments involving the St. Benedict gas catcher, RF carpet, RFQ ion guide and RFQ cooler-buncher will be presented.

beta decay↗

DOE OSTI2025

Emergence and disruption of cooperativity in a denitrifying microbial community

Anthropogenic perturbations to the nitrogen cycle, primarily through use of synthetic fertilizers, is driving an unprecedented increase in the emission of nitrous oxide (N 2 O), a potent greenhouse gas and an ozone depleting substance, causing urgency in identifying the sources and sinks of N 2 O. Microbial denitrification is a primary contributor to biotic production of N 2 O in anoxic regions of soil, marine systems, and wastewater treatment facilities. Here, through comprehensive genome analysis, we show that pathway partitioning is a ubiquitous mechanism of complete denitrification within microbial communities. We have investigated mechanisms and consequences of process partitioning of denitrification through detailed physiological characterization and kinetic modeling of a synthetic community of Rhodanobacter thiooxydans FW510-R12 and Acidovorax sp. GW101-3H11. We have discovered that these two bacterial isolates, from a heavily nitrate (NO 3 − ) contaminated superfund site, complete denitrification through the exchange of nitrite (NO 2 − ) and nitric oxide (NO). The process partitioning of denitrification and other processes, including amino acid metabolism, contribute to increased cooperativity within this denitrifying community. We demonstrate that certain contexts, such as high NO 3 − , cause unbalanced growth of community members, due to differences in their substrate utilization kinetics. The altered growth characteristics of community members drives accumulation of toxic NO 2 − , which disrupts denitrification causing N 2 O off gassing.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Reducing model error using optimized galaxy selection: weak lensing cluster mass estimation

Galaxy clusters are one of the most powerful probes to study extensions of General Relativity and the Standard Cosmological Model. Upcoming surveys like the Vera Rubin Observatory’s Legacy Survey of Space and Time are expected to revolutionise the field, by enabling the analysis of cluster samples of unprecedented size and quality. To reach this era of high-precision cluster cosmology, the mitigation of sources of systematic error is crucial. A particularly important challenge is bias in cluster mass measurements induced by the mismodelling of photometric redshift estimates of source galaxies. This work proposes a method to optimise the source sample selection in cluster weak lensing analyses drawn from wide-field survey lensing catalogs to reduce the bias on reconstructed cluster masses. We use a combinatorial optimisation scheme and methods from variational inference to select galaxies in latent space to produce a probabilistic galaxy source sample catalog for highly accurate cluster mass estimation. We show that our method reduces the critical surface mass density Σ crit relative modelling bias on the 60-70% level, while maintaining up to 90% of galaxies. We highlight that our methodology has applications beyond cluster mass estimation as an approach to jointly combine galaxy selection and model inference under sources of systematics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Lifetimes of two- and four-neutron-hole states in 204 Pb and 206 Pb

The centroid-shift method was used to determine half-lives of excited levels in 204 Pb and 206 Pb. For the states with 𝐼 𝜋 = 12 + , 17 − , and 19 − in 204 Pb, values of 𝑇 1/2 = 4.9⁢(3), 3.9(3), and 0.9(3) ns, respectively, were determined. The half-life of the 𝐼 𝜋 = 9 − level in 206 Pb was measured to be 1.1(3) ns. In 204 Pb, the 12 + state arises from two neutron holes in the 𝑖 13/2 subshell, while the 17 − and 19 − levels are built on 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑝$^{−1}_{1/2}$) and 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑓$^{−1}_{5/2}$) four-neutron-hole configurations, respectively. The 9 − state in 206 Pb results from a 𝜈⁡(𝑖$^{−1}_{13/2}$, 𝑓${−1}_{5/2}$) excitation. As a result, the transition probabilities for the decay of these isomeric levels are reasonably consistent with the expected single-particle values, and these states are described well by shell-model calculations using the kshell code with the KHHE interaction.

190 ≤ A ≤ 219↗

DOE OSTI2025

Core-excited microsecond isomer in 207 Pb

An isomeric state with T 1/2 = 1.46(10) μs and E x = 8835 keV has been established in 207 Pb. Here, this state likely arises from the occupation of neutron holes in the high-j i 13/2 and j 15/2 subshells following the excitation of the 208 Pb core. The high-spin excited level structure of 207 Pb has been extended beyond the previously known states up to E x ≈ 6 MeV. Additionally, a more precise value of the half-life of the I π = 10 + isomer in 208 Pb has been obtained. Shell-model calculations provide a good account of the structure up to the 31/2 + level in 207 Pb, while the energies of higher-spin states built from exclusively high-j nucleons are considerably underestimated.

190 ≤ A ≤ 219↗

DOE OSTI2025

New scaling and nuclear structure aspects in heavy-ion fusion reactions

Three new behaviors have been found in comparisons of fusion cross sections for different collision systems. root (1) Replacing the energy E with a scaling one, E scal = (E-V g )/($\sqrt{2}$W g ), is successful for washing out the Coulomb interaction in the spectra of fusion cross sections, where V g and W g are barrier height and width of the single-Gaussian barrier distribution model. (2) In a representation of σE vs the scaling energy, E scal , all data sets display in parallel. Here, the ratio for sigma E from any two fusion systems over the whole range is a constant value. That behavior is also studied in another representation, in which the data sets display as parallel horizontal lines for any heavy-ion fusion system. (3) The constant ratio value is the ratio of parameter products, $R^2_gW_g$, of the two systems; where R g is the barrier radius obtained in the single-Gaussian barrier distribution model. Moreover, when comparing neighboring collision systems at the same E scal , the ratio of sigma is near a constant value within a few percent over the whole range. Thus a quantitative comparison for the fusion enhancement for neighboring systems is developed. The present finding could be beneficial for predicting unmeasured fusion cross sections.

Jiang, C. L. [Argonne National Laboratory (ANL), A↗