SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Experimentation in the low-density plume of a simulated electrothermal thruster for computer code validation

Pressures and flow angles are measured in the plume of a 20 deg half-angle, conical nozzle in vacuum with Pitot tubes and conical probes. The area of measurement in the plume ranges from the nozzle exit plane to 480 mm axially downstream and from the plume centerline to 60 mm radially. The nozzle has an exit-to-throat area ratio of 100:1 and a throat diameter of 3.2 mm. The nozzle flow exhausts to a vacuum of order 10(exp -2) Pa to simulate a resistojet (an electrothermal rocket of less than 1 N of thrust) operating in space. Experimental data are given for flows of nitrogen at 55 and 68 mg/s, stagnation temperatures between 695 and 921 K, and stagnation pressures ranging from 5600 to 7100 Pa. Data are also given for argon at a rate of 68 mg/s, a stagnation temperature of 648 K, and stagnation pressures of 4500, 4750, and 4770 Pa. Measurements in the nitrogen plume are compared with computational results from a direct-simulation Monte Carlo method.

Meissner, Dana L.↗

NASA NTRS1993

Measurement and analysis of a small nozzle plume in vacuum

Pitot pressures and flow angles are measured in the plume of a nozzle flowing nitrogen and exhausting to a vacuum. Total pressures are measured with Pitot tubes sized for specific regions of the plume and flow angles measured with a conical probe. The measurement area for total pressure extends 480 mm (16 exit diameters) downstream of the nozzle exit plane and radially to 60 mm (1.9 exit diameters) off the plume axis. The measurement area for flow angle extends to 160 mm (5 exit diameters) downstream and radially to 60 mm. The measurements are compared to results from a numerical simulation of the flow that is based on kinetic theory and uses the direct-simulation Monte Carlo (DSMC) method. Comparisons of computed results from the DSMC method with measurements of flow angle display good agreement in the far-field of the plume and improve with increasing distance from the exit plane. Pitot pressures computed from the DSMC method are in reasonably good agreement with experimental results over the entire measurement area.

Penko, P. F.↗

NASA NTRS1992

Parallel implementation of a particle simulation for modeling rarefied gas dynamic flow

When the conditions of flow are rarefied and hypersonic, a more suitable alternative to the use of the Navier-Stokes equations for developing a numerical solution is the Direct Simulation Monte Carlo method (DSMC), a method of simulation which employs a large number of particles in modeling a rarefied gas. The performance of a parallel DSMC code developed for the Intel iPSC/860 Touchstone Gamma prototype computer is studied and the scaleup is found to be very nearly over the range of 16-128 processors.

Fallavollita, M. A.↗

NASA NTRS1993

Near wake structure for a generic ASTV configuration

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using the direct simulation Monte Carlo method. Particular emphasis is given to the near wake flow and its sensitivity to rarefaction and other parametric variations. The flow conditions simulated are attainable in existing low-density hypersonic wind tunnels; that is, Mach 20 nitrogen flow encompassing freestream Knudsen numbers of 0.03 to 0.001. A stable vortex forms in the near wake at and below a freestream Knudsen number of 0.01 and the size of the vortex increases with decreasing freestream Knudsen number. The base region of the flow remains in thermal nonequilibrium for all cases. There is no formation of a lip separation shock or a distinct wake shock at these rarefied conditions.

Dogra, Virendra K.↗

NASA NTRS1993

Temperature dependence of rotational relaxation in shock waves of nitrogen

Computations are carried out for one-dimensional shock waves of diatomic nitrogen, using the direct simulation Monte Carlo method. It was found that the reciprocal shock thickness varies with the upstream temperature condition. This variation was also observed in the experimental data and was simulated numerically by using the temperature-dependent expression for the rotational collision number.

Boyd, I. D.↗

NASA NTRS1993

Return flux contamination of an outgassing spacecraft in low earth orbit

The recontamination of a spacecraft surface directly exposed to the ram flux is investigated. It is assumed that the surface has been contaminated by a rocket plume contact and that, subsequent to the exposure, the contaminated surface is placed in the ram flux. A BGK-simulation was compared to both a Direct Simulation Monte Carlo (DSMC) and a Full Flow Monte Carlo (FFMC) using sensitivity analysis. Results indicate that the BGK-based simulation underpredicted the return flux contamination as compared to both the DSMC and FFMC techniques. Results from the full flow model suggest that the BGK model is inadequate to represent this physical problem.

Justiz, Charles R.↗

NASA NTRS1993

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗

NASA NTRS1993

Computation of axisymmetric and ionized flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo method (DSMC), and the continuum approach solves the Euler equations. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code which is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuation of a previous study which considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry models employed.

Boyd, Iain D.↗

NASA NTRS1993

Stagnation flowfield analysis for an aeroassisted vehicle

Hypersonic reactive flowfield results are presented for the flow about an axisymmetric representation of a three-dimensional aeroassisted vehicle using the viscous shock-layer (VSL) and Navier-Stokes (NS) equations for low- and high-altitude entry conditions, respectively. Solutions are limited to the stagnation region. The NS results have been obtained by employing the surface recombination boundary condition for an 11-species air model, whereas the VSL calculations are carried out by using a 5-species model. Three widely used surface recombination rate models are implemented with these two flowfield calculation procedures. The low-density (or high altitude) NS solutions are obtained with surface slip equations and compare quite well with the direct simulation Monte Carlo predictions. For higher densities (or lower altitudes), the no-slip VSL solutions are considered adequate. Results from this study highlight the effect of surface catalyticity on surface heating and electron number density.

Gupta, Roop N.↗

NASA NTRS1993

Kinetic relaxation of a non-Maxwellian monatomic gas in a state of gross rest

The kinetic relaxation of a monatomic gas in a state of gross rest is studied by use of Maxwell's second-order moment equation. The molecular models considered range from the hard sphere to the Maxwell molecule. The development is exact for an ellipsoidal distribution and becomes a good approximation for more general distributions. The results show that the relaxation of second-order moments is nearly exponential for power-law molecules and that the characteristic time mu/p appears as the principal controlling parameter while the power-law constant plays a secondary role. Comparisons between results of numerical simulations using Bird's direct simulation Monte Carlo method and theory show excellent agreement for both hard-sphere and Maxwell molecules.

Baganoff, D.↗

NASA NTRS1993

Higher-order viscous shock-layer solutions for high altitude flows

A higher-order viscous shock-layer (HVSL) method has been developed and is used to obtain physically consistent results under varying degrees of low-density conditions for perfect-gas and nonequilibrium flows past long slender bodies. This method is implemented with higher-order body and shock slip conditions. Detailed comparisons with the direct simulation Monte Carlo (DSMC) method and Navier-Stokes (NS) calculations establish HVSL as a computationally efficient code for the low-density calculations. Further, it is shown that the deficiencies in the standard viscous shock-layer (VSL) method when predicting low-density flows lie in the governing equations and not with the slip boundary conditions as considered by earlier researchers.

Gupta, Roop N.↗

NASA NTRS1993

Zonally-decoupled DSMC solutions of hypersonic blunt body wake flows

Direct simulation Monte Carlo (DSMC) solutions are presented for the hypersonic flow behind a blunt body in which the wake region is solved in a zonally-decoupled manner. The forebody flow is solved separately using either a DSMC or a Navier-Stokes method, and the forebody exit plane solution is specified as the inflow condition to the decoupled DSMC solution of the wake region. Results are presented for a 70-deg, blunted cone at flow conditions that can be accommodated in existing low-density wind tunnels with the Knudsen number based on base diameter ranging from 0.03 to 0.001. The zonally-decoupled solutions show good agreement with fully-coupled DSMC solutions of the wake flow densities and velocities. The wake closure predicted by the zonally-decoupled solutions is in better agreement with fully-coupled results than that predicted by a fully-coupled Navier-Stokes method indicating the need to account for rarefaction in the wake for the cases considered. The combined use of Navier-Stokes for the forebody with a decoupled DSMC solution for the wake provides an efficient method for solving transitional blunt-body flows where the forebody flow is continuum and the wake is rarefied.

Wilmoth, Richard G.↗

NASA NTRS1993

Models for dynamic surface temperatures during rarefied aeropass maneuvers

New models for calculating dynamic surface temperatures in a direct simulation Monte Carlo particle method during atmospheric entry are described. The models take into account the radiative emissivity, thermal conductivity, and heat capacity of the material. They are applied to thin panels in a highly rarefied flows and yield the transient thermal response of individual surface facets. Simulation results for radiation panels with internal conductivity are verified by comparing them with free molecule theory.

Haas, Brian L.↗

NASA NTRS1993

Flow resolution and domain of influence in rarefied hypersonic blunt-body flows

The study assesses the effects of the upstream domain size and grid resolution upon flow properties and body aerodynamics computed for rarefied flows over cold blunt bodies with a direct simulation Monte Carlo (DSMC) particle method. Empirical correlations are suggested for aerodynamic coefficients for two-dimensional flows past a perpendicular flat plate. Free-stream parameters which were varied in the study include the Mach number, Knudsen number, surface temperature, and intermolecular potential. Insufficient grid resolution leads to overprediction of aerodynamic heating and forces in the DSMC method. Solution accuracy correlates well with the Reynolds number defined at the wall temperature and the stagnation mean free path relative to the cell dimension. Insufficient upstream domain size in the DSMC method leads to overprediction of heating and drag. Errors in aerodynamic coefficients correlate well with the distance ahead of the body where flow temperature reaches half of its peak value. Simulation of a hard-sphere gas is more sensitive to grid resolution, while simulation of a Maxwell gas is more sensitive to upstream domain size.

Haas, Brian L.↗

NASA NTRS1993

Hypersonic blunt body wake computations using DSMC and Navier-Stokes solvers

Numerical results obtained with direct simulation Monte Carlo (DSMC) and Navier-Stokes methods are presented for Mach 20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how does the wake structure change as a function of rarefaction, what are the afterbody levels of heating, and to what limits are continuum models realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗

NASA NTRS1993

Aerodynamics of Shuttle Orbiter at high altitudes

The high-altitude/high-Knudsen number aerodynamics of the Shuttle Orbiter are computed from Low-Earth Orbit down to 100 km using three-dimensional direct simulation Monte Carlo and free molecule codes. Results are compared with Blanchard's latest Shuttle aerodynamic model, which is based on in-flight accelerometer measurements, and bridging formula models. Good comparison is observed, except for the normal force and pitching moment coefficients. The present results were obtained for a generic Shuttle geometry configuration corresponding to a zero deflection for all control surfaces.

Rault, Didier F. G.↗

NASA NTRS1993

DSMC simulation of the interaction between rarefied free jets

This paper presents a direct simulation Monte Carlo (DSMC) calculation of two interacting free jets exhausting into vacuum. The computed flow field is compared against available experimental data and shows excellent agreement everywhere except in the very near field (less than one orifice diameter downstream of the jet exhaust plane). The lack of agreement in this region is attributed to having assumed an inviscid boundary condition for the orifice lip. The results serve both to validate the DSMC code for a very complex, three dimensional non-equilibrium flow field, and to provide some insight as to the complicated nature of this flow.

Dagum, Leonardo↗

NASA NTRS1993

Simulated rarefied aerodynamics of the Magellan spacecraft during aerobraking

Aerodynamic loads upon the Magellan spacecraft during aerobraking through the atmosphere of Venus are computed at off-design attitudes with a direct simulation Monte Carlo (DSMC) particle method. Simulated rarefied flows at nominal altitudes near 140 km and an entry speed of 8.6 km/s were compared to simulated and analytic free molecular results. Aerodynamic moments, forces, and heating for rarefied entry at all attitudes were 7-10 percent below free molecular results. All moments acted to restore the vehicle to its nominal zero-pitch, zero-yaw attitude. Suggested canting of the solar panels is an innovative configuration to assess gas-surface interaction during aerobraking. The resulting roll torques about the central body-axis as predicted in rarefied flow simulations were nearly twice that predicted for free molecular flow, although differences became less distinct for thermal accommodation coefficients well below unity. Roll torques increased dramatically with reduced accommodation coefficients employed in the simulation. In the DSMC code, periodic free-molecule boundary conditions and a coarse computational grid and body resolution served to minimize the simulation size and cost while retaining solution validity.

Haas, Brian L.↗