Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10

Comparison of shock structure solutions using independent continuum and kinetic theory approaches

A vehicle traversing the atmosphere will experience flight regimes at high altitudes in which the thickness of a hypersonic shock wave is not small compared to the shock standoff distance from the hard body. When this occurs, it is essential to compute accurate flow field solutions within the shock structure. In this paper, one-dimensional shock structure is investigated for various monatomic gases from Mach 1.4 to Mach 35. Kinetic theory solutions are computed using the Direct Simulation Monte Carlo method. Steady-state solutions of the Navier-Stokes equations and of a slightly truncated form of the Burnett equations are determined by relaxation to a steady state of the time-dependent continuum equations. Monte Carlo results are in excellent agreement with published experimental data and are used as bases of comparison for continuum solutions. For a Maxwellian gas, the truncated Burnett equations are shown to produce far more accurate solutions of shock structure than the Navier-Stokes equations.

Fiscko, Kurt A.↗

Three-dimensional flow simulation about the AFE vehicle in the transitional regime

The direct-simulation Monte Carlo technique is used to analyze the hypersonic rarefied flow about the three-dimensional NASA Aeroassist Flight Experiment vehicle. Results are given for typical transitional flows encountered during the vehicle's atmospheric entry from altitudes of 200-100 km with an entry velocity of 9.9 km/s. It is found that dissociation is important at altitudes of 110 km and below, and that transitional effects are significant even at an altitude of 200 km.

Celenligil, M. Cevdet↗

Influence of afterbodies on AOTV flows

The role of vibration-dissociation coupling in flows representative of AOTV's with afterbodies is discussed. A physical model designed to reflect the behavior of anharmonic diatomic molecules at high temperatures was developed and incorporated into the direct simulation Monte Carlo method of Bird (1976). The model leads to reduced dissociation and increased convective heating for flows with finite catalytic surfaces.

Olynick, David P.↗

Direct simulation of rotational and vibrational nonequilibrium

The ways in which energy transfer is calculated in the Direct Simulation Monte Carlo method are presented. An energy exchange model that deals with translational and rotational modes is described. A model for simulating the transfer of energy between the translational and vibrational modes is presented as well.

Boyd, Iain D.↗

Computation of electron density in high altitude re-entry flows

The direct simulation Monte Carlo (DSMC) method has been modified to deal with the very small degree of ionization that occurs in the early stages of the 'communications blackout' phase of re-entry. The simplified model of the process assumes that the ionization is entirely a consequence of the reaction N + O yields NO(+) + e-. The measurements from the RAM C-II test vehicle at an altitude of 81 km have been chosen as the test case. It was found that the form of the electron distribution along the shock layer was in agreement with the measurements. That normal to the surface was also in agreement with the measurements if the surface is regarded as fully catalytic. However, quantitative agreement for the fully catalytic case requires a reaction rate about an order of magnitude greater than those generally associated with this reaction.

Bird, G. A.↗

'Hypersonic slip flows' and issues on extending continuum model beyond the Navier-Stokes level

Results of a direct simulation Monte Carlo method for a hypersonic flow about a flat plate at a 40 deg angle of attack have been compared with corresponding results from a theory for fully viscous shock layers (FVSLs). Using the 13-moment equations for a Maxwell gas, it is demonstrated that nonequilibrium and equilibrium FVSL flows can be correlated. With the exception of the pressure density, most of the flow properties along a streamline (including heat flux and shear and normal stresses) are correctly predicted to leading order by the Navier-Stokes model.

Cheng, H. K.↗

Hypersonic rarefied flow about plates at incidence

The direct-simulation Monte Carlo method has been used in a numerical study of the transitional flow about two plate configurations at incidence; one of the two plates, both of which are 12 m long, has zero thickness, while the other has a thickness of 0.5 m and a node radius of 0.5 m. The flow conditions simulated are those of the Space Shuttle Orbiter during 7.5 km/hr reentry, in the 200-100 km altitude range encompassing most of the transitional flow for this vehicle. The results obtained clearly demonstrate that transitional effects are significant even at those altitudes where the flow about a typical space vehicle has been considered free-molecular.

Dogra, Virendra K.↗

Study of hypersonic flow past sharp cones

Calculations using the direct simulation Monte Carlo (DSMC) method of Bird for flow past sharp cones in the near continuum to free molecule flow regime are presented and compared with experiment. It is found that results are sensitive to the grid and the interaction potential. Moreover, the time counter method was found to be as accurate as other methods when the solution is grid independent. Finally, the results show that the effects of the wake on the forebody surface properties are minimal.

Taylor, Jeff C.↗

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