The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non-light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency.This paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

In this study, an integrated coupling method has been developed for solving multiscale fluid-fluid coupling problems in plant-scale safety analysis models in SAM (System Analysis Module) code. In this method, a higher-fidelity multi-dimensional (3D) flow module is used for reactor components of complex flow features (e.g., reactor core) and a lumped parameter one-dimensional (1D) flow module for plant-scale flow loops (e.g., primary loop pipe network), respectively. In this method, the 3D fluid equation/domain and 1D fluid equation/domain are tightly coupled at the residual level and solved simultaneously using the Newton's method to overcome the convergence issues typically seen in existing approaches like separate domain approach, where the 3D fluid equation and 1D fluid equation are solved separately. Extensive and successful code verifications and demonstrations have been performed for this newly developed method. This new modeling approach significantly simplify the work flow in developing high-fidelity plant-scale safety analysis model, e.g. for pool-type reactors and pebble-bed reactors.

One key aspect in analysis of heat pipe microreactors is the efficient modeling of heat pipes and its coupling with the solid reactor core. Various options exist for modeling of heat pipes. Most models require an explicit coupling between the vapor core, which brings in an additional layer of coupling when the heat pipe model is integrated into a system-level safety analysis model. This additional layer of coupling causes both convergence concern and computational burden in practice. This article aims at developing a new heat pipe modeling algorithm, where the heat pipe wall, wick, and vapor core are discretized and coupled in a monolithic fully-implicit manner. The vapor core will be modeled as a one-dimensional compressible flow with the capability of predicting sonic limit inherently; a two-dimensional axisymmetric heat conduction model will be used to model the heat pipe wall and wick region. The heat pipe wick and vapor core are coupled through a conjugate heat transfer interface. Eventually, the coupled system will be solved using the Jacobian-Free Newton-Krylov (JFNK) method. It is demonstrated that the new coupled system works well. Consideration of the vapor compressibility in the two-equation model allows more detailed representation of the vapor core dynamics while remains light-weight in terms of computational complexity. The new model is verified by an approximate analytical solution to the heat pipe vapor core and is validated by a sodium heat pipe experiment.

There has been continuous advancement in numerical techniques and software engineering in the field of system thermal-hydraulics over the past decade. The current state-of-the-art involves the use of fully-implicit high-resolution schemes in conjunction with the Jacobian-Free Newton-Krylov (JFNK) method for both single-phase and two-phase flows. However, the application of these advanced schemes in system thermal-hydraulics codes is not a common practice. Given the opportunity and the need to develop a new modern system thermal-hydraulic code from scratch, there is no doubt that these recent advancements should be leveraged. In this work, we introduce a new modern object-oriented system thermal-hydraulics code named RETA. RETA is characterized by several key features, including: second-order temporal and spatial discretization schemes for both incompressible and compressible flows, as well as single-phase and two-phase flows; fully-implicit solution schemes utilizing both the Newton and PJFNK methods; and an object-oriented design of major fluid components and physical models, which allows for easy extension. A series of verification and demonstration studies are conducted to verify the implementation, demonstrate the accuracy improvement, and evaluate the performance of the new code. It is concluded that results from the second-order scheme are generally more reliable. Without the concern of a code stability issue, the second-order scheme should be the preferred and default choice for a code user or an analyst.

Advanced discretization methods are pursued to improve accuracy of numerical solver for two-phase two-fluid six-equation model. The Weighted Essentially-Non-Oscillatory (WENO) method was a popular high-order discretization one and was successfully applied to many applications. However, development of the WENO-type numerical solver for two-phase two-fluid six-equation model was limited, which was partly due to the lack of analytical eigenvalues and eigenvectors. In previous work, the author derived an approximate analytical eigenvalues and eigenvectors for the two-phase flow model. The analytical eigen-values and eigenvectors were formulated in a compact and structured way and were valid for arbitrary form of equation of state. The analytical eigenvalues and eigenvectors enable the development of a new WENO-type numerical solver for the two-phase flow model. In this work, a new WENO-type numerical solver is developed. Numerical tests show that the newly developed WENO-type solver works very well and is capable of capturing all the characteristic waves and sharp discontinuities. The order of accuracy study shows that the accuracy of WENO is at least third-order for a smooth solution.(c) 2021 Elsevier Ltd. All rights reserved.

In this study, a multiscale model has been developed to include an explicit pebble-temperature model nested in the porous-media model for pebble-bed reactor applications. The multiscale solid-phase energy balance model, including the pebble surface energy balance equation and an explicit modeling of pebble temperature, can predict the macroscopic (pebble bed) and microscopic (pebble) temperature distributions under both steady-state and transient conditions. The proposed multiscale model is solved in a fully coupled manner using the Newton-Krylov method, and therefore iterations between the macroscopic (pebble-bed-scale) and microscopic (pebble-scale) model are avoided. Extensive code verifications, validation, and demonstrations have been performed for this newly developed model. By explicitly modeling pebble temperatures, this new model addresses a major deficiency of the basic porous-media model, which assumes homogeneous solid-phase temperature and is not appropriate for pebble-bed reactor design and safety analyses.

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. The molten-salt-cooled pebble-bed reactor, or pebble-bed FHR (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. To support the development and utilization of the SAM code for PB-FHR safety analysis, activities on SAM code enhancements, reference plant model developments, and code validations have been performed in the past a few years to support near-term industry and NRC needs. This report summarizes recent progress under DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program in SAM code development and demonstration for transient safety analysis of Fluoride-salt-cooled High-temperature Reactors. SAM capabilities has been significantly enhanced over the years to add FHR specific modeling features, including salt freezing and thawing, spherical core channel and pebble bed core modeling, solid-fluid thermal radiation, tritium transport and general species transport in fluids and solids, and the general code enhancements on solver schemes of point kinetics module and reactivity feedback models. A reference PB-FHR model is developed, based on publicly available information from Kairos Power’s generic FHR design and the University of California, Berkeley (UCB) Mk1 design. A reference reactor model is foundational to the methodologies employed by NRC to verify the adequacy of computer codes and evaluation models. The reference FHR model was utilized for a number of selected FHR design basis accidents, including station blackout, loss of heat sink, loss of flow, transient overpower, and overcooling events.

use in safety analyses that support licensing application. Argonne will primarily be responsible for software design, development, and testing; while Kairos will be responsible for software requirements, assessment of software capabilities, needs, gaps, and priorities, and development of proprietary models. The research and development activities for the joint SAM development project include: a series of identification and prioritization studies on design characteristics, event sequences, relevant phenomena, and software capabilities; SAM capability enhancements for specific KP-FHR systems and components; performing code verification and validations; integrating uncertainty quantification (UQ), model calibration, and sensitivity analysis (SA) techniques in safety analyses; and raising the software quality rigor level for commercial-grade applications.

This project aims to raise the SAM code's technical and commercial maturity level to enable the Kairos Power ("Kairos") to use SAM to support its KP-FHR design analysis and licensing application. Argonne National Laboratory ("Argonne") has actively developed and maintained SAM, a modern system-level analysis tool for advanced nonlight water reactors safety analysis. Kairos is currently actively pursuing development of an FHR design and associated technology, and consequently requires a pedigreed safety analysis tool. The purpose of this project is to increase the maturity of the SAM code for the modeling and simulation of the KP-FHR design, thus enabling its use in safety analyses that support licensing application. Argonne will primarily be responsible for software design, development, and testing; while Kairos will be responsible for software requirements, assessment of software capabilities, needs, gaps, and priorities, and development of proprietary models. The research and development activities for the joint SAM development project include: a series of identification and prioritization studies on design characteristics, event sequences, relevant phenomena, and software capabilities; SAM capability enhancements for specific KP-FHR systems and components; performing code verification and validations; integrating uncertainty quantification (UQ), model calibration, and sensitivity analysis (SA) techniques in safety analyses; and raising the software quality rigor level for commercial-grade applications.

The primary objective of this study is to validate the system analysis code, SAM, using experimental data from the Compact Integral Effects Test (CIET) experimental loop. SAM is a modern system analysis code being developed at Argonne National Laboratory for safety analysis of designs for advanced non-light water reactors (non-LWRs), such as sodium-cooled fast reactors, high-temperature gas-cooled reactors, and fluoride salt-cooled high-temperature reactors (FHRs). To support SAM code development for the wide range of non-LWR applications, it is of paramount importance to validate the code against experiments highly relevant to these reactor concepts. The CIET facility, which was designed to reproduce the thermal-hydraulics response of FHRs under both forced- and natural-circulation conditions, has been identified and selected as one of the benchmark test facilities for SAM code validation. In this study, two sets of available CIET tests were selected for SAM code validation purposes, namely, forced convection cooling transient test and steady-state natural-circulation tests. For all selected tests, SAM-predicted results show very good agreement with experimental data. The successful validation of SAM against these selected CIET experiments demonstrates that the computer code is well suited for thermal-hydraulics analysis of FHR designs.

An advanced system analysis tool, SAM, is under development for advanced non-LWR reactor safety analysis, including molten salt reactors (MSR). To support the development and utilization of the SAM code for MSR safety analysis and licensing, continuous efforts have been devoted to enhancing code capabilities and updating reference models for the MSRs. This report documents the FY21 progress in SAM code development, capability enhancements, and reference model development to support transient safety analysis of MSRs, including code enhancements in reactor kinetics and reactivity feedback modeling for liquid fuel reactors, updates of the Molten Salt Reactor Experiment (MSRE) primary system model, developments and updates of the molten salt fast reactor (MSFR) model based on EVOL design, and implementation of a drift flux model for modeling gas transport in MSR systems.

In this study, an explicit model for pebble temperature calculations in the framework of a porous-medium approach has been developed and implemented in the SAM code. The solid-phase energy conservation equation, including the explicit modeling of pebble temperature, is a multi-scale model, which can predict the macroscopic (pebble-bed) and microscopic (pebble) temperature distributions under both steady-state and transient conditions. Extensive and successful code verifications and demonstrations have been performed for this newly developed model. By explicitly modeling pebble temperatures, this new model addresses a major deficiency of the previously implemented model, which assumed a local solid-phase thermal equilibrium that is not appropriate for pebble-bed reactor design and safety analysis applications.

To support the development and utilization of the SAM code for fluoride-salt-cooled high-temperature pebble-bed reactor (PB-FHR) safety analysis and licensing, an effort was devoted to developing reference models for PB-FHR primary loop and the reactor cavity cooling system (RCCS). A reference standard problem of a prototypical reactor design is foundational to NRC and industry to verify the adequacy of computer codes and evaluation models for a specific reactor type. The SAM code was enhanced for the reference PB-FHR model development, including a 1D-3D flow coupling scheme and conjugate heat transfer between porous media and solid structures. The reference FHR primary loop model was developed based on the UC Berkeley Mk1 FHR design with additional design features from the Kairos Power’s KP-FHR core. The PB-FHR core is modeled by a two-dimensional porous medium while the rest of the primary loop is represented by a SAM one-dimensional model. A SAM one-dimensional primary loop model was also developed for comparison using a single channel approach for the reactor core. Furthermore, the water-based NSTF cavity was modeled using SAM multi-dimensional flow module, to assist in evaluating SAM capabilities of modeling the emergency heat removal systems relying on RCCS concepts.

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report summarizes major progress in SAM code development, capability enhancements, demonstration, and validation to support transient safety analysis of advanced non-LWRs. Rapid developments continued in fiscal year 2020 (FY20) to support various needs of the advanced reactor community, especially the NRC and industry on the licensing safety analysis of advanced reactor designs. Significant code changes were made to provide various capability enhancements, bug fixes, and user friendliness improvements. Major code updates are summarized in Section 1, while four important enhancements are detailed in Sections 2-5, including a multi-dimension flow model; reactivity feedback and decay heat models; control and trip system modeling, and additional fluid and solid thermophysical property models. Code validation activities in FY20 include using test data from the Fast Flux Test Facility (FFTF), the High Temperature Test Facility (HTTF), and several separate effects test facilities for pebble-bed modeling.

To support the development and utilization of the SAM code for molten-salt-fueled reactor (MSR) safety analysis and licensing, an effort was devoted to enhancing code capabilities and developing reference models for the MSR primary loop. A reference standard problem of a prototypical reactor design is foundational to NRC to verify the adequacy of computer codes and evaluation models for a specific reactor type. A thermo-fluid model of the Molten Salt Reactor Experiment (MSRE) has been developed. The MSRE primary loop model consists of a 2-D core region and external core components in 1-D or 0-D. Both the steady-state and a transient scenario are simulated. The porous medium model is utilized for both the molten salt fluid channels and the moderator matrix in the MSRE core. The use of the porous medium model for the MSRE core is verified with the higher-fidelity simulation results using 1-D representations of the fluid channels and 3-D modeling of core structures. To further enhance SAM multi-scale simulation capabilities for MSRs, a mass transport model is developed and implemented in SAM multi-dimension flow module for modeling of species transport such as delayed neutron precursors in MSRs.