Search NASA⌕ Search

SEARCH · Search NASA

Results for “monte carlo”

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 325 records · Page 18

A Monte Carlo Sensitivity Analysis of CF2 and CF Radical Densities in a c-C4F8 Plasma

A Monte Carlo sensitivity analysis is used to build a plasma chemistry model for octacyclofluorobutane (c-C4F8) which is commonly used in dielectric etch. Experimental data are used both quantitatively and quantitatively to analyze the gas phase and gas surface reactions for neutral radical chemistry. The sensitivity data of the resulting model identifies a few critical gas phase and surface aided reactions that account for most of the uncertainty in the CF2 and CF radical densities. Electron impact dissociation of small radicals (CF2 and CF) and their surface recombination reactions are found to be the rate-limiting steps in the neutral radical chemistry. The relative rates for these electron impact dissociation and surface recombination reactions are also suggested. The resulting mechanism is able to explain the measurements of CF2 and CF densities available in the literature and also their hollow spatial density profiles.

Bose, Deepak↗

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simuulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rarefaction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the after body sting are emphasizes in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗

Monte Carlo Modeling of Sodium in Mercury's Exosphere During the First Two MESSENGER Flybys

We present a Monte Carlo model of the distribution of neutral sodium in Mercury's exosphere and tail using data from the Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft during the first two flybys of the planet in January and September 2008. We show that the dominant source mechanism for ejecting sodium from the surface is photon-stimulated desorption (PSD) and that the desorption rate is limited by the diffusion rate of sodium from the interior of grains in the regolith to the topmost few monolayers where PSD is effective. In the absence of ion precipitation, we find that the sodium source rate is limited to approximately 10(exp 6) - 10(exp 7) per square centimeter per second, depending on the sticking efficiency of exospheric sodium that returns to the surface. The diffusion rate must be at least a factor of 5 higher in regions of ion precipitation to explain the MASCS observations during the second MESSENGER f1yby. We estimate that impact vaporization of micrometeoroids may provide up to 15% of the total sodium source rate in the regions observed. Although sputtering by precipitating ions was found not to be a significant source of sodium during the MESSENGER flybys, ion precipitation is responsible for increasing the source rate at high latitudes through ion-enhanced diffusion.

Burger, Matthew H.↗

Monte Carlo Simulation of a Knudsen Effusion Mass Spectrometer Sampling System

Knudsen flow is easily simulated with a Monte Carlo method. In this study we develop a Visual Basic for Excel (VBA) code to simulate the molecular beam from a vaporizing solid. The system at NASA Glenn uses the restricted collimation method of Chatillon and colleagues, which consists of two apertures between the effusion cell and the ionizer. The diameter of the first aperture is smaller than the diameter of the effusion cell orifice, so the ionizer effectively only sees inside the effusion cell. The code is able to calculate the transmission coefficient through the cell orifice, through the cell orifice and the first aperture, and through the cell orifice and first and second aperatures. Calculated transmission coefficients through the cell orifice are compared to tabulated values to validate the code. Then transmission coefficients are calculated through the cell orifice and both apertures to the ionizer. This allows the geometry (aperture spacing and diameters) of the sampling system the sampling system to be optimized. Calculated transmission factors are also compared to literature values calculated via an analytic method.

Radke, Michael J.↗

Monte Carlo predictions of DNA fragment-size distributions for large sizes after HZE particle irradiation

DSBs (double-strand breaks) produced by densely ionizing space radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. DSB clustering at large scales, from >100 Mbp down to approximately 2 kbp, is modeled using a Monte-Carlo algorithm. A random-walk model of chromatin is combined with a track model, that predicts the radial distribution of energy from an ion, and the RLC (randomly-located-clusters) formalism, in software called DNAbreak. This model generalizes the random-breakage model, whose broken-stick fragment-size distribution is applicable to low-LET radiation. DSB induction due to track interaction with the DNA volume depends on the radiation quality parameter Q. This dose-independent parameter depends only weakly on LET. Multi-track, high-dose effects depend on the cluster intensity parameter lambda, proportional to fluence as defined by the RLC formalism. After lambda is determined by a numerical experiment, the model reduces to one adjustable parameter Q. The best numerical fits to the experimental data, determining Q, are obtained. The knowledge of lambda and Q allows us to give biophysically based extrapolations of high-dose DNA fragment-size data to low doses or to high LETs.

NASA Center JSC↗

Three-dimensional random earth atmospheres for Monte Carlo trajectory analyses

A set of four computer tapes containing random three dimensional Earth atmospheres is available for Monte Carlo trajectory analyses. The tapes contain sufficient atmospheric tables to allow replications of any trajectory below an altitude of 99 km. The atmospheres were provided by an empirical model designed to generate random atmospheres whose distributions match those in a data base of sounding rocket measurements. A readily implementable means of linking the tapes to any existing trajectory simulation computer program is described involving the addition of three subroutines which are listed in an appendix.

Campbell, J. W.↗

Result of Monte-Carlo simulation of electron-photon cascades in lead and layers of lead-scintillator

Results of Monte-Carlo simulation of electromagnetic cascade development in lead and lead-scintillator sandwiches are analyzed. It is demonstrated that the structure function for core approximation is not applicable in the case in which the primary energy is higher than 100 GeV. The simulation data has shown that introducing an inhomogeneous chamber structure results in subsequent reduction of secondary particles.

Wasilewski, A.↗

Regularized Monte-Carlo calculations of the optical polarization induced by scattering in X-ray binaries

The results of regularized Monte Carlo calculations of the optical polarization Stokes parameters expected as a result of scattering from gas streams or accretion wakes in close binary systems whose properties resemble those of X-ray binaries are presented. The physical processes occurring in the system are represented by a model which follows a selected sample of individual photon paths from each emitting element on the surface of the primary to the photon's detection at earth. The result represents the variations in the intensity from each individual emitting element caused by the effects of limb darkening, gravity darkening, and variation in the projected area of the emitting element as seen from the scattering element.

Dolan, J. F.↗

Vectorization of a Monte Carlo simulation scheme for nonequilibrium gas dynamics

Significant improvement has been obtained in the numerical performance of a Monte Carlo scheme for the analysis of nonequilibrium gas dynamics through an implementation of the algorithm which takes advantage of vector hardware, as presently demonstrated through application to three different problems. These are (1) a 1D standing-shock wave; (2) the flow of an expanding gas through an axisymmetric nozzle; and (3) the hypersonic flow of Ar gas over a 3D wedge. Problem (3) is illustrative of the greatly increased number of molecules which the simulation may involve, thanks to improved algorithm performance.

Boyd, Iain D.↗

Monte Carlo Tree Search for Integrated Planning, Learning, and Execution in Nondeterministic Python

We present a novel use of Monte Carlo Tree Search (MCTS),adapted to explore a search space produced by the choice points embedded in Python code. The choice points are non-deterministic assignment statements and subroutine calls. We present MCTS extensions required for doing tree search in this context which includes control constructs like hierarchical decomposition (subroutine calls), iterative while loops and conditional statements. We demonstrate how the system works in a simulated rideshare scenario in an urban setting, and present preliminary experiments as a proof of concept.

Automatic planning↗

Monte Carlo mixture model of lifetime cancer incidence risk from radiation exposure on shuttle and international space station

Estimating uncertainty in lifetime cancer risk for human exposure to space radiation is a unique challenge. Conventional risk assessment with low-linear-energy-transfer (LET)-based risk from Japanese atomic bomb survivor studies may be inappropriate for relativistic protons and nuclei in space due to track structure effects. This paper develops a Monte Carlo mixture model (MCMM) for transferring additive, National Institutes of Health multiplicative, and multiplicative excess cancer incidence risks based on Japanese atomic bomb survivor data to determine excess incidence risk for various US astronaut exposure profiles. The MCMM serves as an anchor point for future risk projection methods involving biophysical models of DNA damage from space radiation. Lifetime incidence risks of radiation-induced cancer for the MCMM based on low-LET Japanese data for nonleukemia (all cancers except leukemia) were 2.77 (90% confidence limit, 0.75-11.34) for males exposed to 1 Sv at age 45 and 2.20 (90% confidence limit, 0.59-10.12) for males exposed at age 55. For females, mixture model risks for nonleukemia exposed separately to 1 Sv at ages of 45 and 55 were 2.98 (90% confidence limit, 0.90-11.70) and 2.44 (90% confidence limit, 0.70-10.30), respectively. Risks for high-LET 200 MeV protons (LET=0.45 keV/micrometer), 1 MeV alpha-particles (LET=100 keV/micrometer), and 600 MeV iron particles (LET=180 keV/micrometer) were scored on a per particle basis by determining the particle fluence required for an average of one particle per cell nucleus of area 100 micrometer(2). Lifetime risk per proton was 2.68x10(-2)% (90% confidence limit, 0.79x10(-3)%-0. 514x10(-2)%). For alpha-particles, lifetime risk was 14.2% (90% confidence limit, 2.5%-31.2%). Conversely, lifetime risk per iron particle was 23.7% (90% confidence limit, 4.5%-53.0%). Uncertainty in the DDREF for high-LET particles may be less than that for low-LET radiation because typically there is very little dose-rate dependence. Probability density functions for high-LET radiation quality and dose-rate may be preferable to conventional risk assessment approaches. Nuclear reactions and track structure effects in tissue may not be properly estimated by existing data using in vitro models for estimating RBEs. The method used here is being extended to estimate uncertainty in spacecraft shielding effectiveness in various space radiation environments.

NASA Center LaRC↗

A highly optimized vectorized code for Monte Carlo simulations of SU(3) lattice gauge theories

New methods are introduced for improving the performance of the vectorized Monte Carlo SU(3) lattice gauge theory algorithm using the CDC CYBER 205. Structure, algorithm and programming considerations are discussed. The performance achieved for a 16(4) lattice on a 2-pipe system may be phrased in terms of the link update time or overall MFLOPS rates. For 32-bit arithmetic, it is 36.3 microsecond/link for 8 hits per iteration (40.9 microsecond for 10 hits) or 101.5 MFLOPS.

Barkai, D.↗