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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Monte Carlo analysis of satellite debris footprint dispersion

A comprehensive study is performed to investigate satellite debris impact point dispersion using a combination of Monte Carlo statistical analysis and parametric methods. The Monte Carlo technique accounts for nonlinearities in the entry point dispersion, which is represented by a covariance matrix of position and velocity errors. Because downrange distance of impact is a monotonic function of debris ballistic coefficient, a parametric method is useful for determining dispersion boundaries. The scheme is applied in the present analysis to estimate the Skylab footprint dispersions for a controlled reentry. A significant increase in the footprint dispersion is noticed for satellite breakup above a 200,000-ft altitude. A general discussion of the method used for analysis is presented together with some typical results obtained for the Skylab deboost mission, which was designed before NASA abandoned plans for a Skylab controlled reentry.

Rao, P. P.↗
NASA NTRS1994

Monte Carlo simulations of particle acceleration at oblique shocks

The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.

Baring, Matthew G.↗
NASA NTRS1992

Monte Carlo modeling of atomic oxygen attack of polymers with protective coatings on LDEF

Characterization of the behavior of atomic oxygen interaction with materials on the Long Duration Exposure Facility (LDEF) will assist in understanding the mechanisms involved, and will lead to improved reliability in predicting in-space durability of materials based on ground laboratory testing. A computational simulation of atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of assumed mechanistic behavior of atomic oxygen and results of both ground laboratory and LDEF data, a predictive Monte Carlo model was developed which simulates the oxidation processes that occur on polymers with applied protective coatings that have defects. The use of high atomic oxygen fluence-directed ram LDEF results has enabled mechanistic implications to be made by adjusting Monte Carlo modeling assumptions to match observed results based on scanning electron microscopy. Modeling assumptions, implications, and predictions are presented, along with comparison of observed ground laboratory and LDEF results.

Banks, Bruce A.↗
NASA NTRS1989

Monte Carlo simulation of an expanding gas

By application of simple computer graphics techniques, the statistical performance of two Monte Carlo methods used in the simulation of rarefied gas flows are assessed. Specifically, two direct simulation Monte Carlo (DSMC) methods developed by Bird and Nanbu are considered. The graphics techniques are found to be of great benefit in the reduction and interpretation of the large volume of data generated, thus enabling important conclusions to be drawn about the simulation results. Hence, it is discovered that the method of Nanbu suffers from increased statistical fluctuations, thereby prohibiting its use in the solution of practical problems.

Boyd, Iain D.↗
NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗
NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Computational modeling of the erosion of polymers caused by atomic oxygen in low Earth orbit (LEO) is useful for determining areas of concern for spacecraft environment durability. Successful modeling requires that the characteristics of the environment such as atomic oxygen energy distribution, flux, and angular distribution be properly represented in the model. Thus whether the atomic oxygen is arriving normal to or inclined to a surface and whether it arrives in a consistent direction or is sweeping across the surface such as in the case of polymeric solar array blankets is important to determine durability. When atomic oxygen impacts a polymer surface it can react removing a certain volume per incident atom (called the erosion yield), recombine, or be ejected as an active oxygen atom to potentially either react with other polymer atoms or exit into space. Scattered atoms can also have a lower energy as a result of partial or total thermal accommodation. Many solutions to polymer durability in LEO involve protective thin films of metal oxides such as SiO2 to prevent atomic oxygen erosion. Such protective films also have their own interaction characteristics. A Monte Carlo computational model has been developed which takes into account the various types of atomic oxygen arrival and how it reacts with a representative polymer (polyimide Kapton H) and how it reacts at defect sites in an oxide protective coating, such as SiO2 on that polymer. Although this model was initially intended to determine atomic oxygen erosion behavior at defect sites for the International Space Station solar arrays, it has been used to predict atomic oxygen erosion or oxidation behavior on many other spacecraft components including erosion of polymeric joints, durability of solar array blanket box covers, and scattering of atomic oxygen into telescopes and microwave cavities where oxidation of critical component surfaces can take place. The computational model is a two dimensional model which has the capability to tune the interactions of how the atomic oxygen reacts, scatters, or recombines on polymer or nonreactive surfaces. In addition to the specification of atomic oxygen arrival details, a total of 15 atomic oxygen interaction parameters have been identified as necessary to properly simulate observed interactions and resulting polymer erosion that have been observed in LEO. The tuning of the Monte Carlo model has been accomplished by adjusting interaction parameters so the erosion patterns produced by the model match those from several actual LEO space experiments. Surface texturing in LEO can also be predicted by the model. Such comparison of space tests with ground laboratory experiments have enabled confidence in ground laboratory lifetime prediction of protected polymers. Results of Monte Carlo tuning, examples of surface texturing and undercutting erosion prediction, and several examples of how the model can be used to predict other LEO and Mars orbital space results are presented.

Monte Carlo Modeling↗
NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Computational modeling of the erosion of polymers caused by atomic oxygen in low Earth orbit (LEO) is useful for determining areas of concern for spacecraft environment durability. Successful modeling requires that the characteristics of the environment such as atomic oxygen energy distribution, flux, and angular distribution be properly represented in the model. Thus whether the atomic oxygen is arriving normal to or inclined to a surface and whether it arrives in a consistent direction or is sweeping across the surface such as in the case of polymeric solar array blankets is important to determine durability. When atomic oxygen impacts a polymer surface it can react removing a certain volume per incident atom (called the erosion yield), recombine, or be ejected as an active oxygen atom to potentially either react with other polymer atoms or exit into space. Scattered atoms can also have a lower energy as a result of partial or total thermal accommodation. Many solutions to polymer durability in LEO involve protective thin films of metal oxides such as SiO2 to prevent atomic oxygen erosion. Such protective films also have their own interaction characteristics. A Monte Carlo computational model has been developed which takes into account the various types of atomic oxygen arrival and how it reacts with a representative polymer (polyimide Kapton H) and how it reacts at defect sites in an oxide protective coating, such as SiO2 on that polymer. Although this model was initially intended to determine atomic oxygen erosion behavior at defect sites for the International Space Station solar arrays, it has been used to predict atomic oxygen erosion or oxidation behavior on many other spacecraft components including erosion of polymeric joints, durability of solar array blanket box covers, and scattering of atomic oxygen into telescopes and microwave cavities where oxidation of critical component surfaces can take place. The computational model is a two dimensional model which has the capability to tune the interactions of how the atomic oxygen reacts, scatters, or recombines on polymer or nonreactive surfaces. In addition to the specification of atomic oxygen arrival details, a total of 15 atomic oxygen interaction parameters have been identified as necessary to properly simulate observed interactions and resulting polymer erosion that have been observed in LEO. The tuning of the Monte Carlo model has been accomplished by adjusting interaction parameters so the erosion patterns produced by the model match those from several actual LEO space experiments. Surface texturing in LEO can also be predicted by the model. Such comparison of space tests with ground laboratory experiments have enabled confidence in ground laboratory lifetime prediction of protected polymers. Results of Monte Carlo tuning, examples of surface texturing and undercutting erosion prediction, and several examples of how the model can be used to predict other LEO and Mars orbital space results are presented.

Atomic Oxygen↗
NASA NTRS1994

A Monte Carlo study of photochemical model uncertainties affecting the calculation of subsonic perturbations in the upper troposphere

Monte Carlo methods are frequently applied to the evaluation of uncertainties in models with multiple inputs that themselves have associated imprecisions. In the case of photochemical models used to evaluate changes in O3 or OH, inputs analyzed include rate coefficients measured in the laboratory, and chemical and physical constituents measured in the atmosphere. The Monte Carlo method was used with the 1-dimensional GSFC tropospheric photochemical model to examine uncertainty propagation to calculation of NO(x) and the major oxidants (O3, OH) in the upper troposphere. In all cases chemical kinetics inputs are varied and the NO(x) perturbation of a subsonic fleet is simulated. A series of model runs is used to explore sensitivities of model-computed parameters to other parameters: heterogeneous processes, the uncertainty in upper tropospheric H2O vapor measurements, aircraft emissions at different latitudes.

Thompson, Anne M.↗
NASA NTRS1978

Monte Carlo modeling of exospheric bodies - Mercury

In order to study the interaction with the surface, a Monte Carlo program is developed to determine the distribution with altitude as well as the global distribution of density at the surface in a single operation. The analysis presented shows that the appropriate source distribution should be Maxwell-Boltzmann flux if the particles in the distribution are to be treated as components of flux. Monte Carlo calculations with a Maxwell-Boltzmann flux source are compared with Mariner 10 UV spectrometer data. Results indicate that the presently operating models are not capable of fitting the observed Mercury exosphere. It is suggested that an atmosphere calculated with a barometric source distribution is suitable for more realistic future exospheric models.

Smith, G. R.↗
NASA NTRS1973

Monte Carlo applications to acoustical field solutions

The Monte Carlo technique is proposed for the determination of the acoustical pressure-time history at chosen points in a partial enclosure, the central idea of this technique being the tracing of acoustical rays. A statistical model is formulated and an algorithm for pressure is developed, the conformity of which is examined by two approaches and is shown to give the known results. The concepts that are developed are applied to the determination of the transient field due to a sound source in a homogeneous medium in a rectangular enclosure with perfect reflecting walls, and the results are compared with those presented by Mintzer based on the Laplace transform approach, as well as with a normal mode solution.

Haviland, J. K.↗
NASA NTRS1975

Catastrophic rupture of lunar rocks - A Monte Carlo simulation

A computer model based on Monte Carlo techniques was developed to simulate the destruction of lunar rocks by 'catastrophic rupture' due to meteoroid impact. Energies necessary to accomplish catastrophic rupture were derived from laboratory experiments. A crater-production rate derived from lunar rocks was utilized to calculate absolute time scales. Calculated median survival times for crystalline lunar rocks are 1.9, 4.6, 10.3, and 22 m.y. for rock masses of 10, 100, 1000, and 10,000 g, respectively. Corresponding times of 6, 14.5, 32, and 68 million years are required before the probability of destruction reaches 0.99. These results are consistent with absolute exposure ages measured on returned rocks. Some results also substantiate previous conclusions that the catastrophic-rupture process is significantly more effective in obliterating lunar rocks than mass wasting by single-particle abrasion. The view is also corroborated that most rocks presently on the lunar surface either are exhumed from the regolith or are fragments of much larger boulders rather than primary ejecta excavated from pristine bedrock.

Hoerz, F.↗
NASA NTRS1985

Monte Carlo Modeling of Lunar Megaregolith Development

A Monte Carlo model of lunar megaregolith development is proposed. Minimum megaregolith depths are obtained in the preliminary model. The model impacts a gridded, initially flat surface area randomly, with points not constrained to grid intersections. Information is recorded at grid points on evolving topography, deepest depth impacted to, and number of hits at that location. Two types of lunar surfaces ar simulated, maria and highlands. Crater production functions are after run to densities that are multiples of observed lunar crater densities for craters of 8 km in diameter and larger. all craters are considered to have formed in a gravity scaled regime. Edge effect problems are corrected in this model by allowing craters forming outside the grid area or along its edges to effect the gridded area.

Cashore, J.↗
NASA NTRS1986

Monte Carlo study of the atmospheric spread function

Monte Carlo radiative transfer simulations are used to study the atmospheric spread function appropriate to satellite-based sensing of the earth's surface. The parameters which are explored include the nadir angle of view, the size distribution of the atmospheric aerosol, and the aerosol vertical profile.

Pearce, W. A.↗
NASA NTRS1988

Monte-Carlo histories of refractory interstellar dust

Monte-carlo histories of 6 x 10 to the 6th individual dust particles injected uniformly from stars into the interstellar medium during a 6 x 10 to the 9th year history are calculated. The particles are given a two-phase internal structure of successive thermal condensates, and are distributed in initial radius as 1/a-cubed over the value of a between 0.01 and 0.1 micron. The evolution of this system illustrates the distinction between several different lifetimes for interstellar dust. Most are destroyed, but some grow in size. Several important consequences for interstellar dust are described.

Clayton, D. D.↗