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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1995

Effects of Chemistry on Blunt-Body Wake Structure

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using direct simulation Monte Carlo (DSMC) and Navier-Stokes calculations. Particular emphasis is given to the effects of chemistry on the near-wake structure and on the surface quantities and the comparison of the DSMC results with the Navier-Stokes calculations. The flow conditions simulated are those experienced by a space vehicle at an altitude of 85 km and a velocity of 7 km/s during Earth entry. A steady vortex forms in the near wake for these freestream conditions for both chemically reactive and nonreactive air gas models. The size (axial length) of the vortex for the reactive air calculations is 25% larger than that of the nonreactive air calculations. The forebody surface quantities are less sensitive to the chemistry than the base surface quantities. The presence of the afterbody has no effect on the forebody flow structure or the surface quantities. The comparisons of DSMC and Navier-Stokes calculations show good agreement for the wake structure and the forebody surface quantities.

Dogra, Virendra K.↗
NASA NTRS1997

Rarefied Flows of Planetary Entry Capsules

The impact of rarefaction on entry capsules and spacecraft aerothermodynamics is emphasized for various aeroassist missions. The capability of the direct simulation Monte Carlo (DSMC) method to simulate such flows is demonstrated through examples of validation studies and applications. Several space flight projects and ground-based experiments are reviewed for which rarefaction effects have significant effect on spacecraft performance or experimental measurements. This review clearly demonstrates the significant role that the DSMC method plays in characterizing such flows.

James N Moss↗
NASA NTRS1998

Mars Pathfider Rarefied Aerodynamics: Computations and Measurements

On July 4, 1997, after traveling close to 500 million km, the Pathfinder spacecraft successfully completed entry, descent, and landing at Mars. In the present paper, the focus is on the hypersonic rarefied portion of Pathfinder's atmospheric entry where the synergy of flight measurements, aerodynamic calculations, and atmospheric modeling tools are used to extract Pathfinder's attitude and the freestream density. Accuracy of the capsule aerodynamics directly impacts the inferred atmospheric properties extracted from deceleration measurements made by on-board accelerometers. The range of rarefaction considered in this study extends from the free molecular to continuum conditions and angles of attack from O to 30 deg. The aerodynamic computations are made with free-molecular and direct simulation Monte Carlo codes. The calculations show that Pathfinder is statically unstable for much of the transitional rarefied regime. Due to the relatively modest forces and the gyroscopic motion of the spacecraft, the angle of attack excursions were less than 5 deg as inferred from force measurements for the rarefied portion of entry and approached a nominal zero degree trim angle near hypersonic continuum conditions.

Moss, James N.↗
NASA NTRS1998

Aerothermodynamics of the Mars Global Surveyor Spacecraft

The aerothermodynamics characteristics of the Mars Global Surveyor spacecraft are investigated and reported. These results have been used by the Mars Global Surveyor mission planners to design the aerobraking phase of the mission. Analytical and Direct Simulation Monte Carlo computer codes were used with a detailed, three dimensional model of the spacecraft to evaluate spacecraft aerobraking characteristics for flight in free molecular and transitional flow regimes. The spacecraft is found to be aerodynamically stable in aerobraking and planned contingency configurations. Aerodynamic forces, moments, and heating are found to be highly dependent on atmospheric density. Accommodation coefficient. is seen to strongly influence drag coefficient. Transitional flow effects are found to reduce overall solar panel heating. Attitude control thruster plumes are shown to interact with the freestream, diminishing the effectiveness of the attitude control system and even leading to thrust reversal. These plume-freestream interaction effects are found to be highly dependent on freestream density.

Shane, Russell W.↗
NASA NTRS1998

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗
NASA NTRS2000

DSMC Simulations of Shock Interactions About Sharp Double Cones

This paper presents the results of a numerical study of shock interactions resulting from Mach 10 flow about sharp double cones. Computations are made by using the direct simulation Monte Carlo (DSMC) method of Bird. The sensitivity and characteristics of the interactions are examined by varying flow conditions, model size, and configuration. The range of conditions investigated includes those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel.

Moss, James N.↗
NASA NTRS2000

A 3-D Coupled CFD-DSMC Solution Method With Application to the Mars Sample Return Orbiter

A method to obtain coupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC), 3-D flow field solutions for highly blunt bodies at low incidence is presented and applied to one concept of the Mars Sample Return Orbiter vehicle as a demonstration of the technique. CFD is used to solve the high-density blunt forebody flow defining an inflow boundary condition for a DSMC solution of the afterbody wake flow. By combining the two techniques in flow regions where most applicable, the entire mixed flow field is modeled in an appropriate manner.

Glass, Christopher E.↗