Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

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

Direct Monte Carlo Simulations of Hypersonic Viscous Interactions Including Separation

Results of calculations obtained using the direct simulation Monte Carlo method for Mach 25 flow over a control surface are presented. The numerical simulations are for a 35-deg compression ramp at a low-density wind-tunnel test condition. Calculations obtained using both two- and three-dimensional solutions are reviewed, and a qualitative comparison is made with the oil flow pictures highlight separation and three-dimensional flow structure.

Moss, James N.↗
NASA NTRS1996

Modeling of micro thrusters for gravity probe B

The concept of testing Einstein's general theory of relativity by means of orbiting gyroscopes was first proposed in 1959, which lead to the development of the Gravity Probe B experiment. Einstein's theory concerns the predictions of the relativistic precession of a gyroscope in orbit around earth. According to his theory, there will be two precessions due to the warping of space-time by the earth's gravitational field: the geodetic precession in the plane of the orbit, and the frame-dragging effect, in the direction of earth rotation. For a polar orbit, these components are orthogonal. In order to simplify the measurement of the precessions, Gravity Probe B (GP-B) will be placed in a circular polar orbit at 650 km, for which the predicted precessions will be 6.6 arcsec/year (geodetic) and 42 milli-arcsec/year (frame-dragging). As the gyroscope precesses, the orientation of its spin-axis will be measured with respect to the line-of-sight to Rigel, a star whose proper motion is known to be within the required accuracy. The line-of-sight to Rigel will be established using a telescope, and the orientation of the gyroscope spin axis will be measured using very sensitive SQUID (Superconducting Quantum Interference Device) magnetometers. The four gyroscopes will be coated with niobium. Below 2K, the niobium becomes superconducting and a dipole field will be generated which is precisely aligned with the gyroscope spin-axis. The change in orientation of these fields, as well as the spin-axis, is sensed by the SQUID magnetometers. In order to attain the superconducting temperatures for the gyroscopes and the SQUID's, the experiment package will be housed in a dewar filled with liquid helium. The helium flow through a GP-B micro thruster and into a vacuum is investigated using the Direct Simulation Monte Carlo method.

Jones, Kenneth M.↗
NASA NTRS1997

Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices

A study is made of the issues surrounding prediction of microchannel flows using the direct simulation Monte Carlo method. This investigation includes the introduction and use of new inflow and outflow boundary conditions suitable for subsonic flows. A series of test simulations for a moderate-size microchannel indicates that a high degree of grid under-resolution in the streamwise direction may be tolerated without loss of accuracy. In addition, the results demonstrate the importance of physically correct boundary conditions, as well as possibilities for reducing the time associated with the transient phase of a simulation. These results imply that simulations of longer ducts may be more feasible than previously envisioned.

Nance, Robert P.↗
NASA NTRS1998

Magellan Aerodynamic Characteristics During the Termination Experiment Including Thruster Plume-Free Stream Interaction

Results are presented on the aerodynamic characteristics of the Magellan spacecraft during the October 1994 Termination Experiment, including the effects of the thruster engine exhaust plumes upon the molecular free stream around the spacecraft and upon the aerodynamics coefficients. As Magellan passed through the Venusian atmosphere, the solar arrays were turned in opposite directions relative to the free stream creating a torque on the spacecraft. The spacecraft control system was programmed to counter the effects of this torque with attitude control engines to maintain an inertially fixed attitude. The orientation and reaction engine telemetry returned from Magellan are used to create a model of the aerodynamic torques. Geometric models of the Magellan spacecraft are analyzed with the aid of both free molecular and Direct Simulation Monte Carlo codes. The simulated aerodynamic torques determined are compared to the measured torques. The Direct Simulation Monte Carlo method is also used to model the attitude engine exhaust plumes, the free stream disturbance caused by these plumes, and the resulting torques acting on the spacecraft compared to no-exhaust plume cases. The effect of the exhaust plumes was found to be sufficiently large that thrust reversal is possible.

Cestero, Francisco J.↗
NASA NTRS1998

Calculation of Area Loss Due to Meteoroid Penetration

Simple Monte Carlo methods have been applied to the calculation of meteoroid penetration of spacecraft surfaces to account for the distributions of impact angle and velocity on the impact mechanics. Both single sheet and double sheet surfaces have been modeled. The model has been applied to several spacecraft and spacecraft concepts to evaluate light leakage into optical instruments or through sunshades. Results indicate that surfaces are less likely to be penetrated - or will be penetrated less often - than indicated by simply assuming an "average" velocity and normal impact. The reduction is typically about 40 percent. This effect is caused by the angular dependence of the penetration mechanics. However, the probability distribution of penetrating particles is relatively broad and indicates a peak above the ballistic limit particle size associated with the average velocity case.

Cooke, William J., Jr.↗
NASA NTRS2000

Forward Monte Carlo Computations of Polarized Microwave Radiation

Microwave radiative transfer computations continue to acquire greater importance as the emphasis in remote sensing shifts towards the understanding of microphysical properties of clouds and with these to better understand the non linear relation between rainfall rates and satellite-observed radiance. A first step toward realistic radiative simulations has been the introduction of techniques capable of treating 3-dimensional geometry being generated by ever more sophisticated cloud resolving models. To date, a series of numerical codes have been developed to treat spherical and randomly oriented axisymmetric particles. Backward and backward-forward Monte Carlo methods are, indeed, efficient in this field. These methods, however, cannot deal properly with oriented particles, which seem to play an important role in polarization signatures over stratiform precipitation. Moreover, beyond the polarization channel, the next generation of fully polarimetric radiometers challenges us to better understand the behavior of the last two Stokes parameters as well. In order to solve the vector radiative transfer equation, one-dimensional numerical models have been developed, These codes, unfortunately, consider the atmosphere as horizontally homogeneous with horizontally infinite plane parallel layers. The next development step for microwave radiative transfer codes must be fully polarized 3-D methods. Recently a 3-D polarized radiative transfer model based on the discrete ordinate method was presented. A forward MC code was developed that treats oriented nonspherical hydrometeors, but only for plane-parallel situations.

Battaglia, A.↗
NASA NTRS1999

Fluid/Solid Boundary Conditions in Non-Isothermal Systems

The existing theoretical research concerned with thermal creep at fluid/solid interfaces is briefly reviewed, and the importance of microgravity-based experimental data is then discussed. It is noted that the ultimate goal of this research is a rational molecular level theory that predicts the dependence of a dimensionless thermal creep coefficient, Ctc, on relevant dimensionless parameters describing the way fluid molecules interact with the solid surface and how they interact among themselves. The discussion covers thermophoresis of isolated solid spheres and aggregates in gases; solid sphere thermophoresis in liquids and dense vapors; thermophoresis of small immiscible liquid droplets; and applications of the direct simulation Monte Carlo method.

Rosner, Daniel E.↗
NASA NTRS2001

Transport Phenomena of Off-Axis Sputtering Deposition

Various high quality epitaxial films, especially oxides, have been synthesized using off-axis sputtering deposition. In this presentation, we report the experiment results of ZnO films grown by the off-axis sputtering deposition. Films were synthesized in temperatures ranged from room temperature to 600 C, and pressures from 5 mTorr to 150 mTorr. Film growth rate was measured by surface profilometer, ellipsometer, and wavelength dispersive spectrometry. Due to the collisions between the sputtered species and the residue gases, the kinetic energy of species was reduced and the transport of depositing species changed from a ballistic movement for low pressure to a diffusive drift for high pressure in which the transport species were almost thermalized. The measurements show an increase of growth rates along the gravity vector when the Knudsen number of transport species is less than 0.05, which suggests that gravity affected the transport characterization in off-axis sputtering deposition. Because the product of pressure (p) and travel distance (d) of sputtered species, pd, exceeds several mTorr-cm during film deposition, the classical simulations for sputtering process in high vacuum system may not be applied. Based on these experimental measurements, a transport process of the off-axis sputtering deposition is proposed. Several models including the Monte Carlo method and gravity-driven flow dynamics simulation will be discussed.

Zhu, S.↗
NASA NTRS1999

Characterization of Showerhead Performance At Low Pressure

The overall objective of this work is to characterize the flow to rough showerheads by deriving pressure drop versus velocity correlations to at can be then used in reactor scale simulations where the showerhead is approximated as a porous medium. At relatively low Reynolds numbers (less than 1-10 based on the hole length scale) and in the absence of slip flow, Darcy's Law, grad P = mu U/k, can be used to express the relation between the pressure drop and velocity where @mu@ is the fluid viscosity and it is the permeability that can be theoretically predicted as k= e R^2 /8, where e is the porosity. However, at sufficiently small hole diameters and decreased pressures (less than 5 Torr), the Knudsen number based on showerhead tube radius increases, and the flow may be in a transition regime. Different expressions have been proposed to account for this effect in the permeability by expressing k as a function of either pressure or Knudsen number. But at even higher Knudsen numbers, the pressure drop - velocity dependence is non-linear, and Darcy's Law no longer holds such that a permeability cannot be defined. The direct simulation Monte Carlo method is used along side conventional CFD techniques to determine the extent to which the CFD technique is appropriate and helps to derive correlations for the more rarefied cases of interest in these showerhead flows.

Hash, David B.↗
NASA NTRS2002

The Frequency Detuning Correction and the Asymmetry of Line Shapes: The Far Wings of H2O-H2O

A far-wing line shape theory which satisfies the detailed balance principle is applied to the H2O-H2O system. Within this formalism, two line shapes are introduced, corresponding to band-averages over the positive and negative resonance lines, respectively. Using the coordinate representation, the two line shapes can be obtained by evaluating 11-dimensional integrations whose integrands are a product of two factors. One depends on the interaction between the two molecules and is easy to evaluate. The other contains the density matrix of the system and is expressed as a product of two 3-dimensional distributions associated with the density matrices of the absorber and the perturber molecule, respectively. If most of the populated states are included in the averaging process, to obtain these distributions requires extensive computer CPU time, but only have to be computed once for a given temperature. The 11-dimensional integrations are evaluated using the Monte Carlo method, and in order to reduce the variance, the integration variables are chosen such that the sensitivity of the integrands on them is clearly distinguished.

Ma, Q.↗
NASA NTRS1997

A DSMC Study of Low Pressure Argon Discharge

Work toward a self-consistent plasma simulation using the DSMC (Direct Simulation Monte Carlo) method for examination of the flowfields of low-pressure high density plasma reactors is presented. Presently, DSMC simulations for these applications involve either treating the electrons as a fluid or imposing experimentally determined values for the electron number density profile. In either approach, the electrons themselves are not physically simulated. Self-consistent plasma DSMC simulations have been conducted for aerospace applications but at a severe computational cost due in part to the scalar architectures on which the codes were employed. The present work attempts to conduct such simulations at a more reasonable cost using a plasma version of the object-oriented parallel Cornell DSMC code, MONACO, on an IBM SP-2. Due to availability of experimental data, the GEC reference cell is chosen to conduct preliminary investigations. An argon discharge is chosen to conduct preliminary investigations. An argon discharge is examined thus affording a simple chemistry set with eight gas-phase reactions and five species: Ar, Ar(+), Ar(*), Ar(sub 2), and e where Ar(*) is a metastable.

Hash, David B.↗
NASA NTRS2002

High Resolution Soil Water from Regional Databases and Satellite Images

This viewgraph presentation provides information on the ways in which plant growth can be inferred from satellite data and can then be used to infer soil water. There are several steps in this process, the first of which is the acquisition of data from satellite observations and relevant information databases such as the State Soil Geographic Database (STATSGO). Then probabilistic analysis and inversion with the Bayes' theorem reveals sources of uncertainty. The Markov chain Monte Carlo method is also used.

Morris, Robin D.↗
NASA NTRS2002

Differential Cratering of Synchronously Rotating Satellites by Ecliptic Comets

We use Monte Carlo methods to simulate impacts of ecliptic comets on the synchronously rotating satellites of giant planets. We reconfirm the long-standing prediction that the cratering rate should be much higher on the leading hemispheres than on the trailing hemisphere; indeed we find that previously published analytical formulations modestly underestimate the degree of apex-antapex asymmetry to be expected. We then compare our results to new mapping of impact craters on Ganymede, Callisto, and Triton. Ganymede reveals a pronounced apex-antapex asymmetry that is nonetheless much less than predicted. All of Triton's confirmed impact craters are clustered toward the apex of motion, far exceeding the predicted asymmetry. No asymmetry is observed on Callisto. In each case at least one of our basic assumptions must be wrong. Likely candidates include the following: (i) the surfaces of all but the most sparsely cratered satellites are saturated or nearly saturated with impact craters; (ii) these satellites have rotated nonsynchronously over geological time; (iii) most of the craters are made not by heliocentric (Sun-orbiting) comets and asteroids but rather by planetocentric (planet-orbiting) debris of indeterminate origin; or (iv) pathological endogenic resurfacing has created illusions of structure. Callisto's surface is readily classified as nearly saturated. Ganymede's bright terrains, although less heavily cratered than those of Callisto, can also be explained by crater densities approaching saturation on a world where endogenic processes were active. The leading alternative is nonsynchronous rotation, an explanation supported by the distribution of catenae (crater chains produced by impact of tidally disrupted comets). Triton's craters can be explained by planetocentric debris or by capricious resurfacing, but both hypotheses are inherently improbable.

Zahnle, Kevin↗
NASA NTRS1999

Investigating Atmospheric Effects on Impact Ejecta Morphology: Possible Tool for Determining Past Climate Conditions on Mars?

The combined use of impact crater morphology and mechanics provides important information on the physical conditions of both planetary atmospheres and planetary and asteroid surfaces present during crater formation, while an understanding of the rate of crater production on the surface of asteroids provides information of their surface and spin rate evolution. The research performed with support from this project improves our understanding of (1) the mechanics of impact cratering in order to gain insights on the evolution of these physical surface conditions on planets with atmospheres and asteroids, and (2) how impact flux across an asteroid surface may vary due to anisotropic distribution of impactors in the solar system. As part of this project, we have undertaken three studies. In the first study, we investigate atmospheric effects on the morphology of ejecta excavated during a cratering event in order to determine the atmospheric and target conditions from observed crater morphologies. In the second study, we use the physical and morphological consequences of oblique impacts on an asteroid to understand how the asteroid Mathilde (recently imaged by the Near Earth Asteroid Rendezvous - NEAR- spacecraft) could have survived the formation of five giant craters. In a third study, we use a Monte Carlo method to calculate the impact flux on an asteroid given a distribution of impactors on elliptical orbits. In the following section, we present the result obtained from all three studies.

Appleby, John F.↗
NASA NTRS2000

Filtering in Hybrid Dynamic Bayesian Networks

We implement a 2-time slice dynamic Bayesian network (2T-DBN) framework and make a 1-D state estimation simulation, an extension of the experiment in (v.d. Merwe et al., 2000) and compare different filtering techniques. Furthermore, we demonstrate experimentally that inference in a complex hybrid DBN is possible by simulating fault detection in a watertank system, an extension of the experiment in (Koller & Lerner, 2000) using a hybrid 2T-DBN. In both experiments, we perform approximate inference using standard filtering techniques, Monte Carlo methods and combinations of these. In the watertank simulation, we also demonstrate the use of 'non-strict' Rao-Blackwellisation. We show that the unscented Kalman filter (UKF) and UKF in a particle filtering framework outperform the generic particle filter, the extended Kalman filter (EKF) and EKF in a particle filtering framework with respect to accuracy in terms of estimation RMSE and sensitivity with respect to choice of network structure. Especially we demonstrate the superiority of UKF in a PF framework when our beliefs of how data was generated are wrong. Furthermore, we investigate the influence of data noise in the watertank simulation using UKF and PFUKD and show that the algorithms are more sensitive to changes in the measurement noise level that the process noise level. Theory and implementation is based on (v.d. Merwe et al., 2000).

Andersen, Morten Nonboe↗
NASA NTRS1997

DSMC Simulations of Blunt Body Flows for Mars Entries: Mars Pathfinder and Mars Microprobe Capsules

The hypersonic transitional flow aerodynamics of the Mars Pathfinder and Mars Microprobe capsules are simulated with the direct simulation Monte Carlo method. Calculations of axial, normal, and static pitching coefficients were obtained over an angle of attack range comparable to actual flight requirements. Comparisons are made with modified Newtonian and free-molecular-flow calculations. Aerothermal results were also obtained for zero incidence entry conditions.

Moss, James N.↗