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 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Analysis of plume backflow around a nozzle lip in a nuclear rocket

The structure of the flow around a nuclear thermal rocket nozzle lip has been investigated using the direct simulation Monte Carlo method. Special attention has been paid to the behavior of a small amount of harmful particles that may be present in the rocket exhaust gas. The harmful fission product particles are modeled by four inert gases whose molecular weights are in a range of 4 131. Atomic hydrogen, which exists in the flow due to the extremely high nuclear fuel temperature in the reactor, is also included. It is shown that the plume backflow is primarily determined by the thin subsonic fluid layer adjacent to the surface of the nozzle lip, and that the inflow boundary in the plume region has negligible effect on the backflow. It is also shown that a relatively large amount of the lighter species is scattered into the backflow region while the amount of the heavier species becomes negligible in this region due to extreme separation between the species. Results indicate that the backscattered molecules are very energetic and are fast-moving along the surface in the backflow region near the nozzle lip.

Chung, Chan H.↗

NASA NTRS1993

DSMC simulations of Mach 20 nitrogen flows about a 70 degree blunted cone and its wake

Numerical results obtained with the direct simulation Monte Carlo (DSMC) method 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 simulated with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is to characterize the wake flow under rarefied conditions. This is accomplished by calculating the influence of rarefaction on wake structure along with the impact that an afterbody has on flow features. This data report presents extensive information concerning flowfield features and surface quantities.

Moss, James N.↗

NASA NTRS1994

Implementation of a vibrationally linked chemical reaction model for DSMC

A new procedure closely linking dissociation and exchange reactions in air to the vibrational levels of the diatomic molecules has been implemented in both one- and two-dimensional versions of Direct Simulation Monte Carlo (DSMC) programs. The previous modeling of chemical reactions with DSMC was based on the continuum reaction rates for the various possible reactions. The new method is more closely related to the actual physics of dissociation and is more appropriate to the particle nature of DSMC. Two cases are presented: the relaxation to equilibrium of undissociated air initially at 10,000 K, and the axisymmetric calculation of shuttle forebody heating during reentry at 92.35 km and 7500 m/s. Although reaction rates are not used in determining the dissociations or exchange reactions, the new method produces rates which agree astonishingly well with the published rates derived from experiment. The results for gas properties and surface properties also agree well with the results produced by earlier DSMC models, equilibrium air calculations, and experiment.

Carlson, A. B.↗

NASA NTRS1992

Computation of flowfields for hypersonic flight at high altitudes

Recent direct simulation Monte Carlo methods are reviewed focusing on the nonequilibrium aspects of rarefied hypersonic flows typical of reentry conditions. The calculations identify the altitudes for which the onset of dissociation occurs, the magnitude of the surface temperature variations, the effect of the shock wave structure on the chemistry of a dissociating gas, and the radiation emisssion for reentry flowfields.

Moss, James N.↗

NASA NTRS1995

DSMC calculations for 70-deg blunted cone at 3.2 km/s in nitrogen

Numerical results obtained with the direct simulation Monte Carlo (DSMC) method are presented for Mach 15.6 nitrogen flow about a 70-deg spherically blunted cone at zero incidence. This flow condition is one of several generated in the Large Energy National Shock (LENS) tunnel during tests of a 15.24 cm diameter model with an afterbody sting. The freestream Knudsen number, based on model diameter, is 0.0023. The focus of the DSMC calculations is to characterize the near wake flow under conditions where rarefaction effects may influence afterbody aerothermal loads. This report provides information concerning computational details along with flowfield and surface quantities. Calculations show that the flow enveloping the test model is in thermal nonequilibrium and a sizable vortex develops in the near wake. Along the model baseplane the heating rates are about 0.6 percent of the forebody stagnation value while the maximum heating along the sting is about 4.2 percent of the forebody stagnation value. Comparison of a Navier-Stokes solution with the present calculations show good agreement for surface heating, pressure, and skin friction results.

Moss, J. N.↗

NASA NTRS1994

Computation of axisymmetric and ionized hypersonic 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 (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that 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 continuum of a previous study that 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 model employed.

Boyd, Iain D.↗

NASA NTRS1995

Some Developments of the Equilibrium Particle Simulation Method for the Direct Simulation of Compressible Flows

The direct simulation Monte Carlo (DSMC) method is the established technique for the simulation of rarefied gas flows. In some flows of engineering interest, such as occur for aero-braking spacecraft in the upper atmosphere, DSMC can become prohibitively expensive in CPU time because some regions of the flow, particularly on the windward side of blunt bodies, become collision dominated. As an alternative to using a hybrid DSMC and continuum gas solver (Euler or Navier-Stokes solver) this work is aimed at making the particle simulation method efficient in the high density regions of the flow. A high density, infinite collision rate limit of DSMC, the Equilibrium Particle Simulation method (EPSM) was proposed some 15 years ago. EPSM is developed here for the flow of a gas consisting of many different species of molecules and is shown to be computationally efficient (compared to DSMC) for high collision rate flows. It thus offers great potential as part of a hybrid DSMC/EPSM code which could handle flows in the transition regime between rarefied gas flows and fully continuum flows. As a first step towards this goal a pure EPSM code is described. The next step of combining DSMC and EPSM is not attempted here but should be straightforward. EPSM and DSMC are applied to Taylor-Couette flow with Kn = 0.02 and 0.0133 and S(omega) = 3). Toroidal vortices develop for both methods but some differences are found, as might be expected for the given flow conditions. EPSM appears to be less sensitive to the sequence of random numbers used in the simulation than is DSMC and may also be more dissipative. The question of the origin and the magnitude of the dissipation in EPSM is addressed. It is suggested that this analysis is also relevant to DSMC when the usual accuracy requirements on the cell size and decoupling time step are relaxed in the interests of computational efficiency.

Macrossan, M. N.↗

NASA NTRS1995

Numerical investigations in the backflow region of a vacuum plume

Four tasks were completed in this period and results were published in AIAA papers. First, a Boltzmann-2D code, was developed and applied to compute MSFC-A2 nozzle/plume flow field. It solved the two-dimensional Boltzmann-BGK equation using the Finite Difference Discrete Ordinate (FDDO) numerical technique. The code was validated by experimental data for one-dimensional shock structure predictions, paper 95-2056. Successful results for nozzle/plume flow simulation using the developed Boltzmann-2D code were presented at the 1995 AIAA Aerospace Science Conference, paper 95-0627. Second, a computer code solving two-dimensional Burnett equations was developed and applied to low-density nozzle flow field calculation. Results were also published at the 1994 AIAA Thermophysics Conference, paper 94-2055. Third, the developed two-dimensional Burnett code was extended to compute axisymmetric flow field inside MSFC-A2 nozzle, paper 95-2008. The computed nozzle exit conditions are used as input data for Direct Simulation Monte Carlo (DSMC) plume calculation. Fourth, a DSMC code was modified to compute the exhausted plume near the nozzle exit and in the backflow region.

Liaw, Goang-Shin↗

NASA NTRS1995

Comparisons of the Maxwell and CLL gas/surface interaction models using DSMC

The behavior of two different models of gas-surface interactions is studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic forces and heat transfer between the Maxwell and the Cercignani-Lampis-Lord (CLL) models for flat plate configurations at freestream conditions corresponding to a 140 km orbit around Venus. The size of the flat plate represents one of the solar panels on the Magellan spacecraft, and the freestream conditions correspond to those experienced during aerobraking maneuvers. Results are presented for both a single flat plate and a two-plate configuration as a function of angle of attack and gas-surface accommodation coefficients. The two-plate system is not representative of the Magellan geometry but is studied to explore possible experiments that might be used to differentiate between the two gas-surface interaction models. The Maxwell and CLL models produce qualitatively similar results for the aerodynamic forces and heat transfer on a single flat plate. However, the flow fields produced with the two models are qualitatively different for both the single-plate and two-plate calculations. These differences in the flowfield lead to predictions of the angle of attack for maximum heat transfer in a two plate configuration that are distinctly different for the two gas-surface interactions models.

Hedahl, Marc O.↗

NASA NTRS1994

Shuttle Upper Atmosphere Mass Spectrometer Experimental Flight Results

Calibrated pressure measurements for species with mass-to-charge ratios up to 50 amu/e(-) were obtained trom the shuttle upper atmosphere mass spectrometer experiment during re-entry on the STS-35 mission. The principal experimental objective is to obtain measurements of freestream density in the hypersonic rarefied flow flight regime. Data were collected from 180 to about 87 km. However, data above 115 km were contaminated from a source of gas emanating from pressure transdueers connected in parallel to the mass spectrometer. At lower altitudes, the pressure transducer data are compared to the mass spectrometer total pressure with excellent agreement. Near the orifice entrance, a significant amount of CO2 was generated from chemical reactions. The freestream 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 wavelike nature seen on previous flights using accelerometry. Selected spectra are presented at higher altitudes (320 km) showing the effects of the ingestion of gases from a forward fuselage fuel dump.

Blanchard, R. C.↗

NASA NTRS1996

Modeling of micro thrusters for gravity probe B

The concept of testing Einstein's general theory of relativity by means of orbiting gyroscopes was first proposed in 1959, which lead to the development of the Gravity Probe B experiment. Einstein's theory concerns the predictions of the relativistic precession of a gyroscope in orbit around earth. According to his theory, there will be two precessions due to the warping of space-time by the earth's gravitational field: the geodetic precession in the plane of the orbit, and the frame-dragging effect, in the direction of earth rotation. For a polar orbit, these components are orthogonal. In order to simplify the measurement of the precessions, Gravity Probe B (GP-B) will be placed in a circular polar orbit at 650 km, for which the predicted precessions will be 6.6 arcsec/year (geodetic) and 42 milli-arcsec/year (frame-dragging). As the gyroscope precesses, the orientation of its spin-axis will be measured with respect to the line-of-sight to Rigel, a star whose proper motion is known to be within the required accuracy. The line-of-sight to Rigel will be established using a telescope, and the orientation of the gyroscope spin axis will be measured using very sensitive SQUID (Superconducting Quantum Interference Device) magnetometers. The four gyroscopes will be coated with niobium. Below 2K, the niobium becomes superconducting and a dipole field will be generated which is precisely aligned with the gyroscope spin-axis. The change in orientation of these fields, as well as the spin-axis, is sensed by the SQUID magnetometers. In order to attain the superconducting temperatures for the gyroscopes and the SQUID's, the experiment package will be housed in a dewar filled with liquid helium. The helium flow through a GP-B micro thruster and into a vacuum is investigated using the Direct Simulation Monte Carlo method.

Jones, Kenneth M.↗

NASA NTRS1996

Assessment of thermochemical nonequilibrium and slip effects for Orbital Reentry Experiment (OREX)

Results are provided from a viscous shock layer (VSL) analysis of the reentry flowfield around the forebody of the Japanese Orbital Reentry Experiment (OREX) vehicle. This vehicle is a 50 deg. spherically blunted cone with a nose radius of 1.35 m and a base diameter of 3.4 m. Calculations are done for the OREX trajectory from 105 to 48.4 km altitude range. A 7-species chemical model is found adequate for the flowfield analysis. However, for altitudes greater than 84 km, the low density effects (such as thermal nonequilibrium and slip) are to be implemented for good agreement between the predictions and flight inferred heat-transfer rate data. Further, at altitudes lower than 84 km, a finite surface recombination probability is to be employed in place of a non-catalytic surface for better comparison between the calculations and data. VSL results are also compared with the direct simulation Monte Carlo (DSMC) predictions at high altitudes (greater than 80 km) and the electron number density data for three altitudes in the OREX trajectory. Overall, there is a good comparison between the flight data and calculated results. With the ongoing refinements in data extraction procedures, the OREX data should prove valuable for validating theoretical models employed in flowfield codes for calculation of reacting-gas flowfields.

Gupta, Roop N.↗

NASA NTRS1996

Review of blunt body wake flows at hypersonic low density conditions

Recent results of experimental and computational studies concerning hypersonic flows about blunted cones including their near wake are reviewed. Attention is focused on conditions where rarefaction effects are present, particularly in the wake. The experiments have been performed for a common model configuration (70 deg spherically-blunted cone) in five hypersonic facilities that encompass a significant range of rarefaction and nonequilibrium effects. Computational studies using direct simulation Monte Carlo (DSMC) and Navier-Stokes solvers have been applied to selected experiments performed in each of the facilities. In addition, computations have been made for typical flight conditions in both Earth and Mars atmospheres, hence more energetic flows than produced in the ground-based tests. Also, comparisons of DSMC calculations and forebody measurements made for the Japanese Orbital Reentry Experiment (OREX) vehicle (a 50 deg spherically-blunted cone) are presented to bridge the spectrum of ground to flight conditions.

Moss, J. N.↗

NASA NTRS1996

DSMC simulations of OREX entry conditions

Results of direct simulation Monte Carlo (DSMC) solutions are presented for the Japanese Orbital Reentry Experiment (OREX) vehicle, a 50 deg half-angle spherically blunted cone with a nose radius of 1.35 m and a base diameter of 3.4 m. The flow conditions simulated are those for entry into the Earth's atmosphere at a nominal velocity of about 7.4 km/s and zero incidence. Calculations are made for the higher altitude portion of entry, encompassing the transitional flow regime (altitudes of 200 to 80 km). Comparisons with flight measured values are made for axial acceleration, surface pressure, and stagnation point heating.

Moss, James N.↗

NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗

NASA NTRS1995

Comparisons of the Maxwell and CLL Gas/Surface Interaction Models Using DSMC

Two contrasting models of gas-surface interactions are studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic forces and heat transfer between the Maxwell and Cercignani-Lampis-Lord (CLL) models for flat plate configurations at freestream conditions corresponding to a 140 km orbit around Venus. The size of the flat plate is that of one of the solar panels on the Magellan spacecraft, and the freestream conditions are one of those experienced during aerobraking maneuvers. Results are presented for both a single flat plate and a two-plate configuration as a function of angle of attack and gas-surface accommodation coefficients. The two plate system is not representative of the Magellan geometry, but is studied to explore possible experiments that might be used to differentiate between the two gas surface interaction models.

Hedahl, Marc O.↗

NASA NTRS1996

Particle/Continuum Hybrid Simulation in a Parallel Computing Environment

The objective of this study was to modify an existing parallel particle code based on the direct simulation Monte Carlo (DSMC) method to include a Navier-Stokes (NS) calculation so that a hybrid solution could be developed. In carrying out this work, it was determined that the following five issues had to be addressed before extensive program development of a three dimensional capability was pursued: (1) find a set of one-sided kinetic fluxes that are fully compatible with the DSMC method, (2) develop a finite volume scheme to make use of these one-sided kinetic fluxes, (3) make use of the one-sided kinetic fluxes together with DSMC type boundary conditions at a material surface so that velocity slip and temperature slip arise naturally for near-continuum conditions, (4) find a suitable sampling scheme so that the values of the one-sided fluxes predicted by the NS solution at an interface between the two domains can be converted into the correct distribution of particles to be introduced into the DSMC domain, (5) carry out a suitable number of tests to confirm that the developed concepts are valid, individually and in concert for a hybrid scheme.

Baganoff, Donald↗

NASA NTRS1995

NAS Experiences of Porting CM Fortran Codes to HPF on IBM SP2 and SGI Power Challenge

Current Connection Machine (CM) Fortran codes developed for the CM-2 and the CM-5 represent an important class of parallel applications. Several users have employed CM Fortran codes in production mode on the CM-2 and the CM-5 for the last five to six years, constituting a heavy investment in terms of cost and time. With Thinking Machines Corporation's decision to withdraw from the hardware business and with the decommissioning of many CM-2 and CM-5 machines, the best way to protect the substantial investment in CM Fortran codes is to port the codes to High Performance Fortran (HPF) on highly parallel systems. HPF is very similar to CM Fortran and thus represents a natural transition. Conversion issues involved in porting CM Fortran codes on the CM-5 to HPF are presented. In particular, the differences between data distribution directives and the CM Fortran Utility Routines Library, as well as the equivalent functionality in the HPF Library are discussed. Several CM Fortran codes (Cannon algorithm for matrix-matrix multiplication, Linear solver Ax=b, 1-D convolution for 2-D datasets, Laplace's Equation solver, and Direct Simulation Monte Carlo (DSMC) codes have been ported to Subset HPF on the IBM SP2 and the SGI Power Challenge. Speedup ratios versus number of processors for the Linear solver and DSMC code are presented.

Saini, Subhash↗