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

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

Enhancements and applications of DSMC for hypersonic rarefied flows

A review of the direct simulation Monte Carlo (DSMC) method of Bird is presented. The DSMC method provides the capability of simulating real gas flows in the rarefield flow regime. Recent developments and applications of the method for hypersonic flows are reported for both ground-based tests and during entry. The flows considered were three-dimensional, and both an existing general code and a new code under development were used.

Moss, James N.↗

Hypersonic flow - Physical issues and computational simulation

Computational results are presented on hypersonic flows, with emphasis on thermal and chemical nonequilibrium effects and the role played by grids. The continuum equations and their boundary conditions are presented, and the direct simulation Monte Carlo method is described. The importance of using the correct computational and physical models in describing hypersonic nonequilibrium flows is demonstrated.

Hassan, H. A.↗

Aeroassist flight experiment aerodynamics and aerothermodynamics

The problem is to determine the transitional flow aerodynamics and aerothermodynamics, including the base flow characteristics, of the Aeroassist Flight Experiment (AFE). The justification for the computational fluid dynamic (CFD) Application stems from MSFC's system integration responsibility for the AFE. To insure that the AFE objectives are met, MSFC must understand the limitations and uncertainties of the design data. Perhaps the only method capable of handling the complex physics of the rarefied high energy AFE trajectory is Bird's Direct Simulation Monte Carlo (DSMC) technique. The 3-D code used in this analysis is applicable only to the AFE geometry. It uses the Variable Hard Sphere (VHS) collision model and five specie chemistry model available from Langley Research Center. The code is benchmarked against the AFE flight data and used as an Aeroassisted Space Transfer Vehicle (ASTV) design tool. The code is being used to understand the AFE flow field and verify or modify existing design data. Continued application to lower altitudes is testing the capability of the Numerical Aerodynamic Simulation Facility (NASF) to handle 3-D DSMC and its practicality as an ASTV/AFE design tool.

Brewer, Edwin B.↗

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗

Analysis of rotational nonequilibrium in standing shock waves of nitrogen

The one-dimensional standing shock wave is the simplest flow in which nonequilibrium effects may be considered. The direct-simulation Monte Carlo method (DSMC) has yielded excellent agreement with results reported for the case of nitrogen flow at Mach 1.7. The DSMC technique is presently used in conjunction with a variable energy transfer probability model in which energy is transferred between the translational and rotational modes via the Borgnakke-Larsen (1975) phenomenological model.

Boyd, Iain D.↗

Local and overall aerodynamic coefficients for bodies in hypersonic, rarefied flow

A computational method is given for the prediction of local pressure and viscous shear stress on windward surfaces of bluff, convex, axisymmetric or quasi-axisymmetric, hypersonic bodies in the transitional, rarefied flow regime. Overall aerodynamic forces and moments are computed by integration of the local quantities. The method is based upon a correlation of local pressure and shear stress computed by the direct simulation Monte Carlo (DSMC) numerical technique for cold wall, real gas conditions and some supplemental data from low-density, hypersonic wind tunnels. The relative simplicity of the method makes it feasible to do the necessary calculations with a personal computer. Two-dimensional shapes and leeward surfaces are not included in the scope of the method as it is presented here. Results are compared with DSMC computations for both local and overall coefficients. The latter includes sphere and blunt cone drag as well as lift and pitching moment coefficients for the NASA AFE vehicle at various angles of attack. Very satisfactory agreement is shown.

Potter, J. Leith↗