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

Shuttle entry aerothermodynamic flight research - The Orbiter Experiments (OEX) Program

Results of the OEX program are summarized with emphasis on the information on entry aerothermodynamic phenomena derived from Space Shuttle operations. The discussion focuses on OEX experiment complement and operational history, freestream environment and vehicle attitude data, aerodynamic force and moment data, aerodynamic surface data, and vehicle configuration data. Attention is also given to orbiter aerodynamic performance, stability and control, high-altitude atmospheric density variability, direct simulation Monte Carlo validation, orbital drag variation, and computational fluid dynamic technique validation.

Throckmorton, David A.↗

Hypersonic rarefied wake characterization

Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for hypersonic rarefied flow over an aeroassisted space transfer vehicle (ASTV). The emphasis of the study is the characterization of the near wake region which includes the ASTV payload. The study covered the transitional flow regime from near continuum to free molecular. Calculations show that the character of the near wake is significantly affected by the presence of the payload. Flow separation occurs when an afterbody is present throughout the transitional flow regime. In contrast, when no afterbody is present, no separation is observed until the flow approaches continuum.

Brewer, E. B.↗

SUMS experiment flight results on STS-35

Calibrated pressure measurements for species with mass to charge ratios up to 50 amu/e(-) were obtained from the Shuttle Upper Atmosphere Mass Spectrometer (SUMS) experiment during reentry on the STS-35 mission. Data were collected from 180 km, when the signal rose above the background, to about 87 km, when the SUMS system automatically closed the gas inlet value. However, data above 115 km was contaminated from a source of gas emanating from pressure transducers connected in parallel to the mass spectrometer. At lower altitudes, the pressure transducer data is compared with the mass spectrometer total pressure with excellent agreement. The free-stream density in the rarefied flow flight regime is calculated using an orifice pressure coefficient model based upon direct simulation Monte Carlo results. This density, when compared with the 1976 U.S. standard atmosphere model, exhibits the wave-like nature seen on previous flights using accelerometry. In addition, selected spectra are presented at higher altitudes (320 km) showing the effects of the ingestion of gases from a forward fuselage fuel dump. An analysis of the spectra data from this event is presented to show that no significant permanent changes occurred which affected the data interpretation at lower altitudes. Further, the localized chemistry from the individual species during the onset of aerodynamic heating is examined to the extent possible for a closed source system, such as SUMS. Near the orifice entrance, a significant amount of CO2 was generated from chemical reactions with the carbon panels of the Orbiter and absorbed oxygen on the system tubing.

Blanchard, Robert C.↗

Development and application of computational aerothermodynamics flowfield computer codes

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second-approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered. The cases are governed, respectively, by the following: vibrational relaxation; weak dissociation; strong dissociation; and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Venkatapathy, Ethiraj↗

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

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

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

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

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

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

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↗

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

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

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

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

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

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

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