Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Serpent”

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 145 records · Page 8

Studsvik R2 materials test reactor Ad Hoc depletion strategy for the derivation of the fuel isotopic composition of the MPCMIV benchmark

The Ad Hoc Depletion Strategy elaborated by the NINE company, developed in support of the organization of the MPCMIV (Multi-physics Pellet Cladding Mechanical Interaction Validation) benchmark input and output specifications, is presented. This work aims at illustrating the strategy itself and then showing the results obtained with its application over the Studsvik R2 Testing Reactor, which is analyzed in the benchmark. The objective of the application of the strategy is to compute the fuel elements isotopic compositions at the beginning of some core loadings of interest for the benchmark. To this objective, it is necessary to implement first the simulation model of the three single assembly types and perform the infinite lattice depletions, then, to build the full core model and to perform the simulation of the core cycle. All the models and simulations were carried out with the use of the Monte Carlo particle transport code Serpent 2. Finally, the simulations results are assessed against Studsvik isotopic compositions of the fuel elements discharged from the R2 Reactor at the end of the core loading. Several assumptions were necessary during all the steps of the strategy, to overcome the lack of information regarding the core management. For this reason, the solution found at the end of the current analysis may not be completely optimized and further improvements regarding the model assumptions will be tested in a future work. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Simulation of CEFR neutronic start-up tests with FENNECS

This paper presents simulation results of selected Neutronic Start-up Tests of the China Experimental Fast Reactor (CEFR) obtained by the neutronics code FENNECS that have been performed within the frame of the IAEA Coordinated Research Program I31032. The Finite Element Neutronics code FENNECS is developed at GRS and solves the few-group steady-state and transient diffusion equation using a Galerkin-based finite element approach. Its main purpose is the safety assessment of Small Modular Reactors and Micro Reactors with complex geometry (e.g., rotating control drums) which gain increased interest internationally. Serpent has been applied to create reference models of the CEFR and cross-section libraries in 10 energy groups for FENNECS. Using these libraries, the following experiments have been simulated with FENNECS: net criticality, control rod integral and differential worth, void reactivity effect, subassembly exchange reactivity effects and reaction rate distribution. The obtained satisfactory agreements with the measurements represent a valuable contribution to the validation of FENNECS. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and verification of the pyIsoDep depletion package for nuclear thermal propulsion reactors

Nuclear thermal propulsion technology is currently being developed by the National Air and Space Administration for future manned Mars missions. There are multiple on-going activities to develop and enhance computational capabilities, single- and multi-physics, to advance the readiness level of this technology. The prediction of short-lived fission products is crucial to understand the reactivity penalty and the decay heat production during the different burnout stages of the vehicle. Generally, the solution entails coupled neutronic-depletion analysis. However, such simulations are somewhat expensive, especially if the duration of each engine burn is to be optimized. Therefore, there is a need to complement such high-fidelity techniques with time-efficient alternatives that do not require the sequential neutronic and depletion analysis. This paper details the methods and capabilities of pyIsoDep, an open-source code used for depletion analysis of nuclear thermal propulsion reactors. The package is general, in principle, and thus can be applied to any fission system, and can potentially be extended for fusion systems. The package was developed to provide a flexible and efficient alternative to pre-existing higher order depletion codes that could integrate well into mission and transient analysis of nuclear thermal propulsion reactors. Finally, this paper presents a validation and verification of pyIsoDep against analytical and Serpent-generated solutions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sensitivity analyses of a homogeneous model and a RZ model of the MYRRHA reactor in its critical configuration

Two homogenized models of the MYRRHA (Multi-purpose hybrid Research Reactor for High-tech Applications) reactor in its critical configuration are presented in this study. Both models will be verified by comparing the main reactor parameters (K{sub eff}, β{sub eff} and Λ) with those of the heterogeneous model and a sensitivity analysis will be performed. A sensitivity analysis of some important parameters related to reactivity, namely the vacuum coefficient, the Doppler and the power peaking factors, is also presented. The study of the vacuum coefficient has been carried out by analyzing two possible scenarios: 50% vacuum in the coolant (lead bismuth eutectic) and 100%, obtaining a list of the 10 reactions and nuclides that most affect the k{sub eff} value. With respect to the Doppler coefficient, it will be seen that a 300 degree increase in fuel temperature (starting at 800 K) results in a reduction of the total reactivity of the system by 134 pcm.. Finally, the map of the power peaking factors will be shown, as well as a sensitivity analysis. For making the sensitivity calculations Serpent-2 reactor physics Monte Carlo code has been used and JEFF-3.3 library has been chosen in order to continue with the OECD/NEA WPEC SG46 work evaluating the nuclear data of this library. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sensitivity and uncertainty analysis of PWR spent fuel observables to operational and model parameters

Sensitivity and uncertainty analyses of spent nuclear fuel (SNF) observables - decay heat, neutron and γ-ray emission rate - to operational and model parameters have been performed. A 2D model representing a typical PWR 17x17 UO{sub 2} fuel assembly has been taken as reference. The Serpent code and ENDF/B-VII.1 evaluated nuclear data library have been used for the analyses. Relative uncertainty of decay heat, neutron and γ-ray emission rates due to Monte Carlo counting statistics, as well as the number of fuel pin radial regions required for an accurate characterization of SNF neutron emission, have been estimated. A parametric study has been carried out to investigate the impact of the homogenization to an average value of the fuel and moderator temperature, boron concentration power and burnup on the prediction of the selected observables. In addition, linear sensitivity coefficients have been calculated. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Genetic algorithm-based optimisation of the few-group structure for lead fast reactors analysis

The optimal choice of the few-group structure for full-core transient analyses is still an open issue in reactor physics, especially for fast system like the lead fast reactor. One possible approach to select the group boundaries is represented by heuristic search algorithms, such as evolutionary ones. In this paper, a genetic algorithm coupled with the SIMMER code is employed to determine optimized six-group boundaries for the analysis of the ALFRED reactor. The Serpent Monte Carlo code is adopted to produce both the fine-group cross section library and the fine-group flux, used as a figure of merit to drive the genetic optimisation. The results show that the algorithm is indeed able to find satisfactory solutions that comply with the set objectives and can be reasonably interpreted in light of the underlying physics of the considered core. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Verification of Flux Sensitivity Estimates Using the MCNP Tally Perturbation Tool

Nuclear data is commonly used in applications such as nuclear nonproliferation, safeguards, and criticality safety. More specifically, nuclear data is used in predictive simulation codes like the Monte-Carlo N-Particle (MCNP ® ) transport code, Serpent, and similar radiation transport codes. The improvement of nuclear data enables more precise and accurate simulations, which result in higher fidelity designs and reduced operational/procedural costs. Therefore, the improvement of nuclear data is of paramount importance across the nuclear community. Nuclear data is improved and validated through integral benchmark experiments. The design of benchmark experiments is an extensive process; therefore, these experiments are often optimized on multiple characteristics, including sensitivity to the nuclear data, during the design process. Sensitivity is a measure of how much a quantity changes due to changes in independent variables such as experimental configuration. An experimental design that has a larger sensitivity to the nuclear data of interest will have a larger impact on the accuracy and precision of the validated data. Past integral benchmark experiments have primarily used the effective multiplication factor ($k_{eff}$) as the predominant measured quantity; however, experiments designed with other quantities in mind would be able to optimize on validating different areas of the nuclear data. A primary goal of the EUCLID project is to design, constrain, and reduce compensating errors in experiments focused on quantities other than $k_{eff}$ to better validate nuclear data across the board. Currently, there is a capability in MCNP to easily calculate the sensitivity of $k_{eff}$ to specific nuclear data of numerous reactions types and isotopes (KSEN card); however, the sensitivity of other quantities must be estimated in more strenuous manners. For example, the perturbation feature (PERT card) of MCNP can be used to estimate first-order sensitivities of some response in fixed source simulations. A recent announcement revealed that the first- and second-order perturbation features in previous releases of MCNP contained a bug. It was identified that particles were being scored into the wrong energy bin. The bug is in the most recent public release (MCNP6.2); however, a patch has been added to the most up to date version (MCNP6.2.2) that has not been released publicly. A direct comparison of the PERT card results for an F4 (neutron flux averaged over a cell) tally before and after the patch are shown in figure 1. All simulations used in the sensitivity estimates in this report were performed with MCNP6.2.2. This work verifies the patched MCNP perturbation tool by comparing first order sensitivities made using the PERT card to estimates made using manual perturbation of the compact ENDF (ACE) files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Depletion Analysis of Modular High Temperature Gas-cooled Reactor Loaded with LEU/Thorium Fuel

Thorium based fuel has been considered as an option to uranium-based fuel, based on considerations of resource utilization (Thorium is more widely available when compared to Uranium). The fertile isotope of Thorium (Th-232) can be converted to fissile isotope U-233 by neutron capture during the operation of a suitable nuclear reactor such as High Temperature Gas-cooled Reactor (HTGR). However, the fertile Thorium needs a fissile supporter to start and maintain the conversion process such as U-235 or Pu-239. This report presents the results of a study that analyzed the thorium utilization in a prismatic HTGR, namely Modular High Temperature Gas-Cooled Reactor (MHTGR) that was designed by General Atomics (GA). The collected for the modeling of this design come from Chapter 4 of MHTGR Preliminary Safety Information Document that GA sent to Department of Energy (DOE) on 1995. Both full core and unit cell models were used to perform this analysis using SCALE 6.1 and Serpent 1.1.18. Because of the long mean free paths (and migration lengths) of neutrons in HTRs, using a unit cell to represent a whole core can be non-trivial. The sizes of these cells were set to match the spectral index between unit cell and full core domains. It was found that for the purposes of this study an adjusted unit cell model is adequate. Discharge isotopics and one-group cross-sections were delivered to the transmutation analysis team. This report provides documentation for these calculations

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Neutronics and Fuel Performance Evaluation of Accident Tolerant Fuel under Normal Operation Conditions

This report details the analysis of neutronics and fuel performance analysis for enhanced accident tolerance fuel, with Monte Carlo reactor physics code Serpent and INL’s fuel performance code BISON, respectively. The purpose is to evaluate two of the most promising candidate materials, FeCrAl and Silicon Carbide (SiC), as the fuel cladding under normal operating conditions. Substantial neutron penalty is identified when FeCrAl is used as monolithic cladding for current oxide fuel. From the reactor physics standpoint, application of the FeCrAl alloy as coating layer on surface of zircaloy cladding is possible without increasing fuel enrichment. Meanwhile, SiC brings extra reactivity and the neutron penalty is of no concern. Application of either FeCrAl or SiC could be favorable from the fuel performance standpoint. Detailed comparison between monolithic cladding and hybrid cladding (cladding + coating) is discussed. Hybrid cladding is more practical based on the economics evaluation during the transition from current UO 2 /zircaloy to Accident Tolerant Fuel (ATF) system. However, a few issues remain to be resolved, such as the creep behavior of FeCrAl, coating spallation, inter diffusion with zirconium, etc. For SiC, its high thermal conductivity, excellent creep resistance, low thermal neutron absorption cross section, irradiation stability (minimal swelling) make it an excellent candidate materials for future nuclear fuel/cladding system.

99 GENERAL AND MISCELLANEOUS↗

Results for Phase I of the IAEA Coordinated Research Program on HTGR Uncertainties

The quantification of uncertainties in design and safety analysis of reactors is today not only broadly accepted, but in many cases became the preferred way to replace traditional conservative analysis for safety and licensing analysis. The use of a more fundamental methodology is also consistent with the reliable high fidelity physics models and robust, efficient, and accurate codes available today. To facilitate uncertainty analysis applications a comprehensive approach and methodology must be developed and applied. High Temperature Gas-cooled Reactors (HTGR) has its own peculiarities, coated particle design, large graphite quantities, different materials and high temperatures that also require other simulation requirements. The IAEA has therefore launched a Coordinated Research Project (CRP) on the HTGR Uncertainty Analysis in Modeling (UAM) in 2013 to study uncertainty propagation specifically in the HTGR analysis chain. Two benchmark problems are defined, with the prismatic design represented by the General Atomics (GA) MHTGR-350 and a 250 MW modular pebble bed design similar to the HTR-PM (INET, China). This report summarizes the contributions of the HTGR Methods Simulation group at Idaho National Laboratory (INL) up to this point of the CRP. The activities at INL have been focused so far on creating the problem specifications for the prismatic design, as well as providing reference solutions for the exercises defined for Phase I. An overview is provided of the HTGR UAM objectives and scope, and the detailed specifications for Exercises I-1, I-2, I-3 and I-4 are also included here for completeness. The main focus of the report is the compilation and discussion of reference results for Phase I (i.e. for input parameters at their nominal or best-estimate values), which is defined as the first step of the uncertainty quantification process. These reference results can be used by other CRP participants for comparison with other codes or their own reference results. The status on the Monte Carlo modeling of the experimental VHTRC facility is also discussed. Reference results were obtained for the neutronics stand-alone cases (Ex. I-1 and Ex. I-2) using the (relatively new) Monte Carlo code Serpent, and comparisons were performed with the more established Monte Carlo codes MCNP and KENO-VI. For the thermal-fluids stand-alone cases (Ex. I-3 and I-4) the commercial CFD code CFX was utilized to obtain reference results that can be compared with lower fidelity tools.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

FY 2016 Status Report on the Modeling of the M8 Calibration Series using MAMMOTH

This report provides a summary of the progress made towards validating the multi-physics reactor analysis application MAMMOTH using data from measurements performed at the Transient Reactor Test facility, TREAT. The work completed consists of a series of comparisons of TREAT element types (standard and control rod assemblies) in small geometries as well as slotted mini-cores to reference Monte Carlo simulations to ascertain the accuracy of cross section preparation techniques. After the successful completion of these smaller problems, a full core model of the half slotted core used in the M8 Calibration series was assembled. Full core MAMMOTH simulations were compared to Serpent reference calculations to assess the cross section preparation process for this larger configuration. As part of the validation process the M8 Calibration series included a steady state wire irradiation experiment and coupling factors for the experiment region. The shape of the power distribution obtained from the MAMMOTH simulation shows excellent agreement with the experiment. Larger differences were encountered in the calculation of the coupling factors, but there is also great uncertainty on how the experimental values were obtained. Future work will focus on resolving some of these differences.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerator-Driven Subcritical System for Disposing of the U.S. Spent Nuclear Fuel Inventory

The current United States inventory of the spent nuclear fuel (SNF) is ~80,000 metric tons of heavy metal (MTHM), including ~131 tons of minor actinides (MAs) and ~669 tons of plutonium. This study describes a conceptual design of an accelerator-driven subcritical (ADS) system for disposing of this SNF inventory by utilizing the 131 tons of MAs inventory and a fraction of the plutonium inventory for energy production, and transmuting some long-lived fission products. An ADS system with a homogeneous subcritical fission blanket was first examined. A spallation neutron source is used to drive the blanket and it is produced from the interaction of a 1-GeV proton beam with a lead-bismuth eutectic (LBE) target. The blanket has a liquid mobile fuel using LBE as the fuel carrier. The fuel materials are dissolved, mixed, or suspended in the liquid fuel carrier. Monte Carlo analyses were performed to determine the overall parameters of the concept. Steady-state Monte Carlo simulations were performed for three similar fission blankets. Except for, the loaded amount of actinide materials in the LBE is either 5, 7, or 10% of the total volume of the blanket, respectively. The neutron multiplication factors of the three blankets are ~0.98 and the initial MAs blanket inventories are ~10 tons. In addition, Monte Carlo burnup simulations using the MCB5 code were performed to analyze the performance of the three conceptual ADS systems. During operation, fresh fuel was fed into the fission blanket to adjust its reactivity and to control the system power. The burnup analysis shows that the three ADS concepts consume about 1.2 tons of actinides per full power year and produce 3 GW thermal power, with a proton beam power of 25 MW. For the blankets with 5, 7, or 10% actinide fuel particles loaded in the LBE, assuming that the ADS systems can be operated for 35 full-power years, the total MA materials consumed in the three ADS systems are about 30.6, 35.3, and 37.2 tons, respectively. Thus, the corresponding numbers of ADS systems to utilize the 131 tons of MA materials of the SNF inventory are 4.3, 3.7, or 3.5, respectively. ADS concepts with tube bundles inserted in the fission blanket were analyzed to overcome the disadvantages of the homogeneous blanket concept. The liquid lead is used as the target material, the mobile fuel carrier, and the primary coolant to avoid the polonium production from bismuth. Reactor physics and thermal-hydraulic analyses were coupled to determine the parameters of the heterogeneous fission blanket. The engineering requirements for a satisfactory operation performance of the HT-9 ferritic steel structure material have been realized. Two heterogeneous concepts of the subcritical fission blanket with the liquid lead mobile fuel inside or outside the tube bundles were considered. The heterogeneous configuration with the mobile fuel inside the tubes showed better performance than the configuration with mobile fuel outside the bundle tubes. The Monte Carlo burnup codes, MCB5 and SERPENT were both used to simulate the fuel burnup in the ADS concepts with the mobile fuels inside the tubes. The burnup analyses were carried out for 35 full power years. The results show that 5 ADS systems can dispose of the total United States inventory of the spent nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Generation of the Cross Section Library for PROTEUS

The cross sections for a high-fidelity neutron transport code PROTEUS can be generated in multiple ways. Multigroup cross sections in the ISOTXS format can be generated using the external codes such as MC 2 -3 and Monte Carlo codes such as Serpent and OpenMC. For the Monte Carlo codes, their multigroup cross section outputs should be converted to the ISOTXS format using a post processor named GenISOTXS developed by ANL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Transient Reactor Physics Experiment with High-Fidelity, 3-D Flux Measurements for Validation and Verification

This project led to the design and construction of a system of novel micro-pocket fission detectors (MPFDs) and their deployment as part of a series of steady-state and transient experiments at the University of Wisconsin Nuclear Reactor (UWNR). The data acquired during these experiments is intended for use in validation of reactor-physics models, including those developed as part of the NEAMS Toolkit (e.g., PROTEUS) and others. Part of the work led to a completely new computational model of UWNR based on Python-automated input generation for both MCNP and Serpent, and preliminary efforts have been made to develop a mesh for models in PROTEUS and other deterministic tools.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Benchmark Evaluation of TREAT Reactor M2 and M3 Transient Experiments

The objective of this project is to evaluate and produce TREAT reactor benchmarks under steady state and transient conditions. The benchmarks support the TREAT restart effort and are to be contributed to the International Reactor Physics Experiment Evaluation Project (IRPhEP) for inclusion in the International Handbook of Evaluated Reactor Physics Benchmark Experiments. The benchmark evaluation is performed using data generated during the “Power Calibration Experiment for test M2” and “IFR Safety Tests M2 and M3”. The evaluations are conducted using the high fidelity Serpent-OpenFOAM code system. Both the steady state and transient conditions of the M2 and M3 tests were investigated. Furthermore, TREAT has a graphitic core. Historical work has shown that the predictions of the outcomes of TREAT tests are highly sensitive to the graphite thermal neutron scattering cross sections. Therefore, in this project, thermal neutron scattering (i.e., S(α,β)) libraries of graphite, tailored to the specifics of the TREAT core, were generated and utilized in the evaluation

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

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↗

Micro Reactor Simulation Using the PROTEUS Suite in FY19

Two micro reactor cores were simulated using the MOC solver of PROTEUS: MegaPower for a heat pipe micro reactor and Holos-Quad for a high-temperature gas-cooled micro reactor. Verification of PROTEUS standalone was performed by simulating 2D and 3D core benchmark problems developed based on MegaPower and Holos-Quad. Simulation results indicated that PROTEUS successfully modelled 2D and 3D cores of the MegaPower and Holos-Quad micro reactors and PROTEUS eigenvalues agreed well with the SERPENT Monte Carlo solutions. For a heat pipe cooled micro reactor, multiphysics simulations were performed using PROTEUS / FLUENT / ANLHTP. The coupling system of the three codes were developed using Python-based drivers, which exchange the data (power, temperature, and heat rate) required for the coupling and control the boundary conditions of FLUENT and ANLHTP. The coupled system was qualitatively verified using a 3D unit assembly problem composed of six fuel rods and seven heat pipes, which was developed based on the MegaPower compositions. The steady-state problem was solved using FLUENT / ANLHTP with given asymmetric fuel powers, demonstrating the reasonable convergence of temperatures of 7 heat pipes. Transient simulations were also performed by making one out of 7 heat pipes fail. For the transient, the coupled system of PROTEUS / FLUENT / ANLHTP demonstrated reasonable changes of total power and heat pipe temperatures with time, accounting for temperature feedback effects. In this simulation, no effect of structural change was accounted for, which will be performed using ANSYS-mechanical in future along with further verification with larger or whole-core 3D problems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗