SEARCH · Engineering Papers

Results for “NEAMS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Development of Critical Experiments to Provide Validation Data for Multiphysics Coupling Methods (Final Report)

Computational simulation has become one of the most important tools in the modern toolbox for nuclear reactor design. Reactor simulations, like the reactors themselves, involve multiple physical processes including neutronics, fluid flow and heat transfer (thermal-hydraulics), thermo-mechanics, fuel behavior, chemistry, and balance of plant. Each process can affect the behavior of the others, and advanced computational tools have been developed that couple separate physics models together to simulate the feedback effects. The objective of this work is to generate experimental data to support the validation of multi-physics reactor codes. It specifically targets the SHARP toolset from the NEAMS Reactor Product Line, but it is intended to be useful in any set of coded that couples thermal-hydraulics with neutronics simulations. The project takes advantage of the inherent flexibility of the RPI Reactor Critical Facility (RCF), a low power open-pool research reactor. The standard configuration of the RCF is used to explore the effects of moderator temperature on neutronic behavior and reactivity. Because the low power of the reactor generates no appreciable temperature change in the moderator, the moderator temperature is artificially changed with two 18 kW electric immersion heaters. The temperature feedback effects can be observed by recording the change in the reactor power as the reactor goes from a supercritical state, to critical, to a subcritical state. To expand the temperature range and temperature rate of change that can be investigated, an experimental apparatus consisting of a heated water loop was designed and installed in the RCF. In this setup, water from a heated reservoir tank is pumped through a pipe placed through the center of the reactor core. The water in the reservoir can be heated to around 70 °C and pumped through the system. The pump has a variable flow rate and electronic control, so a number of scenarios can be developed. A slug of water can be pushed through the system to generate a reasonably fast change in reactivity. The pump speed can be increased or decreased over the course of a measurement, or even oscillated at different periods, to observe the response of the reactor. To complement the experimental data generated by the RCF, a set of computational models have been developed. They are intended to illustrate the usefulness of the experimental data generated in validating coupled multi-physics simulations. A model previously constructed in MCNP has been used in conjunction with safety analysis at the RCF for many years and was modified to support the design of the heated water loop experiment. This model was translated into Serpent 2, a similar Monte Carlo code, which provided the added advantage of being able to produce collapsed cross-section data files for use in the deterministic neutron transport codes that are part of the SHARP toolset. A complete reconstruction of the RCF model, both with and without the heated fluid loop experiment, were necessary for use with the SHARP codes, PROTEUS (neutronics) and Nek5000 (thermal hydraulics), as these require an unstructured mesh format to operate. The computational requirements for fully coupled three dimensional simulations are significant, and challenges to full implementation and comparison with RCF experiments remain. A wide variety of experiments have been conducted for eventual benchmarking and validation of multi-physics reactor codes. A complete database of detailed time- and temperature-dependent measurements have been generated and stored. Several examples of the experimental results are presented here, with refinement and analysis of experiments ongoing. While these experiments still only address a finite range of operating conditions, in particular being limited to water in the sub-boiling regime, they represent a unique (to date) look at the feedback between temperature and neutronic behavior suited for validation of the coupling routines in multi-physics codes. It is hoped that these will be the foundation of future benchmarks that will provide confidence in the tools that will be used to develop the next generation of nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

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 develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The DSA scheme has been implemented in PROTEUS to speed up SN transport calculations. The existing DSA scheme employs a consistent finite element formulation on the same fine mesh structure as the SN solution While the use of consistent discretizations makes the DSA method effective, it does not necessarily make it efficient because the numerical solution of a discretized elliptic diffusion problem itself can be costly, particularly when the problem size becomes large. We developed a new discontinuous Galerkin (DG) discretization of the diffusion equation, called DG-DSA, which can effectively and efficiently accelerate the SN transport iterations. As compared with the previous work, the novelty of our method is that the diffusion equation is solved on a coarse-mesh grid using the DG methods, and the DG diffusion discretization incorporates local hp adaptation, i.e., local adaptation of mesh size and/or polynomial degree, based on local total cross section (or optical thickness). Therefore, the resulting number of degrees of freedom (DOF) of the DG discretization is much less than the conventional consistent DSA discretizations, and thus DG-DSA can achieve significant improvement in computational efficiency. We implemented this scheme in PROTEUS-SN. In addition, we developed a new stabilization scheme named linear prolongation CMFD (lpCMFD). A novel feature of this scheme is that the conventional flat flux ratio–based scaling approach is replaced with a linear interpolation of the scalar flux differences at the coarse-mesh cell edges between the neutron transport and CMFD calculations. Fourier convergence analysis and numerical results show that lpCMFD is unconditionally stable and more effective than other CMFD based on acceleration schemes such as pCMFD, odCMFD. We also developed a new nonlinear diffusion acceleration scheme for solving neutron transport equations. This scheme, called LR-NDA, employs a local refinement approach on the framework of CMFD by solving a local boundary value problem of the scalar flux on the coarse-mesh structure to replace the piecewise constant scalar flux obtained by CMFD. The refined flux is then used to update the scalar flux in the neutron transport source iteration. It has been demonstrated that that LR-NDA is much more effective and stable than CMFD for a wide range of optical thicknesses. LR-NDA is a local adaptive method, which means LR-NDA does not necessarily require local refinement for all the coarse-mesh cells on the problem domain, i.e., it can be used only for relatively optically thick regions where the standard CMFD scheme would encounter the convergence problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Development of Accident Tolerant Fuel Options For Near Term Applications

In this section, the basis for the proposal followed by description of the work scope and methodologies are discussed. Following the Fukushima disaster in 2011, US congress mandated the department of energy to start a R&D program on accident tolerant fuel (ATF) concepts for the existing light water reactor fleet. By 2015, the start date of this IRP, several ATF concepts were being pursued by various entities. In general, the ATF concepts were divided into two categories: near term and long term. The near term concepts included: Coated Zircaloy clad, fuels with additives and dopants and FeCrAl/steel based claddings. The long term concepts included: SiC composite cladding, high density fuels (U 3 Si 2 , UN) and TRISO type fuel forms (e.g. FCM). It is of particular importance for any DOE funded effort to give background information and original Funding Opportunity Announcement (FOA) goals in its final report to evaluate whether the FOA goals and the proposed work were actually addressed. The FOA for this IRP sought for modeling and simulation capability to predict various ATF performances during normal operation, design basis accident and severe accident conditions. Particularly, the FOA was concern with the state of the core performance during a severe accident. It was postulated that if an ATF retains its structural integrity, the performance of non-fuel structures in the core needs to be addressed. The IRP team focused the proposal on Pressurized Water Reactor (PWR) performance for near term ATF concepts: Coated cladding, Fuel with additives and dopants to accelerate adoption of such concepts and address the unanswered questions regarding their performance. The project also gave limited focus on another near term concept, FeCrAl cladding, by restricting its scope to time-to-failure prediction since substantial effort supported by other programs including ORNL ATF program, GEH ATF program and NEAMS High Impact Problem were already made or underway during the IRP tenure. Despite two previous IRPs funded in 2012 with focus on experimental study of metallic and ceramic Zircaloy coated claddings, this IRP still included a coating program. This was a testament to the project leadership’s belief that the coated clad concept will be the primarily initial pathway for near term industrialization of ATF technology due to presence of enrichment penalty and licensing hurdles of the alternative clad concepts. In doing so, the IRP was positioned to provide key findings on coated Zircaloy as a uniquely funded external program to DOE’s industrial ATF campaign. Currently, in 2019, the coated clad concept is being pursued by all three vendors for licensed application and Lead Test Rods and Lead Test Assemblies (LTRs/LTAs) have already been inserted or planned to be inserted in the existing fleet. At the start of the IRP project, based on the meeting with the projects technical and federal point of contacts, it was also agreed upon if the contractual obligations are met, then the scope of the simulation work will be extended to Boiling Water Reactors as included in Section 2.10.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

ATF material model development and validation for priority fuel concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the Bison fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on the focus was on U 3 Si 2 fuel and FeCrAl cladding as part of a high impact problem through the Nuclear Energy and Advanced Modeling Simulation (NEAMS) program. Then, developments for Cr 2 O 3 -doped UO 2 fuel, and SiC-SiC and Cr-coated zirconium-based claddings began based upon industry interests. In late fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors (CASL) took over further ATF work in Bison in support of the Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC identified their list of priority fuel and cladding concepts, which included Cr 2 O 3 -doped UO 2 and U 3 Si 2 fuels, and Cr-coated zirconium-based and FeCrAl claddings. In particular, the NRC suggested that reports similar in form to NUREG/CR-7024 [1] that was developed for traditional LWR materials UO 2 and zirconium-based claddings (i.e., Zircaloy-4, M5®, ZIRLO) be created for the priority ATF concepts. The approach to ATF capability development in Bison since the beginning has been two-fold: (1) empirical correlations and (2) multiscale model development. Both approaches have uncertainty inherent to them. Uncertainty in empirical correlations is bounded by the experimental data upon which with the correlation was developed. Models developed through a multiscale approach have uncertainty associated with the lower length scale calculations and input parameters that must be propagated to the engineering scale model in Bison. In this report, the recommended models, their range of applicability (e.g., temperature, burnup), and associated uncertainty for the NRC priority fuel concepts Cr 2 O 3 -doped UO 2 and U 3 Si 2 are presented in a manner similar to the aforementioned NUREG. In addition, the Cr 2 O 3 -doped UO 2 models are validated to the Halden Reactor tess IFA-677.1 rods 1 and 5 and IFA-716.1 rod 1. For U 3 Si 2 the models are validated to the recent post-irradiation examination (PIE) data ATF-13 and ATF-15 rods that were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory. Finally, an uncertainty quantification (UQ) and sensitivity analysis (SA) is completed on a single rod for each fuel concept that takes into account the defined uncertainty in select models. The results of the UQ and SA analyses indicate that given the large uncertainty in some of the input ATF models that the uncertainty on the Bison predictions of fuel performance metrics of interest bound the available experimental data. Correlation coefficients are also reported that identified the uncertain inputs with the strongest correlation (positive or negative) with the outputs of interest. Fuel performance metrics investigated include fuel elongation, rod internal pressure, fuel centerline temperature, and fission gas release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2018

Coupled Calculations of the Nuscale FIV Tests: Pre-test Simulations and Sensitivity Analysis

42 ENGINEERING↗

DOE OSTI2019

Demonstration of Fully Coupled Calculations in Helical Steam Generator: Toward Predictions of Fluid Elastic Instability

42 ENGINEERING↗

DOE OSTI2018

Demonstration of Fully Coupled Calculations in Helical Steam Generator: Toward Predictions of Fluid Elastic Instability

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. In addition to being more cost- and space-efficient, such steam generators need to be more reliable, since failure of the primary vessel represents a potential loss of coolant and a safety concern. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict their occurrence based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the design stage. Over the past seven years, the Reactor Product Line has developed the integrated multi-physics code suite SHARP. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside the full reactor core or portions of it with a user-specified fidelity. In particular SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. These tools form a strong basis on which to build a flow-induced vibration modeling capability. In the last two years a significant effort was devoted to 1) develop and validate a fully coupled capability in SHARP for Fluid-Structure calculations (FSI) and 2) validate the use of SHARP for one-way coupled calculations in helical steam generators on legacy data available as well as novel experiments. In this report we discuss an extension of previous one-way coupled calculations performed with Nek5000 and Diablo aimed at simulating available FIV experiments in helical steam generators in the turbulent buffeting regime. It was demonstrated that in the buffeting regime one-way coupling was judged sufficient because the pressure loads do not cause substantial displacements. However, higher speeds led to a reduction in accuracy and the need for a fully coupled capability for velocity or approaching the critical velocity and fluid elastic instability. In this report we demonstrate such fully coupled capability on the Argonne test case discussed in previous reports.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Coupled Calculations of the NuScale FIV Tests – Pre-test Simulations and Sensitivity Analysis

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict FIV based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the conceptual design stage. Over the past several years, the Reactor Product Line has developed a FIV capability in SHARP. SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. In previous reports we have discussed the implementation and validation of an FIV capability for helical steam generators and fuel assemblies. In this report we discuss one-way coupled calculations performed with Nek5000 and Diablo for the NuScale FIV tests. In the turbulent fretting regime one-way coupling is judged to be sufficient because the pressure loads do not cause substantial displacements. It is also the most common source of vibration in helical steam generators per the low flows expected in integral PWRs. The simulations precede the tests, as experimental data was not available at the delivery of this report. Several assumptions had to be made on the support-tube interaction and multiple tests and sensitivity analysis have been performed. Additional work remains to be done and it is discussed in the future work section.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Status Assessment of the PROTEUS-ORIGEN Coupling

As part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign of the U.S. Department of Energy, Office of Nuclear Energy, an initiative was started to provide time dependent isotopic modeling capabilities to the PROTEUS neutronics module of the Simulation based High-efficiency Advanced Reactor Prototyping (SHARP) toolkit. This capability is to be provided through a coupling with the Oak Ridge Isotope Generation (ORIGEN) code, which is part of the SCALE code package developed and maintained at Oak Ridge National Laboratory. This development has been a joint effort between the developers of PROTEUS at Argonne National Laboratory (ANL) and Oak Ridge National Laboratory (ORNL).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Multi-Physics Simulations of Heat Pipe Micro Reactor

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 presents the development of new capabilities for modeling the heat pipe type reactor systems. The need for power at remote locations away from a reliable electrical grid is an important niche for nuclear energy. Heat pipe-cooled fast-spectrum nuclear reactors are well suited for these applications. The key feature of the heat pipe reactors is the use of heat pipes for heat removal from the reactor core. The heat pipe makes use of the phase change of the working fluid and transports a large amount of heat from the evaporator to the condensation end with very small temperature drops. In contrast to the traditional nuclear reactor system that makes use of pumped loop for extracting the thermal power, the heat pipe reactors make use of hundreds of heat pipes for removing the thermal power (including the decay heat) passively. This could potentially significantly improve the reliability and safety of the reactor systems. The essential part in the analysis of a heat pipe type reactor is the modeling of heat transport inside the heat pipe. The capability of SAM was extended to enable the modeling of the conventional heat pipe. The full suite of non-LWR codes for confirmatory analysis at NRC is known as the Comprehensive Reactor Analysis Bundle (CRAB). It makes use of existing NRC codes, and integrates them with several codes developed through the DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this report, the heat pipe micro reactor is simulated using three submodules under CRAB. The MAMMOTH module is used to simulate the reactor kinetics behavior of the micro reactor; the SAM module is used to simulate the heat conduction in the reactor core and heat removal through the heat pipe heat exchangers and Reactor Cavity Cooling System (RCCS); and the MOOSE Tensor Mechanics module is used to simulate the thermal expansion of the reactor cores. The different sub-models are coupled together using MOOSE’s MultiApp system and executed using the BlueCrab application. The new capability has been demonstrated by a steady-state operation analysis, a failure of a single central heat pipe transient analysis, and a loss of heat sink transient analysis. The fully coupled model is shown to work well.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

Development and Validation of a Conjugate Heat Transfer Model for the Two-Phase CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance CFD code NEK5000. The goal of the project is to develop a new two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) [2-3] previously developed at Argonne, but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2019. This report focuses on the FY2019 development and validation of the Conjugate Heat Transfer (CHT) model for the NEK-2P Two-Phase, Computational Fluid Dynamics (CFD) code. The NEK-2P CHT model calculates the coupled behavior of the fluid and solid domains. In the fluid domain NEK-2P calculates the vapor and liquid phase temperatures, velocities and mass fractions. In the solid domain only the temperatures are calculated. The fluid and solid domains are coupled with a heat flux boundary condition at the fluid-solid interface. The heat flux and the wall temperatures at the fluid-solid interface are updated using the NEK-2P wall heat-flux partitioning model. A wall heat flux boundary condition is also applied to the exterior surface of the solid domain. The CHT model validation was illustrated through the analysis of three of the Becker Critical Heat Flux (CHF) tests. Reasonably good agreement with measured data was obtained in predicting the CHF location and post CHF wall temperature behavior illustrating the ability of the NEK-2P CHT model to simulate the CHF phenomena for a wide range of thermal-hydraulic conditions. Apart from the CHT model development and validation, initial simulations were performed for the Virginia Tech. (VT) 3x3 bare rod bundle experiments including both cold and boiling flow simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

SAM Code Validation using the Compact Integral Effects Test (CIET) Experimental Data

The primary objective of this study is to validate the system analysis code, SAM, using experimental data from the Compact Integral Effects Test (CIET) facility. The System Analysis Module (SAM) is an advanced and modern system analysis tool being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Originally aiming at supporting sodium fast reactor (SFR) design, it is now extended to cover other non-light water reactor (non-LWR) concepts, such as high-temperature gas-cooled reactor (HTGR) designs, molten salt-fueled, and molten salt-cooled reactor designs. To support SAM code development for the wide range of non-LWR applications, it is of paramount importance to validate the code against experiments relevant to these reactor concepts. The Compact Integral Effects Test (CIET) facility, which was designed to reproduce the thermal-hydraulics response of fluoride salt-cooled high-temperature reactors (FHR) under both forced- and natural-circulation conditions, has been identified and selected as one of benchmark test facilities for SAM code validation. Among available CIET experiments, three sets of them with distinctive characteristics were selected for SAM code validation purpose. The three sets experiments are: 1) Power step change transient tests; 2) DHX-DRACS natural circulation tests; and 3) heater frequency response tests. For all tests, steady-state or transient, 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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Multi-Scale Modeling of Thermal-Fluid Phenomena Related to Loss of Forced Circulation Transient in HTGRs

Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS), an effort had been pursued to support High Temperature Gas-cooled Reactors (HTGR) technology development and its modeling and simulation needs. There was a particular need for advanced modeling & simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially the core and lower and upper plena, during safety-related transients. In this report, two main such activities were presented: 1) system level HTGR simulations using system analysis code SAM for both normal operating and accident conditions; and 2) three-dimensional (3D) computational fluid dynamics (CFD) simulations on Deteriorated Turbulent Heat Transfer (DTHT) of interest to HTGR transients, including code benchmark and closure correlation development. In this work, Modular High Temperature Gas-cooled Reactor (MHTGR) full loop simulations were performed using System Analysis Module (SAM) code. Two approaches were used to model the reactor core; (i) unit-cell or lumped parameter model, and (ii) detailed ring model, to address the different needs under different reactor conditions. Aside from the core model changes other reactor components in the primary system, such as upper/lower plenum, heat exchanger and blower etc, were kept the same. In the unit-cell core model, the lumped parameter approach was used to model the coolant channels and the heat structures using 1-D and 2-D components, respectively. The fuel assembly heat structure was modeled based on unit-cell thickness approximation. In normal and controlled operating conditions, the calculated maximum temperatures for both coolant and heat structures were within the designed limits. The unit-cell model is, however, not able to capture the radial heat conduction through fuel assemblies, which is the dominant heat transfer mechanism during accident transients, such as pressurized conduction cooldown (PCC) event. To address this modeling challenge, ring models were developed. A simplified ring model was developed first. In this simplified model, the active core was simulated with 9 circular rings consists of 6 rings for homogenized fuel heat structure and 3 rings for gas coolant. It was quickly found that this model gives rather large simulation error as it was not able to correctly capture the solid structure representative geometries. A detailed ring model was developed to address this issue. The active core was simulated with 99 circular rings consists of 66 rings for homogenized fuel heat structure and 33 rings for gas coolant. In both ring models, six additional rings were included to represent inner reflector, outer reflector, core barrel, Reactor Pressure Vessel (RPV), RPV coolant channel, and Reactive Cavity Cooling System (RCCS). The detailed ring model was applied in two different SAM simulations; (i) normal operating conditions, and (ii) accident scenarios (more specifically, PCC event), for MHTGR design. In both simulations, the SAM calculated results for coolant and heat structure temperatures were within the reactor designed limits. The newly proposed detailed ring model approach in this work has proved to be a promising method to address the challenging modeling and simulation needs in the MHTGR design. As demonstrated in this work, this model was suitable for MHTGR simulations under both normal and accident conditions. Apart from the SAM model development for gas-cooled reactors, 3D Computational Fluid Dynamics (CFD) models were developed and validated with relaminarized Deteriorated Turbulent Heat Transfer (DTHT) MIT benchmark test. Two different CFD codes, Nek5000 and STAR-CCM+, were used for this benchmark activity. Compared with MIT data, Nek5000 shows very good agreement between its simulated and experimentally measured wall temperatures. In contrast, STAR-CCM+ under-predicted wall temperature due to over predicted turbulence in the wall region. The Nek5000 code was also used to develop wall-heat transfer correlations for laminar flow in cylindrical tube. Several simulations were performed for various Reynolds flow and wall heat fluxes. Significant deviations were found when compared the predicted heat transfer coefficient and the Sieder-Tate correlation. A new set of heat transfer correlation was proposed for laminar flow which is valid for 400<Re<2000.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Updates and Verifications of the PROTEUS Suite in FY19

In FY19 efforts were made on updating and verifying the PROTEUS suite in terms of depletion and fuel cycle capabilities, updated MSR capabilities, improving code usability, and supporting simulation of micro reactors. The fuel reloading and shuffling schemes as well as the built-in depletion solver were implemented to add the fuel cycle capability to PROTEUS-NODAL. PERSENT can run with PROTEUS-NODAL which now produce the CCCC interface files. The MSR simulation capability was further updated by making PROTEUS-NODAL coupled with SAM via the MOOSE wrapper. The high fidelity solvers of PROTEUS were updated to support the multiphysics simulation of SFRs and micro reactors. The coupled system of PROTEUS / FLUENT / ANLHTP was developed using Python-based external drivers, demonstrating the coupled simulation for small 3D test problems with heat pipe failure transients. In addition, the gamma transport capability was implemented to compute the detailed gamma distributions that determines the primary heat sources for the fuel assembly structure and non-fueled assemblies. Verification tests using fuel assembly and 2D core problems indicated very good agreement in heating distributions between PROTEUS and VARIANT / GAMSOR solutions. For user support, the NEAMS Workbench was updated to include PROTEUS-NODAL and further extended to handle PROTEUS-MOC for which the support of Workbench is limited, instead requiring the off-line generation of mesh files and corresponding cross section sets. The utility tools of PROTEUS were also updated to support users to model and simulate complicated and non-conventional geometry cores.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Summary of Grizzly Development for Advanced Reactor Structural Materials

This summary report contains an overview of work performed during Fiscal Year 2019 under the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) Fuels Product Line work package entitled "MS-19IN020106 - GRIZZLY - INL". Grizzly is a code based on the MOOSE framework that has been under development for several years funded by the US Department of Energy’s Light Water Reactor Sustainability (LWRS) program. Prior to 2019, Grizzly development was focused on development of capabilities needed to simulate the effects of aging mechanisms in light water reactor (LWR) components. The goal of this LWR-focused work has been to develop tools that can be used to assess the ability of existing LWRs to perform safety in long term operation (LTO) scenarios. This LWR-focused work has primarily focused on the effects of reactor pressure vessel (RPV) embrittlement and aging of reinforced concrete structures. Grizzly development to support LWR degradation needs is still ongoing, but starting in 2019, Grizzly’s capabilities are being expanded to also meet growing needs for a tool to address structural component integrity and degradation in advanced reactors. A simulation capability for advanced reactor components has many of the same basic requirements that one for LWRs does, such as the ability to solve coupled systems of partial differential equations that arise from the fact that many of the problems of interest involve multiple coupled physics. The materials issues also inherently involve multiple length scales, with the response at the engineering scale being strongly affected by mechanisms occurring at atomistic and mesoscopic scales.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

One-Way Coupled Flow-Induced Vibration Analysis for a Twisted-Tube Heat Exchanger

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

A Coupled Experimental and Simulation Approach to Investigate the Impact of Grain Growth, Amorphization, and Grain Subdivision in Accident Tolerant U 3 Si 2 Light Water Reactor Fuel

Triuranium disilicide (U 3 Si 2 ) is being considered as an alternative accident tolerant fuel (ATF) due to its higher thermal conductivity. However, there is uncertainty in its use due to a lack of irradiation data at light water reactor (LWR) conditions. In this project, we used an integrated experimental and simulation approach to answer two significant questions pertaining to U 3 Si 2 that were not investigated by the Fuel Cycle Research and Development (FCRD) and Nuclear Energy Advanced Modeling and Simulation (NEAMS) programs: (1) Will grain growth in the hotter portions of the fuel significantly impact U 3 Si 2 LWR fuel behavior? (2) Under what conditions do grain subdivision and amorphization occur in U 3 Si 2 and will either occur at LWR conditions? Our approach to answer these questions used in situ ion irradiation and annealing to investigate grain growth, grain subdivision, and amorphization, along with mesoscale simulations using the MARMOT tool. We found that while the grain boundary mobility of U 3 Si 2 is larger than that of UO 2 , grain growth in U 3 Si 2 fuel pellets will be less than in UO 2 and may be neglected due to the lower thermal conductivity and thus lower centerline temperature. We also found that amorphization will not occur in U 3 Si 2 above 600 K, such that it is not likely to occur at typical LWR operating conditions. Grain subdivision will occur at LWR conditions and will likely occur at lower burnups than in UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Demonstration of MSR salt property database operation with a reactor analysis tool

Successful modeling and simulation of molten salt reactors (MSRs) requires quality thermodynamic and thermophysical data for the working fluid. A goal of the Molten Salt Reactor Campaign is to develop a quality-controlled database of these properties that will be available for a multitude of different salt types. Furthermore, properties shall be defined in a standardized format that can then be used by reactor modeling and analysis tools. This report details work done to develop an approach for integrating salt properties into a MSR analysis tool. This activity required defining a format for thermophysical properties that could be produced from salt data delivered by Argonne National Laboratory (ANL) and which could be read into a reactor analysis tool. The reactor analysis tool used for this activity was the Consortium for Advanced Simulation of Light Water Reactors-developed tool, VERA. The changes implemented in VERA to accept the thermophysical properties file are discussed and tested and the process of reading thermophysical properties into VERA for a MSR simulation is demonstrated. In addition to being able to get fluid properties into the simulation tool, it is also necessary that a process exists for adding new salt data to the salt database. The NEAMS Workbench is being expanded with a new package called ThermPhys that will provide a user interface for adding new salts to the database while performing input validation and ensuring data quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗