Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct simulation Monte Carlo (DSMC)”

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 127 records · Page 7

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

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

Rarefied-continuum gas dynamics transition for SUMS project

This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.

Cheng, Sin-I↗

Monte Carlo simulation of an expanding gas

By application of simple computer graphics techniques, the statistical performance of two Monte Carlo methods used in the simulation of rarefied gas flows are assessed. Specifically, two direct simulation Monte Carlo (DSMC) methods developed by Bird and Nanbu are considered. The graphics techniques are found to be of great benefit in the reduction and interpretation of the large volume of data generated, thus enabling important conclusions to be drawn about the simulation results. Hence, it is discovered that the method of Nanbu suffers from increased statistical fluctuations, thereby prohibiting its use in the solution of practical problems.

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

Direct simulation Monte Carlo analysis on parallel processors

A method is presented for executing a direct simulation Monte Carlo (DSMC) analysis using parallel processing. The method is based on using domain decomposition to distribute the work load among multiple processors, and the DSMC analysis is performed completely in parallel. Message passing is used to transfer molecules between processors and to provide the synchronization necessary for the correct physical simulation. Benchmark problems are described for testing the method and results are presented which demonstrate the performance on two commercially available multicomputers. The results show that reasonable parallel speedup and efficiency can be obtained if the problem is properly sized to the number of processors. It is projected that with a massively parallel system, performance exceeding that of current supercomputers is possible.

Wilmoth, Richard G.↗

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↗

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

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

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

Adaptive domain decomposition for Monte Carlo simulations on parallel processors

A method is described for performing direct simulation Monte Carlo (DSMC) calculations on parallel processors using adaptive domain decomposition to distribute the computational work load. The method has been implemented on a commercially available hypercube and benchmark results are presented which show the performance of the method relative to current supercomputers. The problems studied were simulations of equilibrium conditions in a closed, stationary box, a two-dimensional vortex flow, and the hypersonic, rarefield flow in a two-dimensional channel. For these problems, the parallel DSMC method ran 5 to 13 times faster than on a single processor of a Cray-2. The adaptive decomposition method worked well in uniformly distributing the computational work over an arbitrary number of processors and reduced the average computational time by over a factor of two in certain cases.

Wilmoth, Richard G.↗

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