Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear transport codes”

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

SAM Code Improvements for Molten Salt Reactor Transient and Species Transport Modeling

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes FY24 efforts in addressing the modeling gaps in SAM for molten salt reactor (MSR) applications. These efforts focused on the solver performance related to species transport, usability of the flowing decay heat and delayed neutron precursor capabilities, and modeling gaps related to fission product transport throughout the system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal Neutron Scattering Law (TNSL) Implementation and Testing in FUDGE

Thermal neutron scattering law processing capabilities have been recently implemented in the LLNL code FUDGE (For Updating Data and Generating Evaluations). FUDGE is now capable of producing processed TNSL data for use in Monte Carlo or deterministic transport. To test this new capability, LLNL scientists performed an extensive intercomparison between the Mercury, Ardra, COG and MCNP transport codes. This intercomparison helped probe differences between how TNSL data are handled by the transport codes and by two different processing codes (LLNL’s FUDGE and LANL’s NJOY). This report summarizes recent improvements in FUDGE TNSL capabilities as well as results from the code intercomparison.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SFCOMPO Database of Spent Nuclear Fuel Assay Data – the Next Frontier

SFCOMPO is the world’s largest database for measured spent nuclear fuel assay data. An international effort coordinated by the Nuclear Energy Agency (NEA) resulted in a significant expansion of the database and its release online in 2017 as a downloadable application. The SFCOMPO Technical Review Group (TRG) was recently formed under the direction of NEA’s Nuclear Science Committee/Working Party on Nuclear Criticality Safety and was mandated to maintain and further coordinate the development of SFCOMPO. This TRG is currently focused on (1) critical evaluation of the experimental assay data by independent experts and (2) development of benchmarks and benchmark models that can be applied to validate burnup codes. This will improve the quality and documentation of the experimental datasets and enable their use by the international community to support code validation for design and safety analysis of spent nuclear fuel transportation, storage, and repository applications. It follows the precedent and draws on the experience gained from similar NEA efforts in the International Reactor Physics Experiment Evaluation Project and the International Criticality Safety Benchmark Experiment Project. Ongoing SFCOMPO evaluations have served as a test bed to develop templates for documenting evaluations, develop review guidance, improve approaches for a global uncertainty analysis, and devise a strategy focused on providing practical information of highest value to the user community. The current effort, status, and associated challenges are discussed.

Ilas, Germina↗

Galactic cosmic ray exposure estimates for SAGE-3 mission in polar orbit

An analysis of the effects of galactic cosmic ray (GCR) exposures on charge-coupled devices (CCDs) was performed for the SAGE-III 5-year mission in sun-synchronous orbit between 1996 and 2001. A detailed environment model used in conjunction with a geomagnetic vertical cut-off code provides the predicted 5-year fluence of GCR ions. A computerized solid model of the spacecraft was used to define the effective shield thickness distribution around the CCD detector. The particle fluences at the detector location are calculated with the Langley heavy-ion transport code, and these fluences are used in conjunction with estimated nuclear stopping powers to evaluate dosimetric quantities related to the detector degradation. A previous study analyzing effects of trapped particle and solar flare protons indicated an approximate 20 percent reduction in detector sensitivity for the mission. The galactic cosmic ray contribution was thought to be relatively small and therefore was not previously analyzed. The present study provides quantification of the GCR effects, which are found to contribute less than 1 percent of the total environment degradation.

Nealy, John E.↗

Parameterizations of Pion Energy Spectrum in Nucleon-Nucleon Collisions

The effects of pion (PI) production are expected to play an important role in radiation exposures in the upper atmosphere or on the Martian surface. Nuclear databases for describing pion production are developed for radiation transport codes to support these studies. We analyze the secondary energy spectrum of pions produced in nucleon-nucleon (NN) collisions in the relativistic one-pion exchange model. Parametric formulas of the isospin cross sections for one-pion production channels are discussed and are used to renormalize the model spectrum. Energy spectra for the deuteron related channels (NN yields dPi) are also described.

Cucinotta, Franics A.↗

In vitro cell irradiation systems based on 210Po alpha source: construction and characterisation

One way of studying the risk to human health of low-level radiation exposure is to make biological experiments on living cell cultures. Two 210Po alpha-particle emitting devices, with 0.5 and 100 MBq activity, were designed and constructed to perform such experiments irradiating monolayers of cells. Estimates of dose rate at the cell surface were obtained from measurements by a PIPS alpha-particle spectrometer and from calculations by the SRIM 2000, Monte Carlo charged particle transport code. Particle fluence area distributions were measured by solid state nuclear track detectors. The design and dosimetric characterisation of the devices are discussed. c2002 Elsevier Science Ltd. All rights reserved.

Alpha Particles↗

Isotopic Dependence of GCR Fluence behind Shielding

In this paper we consider the effects of the isotopic composition of the primary galactic cosmic rays (GCR), nuclear fragmentation cross-sections, and isotopic-grid on the solution to transport models used for shielding studies. Satellite measurements are used to describe the isotopic composition of the GCR. For the nuclear interaction data-base and transport solution, we use the quantum multiple-scattering theory of nuclear fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, respectively. The QMSFRG model is shown to accurately describe existing fragmentation data including proper description of the odd-even effects as function of the iso-spin dependence on the projectile nucleus. The principle finding of this study is that large errors (+/-100%) will occur in the mass-fluence spectra when comparing transport models that use a complete isotopic-grid (approx.170 ions) to ones that use a reduced isotopic-grid, for example the 59 ion-grid used in the HZETRN code in the past, however less significant errors (<+/-20%) occur in the elemental-fluence spectra. Because a complete isotopic-grid is readily handled on small computer workstations and is needed for several applications studying GCR propagation and scattering, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

MCNP ® Code Version 6.3.1 Theory & User Manual

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.1, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of CASL VERA for BWR analysis and application to accident tolerant SiC/SiC channel box

Application of the Virtual Environment for Reactor Applications (VERA) to BWR analysis is assessed in this paper by comparing results to those calculated using other widely-used modeling tools, namely the U.S. Nuclear Regulatory Commission’s PARCS/PATHS and the Serpent Monte Carlo particle transport code. Additionally, VERA is used to calculate 3-D temperature and fast neutron flux distributions in silicon carbide (SiC) fiber-reinforced, SiC matrix composite (SiC/SiC) BWR channel boxes, which are being studied as an Accident Tolerant Fuel core structural material concept. We found that the code-to-code comparisons were favorable, and the SiC/SiC channel box evaluation demonstrates the many advanced modeling features VERA provides while also highlighting the non-uniformity in fast neutron flux distributions that can play a role in potential SiC/SiC channel box deformation. Traditional BWR analysis tools do not have the calculation fidelity necessary for coupled assessment of flux and temperature gradients in a SiC/SiC channel box. VERA is a state-of-the-art modeling environment that was developed to increase the safety and economic competitiveness of nuclear power through improved modeling accuracy. While VERA has already been deployed in the nuclear industry for PWR applications, the current study is a vital initial step in the extensive development, validation, and verification that VERA must go through to be useful for BWR applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extension of OpenMC for Fixed Source Transmutation Calculations

This report documents work performed under a Strategic Partnership Project between Argonne National Laboratory (ANL) and the United Kingdom Atomic Energy Authority (UKAEA). The overall goal of this project is to extend OpenMC [1], a community-developed Monte Carlo particle transport code, to be able to perform fixed-source transmutation calculations. In fission and fusion reactors, the high flux of energetic neutrons causes materials within the reactor to “transmute,” or undergo a nuclear reaction that results in the addition/removal of neutrons and protons from the nucleus of an atom. If subjected to these reactions for long enough, the overall composition and physical properties of the material itself begin to change as a result of transmutation. Such a feature is vital for predicting the decrease in tritium production rate within a breeder blanket during the lifetime of a fusion reactor. OpenMC is capable of simulating neutron transport in fission/fusion systems, thereby allowing it to estimate the flux that causes transmutation. It is also capable of solving the transmutation equations, which determine how the composition of a material changes over time due to neutron irradiation and radioactive decay. However, solving the transmutation equations was previously only possible when the source of neutrons came from a fission system. In a fusion system, the source of neutrons is typically determined by a separate code and then given as an input to the particle transport simulation. This is known as a fixed source calculation. Through this project, we have extended OpenMC to solve the transmutation equations for a fixed source calculation. Evaluating the change in material compositions due to transmutation and its effect on physical properties is of key importance to a range of engineering analyses for fission and fusion systems. For example, in a fusion reactor, estimating the dose rate at different physical locations resulting from irradiated materials in the reactor allows designers to ensure that workers are not exposed to doses beyond applicable regulations. In order to properly dispose of irradiated materials, designers also need to estimate the radiotoxicity, which again relies on knowledge of the material composition at some future time. The specific tasks for this project that were agreed to between ANL and UKAEA were as follows: 1. Make changes and additions in the openmc.deplete and related modules in OpenMC to support transmutation calculations following a fixed source transport simulation. 2. Make necessary changes to OpenMC to model transmutation due to an arbitrary set of reactions needed for fusion applications. Use this new capability to generate a depletion chain file based on the TENDL nuclear data library. 3. Improve the openmc.deplete module in OpenMC to keep track of gases produced as a by-product of nuclear reactions during transmutation calculations. 4. Validate the new capabilities by carrying out fixed-source transmutation calculations on a suitable benchmark problem using OpenMC and a comparable Monte Carlo neutron transport code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗

Nonperturbative methods in HZE ion transport

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport. The code is established to operate on the Langley Research Center nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code is highly efficient and compares well with the perturbation approximations.

Wilson, John W.↗

The FLUKA code for space applications: recent developments

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Non-NASA Center↗

Advances in Space Radiation Physics and Transport

The space radiation environment is a complex mixture of particle types and energies originating from sources inside and outside of the galaxy. These environments may be modified by the heliospheric and geomagnetic conditions as well as planetary bodies and vehicle or habitat mass shielding. In low Earth orbit (LEO), the geomagnetic field deflects a portion of the galactic cosmic rays (GCR) and all but the most intense solar particle events (SPE). There are also dynamic belts of trapped electrons and protons with low to medium energy and intense particle count rates. In deep space, the GCR exposure is more severe than in LEO and varies inversely with solar activity. Unpredictable solar storms also present an acute risk to astronauts if adequate shielding is not provided. Near planetary surfaces such as the Earth, moon or Mars, secondary particles are produced when the ambient deep space radiation environment interacts with these surfaces and/or atmospheres. These secondary particles further complicate the local radiation environment and modify the associated health risks. Characterizing the radiation fields in this vast array of scenarios and environments is a challenging task and is currently accomplished with a combination of computational models and dosimetry. The computational tools include models for the ambient space radiation environment, mass shielding geometry, and atomic and nuclear interaction parameters. These models are then coupled to a radiation transport code to describe the radiation field at the location of interest within a vehicle or habitat. Many new advances in these models have been made in the last decade, and the present review article focuses on the progress and contributions made by workers and collaborators at NASA in the same time frame. Although great progress has been made, and models continue to improve, significant gaps remain and are discussed in the context of planned future missions. Of particular interest is the juxtaposition of various review committee findings regarding the accuracy and gaps of combined space radiation environment, physics, and transport models with the progress achieved over the past decade. While current models are now fully capable of characterizing radiation environments in the broad range of forecasted mission scenarios, it should be remembered that uncertainties still remain and need to be addressed.

Space radiation↗

Evaluation of Attila and MCNP computational methods for dose and exposure estimation

Radiation transport calculations are often used to estimate dose or exposure to components and personnel surrounding a radiation source. The sources for these calculations are decaying radionuclides within various nuclear materials. Historically, dose calculations use MCNP (Monte Carlo N-Particle) transport code as the primary particle transport tool without a secondary computational tool to validate the results from the MCNP simulations [1]. The goal of this study is to make an independent check of the Monte Carlo solution from MCNP6 Version 6.2.1 with the discrete ordinates solution from Attila 10.2.0 Beta 3. As an example problem for this study, water-filled, stainless-steel vessels, modeled with an unstructured mesh (UM) with both MCNP and Attila [2], are exposed to 252Cf and 60Co point sources. This report also includes a discussion of the limitations of unstructured mesh in a MCNP calculation.

61 RADIATION PROTECTION AND DOSIMETRY↗

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗