Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear systems code”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Radiation Damage Analysis of FNSF Components Using McCad and MCNP

The Fusion Energy System Studies Fusion Nuclear Science Facility (FESS-FNSF) concept represents a transitional step between ITER and a commercial fusion power plant. The FNSF is a conceptualized D-T fueled tokamak with 518 MW of fusion power that has been extensively used to explore and optimize design features. The energetic 14.1-MeV neutrons can produce significant localized heating and activations, and can cause damage to plasma-facing components, which can determine maintenance/outage scheduling needs and also impact the lifetime of the device as a whole. This study illustrates a neutronics analysis that was conducted on a 22.5-degree symmetric sector of the FNSF with the goal of understanding the neutron heating and radiation damage that can be characterized by quantifying the displacements per atom (dpa). Concurrently, this study also focused on the development of analysis capabilities by converting a three-dimensional computer-aided design model of the FNSF into MCNP6.2 input using the McCad code. Accordingly, some confirmatory results on tritium production and the tritium breeding ratio (TBR) are provided to support model validation. The results produced by MCNP6.2 simulations showed that the highest heating and damage occurred in the outboard region, which concentrated approximately 290 MW of the total nuclear heating, in contrast to 97 MW within the inboard region. These results are consistent with previous studies that employed earlier versions of the FNSF concept and different modeling approaches. Furthermore, this study also provides additional details on neutron wall loading, as well as total heating from neutrons and gammas, results which show the total heating of the device (16 sectors) is approximately 477.83 ± 0.80% MW, indicating a neutron energy multiplication factor of 1.15. Additionally, the capability to calculate hydrogen and helium production, as well as dpa, is illustrated. Finally, the neutronics effects of using alternative materials to tungsten carbide were evaluated for the vacuum vessel, low-temperature shield, and structural ring components, which showed that compounds like YH 2 , Mg(BH 4 ) 2 , and ZrH 2 could reduce the total heating on the magnet and also reduce the TBR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transient analysis of a micro-reactor using the DireWolf code suite

Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory's DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf's ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete rotation of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

CTF Theory Manual: Version 4.3

Coolant-Boiling in Rod Arrays—Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and three-dimensional Cartesian forms of its governing equations are available for the solution. The code was originally developed by Pacific Northwest Laboratory in 1980, and had been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently rebranded as CTF. CTF was later adopted in the early 2010s by Oak Ridge National Laboratory (ORNL) for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, which led to a significant advancement of the code software quality, modeling accuracy, testing systems, and capabilities for improved support of modeling of common LWR nominal and transient behavior. As part of the improvement process, it was necessary to generate sufficient documentation for the public domain code which had lacked such material upon being adopted by RDFMG. This document serves as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code, as well as the numerical solution process utilized. Additional documents available in the CTF documentation suite include the user manual and verification and validation manual.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Fermi-Hubbard Transport Solver v. 1.0

SAND2021-1511 O. The code takes a description of a system comprised of two leads that are sources and sinks of electrons, a set of abstract sites that those electrons can occupy, and their connectivity. It then predicts the current passing between the two leads through the set of abstract sites. There are numerical parameters that define the rate at which electrons tunnel between sites and another parameter that determines the energy of having two electrons on the same site. When the ratio of these energies and biases are set just right, they can be used to predict the existence of Coulomb blockade in which no current flows because the energetic penalty would be too high. That is, the code predicts a current and the absence of a current which is predictive of Coulomb blockade. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Brickson, MitchellI↗
DOE OSTI2020

BISON Capability, Validation and Demonstration for Reactivity-Initiated Accidents

Reactivity Initiated Accidents (RIA) are design basis accidents that could have adverse impacts on the core coolability. In the unlikely event that sufficient reactivity is inserted into the reactor core by the ejected/dropped control rod, prompt energy deposition into the fuel can occur, which when sufficiently high can lead to fuel rod failure or, at large energy deposition levels, expulsion of UO2 fragments or molten UO2 material from the fuel rod. This results in a release of fission product and fuel into the coolant and has the potential to compromise core coolability and threatening the pressure boundary of the primary coolant system. The design basis RIA is one of the industry challenging problems that the CASL aimed to address; the CASL RIA Challenge Problem Charter [1] states, ?The Pressurized Water Reactor (PWR) Rod Ejection Accident (REA) and Boiling Water Reactor (BWR) Control Rod Drop Accident (CRDA) are postulated accidents with consequences that are important to nuclear safety (fuel rod integrity and core coolability). Currently each reload core design must be analyzed to meet regulatory acceptance criteria. The goal of CASLs ModSim capability for RIA is to model the event at a higher fidelity, with validation to existing tests, to better model the transient neutronics and the progression of the fuel and cladding thermal-mechanical behavior. These improved analytical capabilities can be used to better inform reload core design, limits on fuel assembly discharge burnup, restrictions on placement of fuel in the reactor, control rod insertion limits, operating margin, and performance sensitivities.? In support of that charter, BISON, the fuel performance code has been used to demonstrate the simulation of thermal-mechanical behavior of LWR fuels during a reactivity-initiated accident. The combination of mechanical, thermal, and thermal-hydraulic phenomena all present during an RIA makes a multi-physics code such as BISON a valuable tool for modeling these scenarios. This paper will highlight many of the activities associated with BISON relevant to reactivity-initiated accident capability development and validation efforts. These efforts have been performed under both CASL and the Nuclear Energy Advanced Modeling and Simulation (NEAMS) programs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Validation of Thermohydraulic Simulations using RELAP for Critical Dual Purpose Canisters

Simulations using the RELAP5-3D systems analysis code and STAR-CCM+, a high-fidelity CFD code, were performed on a model of a real-scale MPC-32 dual purpose canister undergoing a criticality event. These simulations were performed to provide code-to-code verification of the natural convection capabilities in RELAP5. The STAR-CCM+ simulations revealed various natural convection flow patterns that aided in the cooling of fission heat produced in the fuel rods in the canister. Negligible spatial temperature variations were observed across the canister in the STAR-CCM+ results, supporting the use of RELAP5 in predicting the sub-cooled regime. The surface-averaged rod temperature predicted by RELAP5 compared well with STAR-CCM+, verifying the predictive capability of RELAP5 for sub-cooled conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

MOSCATO Solver Development and Integration Plan

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. Although capable of large-scale simulations, MOSCATO is not well-suited to parametric studies of complete reactor geometries. These types of simulations are instead better handled by reduced-order modeling codes such as ORNL’s Mole code. Reduced-order simulation tools like Mole, however, are dependent on high fidelity correlations to account for complex, coupled three-dimensional phenomena that they do not directly simulate. Tools such as MOSCATO must therefore be used to create these correlations, as suitable empirical relationships are not available for most molten salt systems. Toward that end, we used the Nek-derived version of MOSCATO to create new mass transfer correlations for three relevant cases including tubular, tube bank, and subchannel geometries. These new correlations are more accurate than any existing ones and can be readily integrated into any reduced-order modeling tools that are targeting full-scale MSR simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of fast-running eighth core nodal solution

A fast-running, eighth-core solution technique has been developed for the quick core calculations in ANC, the Westinghouse core analysis code based on the nodal expansion method, and the BEACON{sup TM} Core Monitoring System (BEACON). This technique, referred to as 'Fast Eighth' (Fast8), enables users to obtain satisfactory macroscopic core calculation results such as reactivity, axial offset and assembly power much more quickly than calculations performed in the standard full-core or quarter-core geometry. Fast8 uses three techniques: folding the radial geometry into eighth-core geometry, reducing the axial geometry into a coarse axial mesh geometry, and smearing the spacer grids into nodes on-the-fly after the first calculation in the standard geometry. These geometry adjustments in Fast8 contribute to the reduction of the number of neutronic and material nodes while maintaining the reaction rates in each node. The axial mesh condition dependency in Fast8 is investigated in a xenon transient calculation in single assembly geometry, and the Fast8 maximum axial mesh size is determined based on these results. Fast8 is then applied to the load swing calculation, and the results are compared with the reference, full-core geometry, solution. The critical boron concentration and axial offset of Fast8 show good agreement with the reference solutions, and that Fast8 is a good technique for obtaining macroscopic core calculation results in significantly shorter time. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Demonstration of MELCOR and MACCS Capabilities for Molten Salt Reactor Decay Heat Removal During both Normal Operations and Salt Spill Scenarios

This report provides a demonstration of MELCOR and MELCOR Accident Consequence Code System (MACCS) capabilities to perform a dose assessment for a Molten Salt Reactor (MSR) off-gas system. A primary containment system salt spill is used as the off-normal scenario, along with a normal operation dose assessment for comparison. This report discusses the tools, methods, and information used in this assessment so that it may be utilized as a starting point for future advanced reactor consequence analyses. This report also highlights several gaps, to include the need for reactor inventory information specific to advanced reactors, and the need for specific atmospheric transport models that take into account the unique deposition behaviors of tritium and carbon-14, and makes recommendations for closing these gaps. This report satisfies the DOE NE Milestone M4RD-21SN0601062.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Relativistic Configuration-Interaction and Perturbation Theory Calculations for Heavy Atoms

Heavy atoms present challenges to atomic theory calculations due to the large number of electrons and their complicated interactions. Conventional approaches such as calculations based on Cowan’s code are limited and require a large number of parameters for energy agreement. One promising approach is relativistic configuration-interaction and many-body perturbation theory (CI-MBPT) methods. We present CI-MBPT results for various atomic systems where this approach can lead to reasonable agreement: La I, La II, Th I, Th II, U I, Pu II. Among atomic properties, energies, g-factors, electric dipole moments, lifetimes, hyperfine structure constants, and isotopic shifts are discussed. While in La I and La II accuracy for transitions is better than that obtained with other methods, more work is needed for actinides.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Cross section of neutrons from the H 2 ( n , 2 n ) reaction at E n = 15 MeV

In this work, the double-differential cross section of the deuteron breakup reaction 2 H(n, 2n) has been studied experimentally with a neutron beam energy of 15 MeV. Special attention has been devoted to estimation of background condition and multiple scattering effect in the scattering sample. Experimental data have been compared with models based on phase-space approximation used in the ENDF/B-VIII.0 data library and in the MCNP neutron transport code, as well as with rigorous model based on Faddeev equations used for cross section evaluations in JENDL data library. It was found that experimental data are better reproduced by Faddeev model, however, the model overestimates data in the low-energy region of the neutron spectrum (<4 MeV)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

MRP and NEAMS Meeting - Current Modeling to Support Experiments

Heat-pipe cooled micro-reactor is currently being actively investigated in the NS&T department at INL. Given the high operation temperature and large temperature gradients in the system, the induced thermal stresses need to be well analyzed and kept within the safety margin. Two commercial FEM software (ABAQUS and STAR-CCM+) have been proposed to simulate and compare for both the temperature and stress distributions for the single heat pipe experiments, with the middle single heat pipe and 6 heaters surrounded. This presentation file is aimed to provide the latest updates for the researchers who is under NEAMS program and involved with the DOE micro-reactor program

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Code Standard Evaluation for ASME AG-1-2019, Section FL

ASME AG-1 Section FL Deep Bed Sand Filters includes a requirement to use Special Test Sections to provide a representative assessment of the filter efficiency. The SRPPF project design basis is to test the entire filter bed (versus a representative sample) following the same basis methodology to test the existing filters at the Savannah River Site. Using a combination of: 1) the filter testing experience at SRS since the early 1970’s, 2) the methodology found in ASME AG-1 for testing high capacity HEPA filter housings and 3) requirements for stack or duct sampling systems published by the Environmental Protection Agency (EPA), the Special Test Sections methodology for Deep Bed Sand (DBS) filters in ASME AG-1 Section FL has been replaced with a methodology to test the entire DBS filter system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Remote Radiation Sensing Using Aerial and Ground Platforms

Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Comprehensive Analysis of Streaming and Shutdown Dose Rate Experiments at JET with ORNL Fusion Neutronics Workflows

Current experimental fusion systems and conceptual designs of fusion pilot plants (FPPs) are growing in complexity and size. Several radiation metrics are crucial to the safe operation of fusion machines, including neutron flux streaming through openings and the shutdown dose rate (SDDR). Most current designs of advanced experimental fusion systems—and the most probable candidates for FPPs—are based on the tokamak concept, which is prone to neutron streaming through the myriad openings needed for diagnostic and support systems. SDDR is caused by decay gamma rays from radionuclides that become activated by neutrons during the operation of a fusion system that use deuterium-deuterium (DD), tritium-tritium, or deuterium-tritium plasma. Because computational tools have become essential for determining these radiation metrics, they must be validated against reliable and applicable experimental data. Experiments at the Joint European Torus (JET) provide a unique source of experimental data for validating computational tools and nuclear data used to determine SDDR and neutron fluxes in streaming-dominated geometries. Here, this paper presents the comprehensive analysis of the high-performance DD JET SDDR, and streaming experiments performed using Oak Ridge National Laboratory (ORNL) fusion workflows. The computational results were compared with experimental results that consist of online SDDR measurements with ionization chambers and neutron fluence streaming measurements using thermoluminescent detectors. The ratio of calculated-to-experimental SDDR values ranges from 0.6 to 2.5, and the streaming results range from 0.5 to 8.0. Future work will include analyzing the JET 2021 DTE2 campaign alongside the integration of the Shift Monte Carlo transport code into all ORNL fusion neutronics workflows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗