Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo codes”

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

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

The Energy Spectra of Heavy Nuclei Measured by the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon-borne experiment to measure the spectra and composition of primary cosmic rays in the region of total energy from 100 GeV to near 100 TeV for Z from 1 to 26. ATIC consists of a pixelated silicon matrix detector to measure charge plus a fully active BGO calorimeter, to measure energy, located below a carbon target interleaved with three layers of scintillator hodoscope. The ATIC instrument had a second (scientific) flight from McMurdo, Antarctica from 12/29/02 to 1/18/03, yielding 20 days of good data. The GEANT 3.21 Monte Carlo code with the QGSM event generator and the FLUKA code with the DPMJET-II event generator were used to convert energy deposition measurements to primary energy. We present the preliminary energy spectra for the abundant elements C, O, Ne, Mg, Si and Fe and compare them with the results of the first (test) flight of ATIC in 2000-01 and with results from the HEAO-3 and CRN experiments.

Panov, A. D.↗
NASA NTRS2004

MGGPOD: a Monte Carlo Suite for Modeling Instrumental Line and Continuum Backgrounds in Gamma-Ray Astronomy

Intense and complex instrumental backgrounds, against which the much smaller signals from celestial sources have to be discerned, are a notorious problem for low and intermediate energy gamma-ray astronomy (approximately 50 keV - 10 MeV). Therefore a detailed qualitative and quantitative understanding of instrumental line and continuum backgrounds is crucial for most stages of gamma-ray astronomy missions, ranging from the design and development of new instrumentation through performance prediction to data reduction. We have developed MGGPOD, a user-friendly suite of Monte Carlo codes built around the widely used GEANT (Version 3.21) package, to simulate ab initio the physical processes relevant for the production of instrumental backgrounds. These include the build-up and delayed decay of radioactive isotopes as well as the prompt de-excitation of excited nuclei, both of which give rise to a plethora of instrumental gamma-ray background lines in addition t o continuum backgrounds. The MGGPOD package and documentation are publicly available for download. We demonstrate the capabilities of the MGGPOD suite by modeling high resolution gamma-ray spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board Wind during 1995. The TGRS is a Ge spectrometer operating in the 40 keV to 8 MeV range. Due to its fine energy resolution, these spectra reveal the complex instrumental background in formidable detail, particularly the many prompt and delayed gamma-ray lines. We evaluate the successes and failures of the MGGPOD package in reproducing TGRS data, and provide identifications for the numerous instrumental lines.

Weidenspointner, G.↗
NASA NTRS1996

Orbiter Aerodynamic Acceleration Flight Measurements in the Rarefied-Flow Transition Regime

Acceleration data taken from the Orbital Acceleration Research Experiment (OARE) during reentry on STS-62 have been analyzed using calibration factors taken on orbit. This is the first Orbiter mission which collected OARE data during the Orbiter reentry phase. The data examined include the flight regime from orbital altitudes down to about 90 km which covers the free-molecule-flow regime and the upper altitude fringes of the rarefied-flow transition into the hypersonic continuum. Ancillary flight data on Orbiter position, orientation, velocity, and rotation rates have been used in models to transform the measured accelerations to the Orbiter center-of-gravity, from which aerodynamic accelerations along the Orbiter body axes have been calculated. Residual offsets introduced in the measurements by unmodeled Orbiter forces are identified and discussed. Direct comparisons are made between the OARE flight data and an independent micro-gravity accelerometer experiment, the High Resolution Accelerometer Package (HiRAP), which also obtained flight data on reentry during the mission down to about 95 km. The resulting OARE aerodynamic acceleration measurements along the Orbiter's body axis, aid the normal to axial acceleration ratio in the free-molecule-flow and transition-flow regimes are presented and compared with numerical simulations from three direct simulation Monte Carlo codes.

Blanchard, Robert C.↗
NASA NTRS2004

Density, Velocity and Ionization Structure in Accretion-Disc Winds

This was a project to exploit the unique capabilities of FUSE to monitor variations in the wind- formed spectral lines of the luminous, low-inclination, cataclysmic variables(CV) -- RW Sex. (The original proposal contained two additional objects but these were not approved.) These observations were intended to allow us to determine the relative roles of density and ionization state changes in the outflow and to search for spectroscopic signatures of stochastic small-scale structure and shocked gas. By monitoring the temporal behavior of blue-ward extended absorption lines with a wide range of ionization potentials and excitation energies, we proposed to track the changing physical conditions in the outflow. We planned to use a new Monte Carlo code to calculate the ionization structure of and radiative transfer through the CV wind. The analysis therefore was intended to establish the wind geometry, kinematics and ionization state, both in a time-averaged sense and as a function of time.

Sonneborn, George↗
NASA NTRS2002

Benchmark study for charge deposition by high energy electrons in thick slabs

The charge deposition profiles created when highenergy (1, 10, and 100 MeV) electrons impinge ona thick slab of elemental aluminum, copper, andtungsten are presented in this paper. The chargedeposition profiles were computed using existing representative Monte Carlo codes: TIGER3.0 (1D module of ITS3.0) and MCNP version 4B. The results showed that TIGER3.0 and MCNP4B agree very well (within ~20% of each other) in the majority of the problem geometry. The TIGER results were considered to be accurate based on previous studies. Thus, it was demonstrated that MCNP, with its powerful geometry capability and flexible source and tally options, could be used in calculations of electron charging in high energy electron-rich space radiation environments.

charging high energy electrons ITS MCNP↗
NASA NTRS2002

Benchmark study for total enery electrons in thick slabs

The total energy deposition profiles when highenergy electrons impinge on a thick slab of elemental aluminum, copper, and tungsten have been computed using representative Monte Carlo codes (NOVICE, TIGER, MCNP), and compared in this paper.

total enery deposition MCNP NOVICE TIGER↗
NASA NTRS2007

Comparison of Measured Dark Current Distributions with Calculated Damage Energy Distributions in HgCdTe

This paper presents a combined Monte Carlo and analytic approach to the calculation of the pixel-to-pixel distribution of proton-induced damage in a HgCdTe sensor array and compares the results to measured dark current distributions after damage by 63 MeV protons. The moments of the Coulombic, nuclear elastic and nuclear inelastic damage distributions were extracted from Monte Carlo simulations and combined to form a damage distribution using the analytic techniques first described in [1]. The calculations show that the high energy recoils from the nuclear inelastic reactions (calculated using the Monte Carlo code MCNPX [2]) produce a pronounced skewing of the damage energy distribution. While the nuclear elastic component (also calculated using the MCNPX) contributes only a small fraction of the total nonionizing damage energy, its inclusion in the shape of the damage across the array is significant. The Coulombic contribution was calculated using MRED [3-5], a Geant4 [4,6] application. The comparison with the dark current distribution strongly suggests that mechanisms which are not linearly correlated with nonionizing damage produced according to collision kinematics are responsible for the observed dark current increases. This has important implications for the process of predicting the on-orbit dark current response of the HgCdTe sensor array.

Marshall, C. J.↗
NASA NTRS2009

3D Visualization of Monte-Carlo Simulation's of HZE Track Structure and Initial Chemical Species

Heavy ions biophysics is important for space radiation risk assessment [1] and hadron-therapy [2]. The characteristic of heavy ions tracks include a very high energy deposition region close to the track (<20 nm) denoted as the track core, and an outer penumbra region consisting of individual secondary electrons (6-rays). A still open question is the radiobiological effects of 6- rays relative to the track core. Of importance is the induction of double-strand breaks (DSB) [3] and oxidative damage to the biomolecules and the tissue matrix, considered the most important lesions for acute and long term effects of radiation. In this work, we have simulated a 56Fe26+ ion track of 1 GeV/amu with our Monte-Carlo code RITRACKS [4]. The simulation results have been used to calculate the energy depiction and initial chemical species in a "voxelized" space, which is then visualized in 3D. Several voxels with dose >1000 Gy are found in the penumbra, some located ~0.1 mm from the track core. In computational models, the DSB induction probability is calculated with radial dose [6], which may not take into account the higher RBE of electron track ends for DSB induction. Therefore, these simulations should help improve models of DSB induction and our understanding of heavy ions biophysics.

Plante, Ianik↗
NASA NTRS2010

Meeting Radiation Protection Requirements and Reducing Spacecraft Mass - A Multifunctional Materials Approach

Both crew and radio-sensitive systems, especially electronics must be protected from the effects of the space radiation environment. One method of mitigating this radiation exposure is to use passive-shielding materials. In previous vehicle designs such as the International Space Station (ISS), materials such as aluminum and polyethylene have been used as parasitic shielding to protect crew and electronics from exposure, but these designs add mass and decrease the amount of usable volume inside the vehicle. Thus, it is of interest to understand whether structural materials can also be designed to provide the radiation shielding capability needed for crew and electronics, while still providing weight savings and increased useable volume when compared against previous vehicle shielding designs. In this paper, we present calculations and analysis using the HZETRN (deterministic) and FLUKA (Monte Carlo) codes to investigate the radiation mitigation properties of these structural shielding materials, which includes graded-Z and composite materials. This work is also a follow-on to an earlier paper, that compared computational results for three radiation transport codes, HZETRN, HETC, and FLUKA, using the Feb. 1956 solar particle event (SPE) spectrum. In the following analysis, we consider the October 1989 Ground Level Enhanced (GLE) SPE as the input source term based on the Band function fitting method. Using HZETRN and FLUKA, parametric absorbed doses at the center of a hemispherical structure on the lunar surface are calculated for various thicknesses of graded-Z layups and an all-aluminum structure. HZETRN and FLUKA calculations are compared and are in reasonable (18% to 27%) agreement. Both codes are in agreement with respect to the predicted shielding material performance trends. The results from both HZETRN and FLUKA are analyzed and the radiation protection properties and potential weight savings of various materials and materials lay-ups are compared.

Atwell, William↗
NASA NTRS2011

The Effects of Surface Roughness on the NEAR XRS Elemental Results: Monte-Carlo Modeling

The objective of the NEAR-Shoemaker X-ray Gamma-Ray Spec1roscopy ("XGRS") investigation was to determine the elemental composition of the near-Earth asteroid 433 Eros. The X-ray Spectrometer (XRS) system measured the characteristic fluorescence of six major elements (Mg, Al, Si, S, Ca, Fe) in the 1-10 keV energy range excited by the interaction of solar X-rays with the upper 100 microns of the surface of 433 Eros. Various investigators, using both laboratory experiments and computer simulations have established that X-ray fluorescent line ratios can be influenced by small-scale surface roughness at high incidence or emission angles. The effect on the line ratio is specific to the geometry, excitation spectrum, and composition involved, In general, however, the effect is only substantial for ratios of lines with a significant energy difference between them: Fe/Si and Ca/Si are much more likely to be affected than AI/Si or Mg/Si. We apply a Monte-Carlo code to the specific geometry and spectrum of a major NEAR XRS solar flare observation, using an H chondrite composition as the substrate. The seventeen most abundant elements were included in the composition model, from oxygen to titanium.

Lin, Lucy F.↗
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 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↗