Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Parallel implementation of a particle simulation for modeling rarefied gas dynamic flow

When the conditions of flow are rarefied and hypersonic, a more suitable alternative to the use of the Navier-Stokes equations for developing a numerical solution is the Direct Simulation Monte Carlo method (DSMC), a method of simulation which employs a large number of particles in modeling a rarefied gas. The performance of a parallel DSMC code developed for the Intel iPSC/860 Touchstone Gamma prototype computer is studied and the scaleup is found to be very nearly over the range of 16-128 processors.

Fallavollita, M. A.↗
NASA NTRS1993

Near wake structure for a generic ASTV configuration

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using the direct simulation Monte Carlo method. Particular emphasis is given to the near wake flow and its sensitivity to rarefaction and other parametric variations. The flow conditions simulated are attainable in existing low-density hypersonic wind tunnels; that is, Mach 20 nitrogen flow encompassing freestream Knudsen numbers of 0.03 to 0.001. A stable vortex forms in the near wake at and below a freestream Knudsen number of 0.01 and the size of the vortex increases with decreasing freestream Knudsen number. The base region of the flow remains in thermal nonequilibrium for all cases. There is no formation of a lip separation shock or a distinct wake shock at these rarefied conditions.

Dogra, Virendra K.↗
NASA NTRS1993

Model scattering envelopes of young stellar objects. II - Infalling envelopes

We present scattered light images for models of young stellar objects surrounded by dusty envelopes. The envelopes are assumed to have finite angular momentum and are falling in steady flow onto a disk. The model envelopes include holes, such as might be created by energetic bipolar flows. We calculate images using the Monte Carlo method to follow the light scattered in the dusty envelope and circumstellar disk, assuming that the photons originate from the central source. Adopting typical interstellar medium dust opacities and expected mass infall rates for protostars of about 10 exp -6 solar mass/yr, we find that detectable amounts of optical radiation can escape from envelopes falling into a disk as small as about 10-100 AU, depending upon the viewing angle and the size of the bipolar flow cavity. We suggest that the extended optical and near-IR light observed around several young stars is scattered by dusty infalling envelopes rather than disks.

Whitney, Barbara A.↗
NASA NTRS1993

Temperature dependence of rotational relaxation in shock waves of nitrogen

Computations are carried out for one-dimensional shock waves of diatomic nitrogen, using the direct simulation Monte Carlo method. It was found that the reciprocal shock thickness varies with the upstream temperature condition. This variation was also observed in the experimental data and was simulated numerically by using the temperature-dependent expression for the rotational collision number.

Boyd, I. D.↗
NASA NTRS1993

Spacecraft contamination investigation by direct simulation Monte Carlo - Contamination on UARS/HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to study the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite, with the ultimate goal of estimating the cumulative level of molecular deposits on critical optical surfaces of the Halogen Occultation Experiment. Emphasis is placed on describing the development of preprocessing codes that simplify the incorporation of geometry and outgassing data into the main run code. A novel use of parallel processing software allows use of multiple processors to solve different regions of the overall flowfield. Initial results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the Halogen Occultation telescope's aperture plane are presented, along with data related to code performance.

Rault, Didier F. G.↗
NASA NTRS1993

A generalized hard-sphere model for Monte Carlo simulation

A new molecular model, called the generalized hard-sphere, or GHS model, is introduced. This model contains, as a special case, the variable hard-sphere model of Bird (1981) and is capable of reproducing all of the analytic viscosity coefficients available in the literature that are derived for a variety of interaction potentials incorporating attraction and repulsion. In addition, a new procedure for determining interaction potentials in a gas mixture is outlined. Expressions needed for implementing the new model in the direct simulation Monte Carlo methods are derived. This development makes it possible to employ interaction models that have the same level of complexity as used in Navier-Stokes calculations.

Hassan, H. A.↗
NASA NTRS1993

Transient motion of a confined rarefied gas due to wall heating or cooling

The direct Monte Carlo method is used to simulate the transient motion of a confined rarefied gas as a container wall is rapidly heated or cooled. Results show that, for the rapid wall temperature changes, the gas response is better measured by an acoustic, rather than conduction, timescale, with steady state conditions being achieved after about ten wave crossing times.

Wadsworth, Dean C.↗
NASA NTRS1993

Kinetic relaxation of a non-Maxwellian monatomic gas in a state of gross rest

The kinetic relaxation of a monatomic gas in a state of gross rest is studied by use of Maxwell's second-order moment equation. The molecular models considered range from the hard sphere to the Maxwell molecule. The development is exact for an ellipsoidal distribution and becomes a good approximation for more general distributions. The results show that the relaxation of second-order moments is nearly exponential for power-law molecules and that the characteristic time mu/p appears as the principal controlling parameter while the power-law constant plays a secondary role. Comparisons between results of numerical simulations using Bird's direct simulation Monte Carlo method and theory show excellent agreement for both hard-sphere and Maxwell molecules.

Baganoff, D.↗
NASA NTRS1993

Particulate multi-phase flowfield analysis for advanced solid rocket motor

Particulate multi-phase flowfield with chemical reaction for a 2D advanced solid rocket motor (ASRM) is analyzed using the finite difference Navier-Stokes (FDNS) code. The flowfield in the aft dome cavity of the ASRM is examined and its significant impact on the motor operation and performance is demonstrated. Chemical reaction analysis is performed for H2O, O2, H2, O, H, OH, CO, CO2, Cl, Cl2, HCl, and N2. The turbulent dispersion effect is calculated with the Monte Carlo method. Result show that a recirculation zone exists at the entry of the aft-dome cavity. The particle impingement could cause the erosion and damage nozzle wall. Accumulating in the impingement area the particles change the wall shape and affect the motor performance.

Liaw, Paul↗
NASA NTRS1993

On predicting contamination levels of HALOE optics aboard UARS using direct simulation Monte Carlo

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 flowfield and surface conditions and geometric orientations in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. Problems resolving species outgassing and vent flux rates that varied over many orders of magnitude were handled using species weighting factors. Results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface are presented, along with data related to code performance. Using procedures developed in standard contamination analyses, the cumulative level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated to be about 2700A.

Woronowicz, Michael S.↗
NASA NTRS1993

Analysis of plume backflow around a nozzle lip in a nuclear rocket

The structure of the flow around a nuclear thermal rocket nozzle lip has been investigated using the direct simulation Monte Carlo method. Special attention has been paid to the behavior of a small amount of harmful particles that may be present in the rocket exhaust gas. The harmful fission product particles are modeled by four inert gases whose molecular weights are in a range of 4 131. Atomic hydrogen, which exists in the flow due to the extremely high nuclear fuel temperature in the reactor, is also included. It is shown that the plume backflow is primarily determined by the thin subsonic fluid layer adjacent to the surface of the nozzle lip, and that the inflow boundary in the plume region has negligible effect on the backflow. It is also shown that a relatively large amount of the lighter species is scattered into the backflow region while the amount of the heavier species becomes negligible in this region due to extreme separation between the species. Results indicate that the backscattered molecules are very energetic and are fast-moving along the surface in the backflow region near the nozzle lip.

Chung, Chan H.↗
NASA NTRS1992

On the decay of the magnetic fields of single radio pulsars

We investigate the statistical evidence for the decay of the magnetic field of single radio pulsars. We perform population syntheses for different assumed values for the time scale of field decay using a Monte Carlo method. We allow for the selection effects in pulsar surveys and compare the synthesized populations with the observed pulsars. We take account of the finite scale height of the distribution in the Galaxy of free electrons, which determine the dispersion measure and hence the apparent distance of radio pulsars. Our simulations give much better agreement with the observations if the time scale for the field decay is assumed to be longer than the typical active life time of a radio pulsar. This indicates that no significant field decay occurs in single radio pulsars.

Bhattacharya, Dipankar↗
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.↗