Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Full core”

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 109 records · Page 6

Neutronics and Thermo-Fluids Simulation of Generic Pebble-Bed Fluoride-Salt-Cooled High-Temperature Reactor

The fluoride-salt-cooled high-temperature reactor (FHR) is one of the advanced reactors that has been attracting considerable interest from both the research community and the nuclear industry. To help facilitate the nuclear community's familiarity with the FHR, Kairos Power has developed a generic FHR (gFHR) benchmark. In the research performed here, this benchmark was used to assess innovative modeling methods that combine stochastic and deterministic computer codes to perform the design and analysis of the gFHR. Further, the Monte Carlo code Serpent 2 was used to generate few-group cross sections that were then used in the neutron diffusion and thermal-fluids code AGREE to perform full-core neutronics and thermal-fluids steady-state and transient core analysis. The Argonne National Laboratory code SAM was then used to model the gFHR system and to simulate the load-follow operation of the gFHR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DUCTED ASSEMBLY STEADY STATE HEAT TRANSFER SOFTWARE

The Ducted Assembly Steady State Heat transfer software (DASSH) performs full core, steady state, subchannel temperature calculations for hexagonal, ducted assemblies. DASSH is targeted for use in the early stages of the reactor design process when assembly characteristics may not be defined or are undergoing considerable changes. DASSH provides a relatively rapid assessment of the temperature and flow distribution to inform reactor design updates and analysis.

ATZ, MILOSIVO↗

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

VERA BWR progression problems

During the first phase of the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, the Virtual Environment for Reactor Applications (VERA) was developed with a focus on capabilities for high-fidelity, multiphysics simulation of pressurized water reactors (PWRs). During this development effort, a set of progression problems was created ranging from smaller pin cell calculations to larger 3D full-core calculations. These progression problems helped to guide the development of the software and served as benchmarks against which to test VERA. Since 2019, efforts have been made to extend the capabilities of VERA to model boiling water reactors (BWRs). Because BWR simulations come with many unique challenges, a set of BWR progression problems was developed to aid in this new effort. The BWR progression problems range from 2D lattice calculations to 3D mini-core problems, and reference neutronic solutions were computed using continuous-energy Monte Carlo codes. MPACT, one of the neutronics code in VERA, was benchmarked using the BWR progression problems. The code is capable of computing solutions to all problems. The eigenvalues computed by MPACT agree well with the Monte Carlo reference solutions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Overview of recent SCALE activities for Non-LWR inventory and decay heat analysis

In 2019, the US Nuclear Regulatory Commission initiated a project for the development and assessments of non-light-water reactor (non-LWR) accident progression using the SCALE and MELCOR simulation tools. SCALE simulations are used to generate nuclide inventories, full-core power distributions, decay heat, and kinetics parameters to initialize MELCOR simulations of severe accident scenarios. Five non-LWR concepts were studied: high-temperature gas-cooled reactor (HTGR), heat pipe reactor (HPR), high-temperature fluoride salt-cooled reactor (FHR), molten salt-fueled reactor (MSR), and sodium-cooled fast reactor (SFR). This paper summarizes the SCALE results obtained in 2021 for the first three non-LWR concepts, compares characteristics and results to common LWRs, and provides the strategy for the analysis of the remaining two non-LWRs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

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↗

As-Run Neutronics Evaluation for the CSM-10584 Experiment in the ATR

This engineering calculations and analysis report (ECAR) documents the results of the Advanced Test Reactor (ATR) detailed Monte Carlo N-Particle (MCNP) code full-core model as-run physics analysis performed to support the Colorado School of Mines (CSM) experiment in the B-5 position.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ECAR 4951 BSU-8242 3 DPA As-Run Physics Analysis

The purpose of this report is to document the as-run physics analysis for the Boise State University (BSU-8242) experiment in the A6, A7 and A8 positions, correlating with the 3 DPA specimens irradiated in Cycles 164A, 164B, and 166A. In this paper the heat rate, flux, fluence, and DPA for all specimens are scaled to the appropriate as-run lobe powers for an accurate assessment and description of the irradiated materials. Source terms per capsule are provided, and a bounded source term per specimen type is given on a per gram basis. Verified and validated MCNP full-core models were used for the neutronics analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The STAT7 Code for Statistical Propagation of Uncertainties In Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

STAT7 was written to automate many of the steady-state thermal hydraulic safety calculations for the MIT research reactor, both for conversion of the reactor from highly enriched uranium fuel to low-enriched uranium fuel and for future fuel re-loads after the conversion. A Monte-Carlo statistical propagation approach is used to treat uncertainties in important parameters in the analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to ONB for each axial node of each stripe of each plate of each element in the core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Detailed analyses of a TRISO-fueled microreactor: Modeling of a Micro-Reactor System using NEAMS Tools

Our analysis provides preliminary confirmation of the feasibility of applying the NEAMS codes for multiphysics simulations to model load following transients and accidental scenarios. Detailed demonstration was performed on a simple unit-cell model to perform transient simulations, while verifying accuracy of the multiphysics solution obtained with the NEAMS codes. Results on the full core focused on showcasing the multiphysics workflow rather than detailed verification of the results, at least for this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Fabricating Fuel for the Versatile Test Reactor

A metal driver fuel has been proposed for the Versatile Test Reactor (VTR). About 30 years ago, the Experimental Breeder Reactor-II (EBR-II) was the last reactor in the U.S. to utilize a full core of metal driver fuel. While the necessary knowledge to make metal fuels is well preserved and is practiced for research activities today, re-establishing a production line to support the fuel needs of a 300 MWth reactor has unique technical and engineering challenges. These challenges are the focus of a multi-laboratory and private sector team that has been tasked with the responsibility to fabricate fuel for the VTR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Implementation of an Orificing Optimization Algorithm in the DASSH Subchannel Analysis Code

The Ducted Assembly Steady-State Heat transfer code (DASSH) performs full-core subchannel thermal hydraulics calculations in liquid metal fast reactors. One of the applications of subchannel codes is to optimize coolant flow orificing. As a design activity, the primary task is to determine the best way to divide assemblies into groups and distribute coolant flow rates among them. This report documents an algorithm implemented in DASSH to automatically optimize coolant orificing. Over the course of multiple iterations, DASSH determines the orifice grouping and flow distribution that minimizes peak coolant, clad, or fuel temperatures across all timesteps for a user-specified number of assembly groups. The total coolant flow rate in the reactor is constrained to achieve the specified core-average outlet temperature. The flow rate to each orifice group may also be constrained by the allowable pressure drop. The distribution of coolant flow among groups is accelerated using a predictor-corrector algorithm based on interpolated results from single-assembly parametric calculations. The assembly orificing grouping is initially predicted based on assembly power but can be refined if results demonstrate that an assembly would fit better in another group. The algorithm is demonstrated with two case studies. The first is a simple model for a reactor core consisting of just fuel assemblies; the pin power distributions are specified to create a situation where the initial assembly grouping prediction is suboptimal. This example is used to describe the initial grouping, demonstrate convergence over multiple iterations, and highlight the impact of regrouping. Then, the algorithm is applied to minimize peak clad and fuel temperatures in an example sodium-cooled fast reactor, the Versatile Test Reactor. The multicycle optimization confirms prior calculations for the reference core design. The example highlights how optimizing for different peak temperatures affects the results and demonstrates the use of the pressure drop constraint to limit the maximum flow rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SCALE 6.2.4: Validation Overview

SCALE is a multi-application code system with tools for reactor physics, criticality safety, radiation shielding, and spent fuel characterization for a range of nuclear systems. The last of the SCALE 6.2 series, version 6.2.4, was released in June 2020. Validation studies in nuclear criticality safety, reactor physics, and shielding applications were performed to confirm results, assess accuracy, and identify biases from SCALE simulations. Nuclear criticality safety applications yield small code biases for a wide range of systems, generally overestimating cross section uncertainty. Reactor physics applications improvements in nuclide inventory and full-core eigenvalue predictions, with similar performance in decay heat predictions. Radiation shielding applications yield sufficient agreement in comparisons to measurements and alternate simulations for eight shielding experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The STAT7 Code for Statistical Propagation of Uncertainties in Steady State Thermal Hydraulics Analysis of Plate Fueled Reactors

The STAT7 software was developed to perform steady-state thermal hydraulic analyses. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as MITR-II. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model common fabrication parameters and other key reactor analysis uncertainties required for research and test reactor thermal hydraulic analyses. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one lateral stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each element in the core.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of Lateral Conduction in a Rectangular Fuel Plate and Optimum Stripe Widths for the OPAL Reactor

The impetus for this work is the Proliferation Resistance Optimization (PRO-X) program, which was created by the National Nuclear Security Administration (NNSA) of the US Department of Energy (DOE) to provide a framework for developing reactor designs that minimize the production of special nuclear materials while maximizing performance for peaceful uses. The Argentine company INVAP and Argonne National Laboratory collaborated to share knowledge and expertise regarding the thermal analysis of research reactors. Working together, they designed a series of test problems based on the Open Pool Australian Light Water (OPAL) reactor, whose fuel is in the form of flat plates cooled by water flowing through narrow rectangular channels. ANL used their PLTEMP/ANL code. INVAP used their TERMIC code (and its new multi-plate version, TERMIC-MP). Then, INVAP and ANL compared the two sets of results. 1-D models were used during the collaboration since for research reactor analysis, largely 1-D thermal-hydraulics models are commonly used because they are simple to apply and technically defendable. The collaboration between the two organizations aims for improvements in thermal-hydraulic models. Improvements to methods and models can reduce modeling uncertainties and excessive conservatism, allowing greater reactor performance without reducing predicted safety margins. Therefore, these improvements potentially expand the design space, making additional designs feasible, which, in turn, may achieve better fuel utilization and proliferation resistance in the designs. This report provides justification for the lateral node size (stripe width) chosen for the 1-D computer models during the collaboration. The lateral heat conduction in OPAL fuel plates is evaluated using the two-dimensional conduction capability of PLTEMP/ANL to determine a representative hot-stripe width that can be conservatively used for one-dimensional steady-state thermal-hydraulic analysis. Basically, a representative hot-stripe width needs to be selected so that the actual heat flux averaged over this stripe width is greater than or equal to the maximum heat flux obtained by the two-dimensional analysis. Detailed results in this report show that a hot-stripe width of 8.125 mm can be conservatively used for the one-dimensional thermal-hydraulics analysis of all plates. Therefore, eight equal stripes (8.125 mm each) are recommended to be modeled in the MCNP and PLTEMP/ANL full core models of the OPAL reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multiphysics Simulations of MSRE with NEAMS Thermal Hydraulics Tools

This report documents the benchmarks being developed and simulations performed using tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing MSRE experimental data. In FY23, three main work scopes were investigated under the NEAMS MSR work package at ANL. The first scope investigated the Griffin-SAM coupling model for simulating the pump startup transient experiment of MSRE. The analyses start with a simple model (single-channel, single-lattice), gradually adding more details (multi-channel, full-core) into the model. The results show that the reactivity loss curve is very sensitive to the axial boundary conditions and the radial core discretization. The simple model can predict a similar reactivity trend as that of the more sophisticated model, which is likely due to error cancellation. Accurately modeling the axial boundary condition may further improve the reactivity trend but would require significant efforts to generate the mesh of the MSRE inlet and upper plenum. The core channel radial discretization for the Griffin-SAM coupled model also depends on the flow distribution. Given the complex geometry in the inlet plenum, the flow distribution needed to be calculated from CFD analysis, which was performed using the NekRS code. This analysis employed a MSRE CAD model developed by Copenhagen Atomics. The CAD model was disassembled to keep the inlet plenum region only, which was subsequently cleaned and modified so that the mesh generated is under the memory limit. The results are merged to a few radial regions to show that the flow rate is highest in the central region. This would be useful for future improvement of the Griffin-SAM coupling model of the MSRE core. The last task investigated is tritium transport modeling using the standalone SAM code. This task aimed to initiate the effort to demonstrate and validate the tritium transport model implemented in SAM. The preliminary investigation employed an MSRE model consisting of the primary loop. Three tritium transport pathways were examined including the retention in the graphite, the permeation through the HX tube wall, and the removal from the off-gas system. The results compare well with the MSRE data, but improvements are still needed on the initial conditions (i.e., the present state may not have reached equilibrium), the boundary conditions, the off-gas system modeling, and a better numerical strategy to reach the equilibrium state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

STAT7 v1.2 User Guide: The STAT7 Code for Statistical Propagation of Uncertainties in Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for thermal hydraulics analyses of research and test reactors. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to the fit functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one latera stripe of a fuel element. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗