SEARCH · Search NASA

Results for “Monte Carlo neutronics 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Coupled Reactor and Engine Nuclear Thermal Propulsion Modeling Methodology

Jacob Stonehill↗

NASA NTRSDate not supplied

Coupled Reactor and Engine Nuclear Thermal Propulsion Modeling Methodology

The design and development process of a Nuclear Thermal Propulsion (NTP) system requires extensive multiphysics modeling to couple the neutron physics and thermal feedback effects to determine the reactor’s power shape. Propulsion system performance codes utilize this power shape to determine NTP key performance parameters. While the power shape is heavily dependent on the temperature profile and geometry of the reactor, many analyses either assume a constant power shape, or use neutronics analysis to determine a power shape for a specific reactor configuration. The development of a coupling interface for a propulsion system performance code and a Monte Carlo neutron transport code (OpenMC) allows for the reactor power shape to be calculated in an iteration loop. The interface utilizes a file share system to transfer geometry dimensions, temperatures, and material identifiers to OpenMC, which is used to perform a neutron transport simulation of a design like the government Testing Reference Design reactor. The interface is then able to post-process the results from OpenMC and use the same file share system to share a power shape and other important neutron transport parameters to the system performance code. Initial results show that neglecting the changes to power shape when comparing reactor configurations can yield inaccurate results. Furthermore, utilizing propellants other than hydrogen gas can cause significant changes to the power shape, and thus, the thermal performance of a specific reactor design. This methodology is being expanded to allow for multiple families of NTP reactors to be analyzed, including block moderator, particle bed, and NERVA-derived reactors.

multiphysics coupling↗

NASA NTRSDate not supplied

Coupled Reactor and Engine Nuclear Thermal Propulsion Modeling Methodology

multiphysics coupling↗

NASA NTRS1962

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗

NASA NTRS1969

Fast neutron streaming in helical ducts

Fast neutron transport characteristics of helical ducts using Monte Carlo code FASTER, computing neutron fluxes and attenuation

Jordan, T. M.↗

NASA NTRS2004

High Altitude Radiations Relevant to the High Speed Civil Transport (HSCT)

The Langley Research Center (LaRC) performed atmospheric radiation studies under the SST development program in which important ionizing radiation components were measured and extended by calculations to develop the existing atmospheric ionizing radiation (AIR) model. In that program the measured neutron spectrum was limited to less than 10 MeV by the available 1960-1970 instrumentation. Extension of the neutron spectrum to high energies was made using the LaRC PROPER-3C monte carlo code. It was found that the atmospheric neutrons contributed about half of the dose equivalent and approximately half of the neutron contribution was from high energy neutrons above 10 MeV. Furthermore, monte carlo calculations of solar particle events showed that potential exposures as large as 10-100 mSv/hr may occur on important high latitude routes but acceptable levels of exposure could be obtained if timely descent to subsonic altitudes could be made. The principal concern was for pregnant occupants onboard the aircraft. As a result of these studies the FAA Advisory Committee on the Radiobiological Aspects of the SST recommended: 1. Crew members will have to be informed of their exposure levels 2. Maximum exposures on any flight to be limited to 5 mSv 3. Airborne radiation detection devices for total exposure and exposure rates 4. Satellite monitoring system to provide SST aircraft real-time information on atmospheric radiation levels for exposure mitigation 5. A solar forecasting system to warn flight operations of an impending solar event for flight scheduling and alert status. These recommendations are a reasonable starting point to requirements for the HSCT with some modification reflecting new standards of protection as a result of changing risk coefficients.

Wilson, J. W.↗

NASA NTRS1972

MORSE Monte Carlo shielding calculations for the zirconium hydride reference reactor

Verification of DOT-SPACETRAN transport calculations of a lithium hydride and tungsten shield for a SNAP reactor was performed using the MORSE (Monte Carlo) code. Transport of both neutrons and gamma rays was considered. Importance sampling was utilized in the MORSE calculations. Several quantities internal to the shield, as well as dose at several points outside of the configuration, were in satisfactory agreement with the DOT calculations of the same.

Burgart, C. E.↗

NASA NTRS1972

Radiation transport codes for potential applications related to radiobiology and radiotherapy using protons, neutrons, and negatively charged pions

Several Monte Carlo radiation transport computer codes are used to predict quantities of interest in the fields of radiotherapy and radiobiology. The calculational methods are described and comparisions of calculated and experimental results are presented for dose distributions produced by protons, neutrons, and negatively charged pions. Comparisons of calculated and experimental cell survival probabilities are also presented.

Armstrong, T. W.↗

NASA NTRS1997

Earth and Planetary Science Letters

Systematic measurements of the concentrations of cosmogen Ca-41 (half-life = 1.04 x 10(exp 5) yr) in the Apollo 15 long core 15001-15006 were performed by accelerator mass spectroscopy. Earlier measurements of cosmogenic Be-10, C-14, Al-26, Cl-36, and Mn-53 in the same core have provided confirmation and improvement of theoretical models for predicting production profiles of nuclides by cosmic ray induced spallation in the Moon and large meteorites. Unlike these nuclides, Ca-40 in the lunar surface is produced mainly by thermal neutron capture reactions on Ca-40. The maximum production of Ca-41, about 1 dpm/g Ca, was observed at a depth in the Moon of about 150 g/sq cm. For depths below about 300 g/sq cm, Ca-41 production falls off exponentially with an e-folding length of 175 g/sq cm. Neutron production in the Moon was modeled with the Los Alamos High Energy Transport Code System, and yields of nuclei produced by low-energy thermal and epithermal neutrons were calculated with the Monte Carlo N-Particle code. The new theoretical calculations using these codes are in good agreement with our measured Ca-41 concentrations as well as with Co-60 and direct neutron fluence measurements in the Moon.

K. Nishiizumi↗

NASA NTRS2011

Use of Apollo 17 Epoch Neutron Spectrum as a Benchmark in Testing LEND Collimated Sensor

The Apollo 17 neutron experiment LPNE provided a unique set of data on production of neutrons in the Lunar soil bombarded by Galactic Cosmic Rays (GCR). It serves as valuable "ground-truth" in the age of orbital remote sensing. We used the neutron data attributed to Apollo 17 epoch as a benchmark for testing the LEND's collimated sensor, as introduced by the geometry of collimator and efficiency of He3 counters. The latter is defined by the size of gas counter and pressure inside it. The intensity and energy spectrum of neutrons escaping the lunar surface are dependent on incident flux of Galactic Cosmic Rays (GCR) whose variability is associated with Solar Cycle and its peculiarities. We obtain first the share of neutrons entering through the field of view of collimator as a fraction of the total neutron flux by using the angular distribution of neutron exiting the Moon described by our Monte Carlo code. We computed next the count rate of the 3He sensor by using the neutron energy spectrum from McKinney et al. [JGR, 2006] and by consider geometry and gas pressure of the LEND sensor. Finally the neutron count rate obtained for the Apollo 17 epoch characterized by intermediate solar activity was adjusted to the LRO epoch characterized by low solar activity. It has been done by taking into account solar modulation potential, which affects the GCR flux, and in turn changes the neutron albedo flux.

Chin, Gordon↗

NASA NTRS2015

Using Neutron Spectroscopy to Constrain the Composition and Provenance of Phobos and Deimos

The origin of the Martian moons Phobos and Deimos is obscure and enigmatic. Hypotheses include the capture of asteroids originally from the outer main belt or beyond, residual material left over from Mars' formation, and accreted ejecta from a large impact on Mars, among others. Measurements of reflectance spectra indicate a similarity to dark, red D-type asteroids, but could indicate a highly space-weathered veneer. Here we suggest a way of constraining the near-surface composition of the two moons, for comparison to known meteoritic compositions. Neutron spectroscopy, particularly the thermal and epithermal neutron flux, distinguishes clearly between various classes of meteorites and varying hydrogen (water) abundances. Perhaps most surprising of all, a rendezvous with Phobos or Deimos is not necessary to achieve this. A low-cost mission based on the LADEE spacecraft design in an eccentric orbit around Mars can encounter Phobos every 2 weeks. As few as five flyby encounters at speeds of 2.3 kilometers per second and closest-approach distance of 3 kilometers provide sufficient data to distinguish between ordinary chondrite, water-bearing carbonaceous chondrite, ureilite, Mars surface, and aubrite compositions. A one-Earth year mission design includes many more flybys at lower speeds and closer approach distances, as well as similar multiple flybys at Deimos in the second mission phase, as described in the Phobos And Deimos Mars Environment (PADME) mission concept. This presentation will describe the expected thermal and epithermal neutron fluxes based on MCNP6 (Monte Carlo N (i.e. Neutron)-Particle transport code (version 6) simulations of different meteorite compositions and their uncertainties.

thtermal and epithermal neutron fluxes↗

NASA NTRS2017

Neutron Angular Scatter Effects in 3DHZETRN: Quasi-Elastic

The current 3DHZETRN code has a detailed three dimensional (3D) treatment of neutron transport based on a forward/isotropic assumption and has been compared to Monte Carlo (MC) simulation codes in various geometries. In most cases, it has been found that 3DHZETRN agrees with the MC codes to the extent they agree with each other. However, a recent study of neutron leakage from finite geometries revealed that further improvements to the 3DHZETRN formalism are needed. In the present report, angular scattering corrections to the neutron fluence are provided in an attempt to improve fluence estimates from a uniform sphere. It is found that further developments in the nuclear production models are required to fully evaluate the impact of transport model updates. A model for the quasi-elastic neutron production spectra is therefore developed and implemented into 3DHZETRN.

Wilson, John W.↗

NASA NTRS2015

Space Radiation Transport Code Development: 3DHZETRN

The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.

Wilson, John W.↗

NASA NTRS1984

Vectorized Monte Carlo methods for reactor lattice analysis

Some of the new computational methods and equivalent mathematical representations of physics models used in the MCV code, a vectorized continuous-enery Monte Carlo code for use on the CYBER-205 computer are discussed. While the principal application of MCV is the neutronics analysis of repeating reactor lattices, the new methods used in MCV should be generally useful for vectorizing Monte Carlo for other applications. For background, a brief overview of the vector processing features of the CYBER-205 is included, followed by a discussion of the fundamentals of Monte Carlo vectorization. The physics models used in the MCV vectorized Monte Carlo code are then summarized. The new methods used in scattering analysis are presented along with details of several key, highly specialized computational routines. Finally, speedups relative to CDC-7600 scalar Monte Carlo are discussed.

Brown, F. B.↗

NASA NTRS2018

Electromagnetic Chirps from Neutron Star-Black Hole Mergers

We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The gamma-ray light curve is modulated by a number of relativistic effects, including Doppler beaming and gravitational lensing. Because the photons originate from the inspiraling neutron star, the light curve closely resembles the corresponding gravitational waveform: a chirp signal characterized by a steadily increasing frequency and amplitude. We propose to search for these electromagnetic chirps using matched filtering algorithms similar to those used in LIGO data analysis.

gamma-ray burst↗

NASA NTRS2000

Materials for Low-Energy Neutron Radiation Shielding

Various candidate aircraft and spacecraft materials were analyzed and compared in a low-energy neutron environment using the Monte Carlo N-Particle (MCNP) transport code with an energy range up to 20 MeV. Some candidate materials have been tested in particle beams, and others seemed reasonable to analyze in this manner before deciding to test them. The two metal alloys analyzed are actual materials being designed into or used in aircraft and spacecraft today. This analysis shows that hydrogen-bearing materials have the best shielding characteristics over the metal alloys. It also shows that neutrons above 1 MeV are reflected out of the face of the slab better by larger quantities of carbon in the material. If a low-energy absorber is added to the material, fewer neutrons are transmitted through the material. Future analyses should focus on combinations of scatterers and absorbers to optimize these reaction channels and on the higher energy neutron component (above 50 MeV).

Singleterry, Robert C., Jr.↗

NASA NTRS2024

The Role of a Neutron Component in the Photospheric Emission of Long Duration Gamma-Ray Burst Jets

Long-duration gamma-ray bursts (LGRBs), thought to be produced during core-collapse supernovæ, may have a prominent neutron component in the outflow material. If present, neutrons can change how photons scatter in the outflow by reducing its opacity, thereby allowing the photons to decouple sooner than if there were no neutrons present. Understanding the details of this process could therefore allow us to probe the central engine of LGRBs, which is otherwise hidden. Here, we present results of the photospheric emission from an LGRB jet, using a combination of relativistic hydrodynamic simulations and radiative transfer post-processing using the Monte Carlo Radiation Transfer (MCRaT) code. We control the size of the neutron component in the jet material by varying the equilibrium electron fraction Y e , and we find that the presence of neutrons in the GRB fireball affects the Band parameters α and E 0 , while the picture with the β parameter is less clear. In particular, the break energy E 0 is shifted to higher energies. Additionally, we find that increasing the size of the neutron component also increases the total radiated energy of the outflow across multiple viewing angles. Our results not only shed light on LGRBs, but are also relevant to short-duration gamma-ray bursts associated with binary neutron star mergers, due to the likelihood of a prominent neutron component in such systems.

Gamma-ray bursts↗

NASA NTRS2014

Path Toward a Unifid Geometry for Radiation Transport

The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry↗