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.
A multi-modal approach combining in situ X-ray absorption spectroscopy, optical spectroscopy, and ab initio molecular dynamics simulations was used to investigate and understand nickel metal speciation in molten salt systems.
Global markets for energy storage are growing rapidly, with some applications transitioning from traditional LiNi x Mn y Co 1−x−y O 2 (NMC) toward LiFePO 4 (LFP) due to cost, safety, and performance advantages. Battery growth has seeded interest in critical material supplies such as high-purity lithium precursors. In contrast, challenges in securing high-purity iron and phosphorus, historically not considered critical materials, are often overlooked. Precursors must remain inexpensive to maintain LFP's current cost advantage, which leverages the low-cost (∼$\$100$ per t) FeSO 4 byproduct from titanium dioxide manufacturing and will not be available as LFP manufacturing expands. As an alternative, iron is mined primarily for steel manufacturing, for which the existing supply chain and beneficiation process is optimized. LFP batteries require small iron volumes compared to steel (0.11%), but profit margins associated with existing low-cost iron ores and high costs (∼$27 000 per t) associated with low volume high-purity iron oxides may limit interest in manufacturing small volumes of high-purity, specialized iron battery precursors. The complementary LFP precursor, phosphoric acid (H 3 PO 4 ), is primarily utilized in fertilizers. Of the current phosphate ore demand for fertilizers, 4–22.8% would be required to meet projected 2045 LFP H 3 PO 4 demand, suggesting significant supply chain planning is needed to achieve projected demand. If only higher-grade material is considered, demand jumps to 15–77% of current world production. While lithium precursor purity requirements have been evaluated (Li 2 CO 3 is commonly defined as ≥99.5%), “battery grade” iron and phosphorus precursors remain poorly defined, with no internationally accessible and widely adopted standard, further challenging expanding industry by creating manufacturing uncertainty and increasing potential costs.
Current-mode Cherenkov aerogel detectors have been used to measure time resolved spectra of MeV radiographic x-ray sources and have been considered to replace gas Cherenkov systems. The simulations of the Cherenkov response of aerogel detectors have uncertainties in how to model the complex structure of aerogels. Data from three sites—the Idaho electron beam, the Los Alamos Site electron beam and microtron sources at Los Alamos—have been used to calibrate the relative energy response of 15 aerogels, varying in density from 10 to 430 mg/cm 3 , from incident electrons with energies from 1 to 12 MeV. Detector simulations need to introduce density variation in the aerogels to match the collected data.
Liquid-fueled, thermal-spectrum molten salt breeder reactors (TS-MSBRs) offer the potential for affordable, safe, inexhaustible energy with minimal potential for nuclear material misuse and without significant actinide waste generation. Realizing the full set of TS-MSBR capabilities is only now becoming possible with the advent of advanced fuel-salt processing techniques, improved materials, and a more detailed understanding of fuel-salt properties. Additionally, modern higher-fidelity modeling and simulation methods enable a more detailed evaluation of TS-MSBR design options. TS-MSBRs, however, remain immature and will require substantial, sustained development resources.
The ability of nuclear reactors to operate their power conversion cycles more flexibly will enhance their value to energy grids with variable pricing. Current nuclear control systems are typically classical controllers that are often based on proportional-integral-derivative (PID) control. This paper presents a method of augmenting the existing PID control for difficult transient operations in nuclear power plants using a reinforcement learning–derived feedforward signal applied in real time. The agents, which are trained on a test thermal load-following problem, are designed to improve steam generator outlet temperature control for a range of fast load-following scenarios covering ramp rates from 9%/min to 15%/min. Several reinforcement learning algorithms were initially investigated for the training of the feedforward agents with deep Q-learning (DQN) and proximal policy optimization (PPO) networks, which were found to be the most promising. The DQN controllers utilize discrete actions, giving them a better disturbance rejection at steady state but inconsistent response to initial temperature deviations. In contrast, PPO-trained agents, which take continuous actions except for a dead zone around zero, were shown to have the best combination of high disturbance rejection at steady state and good tracking of the desired temperature value. The ability of the PPO agent was also examined, with the average time of decision making found to be on the order of 1 ms. The fault properties of the controller under the loss of the reinforcement learning agent feedforward signal were also examined. The controller showed strong performance in situations of “no-signal” faults. but was less good at handling “stuck-at” faults, where the feedforward signal remains at a set value. In both cases, however, the PID was able to successfully maintain stability, eventually returning the system to a steady state. It is hoped that this work will allow for the proposed control architecture to be examined for more difficult control problems such that it may eventually be used to adapt existing nuclear plants for more aggressive load-following on grids of the future.
Molten salt reactors (MSRs) offer potential operability and safety advantages when compared to commercial light water reactors (LWRs). However, operating experience with MSRs is sparse in comparison to what exists for LWRs. Further, the chemical and isotopic composition of the fuel and/or coolant salt is dynamic and difficult to characterize continuously, posing potential safety, operability, and safeguards unknowns that need to be addressed. A molten salt sampling system (MSSS) is regarded as a necessary subsystem within first generation MSRs used to obtain samples of salt for chemical and isotopic analysis in support of the need to monitor and control salt composition during operation. The MSSS is being developed using the Safety-in-Design (SiD) methodology, which incorporates incremental integration of safety analysis into the design process. The MSSS conceptual design emerging from the application of the early stages of the SiD methodology consists of a sample collection system and its housing, a freeze port, and inert gas control and delivery systems. This article describes the prototypes developed to test the functions of these MSSS subsystems, presents the results of testing in both dry and molten salt environments (including reliability data collection performed in accordance with the principles of SiD and the development of a semiquantitative fault tree model), and summarizes the opportunities for future design and testing enhancements based on the results of prototype testing.
RELAP5-3D is the premier systems thermal-hydraulic code in the world. The code has been used significantly for high-temperature gas-cooled reactor (HTGR) analyses over the last 10 years, but there has not been a comprehensive summary of RELAP5-3D capabilities and methods for HTGR modeling. In this paper, we provide a guide for how novice users can build models for prismatic HTGRs. This guide includes best practices for developing heat structures and setting up radiation and conduction heat transfer. We also provide a summary of relevant validation studies. These studies are divided into three categories: separate effects tests, integral effects tests, and nuclear reactor studies. We present a comprehensive review of RELAP5-3D HTGR applications for reactor modeling over the last 10 years. These applications include work for the Modular High-Temperature Gas-Cooled Reactor benchmark, design activities for a prismatic HTGR as part of the U.S. Department of Energy’s Advanced Demonstration and Test Reactor Options Study, and studies performed to design gas-cooled microreactors. Following the comprehensive literature review, we present a gap analysis highlighting capabilities of RELAP5-3D that are available and missing for a complete analysis of HTGR thermal hydraulics. We also discuss the limitations of existing code capabilities and potential improvements that could improve analyses. Current research efforts have yet to fully validate RELAP5-3D for a prismatic HTGR analysis; however, some of this is attributable to confounding factors in the available experimental data as opposed to code capabilities. In conclusion, this paper presents the most comprehensive review to date on RELAP5-3D capabilities, limitations, and applications for HTGR analysis.
The accurate thermal-hydraulic assessment of thermal storage systems is crucial for addressing the integrity and performance of the storage system design. This is particularly important for thermal storage systems using molten salts as a heat transport and storage medium, as the chemical corrosion and erosion by the high-temperature molten salts can cause serious damage to the structural components of the storage system. To understand these effects for TerraPower’s Natrium® Demonstration Reactor design, the system thermal-hydraulic analysis code, RELAP5-3D, is utilized to analyze the thermal storage system. This study evaluates two heat transfer correlations implemented in the RELAP5-3D code—the Dittus-Boelter and Gnielinski correlations—using selected benchmark cases from two molten salt heat transfer experiments conducted at Xi’an Jiao Tong University and the German Aerospace Center. The Nusselt numbers calculated using the Gnielinski correlation agree with the experimental data within ±5% of the relative deviations at 300°C and 400°C. With the increasing Reynolds numbers, the Dittus-Boelter correlation underestimates the Nusselt numbers by up to 22% compared to the Gnielinski correlation. Comparison of the RELAP5-3D calculation results with experimental data at a high temperature of 550°C revealed that the thermo-chemical characteristics of solar salt at high temperatures surpassing the decomposition temperatures of nitrate salts reduce the accuracy of heat transfer model of the RELAP5-3D code. Specifically, it has been observed that the temperatures of the wall surface and liquid near the wall can surpass the decomposition temperatures, even though the bulk temperature of molten salts remains significantly lower than these decomposition temperatures. In summary, this study recommends the use of the Gnielinski correlation for the heat transfer analysis of molten salt, rather than the Dittus-Boelter correlation.
The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.
Here, analyses are completed for experiments prior to experiment irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Various codes are used to qualify all experiments planned for insertion in the reactor, thereby ensuring that all safety and programmatic requirements are satisfied preirradiation. Among the common experiment analysis tools at INL are MCNP5 coupled to ORIGEN2 (MOPY) and MC21. MOPY uses MCNP5 for transport calculations along with calculations for fluxes and select reaction rates, and then ORIGEN2 handles the step-by-step and postirradiation depletion. MC21 handles all in-reactor transport and step-by-step, during-irradiation, depletion calculations, and then ORIGEN (SCALE 6.2.3) is used for decay and dose calculations postirradiation. The MOPY results, along with those obtained via two variations of the MC21 model, were compared in terms of 238 Pu production in the ATR’s H10 position. For the MOPY model, the MC21 model utilizing the HELIOS-based fission product (FP) library, and the MC21 model utilizing the expanded 1300 FP library, the during-cycle irradiation in-core heating results were sufficiently equivalent; however, the MOPY model and the MC21 model with the HELIOS library showed some differences relating to the respective FP libraries. Ultimately, the MC21 model with a 1300 FP library produced the most consistent results throughout the cycle, whereas the MC21 model that utilized the (smaller) HELIOS library was able to handle during-irradiation analysis but lacked certain short-lived FPs that significantly contributed to the total decay heat at shutdown. MOPY, on the other hand, was found to overpredict fission gas production, as a result of limitations in the ORIGEN2 code.
The pebble tracking transport (PTT) algorithm offers a high-fidelity deterministic approach for neutron transport for pebble bed reactors (PBRs). This approach requires the mesh for the active-core region to consist exclusively of tetrahedral elements, where each node in the pebble-packing region represents a pebble centroid. This paper investigates the application of PTT for full-scale PBRs, considering both the isothermal and the temperature-dependent core conditions. Macroscopic cross sections are generated using Serpent 2 full-core eigenvalue simulations where pebbles are grouped into disjoint subsets using machine learning. To minimize the need for individual cross-section sets for each pebble in the core, K-means clustering is used to group pebbles by temperature and neutronic environment parameters. Here, we compare the multiplication factor and power rate distributions between PTT simulations using the Griffin reactor physics software and reference solutions from Serpent 2. Our analysis shows that a full-core, high-fidelity PTT calculation produces accurate results with minimal local (pebblewise) errors. Additionally, timing results indicate that PTT simulations converge rapidly on modern supercomputing platforms.
Spent tri-structural isotropic (TRISO)–based fuels have a strong track record in storage and transportation without documented incidents. This work seeks to reduce uncertainty to aid in more informed spent fuel management of TRISO-based fuels by modeling both fresh and spent pebble bed reactor (PBR) fuel and comparing the results to the regulatory standards from 10 CFR 71. SCALE was used for all modeling due to it having fast and accurate methods for handling PBR fuel modeling, as well as having an efficient method for shielding calculations in monaco with automated variance reduction using importance calculations (MAVRIC), which utilizes the consistent adjoint-driven importance sampling (CADIS) and the forward-weighted consistent adjoint-driven importance sampling (FW-CADIS) methods. KENO-VI was used for all criticality calculations, TSUNAMI was used for uncertainty quantification on k-effective, TRITON and the Oak Ridge isotope generation code (ORIGEN) were both used for depletion of the fuel, and MAVRIC was used for shielding calculations. For criticality assessments, this study focused on the requirement that the value of the neutron multiplication factor, k-effective (k-eff), would not exceed a peak value of 0.95, including uncertainty, with 95% confidence. Criticality was initially examined by modeling fresh fuel from three different designs—HTR-10 fuel, PBMR-400 fuel, and demonstration fuel representative of a TRISO-fueled modern high-temperature gas reactor (HTGR) design, henceforth referred to as Demo HTGR—and placing them into various sized containers with conditions described in 10 CFR 71 to quantify the peak k-eff state. When the peak value of 0.95 k-eff was exceeded, mitigation methods were examined in those scenarios. Burnup credit, pebble displacement in areas of strong neutron multiplication, and random pebble replacement using pebbles of various compositions and replacement fractions were examined. In summary, the criticality of PBR fuels can be well accounted for by restricting container size, taking credit for burnup, or by displacing/replacing pebbles. Uncertainty of the k-eff due to nuclear data uncertainties was recorded at ~0.6644%Δk/k, or roughly 664% mil (pcm). The nuclear data–induced uncertainty was relatively small and should not require significant modification in the design to be accounted for. Revisions to the evaluated nuclear data file values have been shown to have a larger impact than nuclear data–induced uncertainty. For dose rate aspects, U.S. Nuclear Regulatory Commission regulations require a maximum dose rate of 10 millirem per hour (mrem/h) at 2 meters. In examining the dose rate behavior of spent PBR fuel, the representative Demo HTGR fuel was modeled exclusively due to it possessing the highest target burnup of the examined fuels. Equilibrium cycle modeling methods were used to produce a higher-fidelity discharge isotopic composition than simple assumptions, such as reflected pebbles. The discharge composition was used as a source term in the fixed-source transport shielding calculations, and dose rates were calculated at 2 m for the shortest possible cooling time. The low concentration of fuel material led to dose rates that were in line with regulatory limits, despite the high burnup when compared to traditional light water reactor fuels. In conclusion, the methods employed in this study would require more work to further verify and validate and are limited to the criticality and dose rate analyses performed.
Nuclear power plant (NPP) monitoring and diagnostic centers are actively investigating and implementing automated anomaly detection algorithms to help plants catch anomalies sooner, thereby preventing or reducing the duration of unexpected shutdowns. Current machine learning-based anomaly detection methods are expected to be highly effective during stable, full-power operations because NPPs typically operate as baseload power generators, meaning there are extensive operating data available from plant equipment. However, it is expected that anomaly detection methods will face significant challenges during transient conditions (i.e., when power output falls below full power) because plants only occasionally operate at these lower power levels, generating sparse transient operational data, and resulting in false alarms or missed detections. Here, to address this issue, transfer learning is used, which for this problem leverages knowledge (in the form of learned features) from stable, full-power operations to improve detection accuracy during transient conditions, even with limited data. In this effort, a novel subspace approach is developed to transfer a subset of the data features from full power operation to transients. This approach is validated through experiments using synthetic data and was found to outperform two baseline transfer learning approaches in anomaly detection performance across a range of amounts of transient data used in the training process.
The U.S. nuclear industry is pursuing extensions of light water reactor (LWR) fuel burnup and enrichment limits to approximately 75 GWd/t and 10 wt.% 235 U to achieve economic and operational benefits. A central safety consideration in this effort is the behavior of high burnup (HBu) fuel during loss-of-coolant accidents (LOCAs), particularly fuel fragmentation, relocation, and dispersal (FFRD). Here, this work provides a historical and technical review of U.S. LOCA regulation and experimentation, clarifying how the evolution of Emergency Core Cooling System (ECCS) acceptance criteria in 10 CFR 50.46 has shaped both testing approaches and interpretations of fuel safety. The study revisits the original intent of the ECCS criteria, showing that the peak cladding temperature and equivalent cladding reacted limits were developed as surrogates to preserve a coolable geometry. The explicit inclusion of the coolable geometry criterion in the regulation was intended to emphasize the underlying safety philosophy and as a safeguard against unforeseen failure modes, an intent that remains directly relevant to modern concerns regarding FFRD. The review traces the lineage of HBu LOCA experiments to the Argonne National Laboratory furnace tests, from which subsequent programs at Studsvik, Halden, and Oak Ridge National Laboratory were derived. These tests employed a 5 °C/s heating rate inherited from early embrittlement studies, a stylized temperature history that does not represent actual LWR LOCA thermal-hydraulics. Comparison of these test conditions to pressurized water reactor large break LOCAs and separate effects data indicates that the existing HBu LOCA database may not be fully applicable to all LWR LOCA scenarios, from which a qualitative framework for applicability is proposed.
An irradiation experiment analysis can be informed by high-fidelity reactor engineering depletion results, but this comes at a computational cost. Applying depletion chain simplification to the advanced test reactor driver fuel before performing experiment depletions permits their programmatic parameters to be calculated faster, with a small penalty to accuracy. Here, this work contrasts the results of two irradiation experiments with different neutronic characteristics. Overall, the simplified nuclide library produced using a simple one-group microscopic cross-section library for a pressurized water reactor in the depletion chain simplification process performed comparably in terms of accuracy and runtime to the simplified nuclide library produced using a three-group microscopic cross-section library generated specifically for the advanced test reactor experiments being modeled. This is attributed to the additional nuclides and transmutation pathways preserved in the one-group cross-section library, which has data for 297 nuclides, compared to the three-group cross-section library, which has data for 217 nuclides. This indicates that a cross-section library with more nuclides is better than a cross-section library with fewer nuclides for the depletion chain simplification process, even if the cross-section library with fewer nuclides better represents the flux spectrum of the system being considered.
Tritium, a radionuclide produced through neutron capture by lithium and other elements (beryllium and fluoride) in molten salts, presents unique challenges to radionuclide release. This is true for both fusion energy breeder blankets and molten salt fission reactors. The fundamental understanding of tritium transport is crucial to the safe design and operation of these reactors. The Molten Salt Tritium Transport Experiment (MSTTE), currently under construction at Idaho National Laboratory, aims to investigate tritium transport phenomena using a forced-convection fluoride salt loop. This loop is designed to study various transport mechanisms, such as permeation through metals and gas-liquid interactions, and is intended to support future research on tritium extraction units. A critical aspect of the MSTTE loop design is ensuring a fully developed velocity profile before the fluid reaches the permeation test section where measurements are made. This study employs computational fluid dynamics to model the salt flow behavior within the MSTTE permeation test section. A realizable k-ε turbulent model with enhanced wall treatment is used to simulate the single-phase, vertical upward flow of molten salt FLiNaK under isothermal conditions. The simulation results indicated flow distortion and underdeveloped profiles at all planned flow rates within the test section due to the 85-deg sharp bend. To address this issue, a reduced diameter with a reducer and expander and a flow conditioner are investigated to achieve fully developed flow. The analysis showed that the flow conditioner successfully corrected the flow profile, achieving fully developed behavior at a flow rate of 50 liters per minute (LPM). In conclusion, this research enhances our understanding of flow dynamics in molten salt systems and contributes to optimizing tritium transport control technologies.
Dithiophosphinic acids undergo radiolytic degradation during the extraction of actinides in used nuclear fuel. These will degrade into monothiophosphinic acids and then to phosphinic acids. To elucidate how the complexes that are formed during these radioactive separations change as the ligand degrades, the mixed donor ligand O,O′- diethylmonothiophosphate is chosen as an analog for the monothiophosphinic intermediate. Herein, the monothiophosphate complexes Ln 2 (OPS(OEt) 2 ) 6 (H 2 O) 8 (Ln = La) (La 2 L 6 H 2 O), Ln 2 (OPS(OEt) 2 ) 6 (EtOH) 4 (Ln = La) (La 2 L 6 EtOH), K 2 [Ln(OPS(OEt) 2 ) 5 (H 2 O) 2 ]·H 2 O·CH 2 Cl 2 , (Ln = Ce) (CeL 5 -α), K 2 [Ln(OPS(OEt) 2 ) 5 (H 2 O) 2 ]·H 2 O·CH 2 Cl 2 , (Ln = Pr) (PrL 5 -β), K[Ln(OPS(OEt) 2 ) 4 (H 2 O) 3 ], (Ln = Pr, Sm-Er) (ML 4 ), and K 3 [Ln(OPS(OEt) 2 ) 6 ], (Ln = Dy) (DyL 6 ) were synthesized and characterized using single-crystal X-ray diffraction and optical spectroscopy. Although the lanthanides contract in a nearly linear fashion, the structural changes observed as the f-block is traversed in these compounds are not necessarily a hard line but more so a blend of different structure types possible for each f-element. Furthermore, comparison of the Ln−O bond lengths shows a nearly linear contraction, but the Ln−S bond lengths do not monotonically decrease because of the hard Lewis acidity of the Ln 3+ cations.
Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.