Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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 217 records · Page 12

Computation of rarefied hypersonic flows

Numerical techniques for the simulation of hypersonic flows of rarefied gases are examined in an analytical review. The direct-simulation Monte Carlo (DSMC) method developed to interpret measurement data obtained by the Space Shuttle in the SUMS project is described; the fundamental limitations of the DMSC approach are discussed; and modifications to improve the physical plausibility of DSMC predictions are proposed. Particular attention is given to the use of a double-peaked molecular distribution function for the internal flow in the SUMS probe, a downstream vacuum-reservoir condition as a simplifying assumption, and an explicit forward-time centered-space differencing scheme for the discretization of the SUMS problem. Typical simulation results are presented in extensive graphs and briefly characterized.

Cheng, Sin-I↗

Nonequilibrium molecular motion in a hypersonic shock wave

Molecular velocities have been measured inside a hypersonic, normal shock wave, where the gas experiences rapid changes in its macroscopic properties. As first hypothesized by Mott-Smith, but never directly observed, the molecular velocity distribution exhibits a qualitatively bimodal character that is derived from the distribution functions on either side of the shock. Quantitatively correct forms of the molecular velocity distribution function in highly nonequilibrium flows can be calculated, by means of the Direct Simulation Monte Carlo technique.

Pham-Van-diep, G.↗

Rotational and vibrational nonequilibrium effects in rarefied, hypersonic flow

Results are reported for an investigation into the methods by which energy transfer is calculated in the Direct Simulation Monte Carlo method. Description is made of a recently developed energy exchange model that deals with the translational and rotational modes. A new model for simulating the transfer of energy between the translational and vibrational modes is also explained. This model allows the vibrational relaxation time to follow the temperature dependence predicted by the Landau-Teller theory at moderate temperatures. For temperatures in excess of about 8000K the vibrational model is extended to include an empirical result for the relaxation time. The effect of introducing these temperature dependent collision numbers into the DSMC technique is assessed by making calculations representative of the stagnation streamline of a hypersonic space vehicle. Both thermal and chemical nonequilibrium effects are included while the flow conditions have been chosen such that ionization and radiation may be neglected. The introduction of these new models is found to significantly affect the degree of thermal nonequilibrium observed in the flowfield. Larger, and more widely ranging, differences in the results obtained with the different energy exchange probabilities are found when a significant amount of internal energy is included in the calculation of chemical nonequilibrium.

Boyd, Iain D.↗

Plume flowfield analysis of the shuttle primary Reaction Control System (RCS) rocket engine

A solution was generated for the physical properties of the Shuttle RCS 4000 N (900 lb) rocket engine exhaust plume flowfield. The modeled exhaust gas consists of the five most abundant molecular species, H2, N2, H2O, CO, and CO2. The solution is for a bare RCS engine firing into a vacuum; the only additional hardware surface in the flowfield is a cylinder (=engine mount) which coincides with the nozzle lip outer corner at X = 0, extends to the flowfield outer boundary at X = -137 m and is coaxial with the negative symmetry axis. Continuum gas dynamic methods and the Direct Simulation Monte Carlo (DSMC) method were combined in an iterative procedure to produce a selfconsistent solution. Continuum methods were used in the RCS nozzle and in the plume as far as the P = 0.03 breakdown contour; the DSMC method was used downstream of this continuum flow boundary. The DSMC flowfield extends beyond 100 m from the nozzle exit and thus the solution includes the farfield flow properties, but substantial information is developed on lip flow dynamics and thus results are also presented for the flow properties in the vicinity of the nozzle lip.

Hueser, J. E.↗

Direct simulation of hypersonic rarefied flow about a delta wing

Three-dimensional hypersonic rarefied flow about a delta wing at angle of incidence has been studied using the direct simulation Monte Carlo technique. Results are obtained for a transitional flow case which has been investigated in a nitrogen wind-tunnel experiment. The simulations are performed using a nonreacting, single-species gas model that accounts for rotational and vibrational internal energies. The computations yield an attached leeside flow associated with supersonic expansion over the leeside of the wing. Results are presented for the computed flowfield and surface quantities and overall aerodynamic coefficients.

Celenligil, M. Cevdet↗

Direct simulation of reentry flows with ionization

The Direct Simulation Monte Carlo (DSMC) method is applied in this paper to the study of rarefied hypersonic reentry flows. The assumptions and simplifications involved with the treatment of ionization, free electrons and the electric field are investigated. A new method is presented for the calculation of the electric field and handling of charged particles with DSMC. In addition, a two-step model for electron impact ionization is implemented. The flowfield representing a 10 km/sec shock at an altitude of 65 km is calculated. The effects of the new modeling techniques on the calculation results are presented and discussed.

Carlson, Ann B.↗

Assessment of chemical nonequilibrium in rarefied hypersonic flow

Calculations are presented of chemical nonequilibrium in the rarefied hypersonic flowfields surrounding reentering spacecraft in the upper atmosphere. The direct simulation Monte Carlo method is the solution technique employed. Methods by which chemical reactions are included in this simulation scheme are focused upon. In particular, the coupling of vibrational and dissociative nonequilibrium is investigated. The assumption of a fully excited mode has very little effect on the calculated results when the vibrational mode does not contribute to the energy available for reaction. A steric factor for dissociative reactions is developed which includes the vibrational energy of the dissociating molecule. Calculations are presented for the flow along the stagnation streamline of a blunt body for conditions representative of the entry and exit phases of the Aeroassisted Flight Experiment vehicle at an altitude of 90 km. It is found that both the flow properties and surface quantities are significantly affected by the introduction of the steric factor.

Boyd, Iain D.↗

Influence of local configuration on the backflow from small rocket thrusters

A numerical study has been made of the low density plume flow from a small thruster motor. For the case studied, the addition of a thin plume shield with a diameter 1.67 times the nozzle exit diameter was found to lead to a reduction in the backflow by several orders of magnitude. Calculations were also made for the three-dimensional flow produced by a pair of motors, in order to determine the increase in the backflow due to the interaction of the plumes. The calculations were made with the direct simulation Monte Carlo (DSMC) method, and one of the purposes of the work was to test a new and faster set of algorithms for this method.

Bird, G. A.↗

Assessment of impingement effects in the isentropic core of a small satellite control thruster plume

Two computational techniques commonly employed in the calculation of rocket and thruster expansion plumes are assessed. These are the method of characteristics (MOC), which is derived from the continuum Euler equations, and the direct simulation Monte Carlo (DSMC) method, which adopts a discrete particle approach. These techniques vary both in the computational expense and in the accuracy and detail of the solutions that they provide, depending upon the regime of application. The assessment is made with reference to the plume expanding from a small monopropellant hydrazine thruster and concentrates on the isentropic core of the jet for the flow regime lying between the continuum and free molecular limits. It is found that the more numerically intensive DSMC method offers the better correspondence to the available experimental data. In addition, large differences in typical impingement effects such as drag force and heat transfer are found at the free molecular limit of the plume expansion for the two predictive techniques. It is concluded that accurate estimation of impingement potential may only be achieved through application of the discrete particle method.

Boyd, I. D.↗

Grids for DSMC computation of two and three-dimensional flows

Contemporary grid schemes for use in Direct Simulation Monte Carlo (DSMC) method of computing two- and three-dimensional flows are briefly surveyed. The grids described include: (1) simple, constant area cells defined by a rectangular grid; (2) multilevel rectangular grids; (3) adaptive rectangular cells; (4) 'point reference' scheme; (5) body-fitted 'analytical' cells; (6) 'general' schemes; (7) 'general' schemes with subcells; (8) orthogonal transformation to rectangular grid; (9) structural tetrahedral grids; (10) rectangular multilevel subcells with adaptive body-fitted cells.

Bird, G. A.↗

Perception of numerical methods in rarefied gasdynamics

The relationships between various numerical methods applied to problems in rarefied gasdynamics are discussed, with emphasis on conflicting viewpoints and computational requirements associated with physical simulation versus the numerical solution of the Boltzmann equation. The basic differences between the molecular dynamics and direct simulation methods are shown to affect their applicability to dense and rarefied flows. Methods for the probabilistic selection of representative collision in the direct simulation Monte Carlo method are reviewed. A method combining the most desirable features of the earlier methods is presented.

Bird, G. A.↗

Validation of MCDS by comparison of predicted with experimental velocity distribution functions in rarefied normal shocks

Velocity distribution functions in normal shock waves in argon and helium are calculated using Monte Carlo direct simulation. These are compared with experimental results for argon at M = 7.18 and for helium at M = 1.59 and 20. For both argon and helium, the variable-hard-sphere (VHS) model is used for the elastic scattering cross section, with the velocity dependence derived from a viscosity-temperature power-law relationship in the way normally used by Bird (1976).

Pham-Van-diep, Gerald C.↗

Direct simulation of AFE forebody and wake flow with thermal radiation

Calculated results for the flowfield structure and surface quantities are presented for an axisymmetric representation of an aeroassist flight experiment vehicle. The direct simulation Monte Carlo (DSMC) method is used to perform the calculations, since the flow is highly nonequilibrium about the vehicle during both the compression and expansion phases. The body configuration is an elliptically blunt nose followed by a skirt with a circular radius and an afterbody. Freestream conditions correspond to a single point along the entry trajectory at an altitude of 90 km and a velocity of 9.9 km/s. The calculations account for nonequilibrium in the translational and internal modes, dissociation, ionization, and thermal radiation. The degree of dissociation is large, but the maximum ionization is only about 2 percent by mole fraction. The blunt forebody flow experiences a high degree of thermal nonequilibrium in which the translational temperature is generally greater than the internal temperature. However, as the flow expands about the aerobrake skirt and afterbody, the internal temperature is generally greater than the translational temperature.

Moss, James N.↗

Direct simulation of three-dimensional flow about the AFE vehicle at high altitudes

Three-dimensional hypersonic rarefied flow about the Aeroassist Flight Experiment (AFE) vehicle was studied using the direct simulation Monte Carlo (DSMC) technique. Results are presented for the transitional flow regime encountered between 120 and 200 km altitudes with a reentry velocity of 9.92 km/s. In the simulations, a five-species reacting real-gas model that accounts for internal energies (rotational and vibrational) is used. The results indicate that the transitional effects are significant even at an altitude of 200 km and influence the overall vehicle aerodynamics. For the cases considered, the aerodynamic coefficients, surface pressures, convective heating, and flow field structure variations with rarefaction effects are presented.

Celenligil, M. Cevdet↗

Rarefied flow past a flat plate at incidence

Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for the transitional flow about a flat plate at 40 deg incidence. The plate has zero thickness and a length of 1.0 m. The flow conditions simulated are those experienced by the Shuttle Orbiter during reentry at 7.5 km/s. The range of freestream conditions are such that the freestream Knudsen number values are between 0.02 and 8.4, i.e., conditions that encompass most of the transitional flow regime. The DSMC simulations show that transitional effects are evident when compared with free molecule results for all cases considered. The calculated results demonstrate clearly the necessity of having a means of identifying the effects of transitional flow when making aerodynamic flight measurements as are currently being made with the Space Shuttle Orbiter vehicles. Previous flight data analyses have relied exclusively on adjustments in the gas-surface interaction models without accounting for the transitional effect which can be comparable in magnitude. The present calculations show that the transitional effect at 175 km would increase the Space Shuttle Orbiter lift-drag ratio by 90 percent over the free molecule value.

Dogra, Virendra K.↗

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗

Shock waves in mixtures - A re-examination

The structure of normal shock waves in binary mixtures is analyzed, using Monte Carlo direct simulation, and compared with previous experimental data. Previous work by Erwin et al. used a collision method based directly on the intermolecular potential function; the present work extends this approach to binary mixtures, obtaining heteromolecular collision cross sections from the corresponding potentials. Results are presented for helium-argon mixtures, with conditions corresponding to the experiments of Harnett and Muntz (1972).

Erwin, Daniel A.↗

Bipropellant engine plume analyses for the Mars Observer scientific spacecraft

This paper discusses the methodology utilized to characterize exhaust plume expansion of the 490 N bipropellant engine and provides an analysis of exhaust plune impingement effects on the Mars Observer oxidizer tank structure. Nozzle and plume continuum flowfields are computed with a method of characteristics solution while transitional and rarefied regions of the exhaust gas flow were defined by a direct simulation Monte Carlo method. Results of this study provide sufficient data for an assessment of spacecraft thermal impacts due to plume heating effects and thrust losses from exhaust plume impingement.

Rattenni, L., Jr.↗