Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models & methods for nuclear reactions”

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 OSTI2024

Predicting Radiation-Induced Plutonium Redox Chemistry using Multi-scale Modeling Methods

Over the the last 70 years plutonium (Pu) has been integral in the development of several technologies that have changed the world, yet our fundamental understanding of its chemistry is still far from complete. This is a testament to this element’s unique and complex properties, such as its ability to coexist as multiple oxidation states in aqueous solution. Careful manipulation of plutonium oxidation states is essential in the study and utilization of its rich chemistry. To achieve this level of control, a comprehensive mechanistic understanding of radiation-induced plutonium redox chemistry is critical due to the unavoidable exposure of plutonium to ionizing radiation fields, both inherent and from in-process applications. For this reason, we have developed an experimentally evaluated multi-scale computer model for the prediction of gamma radiation-induced Pu(IV) redox chemistry in concentrated nitric acid solutions (1.0, 3.0, and 6.0 M). Under these acidic, aqueous solution conditions, cobalt-60 gamma irradiation afforded negligible net change in the steady-state oxidation state distribution of Pu(IV). Multi-scale calculations, which are in excellent agreement with experimental data, indicate that this observation is due to radiation-induced redox cycling between Pu(IV) and Pu(III), as achieved by the reduction of Pu(IV) by radiolytic nitrous acid and hydrogen peroxide, and the oxidation of Pu(III) by nitrate and hydroxyl radicals. These radiation-induced redox processes are augmented by plutonium’s inherent disproportionation reactions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE

Analyses of cross section sensitivity data from systems with fissile material allow analysts to associate an importance for each material, nuclide, reaction, and neutron energy by simulating real world criticality scenarios. Although criticality safety validation efforts can be guided by the cross-section sensitivity and uncertainty data generated for a particular system, these calculations can often be computationally expensive and sometimes cumbersome without proper guidance. The TSUNAMI suite within the SCALE code package has several methods for generating sensitivity data, including multigroup and continuous energy (CE) capabilities. The release of SCALE 6.3 has three different CE methods for generating cross section sensitivity data: (1) the Iterated Fission Probability (IFP) method with the KENO Monte Carlo transport solver, (2) the IFP method with the Shift Monte Carlo transport solver, and (3) the Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method with the KENO Monte Carlo transport solver. Although the CLUTCH method has additional parameters for generating sensitivity data files relative to the IFP method, all three methods use latent generations, which are the generations between an event (i.e., fission) and the assessment of importance based on the asymptotic population of progeny neutrons. Increasing the number of latent generations in a calculation leads to increased discrimination of the sensitivity coefficients but at the cost of the increased uncertainty associated with those generated values. Analysts must balance the accuracy of the sensitivity calculations and its uncertainty with the associated computational cost involved in generating the values. This paper discusses the impact of adjusting the latent generation parameter for a range of sensitivity values and how these changes compare with the direct perturbation values obtained from a change of ±0.5% Δ k in both benchmark and safety application models. Two benchmarks from the International Handbook of Evaluated Criticality Safety Benchmark Experiments and the MPC-32 dual purpose canister for spent nuclear fuel are used for analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Optical follow-up observation for GW event S190510g using Subaru/Hyper Suprime-Cam

Abstract A gravitational wave event, S190510g, which was classified as a binary-neutron-star coalescence at the time of preliminary alert, was detected by LIGO/Virgo collaboration on 2019 May 10. At 1.7 hours after the issue of its preliminary alert, we started a target-of-opportunity imaging observation in the Y band to search for its optical counterpart using the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The observation covers a 118.8 deg2 sky area corresponding to $11.6\%$ confidence in the localization skymap released in the preliminary alert and $1.2\%$ in the updated skymap. We divided the observed area into two fields based on the availability of HSC reference images. For the fields with the HSC reference images, we applied an image subtraction technique; for the fields without the HSC reference images, we sought individual HSC images by matching a catalog of observed objects with the PS1 catalog. The search depth is 22.28 mag in the former method and the limit of search depth is 21.3 mag in the latter method. Subsequently, we performed visual inspection and obtained 83 candidates using the former method and 50 candidates using the latter method. Since we only have the one-day photometric data, we evaluated the probability of candidates being located inside the 3D skymap by estimating their distances with photometry of associated extended objects. We found three candidates are likely located inside the 3D skymap and concluded they could be a counterpart of S190510g, while most of the 133 candidates were likely to be supernovae because the number density of candidates was consistent with the expected number of supernova detections. By comparing our observational depth with a light-curve model of such a kilonova reproducing AT2017gfo, we show that early deep observations with the Subaru/HSC can capture the rising phase of the blue component of a kilonova at the estimated distance of S190510g (∼230 Mpc).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2026

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗
DOE OSTI2022

Checking nonflow assumptions and results via PHENIX published correlations in p+p, p+Au, d+Au, and 3 He+Au at √ s NN=200 GeV

Recently the PHENIX Collaboration has made available two-particle correlation Fourier coefficients for multiple detector combinations in minimum bias p+p and 0–5% central p+Au, d+Au, and 3 He+Au collisions at √ s NN = 200 GeV [Phys. Rev. C 105, 024901 (2022)]. Using these coefficients for three sets of two-particle correlations, azimuthal anisotropy coefficients v 2 and v 3 are extracted for midrapidity charged hadrons as a function of transverse momentum. In this paper, we use the available coefficients to explore various nonflow hypotheses as well as to compare the results with theoretical model calculations. The nonflow methods fail basic closure tests with ampt and pythia/angantyr, particularly when including correlations with particles in the low multiplicity light-projectile going direction. In data, the nonflow adjusted v 2 results are modestly lower in p+Au and the adjusted v 3 results are more significantly higher in p+Au and d+Au. However, the resulting higher values for the ratio v 3 /v 2 in p+Au at RHIC compared to p+Pb at the LHC is additional evidence for a significant overcorrection. Incorporating these additional checks, the conclusion that these flow coefficients are dominated by initial geometry coupled with final-state interactions (e.g., hydrodynamic expansion of quark-gluon plasma) remains true, and explanations based solely on initial-state glasma are ruled out. The detailed balance between intrinsic and fluctuation-driven geometry and the exact role of weakly versus strongly coupled prehydrodynamic evolution remains an open question for triangular flow, requiring further theoretical and experimental investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Neutron-neutron scattering length from the 6 He(p,pα)nn reaction

We propose a novel method to measure the neutron-neutron scattering length using the 6 He(p,pα)nn reaction in inverse kinematics at high energies. The method is based on the final-state interaction (FSI) between the neutrons after the sudden knockout of the α particle. We show that the details of the neutron-neutron relative-energy distribution allow for a precise extraction of the s-wave scattering length. Here, we present the state of the art in regard to the theory of this distribution. The distribution is calculated in two steps. First, we calculate the ground-state wave function of 6 He as a αnn three-body system. For this purpose we use Halo effective field theory, which also provides uncertainty estimates for the results. We compare our results at this stage to model calculations done with the computer code face. In a second step we determine the effects of the nn FSI using the nn t-matrix. We compare these FSI results to approximate FSI approaches based on standard FSI enhancement factors. While the final distribution is sensitive to the nn scattering length, it depends only weakly on the effective range. Throughout we emphasize the impact of theoretical uncertainties on the neutron-neutron relative-energy distribution, and discuss the extent to which those uncertainties limit the extraction of the neutron-neutron scattering length from the reaction 6 He(p,pα)nn.

6 ≤ A ≤ 19↗
DOE OSTI2020

Quasielastic lepton scattering and back-to-back nucleons in the short-time approximation

Understanding quasielastic electron- and neutrino-scattering from nuclei has taken on new urgency with current and planned neutrino oscillation experiments, and with electron scattering experiments measuring specific final states, such as those involving nucleon pairs in “back-to-back” configurations. Accurate many-body methods are available for calculating the response of light (A ≤ 12) nuclei to electromagnetic and weak probes, but they are computationally intensive and only applicable to the inclusive response. In the present work we introduce a novel approach, based on realistic models of nuclear interactions and currents, to evaluate the short-time (high-energy) inclusive and exclusive response of nuclei. The approach accounts reliably for crucial two-nucleon dynamics, including correlations and currents, and provides information on back-to-back nucleons observed in electron and neutrino scattering experiments. We demonstrate that in the quasielastic regime and at moderate momentum transfers both initial- and final-state correlations, and two-nucleon currents are important for a quantitatively successful description of the inclusive response and final state nucleons. Finally, the approach can be extended to include relativistic—kinematical and dynamical—effects, at least approximately in the two-nucleon sector, and to describe the response in the resonance-excitation region.

6 ≤ A ≤ 19↗
DOE OSTI2022

Novel fabrication tools for dynamic compression targets with engineered voids using photolithography methods

Mesoscale imperfections, such as pores and voids, can strongly modify the properties and the mechanical response of materials under extreme conditions. Tracking the material response and microstructure evolution during void collapse is crucial for understanding its performance. In particular, imperfections in the ablator materials, such as voids, can limit the efficiency of the fusion reaction and ultimately hinder ignition. To characterize how voids influence the response of materials during dynamic loading and seed hydrodynamic instabilities, we, in this paper, have developed a tailored fabrication procedure for designer targets with voids at specific locations. Our procedure uses SU-8 as a proxy for the ablator materials and hollow silica microspheres as a proxy for voids and pores. By using photolithography to design the targets’ geometry, we demonstrate precise and highly reproducible placement of a single void within the sample, which is key for a detailed understanding of its behavior under shock compression. This fabrication technique will benefit high-repetition rate experiments at x-ray and laser facilities. Insight from shock compression experiments will provide benchmarks for the next generation of microphysics modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Production of neutron-rich isotopes for 𝑍 ≥ 98 in the 238 U + 248 Cm reaction

Multinucleon transfer (MNT) reactions in actinide systems are a promising method to synthesize transuranium neutron-rich elements. Appropriate theoretical approaches are needed to understand the mechanisms behind MNT. We employ a microscopic approach to calculate neutron-rich isotope production in the reaction 238 U + 248 Cm system. Here, the stochastic mean-field (SMF) approach is used to calculate the primary cross sections in MNT reactions based on the quasifission and inverse quasifission processes, and a statistical de-excitation model with GEMINI ++ code to calculate the secondary fragment cross sections. The calculated cross sections using SMF and GEMINI ++ explain available experimental results for the 238 U + 248 Cm system at 𝐸 c.m. = 898.7 MeV energy. This shows the effectiveness and applicability of the quantal diffusion approach, based on the SMF theory, in heavy-ion collisions. Production of transuranium neutron-rich elements with a proton number up to 𝑍=101 is obtained with sizable cross sections. Theoretical results calculated for the 𝑍=102–105 region, for which there are no experimental data, show that the cross-section values would be lower than the microbarn level. SMF theory does not contain any adjustable parameters other than the standard parameters of the energy density functional used in the TDHF theory and is an important approach for the microscopic understanding of reaction mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

𝐵⁡(𝐸⁢2) measurements in the yrast band of 28 Mg: Implications for the 𝑁 = 20 island of inversion

High-precision lifetime measurements in 28 Mg were performed to study neutron shell evolution in Mg isotopes and the onset of the N = 20 island of inversion. Here, using both the recoil distance and Doppler shift attenuation methods, five lifetimes were measured in addition to six upper limits. The observation of two long-lived, negative-parity states demonstrate the importance of studying Mg isotopes for the contribution of intruder configurations to sd-shell nuclei. Lifetimes of the $2^+_1$ and $4^+_1$ states of 1.81(5) ps and 172⁢($^{+11}_{–10}$) stat.⁢ (4) stop. ⁢(8) feed.⁢ (4) targ. fs, respectively, demonstrate a loss of collectivity with increasing spin in the yrast band, permitting for distinguishing between current theoretical models. These measurements also highlight the progression of yrast structure across the Mg isotopic chain from rotational at N = 12 to large shape mixing at N = 16 and back to collective behavior at N = 20 but with dominating intruder configurations.

20 ≤ A ≤ 38↗
DOE OSTI2020

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

In-situ irradiation-induced studies of grain growth kinetics of nanocrystalline UO 2

The thermal conductivity of UO 2 fuel needs to be high enough to dissipate the heat generated from the fission reaction. Since grain size affects thermal conductivity and grain size can evolve with irradiation, it is critical to understand in-reactor UO 2 grain growth. Most studies of grain growth in UO 2 are based on thermally driven processes at elevated temperatures. However, studies have shown that grain growth can occur even at cryogenic temperatures by ballistic processes. Such irradiation-induced grain growth in UO 2 is yet to be studied. Advanced in-situ Kr ion irradiation and transmission electron microscopy were systematically performed on nanocrystalline UO 2 thin films at temperatures ranging from 50 K to 1073 K; grain growth was observed at all temperatures. A combination of manual and machine learning techniques was used to measure and plot grain size evolution against irradiation fluence at various irradiation temperatures. The machine learning method has significantly improved the analysis efficiency and reduced human labors. The grain diameter data were fitted using classical grain growth and thermal spike models to describe grain growth kinetics with and without irradiation effect. Grain growth during low temperature irradiation (≤ 475 K) can be well described by the thermal spike model. Above 475 K, there were additional thermally assisted processes that further accelerate the grain growth. At the highest irradiation temperature about 1075 K, both irradiation-induced dislocation loops and cavities/bubbles were observed to form in the UO 2 . In this report, the effects of irradiation-induced defects on grain growth kinetics are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2011

Z-Pinch Magneto-Inertial Fusion Propulsion Engine Design Concept

Fusion-based nuclear propulsion has the potential to enable fast interplanetary transportation. Due to the great distances between the planets of our solar system and the harmful radiation environment of interplanetary space, high specific impulse (Isp) propulsion in vehicles with high payload mass fractions must be developed to provide practical and safe vehicles for human spaceflight missions. Magneto-Inertial Fusion (MIF) is an approach which has been shown to potentially lead to a low cost, small fusion reactor/engine assembly (1). The Z-Pinch dense plasma focus method is an MIF concept in which a column of gas is compressed to thermonuclear conditions by an estimated axial current of approximately 100 MA. Recent advancements in experiments and the theoretical understanding of this concept suggest favorable scaling of fusion power output yield as I(sup 4) (2). The magnetic field resulting from the large current compresses the plasma to fusion conditions, and this is repeated over short timescales (10(exp -6) sec). This plasma formation is widely used in the field of Nuclear Weapons Effects (NWE) testing in the defense industry, as well as in fusion energy research. There is a wealth of literature characterizing Z-Pinch physics and existing models (3-5). In order to be useful in engineering analysis, a simplified Z-Pinch fusion thermodynamic model was developed to determine the quantity of plasma, plasma temperature, rate of expansion, energy production, etc. to calculate the parameters that characterize a propulsion system. The amount of nuclear fuel per pulse, mixture ratio of the D-T and nozzle liner propellant, and assumptions about the efficiency of the engine, enabled the sizing of the propulsion system and resulted in an estimate of the thrust and Isp of a Z-Pinch fusion propulsion system for the concept vehicle. MIF requires a magnetic nozzle to contain and direct the nuclear pulses, as well as a robust structure and radiation shielding. The structure, configuration, and materials of the nozzle must meet many severe requirements. The configuration would focus, in a conical manner, the Deuterium-Tritium (D-T) fuel and Lithium-6/7 liner fluid to meet at a specific point that acts as a cathode so the Li-6 can serve as a current return path to complete the circuit. In addition to serving as a current return path, the Li liner also serves as a radiation shield. The advantage to this configuration is the reaction between neutrons and Li-6 results in the production of additional Tritium, thus adding further fuel to the fusion reaction and boosting the energy output. To understand the applicability of Z-Pinch propulsion to interplanetary travel, it is necessary to design a concept vehicle that uses it. The propulsion system significantly impacts the design of the electrical, thermal control, avionics, radiation shielding, and structural subsystems of a vehicle. The design reference mission is the transport of crew and cargo to Mars and back, with the intention that the vehicle be reused for other missions. Several aspects of this vehicle are based on a previous crewed fusion vehicle study called Human Outer Planet Exploration (HOPE), which employed a Magnetized Target Fusion (MTF) propulsion concept. Analysis of this propulsion system concludes that a 40-fold increase of Isp over chemical propulsion is predicted. This along with a greater than 30% predicted payload mass fraction certainly warrants further development of enabling technologies. The vehicle is designed for multiple interplanetary missions and conceivably may be suited for an automated one-way interstellar voyage.

Miernik, Janie H.↗
DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Developing a microwave-driven reactor for ammonia synthesis: insights into the unique challenges of microwave catalysis

Rapid development of new ammonia (NH 3 ) synthesis techniques that enable modular, intermittent production is essential to actualizing NH 3 's potential as a clean energy carrier, since contemporary methods are configured to centralized, continuous production methods with high emissions and are incompatible with renewable sources. Here, in this mission, microwave-driven catalysis is promising for its ability to enhance reaction kinetics and apply targeted heating for efficient energy use. However, owing to an incomplete understanding of the interaction between microwave fields and catalyst beds, the development of such microwave-catalysis systems remains underexplored and challenging. This paper investigates the 10× scale-up of a microwave-based NH 3 synthesis reactor using numerical and experimental approaches, achieving the largest reported microwave-driven NH 3 reactor to date. Results elucidate phenomena unique to microwave processes that are only predictable through numerical modeling, including how a catalyst's dielectric properties influence microwave field distribution by affecting penetration depth and how energy utilization can be poor even with sufficient attenuation. These dynamics change with scale, constrain reactor geometry, and potentially hamper performance. Nonetheless, we demonstrate a production rate of 56.6 g NH$_3$ per day, the highest reported NH 3 synthesis rate for laboratory-scale alternative techniques; correspondingly, the benchmark energy efficiency achieved in this paper (45.6 g NH$_3$ kW h -1 ) is the highest reported for such reactors of sufficient scale. Even with this exemplary energy efficiency, energy losses were found in excess of 50%, an issue resolvable through scale and reactor design. The efficiencies imparted by microwaves were key in these achievements, warranting further investigation toward development of microwave-driven NH 3 systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗