Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transient Testing Reactor”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Improvement and Verification of Online Cross Section Generation Capability of Griffin for TRISO-fueled Reactors

Griffin, a MOOSE-based reactor multiphysics code jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE Office of Nuclear Energy’s NEAMS program, has pursued the development of an online multigroup cross section generation capability for a few years to enable high-fidelity, problem-dependent neutronics analyses of advanced thermal reactors. Recent advancements in Griffin’s online multigroup cross section generation capability have significantly improved the accuracy, robustness, and efficiency of self-shielding calculations for both prismatic and pebble-bed TRISO-fueled reactor applications. Key developments include a unified fuel self-shielding method applicable to both TRISO and annular compact/spherical shell fuel zone geometries; an advanced Dancoff Category-based Equivalence Theory using a bell function for non-fuel resonance treatment, achieving more than an order-of-magnitude speedup compared to the Tone method; an on-the-fly multigroup equivalence approach to mitigate group condensation errors; and a streaming correction method for pebble-bed homogenization. A proof-of-concept demonstration of on-the-fly group condensation with consistent P0 transport correction was also achieved. The method reproduced direct fine-group solutions with excellent accuracy (eigenvalue errors within 10 pcm and pin-power differences within 0.5%), but due to performance limitations of the current fixed-source solver, improvements to solver efficiency will be addressed in future work. Verification tests were performed on graphite-moderated TRISO-fueled two-dimensional core benchmark problems representing gas-cooled microreactors, heat pipe-cooled microreactors, gas-cooled pebble-bed reactors, and fluoride salt-cooled high-temperature reactors. Across all cases, Griffin showed excellent agreement with Serpent2 continuous energy Monte Carlo solutions: eigenvalue errors within 200 pcm, pin-power root-mean-square errors within 2%, and control rod and drum worth errors less than 2%. It should be noted that, for the benchmark problem, cross section generation contributed less than 3% of the total simulation times. These results demonstrate that Griffin’s online cross section generation capability delivers accurate and efficient reactor physics solutions across a wide spectrum of TRISO-fueled advanced reactor designs. With further improvements to the fine-group fixed-source solver and planned extensions to depletion, transients, and coupled neutron–gamma transport, Griffin will be well-positioned to become a powerful and comprehensive tool for advanced reactor analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The Monolithic Heat Pipe Microreactor Reference Plant Model

This work introduces a reference plant model for a generic monolithic heat-pipe-cooled microreactor. The model will serve as a springboard to develop future evaluation models in the licensing process of similar microreactor designs at the U.S. Nuclear Regulatory Commission. This model has been developed with the Comprehensive Reactor Analysis Bundle and its specifications are based on open literature publications for the eVinci TM design. BlueCRAB is the U.S. Nu- clear Regulatory Commission non-light-water reactor analysis system based on the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications includes tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C. No moderator or burnable poison pins are used in the design. The fuel enrichment is reduced to control excess reactivity in the core. This core design is not optimized and only serves for testing purposes, since the primary objective of this work is to exercise the multiphysics coupling for this type of reactor system. A three dimensional (3D) core heterogeneous Griffin discrete ordinates (SN) transport model allows the precise calculation of the flux distribution and pin powers. Griffin transfers the power density distribution and obtains a temperature distribution to and from BISON. The BISON model com- putes the 3D core temperature distribution and is coupled to 876 Sockeye subapplications running a heat pipe model. This 3D conduction model is coupled to the various heat pipes via heat flux boundary conditions. The model includes a small gap between the heat pipe and the monolith. Convective heat transfer boundaries with either ambient temperature or condenser temperature as heat sinks are imposed at the model boundaries. The 2D Sockeye heat pipe model uses a vapor- only methodology, which provides the needed resolution for transient calculations and allows the determination of various heat pipe limits. This approach is superior to the superconductor model traditionally used in steady-state calculations. BlueCRAB computes steady-state power and temperature distributions that serve as the initial condition for a loss-of-heat-sink transient simulation. The steady-state results show significant peaking due to the position of the control drum, but this is a characteristic of the particular design used, which is not optimized at this stage. The transient results show the reactor power slowly stabilizing towards a 3% power level after the partial loss of secondary heat removal. Several recriticalities are observed due to cooling through the secondary system but the reactor is self-stabilizing and behaves as expected.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Listening to Radiation Damage In Situ: Passive and Active Acoustic Techniques

Knowing when, why, and how materials evolve, degrade, or fail in radiation environments is pivotal to a wide range of fields from semiconductor processing to advanced nuclear reactor design. A variety of methods, including optical and electron microscopy, mechanical testing, and thermal techniques, have been used in the past to successfully monitor the microstructural and property evolution of materials exposed to extreme radiation environments. Acoustic techniques have also been used in the past for this purpose, although most methodologies have not achieved widespread adoption. However, with an increasing desire to understand microstructure and property evolution in situ, acoustic methods provide a promising pathway to uncover information not accessible to more traditional characterization techniques. This work highlights how two different classes of acoustic techniques may be used to monitor material evolution during in situ ion beam irradiation. The passive listening technique of acoustic emission is demonstrated on two model systems, quartz and palladium, and shown to be a useful tool in identifying the onset of damage events such as microcracking. An active acoustic technique in the form of transient grating spectroscopy is used to indirectly monitor the formation of small defect clusters in copper irradiated with self-ions at high temperature through the evolution of surface acoustic wave speeds. Here, these studies together demonstrate the large potential for using acoustic techniques as in situ diagnostics. Such tools could be used to optimize ion beam processing techniques or identify modes and kinetics of materials degradation in extreme radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Integral Reactor Containment Condensation Model and Experimental Validation

This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). In the three years of investigation duration, following the original proposal, the planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). A film flow model was developed for the scaling analysis, and the results suggest that the 1/3 scaled test facility covers large portion of laminar film flow, leading to a lower average heat transfer coefficient comparing to the prototypic value. Although it is conservative in reactor safety analysis, the significant reduction of heat transfer coefficient (50%) could under estimate the prototypic condensation heat transfer rate, resulting in inaccurate prediction of the decay heat removal capability. Further investigation is thus needed to quantify the scaling distortion for safety analysis code validation. Experimental investigations were performed in the existing MASLWR test facility at OrST with minor modifications. A total of 13 containment condensation tests were conducted for pressure ranging from 4 to 21 bar with three different static inventories of non-condensable gas. Condensation and heat transfer rates were evaluated employing several methods, notably from measured temperature gradients in the HTP as well as measured condensate formation rates. A detailed mass and energy accounting was used to assess the various measurement methods and to support simplifying assumptions required for the analysis. Condensation heat fluxes and heat transfer coefficients are calculated and presented as a function of pressure to satisfy the objectives of this investigation. The major conclusions for those tests are summarized below: (1) In the steam blow-down tests, the initial condensation heat transfer process involves the heating-up of the containment heat transfer plate. An inverse heat conduction model was developed to capture the rapid transient transfer characteristics, and the analysis method is applicable to SMR safety analysis. (2) The average condensation heat transfer coefficients for different pressure conditions and non-condensable gas mass fractions were obtained from the integral test facility, through the measurements of the heat conduction rate across the containment heat transfer plate, and from the water condensation rates measurement based on the total energy balance equation. 15 (3) The test results using the measured HTP wall temperatures are considerably lower than popular condensation models would predict mainly due to the side wall conduction effects in the existing MASLWR integral test facility. The data revealed the detailed heat transfer characteristics of the model containment, important to the SMR safety analysis and the validation of associated evaluation model. However this approach, unlike separate effect tests, cannot isolate the condensation heat transfer coefficient over the containment wall, and therefore is not suitable for the assessment of the condensation heat transfer coefficient against system pressure and noncondensable gas mass fraction. (4) The average condensation heat transfer coefficients measured from the water condensation rates through energy balance analysis are appropriate, however, with considerable uncertainties due to the heat loss and temperature distribution on the containment wall. With the consideration of the side wall conduction effects, the results indicate that the measured heat transfer coefficients in the tests is about 20% lower than the prediction of Dehbi’s correlation, mainly due to the side wall conduction effects. The investigation also indicates an increase in the condensation heat transfer coefficient at high containment pressure conditions, but the uncertainties invoked with this method appear to be substantial. (5) Non-condensable gas in the tests has little effects on the condensation heat transfer at high elevation measurement ports. It does affect the bottom measurements near the water level position. The results suggest that the heavier non-condensable gas is accumulated in the lower portion of the containment due to stratification in the narrow containment space. The overall effects of the non-condensable gas on the heat transfer process should thus be negligible for tall containments of narrow condensation spaces in most SMR designs. Therefore, the previous correlations with noncondensable gas effects are not appropriate to those small SMR containments due to the very poor mixing of steam and non-condensable gas. The MELCOR simulation results agree with the experimental data reasonably well. However, it is observed that the MELCOR overpredicts the heat flux for all analyzed tests. The MELCOR predicts that the heat fluxes for CCT’s approximately range from 30 to 45 kW/m2 whereas the experimental data (averaged) ranges from about 25 to 40 kW/m2. This may be due to the limited availability of liquid film models included in MELCOR. Also, it is believed that due to complex test geometry, measured temperature gradients across the heat transfer plate may have been underestimated and thus the heat flux had been underestimated. The MELCOR model predicts a film thickness on the order of 100 microns, which agrees very well with film flow model developed in this study for scaling analysis. However, the expected differences in film thicknesses for near vacuum and near atmospheric test conditions are not significant. Further study on the behavior of condensate film is expected to refine the simulation results. Possible refinements include but are not limited to, the followings: CFD simulation focusing on the liquid film behavior and benchmarking with experimental analyses for simpler geometries. 16 1 INTRODUCTION This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). The experimental results are employed to validate the containment condensation model in reactor containment system safety analysis code for integral SMRs. Such a containment condensation model is important to demonstrate the adequate cooling. In the three years of investigation, following the original proposal, the following planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). The results are applicable to integral Small Modular Reactor (SMR) designs, including NuScale, mPower, Westinghouse SMR, Holtec-160 and other integral reactors with small containments of relatively high pressures under accidental conditions. Testing has been conducted at the OrSU laboratory in the existing MASLWR (Multi-Application Small Light Water Reactor) integral test facility sponsored by the US Department of Energy. Its highpressure stainless steel containment model (~2 MPa) is scaled to the NuScale SMR currently under development at NuScale Power, Inc.. Minor modifications to the model containment have been made to control the non-condensable gas fraction and to utilize the secondary loop stable steam flow for condensation testing. UW-Madison has developed a containment condensation model, which leveraged previous validated containment heat transfer work carried out at UW-Madison, and extended the range of applicability of the model to integral SMR designs that utilize containment vessels of high heat transfer efficiencies. In this final report, the research background and literature survey are presented in Chapter 2 and 3, respectively. The test facility description and modifications are summarized in Chapter 4, and the scaling analysis is introduced in Chapter 5. The tests description, procedures, and data analysis are presented in Chapter 6, while the numerical modeling is presented in Chapter 7, followed by a conclusion section in Chapter 8.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Demonstration and Verification of Thermo-Mechanical Bowing in a Limited Free-Bow SFR Concept Using MOOSE

Core bowing due to thermal gradients and irradiation induced swelling and creep introduces significant reactivity feedback effects in liquid metal-cooled fast spectrum reactors. During startup, normal operations, and transient events, differential thermal and flux gradients cause expansion, creep, and swelling, which results in bowing in preferential directions depending on local material properties and load-pad and restraint-ring design. The bowing phenomenon produces negative reactivity during accident conditions provided the restraint system has been properly designed to optimally guide the deformation outwards in the active core region. Under the Department of Energy Nuclear Energy Advanced Modeling and Simulation program, a Multiphysics Object Oriented Simulation Environment (MOOSE)-based multiphysics approach to model core bowing is being developed. The present work expands on previous modeling of simpler International Atomic Energy Agency (IAEA) verification problems with these tools by modeling IAEA Verification Problem 4, which involves a symmetric sector of a reactor core with ducted assemblies undergoing differential thermal expansion due to thermal gradients bowing outward with duct-to-duct contact. This model is available on the National Reactor Innovation Center Virtual Test Bed repository. This example was verified against benchmark participant results, which includes bowing deformation evaluation and duct-to-duct interactions at load pads with mechanical contact.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

MPACT Verification and Validation Manual (rev 1)

The MPACT code provides the capability to calculate the steady-state and transient reactor core neutron distribution for operating nuclear power plant (PWR) conditions across multiple fuel cycles. As required in the CASL SQA plan, it is the responsibility of UM and ORNL as joint owners of MPACT to ensure that verification and validation activities are performed and documented in a V&V manual with supporting publications and CASL technical reports which can be distributed for reference and distribution within VERA. This document provides the first revision of the MPACT verification and validation (V&V) manual and describes the current state of MPACT V&V and updates the framework for future MPACT V&V activities. During FY16, MPACT V&V activities in the areas of code verification have continued to mature and work has progressed on code error evaluation using the Method of Manufactured Solutions (MMS). Several specific tasks have also been completed to improve MPACT code verification to include tasks to improve the unit and regression test coverage. A joint UM/ORNL task was completed to develop standardize coding standards and workflow for the continued collaborative development of MPACT. During the past year, MPACT validation work has continued in both the areas of measured data from critical experiments as well as measured data from operating nuclear power plants. Preliminary work was also performed in FY16 on the use of post-irradiation examinations for evaluation and validation of the isotopic depletion capability in MPACT. Based on the results of V&V efforts this past year, the confidence level has increased in the ability of MPACT to model an operational Pressurized Water Reactor. Areas were identified which require increased emphasis during FY17 to include updating all the V&V problems with the new 51-group library, the completion of cycle 2 of the BEAVRS benchmark, initial validation for BWR applications, the addition of Uncertainty Quantification results into the benchmark results, and the completion of specific validation tasks for the CASL Challenge problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

System Modeling of the HTTR and Economic Dispatch Model of the Secondary System

High Temperature Gas-cooled Reactors (HTGRs) can be used for the generation of electricity and their process heat can be used to improve the efficiency of chemical processes such as hydrogen production. The JAEA-operated High-Temperature engineering Test Reactor (HTTR-GT/H2) is exploring using the reactor for electricity and hydrogen production. A RELAP5-3D model of the HTTR-GT/H2 secondary system has been developed using design information. The various components and heat exchangers in the secondary system were modeled and results were compared to the design conditions. The results for the sole-power generation mode were shown to fit the design conditions very well. The largest temperature difference was on the order of 7 K, and the largest pressure difference was on the order of 0.05 MPa. The results for the hydrogen cogeneration mode did not match the design conditions nearly as well. The largest temperature difference was about 39 K and the largest pressure difference was about 0.27 MPa at the compressor outlet. The larger differences for the hydrogen cogeneration mode are attributed to the various complex components and the flow being split in the secondary loop. A transient reduction in heat removal capability of the secondary system was investigated. Reactor temperatures are anticipated to rise as a result. The core reactivity response due to this increase in temperature is investigated and is expected to add negative reactivity to the reactor. An economic dispatch model was developed for a nuclear-driven iodine-sulfur cycle system to determine hydrogen sale prices that would make such a system profitable. The study focuses on the development of the economic model and the role that input data plays on final calculated values. It was found that the input electricity prices, whether using historical data or a host of synthetic time histories, produce significantly different breakeven hydrogen sale prices. As such, great care should be used in these economic dispatch analyses to select reasonable input assumptions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Using X-ray spectrum imaging to quantify fission gas in high-burnup UO2 fuel

Analysis of fission gas bubbles (FGBs) in light-water reactor (LWR) fuel is needed to improve the understanding and predictive capability of fuel evolution under normal- and transient conditions. One of the most important parameters of a FGB is the pressure of the gas, primarily Xe, inside it. Bubble volume and location (inter- vs. intragranular) are important considerations as well. However, experimental analysis of such bubbles is challenging due to their small size and embedded nature, and usefulness of the data requires large numbers be analyzed. This paper proposes a method to measure the pressure of Xe bubbles using X-ray spectrum imaging (XSI) in scanning transmission electron microscopy (STEM). From X-ray generation and instrumental parameters, the number of Xe atoms yielding a given number of Xe L-series X-ray counts is estimated, and then from the estimated volume of the bubble the Xe density is estimated, which can then be converted to pressure via an equation-of-state. We apply the method to XSI datasets from high burnup (HBu) fuel from North Anna 1 reactor in as-irradiated and post-LOCA-test conditions and find Xe pressures in the FGBs clustered around 1 GPa.

Parish, Chad [ORNL] (ORCID:0000000312097439)↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2017

Development of a Reduced-Order Three-Dimensional Flow Model for Thermal Mixing and Stratification Simulation during Reactor Transients

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the multi-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

Hu, Rui↗
DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Three-dimensional flow model development for thermal mixing and stratification modeling in reactor system transients analyses

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the mull-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Here, based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Integrating PCTRAN with AI-Driven Host-Intrusion Detection and Secured Container Systems for Advanced Malware Analysis (Summer Internship Report)

This study presents a solution for enhancing the security of the Personal Computer Transient Analyzer (PCTRAN) PC-based Nuclear Power Plant Simulator by integrating the software with an artificial intelligence (AI)-driven host-intrusion detection system (HIDS), in addition to a secured container system, for malware analysis. PCTRAN is a Windows XP-based software package that has the ability to simulate a variety of accident and transient conditions for nuclear power plants (NPPs). It offers a high-resolution replica of the Nuclear Steam Supply System (NSSS) and displays the status of important parameters allowing for operator interaction. By including AI-driven HIDS for the NSSS, the framework can identify security threats in real-time, ensuring the integrity of the nuclear simulation environment. Additionally, the secured container system offers the ability to isolate and analyze malware, preventing potential threats from affecting core systems. The integration process involves extensive testing and validation in order to ensure accuracy, reliability, and compliance with security policies. This framework sets a new precedent for secure simulation and training in NPP operations, and offers insight for future advancements in cybersecurity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

RELAP5-3D modelling of the HTTR-GT/H 2 secondary system and turbomachinery

High Temperature Gas-cooled Reactors (HTGRs) can be used for the cogeneration of electricity and hydrogen. The JAEA-operated High-Temperature engineering Test Reactor (HTTR-GT/H 2 ) is exploring using the reactor for electricity and hydrogen production. Here, a RELAP5-3D model of the HTTR-GT/H 2 secondary system has been developed using design information. The various components and heat exchangers in the secondary system were investigated, modified, and improved where applicable. With these changes made, the model was run, and the results were compared to the design conditions. The results for the sole-power generation mode were shown to fit the design conditions very well. The largest temperature difference was on the order of 7 K, and the largest pressure difference was on the order of 0.05 MPa. The results for the hydrogen cogeneration mode did not fit the design conditions as well as the other case. The largest temperature difference was about 39 K and the largest pressure difference was about 0.27 MPa at the compressor outlet. The larger differences for the hydrogen cogeneration mode are attributed to the various complex components and the flow being split in the secondary loop. A transient reduction in heat removal capability of the secondary system was investigated. Reactor temperatures were predicted to rise as a result. The core reactivity response due to this increase in temperature is investigated and is expected to add negative reactivity to the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Thermal Gradient Effects on Redox Evolution and Volatility-Driven Fractionation in Ternary U/Ce/Cs Condensates

Understanding how thermal history influences redox evolution and chemical fractionation is essential for characterizing high-temperature condensation in complex materials, including nuclear debris. Here, we tested the hypothesis that distinct thermal regimes in a plasma flow reactor influence redox pathways and elemental partitioning in ternary U/Ce/Cs systems. A configurable plasma flow reactor was modified with an external tube furnace to impose two distinct thermal gradients: continuous ambient cooling and a furnace-assisted thermal hold-up near 1400 K followed by rapid cooling. Transmission electron microscopy characterized phase identity, morphology, and nanoscale element distributions, while inductively coupled plasma-mass spectrometry quantified bulk elemental ratios. Across both thermal regimes, uranium and cerium condensed as UO 2 and CeO 2 as dominant refractory oxide products. Uranium partially oxidized to α-UO 3 during extended ambient cooling, while furnace-assisted hold-up preserved UO 2 and produced partial reduction of cerium to Ce 2 O 3 . Cesium remained volatile upstream and condensed later in the reactor, forming Cs 2 O and Cs-uranate phases with the highest incorporation after thermal hold-up. Bulk ICP-MS measurements supported these observations. U/Ce ratios remained comparatively stable and Cs displayed delayed and apparent transient enrichment that matched the nanoscale measurements. This integrated approach provides a quantitative method for linking thermal gradients to redox evolution and volatility-driven fractionation. These results show how the plasma flow reactor can identify where equilibrium descriptions remain adequate and where kinetic effects from residence time and temperature history must be considered when interpreting condensation behavior in multicomponent systems.

and nuclear chemistry↗