SEARCH · Engineering Papers

Results for “direct simulation 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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Monte Carlo analysis of dissociation and recombination behind strong shock waves in nitrogen

Computations are presented for the relaxation zone behind strong, 1D shock waves in nitrogen. The analysis is performed with the direct simulation Monte Carlo method (DSMC). The DSMC code is vectorized for efficient use on a supercomputer. The code simulates translational, rotational and vibrational energy exchange and dissociative and recombinative chemical reactions. A model is proposed for the treatment of three body-recombination collisions in the DSMC technique which usually simulates binary collision events. The model improves previous models because it can be employed with a large range of chemical-rate data, does not introduce into the flow field troublesome pairs of atoms which may recombine upon further collision (pseudoparticles) and is compatible with the vectorized code. The computational results are compared with existing experimental data. It is shown that the derivation of chemical-rate coefficients must account for the degree of vibrational nonequilibrium in the flow. A nonequilibrium-chemistry model is employed together with equilibrium-rate data to compute the flow in several different nitrogen shock waves.

Boyd, I. D.↗

NASA NTRS1992

Monte Carlo simulation of reentry flows with ionization

The Direct Simulation Monte Carlo method is applied to a rarefied, weakly ionized, hypersonic flow over a blunt axisymmetric body. An ionization model based on the concept of ambipolar diffusion is used and a model for the sheath is presented. The effects of the new modeling techniques are investigated for flow over the Project Fire II configuration at 11.37 km/s at an altitude of 84.6 km. The calculated results are presented and compared with both experimental data and solutions where ionization effects were not included. In general, the calculated results overpredict the experimental values by about 15-20 percent.

Taylor, Jeff C.↗

NASA NTRS1992

Monte Carlo simulation of entry in the Martian atmosphere

The Direct Simulation Monte Carlo method of Bird is used to investigate the characteristics of low density hypersonic flowfields for typical aerobrakes during Martian atmospheric entry. The method allows for both thermal and chemical nonequilibrium. Results are presented for a sixty-degree spherically blunt cone for various nose radii and altitudes.

Hash, David B.↗

NASA NTRS1992

Hypersonic rarefied flow past spheres including wake structure

Results of a numerical study using the direct simulation Monte Carlo method are presented for hypersonic rarefied flow past spheres. The flow conditions considered are those corresponding to low density wind tunnel test conditions. The set of the experimental conditions for the calculations encompasses the transitional to near-continuum flow regimes. Comparison of the calculated drag with experimental results shows good agreement to well within the experimental error. Particular attention is focused on the wake structure. Calculations show that the wake is very rarefied with considerable thermal nonequilibrium for all the cases considered. No flow separation is observed in the wake for the near-continuum case where a vortex has been predicted by Navier-Stokes type calculations.

Dogra, Virendra K.↗

NASA NTRS1992

Internal structure of shock waves in disparate mass mixtures

The detailed flow structure of a normal shock wave for a gas mixture is investigated using the direct-simulation Monte Carlo method. A variable diameter hard-sphere (VDHS) model is employed to investigate the effect of different viscosity temperature exponents (VTE) for each species in a gas mixture. Special attention is paid to the irregular behavior in the density profiles which was previously observed in a helium-xenon experiment. It is shown that the VTE can have substantial effects in the prediction of the structure of shock waves. The variable hard-sphere model of Bird shows good agreement, but with some limitations, with the experimental data if a common VTE is chosen properly for each case. The VDHS model shows better agreement with the experimental data without adjusting the VTE. The irregular behavior of the light-gas component in shock waves of disparate mass mixtures is observed not only in the density profile, but also in the parallel temperature profile. The strength of the shock wave, the type of molecular interactions, and the mole fraction of heavy species have substantial effects on the existence and structure of the irregularities.

Chung, Chan-Hong↗

NASA NTRS1992

DSMC analysis of species separation in rarefied nozzle flows

The direct-simulation Monte Carlo method has been used to investigate the behavior of a small amount of a harmful species in the plume and the backflow region of nuclear thermal propulsion rockets. Species separation due to pressure diffusion and nonequilibrium effects due to rapid expansion into a surrounding low-density environment are the most important factors in this type of flow. It is shown that a relatively large amount of the lighter species is scattered into the backflow region and the heavier species becomes negligible in this region due to the extreme separation between species. It is also shown that the type of molecular interaction between the species can have a substantial effect on separation of the species.

Chung, Chan-Hong↗

NASA NTRS1992

A new two-temperature dissociation model for reacting flows

A new two-temperature dissociation model for flows undergoing compression is derived from kinetic theory. The model minimizes uncertainties associated with the two-temperature model of Park. The effects of the model on AOTV type flowfields are examined and compared with the Park model. Calculations are carried out for flows with and without ionization. When considering flows with ionization, a four temperature model is employed. For Fire II conditions, the assumption of equilibrium between the vibrational and electron-electronic temperatures is somewhat poor. A similar statement holds for the translational and rotational temperatures. These trends are consistent with results obtained using the direct simulation Monte Carlo method.

Olynick, David R.↗

NASA NTRS1992

Evaluation of thermochemical models for particle and continuum simulations of hypersonic flow

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered that are dominated in turn by vibrational relaxation, weak dissociation, strong dissociation and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Boyd, Iain D.↗

NASA NTRS1991

Monte Carlo study of vibrational relaxation processes

A new model is proposed for the computation of vibrational nonequilibrium in the direct simulation Monte Carlo method (DSMC). This model permits level to level vibrational transitions for the first time in a Monte Carlo flowfield simulation. The model follows the Landau-Teller theory for a harmonic oscillator in which the rates of transition are related to an experimental correlation for the vibrational relaxation time. The usual method for simulating such processes in the DSMC technique applies a constant exchange probability to each collision and the vibrational energy is treated as a continuum. A comparison of these two methods is made for the flow of nitrogen over a wedge. Significant differences exist for the vibrational temperatures computed. These arise as a consequence of the incorrect application of a constant exchange probability in the old method. It is found that the numerical performances of the two vibrational relaxation models are equal.

Boyd, Iain D.↗

NASA NTRS1993

Development and application of computational aerothermodynamics flowfield computer codes

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second-approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered. The cases are governed, respectively, by the following: vibrational relaxation; weak dissociation; strong dissociation; and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Venkatapathy, Ethiraj↗

NASA NTRS1993

Experimentation in the low-density plume of a simulated electrothermal thruster for computer code validation

Pressures and flow angles are measured in the plume of a 20 deg half-angle, conical nozzle in vacuum with Pitot tubes and conical probes. The area of measurement in the plume ranges from the nozzle exit plane to 480 mm axially downstream and from the plume centerline to 60 mm radially. The nozzle has an exit-to-throat area ratio of 100:1 and a throat diameter of 3.2 mm. The nozzle flow exhausts to a vacuum of order 10(exp -2) Pa to simulate a resistojet (an electrothermal rocket of less than 1 N of thrust) operating in space. Experimental data are given for flows of nitrogen at 55 and 68 mg/s, stagnation temperatures between 695 and 921 K, and stagnation pressures ranging from 5600 to 7100 Pa. Data are also given for argon at a rate of 68 mg/s, a stagnation temperature of 648 K, and stagnation pressures of 4500, 4750, and 4770 Pa. Measurements in the nitrogen plume are compared with computational results from a direct-simulation Monte Carlo method.

Meissner, Dana L.↗

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

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

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

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