Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 559 records · Page 31

Investigation of the external flow analysis for density measurements at high altitude

A Monte Carlo program was developed for modeling the flow field around the space shuttle in the vicinity of the shuttle upper atmosphere mass spectrometer experiment. The operation of the EXTERNAL code is summarized. Issues associated with geometric modeling of the shuttle nose region and the modeling of intermolecular collisions including rotational energy exchange are discussed as well as a preliminary analysis of vibrational excitation and dissociation effects. The selection of trial runs is described and the parameters used for them is justified. The original version and the modified INTERNAL code for the entrance problem are reviewed. The code listing is included.

Bienkowski, G. K.↗

Studies of radiative transfer in planetary atmospheres

Schloerb and Claussen continued their analysis of the very high quality data set obtained on the 18 centimeter OH line from the Comet P/Halley with the NRAO 43 meter antenna. The high spectral resolution (0.22 km/sec) and high signal-to-noise of the OH spectra make them ideal for the study of kinematics in the coma. Synthetic profiles were initiated for comparison with the data. A vectorial model was developed using the Monte Carlo techniques originated by Combi and Delsemme. Analysis of the millimeter wavelength observations of HCN emission from P/Halley obtained throughout much of the recent apparition were continued using the University of Massachusetts 14 millimeter-wavelength (FCRAO) antenna. A detailed analysis of the HCN lineshpaes was performed over the last six months. The excitation of HCN in the coma was studied to obtain a detailed match to the observed spectra. The passive millimeter wave radiometer was used to probe the physical and chemical nature of comets from spacecraft. Work was continued on an improved theory of radiative transfer for rough and porous surfaces, such as the regoliths of satellites, asteroids, and comets.

Irvine, W. M.↗

Variability analysis in low count rate sources

A method, based on the chi-square statistics, is described for detecting pulselike time variability in low count rate sources observed with photon-counting instruments. This method can be used even in the presence of observational gaps, takes full advantage of the filtering effect due to binning with different bin sizes, and takes into account the arbitrariness introduced by the binning phase. The procedure developed to limit the dependence of the results on the binning phase and ensure statistically correct results is described along with the application of the proposed procedure to a model of a variable source. Monte Carlo simulations are used to show how the method can be used to derive the characteristic variability time scales and that the method is more sensitive than the nonparametric Kolmogorov-Smirnov test in detecting variability to a given confidence level.

Collura, A.↗

The economics of satellite retrieval

The economics of space operations with and without the Space Station have been studied in terms of the financial performance of a typical communications-satellite business venture. A stochastic Monte-Carlo communications-satellite business model is employed which includes factors such as satellite configuration, random and wearout failures, reliability of launch and space operations, stand-down time resulting from failures, and insurance by operation. Financial performance impacts have been evaluated in terms of the magnitude of investment, net present value, and return on investment.

Price, Kent M.↗

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.↗

Lidar multiple scattering: Dependence on atmospheric parameters

Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).

Wiegner, Matthias↗

Improving the chi-squared approximation for bivariate normal tolerance regions

Let X be a two-dimensional random variable distributed according to N2(mu,Sigma) and let bar-X and S be the respective sample mean and covariance matrix calculated from N observations of X. Given a containment probability beta and a level of confidence gamma, we seek a number c, depending only on N, beta, and gamma such that the ellipsoid R = (x: (x - bar-X)'S(exp -1) (x - bar-X) less than or = c) is a tolerance region of content beta and level gamma; i.e., R has probability gamma of containing at least 100 beta percent of the distribution of X. Various approximations for c exist in the literature, but one of the simplest to compute -- a multiple of the ratio of certain chi-squared percentage points -- is badly biased for small N. For the bivariate normal case, most of the bias can be removed by simple adjustment using a factor A which depends on beta and gamma. This paper provides values of A for various beta and gamma so that the simple approximation for c can be made viable for any reasonable sample size. The methodology provides an illustrative example of how a combination of Monte-Carlo simulation and simple regression modelling can be used to improve an existing approximation.

Feiveson, Alan H.↗

Probabilistic micromechanics for metal matrix composites

A probabilistic micromechanics-based nonlinear analysis procedure is developed to predict and quantify the variability in the properties of high temperature metal matrix composites. Monte Carlo simulation is used to model the probabilistic distributions of the constituent level properties including fiber, matrix, and interphase properties, volume and void ratios, strengths, fiber misalignment, and nonlinear empirical parameters. The procedure predicts the resultant ply properties and quantifies their statistical scatter. Graphite copper and Silicon Carbide Titanlum Aluminide (SCS-6 TI15) unidirectional plies are considered to demonstrate the predictive capabilities. The procedure is believed to have a high potential for use in material characterization and selection to precede and assist in experimental studies of new high temperature metal matrix composites.

Engelstad, S. P.↗

Imaging the sun in hard X-rays - Spatial and rotating modulation collimators

Several approaches to imaging hard X-rays emitted from solar flares have been proposed or are planned for the nineties including the spatial modulation collimator (SMC) and the rotating modulation collimator (RMC). A survey of current solar flare theoretical literature indicates the desirability of spatial resolutions down to 1 arcsecond, field of views greater than the full solar disk (i.e., 32 arcminutes), and temporal resolutions down to 1 second. Although the sun typically provides relatively high flux levels, the requirement for 1 second temporal resolution raises the question as to the viability of Fourier telescopes subject to the aforementioned constraints. A basic photon counting, Monte Carlo 'end-to-end' model telescope was employed using the Astronomical Image Processing System (AIPS) for image reconstruction. The resulting solar flare hard X-ray images compared against typical observations indicated that both telescopes show promise for the future.

Campbell, Jonathan W.↗

Powder agglomeration in a microgravity environment

This is the final report for NASA Grant NAG3-755 entitled 'Powder Agglomeration in a Microgravity Environment.' The research program included both two types of numerical models and two types of experiments. The numerical modeling included the use of Monte Carlo type simulations of agglomerate growth including hydrodynamic screening and molecular dynamics type simulations of the rearrangement of particles within an agglomerate under a gravitational field. Experiments included direct observation of the agglomeration of submicron alumina and indirect observation, using small angle light scattering, of the agglomeration of colloidal silica and aluminum monohydroxide. In the former class of experiments, the powders were constrained to move on a two-dimensional surface oriented to minimize the effect of gravity. In the latter, some experiments involved mixture of suspensions containing particles of opposite charge which resulted in agglomeration on a very short time scale relative to settling under gravity.

Cawley, James D.↗

The global resurfacing of Venus

The impact cratering record on Venus is unique among the terrestrial planets. Fully 84% of the craters are in pristine condition, and only 12% are fractured. Remarkably, only 2.5% of the craters and crater-related features are embayed by lava, although intense volcanism and tectonism have affected the entire planet. Furthermore, the spatial and hypsometric distribution of the craters is consistent with a completely random one, including stochastic variations. Monte Carlo simulations of equilibrium resurfacing models result in a minimum of 17 times more embayed craters than observed, or unobserved nonrandom crater distributions for resurfacing areas between 0.03% and 100% of the planet's surface. These models also are not consistent with the number and nonrandom distribution of volcanoes, and the nonrandom distribution of embayed and heavily fractured craters. The constraints imposed by the cratering record strongly indicate that Venus experienced a global resurfacing event about 300 m.y. ago followed by a dramatic reduction of volcanism and tectonism. This global resurfacing event ended abruptly (less than 10 m.y.). The present crater population has accumulated since then and remains largely intact. Thermal history models suggest that similar global resurfacing events probably occured episodically in the past. We show that neither the present level and style of geologic activity nor anything less than global resurfacing could have produced the observed cratering record. The effects of recent geologic activity are much less than those of the earlier global resurfacing event, when the record of all the early heavy bombardment and much of the later light bombardment was erased from the surface by massive volcanism and tectonic activity. Episodic regional resurfacing events that had global effects also occurred on Earth (e.g., the mid-Cretaceous superplume) and probably on Mars. On Mars they may have triggered the catastrophic releases of water that formed the outflow channels.

Strom, Robert G.↗

System Design under Uncertainty: Evolutionary Optimization of the Gravity Probe-B Spacecraft

This paper discusses the application of evolutionary random-search algorithms (Simulated Annealing and Genetic Algorithms) to the problem of spacecraft design under performance uncertainty. Traditionally, spacecraft performance uncertainty has been measured by reliability. Published algorithms for reliability optimization are seldom used in practice because they oversimplify reality. The algorithm developed here uses random-search optimization to allow us to model the problem more realistically. Monte Carlo simulations are used to evaluate the objective function for each trial design solution. These methods have been applied to the Gravity Probe-B (GP-B) spacecraft being developed at Stanford University for launch in 1999, Results of the algorithm developed here for GP-13 are shown, and their implications for design optimization by evolutionary algorithms are discussed.

Pullen, Samuel P.↗

MeV Gamma Ray Detection Algorithms for Stacked Silicon Detectors

By making use of the signature of a gamma ray event as it appears in N = 5 to 20 lithium-drifted silicon detectors and applying smart selection algorithms, gamma rays in the energy range of 1 to 8 MeV can be detected with good efficiency and selectivity. Examples of the types of algorithms used for different energy regions include the simple sum mode, the sum-coincidence mode used in segmented detectors, unique variations on sum-coincidence for an N-dimensional vector event, and a new and extremely useful mode for double escape peak spectroscopy at pair-production energies. The latter algorithm yields a spectrum similar to that of the pair spectrometer, but without the need of the dual external segments for double escape coincidence, and without the large loss in efficiency of double escape events. Background events due to Compton scattering are largely suppressed. Monte Carlo calculations were used to model the gamma ray interactions in the silicon, in order to enable testing of a wide array of different algorithms on the event N-vectors for a large-N stack.

McMurray, Robert E. Jr.↗

Simulations of Ground and Space-Based Oxygen Atom Experiments

A low-earth orbit (LEO) materials erosion scenario and the ground-based experiment designed to simulate it are compared using the direct-simulation Monte Carlo (DSMC) method. The DSMC model provides a detailed description of the interactions between the hyperthermal gas flow and a normally oriented flat plate for each case. We find that while the general characteristics of the LEO exposure are represented in the ground-based experiment, multi-collision effects can potentially alter the impact energy and directionality of the impinging molecules in the ground-based experiment. Multi-collision phenomena also affect downstream flux measurements.

Finchum, A.↗

Space Radiation Cancer Risks and Uncertainties for Mars Missions

Projecting cancer risks from exposure to space radiation is highly uncertain because of the absence of data for humans and because of the limited radiobiology data available for estimating late effects from the high-energy and charge (HZE) ions present in the galactic cosmic rays (GCR). Cancer risk projections involve many biological and physical factors, each of which has a differential range of uncertainty due to the lack of data and knowledge. We discuss an uncertainty assessment within the linear-additivity model using the approach of Monte Carlo sampling from subjective error distributions that represent the lack of knowledge in each factor to quantify the overall uncertainty in risk projections. Calculations are performed using the space radiation environment and transport codes for several Mars mission scenarios. This approach leads to estimates of the uncertainties in cancer risk projections of 400-600% for a Mars mission. The uncertainties in the quality factors are dominant. Using safety standards developed for low-Earth orbit, long-term space missions (>90 days) outside the Earth's magnetic field are currently unacceptable if the confidence levels in risk projections are considered. Because GCR exposures involve multiple particle or delta-ray tracks per cellular array, our results suggest that the shape of the dose response at low dose rates may be an additional uncertainty for estimating space radiation risks.

NASA Discipline Radiation Health↗

DNA Damage by Ionizing Radiation: Tandem Double Lesions by Charged Particles

Oxidative damages by ionizing radiation are the source of radiation-induced carcinogenesis, damage to the central nervous system, lowering of the immune response, as well as other radiation-induced damages to human health. Monte Carlo track simulations and kinetic modeling of radiation damages to the DNA employ available molecular and cellular data to simulate the biological effect of high and low LET radiation io the DNA. While the simulations predict single and double strand breaks and base damages, so far all complex lesions are the result of stochastic coincidence from independent processes. Tandem double lesions have not yet been taken into account. Unlike the standard double lesions that are produced by two separate attacks by charged particles or radicals, tandem double lesions are produced by one single attack. The standard double lesions dominate at the high dosage regime. On the other hand, tandem double lesions do not depend on stochastic coincidences and become important at the low dosage regime of particular interest to NASA. Tandem double lesions by hydroxyl radical attack of guanine in isolated DNA have been reported at a dosage of radiation as low as 10 Gy. The formation of two tandem base lesions was found to be linear with the applied doses, a characteristic of tandem lesions. However, tandem double lesions from attack by a charged particle have not been reported.

Huo, Winifred M.↗

MCNPX Cosmic Ray Shielding Calculations with the NORMAN Phantom Model

The United States is planning manned lunar and interplanetary missions in the coming years. Shielding from cosmic rays is a critical aspect of manned spaceflight. These ventures will present exposure issues involving the interplanetary Galactic Cosmic Ray (GCR) environment. GCRs are comprised primarily of protons (approx.84.5%) and alpha-particles (approx.14.7%), while the remainder is comprised of massive, highly energetic nuclei. The National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) has commissioned a joint study with Los Alamos National Laboratory (LANL) to investigate the interaction of the GCR environment with humans using high-fidelity, state-of-the-art computer simulations. The simulations involve shielding and dose calculations in order to assess radiation effects in various organs. The simulations are being conducted using high-resolution voxel-phantom models and the MCNPX[1] Monte Carlo radiation-transport code. Recent advances in MCNPX physics packages now enable simulated transport over 2200 types of ions of widely varying energies in large, intricate geometries. We report here initial results obtained using a GCR spectrum and a NORMAN[3] phantom.

James, Michael R.↗

A Flight Dynamics Perspective of the Orion Pad Abort One Flight Test

The Orion Crew Exploration Vehicle is America s next generation of human rated spacecraft. The Orion Launch Abort System will take the astronauts away from the exploration vehicle in the event of an aborted launch. The pad abort mode of the Launch Abort System will be flight-tested in 2009 from the White Sands Missile Range in New Mexico. This paper examines some of the efforts currently underway at the NASA Dryden Flight Research Center by the Controls & Dynamics group in preparation for the flight test. The concept of operation for the pad abort flight is presented along with an overview of the guidance, control and navigation systems. Preparations for the flight test, such as hardware testing and development of the real-time displays, are examined. The results from the validation and verification efforts for the aerodynamic and atmospheric models are shown along with Monte Carlo analysis results.

Idicula, Jinu↗