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 163 records · Page 9

Generalization Equations for Machine Learners Based on Physical and Abstract Laws

The physical and abstract laws derived from the first principles have been recently exploited to customize and sharpen machine learning (ML) methods and also derive their generalization equations. These laws often encapsulate knowledge that complements datasets and ML models. We present a generic framework that uses these laws to provide ML codes that are transferable across multiple areas, including data transport infrastructures and thermal hydraulics analytics of nuclear reactors. By anchoring on datasets from these areas and the statistical generalization theory, we present a rigorous approach to co-develop ML solutions and the generalization equations that characterize them, by exploiting the structure and constraints from the laws. We present illustrative examples using practical problems from existing literature on the performance characterization of data transport infrastructures, and the sensor error and power level estimation in nuclear reactor systems using sensor measurements of primary and secondary coolant systems, respectively.

Rao, Nageswara↗

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis of NCERC Critical Experiments with ENDF/B Nuclear Data Libraries

Nuclear data (ND) libraries are the backbone of the nuclear industry, as they are the collections of tabulated probabilities that define sub-atomic particle interactions with matter. In the areas of criticality safety and experiments, no evaluated nuclear data files (ENDF) are more important than those containing neutron cross section data. It is upon these files, and accompanying radiation transport codes, that practitioners are enabled to safely design both subcritical and critical systems. Likewise, in a symbiotic fashion, it is the same critical assemblies which are used primarily to validate that the cross sections are correct. The ENDF/B library, the United States’ national library maintained by the National Nuclear Data Center (NNDC) at Brookhaven National Laboratory (BNL), is soon releasing a new version, ENDF/B-VIII.1. Prior to the official release, several beta versions of the library were prepared and tested in simulation suites. The work presented here are results comparing the newest ENDF/B beta release (ENDF/B-VIII.1b3) and ENDF/B-VIII.0 with recent experiments done at the National Criticality Experiments Research Center (NCERC) as well as correlated experiments from the Los Alamos Critical Experiments Facility (LACEF). These NCERC and LACEF experiments were performed in part to provide validation for various cross sections that were identified as insufficient in the ENDF/B-VIII.0. In particular, lead, copper, fluorine, and chlorine, as well as the major actinides, were targeted from the last decade of critical experiments.

97 MATHEMATICS AND COMPUTING↗

Integral Nuclear Data and Benchmarking Needs for Fusion Energy Systems

Fusion energy systems are currently being designed and optimized using radiation transport codes. To deal with the unique environment inside a fusion-based system, many of these designs incorporate novel materials able to withstand the high radiation fields, ensure adequate cooling and thermal protection, and produce tritium. Validation plays a vital role in building trust in the predictive power of these models and computational methods. Validation of a code consists of modeling documented real-world experiments and comparing the code-predicted response to the measured response. Adequate validation requires measured responses from real-world experiments, also known as integral data, that mimic the system being designed, including materials, impinging radiation, and temperature, among other variables. The most trusted integral data are experimental responses that have been through a rigorous benchmarking process that develops a recommended computational model and evaluates all experimental uncertainties. Finally, there are a few research groups around the world that have been producing integral data for fusion applications, but a substantial investment is needed to address the unique validation needs of the fusion community.

Fusion↗

Biophysical aspects of heavy ion interactions in matter

The biological effects of high energy, high charge nuclei (HZE particles) occupy a central role in the management of space radiation hazards due to galactic cosmic rays. For the energy range of interest, the mean free path for nuclear interactions of these heavy ions is comparable to the thickness of the material traversed, and a significant fraction of stopping particles will undergo a nuclear reaction with the nuclei of the stopping material. Transport methods for HZE particles are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The basic physical properties of HZE particles will be reviewed and illustrated with the results of nuclear fragmentation experiments performed with 670A MeV neon ions incident on a water absorber and with measurements of multiple Coulomb scattering of uranium beams in copper. Finally, the extent to which physical measurements yield radiobiological predictions is illustrated for the example of neon.

Schimmerling, Walter↗

COBRA-SFS Version 6.2.1 User Guide: A Thermal Hydraulic Analysis Code for Spent Fuel Storage and Transportation Casks

This document is an updated and reduced version of the complete documentation of the COBRA-SFS thermal analysis code (COolant Boiling in Rod Arrays – Spent Fuel Storage) (Michener et al. 1995). This manual provides the analyst with input instructions and fundamental guidelines for contrasting models using COBRA-SFS. COBRA-SFS is designed to analyze the dry storage and transportation casks that house spent or used nuclear fuel. The primary focus and application is on square light water reactor fuel assemblies, however the code also has the capability to model hexagonal assemblies and is will suited to any geometry that can be defined as an extruded 2-D map.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermochemically-Informed Mass Transport Model for Zr in U-Zr Fuel

Recent improvements to the coupled Thermochimica-MOOSE/BISON code system have enabled efficient calculations of species transport based on direct evaluation of composition and temperature dependent chemical potentials of the species. This presents an alternative to the traditional approach to species transport in nuclear fuels, which has been to employ a diffusion formulation that combines concentration-gradient driven Fickian diffusion with a Soret term based on a heat of transport fit to experimental data. Here we describe the application of the coupled code system to the diffusion of Zr in U-Zr metallic fuel. New classes implemented in BISON to solve this problem are documented. The Zr concentration profile after 50 years of diffusion is found to be strongly dependent on the assumptions made pertaining how to mobility is calculated in multi-phase regions of the fuel element. Two assumptions are compared (simple averaging of mobilities and using the majority phase mobility), and good qualitative agreement with experimental measurements is obtained using the majority phase assumption.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Facility Targeting, Protection and Mission Decision Making Using the VISAC Code

The Visual Interactive Site Analysis Code (VISAC) has been used by DTRA and several other agencies to aid in targeting facilities and to predict the associated collateral effects for the go, no go mission decision making process. VISAC integrates the three concepts of target geometric modeling, damage assessment capabilities, and an event/fault tree methodology for evaluating accident/incident consequences. It can analyze a variety of accidents/incidents at nuclear or industrial facilities, ranging from simple component sabotage to an attack with military or terrorist weapons. For nuclear facilities, VISAC predicts the facility damage, estimated downtime, amount and timing of any radionuclides released. Used in conjunction with DTRA's HPAC code, VISAC also can analyze transport and dispersion of the radionuclides, levels of contamination of the surrounding area, and the population at risk. VISAC has also been used by the NRC to aid in the development of protective measures for nuclear facilities that may be subjected to attacks by car/truck bombs.

Morris, Robert H.↗

Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length

A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NucMesh: nuclear reactor geometry creation and mesh generation module in NEMoSys

NucMesh is a parameterized geometry and mesh generator for nuclear reactors developed within the Nuclear Energy Modeling System NEMoSys at Illinois Rocstar. NEMoSys is a platform developed for mesh generation, adaptive refinement, and solution verification. NucMesh is implemented to be generalized and extensible with a robust computer aided design engine and multiple mesh generation algorithms for unstructured triangular, quad-dominant, and structured quadrilateral meshing. In this paper, we present the geometric and meshing features of NucMesh. Geometrically objects are constructed bottom-to-top and overlaps are addressed automatically. A sophisticated object tracking algorithm prevents data from being lost for segmented objects. We discuss the primitive objects of circle and polygons that constitute the module and show how they are used with example inputs. Arrays of primitives and arrays of arrays are utilized to build large assemblies of objects. The concept of saved objects is discussed to demonstrate how repetitive objects can be reused easily and augmented in place. Three dimensional meshes can be obtained through mesh extrusion where all materials and side sets are extended to three dimensions. We show that side sets can be defined nearly anywhere within the geometry and can then be applied to the mesh. Finally, example reactor meshes are demonstrated for the Idaho National Laboratory Advanced Test Reactor and Los Alamos National Lab Empire reactor, both of which use control drums that NucMesh handles easily. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

(U) Transport Corrections Implemented in SENSMG

Truncating the spherical harmonics expansion of the neutron scattering source in the Boltzmann transport equation leads to truncation error. Methods for approximately correcting for this truncation error are called transport corrections. The PARTISN multigroup discrete ordinates neutron transport code has three transport correction options. These options can be specified in the SENSMG multigroup neutron sensitivity code and passed to PARTISN for the neutron transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements of the secondary particle energy spectra in the Space Shuttle

Measurements of the energy spectra of secondary particles produced by galactic cosmic rays and trapped protons due to the nuclear interactions of these particles with the Shuttle shielding provide a powerful tool for validating radiation transport codes. A code validated in this way can be used to better estimate the dose and dose equivalent to body organs, measurements that cannot be made directly. The principal cause of single event upsets in electronic devices in the region of the South Atlantic Anomaly is secondary particles, and even in the region of galactic cosmic radiation a significant fraction is produced by secondary particles. In this paper, we describe the first direct measurements of the energy spectra of secondary protons, deuterons, tritons, 3He and 4He produced by galactic cosmic rays inside the Space Shuttle using a charged particle spectrometer. A comparison of these spectra with radiation transport code HZETRN showed reasonably good agreement for secondary protons. However, the code seriously underestimated the flux of all other light ions. The code has been modified to include pick-up and knock-on processes. The modified code leads to good agreement for deuterons and 3He but not for other light ions. This revised code leads to about 10% higher dose equivalent than the original code under moderate shielding, if we assume that higher charge ion fluxes are correctly predicted by the model.

short duration↗

Projectile and Lab Frame Differential Cross Sections for Electromagnetic Dissociation

Differential cross sections for electromagnetic dissociation in nuclear collisions are calculated for the first time. In order to be useful for three - dimensional transport codes, these cross sections have been calculated in both the projectile and lab frames. The formulas for these cross sections are such that they can be immediately used in space radiation transport codes. Only a limited amount of data exists, but the comparison between theory and experiment is good.

Norbury, John W.↗

Lunar Surface Reactor Shielding Study

Nuclear reactor system could provide power to support a long term human exploration to the moon. Such a system would require shielding to protect astronauts from its emitted radiations. Shielding studies have been performed for a Gas Cooled Reactor (GCR) system because it is considered to be the most suitable nuclear reactor system available for lunar exploration, based on its tolerance of oxidizing lunar regolith and its good conversion efficiency (Wright, 2003). The goals of the shielding studies were to provide optimal material shielding configuration that reduces the dose (rem) to the required level in order to protect astronauts, and to estimate the mass of regolith that would provide an equivalent protective effect if it were used as the shielding material. All calculations were performed using MCNPX code, a Monte Carlo transport code.

Monte Carlo N-Particle eXtended (MCNPX)↗

ANTARES: a new IRSN computational chain for safety analysis

Since some years, the IRSN has developed a new project to simulate accurately the interaction between neutronic and thermal hydraulic phenomena. This project called ANTARES (Advanced Neutronics and Thermal-hydraulic for the Analysis of Reactor Safety) couples the thermal hydraulic code CATHARE-3 with the 3D neutronic nodal code PARCS. The purpose of this project is to take the advantage of the complementary of these 2 codes, by the knowledge of the 3D core power distribution in CATHARE-3 (computed by PARCS) and improving the thermal hydraulic boundaries conditions in PARCS. In this way, the IRSN diversifies and enlarges its available tools to perform safety analysis with improved accuracy. The current technique usually adopted in France for the safety demonstrations for Reactivity Insertion Accident (RIA) is the so-called 'conservative' approach, which consists of reducing all the feedback (Doppler and moderator effects) and in modifying some physical quantities in such a way to increase a power peak. For this reason, these facilities ('penalties') have been implemented in ANTARES. In this paper we will give two examples of accidental transients that can be simulated with ANTARES: a Main Steam Line Break (MSLB) and an inadvertent boron dilution event. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Threat Sources for Creating Synthetic Urban Search Data

Equivalent point source energy emission distributions were computed for various threat sources for use in simulating the detector responses for urban search scenarios. The sources include standard isotopic sources used in detector testing, medical and industrial sources occasionally encountered in urban searches, and several types of special nuclear materials. Most of the equivalent point source distributions represent small sources inside some amount of shielding, but the special nuclear material sources represent volumetrically distributed sources in spheres of metal. Text-based inputs for emission distributions are available for the Monte Carlo transport codes Monte Carlo N-Particle, SCALE/MAVRIC, and Omnibus/Shift, any of which can easily be converted to other formats. These sources were developed for use in the Radiological Anomaly Detection and Identification (RADAI) project and the follow-on project, the RADAI-Extended project, sponsored by the National Nuclear Security Administration Office of Defense Nuclear Nonproliferation Research and Development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Overview of the Graphical User Interface for the GERM Code (GCR Event-Based Risk Model

The descriptions of biophysical events from heavy ions are of interest in radiobiology, cancer therapy, and space exploration. The biophysical description of the passage of heavy ions in tissue and shielding materials is best described by a stochastic approach that includes both ion track structure and nuclear interactions. A new computer model called the GCR Event-based Risk Model (GERM) code was developed for the description of biophysical events from heavy ion beams at the NASA Space Radiation Laboratory (NSRL). The GERM code calculates basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at NSRL for the purpose of simulating space radiobiological effects. For mono-energetic beams, the code evaluates the linear-energy transfer (LET), range (R), and absorption in tissue equivalent material for a given Charge (Z), Mass Number (A) and kinetic energy (E) of an ion. In addition, a set of biophysical properties are evaluated such as the Poisson distribution of ion or delta-ray hits for a specified cellular area, cell survival curves, and mutation and tumor probabilities. The GERM code also calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle. The contributions from primary ion and nuclear secondaries are evaluated. The GERM code accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections, and has been used by the GERM code for application to thick target experiments. The GERM code provides scientists participating in NSRL experiments with the data needed for the interpretation of their experiments, including the ability to model the beam line, the shielding of samples and sample holders, and the estimates of basic physical and biological outputs of the designed experiments. We present an overview of the GERM code GUI, as well as providing training applications.

Kim, Myung-Hee↗

A Prototype Thick-Target Bremsstrahlung Model with Angularly-Dependent Emission in the MCNP6 ® Code

This document summarizes the current thick-target bremsstrahlung (TTB) model in MCNP and provides test results for an alternative implementation to improve the accuracy with reduced cost compared to full electron transport. It has been observed that the current TTB model produces inaccurate results in problems where the medium is thick with respect to electrons, but the photon distribution in the problem has a strong directionality. An example of such a simulation is detectors surrounding a metal target irradiated with a radiographic beam of high energy photons. The primary cause of this discrepancy is the current TTB method emits all bremsstrahlung photons in the same direction as the primary electron produced from each (γ, e ± ) interaction, leading to artificially forward peaked photon distributions for intermediate to high-energy incident photons. To improve the TTB model, we have implemented an angularly-dependent TTB model in a developer version of the MCNP6 ® code; for developers, this was done on the branch prototype/angular_ttb in the mcnp6 repo on bitbucket. The angularly-dependent TTB model accounts for the energy and scattering of electrons as they slow down in the current material, but does not sample the computationally expensive energy straggling, secondary electron events, and tracking electrons; this approach is significantly less computationally expensive than full electron transport and can be comparable to the original TTB method for problems with sufficiently complex materials and geometry. To evaluate the method, we have modeled a simple problem of a beam of 5 MeV photons incident on a sphere of plutonium surrounded by detectors at different deflection angles. For this problem, the angularly-dependent TTB produces a photon flux within 8.0% for a 90 degree deflection angle and 0.8% along the beam axis, as compared to the electron transport solution. This is an improvement compared to a 43% and 81% discrepancy with the original TTB method, respectively. The rest of this work includes the following: the first section details the current TTB treatment in MCNP, which has not been well documented elsewhere. Then, the modified TTB algorithm is detailed and the approximations compared to the condensed history algorithm are compared. Results are given comparing the two TTB methods to the condensed history transport algorithm. The appendix includes details for code developers on relevant electron transport implementation details and potential code improvements for future work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗