SEARCH · Engineering Papers

Results for “direct simulation Monte Carlo methods”

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 109 records · Page 6

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 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 NTRS1994

Direct simulation Monte Carlo prediction of on-orbit contaminant deposit levels for HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flow field and surface conditions and geometric orientations for the satellite in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. A detailed description of the adaptation of this solution method to the study of the satellite's environment is also presented. Results pertaining to the satellite's environment are presented regarding contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface, along with data related to code performance. Using procedures developed in standard contamination analyses, along with many worst-case assumptions, the cumulative upper-limit level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated at about 13,350 A.

Woronowicz, Michael S.↗

NASA NTRS1992

Procedure for Adapting Direct Simulation Monte Carlo Meshes

A technique is presented for adapting computational meshes used in the G2 version of the direct simulation Monte Carlo method. The physical ideas underlying the technique are discussed, and adaptation formulas are developed for use on solutions generated from an initial mesh. The effect of statistical scatter on adaptation is addressed, and results demonstrate the ability of this technique to achieve more accurate results without increasing necessary computational resources.

Woronowicz, Michael S.↗

NASA NTRS1993

Direct Monte Carlo Simulations of Hypersonic Viscous Interactions Including Separation

Results of calculations obtained using the direct simulation Monte Carlo method for Mach 25 flow over a control surface are presented. The numerical simulations are for a 35-deg compression ramp at a low-density wind-tunnel test condition. Calculations obtained using both two- and three-dimensional solutions are reviewed, and a qualitative comparison is made with the oil flow pictures highlight separation and three-dimensional flow structure.

Moss, James N.↗

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 NTRS1998

Magellan Aerodynamic Characteristics During the Termination Experiment Including Thruster Plume-Free Stream Interaction

Results are presented on the aerodynamic characteristics of the Magellan spacecraft during the October 1994 Termination Experiment, including the effects of the thruster engine exhaust plumes upon the molecular free stream around the spacecraft and upon the aerodynamics coefficients. As Magellan passed through the Venusian atmosphere, the solar arrays were turned in opposite directions relative to the free stream creating a torque on the spacecraft. The spacecraft control system was programmed to counter the effects of this torque with attitude control engines to maintain an inertially fixed attitude. The orientation and reaction engine telemetry returned from Magellan are used to create a model of the aerodynamic torques. Geometric models of the Magellan spacecraft are analyzed with the aid of both free molecular and Direct Simulation Monte Carlo codes. The simulated aerodynamic torques determined are compared to the measured torques. The Direct Simulation Monte Carlo method is also used to model the attitude engine exhaust plumes, the free stream disturbance caused by these plumes, and the resulting torques acting on the spacecraft compared to no-exhaust plume cases. The effect of the exhaust plumes was found to be sufficiently large that thrust reversal is possible.

Cestero, Francisco J.↗

NASA NTRS1999

Fluid/Solid Boundary Conditions in Non-Isothermal Systems

The existing theoretical research concerned with thermal creep at fluid/solid interfaces is briefly reviewed, and the importance of microgravity-based experimental data is then discussed. It is noted that the ultimate goal of this research is a rational molecular level theory that predicts the dependence of a dimensionless thermal creep coefficient, Ctc, on relevant dimensionless parameters describing the way fluid molecules interact with the solid surface and how they interact among themselves. The discussion covers thermophoresis of isolated solid spheres and aggregates in gases; solid sphere thermophoresis in liquids and dense vapors; thermophoresis of small immiscible liquid droplets; and applications of the direct simulation Monte Carlo method.

Rosner, Daniel E.↗

NASA NTRS1999

Characterization of Showerhead Performance At Low Pressure

The overall objective of this work is to characterize the flow to rough showerheads by deriving pressure drop versus velocity correlations to at can be then used in reactor scale simulations where the showerhead is approximated as a porous medium. At relatively low Reynolds numbers (less than 1-10 based on the hole length scale) and in the absence of slip flow, Darcy's Law, grad P = mu U/k, can be used to express the relation between the pressure drop and velocity where @mu@ is the fluid viscosity and it is the permeability that can be theoretically predicted as k= e R^2 /8, where e is the porosity. However, at sufficiently small hole diameters and decreased pressures (less than 5 Torr), the Knudsen number based on showerhead tube radius increases, and the flow may be in a transition regime. Different expressions have been proposed to account for this effect in the permeability by expressing k as a function of either pressure or Knudsen number. But at even higher Knudsen numbers, the pressure drop - velocity dependence is non-linear, and Darcy's Law no longer holds such that a permeability cannot be defined. The direct simulation Monte Carlo method is used along side conventional CFD techniques to determine the extent to which the CFD technique is appropriate and helps to derive correlations for the more rarefied cases of interest in these showerhead flows.

Hash, David B.↗

NASA NTRS1997

A DSMC Study of Low Pressure Argon Discharge

Work toward a self-consistent plasma simulation using the DSMC (Direct Simulation Monte Carlo) method for examination of the flowfields of low-pressure high density plasma reactors is presented. Presently, DSMC simulations for these applications involve either treating the electrons as a fluid or imposing experimentally determined values for the electron number density profile. In either approach, the electrons themselves are not physically simulated. Self-consistent plasma DSMC simulations have been conducted for aerospace applications but at a severe computational cost due in part to the scalar architectures on which the codes were employed. The present work attempts to conduct such simulations at a more reasonable cost using a plasma version of the object-oriented parallel Cornell DSMC code, MONACO, on an IBM SP-2. Due to availability of experimental data, the GEC reference cell is chosen to conduct preliminary investigations. An argon discharge is chosen to conduct preliminary investigations. An argon discharge is examined thus affording a simple chemistry set with eight gas-phase reactions and five species: Ar, Ar(+), Ar(*), Ar(sub 2), and e where Ar(*) is a metastable.

Hash, David B.↗

NASA NTRS1997

DSMC Simulations of Blunt Body Flows for Mars Entries: Mars Pathfinder and Mars Microprobe Capsules

The hypersonic transitional flow aerodynamics of the Mars Pathfinder and Mars Microprobe capsules are simulated with the direct simulation Monte Carlo method. Calculations of axial, normal, and static pitching coefficients were obtained over an angle of attack range comparable to actual flight requirements. Comparisons are made with modified Newtonian and free-molecular-flow calculations. Aerothermal results were also obtained for zero incidence entry conditions.

Moss, James N.↗

NASA NTRS2003

Direct Simulation Monte Carlo Calculations in Support of the Columbia Shuttle Orbiter Accident Investigation

The Direct Simulation Monte Carlo method was used to provide 3-D simulations of the early entry phase of the Shuttle Orbiter. Undamaged and damaged scenarios were modeled to provide calibration points for engineering "bridging function" type of analysis. Currently the simulation technology (software and hardware) are mature enough to allow realistic simulations of three dimensional vehicles.

Gallis, Michael A.↗

NASA NTRS2011

Collisional Evolution of the Enceladus Neutral Cloud

Water vapor ejected from Saturn's small moon Enceladus easily escapes its meager gravity to form a Saturn-encircling cloud with a low collision rate. Observations show that the cloud is quite broad in the radial direction, and we show here that collisions, though quite rare, may be largely responsible for this radial spreading. We modeled this cloud using the Direct Simulation Monte Carlo method, as fluid methods would be inappropriate for such a tenuous gas.

magnetic fields↗

NASA NTRS2016

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Direct simulation Monte Carlo analysis on parallel processors

A method is presented for executing a direct simulation Monte Carlo (DSMC) analysis using parallel processing. The method is based on using domain decomposition to distribute the work load among multiple processors, and the DSMC analysis is performed completely in parallel. Message passing is used to transfer molecules between processors and to provide the synchronization necessary for the correct physical simulation. Benchmark problems are described for testing the method and results are presented which demonstrate the performance on two commercially available multicomputers. The results show that reasonable parallel speedup and efficiency can be obtained if the problem is properly sized to the number of processors. It is projected that with a massively parallel system, performance exceeding that of current supercomputers is possible.

Wilmoth, Richard G.↗

NASA NTRS1987

Direct simulation of high-vorticity gas flows

The computational limitations associated with the molecular dynamics (MD) method and the direct simulation Monte Carlo (DSMC) method are reviewed in the context of the computation of dilute gas flows with high vorticity. It is concluded that the MD method is generally limited to the dense gas case in which the molecular diameter is one-tenth or more of the mean free path. It is shown that the cell size in DSMC calculations should be small in comparison with the mean free path, and that this may be facilitated by a new subcell procedure for the selection of collision partners.

Bird, G. A.↗