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 235 records · Page 13

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↗

Application of a parallel DSMC method to hypersonic rarefied flows

This paper describes a method for doing direct simulation Monte Carlo (DSMC) calculations using parallel processing and presents some results of applying the method to several hypersonic, rarefied flow problems. The performance and efficiency of the parallel method are discussed. The applications described are the flow in a channel and the flow about a flat plate at incidence. The results show significant advantages of parallel processing over conventional scalar processing and demonstrate the scalability of the method to large problems.

Wilmoth, Richard G.↗

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

DSMC calculations for the double ellipse

The direct simulation Monte Carlo (DSMC) method involves the simultaneous computation of the trajectories of thousands of simulated molecules in simulated physical space. Rarefied flow about the double ellipse for test case 6.4.1 has been calculated with the DSMC method of Bird. The gas is assumed to be nonreacting nitrogen flowing at a 30 degree incidence with respect to the body axis, and for the surface boundary conditions, the wall is assumed to be diffuse with full thermal accommodation and at a constant wall temperature of 620 K. A parametric study is presented that considers the effect of variations of computational domain, gas model, cell size, and freestream density on surface quantities.

Moss, James N.↗

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↗

A shock-layer theory based on thirteen-moment equations and DSMC calculations of rarefied hypersonic flows

Grad's thirteen-moment equations are applied to the flow behind a bow shock under the formalism of a thin shock layer. Comparison of this version of the theory with Direct Simulation Monte Carlo calculations of flows about a flat plate at finite attack angle has lent support to the approach as a useful extension of the continuum model for studying translational nonequilibrium in the shock layer. This paper reassesses the physical basis and limitations of the development with additional calculations and comparisons. The streamline correlation principle, which allows transformation of the 13-moment based system to one based on the Navier-Stokes equations, is extended to a three-dimensional formulation. The development yields a strip theory for planar lifting surfaces at finite incidences. Examples reveal that the lift-to-drag ratio is little influenced by planform geometry and varies with altitudes according to a 'bridging function' determined by correlated two-dimensional calculations.

Cheng, H. K.↗

Three-dimensional rarefied flow simulations for the Aeroassist Flight Experiment vehicle

Hypersonic rarefied flow about the Aeroassist Flight Experiment vehicle has been investigated using a three-dimensional direct simulation Monte Carlo method. Calculations are performed for the transitional flows encountered during the vehicle's atmospheric entry for altitudes of 110 and 100 km with an entry velocity of 9.9 km/s. The simulations are performed using a five-species reacting gas model that account for rotational and vibrational internal energies. The solutions indicate that dissociation is important at altitudes of 110 km and below. Results are presented for surface pressures, convective heating, flowfield structure, and aerodynamic coefficient variations with altitude.

Celenligil, M. Cevdet↗

DSMC calculations for the delta wing

Results are reported from three-dimensional direct simulation Monte Carlo (DSMC) computations, using a variable-hard-sphere molecular model, of hypersonic flow on a delta wing. The body-fitted grid is made up of deformed hexahedral cells divided into six tetrahedral subcells with well defined triangular faces; the simulation is carried out for 9000 time steps using 150,000 molecules. The uniform freestream conditions include M = 20.2, T = 13.32 K, rho = 0.00001729 kg/cu m, and T(wall) = 620 K, corresponding to lambda = 0.00153 m and Re = 14,000. The results are presented in graphs and briefly discussed. It is found that, as the flow expands supersonically around the leading edge, an attached leeside flow develops around the wing, and the near-surface density distribution has a maximum downstream from the stagnation point. Coefficients calculated include C(H) = 0.067, C(DP) = 0.178, C(DF) = 0.110, C(L) = 0.714, and C(D) = 1.089. The calculations required 56 h of CPU time on the NASA Langley Voyager CRAY-2 supercomputer.

Celenligil, M. Cevdet↗

Flat plate at incidence as a waverider in rarefied hypersonic flow

The physical validity of continuum models and their ability to predict the critical aerothermodynamic properties of a waverider at high altitudes are examined using a flat plate at angle of attack as a generic hypersonic lifting vehicle. For a shock layer far from local translational equilibrium, a theoretical study based on Grad's thirteen-moment equations shows that the Navier-Stokes based solutions can correctly predict the drag, lift, and surface heat transfer rate, with the prediction error comparable to that of the standard shock-layer theory. The conclusion is supported by a comparison with direct simulation Monte Carlo calculations.

Cheng, H. K.↗

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

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

Hypersonic rarefied flow about a compression corner - DSMC simulation and experiment

The phenomena of shock/boundary-layer interactions and flow separation are investigated using both computational and experimental methods for low-density hypersonic flow about two-dimensional compression corners. The numerical calculations are made with the direct simulation Monte Carlo (DSMC) method. Experimental measurements provide information concerning the flowfield structure and surface flow patterns by means of gas glow discharge and oil flow pictures, respectively. Comparison of the two data sets provides a qualitative basis for assessing the ability of the DSMC method to describe such flows.

Moss, James N.↗