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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Celeritas Midterm SciDAC Report

Celeritas is a new Monte Carlo (MC) code that helps satisfy the increasing demand for high energy physics (HEP) detector simulation, using Graphics Processing Unit (GPU) hardware on high performance computing (HPC) systems to model Large Hadron Collider (LHC) experiments and beyond. This report details the project’s progress midway through its SciDAC funding period, highlighting the first complete implementation of standard electromagnetic (EM) physics on GPUs, initial results for performance and scalability on Leadership Computing Facilities (LCFs), and preliminary integration into the CMS and ATLAS experiments. By integrating HEP domain knowledge with expertise in MC transport, Celeritas has catalyzed a shift in the HEP community’s perception of GPU platforms as the future for HPC simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

COG Software Architecture Design Description Document

This COG Software Architecture Design Description Document describes the organization and functionality of the COG Multiparticle Monte Carlo Transport Code for radiation shielding and criticality calculations, at a level of detail suitable for guiding a new code developer in the maintenance and enhancement of COG. The intended audience also includes managers and scientists and engineers who wish to have a general knowledge of how the code works. This Document is not intended for end-users.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Monte Carlo Study of Imaging Plate Response to Laser-Driven Aluminum Ion Beams

We measured the response of BAS-TR imaging plate (IP) to energetic aluminum ions up to 222 MeV, and compared it with predictions from a Monte Carlo simulation code using two different IP response models. Energetic aluminum ions were produced with an intense laser pulse, and the response was evaluated from cross-calibration between CR-39 track detector and IP energy spectrometer. For the first time, we obtained the response function of the BAS-TR IP for aluminum ions with a kinetic energy as high as 222 MeV. On close examination of the two IP response models, we confirm that the exponential model fits our experimental data better. Moreover, we find that the IP sensitivity in the exponential model is nearly constant in this energy range, suggesting that the response function can be determined even with little experimental data.

Won, Junho↗
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modelling Tritium Production and Release at High-Energy Proton Accelerators

Tritium is a well-known byproduct of particle accelerator operations. To keep levels of tritium below regulatory limits, tritium production is actively monitored and managed at Fermilab. We plan to study tritium production in the targets, beamline components, and shielding elements of the Fermilab facilities such as NuMI, BNB, and MI-65. To facilitate the analysis, we construct a simple model and use three Monte Carlo radiation codes, FLUKA, MARS, and PHITS, to estimate the amount of tritium produced in these facilities. The analysis could also serve as an intercomparison between these code results related to tritium production. To assess the actual amounts of tritium that would be released from various materials, we employ a semi-empirical diffusion model. The results of this analysis are compared to experimental data whenever possible. This approach also helps to optimize proposed target materials with respect to the tritium production and release.

Georgobiani, Dali [Fermilab]↗
NASA NTRS1981

Double Compton scatter telescope calibration

Calibration techniques for a medium energy gamma ray telescope are described. Gain calibration using Compton edge spectra involves comparisons of pulse height spectra with spectra simulated by a Monte Carlo computer code which includes Compton scattering and pair production, plural scattering and variable energy resolution, and cell size. The telescope considered comprises 56 cells of liquid scintillator in four size groups, with a total liquid volume of 325 l; each cell has its own photomultiplier tube. Energy and angular resolutions and the PMT gain calibration procedure are verified with double scatter data for monoenergetic gamma rays at a known location. Detection probabilities for any cell combination in the two telescope arrays are calculated per steradian as a function of the scattering for a number of different energies with a Van de Graaff accelerator.

Dayton, B.↗
NASA NTRS1985

Monte Carlo calibration of the SMM gamma ray spectrometer for high energy gamma rays and neutrons

The Gamma Ray Spectrometer (GRS) on the Solar Maximum Mission spacecraft was primarily designed and calibrated for nuclear gamma ray line measurements, but also has a high energy mode which allows the detection of gamma rays at energies above 10 MeV and solar neutrons above 20 MeV. The GRS response has been extrapolated until now for high energy gamma rays from an early design study employing Monte Carlo calculations. The response to 50 to 600 MeV solar neutrons was estimated from a simple model which did not consider secondary charged particles escaping into the veto shields. In view of numerous detections by the GRS of solar flares emitting high energy gamma rays, including at least two emitting directly detectable neutrons, the calibration of the high energy mode in the flight model has been recalculated by the use of more sophisticated Monte Carlo computer codes. New results presented show that the GRS response to gamma rays above 20 MeV and to neutrons above 100 MeV is significantly lower than the earlier estimates.

Cooper, J. F.↗
NASA NTRS1985

The polar ionosphere as a source of the storm time ring current

The transport of ions from the polar ionosphere to the inner magnetosphere during storm time conditions has been computed using a Monte Carlo diffusion code. The effect of the electrostatic turbulence assumed to be present during the substorm expansion phase was simulated by a process that accelerated the ions stochastically perpendicular to the magnetic field with a diffusion coefficient proportional to the rate of energization of the ions by the induced electric field. This diffusion process was continued as the ions were convected from the plasma sheet boundary layer to the double-spiral injection boundary. Inward of the injection boundary the ions were convected adiabatically. By using as input an O(+) flux of 2.8 x 10 to the 8th per sq cm per s (w greater than 10 eV) and an H(+) flux of 5.5 x 10 to the 8th per sq cm per s (w greater than .63 eV) the computed distribution functions of the ions in the ring current were found to be in good agreement, over a wide range in L (4-8), with measurements made with the ISEE 1 satellite during a storm. This O(+) flux and a large part of the H(+) flux appear to be consistent with the DE-1 and DE-2 satellite measurements of the polar ionospheric outflow during disturbed times.

Cladis, J. B.↗
NASA NTRS1985

Transport of ions in presence of induced electric field and electrostatic turbulence - Source of ions injected into ring current

The transport of ions from the polar ionosphere to the inner magnetosphere during stormtime conditions has been computed using a Monte Carlo diffusion code. The effect of the electrostatic turbulence assumed to be present during the substorm expansion phase was simulated by a process that accelerated the ions stochastically perpendicular to the magnetic field with a diffusion coefficient proportional to the energization rate of the ions by the induced electric field. This diffusion process was continued as the ions were convected from the plasma sheet boundary layer to the double-spiral injection boundary. Inward of the injection boundary, the ions were convected adiabatically. By using as input an O(+) flux of 2.8 x 10 to the 8th per sq cm per s (w greater than 10 eV) and an H(+) flux of 5.5 x 10 to the 8th per sq cm per s (w greater than 0.63 eV), the computed distribution functions of the ions in the ring current were found to be in good agreement, over a wide range in L (4 to 8), with measurements made with the ISEE-1 satellite during a storm. This O(+) flux and a large part of the H(+) flux are consistent with the DE satellite measurements of the polar ionospheric outflow during disturbed times.

Cladis, J. B.↗
NASA NTRS1992

X-ray scattering and fluorescence in the wind of a massive X-ray binary

Spectral data from the binary X-ray pulsar 4 U 0900 - 40 obtained with the Ginga satellite and covering a range of orbital phases are presented and interpreted using simulated spectra created with a Monte Carlo scattering code. It is found that scattering and fluorescence in a simple spherically symmetric distribution of matter surrounding the companion star can reproduce the eclipse spectrum, as well as the 'soft X-ray-excess' observed during egress and other orbital phases. Reasonably secure values are found for a number of the parameters that characterize the density profile of the stellar wind and atmosphere of the companion star. The Fe abundance is within a factor of about 1.3 of that in the solar neighborhood. It is shown that ionization zones are not critical to understanding the orbital-phase-resolved spectra in this source. It is also found that the contribution by scattering from interstellar dust grains to the observed spectra during eclipse is negligible, while that from diffuse emission from the 'Galactic ridge' is significant.

Lewis, W.↗
NASA NTRS1994

Rates for neutron-capture reactions on tungsten isotopes in iron meteorites

High-precision W isotopic analyses by Harper and Jacobsen indicate the W-182/W-183 ratio in the Toluca iron meteorite is shifted by -(3.0 +/- 0.9) x 10(exp -4) relative to a terrestrial standard. Possible causes of this shift are neutron-capture reactions on W during Toluca's approximately 600-Ma exposure to cosmic ray particles or radiogenic growth of W-182 from 9-Ma Hf-182 in the silicate portion of the Earth after removal of W to the Earth's core. Calculations for the rates of neutron-capture reactions on W isotopes were done to study the first possibility. The LAHET Code System (LCS) which consists of the Los Alamos High Energy Transport (LAHET) code and the Monte Carlo N-Particle(MCNP) transport code was used to numerically simulate the irradiation of the Toluca iron meteorite by galactic-cosmic-ray (GCR) particles and to calculate the rates of W(n, gamma) reactions. Toluca was modeled as a 3.9-m-radius sphere with the composition of a typical IA iron meteorite. The incident GCR protons and their interactions were modeled with LAHET, which also handled the interactions of neutrons with energies above 20 MeV. The rates for the capture of neutrons by W-182, W-183, and W-186 were calculated using the detailed library of (n, gamma) cross sections in MCNP. For this study of the possible effect of W(n, gamma) reactions on W isotope systematics, we consider the peak rates. The calculated maximum change in the normalized W-182/W-183 ratio due to neutron-capture reactions cannot account for more than 25% of the mass 182 deficit observed in Toluca W.

Masarik, J.↗
NASA NTRS1994

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.↗
NASA NTRS1995

The detector response matrices of the burst and transient source experiment (BATSE) on the Compton Gamma Ray Observatory

The detector response matrices for the Burst And Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) are described, including their creation and operation in data analysis. These response matrices are a detailed abstract representation of the gamma-ray detectors' operating characteristics that are needed for data analysis. They are constructed from an extensive set of calibration data coupled with a complex geometry electromagnetic cascade Monte Carlo simulation code. The calibration tests and simulation algorithm optimization are described. The characteristics of the BATSE detectors in the spacecraft environment are also described.

Pendleton, Geoffrey N.↗
NASA NTRS1997

Development of a Space Radiation Monte Carlo Computer Simulation

The ultimate purpose of this effort is to undertake the development of a computer simulation of the radiation environment encountered in spacecraft which is based upon the Monte Carlo technique. The current plan is to adapt and modify a Monte Carlo calculation code known as FLUKA, which is presently used in high energy and heavy ion physics, to simulate the radiation environment present in spacecraft during missions. The initial effort would be directed towards modeling the MIR and Space Shuttle environments, but the long range goal is to develop a program for the accurate prediction of the radiation environment likely to be encountered on future planned endeavors such as the Space Station, a Lunar Return Mission, or a Mars Mission. The longer the mission, especially those which will not have the shielding protection of the earth's magnetic field, the more critical the radiation threat will be. The ultimate goal of this research is to produce a code that will be useful to mission planners and engineers who need to have detailed projections of radiation exposures at specified locations within the spacecraft and for either specific times during the mission or integrated over the entire mission. In concert with the development of the simulation, it is desired to integrate it with a state-of-the-art interactive 3-D graphics-capable analysis package known as ROOT, to allow easy investigation and visualization of the results. The efforts reported on here include the initial development of the program and the demonstration of the efficacy of the technique through a model simulation of the MIR environment. This information was used to write a proposal to obtain follow-on permanent funding for this project.

Pinsky, Lawrence S.↗
NASA NTRS1998

Aerothermodynamics of the Mars Global Surveyor Spacecraft

The aerothermodynamics characteristics of the Mars Global Surveyor spacecraft are investigated and reported. These results have been used by the Mars Global Surveyor mission planners to design the aerobraking phase of the mission. Analytical and Direct Simulation Monte Carlo computer codes were used with a detailed, three dimensional model of the spacecraft to evaluate spacecraft aerobraking characteristics for flight in free molecular and transitional flow regimes. The spacecraft is found to be aerodynamically stable in aerobraking and planned contingency configurations. Aerodynamic forces, moments, and heating are found to be highly dependent on atmospheric density. Accommodation coefficient. is seen to strongly influence drag coefficient. Transitional flow effects are found to reduce overall solar panel heating. Attitude control thruster plumes are shown to interact with the freestream, diminishing the effectiveness of the attitude control system and even leading to thrust reversal. These plume-freestream interaction effects are found to be highly dependent on freestream density.

Shane, Russell W.↗
NASA NTRS1998

Estimation of 3-D Cloud Effects on TOMS Satellite Retrieval of Surface UV Irradiance

To improve surface UV irradiance retrieval from the Total Ozone Mapping Spectrometer (TOMS) we simulate errors of the TOMS cloud correction algorithm for summertime broken cloud conditions. Cloud scenes (50 km by 50 km) are modeled by a normal random (Gaussian) field with a fixed lower boundary and conservative scattering. The model relates stochastic field characteristics with the cloud amount, mean cloud diameter and aspect ratio. Clouds are embedded into Rayleigh atmosphere with standard ozone profile. Radiative transfer calculations of the radiance at the top of the atmosphere and irradiance at the surface were performed using 3-D Monte Carlo (MC) code. The results are averaged over the satellite field of view on the surface (50 km by 50 km) and compared with TOMS predicted surface irradiance for the same scene reflectance. The TOMS algorithm assumes horizontally homogeneous Cl-type cloud between 3 km and 5.5 km. The effective optical depth is determined by fitting observed (MC) radiance at 380 nm. Having the same radiance at the satellite the homogeneous and broken cloud models predict different average irradiances at the surface. This is due to the differences in Bidirectional Reflection Distribution Function (BRDF) for homogeneous and broken cloud scenes with the same hemispherical albedo. For typical TOMS observational geometry at mid-latitudes the simulated single pixels errors may be as large as +/- 20%. Qualitatively these errors are due to the dominance of the non-horizontal cloud surfaces, which are not accounted for in the homogeneous cloud model. However, due to high variability of the real cloud shapes and types it is unclear how these single pixel errors would affect TOMS time-integrated UV exposure over extended periods (weeks to months) for different regions.

Krotkov, Nickolay A.↗