Engineering Papers⌕ Search

Engineering topics

Royston, Katherine

Publications and source records attributed to Royston, Katherine.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Photonuclear Physics in SCALE

Radiation shielding modeling involves a wide range of applications, such as power reactors, critical experiments, portable sources, linear accelerators, and fusion systems. Particle energy ranges vary from thermal energies to fast energy spectra. MAVRIC, as the radiation shielding sequence of the SCALE code system [1] developed and maintained by Oak Ridge National Laboratory, is capable of modeling both neutron and photon interactions in matter and is widely used by the US Nuclear Regulatory Commission, Department of Energy, the National Nuclear Security Administration, and radiation shielding practitioners all around the world. Although most applications don’t involve highly energetic photons that can cause photonuclear reactions and create coupled chains of high energy, ionized electrons and bremsstrahlung radiation within a material medium, when present, these reactions can have a significant impact upon the downstream particle fluxes from incident photons above 20 MeV. For example, linear accelerators used for medical isotope production are examples of applications in which high-energetic photons are generated and used for driving or controlling other components of a complex facility. Radiation shielding analyses of such a facility would require accounting for any and all significant sources of radiation. Unfortunately, MAVRIC included no data and methods to model either photonuclear reactions or the bremsstrahlung photons until this work, which enabled accounting for both. Photonuclear reaction cross sections were imported from the Evaluated Nuclear Data File libraries, and corresponding methods were implemented in the SCALE code system; see Section 2 for how the nuclear data are imported and used. Bremsstrahlung photons were modeled by assuming thick-target approximation whereby all generated electrons and positrons due to photon interactions are assumed to be absorbed within the same material in which they were born, yielding a quick particle transport simulation method by using fast lookup tables for generating subsequent gammas; see Section 4 for details. A suite of verification models was used to test implemented methods and associated data. Tests yielded satisfactory results for the new features. Enhancements will be released as part of the next major SCALE release (v7.0), which is anticipated to be deployed in 2025, with beta releases available in the interim.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2021

Validation of Light Water Reactor Ex-Core Calculations with VERA

The Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) is a reactor simulation software. It offers unique capabilities by combining high-fidelity in-core radiation transport with temperature feedback by using MPACT (a deterministic neutron transport code) and COBRA-TF (a thermal-hydraulic code) with follow-on, fixed-source transport calculations using the Shift Monte Carlo code to calculate ex-core quantities of interest. In these coupled calculations, MPACT provides Shift with the fission source for follow-on ex-core calculations. These ex-core simulations can be set up to calculate detector responses, as well as the flux and fluence in ex-core regions of interest, such as the reactor pressure vessel, nozzle, and irradiated capsules. A Watts Bar Nuclear Plant Unit 1 (WBN1) ex-core model was developed, as described in this paper, and this model was used to perform coupon calculations. The results for the coupon flux calculations show close agreement with the reference values for cycle 1 produced by the two-dimensional Discrete Ordinates Transport (DORT) code and presented in a BWXT Services Inc. report. However, differences in the results (10%) seen in cycles 2 and 3 and the reasons for these differences are discussed in this paper. The VERA WBN1 model was also used to perform a vessel fluence calculation for cycle 1. Additionally, a collaboration between CASL and Duke Energy led to the first code-to-code validation of VERA for reactor ex-core applications that used a model for the Shearon Harris reactor. Results from this collaboration show excellent agreement between VERA and the Monte Carlo N-Particle Transport Code for the detector response calculations. The work performed under this collaboration is also detailed in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Secondary-Source Core Reload Modeling with VERA

The CASL reactor simulation package VERA has been adapted to provide high-fidelity simulation capabilities for modeling source range detector response during subcritical reactor configurations. New features include the activation and shuffling of secondary-source assemblies, use of burned fuel neutron emission data from the ORIGEN depletion solver to the MPACT deterministic neutron transport solver, allowance of user-defined sources in MPACT based on material composition, ability to solve the subcritical source-driven system with neutron multiplication using the MPACT diffusion solver, and transfer of the calculated fission source from MPACT to the continuous-energy Monte Carlo solver Shift for final detector response evaluation using the CADIS methodology for variance reduction. These new capabilities were benchmarked against Watts Bar Unit 1 plant operating data for the first few fuel loading steps and were found to demonstrate excellent agreement with the measured data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Impact of Radial Reflector Fidelity on Neutronics and Vessel Fluence Simulations

The Consortium for Advanced Simulation of Light Water Reactors is developing the Virtual Environment for Reactor Applications (VERA), and the MPACT code, which is the primary deterministic neutron transport solver in VERA, provides sub-pin level flux and power distributions as part of full-scale cycle depletion and analysis. In such calculations, an important aspect is the radial reflector treatment. To improve the fidelity of the radial reflector treatment, MPACT was extended to approximate the modeling of the reactor’s structural components such as the core shroud, barrel, neutron pads, and vessel. This work explores several modeling configurations with varying levels of fidelity and computational burden and assesses the importance of modeling fidelity on the eigenvalue and pin power distribution. Two two-dimensional (2-D) problems were analyzed to assess the impact on eigenvalue and pin power distributions with low-fidelity, coarse square cell reflector representations: (1) a Watts Bar Nuclear Plant Unit 1 (WBN1) quarter-core slice with depletion and (2) an AP1000 quarter-core slice. In this work, the analyses showed that the effect on eigenvalue is fairly small, but the effect on pin power is more pronounced, especially locally in the assemblies closest to the periphery, where the maximum pin power difference is nearly 3.5% in the AP1000 case. Two additional 2-D problems were used to assess the comparison between the low-fidelity coarse square cell treatment and a high-fidelity geometric representation that uses subpin material specification: (1) the same WBN1 quarter-core slice and (2) a representative model of the NuScale small modular reactor (SMR), which features a solid reflector design with moderator holes. These results demonstrate that even a coarse, low-fidelity representation adequately captures the necessary simulation characteristics. Last, these capabilities were applied to the 2-D WBN1 quarter-core depletion to assess the impact on vessel fluence using VeraShift. From adjoint calculations, pins along the periphery were observed to be of highest importance for fluence calculation, so the impact of the reflector representation in MPACT could theoretically substantially affect the predicted result. However, it was observed that the change in pin powers along the periphery minimally impacts the maximum vessel fluence with a difference within the statistical uncertainty but provides terrific insight on the sensitivity of the peripheral pins.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

NEAMS Milestone Report: M2MS–20OR030102 FW–CADIS PWR Ex-Core Analysis with Shift Through VERA

This report presents the work completed for the NEAMS milestone M2MS–20OR030102 titled "FW–CADIS PWR Ex-Core Analysis with Shift through VERA." The work completed for this milestone includes the implementation, integration, and optimization of memory and performance improvement methods in Shift for fully coupled ex-core calculations through Virtual Environment for Reactor Applications (VERA). Fully coupled in this context means the transfer of moderator boron concentration, pin-wise fission source, depleted compositions, temperatures, and moderator densities from MPACT (with COBRA-TF (CTF)) to Shift. The ability to run ex-core calculations with VERA has been enabled and used for several years by Consortium for Advanced Simulation of Light Water Reactors (CASL) partners. However, this implementation was limited and potentially computationally burdensome. This work has enabled the ability to run higher-fidelity ex-core calculations on moderate computing clusters by focusing on multithreading, domain decomposition, and Forward-Weighted CADIS (FW-CADIS) variance reduction. Tests performed on small cores, a small modular reactor (SMR), and CASL progression problems show very promising memory reduction and computational performance. Recommendations for settings when running fully coupled high-fidelity ex-core calculations with VERA are documented. Without these optimization methods, many processors on a compute node would be left unused for the entire ex-core calculation. Therefore, these methods enable the user to better use the resources available and reduce computation time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Subcritical Multiplication Ex-Core Calculations with Shift

This report details the development and testing of Virtual Environment for Reactor Applications (VERA) as applied to full core subcritical multiplication ex-core calculations with Shift. The limiting factor in being able to run these detector calculations with Monte Carlo is the amount of computational memory used by a fully-loaded full core problem. The main development features added to VERA in this effort were enabling source biasing and separable weight windows and integrating multithreading in Shift. With these capabilities, a subcritical multiplication detector response calculation on a full core Watts Bar Nuclear Plant Unit 1 (WBN1) ex-core model was successfully run on a moderately sized computing cluster.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Demonstration of Comprehensive Ex-Vessel Fluence Capability

Several recent developments in the Consortium for Advanced Simulation of Light Water Reactors’ (CASL) Virtual Environment for Reactor Applications (VERA) were tested to highlight its capability to perform ex-vessel fluence and reaction rate calculations. VERA offers unique capabilities for integrating the deterministic neutronics code, MPACT, with Shift, a Monte Carlo code, to perform high-fidelity in-core and ex-core radiation transport. Applications such as pressure vessel fluence, ex-core detector response, and coupon irradiation analyses take advantage of this coupling. For these applications, MPACT performs the in-core radiation transport with temperature feedback and isotopic depletion through the direct coupling with the CTF subchannel thermal-hydraulics code. Then, MPACT provides the fission source to Shift for a follow-on fixed source radiation transport calculation that tallies all the ex-core responses of interest for each time-dependent statepoint. The variance reduction method, CADIS, which is implemented in Shift, allows for efficient performance of ex-core transport for calculation of ex-core quantities of interest. This milestone report identifies limitations to the current methodology in place and showcases VERA’s current capability to model ex-core quantities of interest.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Source Range Detector Response Modeling using VERA

The CASL reactor simulation package VERA has been adapted to provide high-fidelity simulation capabilities for modeling source range detector response during subcritical reactor configurations. New features include the activation and shuffling of secondary source assemblies, transfer of burned fuel neutron emission data from the ORIGEN depletion solver to the MPACT deterministic neutron transport solver, allowance of user-defined sources in MPACT based on material composition, ability to solve the subcritical source-driven system with neutron multiplication using the MPACT diffusion solver, and transfer of the calculated fission source from MPACT to the continuous energy Monte Carlo solver Shift for final detector response evaluation. These new capabilities were benchmarked against Watts Bar Unit 1 plant operating data and found to be in very good agreement.

Gentry, Cole↗
DOE OSTI2019

Assessment of Activation on Level L3 of the Tokamak Building due to the ITER Tokamak Cooling Water System

The ITER fusion reactor is being built to demonstrate the feasibility of fusion power and will be the largest tokamak in the world. The tokamak cooling water system (TCWS) will extract the heat generated during operations and includes large amounts of piping and equipment such as pumps and heat exchangers (HXs) that are located in a large shielded region on level L3 of the tokamak building. During operation, water in the TCWS will be activated by plasma neutrons and then flow into this shielded region. The activated coolant will in turn activate the steel in the TCWS during operation and result in an activation gamma source and radiation responses that must be assessed to inform equipment selection and maintenance schedules.The activation of materials in the shielded region of level L3 was assessed at several decay times and for different equipment options using the Oak Ridge National Laboratory (ORNL) shutdown dose rate (SDDR) code suite. The ORNL SDDR code suite implements the rigorous two-step method using the Multi-Step Consistent Adjoint-Driven Importance Sampling (MS-CADIS) method to create effective neutron variance reduction parameters for the photon response of interest. In conclusion, two different HX designs, shell and tube and shell and plate, were considered, as well as the impact of cobalt impurities in steel equipment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Increased Plutonium-238 Production via High Flux Isotope Reactor Permanent Beryllium Reflector Redesign

Irradiation of 237Np-bearing targets in the permanent beryllium reflector (PBR) of Oak Ridge National Laboratory’s (ORNL) High Flux Isotope Reactor (HFIR) results in high-purity 238Pu that can be used as a reliable power source for deep space and planetary National Aeronautics and Space Administration missions. However, HFIR’s 238Pu production capability is constrained by the available irradiation volume in its PBR. In preparation for the HFIR beryllium change-out in 2023, ORNL staff have redesigned the PBR to include six additional irradiation sites, be more versatile with respect to irradiation and scattering experiments, and enhance its thermal-structural performance. The new PBR design offers large potential increases in annual 238Pu production at ORNL.

Chandler, David↗
DOE OSTI2018

VERA Infrastructure

The Virtual Environment for Reactor Applications (VERA) Infrastructure provides the infrastructure to build, configure, test, install, and package for release the VERAS suite of codes that is being developed by the Consortium for the Advanced Simulation of Light-water reactors (CASL). It includes the list of codes that can be built within VERA, the location on the configuration-controlled CASL development computer, the version number used in any particular software installation, the location of the third-party libraries used for testing VERA on a suite of the CASL testing and development computers, and the scripts to configure, build, test, install, and package for release a particular version of VERA. The VERA Infrastructure contains a list of scripts, configuration options, and utilities that are specific to CASL, but contain no physics or numeric, nuclear technology, or encryption technology.

Baird, Mark↗