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 343 records · Page 19

Shuttle vertical fin flowfield by the direct simulation Monte Carlo method

The flow properties in a model flowfield, simulating the shuttle vertical fin, determined using the Direct Simulation Monte Carlo method. The case analyzed corresponds to an orbit height of 225 km with the freestream velocity vector orthogonal to the fin surface. Contour plots of the flowfield distributions of density, temperature, velocity and flow angle are presented. The results also include mean molecular collision frequency (which reaches 1/60 sec near the surface), collision frequency density (approaches 7 x 10 to the 18/cu m sec at the surface) and the mean free path (19 m at the surface).

Hueser, J. E.↗

Morphological evolution of growing crystals - A Monte Carlo simulation

The combined effects of nutrient diffusion and surface kinetics on the crystal morphology were investigated using a Monte Carlo model to simulate the evolving morphology of a crystal growing from a two-component gaseous nutrient phase. The model combines nutrient diffusion, based on a modified diffusion-limited aggregation process, with anisotropic surface-attachment kinetics and surface diffusion. A variety of conditions, ranging from kinetic-controlled to diffusion-controlled growth, were examined. Successive transitions from compact faceted (dominant surface kinetics) to open dendritic morphologies (dominant volume diffusion) were obtained.

Xiao, Rong-Fu↗

Nonequilibrium gas flows. I - A detailed validation of Monte Carlo direct simulation for monatomic gases

One-dimensional shock wave properties in helium and argon are predicted using Monte Carlo direct simulation. The collision model is based directly on the interatomic potential taking angular scattering into account. The potential is assumed to be of the Maitland-Smith n(r)-6 form. The detailed validity of the simulation is studied by comparing the predicted macroscopic and miroscopic flow properties in shock waves to a wide range of available data.

Erwin, Daniel A.↗

Direct simulation Monte Carlo with ionization and radiation

Improvements in the modeling of radiation in low density shock waves with Direct Simulation Monte Carlo (DSMC) are the subject of this study. The physical processes which determine the amount of radiation in a shock wave were investigated and the way in which they were modeled with DSMC was evaluated. Three physical processes were identified for which an improvement in the modeling technique could result in improved radiation predictions. New physical modeling schemes are introduced in this report for the three processes. First a method for determining the electric field and its effect on the flow is introduced. Second, a two step reaction process for electron impact ionization reactions is evaluated. Finally, a new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. Each new scheme attempts to move the DSMC method toward more accurate physics or more reliance on experimental data. The new schemes are all compared to the current modeling techniques and the differences in the results are evaluated. In all cases the results agree with the available data as well as, or better than the results from earlier schemes.

Carlson, Ann B.↗

Procedure for Adapting Direct Simulation Monte Carlo Meshes

A technique is presented for adapting computational meshes used in the G2 version of the direct simulation Monte Carlo method. The physical ideas underlying the technique are discussed, and adaptation formulas are developed for use on solutions generated from an initial mesh. The effect of statistical scatter on adaptation is addressed, and results demonstrate the ability of this technique to achieve more accurate results without increasing necessary computational resources.

Woronowicz, Michael S.↗

Monte Carlo modeling of the net effects of coma scattering and thermal reradiation on the energy input to cometary nucleus

A Monte Carlo simulation method is presented that can, to an accuracy of a few percent, calculate the effects of a dusty coma on the total energy input to the cometary nucleus. This method treats nonconservative nonisotropic scattering, as well as the reflection from the nucleus surface. Results are presented as a function of the optical thickness of the dust column in the sun-comet axis. The total energy input to the nucleus appears to be only weakly dependent on the opacity of the coma, the radial distribution of the dust, or the details of the extinction processes.

Salo, H.↗

Monte Carlo simulation of reentry plasmas

Attention is given to the treatment of ionization and plasma effects in the direct simulation Monte Carlo method. The requirements for accurate modeling of reentry plasmas are discussed along with the difficulties these requirements present. The current method for modeling such plasmas is reviewed and an alternative method is presented. Both methods are applied to the flow of a 10 km/s shock wave in air at 0.1 torr; a flowfield directly relevant to the projected aeroassisted orbital transfer vehicle. The results are compared and the differences between the methods are discussed.

Carlson, Ann B.↗

A Monte Carlo simulation technique for low-altitude, wind-shear turbulence

A case is made for including anisotropy in a Monte Carlo flight simulation scheme of low-altitude wind-shear turbulence by means of power spectral density. This study attempts to eliminate all flight simulation-induced deficiencies in the basic turbulence model. A full-scale low-altitude wind-shear turbulence simulation scheme is proposed with particular emphasis on low cost and practicality for near-ground flight. The power spectral density statistic is used to highlight the need for realistic estimates of energy transfer associated with low-altitude wind-shear turbulence. The simulation of a particular anisotropic turbulence model is shown to be a relatively simple extension from that of traditional isotropic (Dryden) turbulence.

Bowles, Roland L.↗

Monte Carlo modeling of the resurfacing of Venus

We have developed a three-dimensional model of venusian resurfacing that employs Monte Carlo simulations of both impact cratering and volcanism. The model simulates the production of craters on Venus by using the observed mass distributions of Earth- and Venus-crossing asteroids and comets. Lava flows are modeled by an energy minimization technique to simulate the effects of local topography on the shape and extent of flows. The model is run under a wide range of assumptions regarding the scale and time evolution of volcanism on Venus. Regions of the parameter space that result in impact crater distributions and modifications that are currently observed will be explored to place limits on the possible volcanic resurfacing history of Venus.

Bullock, M. A.↗

Direct simulation Monte Carlo prediction of on-orbit contaminant deposit levels for HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flow field and surface conditions and geometric orientations for the satellite in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. A detailed description of the adaptation of this solution method to the study of the satellite's environment is also presented. Results pertaining to the satellite's environment are presented regarding contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface, along with data related to code performance. Using procedures developed in standard contamination analyses, along with many worst-case assumptions, the cumulative upper-limit level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated at about 13,350 A.

Woronowicz, Michael S.↗

SUPREM-DSMC: A New Scalable, Parallel, Reacting, Multidimensional Direct Simulation Monte Carlo Flow Code

An AFRL/NRL team has recently been selected to develop a scalable, parallel, reacting, multidimensional (SUPREM) Direct Simulation Monte Carlo (DSMC) code for the DoD user community under the High Performance Computing Modernization Office (HPCMO) Common High Performance Computing Software Support Initiative (CHSSI). This paper will introduce the JANNAF Exhaust Plume community to this three-year development effort and present the overall goals, schedule, and current status of this new code.

Campbell, David↗

CLUST - EVAP Monte Carlo Simulation Applications for Determining Effective Energy Deposition in Silicon by High Energy Protons

The CLUST-EVAP is a Monte Carlo simulation of the interaction of high energy (25 - 400 MeV) protons with silicon nuclei. The initial nuclear cascade stage is modeled using the CLUST model developed by Indiana University over 30 years ago. The second stage, in which the excited nucleus evaporates particles in random directions, is modeled according to the evaporation algorithm provided by H. H. K. Tang of IBM. Using the CLUST-EVAP code to model fragment produ6tion and the Vavilov-Landau theory to model fluctuations in direct ionization in thin silicon layers, we have predicted energy deposition in silicon components for various geometrical configurations. We have compared actual measurements with model predictions for geometry's such as single, thin silicon particle detectors, telescopic particle detectors flown in space to measure the environment, and thin sensitive volumes of modern micro-electronic components. We have recently compared the model predictions with actual measurements made by the DOSTEL spectrometer flown in the Shuttle payload bay on STS-84. The model faithfully reproduces the features and aids in interpretation of flight results of this instrument. We have also applied the CLUST-EVAP model to determine energy deposition in the thin sensitive volumes of modern micro-electronic components. We have accessed the ability of high energy (200 MeV) protons to induce latch-up in certain devices that are known to latch up in heavy ion environments. However, some devices are not nearly as susceptible to proton induced latch-up as expected according to their measured heavy ion latch-up cross sections. The discrepancy is believed to be caused by the limited range of the proton-silicon interaction fragments. The CLUST-EV AP model was used to determine a distribution of these fragments and their range and this is compared to knowledge of the ranges required based on the known device structure. This information is especially useful in accessing the risk to on-orbit perfonnance in a heavy ion environment based on testing performed with only protons.

ONeill, Pat M.↗

Refining Lucy Mission Delta-V During Spacecraft Design Using Trajectory Optimization Within High-Fidelity Monte Carlo Maneuver Analysis

Recent advances linking medium-fidelity trajectory optimization and high-fidelity trajectory propagation/maneuver design software with Monte Carlo maneuver analysis and parallel processing enabled realistic statistical delta-V estimation well before launch. Completing this high-confidence, refined statistical maneuver analysis early enabled release of excess delta-V margin for increased dry mass margin for the Lucy Jupiter Trojan flyby mission. By 3.3 years before launch, 16 of 34 TCMs had 1000 re-optimized trajectory design samples, yielding tens of m/s lower 99%-probability delta-V versus targeting maneuvers to one optimal trajectory. One year later, 1000 re-optimized samples of all deterministic maneuvers and subsequent flybys further lowered estimated delta-V.

Trojan↗

Monte Carlo Analysis of the Commissioning Phase Maneuvers of the Soil Moisture Active Passive (SMAP) Mission

The Soil Moisture Active Passive (SMAP) mission will perform soil moisture content and freeze/thaw state observations from a low-Earth orbit. The observatory is scheduled to launch in October 2014 and will perform observations from a near-polar, frozen, and sun-synchronous Science Orbit for a 3-year data collection mission. At launch, the observatory is delivered to an Injection Orbit that is biased below the Science Orbit; the spacecraft will maneuver to the Science Orbit during the mission Commissioning Phase. The delta V needed to maneuver from the Injection Orbit to the Science Orbit is computed statistically via a Monte Carlo simulation; the 99th percentile delta V (delta V99) is carried as a line item in the mission delta V budget. This paper details the simulation and analysis performed to compute this figure and the delta V99 computed per current mission parameters.

Monte Carlo simulations↗

Monte Carlo transport model comparison with 1A GeV accelerated iron experiment: heavy-ion shielding evaluation of NASA space flight-crew foodstuff

Deep-space manned flight as a reality depends on a viable solution to the radiation problem. Both acute and chronic radiation health threats are known to exist, with solar particle events as an example of the former and galactic cosmic rays (GCR) of the latter. In this experiment Iron ions of 1A GeV are used to simulate GCR and to determine the secondary radiation field created as the GCR-like particles interact with a thick target. A NASA prepared food pantry locker was subjected to the iron beam and the secondary fluence recorded. A modified version of the Monte Carlo heavy ion transport code developed by Zeitlin at LBNL is compared with experimental fluence. The foodstuff is modeled as mixed nuts as defined by the 71st edition of the Chemical Rubber Company (CRC) Handbook of Physics and Chemistry. The results indicate a good agreement between the experimental data and the model. The agreement between model and experiment is determined using a linear fit to ordered pairs of data. The intercept is forced to zero. The slope fit is 0.825 and the R2 value is 0.429 over the resolved fluence region. The removal of an outlier, Z=14, gives values of 0.888 and 0.705 for slope and R2 respectively. c2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

NASA Discipline Radiation Health↗

A Monte Carlo analysis of the Viking lander dynamics at touchdown

The performance of the Viking lander has been evaluated by using a Monte Carlo simulation, and all results are presented in statistical form. The primary objectives of this analysis were as follows: (1) to determine the three sigma design values of maximum rigid body accelerations and the minimum clearance of the lander body during landing; (2) to determine the probability of an unstable landing; and (3) to determine the probability of the lander body striking a rock. Two configurations were analyzed with the only difference being in the ability of the primary landing gear struts to carry tension loads.

Muraca, R. J.↗