Engineering Papers⌕ Search

Engineering topics

Hu, Rui (ORCID:0000000237712920)

Publications and source records attributed to Hu, Rui (ORCID:0000000237712920).

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗

Development of SAM Code Capabilities for Safety Analysis of GCR Air-ingress Events

Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.

Yang, Gang↗

Influence of Surfactants on Noble Metal Liquid-Gas Interface Mass Transfer in Molten Salt Reactors

Noble metals, a group of insoluble fission products (FPs), can accumulate in molten salt reactors (MSRs) and influence system performance and safety through deposition. Upon generation in fission, noble metals are transported by circulating salt and may either deposit on structural surfaces or enter the circulating gas bubble interfaces. While sparging was primarily designed for xenon and krypton removal, the mass transfer of noble metals into bubble interfaces highlights a potential pathway for their transport and partial extraction. The properties of the liquid-gas interface therefore play an important role in noble metal deposition and removal behaviors. Observations from earlier studies suggest that noble metals accumulated at the bubble interface can act as surfactants, making the bubble interface partially immobile. This behavior reduces the mass transfer capacity of circulating bubbles and affects both the overall distribution of noble metals in the loop and their removal efficiency. To investigate the influence of surfactants on noble metal transport, deposition, and extraction behaviors in MSRs, a previously developed and validated species transport model incorporating noble metal interphase mass transfer was used to simulate these behaviors. Two designed cases were studied to predict the distribution of 132Te along the Molten Salt Reactor Experiment (MSRE) loop under different bubble interface conditions. The results show that the mass transfer coefficients at the liquid-gas interface, determined by bubble interface characteristics, significantly influence noble metal distribution within the reactor loop. These findings emphasize the importance of continuously introducing fresh circulating bubbles into the molten salt to improve the removal efficiency of insoluble FPs from the MSR primary loop.

MSR↗

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗

Effect of Inlet Throttling on Thermohydraulic Instability in a Large Scale Waterbased RCCS: A System-Level Analysis with RELAP5-3D

This paper presents results from system -level modeling of a water -based reactor cavity cooling system using RELAP5-3D. The computational model is benchmarked with experimental data from a half -scale RCCS test facility at Argonne National Laboratory. The model prediction is first compared with a two-phase oscillatory baseline experimental case where mixed accuracy is obtained. The model shows reasonable prediction of mass flow rate, pressure, and temperature but significant overprediction of void fraction. The model prediction is then compared with a fault case where the inlet of the risers is gradually reduced using a throttling valve. As the valve is closed, the model is able to predict some major flow phenomena observed in the experiment such as the dampening of oscillations, the reintroduction of oscillations, as well as boiling, flashing, and geysering in the risers. However, the timeline of these events are not well captured by the model. The model is also used to investigate the evolution of flow regime in the chimney. This work highlights that the semiempirical constitutive relations used in RELAP-3D could have a strong influence on the accuracy of the model in two-phase oscillatory flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

FY23 Report on Water NSTF Testing at Two-Phase Conditions: Off-normal Scenarios

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report details the experimental activities and upgrades performed during the program’s fifth year of water-based operation. The theme of this year of testing was examining off-normal scenarios. In practice, that meant examining scenarios with blockages or other design basis conditions. From a maintenance and capability standpoint that meant assessing valving options for implementing blockages, assessing drain and refill capabilities, and ensuring in all cases that operators could safely, efficiently, and repeatably bring the facility to an off-normal condition and return it to the normal operating state. In one notable case, no solution existed, nor was a solution commercially available, that met the facility needs for implementing blockages within the two-phase region of the chimney. A valve and feedthrough were custom-designed for NSTF to operate submerged in saturated water while being safely actuated from outside the tank. The fifth year of testing included eight matrix tests consisting of 162 hours of active heating, 9,735 kWh of electrical heating, with seven tests classified as Accepted per NQA-1 and one as Trending. All matrix tests were performed at two-phase flow conditions. The power parametric series from previous years was extended to include another scenario with a decay load equivalent to 1.75 MWt, full scale, further resolving the system flow oscillation response as a function of input power. A depletion test was performed to extend previous work and examine the system’s approach to stagnation and geysering that were not previously observed. This test began at an initial inventory level of 50%, prototypic power of 2.1 MWt, and employed an accelerated drain of 0.75 gpm plus boiloff until stagnation conditions were achieved and multiple geysering events were recorded. Further exploring the geysering phenomenon was a separate effects test examining static boiling in the risers and chimney. This static boiling test was performed at 48 kWe, a higher power than testing in the previous year, and included enhancements to the drain system to remove inventory in a better-controlled and more reliable manner. Two tests were performed, at prototypic conditions of 1.4 MWt and 2.4 MWt, respectively, examining the system response to throttled conditions in the singlephase region of the piping. A ball valve immediately before he riser inlet header was increasingly throttled, progressing the system through a repeatable pattern of two-phase oscillations and stable regimes before ultimately inducing stagnation and geysering. Concluding the year, two tests were performed at baseline conditions of 70% initial inventory level and 2.1 MWt of prototypic decay heat to examine both repeatability and the response to throttling within the two-phase region, at the inventory tank inlet using the custom valve mentioned above. The first such test repeated transient operating conditions and revealed significant sensitivity t

22 GENERAL STUDIES OF NUCLEAR REACTORS↗