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 199 records · Page 11

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

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

An efficient DSMC algorithm applied to a delta wing

A new algorithm for 3D direct simulation Monte Carlo (DSMC) is tested and numerical results are compared with wind tunnel data and results obtained earlier with a more traditional DSMC code. The test case is the flowfield around a delta wing at incidence at Knudsen number of 0.016 and Mach number of 20.2. The results are shown to compare favorably with both experimental and earlier numerical results. The new algorithm is described with special emphasis placed on its distinctive features: Cartesian/unstructured combination grid, special body surface definition, discretization in physical space.

Rault, D. F. G.↗

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

New approach in direct-simulation of gas mixtures

Results are reported for an investigation of a new direct-simulation Monte Carlo method by which energy transfer and chemical reactions are calculated. The new method, which reduces to the variable cross-section hard sphere model as a special case, allows different viscosity-temperature exponents for each species in a gas mixture when combined with a modified Larsen-Borgnakke phenomenological model. This removes the most serious limitation of the usefulness of the model for engineering simulations. The necessary kinetic theory for the application of the new method to mixtures of monatomic or polyatomic gases is presented, including gas mixtures involving chemical reactions. Calculations are made for the relaxation of a diatomic gas mixture, a plane shock wave in a gas mixture, and a chemically reacting gas flow along the stagnation streamline in front of a hypersonic vehicle. Calculated results show that the introduction of different molecular interactions for each species in a gas mixture produces significant differences in comparison with a common molecular interaction for all species in the mixture. This effect should not be neglected for accurate DSMC simulations in an engineering context.

Chung, Chan-Hong↗

Numerical and experimental investigations of rarefied nozzle and plume flows of nitrogen

Numerical and experimental investigations are performed for the rarefied flow of nitrogen through a small nozzle which is expanded into near-vacuum conditions. Two different numerical studies are undertaken: the first employs a continuum approach in solving the Navier-Stokes equations, and the second employs a particle approach through use of the direct simulation Monte Carlo method (DSMC). The experimental investigation concerns the measurement of pressure, using a Pitot tube, in the nozzle exit plane and near-field of the plume. Comparison of the experimental and numerical data at the nozzle exit reveals that the DSMC technique provides the more accurate description of the expanding flow. It is discovered that the DSMC solutions are quite sensitive to the model employed to simulate the interaction between the gas and the nozzle wall surface. It is concluded that the simplistic fully diffuse model is quite satisfactory for the present application.

Boyd, Iain D.↗

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

Numerical analysis of rarefied slit flows. II - Navier-Stokes simulations

The model problem of pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is analyzed via full Navier-Stokes numerical simulation. Parametric solutions are generated for slit-height based Knudsen number ranging from continuum to transitional flow and for reservoir pressure ratios leading to subsonic and supersonic flow. The change in the structure of the flowfield near the slit as a function of pressure ratio and Knudsen number are quantified from a purely continuum standpoint. The choice of numerical domain size, boundary conditions and treatment of the slit are also discussed. As expected, comparison with a Direct Simulation Monte Carlo solution for a highly rarefied case shows large differences in the predicted mass flow. The cause of these differences can be quantified through detailed comparison of the local flowfield properties. For the larger pressure ratio cases qualitative trends with increasing rarefaction are discussed, including the change in the sonic line shape in the slit and the in total mass flow.

Wadsworth, D. C.↗