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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

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

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↗
NASA NTRS1991

Aerothermodynamics of a 1.6-m-diameter sphere in hypersonic rarefied flow

Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for hypersonic rarefied flow about a 1.6-m-diameter sphere. The flow conditions considered are those experienced by a typical satellite in orbit or by a space vehicle during entry. The altitude range considered is that from 90 to 200 km, which encompasses the near continuum, transitional and free-molecular flow regimes. A freestream velocity of 7.5 km/s is assumed in the simulations. The results show that transitional effects are significant at all altitudes below 200 km, but at 200 km the flow about the sphere attains the free-molecular limit. Very little chemical activity is present above 120 km. Both the stagnation point heat transfer and the sphere drag approach their respective free molecule values at 200 km. Results highlight the thermal and chemical nonequilibrium nature of the flowfield. Nonequilibrium effects on the surface heating and body drag are also investigated.

Dogra, Virendra K.↗
NASA NTRS1991

Accuracy of the Burnett equations for hypersonic real gas flows

Burnett solutions for the structure of normal shock waves for Mach numbers ranging from 1.1 to 50 have been obtained by using the time-dependent equation set in conjunction with modern computational fluid dynamic techniques. The current work seeks to provide corrected and improved high-resolution direct simulation Monte Carlo (DSMC) results for hypersonic shocks in monatomic gases. These DSMC results will serve as a reliable test base for assessing the accuracy of various continuum equations more advanced than Navier-Stokes, such as the Burnett equations. The main conclusion of Fisko and Chapman (1988) that the Burnett equations are more accurate than the Navier-Stokes equations for all gases and at all Mach numbers in matching DSMC shock wave results is reassesed by comparing the improved and corrected DSMC results with continuum Burnett results. Various Burnett results for shocks in a diatomic gas, nitrogen, which include improvements to the work of Lumpkin and Chapman (1989) are presented. It is concluded that the Burnett equations represent a significant improvement over Navier-Stokes equations for flows exhibiting translational nonequilibrium.

Lumpkin, Forrest E., III↗
NASA NTRS1990

Vectorization of a particle code used in the simulation of rarefied hypersonic flow

A limitation of the direct simulation Monte Carlo (DSMC) method is that it does not allow efficient use of vector architectures that predominate in current supercomputers. Consequently, the problems that can be handled are limited to those of one- and two-dimensional flows. This work focuses on a reformulation of the DSMC method with the objective of designing a procedure that is optimized to the vector architectures found on machines such as the Cray-2. In addition, it focuses on finding a better balance between algorithmic complexity and the total number of particles employed in a simulation so that the overall performance of a particle simulation scheme can be greatly improved. Simulations of the flow about a 3D blunt body are performed with 10 to the 7th particles and 4 x 10 to the 5th mesh cells. Good statistics are obtained with time averaging over 800 time steps using 4.5 h of Cray-2 single-processor CPU time.

Baganoff, D.↗
NASA NTRS1991

Shock wave-boundary layer interactions in rarefied gas flows

A numerical study is presented, using the direct simulation Monte Carlo (DSMC) method, of shock wave-boundary layer interactions in low density supersonic flows. Test cases include two-dimensional, axially-symmetric and three-dimensional flows. The effective displacement angle of the boundary layer is calculated for representative flat plate, wedge, and cone flows. The maximum pressure, shear stress, and heat transfer in the shock formation region is determined in each case. The two-dimensional reflection of an oblique shock wave from a flat plate is studied, as is the three-dimensional interaction of such a wave with a sidewall boundary layer.

Bird, G. A.↗
NASA NTRS1991

Hypersonic rarefied flow about a delta wing - Direct simulation and comparison with experiment

Three-dimensional simulations of hypersonic rarefied flow about a delta wing are made using the direct simulation Monte Carlo (DSMC) method of Bird, and the results of the computations are compared with recent experimental data obtained in a vacuum wind tunnel at the DLR in Gottingen, Germany. The present study considers Mach 8.89 nitrogen flow for a range of conditions that include Knudsen numbers of 0.016 to 3.505 for an incidence angle of 30 deg, and angles of incidence of 15 to 60 deg for a constant Knudsen number of 0.389. The calculations provide details concerning the flowfield structure and surface quantities. Comparisons between the calculations and the available experimental measurements are made for aerodynamic and overall heat-transfer coefficients and recovery temperature. The agreement between the measured and calculated data are very good, well within the estimated measurement uncertainty. Comparisons are also made with modified Newtonian and free-molecule theories.

Celenligil, M. C.↗
NASA NTRS1991

Energetics of gas-surface interactions in transitional flows at entry velocities

The direct simulation Monte Carlo (DSMC) method has been used to calculate the molecular velocity and energy distributions of molecules striking a surface after traversing a shock layer in hypersonic transitional flow. The calculations were performed for a 1.6-m-diameter sphere at a nominal velocity for re-entry of 7.5 km/s over an altitude range of 130 to 90 km. Real gas effects and chemical reactions were included in the DSMC simulations. Results are presented for these conditions and the need for gas-surface interaction experiments is discussed.

Wilmoth, R. G.↗
NASA NTRS1991

Estimated accuracy of method of characteristics viscous plume solutions for an orbit plume induced environment prediction

A study has been performed to determine an estimate of the accuracy of rocket exhaust plume induced spacecraft environments calculated with the more widely used computational methods. Calculations of plume flowfields and plume induced environments using Method of Characteristic (MOC) based models are compared with both experimental results and direct simulation Monte Carlo (DSMC) results to determine the applicability and accuracy of Method of Characteristic models.

Smith, Sheldon D.↗
NASA NTRS1991

Radiation modeling with direct simulation Monte Carlo

Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available experimental data. The new method is compared with the current modeling technique and both techniques are compared with available data. The differences in the results are evaluated. The test case is based on an AVCO-Everett shock tube experiment, a 10-km/s standing shock wave in air at 0.1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to different flow conditions.

Carlson, Ann B.↗
NASA NTRS1991

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

Direct simulation Monte Carlo with ionization and radiation

Improvements in the modeling of radiation in low density shock waves with Direct Simulation Monte Carlo (DSMC) are the subject of this study. The physical processes which determine the amount of radiation in a shock wave were investigated and the way in which they were modeled with DSMC was evaluated. Three physical processes were identified for which an improvement in the modeling technique could result in improved radiation predictions. New physical modeling schemes are introduced in this report for the three processes. First a method for determining the electric field and its effect on the flow is introduced. Second, a two step reaction process for electron impact ionization reactions is evaluated. Finally, a new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. Each new scheme attempts to move the DSMC method toward more accurate physics or more reliance on experimental data. The new schemes are all compared to the current modeling techniques and the differences in the results are evaluated. In all cases the results agree with the available data as well as, or better than the results from earlier schemes.

Carlson, Ann B.↗
NASA NTRS1992

Shock interference prediction using direct simulation Monte Carlo

The shock interaction produced when an incident shock impinges on an inlet cowl lip is investigated for a cowl lip radius of 0.1 in. and flight conditions of Mach 15 and 35-km altitude. This problem is of interest in the design of hypersonic flight vehicles because the interference heating at the cowl lip surface from such a shock interaction is expected to be very high and the design of adequate thermal protection in this region will be challenging. Noncontinuum effects may be significant at this combination of altitude and small dimensions. If this is so, a computational solution based on the Navier-Stokes equations could significantly overpredict the surface heating on the cowl lip. Therefore, a kinetic theory solution technique, direct simulation Monte Carlo (DSMC), is used to predict the flowfield and surface heating. The details of modeling this problem and some preliminary results are presented. The same qualitative flow patterns are seen with DSMC and a continuum solution of the shock interference case, but the peak heating predicted with DSMC is somewhat lower.

Carlson, Ann B.↗
NASA NTRS1992

Decoupled predictions of radiative heating in air using a particle simulation method

A particle simulation technique with decoupled radiation is used to estimate the radiative emission along the stagnation streamline and the radiative heating at the stagnation point of a blunt-nosed vehicle entering the earth's atmosphere at hypersonic speed. The direct simulation Monte Carlo (DSMC) method is used to compute the fluid mechanics of the weakly ionized flow. The radiative emission is computed using the NEQAIR computer code with the radiation decoupled from the flowfield solution. It is concluded that the new models make it possible to decrease the predicted total radiative heating at the stagnation point of the vehicle by a factor of 15. A comparison of the DSMC approach with a continuum flow model shows that the total predicted radiative heating estimates agree within a factor of 2.

Boyd, Iain D.↗
NASA NTRS1992

Monte Carlo simulation of a near-continuum shock-shock interaction problem

A complex shock interaction is calculated with direct simulation Monte Carlo (DSMC). The calculation is performed for the near-continuum flow produced when an incident shock impinges on the bow shock of a 0.1 in. radius cowl lip for freestream conditions of approximately Mach 15 and 35 km altitude. Solutions are presented both for a full finite-rate chemistry calculation and for a case with chemical reactions suppressed. In each case, both the undisturbed flow about the cowl lip and the full shock interaction flowfields are calculated. Good agreement has been obtained between the no-chemistry simulation of the undisturbed flow and a perfect gas solution obtained with the viscous shock-layer method. Large differences in calculated surface properties when different chemical models are used demonstrate the necessity of adequately representing the chemistry when making surface property predictions. Preliminary grid refinement studies make it possible to estimate the accuracy of the solutions.

Carlson, Ann B.↗
NASA NTRS1992

Direct simulation with vibration-dissociation coupling

The majority of implementations of the Direct Simulation Monte Carlo (DSMC) method of Bird do not account for vibration-dissociation coupling. Haas and Boyd have proposed the vibrationally-favored dissociation model to accomplish this task. This model requires measurements of induction distance to determine model constants. A more general expression has been derived that does not require any experimental input. The model is used to calculate one-dimensional shock waves in nitrogen and the flow past a lunar transfer vehicle (LTV). For the conditions considered in the simulation, the influence of vibration-dissociation coupling on heat transfer in the stagnation region of the LTV can be significant.

Hash, David B.↗
NASA NTRS1992

Measurement and analysis of a small nozzle plume in vacuum

Measurements of Pitot pressure and flow angle were made in the plume of a nozzle flowing nitrogen and exhausting to a vacuum. The measurements were compared to results from a numerical simulation of the flow that was based on kinetic theory and used the direct-simulation Monte Carlo (DSMC) method. Numerical results were compared with measurements made in the plume at various axial and radial stations. Total pressure measurements were made with Pitot tubes sized for specific regions of the plume. Flow angle measurements were made with a conical probe. The measurement area for flow angle extended to 160 mm (5 exit diameters) downstream of the nozzle exit plane and radially to 60 mm (1.9 exit diameters) from the plume axis. The total pressure measurements extended 480 mm (16 exit diameters) downstream and radially to 60 mm. Comparisons of computed results from the DSMC method with measurements of flow angle displayed improved agreement with increasing distance from the exit plane. Pitot pressures computed from the DSMC method were in reasonably good agreement with experimental results over the entire measurement area.

Penko, Paul F.↗
NASA NTRS1991

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