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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,639 records · Page 91

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Molecular gas dynamics applied to low-thrust propulsion

The Direct Simulation Monte Carlo method is currently being applied to study flowfields of small thrusters, including both the internal nozzle and the external plume flow. The DSMC method is employed because of its inherent ability to capture nonequilibrium effects and proper boundary physics in low-density flow that are not readily obtained by continuum methods. Accurate prediction of both the internal and external nozzle flow is important in determining plume expansion which, in turn, bears directly on impingement and contamination effects.

Zelesnik, Donna↗

NASA NTRS1993

Numerical solutions of the complete Navier-Stokes equations

The objective of this study is to compare the use of assumed pdf (probability density function) approaches for modeling supersonic turbulent reacting flowfields with the more elaborate approach where the pdf evolution equation is solved. Assumed pdf approaches for averaging the chemical source terms require modest increases in CPU time typically of the order of 20 percent above treating the source terms as 'laminar.' However, it is difficult to assume a form for these pdf's a priori that correctly mimics the behavior of the actual pdf governing the flow. Solving the evolution equation for the pdf is a theoretically sound approach, but because of the large dimensionality of this function, its solution requires a Monte Carlo method which is computationally expensive and slow to coverage. Preliminary results show both pdf approaches to yield similar solutions for the mean flow variables.

Hassan, H. A.↗

NASA NTRS1994

Stochastic solution to quantum dynamics

The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. The technique, which is developed and validated in two- and three- dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space.

John, Sarah↗

NASA NTRS1992

Computation of flowfields for hypersonic flight at high altitudes

Recent direct simulation Monte Carlo methods are reviewed focusing on the nonequilibrium aspects of rarefied hypersonic flows typical of reentry conditions. The calculations identify the altitudes for which the onset of dissociation occurs, the magnitude of the surface temperature variations, the effect of the shock wave structure on the chemistry of a dissociating gas, and the radiation emisssion for reentry flowfields.

Moss, James N.↗

NASA NTRS1994

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

NASA NTRS1994

X-ray bumps, iron K-alpha lines, and X-ray suppression by obscuring tori in Seyfert galaxies

We investigate the X-ray spectral properties of unobscured type 1 and obscured type 2 Seyferts as predicted by the unified Seyfert scheme. We consider the reprocessing of X-ray photons by photoelectric absorption, iron fluorescence, and Compton downscattering in the obscuring tori surrounding these active nuclei, and compute by Monte Carlo methods the reprocessed spectra as a function of the viewing angle. Depending on the optical depth and shape of the torus, and on the viewing angle, the X-ray flux can be suppressed by substantial factors when our line of sight is obscured. We show that an immediate consequence of the existence of an obscuring thick torus is the production in the spectra of type 1 Seyfert galaxies of a bump in the continuum above 10-20 keV and an Fe K-alpha line with significant equivalent width. In those type 2 Seyferts for which the hard X-ray spectrum has been substantially suppressed, the equivalent width of the Fe K-alpha line in the transmitted spectrum can be very large.

Krolik, Julian H.↗

NASA NTRS1995

Polarization from scattering in blobs

We describe a Monte Carlo method for calculating polarization due to scattering in various geometries. We apply the method in this paper to the case of scattering in a spherical cloud illuminated by parallel rays. The polarization is largest in small optical depth blobs, where single scattering dominates. The polarization decreases as optical depth increases until a limit is reached where scattering occurs on the surface, and increasing optical depth has little effect on the results. For a given optical depth, decreasing the albedo leads to an increase of the degree of polarization of the scattered light because single scattering increases relative to multiple scattering. As the optical depth of the blob increases, the scattered flux becomes increasingly backward throwing. In the high albedo blobs, the maximum polarization becomes skewed towards the illuminated hemisphere at high optical depths. We discuss the applications to polarization in supergiants and R CrB stars, and scattering in a clumpy interstellar medium.

Code, Arthur D.↗

NASA NTRS1995

Discovery and imaging of a Galactic cirrus cloud with the far ultraviolet space telescope

We present new far-ultraviolet (1400-1800 A) data concerning a Galactic cirrus cloud G251.2+73.3 near the north Galactic pole obtained with the space-borne imaging telescope FAUST (Far Ultraviolet Space Telescope). We obtain a good correlation between the far-ultraviolet (FUV) and IRAS 100 micrometers surface brightnesses, their relation being I(sub FUV) = (128 +/- 3) I(sub 100 micrometers) - (264 +/- 9), where the I(sub FUV) flux is given in units of photon/s/sq cm/A/sr and I(sub 100 micrometers) in MJy/sr. Using uvbyH-beta photometry, we get a distance of 120 pc and a visual extinction in the center of the cloud of 0.39 mag corresponding to an extinction of 1.0 mag at 1565 A. We have performed a multiple scattering calculation for the scattered light using the Monte Carlo method. These calculations provide restrictions on the FUV scattering properties of the interstellar dust.

Haikala, Lauri K.↗

NASA NTRS1995

Computing Interactions Of Free-Space Radiation With Matter

High Charge and Energy Transport (HZETRN) computer program computationally efficient, user-friendly package of software adressing problem of transport of, and shielding against, radiation in free space. Designed as "black box" for design engineers not concerned with physics of underlying atomic and nuclear radiation processes in free-space environment, but rather primarily interested in obtaining fast and accurate dosimetric information for design and construction of modules and devices for use in free space. Computational efficiency achieved by unique algorithm based on deterministic approach to solution of Boltzmann equation rather than computationally intensive statistical Monte Carlo method. Written in FORTRAN.

Wilson, J. W.↗

NASA NTRS1995

System-Availability And Resource-Allocation Program

ACARA analyzes availability, life-cycle cost, and scheduling of resources. Uses statistical Monte Carlo method to simulate capacity states of system as well as failure and repair of components. Failures of components modeled mathematically by use of combination of exponential and Weibull probability distributions. Schedules replacement of components to optimize performance of system. Made to comply with any constraints on production of components, capacities of resupply vehicles, spares kept on site, crews, and/or equipment. Written in APL2.

Viterna, L. A.↗

NASA NTRS1995

HINCOF-1: a Code for Hail Ingestion in Engine Inlets

One of the major concerns during hail ingestion into an engine is the resulting amount and space- and time-wise distribution of hail at the engine face for a given geometry of inlet and set of atmospheric and flight conditions. The appearance of hail in the capture streamtube is invariably random in space and time, with respect to size and momentum. During the motion of a hailstone through an inlet, a hailstone undergoes several processes, namely impact with other hailstones and material surfaces of the inlet and spinner, rolling and rebound following impact; heat and mass transfer; phase change; and shattering, the latter three due to friction and impact. Taking all of these factors into account, a numerical code, designated HINCOF-I, has been developed for determining the motion hailstones from the atmosphere, through an inlet, and up to the engine face. The numerical procedure is based on the Monte-Carlo method. The report presents a description of the code, along with several illustrative cases. The code can be utilized to relate the spinner geometry - conical or, more effective, elliptical - to the possible diversion of hail at the engine face into the bypass stream. The code is also useful for assessing the influence of various hail characteristics on the ingestion and distribution of hailstones over the engine face.

Gopalaswamy, N.↗

NASA NTRS1992

Direct Monte Carlo Simulations of Hypersonic Low-Density Flows about an ASTV Including Wake Structure

Results of a numerical study concerning flow past a 70-deg blunted cone in hypersonic low-density flow environments are presented using the direct simulation Monte-Carlo method. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Results indicate that a stable vortex forms in the near wake at and below a freestream Knudsen number (based on cone diameter) of 0.01 and the size of the vortex increases with decreasing Knudsen number. The base region of the flow remains in thermal nonequilibrium for all cases considered herein.

Dogra, V. K.↗

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