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 55 records · Page 3

Assessment of Modeling and Nuclear Data Needs for Active Neutron Interrogation

This document is the primary deliverable for a scoping study proposed to DOE National Laboratory Announcement Number LAB 19-2114 in the NNSA research area. The study supports user applications employing active neutron interrogation by providing a science plan to improve the modeling capability and the nuclear data that radiation transport codes use. Users rely on the accuracy of the elastic scattering and non-elastic cross-sections spanning thermal energies to 14 MeV (and higher in some cases) for modeling the neutron transport through complex geometries of materials potentially composed of many elements. While the elastic scattering cross-section data are accepted for all commonly occurring elements, the non-elastic cross-section data and the associated emission data include reaction channels that require attention. The study focused on the non-elastic reactions that emit secondary, also known as prompt, gammas with the premise that many users would benefit from improved modeling of these reactions. Many users develop material assay technologies based upon gamma signatures from radiative capture, inelastic scattering, and reactions on low-Z isotopes emitting multiple particles, so the nuclear data gaps, modeling deficiencies, and recommendations for addressing the shortfalls were assessed for these reactions. Fission gammas were excluded from this study because there are other efforts underway to address known shortfalls. Follow-on efforts that successfully execute the recommendations will tangibly improve to the ability to model gamma signatures and backgrounds for user applications, such as controlled substance detection, oil-well logging, and space exploration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

HZETRN: neutron and proton production in quasi-elastic scattering of GCR heavy-ions

The development of transport models for radiation shielding design and evaluation has provided a series of deterministic computer codes that describe galactic cosmic radiation (GCR), solar particle events, and experimental beams at particle accelerators. These codes continue to be modified to accommodate new theory and improvements to the particle interaction database (Cucinotta et al., 1994, NASA Technical Paper 3472, US Government Printing Office, Washington DC). The solution employed by the heavy-ion transport code HZETRN was derived with the assumption that nuclear fragments are emitted with the same velocity as the incident ion through velocity conserving nuclear interactions. This paper presents a version of the HZETRN transport code that provides a more realistic distribution of the energy of protons and neutrons emitted from GCR interactions in shields. This study shows that the expected GCR dose equivalent is lower than previously calculated for water shields that are less than 110 g cm-2 thick. Calculations of neutron energy spectra in low Earth orbit indicate substantial contributions from relativistic neutrons. c2001 Elsevier Science Ltd. All rights reseved.

NASA Center LaRC↗

Update of CEA DES Criticality-Safety Activities and Perspectives [Slides]

Neutronics staff at CEA DES are in charge of a whole set of modeling and simulation tools for neutronics: nuclear data evaluation and processing, transport codes development, calculation sequences edition, verification, validation and uncertainty analysis. Advanced user groups are in charge of specific nuclear analysis in reactor physics, fuel cycle, criticality-safety, radiation shielding and nuclear instrumentation for CEA and its partners, at all stages of nuclear facilities life cycle. CEA is willing to participate/collaborate on the following issues: detector design, experiments, analysis and evaluation, DH validation by combining TAGS measurements with Fission Yield variance-covariance matrices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GCR Transport in the Brain: Assessment of Self-Shielding, Columnar Damage, and Nuclear Reactions on Cell Inactivation Rates

Radiation shield design is driven by the need to limit radiation risks while optimizing risk reduction with launch mass/expense penalties. Both limitation and optimization objectives require the development of accurate and complete means for evaluating the effectiveness of various shield materials and body-self shielding. For galactic cosmic rays (GCR), biophysical response models indicate that track structure effects lead to substantially different assessments of shielding effectiveness relative to assessments based on LET-dependent quality factors. Methods for assessing risk to the central nervous system (CNS) from heavy ions are poorly understood at this time. High-energy and charge (HZE) ion can produce tissue events resulting in damage to clusters of cells in a columnar fashion, especially for stopping heavy ions. Grahn (1973) and Todd (1986) have discussed a microlesion concept or model of stochastic tissue events in analyzing damage from HZE's. Some tissues, including the CNS, maybe sensitive to microlesion's or stochastic tissue events in a manner not illuminated by either conventional dosimetry or fluence-based risk factors. HZE ions may also produce important lateral damage to adjacent cells. Fluences of high-energy proton and alpha particles in the GCR are many times higher than HZE ions. Behind spacecraft and body self-shielding the ratio of protons, alpha particles, and neutrons to HZE ions increases several-fold from free-space values. Models of GCR damage behind shielding have placed large concern on the role of target fragments produced from tissue atoms. The self-shielding of the brain reduces the number of heavy ions reaching the interior regions by a large amount and the remaining light particle environment (protons, neutrons, deuterons. and alpha particles) may be the greatest concern. Tracks of high-energy proton produce nuclear reactions in tissue, which can deposit doses of more than 1 Gv within 5 - 10 cell layers. Information on rates of cell killing from GCR, including patterns of cell killing from single particle tracks. can provide useful information on expected differences between proton and HZE tracks and clinical experiences with photon irradiation. To model effects on cells in the brain, it is important that transport models accurately describe changes in the GCR due to interactions in the cranium and proximate tissues. We describe calculations of the attenuated GCR particle fluxes at three dose-points in the brain and associated patterns of cell killing using biophysical models. The effects of the brain self-shielding and bone-tissue interface of the skull in modulating the GCR environment are considered. For each brain dose-point, the mass distribution in the surrounding 4(pi) solid angle is characterized using the CAM model to trace 512 rays. The CAM model describes the self-shielding by converting the tissue distribution to mass-equivalent aluminum, and nominal values of spacecraft shielding is considered. Particle transport is performed with the proton, neutron, and heavy-ion transport code HZETRN with the nuclear fragmentation model QMSFRG. The distribution of cells killed along the path of individual GCR ions is modeled using in vitro cell inactivation data for cells with varying sensitivity. Monte Carlo simulations of arrays of inactivated cells are considered for protons and heavy ions and used to describe the absolute number of cell killing events of various magnitude in the brain from the GCR. Included are simulations of positions of inactivated cells from stopping heavy ions and nuclear stars produced by high-energy ions most importantly, protons and neutrons.

Shavers, M. R.↗

Critical Simulation Pipeline for COG Suites [Poster]

The CRItical Simulation Pipeline (CRISP) is a Python package for automating validation of reactor criticality benchmarks. CRISP supplies COG—a multi-particle radiation transport code maintained by the Nuclear Criticality Safety Division—with a pipeline to calculate k eff performance for 400+ benchmark experiments with 3,400+ configurations from the International Criticality Safety Benchmark Evaluation Project (ICSBEP). The pipeline includes four stages: materials configuration, input card templating, cluster submission, and results analysis. CRISP includes a command-line interface to facilitate user interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Recent MCNP6 ® Code Developments and Improvements for Nuclear Engineering Applications [Slides]

The Los Alamos MCNP Monte Carlo radiation transport code has been the international gold standard for particle transport applications for over three decades. Many developments to the code have taken place with several significant new feature additions, major improvements, and enhancements to existing features. With significant institutional and programmatic investment in the code since the time of the last public release in 2018, important code development and infrastructure modernization has taken place and remains a high priority for all ongoing efforts across the code development team.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electromagnetic Dissociation Cross Sections using Weisskopf-Ewing Theory

It is important that accurate estimates of crew exposure to radiation are obtained for future long-term space missions. Presently, several space radiation transport codes exist to predict the radiation environment, all of which take as input particle interaction cross sections that describe the nuclear interactions between the particles and the shielding material. The space radiation transport code HZETRN uses the nuclear fragmentation model NUCFRG2 to calculate Electromagnetic Dissociation (EMD) cross sections. Currently, NUCFRG2 employs energy independent branching ratios to calculate these cross sections. Using Weisskopf-Ewing (WE) theory to calculate branching ratios, however, is more advantageous than the method currently employed in NUCFRG2. The WE theory can calculate not only neutron and proton emission, as in the energy independent branching ratio formalism used in NUCFRG2, but also deuteron, triton, helion, and alpha particle emission. These particles can contribute significantly to total exposure estimates. In this work, photonuclear cross sections are calculated using WE theory and the energy independent branching ratios used in NUCFRG2 and then compared to experimental data. It is found that the WE theory gives comparable, but mainly better agreement with data than the energy independent branching ratio. Furthermore, EMD cross sections for single neutron, proton, and alpha particle removal are calculated using WE theory and an energy independent branching ratio used in NUCFRG2 and compared to experimental data.

Adamczyk, Anne M.↗

Comparative study of decay heat calculations with FLUKA and MCNP/CINDER2008

In designing a high-power tungsten target, decay heat driven temperature rise in the spallation volume is a safety concern during maintenance and in loss of coolant accidents. As tungsten hydrates and becomes volatile in steam at above 800 ◦C, it is important to keep the target temperature below this threshold when active cooling is unavailable. Decay heat in a target is calculated with particle transport simulation codes combined with transmutation codes. The calculated decay heat usually differs depending on the nuclear cross sections and the decay particle transport models built in the code architecture. In this paper, we calculated decay heat of a water-cooled tungsten target using popular particle transport codes, FLUKA and MCNP6® paired with CINDER2008 and compared the results. The target-moderator-reflector (TMR) system is modeled with a water-cooled solid tungsten target, water premoderators, liquid hydrogen cold moderators and beryllium reflectors. Water-cooled stainless-steel shielding is modeled around the TMR system. The tungsten volume is clad with a thin layer of erosion/corrosion resistant material. This study provides information about the uncertainty range in decay heat prediction of high-power spallation targets for hazard analysis.

Lee, Yong Joong↗

Recent MCNP® Code Developments and Improvements for Nuclear Energy Applications [Slides]

The MCNP® code is a general-purpose radiation transport code developed at LANL over the past 46+ years. This code is capable of modeling the fundamental physics of particles as they move through and interact with materials. It is used in many LANL and worldwide radiation transport applications (see next slide) for example in the Fundamental nuclear physics and data experiments at LANSCE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Resonance Self-Shielding: Why it is so Important

This paper is one of a series that I am writing to document my 58 years of experience with ENDF and Neutron Transport calculations, beginning when I worked at the National Nuclear Data Center (NNDC), Brookhaven National Laboratory (BNL), from 1967 to 1972. During those years I was the head of the computer unit of NNDC, assigned to develop computer codes to pre-process, view and test ENDF/B data. Since then, I have continued to support the ENDF effort without any official position or monetary compensation, because I realized how important accurate nuclear data is for use in use in our Engineering applications. It is so important to realize that regardless of how accurate or even perfect our application codes may be to transport particles, without accurate nuclear data we are in a “Garbage In = Garbage Out” situation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of MCNP Training Modules for Safeguards Practitioners [Abstract]

The Monte-Carlo N-Particle (MCNP) software developed at LANL is the most widely used neutron transport code in the world. It is an essential tool for a variety of applications including detector development and design, nuclear fuel burnup simulation, criticality safety, and nondestructive assay system optimization. For this reason, it is indispensable within the safeguards and materials control & accountability (MC&A) communities. Multiple MCNP training courses have been created and taught over the last several decades by the MCNP development team at LANL, however there are no existing courses that cover specialized topics considered fundamental to NDA and safeguards models. To fill this gap, the MCNP team and Safeguards Science and Technology group at LANL have co-created a set of training modules customized to meet the specialized needs of the safeguards and MC&A communities. The basic modules cover concepts such as NDA system optimization, He-specific and other capture tallies, and tools for improved theoretical understanding. An advanced module was also created to cover topics including variance reduction for active interrogation simulations, use of the LANL MCNPTools post-processor, PTRAC (particle tracking) and list-mode data simulations, and fuel burnup simulations. The training modules teach to the latest and most state-of-the-art MCNP features and tools released by the development team at LANL and are intended to be taught jointly by the developers and safeguards experts. Ultimately, we hope that creation of these modules will serve to capture and convey the safeguards modeling and MCNP expertise at LANL, and that we will be able to share the modules more broadly with the MC&A and safeguards communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exascale models of stellar explosions: Quintessential multi-physics simulation

The ExaStar project aims to deliver an efficient, versatile, and portable software ecosystem for multi-physics astrophysics simulations run on exascale machines. The code suite is a component-based multi-physics toolkit, built on the capabilities of current simulation codes (in particular Flash-X and Castro), and based on the massively parallel adaptive mesh refinement framework AMReX. It includes modules for hydrodynamics, advanced radiation transport, thermonuclear kinetics, and nuclear microphysics. The code will reach exascale efficiency by building upon current multi- and many-core packages integrated into an orchestration system that uses a combination of configuration tools, code translators, and a domain-specific asynchronous runtime to manage performance across a range of platform architectures. The target science includes multi-physics simulations of astrophysical explosions (such as supernovae and neutron star mergers) to understand the cosmic origin of the elements and the fundamental physics of matter and neutrinos under extreme conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

DRiFT - Release 1.1.1: Organic Scintillators

DRiFT (a Detector Response Function Toolkit) is LANL-developed software that postprocesses output from the extensively validated radiation transport code, MCNP, and generates realistic nuclear instrumentation response. DRiFT is designed to be flexible, enabling users to specify detector type and many experimental settings, as well as accommodating the addition of their own desired features. Although DRiFT development has included scintillator, gas, and semiconductor features, the focus of this release is on organic scintillator and associated capabilities. Organic scintillators are widely used in the areas of nuclear safeguards and nuclear non-proliferation efforts. DRiFT has several diagnostic and detector physics features relevant to detailed scintillator simulations including: tracking source particle information, scintillation light production, the effects of PMT quantum efficiency and gain, and digitizer settings. Users can select responses from many scintillator and PMT types supported natively by DRiFT, or add their own by following the instructions in this document. We acknowledge that DRiFT is under active development, bug reports and general questions and comments should be directed to Madison Andrews, madison@lanl.gov. This manual is divided into four parts: I) An overview of DRiFT, including how to obtain and install the executable, II) A description of the detector physics related to scintillators available, III) a description of more general DRiFT features the user may find useful, and IV) a description of the test suite and examples made available with the code release.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DRiFT - Release 2.1.0: Organic Scintillators and Gas Detectors

DRiFT (a Detector Response Function Toolkit) is LANL-developed software that post-processes output from the extensively validated radiation transport code, MCNP, and generates realistic nuclear instrumentation response. DRiFT is designed to be flexible, enabling users to specify detector type and many experimental settings, as well as accommodating the addition of their own desired features. The focus of this release is on organic scintillator, gas detector, and associated capabilities, while semiconductor features are still under development. Organic scintillators are widely used in the areas of nuclear safeguards and nuclear non-proliferation efforts. DRiFT has several diagnostic and detector physics features relevant to detailed scintillator simulations including: tracking source particle information, scintillation light production, the effects of PMT quantum efficiency and gain, and digitizer settings. Users can select responses from many scintillator and PMT types supported natively by DRiFT, or add their own by following the instructions in this document. DRiFT also has several diagnostic and detector physics features relevant to neutron gas detector simulations including: gas detector wall effects, inactive areas at the end of detector tubes, and effects of a preamplifier. We acknowledge that DRiFT is under active development, bug reports and general questions and comments should be directed to Madison Andrews, madison@lanl.gov. This manual is divided into four parts: I) An overview of DRiFT, including how to obtain and install the executable, II) A description of the detector physics related to scintillators available, III) a description of more general DRiFT features the user may find useful, and IV) a description of the test suite and examples made available with the code release.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Subcritical Multiplication with a Fixed Source

In a subcritical, multiplying medium, the system multiplication describes the expected total number of neutrons created by a single source neutron. Subcriticality plays a large role in criticality safety and thus it is vital for the subcritical multiplication factor be accurate, especially as a system approaches criticality. This work examines the accuracy of calculating the system multiplication using the MCNP6.2 ® k-eigenvalue power iteration (KCODE) method when a fixed-point source is present in a multiplying medium, for near critical systems. This work compares the standard approach for calculating system multiplication, using the fixed-source calculational approach, to a new, single k-eigenvalue power iteration approach that incorporates a fixed-source component and a fission-source component into a single calculation. For the remainder of this paper, some theoretical background and numerical results for an approximate k eigenvalue approach, an accurate fixed-source approach and a new and more accurate k-eigenvalue approach to computing system multiplication are provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interpretation of TEPC Measurements in Space Flights for Radiation Monitoring

For the proper interpretation of radiation data measured in space, the results of integrated radiation transport models were compared with the tissue equivalent proportional counter (TEPC) measurements. TEPC is a simple, time-dependent approach to radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Qave(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Qave(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. The deviation of y from LET is caused by energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. The response distribution functions of the wall-less and walled TEPCs were calculated from Monte-Carlo track simulations. Using an integrated space radiation model (which includes the transport codes HZETRN and BRYNTRN, and the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions from Monte-Carlo track simulations, we compared model calculations with the walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. Good agreement was found for Qave(y) between the model and measured spectra from past NASA missions. The Qave(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Qave(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Qave(y). The GCR values are decreasing with the shield thickness. Our analysis of the measurements of TEPCs can be used for a proper interpretation of observed data of monitoring the space radiation environment.

Kim, Myung-Hee Y.↗

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗