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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

High-Burnup LOCA Burst Susceptibility BISON Analysis in PWRs and BWRs

Accurately assessing high-burnup fuel behavior during loss-of-coolant accidents (LOCAs) is essential for understanding fuel fragmentation, relocation, and dispersal (FFRD) risks across the US light-water reactor fleet. This work updates previous Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program multiphysics LOCA analyses for a pressurized water reactor (PWR) and a boiling water reactor (BWR) by incorporating recent model and material property advancements in the BISON fuel performance code, including a high-burnup structure (HBS) model, revised cladding burst criteria, and updated thermal–mechanical correlations. This update was needed to support ongoing industry initiatives and upcoming regulatory changes. Full-core, rod-resolved operating histories generated using Virtual Environment for Reactor Analysis (VERA) and system-level LOCA conditions obtained from TRACE were applied to statistically representative rod samples in BISON to evaluate burst behavior and FFRD susceptibility. These calculations used two cladding burst correlations and three fuel pulverization models so that the predictions of these models could be compared. The updated PWR simulations show markedly improved numerical stability as the number of crashed simulations decreased by 95% compared to the previous study, and hence higher confidence in results. The updated PWR simulations predicted cladding bursts exclusively among once-burned, high-power rods, with two different cladding burst models identifying the same burst-susceptible population. Resulting FFRD susceptibility estimates are significantly reduced compared with earlier studies, driven by cooler predicted fuel and plenum temperatures, lower hoop strains, and reduced fission gas release in the updated models. In contrast, none of the BWR rods were predicted to burst under either burst criterion, reaffirming minimal BWR FFRD susceptibility even with updated HBS and material models. Comparisons between the PWR and BWR end-of-cycle predictions are made. Comparison with prior work highlights significant shifts in PWR fuel performance metrics and confirmation of earlier BWR conclusions. Overall, the updated results underscore the importance of having high-resolution detailed modeling capability and continuously integrating evolving material models and physics into high-resolution multiphysics simulations. The unified assessment presented here strengthens confidence in predicting high-burnup LOCA behavior by improving agreement between different cladding burst correlations. These results also provide an improved foundation for future BISON model development, FFRD susceptibility calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Update on Radiochemical Assessment of High Burnup Commercially Irradiated Fuel

This work documents an effort to collect burnup measurements on a high burnup rod, designated 6XV, and first cycle accident tolerant fuel (ATF) rod, designated 47I, to enable benchmarking of fuel performance codes and neutronics codes. In addition to measurements, Virtual Environment for Reactor Applications (VERA) full-core-depletion analysis was also performed for the rods that were experimentally analyzed to provide an opportunity for code validation. This effort focuses on collecting data from rods irradiated at Byron Generating Station and shipped to the Oak Ridge National Laboratory (ORNL) hot-cells. This data will also anchor non-destructive examination evaluations of burnup of the various fuel rods undergoing postirradiation examination (PIE) at ORNL. Previous PIE of these fuel rods provides some guidance on the burnup trend across the fuel. Axial gamma spectroscopy scans provide a measure of relative changes in burnup across a fuel pin. Mass spectrometry based burnup measurements performed for this work at specific axial locations in the fuel are fully quantitative. By combining the mass spectrometry data with the gamma scans it is possible to more quantitatively evaluate axial variations in burnup across the entire fuel pin [1]. The combined set of burnup evaluations will be made available to other organizations that have an interest in high burnup radiochemistry data for validation of neutronic simulations and source term evaluation such as the Nuclear Regulatory Commission (NRC).

Harp, Jason [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2021

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗
DOE OSTI2026

Nuclear Data Impact on Key Metrics for a Representative Molten Chloride Fast Reactor Model

Nuclear data are an essential component of the foundation on which all modeling and simulation methods and tools are relying upon, from the front end to the back end of the nuclear fuel cycle. In this study, the impact of uncertainties in nuclear data is investigated for a representative molten chloride fast reactor, for several important metrics, including eigenvalue, reactivity differences, and nuclide inventories in fuel at 5-yr irradiation. Uncertainty of keff for a full core model was found to be similar between the fresh fuel and the irradiated fuel states (1.7-1.8%), with its primary driver being the uncertainty in the 235U (n,γ) cross section. The results obtained for the reactivity differences show large uncertainties, of over 100%, in elastic scattering sensitivities of several nuclides, which led to large uncertainties of temperature reactivity differences for cladding and reflector. These results provide evidence that the currently applied methods may not be sufficiently adequate for ensuring the reliable determination of such metrics.

Procop, Germina [ORNL] (ORCID:0000000342226393)↗
NASA NTRS1994

Probing the galactic disk and halo. 1: The NGC 3783 sight line

We report a study of Galactic disk and halo absorption toward the Seyfert galaxy NGC 3783 which has Galactic coordinates l = 287.46 and b = +22.95. The data were obtained with the Goddard High Resolution Spectrograph operating at medium resolution with the Large Science Aperture, which produces a line spread function having a sharp core (full width at half maximum (FWHM) approximately 20 km/s) and broad wings extending to +/- 70 km/s. Ion species detected in absorption near zero LSR velocity include C IV and N V for high ions, and C I, Mg II, Si II, and S II for low ions. Absorption from a high-velocity cloud (HVC) at a velocity of +240 km/s along the sight line is also detected in the ion species of S II, Si II, and probably C I. This is the first reported case where S II and C I absorption has been detected in a HVC. The S II lines are especially useful since metal abundance estimates based on S are largely unaffected by dust grains. The study is aided by the availability of 21 cm emission data.

Lu, Limin↗
DOE OSTI2020

As-run physics analysis for the EPRI-2 experiment for cycle 153B

The purpose of this ECAR is to document the as-run physics analysis and source term in Curies and decay heat rate (Watts) for the shipment of the EPRI-2 experiment located in the Center Flux Trap (CFT). EPRI-2 was irradiated during cycle 153B. The experiment irradiation ended on April 13, 2013. The as-run heat rates and flux/fluence results were calculated using the MCNP ATR full core model. The heat rates and flux are calculated based on an average center lobe power of 30.8 MW. The fluence is based on the operating time of 13.45 days. ORIGEN2 Version 2.2 is used to determine the decay heat rate, gamma spectrum, and source term for the EPRI-2 experiment. The source term analysis was performed as requested by the project staff to demonstrate compliance with shipping requirements in the GE-2000 cask following irradiation as well as provide an estimated source term to support Post-Irradiation Examination as needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

As-Run Physics Analysis for the EPRI-1 Experiment

The purpose of this ECAR is to document the as-run physics analysis for the Electric Power Research Institute (EPRI)-1 experiment located in the Center Flux Trap (CFT). EPRI-1 was irradiated during cycle 157C. The experiment irradiation ended on February 15, 2015. The as-run heat rates, DPA, and flux/fluence results were calculated using the Monte Carlo N-Particle (MCNP) Advanced Test Reactor (ATR) full core model. The heat rates, DPA, and flux are calculated based on an average center lobe power of 21.1 MW. The fluence and DPA is based on the operating time of 5.4 days. ORIGEN2 Version 2.2 was used to calculate the isotopic concentration in grams for the EPRI-1 experiment. The source term analysis was performed as requested by the project staff. The results of this analysis may be used to support Post-Irradiation Examination (PIE)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

JMOCUP Physics Depletion Calculations for the As-Run AGR-5/6/7 TRISO Particle Experiment in ATR Northeast Flux Trap

This ECAR documents the as-run Jim Sterbentz’s MCNP ORIGEN Coupled Utility Program (JMOCUP) physics depletion calculation and the calculated results for the AGR-5/6/7 irradiation experiment in the Advanced Test Reactor (ATR). AGR 5/6/7 was irradiated for nine power cycles in the northeast flux trap. The depletion calculations were performed to provide input data for a variety of other engineering analyses supporting the AGR-5/6/7 experiment along with post-irradiation characterization of the tri-structural isotropic (TRISO) particle fuel compacts. Detailed full-core MCNP models and Oak Ridge Isotope Generation (ORIGEN2) radionuclide generation models were specifically developed as part of the JMOCUP Monte Carlo depletion calculations. The MCNP ORIGEN2 computer codes were coupled using the well-established JMOCUP utility modules to couple the two codes and run the depletion calculations. The physics calculations were performed in support of the Advanced Gas Reactor (AGR) program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling and Analysis of a Micro-reactor using the DireWolf Code Suite

Micro-reactors, specifically those cooled with heat pipes, are tightly coupled systems which require multi-physics tools to accurately model both steady-state and transient behavior. Idaho National Laboratory’s DireWolf code suite is tailor-built to model heat-pipe reactors. DireWolf’s ability to model the coupled steady-state behavior of a heat-pipe micro-reactor are demonstrated herein. The DireWolf modeling methodology is described, and a micro-reactor design is introduced. A full-core model is developed, and the steady-state coupled neutronics-thermal analysis is performed. Results of the analyses, including flux, power, and temperature distributions are presented and discussed. This is a Westinghouse Electric Company (WEC) led publication.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

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 [1, 2] 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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2022

Native Natural circulation in VERA

This work seeks to introduce a native natural circulation solver to VERA's thermohydraulics capabilities. This work begins with initial investigations into the viability of such a solver by determining the fundamental equations as well as necessary inputs. The method was then implemented in the form of a Python script coupled to CTF which was compared to TRACE. Once viability of the natural circulation solver was established, it was then implemented natively in MPACT for use with previously created simplified thermohydraulics solvers. This new implementation of the native natural circulation solver in MPACT is then tested for single physics problems for both a BWR and PWR single channel model and compared to TRACE where it shows an error of about 30% for BWRs and about 15% for PWRs. The solver is additionally tested for NuScale based models where it demonstrates an error to the reported flow rate of about 10%, resulting in an eigenvalue error of only 138 pcm when comparing fixed flow rate multiphysics runs. Finally, the solver is tested in multiphysics systems for a small PWR assembly, a small BWR assembly, and a full core approximate model of the NuScale design. In these multiphysics runs, it demonstrates no impact on multiphysics convergence stability, as well as negligible impact on runtimes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Development of a high fidelity model of the CROCUS experimental reactor

Measurements of scalar flux distributions with fine spatial and energy resolutions are needed to remedy one of the validation shortcomings of the novel neutronics full core solvers, such as MPACT and nTRACER. While a very detailed resolution of the flux can be calculated with such codes, only a limited experimental data set is available to check their accuracy. Such type of measurements are on-going at the zero power reactor CROCUS, operated at the Laboratory for Reactor Physics and System Behaviours of the EPFL, thanks to the development of advanced miniature neutron detection systems. This kind of experimental data would provide the community with a suitable benchmark for the validation of high fidelity neutronics solvers. In parallel, a multi-physics solver for steady-state and transient analysis of nuclear reactors, named GeN-Foam, has been developed. Based on the finite-volume OpenFOAM library, GeN-Foam provides us with enough flexibility to analyze non-conventional reactor geometries such as that of CROCUS. While CROCUS heterogeneities cannot be modeled by MPACT and nTRACER for the moment, GeN-Foam offers a unique opportunity to build a high fidelity model which mimics these codes' method to reach sub-pin simulation resolution. This document aims at describing the work achieved to get from the existing GeN-Foam model of the CROCUS reactor based on a structured mesh and using the neutron diffusion, the first high-fidelity model using discrete ordinates method as an approximation to neutron transport and an unstructured mesh for inter lattice water gap description and sub-pin heterogeneous modeling. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Scalable algorithms for domain decomposition of Monte Carlo neutral particle transport simulations on unstructured mesh

Monte Carlo simulations are the gold standard in particle transport simulations, for shielding and nuclear reactor simulations Full core simulations can require terabytes of memory for tracking pin-resolved (up to 50k pins), nuclide (200 tracked), axial (20-200 zones) reaction rates Domain decomposition, splitting the calculation and the memory storage, among several nodes is challenging as it involves communication of particles across domain boundaries We created a new MOOSE application, MaCaw, to study this problem Algorithms developed for ray tracing in MOOSE were adapted for simulating neutral particle transport Physics and tally routines in OpenMC, dynamically linked, are called from the application

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Uncertainty Quantification and Sensitivity Analysis of Non-Nuclear Advanced Controls Testbed Reactor Mockup

The research presented in this report describes our progress in applying stochastic methods and uncertainty quantification, parametric study, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated with Griffin, with the goal of developing a reduced order (surrogate) model which can be rapidly sampled while perturbing multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums, which are controlled by a hybrid proportional, integral and derivative (PID) controller, actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time and $\beta$ from a transient Griffin simulation executed with the improved quasi-static method to provide the kinetic parameters) as inputs to functions which control the CD rotation angle. Using a number of stochastic method approaches, we developed a dual purpose training-surrogate model of the NTP system using polynomial regression. The trained model can be rapidly sampled while simultaneously perturbing various input parameters of the model, such as coefficients on the PID control, or temperature (directly affect the neutron cross section). The surrogate model delivers accurate results orders-of-magnitude faster (minutes, not days) than the base model. Once the base model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs and their effect on the output. For example, coefficients used in the PID control system may vary due to some physical interference, or there may be uncertainty in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters and operational boundaries can be determined. The goal of this work is to support development of an advanced control system to operate CDs in a functioning NTP system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Validation of Pronghorn’s Subchannel code using the EBR-II shutdown heat removal tests

A Subchannel application (Pronghorn-SC) is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This work utilizes the EBR-II SHRT tests to validate the subchannel capabilities. Multi-scale and multi-physics coupling is used to improve the fidelity of the subchannel code calculations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Modeling of the Molten Salt Reactor Experiment with SCALE

A SCALE model was developed for the Molten Salt Reactor Experiment (MSRE) benchmark that was recently added to the International Handbook of Evaluated Reactor Physics Benchmark Experiments. This SCALE model served as a basis for criticality calculations and nuclear data sensitivity and uncertainty analyses with the Monte Carlo code Shift and the TSUNAMI computational capabilities in the SCALE code system. The focus of this work is the assessment of the impact of nuclear data on the calculated eigenvalue results in support of the discussion of differences between the calculated and the experimental eigenvalue result. The differences in the eigenvalues obtained using the ENDF/B-VII.0, ENDF/B-VII.1, and ENDF/B-VIII.0 nuclear data libraries cover a relatively small range of ~230 pcm. Since eigenvalue sensitivity of the MSRE is dominated by the neutron multiplicity and neutron capture of 235 U and elastic scattering in graphite, relevant changes in the ENDF/B libraries for nuclear reactions (such as carbon capture) that caused large differences in other graphite-moderated systems did not have a significant impact. Propagation of nuclear data uncertainty results in an eigenvalue uncertainty of ~700 pcm with the major contributors being 235 U neutron multiplicity, graphite elastic scattering, and 7Li neutron capture. All calculations resulted in large differences of ~2000 pcm in eigenvalue compared to the benchmark experimental value. Several potential contributors to this difference—including uncertainties and gaps in the knowledge of the material, geometry, and nuclear data—were identified. Simplified models of the full MSRE core were developed, and similarity assessments were conduced with the full MSRE core model. It was found that simplified models can serve as adequate surrogates of the full-core model such that they can be used for performing selected nuclear data performance assessments with a lower computational burden.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗