Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Effect on the Lunar Exosphere of a CME Passage

It has long been recognized that solar wind bombardment onto exposed surfaces in the solar system will produce an energetic component to the exospheres about those bodies. Laboratory experiments have shown that the sputter yield can be noticeably increased in the case of a good insulating surface. It is now known that the solar wind composition is highly dependent on the origin of the particular plasma. Using the measured composition of the slow wind, fast wind, solar energetic particle (SEP) population, and coronal mass ejection (CME), broken down into its various components, we have estimated the total sputter yield for each type of solar wind. The heavy ion component, especially the He++ component, greatly enhances the total sputter yield during times when the heavy ion population is enhanced, most notably during a coronal mass ejection. To simulate the effect on the lunar exosphere of a CME passage past the Moon, we ran a Monte Carlo code for the species Na, K, Mg and Ca.

Killen, Rosemary M.↗
NASA NTRS2011

Effect on the Lunar Exosphere of a CME Passage

It has long been recognized that solar wind bombardment onto exposed surfaces in the solar system will produce an energetic component to the exospheres about those bodies. Laboratory experiments have shown that the sputter yield can be noticeably increased in the case of a good insulating surface. It is now known that the solar wind composition is highly dependent on the origin of the particular plasma. Using the measured composition of the slow wind. fast wind. solar energetic particle (SEP) population. and coronal mass ejection (CME), broken down into its various components, we have estimated the total sputter yield for each type of solar wind. The heavy ion component, especially the He(++) component, greatly enhances the total sputter yield during times when the heavy ion population is enhanced, most notably during a coronal mass ejection. To simulate the etfect on the lunar exosphere of a CME passage past the Moon, we ran a Monte Carlo code for the species Na, K, Mg and Ca.

Killen, Roseamry M.↗
NASA NTRS2013

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗
NASA NTRS2014

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species with Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for near-equilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion.

Liechty, Derek S.↗
NASA NTRS2014

Deep Mapping of Small Solar System Bodies with Galactic Cosmic Ray Secondary Particle Showers

Galactic cosmic rays rain steadily from all directions onto asteroids and comets. The interaction of these high-energy ions produces a cascade of secondary particles, including muons, which can penetrate the solid interiors of small solar system bodies. Muons, which are produced in abundance in Earth's atmosphere, have been used to image large structures on Earth, including the Great Pyramids and the interior of volcanoes. In this study, we demonstrate that the transmitted flux of muons is sensitive to the interior density structure of asteroids and comets, less than a few hundred meters in diameter. Muonography has the potential to fill a critical gap in our knowledge of the deep interiors of small bodies, providing information needed for planetary defense, in situ resource utilization, and planetary science. We use Monte Carlo codes (MCNPX and FLUKA), which accurately model galactic cosmic ray showers, to explore systematic variations in the production of muons in solid surfaces. Results of these calculations confirm the scaling of muon production in Earth's atmosphere to solid regolith materials, as predicted by a simple, semi-empirical model. Muons are primarily produced in the top meter of the regoliths of asteroids and comets. Their rate of production is over three orders of magnitude lower than in Earth's atmosphere and depends strongly on regolith density. In practice, the use of muonography to characterize the interiors of small solar system bodies must overcome their low rate of production and their dependence on regolith density, which can vary over the surface of asteroids and comets. We show that interior contrast can be resolved using a muon telescope (hodoscope) with about 1 sq m aperture with integration times ranging from hours to weeks. Design concepts for a practical hodoscope that could be deployed in situ or on an orbiting spacecraft, are described. Regolith density within the top meter of an asteroid can be determined from radar observations. A concept for a pilot mission that combines remote radar measurements with in situ muonography of a near-Earth asteroid is presented. Perceived challenges and next steps for the development of the concept are described.

Prettyman, T. H.↗
NASA NTRS2015

The Distribution and Annihilation of Dark Matter Around Black Holes

We use a Monte Carlo code to calculate the geodesic orbits of test particles around Kerr black holes, generating a distribution function of both bound and unbound populations of dark matter (DM) particles. From this distribution function, we calculate annihilation rates and observable gamma-ray spectra for a few simple DM models. The features of these spectra are sensitive to the black hole spin, observer inclination, and detailed properties of the DM annihilation cross-section and density profile. Confirming earlier analytic work, we find that for rapidly spinning black holes, the collisional Penrose process can reach efficiencies exceeding 600%, leading to a high-energy tail in the annihilation spectrum. The high particle density and large proper volume of the region immediately surrounding the horizon ensures that the observed flux from these extreme events is non-negligible.

astroparticle physics↗
NASA NTRS2016

Carbon Back Sputter Modeling for Hall Thruster Testing

Lifetime requirements for electric propulsion devices, including Hall Effect thrusters, are continually increasing, driven in part by NASA's inclusion of this technology in it's exploration architecture. NASA will demonstrate high-power electric propulsion system on the Solar Electric Propulsion Technology Demonstration Mission (SEP TDM). The Asteroid Redirect Robotic mission is one candidate SEP TDM, which is projected to require tens of thousands of thruster life. As thruster life is increased, for example through the use of improved magnetic field designs, the relative influence of facility effects increases. One such effect is the sputtering and redeposition, or back sputter, of facility materials by the high energy thruster plumes. In support of wear testing for the Hall Effect Rocket with Magnetic Shielding (HERMeS) project, the back sputter from a Hall effect thruster plume has been modeled for the NASA Glenn Research Center's Vacuum Facility 5. The predicted wear at a near-worst case condition of 600 V, 12.5 kW was found to be on the order of 1 micron/kh in a fully carbon-lined chamber. A more detailed numerical Monte Carlo code was also modified to estimate back sputter for a detailed facility and pumping configuration. This code demonstrated similar back sputter rate distributions, but is not yet accurately modeling the magnitudes. The modeling has been benchmarked to recent HERMeS wear testing, using multiple microbalance measurements. These recent measurements have yielded values on the order of 1.5 - 2 micron/kh at 600 V and 12.5 kW.

Propulsion↗
NASA NTRS2016

Carbon Back Sputter Modeling for Hall Thruster Testing

In support of wear testing for the Hall Effect Rocket with Magnetic Shielding (HERMeS) program, the back sputter from a Hall effect thruster plume has been modeled for the NASA Glenn Research Centers Vacuum Facility 5. The predicted wear at a near-worst case condition of 600 V, 12.5 kW was found to be on the order of 3 4 mkhour in a fully carbon-lined chamber. A more detailed numerical monte carlo code was also modified to estimate back sputter for a detailed facility and pumping configuration. This code demonstrated similar back sputter rate distributions, but is not yet accurately modeling the magnitudes. The modeling has been benchmarked to recent HERMeS wear testing, using multiple microbalance measurements. These recent measurements have yielded values, on the order of 1.5- 2 microns/khour.

Facility↗
NASA NTRS2017

Materials Degradation in the Jovian Radiation Environment

The radiation environment of Jupiter represents a significant hazard for Europa Lander deorbit stage components, and presents a significant potential mission risk. The radiolytic degradation of ammonium perchlorate (AP) oxidizer in solid propellants may affect its properties and performance. The Monte Carlo code MONSOL was used for modeling of laboratory experiments on the electron irradiation of propellant samples. An approach for flattening dose profiles along the depth of irradiated samples is proposed. Depth-dose distributions produced by Jovian electrons in multi-layer slabs of materials are calculated. It is found that the absorbed dose in a particular slab is significantly affected by backscattered electrons and photons from neighboring slabs. The dose and radiolytic decomposition of AP crystals are investigated and radiation-induced chemical yields and weight percent of radical products are reported.

Miloshevsky, Gennady↗
NASA NTRS2020

Event Selection and Background Rejection in Time Projection Chambers Using Convolutional Neural Networks and a Specific Application to the AdEPT Gamma-ray Polarimeter Mission

The Advanced Energetic Pair Telescope gamma-ray polarimeter uses a time projection chamber for measuring pair production events and is expected to generate a raw instrument data rate four orders of magnitude greater than is transmittable with typical satellite data communications. GammaNet, a convolutional neural network, proposes to solve this problem by performing event classification on-board for pair production and background events, reducing the data rate to a level that can be accommodated by typical satellite communication systems. In order to train GammaNet, a set of 1.1 × 10(exp 6) pair production events and 10(exp 6) background events were simulated for the Advanced Energetic Pair Telescope using the Geant4 Monte Carlo code. An additional set of 10(exp 6) pair production and 10 background events were simulated to test GammaNet’s capability for background discrimination. With optimization, GammaNet has achieved the proposed background rejection requirements for Galactic Cosmic Ray proton events. Given the best case assumption for downlink speeds, signal sensitivity for pair production ranged between 1.1 ± 0.5% to 69 ± 2% for 5 and 250 MeV incident gamma rays. This range became 0.1 ± 0.1% to 17 ± 2% for the worst case scenario of downlink speeds. The application of a feature visualization algorithm to GammaNet demonstrated decreased response to electronic noise and events exiting or entering the frame and increased response to parallel tracks that are close in proximity. GammaNet has been successfully implemented and shows promising results.

Richard L. Garnett↗
NASA NTRSDate not supplied

Simulation of Vacuum Chamber Pressure Distribution with Surrogate Modeling and Uncertainty Quantification

A major challenge in understanding differences in electric propulsion performance in ground tests and in space operations concerns the pressure distribution within the test vacuum chamber. The chamber backpressure is much higher than experienced in space, modifying thruster performance and plume dynamics. Numerical simulation is a key element to determining the background conditions in non-ideal vacuum chamber environments. An important parameter for the accurate simulation of chamber backpressure is the probability that an atom will stick to a cryogenic panel used to pump away the plume gases. This quantity can be used to model vacuum pumps in particle-based kinetic numerical methods. In this work, a three-dimensional direct simulation Monte Carlo code is used to model neutral xenon atoms flowing from the anode of the H9 Hall Effect Thruster within the University of Michigan’s Large Vacuum Test Facility. Simulated pressures are compared with ion gauge pressure measurements to infer the effective sticking coefficient of the chamber’s vacuum pumps. A pressure predicting surrogate model is developed for inference of pump sticking coefficients and for uncertainty quantification. This information enables accurate and useful kinetic simulations of electric propulsion thruster plasma plumes in vacuum chambers.

DSMC↗
DOE OSTI2011

New Improved Nuclear Data for Nuclear Criticality and Safety

The Geel Electron Linear Accelerator (GELINA) was used to measure neutron total and capture cross sections of 182,183,184,186 W and 63,65 Cu in the energy range from 100 eV to {approx}200 keV using the time-of-flight method. GELINA is the only high-power white neutron source with excellent timing resolution and ideally suited for these experiments. Concerns about the use of existing cross-section data in nuclear criticality calculations using Monte Carlo codes and benchmarks were a prime motivator for the new cross-section measurements. To support the Nuclear Criticality Safety Program, neutron cross-section measurements were initiated using GELINA at the EC-JRC-IRMM. Concerns about data deficiencies in some existing cross-section evaluations from libraries such as ENDF/B, JEFF, or JENDL for nuclear criticality calculations were the prime motivator for new cross-section measurements. Over the past years many troubles with existing nuclear data have emerged, such as problems related to proper normalization, neutron sensitivity backgrounds, poorly characterized samples, and use of improper pulse-height weighting functions. These deficiencies may occur in the resolved- and unresolved-resonance region and may lead to erroneous nuclear criticality calculations. An example is the use of the evaluated neutron cross-section data for tungsten in nuclear criticality safety calculations, which exhibit discrepancies in benchmark calculations and show the need for reliable covariance data. We measured the neutron total and capture cross sections of 182,183,184,186 W and 63,65 Cu in the neutron energy range from 100 eV to several hundred keV. This will help to improve the representation of the cross sections since most of the available evaluated data rely only on old measurements. Usually, these measurements were done with poor experimental resolution or only over a very limited energy range, which is insufficient for the current application.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

THE CALCULATION OF BURNABLE POISON CORRECTION FACTORS FOR PWR FRESH FUEL ACTIVE COLLAR MEASUREMENTS

Verification of commercial low enriched uranium light water reactor fuel takes place at the fuel fabrication facility as part of the overall international nuclear safeguards solution to the civilian use of nuclear technology. The fissile mass per unit length is determined nondestructively by active neutron coincidence counting using a neutron collar. A collar comprises four slabs of high density polyethylene that surround the assembly. Three of the slabs contain {sup 3}He filled proportional counters to detect time correlated fission neutrons induced by an AmLi source placed in the fourth slab. Historically, the response of a particular collar design to a particular fuel assembly type has been established by careful cross-calibration to experimental absolute calibrations. Traceability exists to sources and materials held at Los Alamos National Laboratory for over 35 years. This simple yet powerful approach has ensured consistency of application. Since the 1980's there has been a steady improvement in fuel performance. The trend has been to higher burn up. This requires the use of both higher initial enrichment and greater concentrations of burnable poisons. The original analytical relationships to correct for varying fuel composition are consequently being challenged because the experimental basis for them made use of fuels of lower enrichment and lower poison content than is in use today and is envisioned for use in the near term. Thus a reassessment of the correction factors is needed. Experimental reassessment is expensive and time consuming given the great variation between fuel assemblies in circulation. Fortunately current modeling methods enable relative response functions to be calculated with high accuracy. Hence modeling provides a more convenient and cost effective means to derive correction factors which are fit for purpose with confidence. In this work we use the Monte Carlo code MCNPX with neutron coincidence tallies to calculate the influence of Gd{sub 2}O{sub 3} burnable poison on the measurement of fresh pressurized water reactor fuel. To empirically determine the response function over the range of historical and future use we have considered enrichments up to 5 wt% {sup 235}U/{sup tot}U and Gd weight fractions of up to 10 % Gd/UO{sub 2}. Parameterized correction factors are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Features of MCNP6 Relevant to Medical Radiation Physics

MCNP (Monte Carlo N-Particle) is a general-purpose Monte Carlo code for simulating the transport of neutrons, photons, electrons, positrons, and more recently other fundamental particles and heavy ions. Over many years MCNP has found a wide range of applications in many different fields, including medical radiation physics. In this presentation we will describe and illustrate a number of significant recently-developed features in the current version of the code, MCNP6, having particular utility for medical physics. Among these are major extensions of the ability to simulate large, complex geometries, improvement in memory requirements and speed for large lattices, introduction of mesh-based isotopic reaction tallies, advances in radiography simulation, expanded variance-reduction capabilities, especially for pulse-height tallies, and a large number of enhancements in photon/electron transport.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2011

The Rossi Alpha Validation Suite for MCNP

A Rossi α validation suite has been created for the MCNP Monte Carlo code. The suite includes 13 benchmarks based on specifications given in the International Handbook of Evaluated Criticality Safety Benchmark Experiments . The benchmarks are divided into four categories of fuel - 233 U, highly enriched uranium, intermediate enriched uranium, and plutonium. The benchmarks also cover fast, intermediate, and thermal spectra. Succinct descriptions are provided for each benchmark, along with computed values for Rossi α using ENDF/B-VI and ENDF/B-VII.0 nuclear data libraries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

OBJECT KINETIC MONTE CARLO SIMULATIONS OF RADIATION DAMAGE IN TUNGSTEN

We used our recently developed lattice-based object kinetic Monte Carlo code; KSOME [1] to carryout simulations of radiation damage in bulk tungsten at temperatures of 300, and 2050 K for various dose rates. Displacement cascades generated from molecular dynamics (MD) simulations for PKA energies at 60, 75 and 100 keV provided residual point defect distributions. It was found that the number density of vacancies in the simulation box does not change with dose rate while the number density of vacancy clusters slightly decreases with dose rate indicating that bigger clusters are formed at larger dose rates. At 300 K, although the average vacancy cluster size increases slightly, the vast majority of vacancies exist as mono-vacancies. At 2050 K no accumulation of defects was observed during irradiation over a wide range of dose rates for all PKA energies studied in this work.

Nandipati, Giridhar↗
DOE OSTI2016

NJOY2016

NJOY2016 is used to convert evaluations in the Evaluated Nuclear Data Files (ENDF) format into forms useful for practical applications such as fission and fusion reactor analysis, criticality safety, radiation shielding, nuclear waste management, nuclear medicine procedures, and more. It begins with the generation of pointwise libraries, including reaction and resonance reconstruction, Doppler broadening, radiation heating and damage, thermal scattering data, unresolved resonance data, and gas production. It reviews the production of libraries for the continuous-energy Monte Carlo code MCNP, multigroup neutron, photon, and particle cross sections and matrices, and photon interaction data. Uncertainty information for ENDF data is processed, including the capability for calculating covariances of resonance data, angular distributions, energy distributions, and radioactive nuclide production. NJOY can prepare and process thermal scattering data evaluations into application libraries. It can also plot cross sections, angular distributions and secondary emission spectra.

Kahler III, A. C. (Skip)↗
DOE OSTI2016

MARS15 Simulation of Radiation Environment at the ESS Linac

Comprehensive studies with the MARS15(2016) Monte-Carlo code are described on evaluation of prompt and residual radiation levels induced by nominal and accidental beam losses in the 5-MW, 2-GeV European Spallation Source (ESS) Linac. These are to provide a basis for radiation shielding design verification through the accelerator complex. The calculation model is based on the latest engineering design and includes a sophisticated algorithm for particle tracking in the machine RF cavities as well as a well-established model of the beam loss. Substantial efforts were put in solving the deep-penetration problem for the thick shielding around the tunnel with numerous complex penetrations. It allowed us to study in detail not only the prompt dose, but also component and air activation, radiation loads on the soil outside the tunnel, and skyshine studies for the complicated 3-D surface above the machine. Among the other things, the newest features in MARS15 (2016), such as a ROOT-based beamline builder and a TENDL-based event generator for nuclear interactions below 100 MeV, were very useful in this challenging application

Mokhov, N. V.↗