Search NASA⌕ Search

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.

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

At least 73 records · Page 4

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

Monte Carlo Methodology Serves Up a Software Success

Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Source record↗

Monte Carlo simulation of radiating reentry flows

The Direct Simulation Monte Carlo (DSMC) method is applied to a radiating, hypersonic, axisymmetric flow over a blunt body in the near continuum regime. The ability of the method to predict the flowfield radiation and the radiative heating is investigated for flow over the Project Fire II configuration at 11.36 kilometers per second at an altitude of 76.42 kilometers. Two methods that differ in the manner in which they treat ionization and estimate electronic excitation are employed. The calculated results are presented and compared with both experimental data and solutions where radiation effects were not included. Differences in the results are discussed. Both methods ignore self absorption and, as a result, overpredict measured radiative heating.

Taylor, Jeff C.↗

Monte Carlo simulation of reentry flows with ionization

The Direct Simulation Monte Carlo method is applied to a rarefied, weakly ionized, hypersonic flow over a blunt axisymmetric body. An ionization model based on the concept of ambipolar diffusion is used and a model for the sheath is presented. The effects of the new modeling techniques are investigated for flow over the Project Fire II configuration at 11.37 km/s at an altitude of 84.6 km. The calculated results are presented and compared with both experimental data and solutions where ionization effects were not included. In general, the calculated results overpredict the experimental values by about 15-20 percent.

Taylor, Jeff C.↗

Magnetohydrodynamic and Aerodynamic Assessment of Ballistic Entry of a 70deg Spherecone at Mars and Venus

An electrical conductivity database for continuum flow in a CO2 atmosphere over a 70 deg spherecone was created using the Data Parallel Line Relaxation Code computational fluid dynamics software to inform development of future magnetohydrodynamic subsystems at Venus and Mars. Sixteen freestream conditions were considered at Mars with atmospheric relative velocities from 5 to 8 km/s and altitudes between 20 and 80 km. Sixteen freestream conditions were considered at Venus with atmospheric relative velocities from 9 to 12 km/s and altitudes between 85 and 115 km. Results indicate that the total electrical conductivity in the flow volume always increases as velocity increases. At low velocities, the electrical conductivity is higher at high altitudes, while at high velocities, the electrical conductivity is higher at low altitudes. Three of the 80 km altitude computational fluid dynamics solutions show good agreement with Direct Simulation Monte Carlo results. In general, computational fluid dynamics predicts thinner shocks, higher electron number density, and similar vibrational temperatures as direct simulation Monte Carlo. The magnetohydrodynamic force was calculated at both Mars and Venus. Results indicate there may not be sufficient control authority to use magnetohydrodynamics as a trajectory control mechanism at Mars without artificially increasing the electrical conductivity of the flow, but there may be appreciable control authority for a drag-modulated aerocapture at Venus.

Destiny M Fawley↗

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

Wake shield

Progress on the modeling of the flow field around a wake shield using a recently obtained code based on the Monte Carlo method is discussed. The direct simulation Monte Carlo method is a method for solving the Boltzman Equation using an approximation to the collision integral term. The collision integrand is evaluated for randomly selected values of its arguments and the summation will approach the integral for large enough samples. The collision effects may be modeled for either hard sphere or various power law potentials. The convective side of the Boltzman equation is approximated over a time step using a simple trajectory calculation of molecules as they travel through the domain of interest.

Bannister, Tommy↗

Monte Carlo simulation of entry in the Martian atmosphere

The Direct Simulation Monte Carlo method of Bird is used to investigate the characteristics of low density hypersonic flowfields for typical aerobrakes during Martian atmospheric entry. The method allows for both thermal and chemical nonequilibrium. Results are presented for a sixty-degree spherically blunt cone for various nose radii and altitudes.

Hash, David B.↗

Adaptive domain decomposition for Monte Carlo simulations on parallel processors

A method is described for performing direct simulation Monte Carlo (DSMC) calculations on parallel processors using adaptive domain decomposition to distribute the computational work load. The method has been implemented on a commercially available hypercube and benchmark results are presented which show the performance of the method relative to current supercomputers. The problems studied were simulations of equilibrium conditions in a closed, stationary box, a two-dimensional vortex flow, and the hypersonic, rarefield flow in a two-dimensional channel. For these problems, the parallel DSMC method ran 5 to 13 times faster than on a single processor of a Cray-2. The adaptive decomposition method worked well in uniformly distributing the computational work over an arbitrary number of processors and reduced the average computational time by over a factor of two in certain cases.

Wilmoth, Richard G.↗

Adaptive domain decomposition for Monte Carlo simulations on parallel processors

A method is described for performing direct simulation Monte Carlo (DSMC) calculations on parallel processors using adaptive domain decomposition to distribute the computational work load. The method has been implemented on a commercially available hypercube and benchmark results are presented which show the performance of the method relative to current supercomputers. The problems studied were simulations of equilibrium conditions in a closed, stationary box, a two-dimensional vortex flow, and the hypersonic, rarefied flow in a two-dimensional channel. For these problems, the parallel DSMC method ran 5 to 13 times faster than on a single processor of a Cray-2. The adaptive decomposition method worked well in uniformly distributing the computational work over an arbitrary number of processors and reduced the average computational time by over a factor of two in certain cases.

Wilmoth, Richard G.↗

Monte Carlo calculations of diatomic molecule gas flows including rotational mode excitation

The direct simulation Monte Carlo method was used to solve the Boltzmann equation for flows of an internally excited nonequilibrium gas, namely, of rotationally excited homonuclear diatomic nitrogen. The semi-classical transition probability model of Itikawa was investigated for its ability to simulate flow fields far from equilibrium. The behavior of diatomic nitrogen was examined for several different nonequilibrium initial states that are subjected to uniform mean flow without boundary interactions. A sample of 1000 model molecules was observed as the gas relaxed to a steady state starting from three specified initial states. The initial states considered are: (1) complete equilibrium, (2) nonequilibrium, equipartition (all rotational energy states are assigned the mean energy level obtained at equilibrium with a Boltzmann distribution at the translational temperature), and (3) nonequipartition (the mean rotational energy is different from the equilibrium mean value with respect to the translational energy states). In all cases investigated the present model satisfactorily simulated the principal features of the relaxation effects in nonequilibrium flow of diatomic molecules.

Yoshikawa, K. K.↗

Monte Carlo simulation of reentry plasmas

Attention is given to the treatment of ionization and plasma effects in the direct simulation Monte Carlo method. The requirements for accurate modeling of reentry plasmas are discussed along with the difficulties these requirements present. The current method for modeling such plasmas is reviewed and an alternative method is presented. Both methods are applied to the flow of a 10 km/s shock wave in air at 0.1 torr; a flowfield directly relevant to the projected aeroassisted orbital transfer vehicle. The results are compared and the differences between the methods are discussed.

Carlson, Ann B.↗

Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices

A study is made of the issues surrounding prediction of microchannel flows using the direct simulation Monte Carlo method. This investigation includes the introduction and use of new inflow and outflow boundary conditions suitable for subsonic flows. A series of test simulations for a moderate-size microchannel indicates that a high degree of grid under-resolution in the streamwise direction may be tolerated without loss of accuracy. In addition, the results demonstrate the importance of physically correct boundary conditions, as well as possibilities for reducing the time associated with the transient phase of a simulation. These results imply that simulations of longer ducts may be more feasible than previously envisioned.

Nance, Robert P.↗

Monte Carlo simulation of a near-continuum shock-shock interaction problem

A complex shock interaction is calculated with direct simulation Monte Carlo (DSMC). The calculation is performed for the near-continuum flow produced when an incident shock impinges on the bow shock of a 0.1 in. radius cowl lip for freestream conditions of approximately Mach 15 and 35 km altitude. Solutions are presented both for a full finite-rate chemistry calculation and for a case with chemical reactions suppressed. In each case, both the undisturbed flow about the cowl lip and the full shock interaction flowfields are calculated. Good agreement has been obtained between the no-chemistry simulation of the undisturbed flow and a perfect gas solution obtained with the viscous shock-layer method. Large differences in calculated surface properties when different chemical models are used demonstrate the necessity of adequately representing the chemistry when making surface property predictions. Preliminary grid refinement studies make it possible to estimate the accuracy of the solutions.

Carlson, Ann B.↗

Aerothermal Effects of Hypersonic Shock Impingement on a Highly Swept Leading Edge

A numerical study is conducted to investigate aerothermal implications of oblique shock impingement on a highly swept leading edge in hypersonic flow, with flow conditions and geometry relevant to non-propulsive leading edge heating on hypersonic cruise vehicles. In contrast to propulsive inlet leading edge shock impingement, which has been subject to extensive historical study, bow shock impingement on a highly swept wing leading edge has received comparatively little attention in the literature. Distinguishing characteristics of this type of flow include a potentially dominant surface heating contribution from crossflow-induced viscous dissipation, possible wall velocity slip and temperature jump effects, and various modes of shock layer thermal nonequilibrium that can interact with shock impingement flowfield structures. One overall goal of this work is to assess the impact of bow shock impingement and related nonequilibrium phenomena for wing leading edge design, materials selection and thermal management. Another goal is to improve understanding – and thereby reduce associated design margins – for the physical processes and mechanisms that may influence wing leading edge heating in the presence of an impinging shock. A combination of direct simulation Monte Carlo and Reynolds-averaged Navier-Stokes analysis is used to assess surface and flowfield characteristics for a series of hypersonic shock impingement problems, based on historical ground tests and related flight conditions. Favorable comparisons are performed with experimental data, and various thermal nonequilibrium effects are quantified.

direct simulation Monte Carlo↗

Monte Carlo simulation of an expanding gas

By application of simple computer graphics techniques, the statistical performance of two Monte Carlo methods used in the simulation of rarefied gas flows are assessed. Specifically, two direct simulation Monte Carlo (DSMC) methods developed by Bird and Nanbu are considered. The graphics techniques are found to be of great benefit in the reduction and interpretation of the large volume of data generated, thus enabling important conclusions to be drawn about the simulation results. Hence, it is discovered that the method of Nanbu suffers from increased statistical fluctuations, thereby prohibiting its use in the solution of practical problems.

Boyd, Iain D.↗

Monte Carlo simulations in support of the Shuttle upper atmospheric mass spectrometer experiment

This paper presents the results of calculations obtained with a Direct Simulation Monte Carlo (DSMC) method that describes both the external flow about the nose region of the Shuttle Orbiter and the internal flow in an inlet tube that provides the initial path between the shock-processed gases and a mass spectrometer mounted inside the Orbiter. A dedicated, three-dimensional version of the DSMC was developed during this study for the internal flow simulation along with the procedure for interfacing the external and internal flows. The calculations span the 140 to 95 km altitude range, which includes most of the hypersonic transitional flow encountered during reentry. All calculations were for a multicomponent gas mixture consisting of five chemical species while simulating the effects of transitional, rotational, vibrational, and chemical nonequilibrium. The results of the calculations show that within the entry region of the inlet tube where the gas has equilibrated with the sidewall tube temperature, the pressure is substantially less than the pressure at the external surface. This pressure correction for the entry region is significant for all conditions investigated and increases with altitude. The results highlight the structure of both the external and internal flow. Results of parametric studies show the sensitivity of the equilibrated inlet sidewall pressure to mass flow rates, gas-surface reflection model, tube side-wall temperature variations, and surface recombinations.

Moss, J. N.↗