Engineering Papers⌕ Search

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.

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

At least 73 records · Page 4

Direct simulation Monte Carlo technique for modeling of the environment in the vicinity of the Space Shuttle Orbiter

A program for analyzing the flowfield parameters in the neighborhood of the Space Shuttle Orbiter has been developed. The program uses the direct simulation Monte Carlo method, which is a completely probabilistic Monte Carlo technique capable of analyzing 3-dimensional steady or unsteady flow with prescribed internal and external boundary conditions. The freestream flux densities incident on the external flowfield boundaries are calculated from the drifting Maxwellian gas properties of the freestream. The flux entering the flowfield from the internal boundary is calculated from the outgassing flux density distribution over the Shuttle external surface and the prescribed discrete source fluxes. This technique produces a numerical flowfield solution which is the probabilistic equivalent of a complete solution of the time-dependent, 3-dimensional Boltzmann equation. Flowfield results are presented for the following configurations: (1) Shuttle angles of attack of 0 and 90 deg, (2) freestream density 10 to the 9th - 10 to the 11th per cu cm, (3) Shuttle outgassing flux density, and (4) operation of the aft downfiring vernier Reaction Control System (RCS) engine. Results are presented for column density distribution of outgassed and engine species and for the flux density of outgassed and engine species incident on the Shuttle bay and the windshield.

Melfi, L. T., Jr.↗

Direct simulation of typical AOTV entry flows

The direct simulation Monte Carlo method is applied to two flows that are relevant to the projected aeroassisted orbital transfer vehicle. The first is the axially symmetric flow about a 70 degree blunted cone shield with a cylindrical afterbody, and the second is the one-dimensional flow along the stagnation streamline of the first flow. The shock standoff distance for the stagnation streamline study is obtained from the full axially symmetric calculation. Advantage is then taken of the speed of the one-dimensional calculation to study the sensitivity of the stagnation point heat transfer to a number of physical variables. These include the chemical reaction rate coefficients, the surface catalycity, the electric field associated with ambipolar diffusion, and the nonequilibrium thermal radiation from the shock front.

Bird, G. A.↗

Bipropellant rocket exhaust plume analysis on the Galileo spacecraft

This paper describes efforts to quantify the contaminant flow field produced by 10 N thrust bipropellant rocket engines used on the Galileo spacecraft. The prediction of the composition of the rocket exhaust by conventional techniques is found to be inadequate to explain experimental observations of contaminant deposition on moderately cold (200 K) surfaces. It is hypothesized that low volatility contaminants are formed by chemical reactions which occur on the surfaces. The flow field calculations performed using the direct simulation Monte Carlo method give the expected result that the use of line-of-sight plume shields may have very little effect on the flux of vapor phase contaminant species to a surface, especially if the plume shields are located so close to the engine that the interaction of the plume with the shield is in the transition flow regime. It is shown that significant variations in the exhaust plume composition caused by nonequilibrium effects in the flow field lead to very low concentrations of species which have high molecular weights in the more rarefied regions of the flow field. Recommendations for the design of spacecraft plume shields and further work are made.

Guernsey, C. S.↗

Nonequilibrium thermal radiation for an aeroassist flight experiment vehicle

The direct-simulation Monte Carlo method incorporating a dissociating and ionizing gas model for air with thermal radiation is used to characterize the hypersonic flow about an axisymmetric representation of an aeroassist flight experiment (AFE) vehicle, whose freestream conditions correspond to selected points along the entry, aerobraking, and exit phases of the trajectory. Calculations for two trajectory conditions indicate that the radiative heating of the AFE forebody is lower than the convective heating, but becomes significant as the maximum convective heating rate condition is approached.

Moss, James N.↗

Interference effects on the hypersonic, rarefied flow about a flat plate

The Direct Simulation Monte Carlo method is used to study the hypersonic, rarified flow interference effects on a flat plate caused by nearby surfaces. Calculations focus on shock-boundary-layer and shock-lip interactions in hypersonic inlets. Results are presented for geometries consisting of a flat plate with different leading-edge shapes over a flat lower wall and a blunt-edge flat plate over a 5-degree wedge. The problems simulated correspond to a typical entry flight condition of 7.5 km/s at altitudes of 75 to 90 km. The results show increases in predicted local heating rates for shock-boundary-layer and shock-lip interactions that are quantitatively similar to those observed experimentally at much higher densities.

Wilmoth, Richard G.↗

Comparison of shock structure solutions using independent continuum and kinetic theory approaches

A vehicle traversing the atmosphere will experience flight regimes at high altitudes in which the thickness of a hypersonic shock wave is not small compared to the shock standoff distance from the hard body. When this occurs, it is essential to compute accurate flow field solutions within the shock structure. In this paper, one-dimensional shock structure is investigated for various monatomic gases from Mach 1.4 to Mach 35. Kinetic theory solutions are computed using the Direct Simulation Monte Carlo method. Steady-state solutions of the Navier-Stokes equations and of a slightly truncated form of the Burnett equations are determined by relaxation to a steady state of the time-dependent continuum equations. Monte Carlo results are in excellent agreement with published experimental data and are used as bases of comparison for continuum solutions. For a Maxwellian gas, the truncated Burnett equations are shown to produce far more accurate solutions of shock structure than the Navier-Stokes equations.

Fiscko, Kurt A.↗

Influence of afterbodies on AOTV flows

The role of vibration-dissociation coupling in flows representative of AOTV's with afterbodies is discussed. A physical model designed to reflect the behavior of anharmonic diatomic molecules at high temperatures was developed and incorporated into the direct simulation Monte Carlo method of Bird (1976). The model leads to reduced dissociation and increased convective heating for flows with finite catalytic surfaces.

Olynick, David P.↗

Direct simulation of rotational and vibrational nonequilibrium

The ways in which energy transfer is calculated in the Direct Simulation Monte Carlo method are presented. An energy exchange model that deals with translational and rotational modes is described. A model for simulating the transfer of energy between the translational and vibrational modes is presented as well.

Boyd, Iain D.↗

'Hypersonic slip flows' and issues on extending continuum model beyond the Navier-Stokes level

Results of a direct simulation Monte Carlo method for a hypersonic flow about a flat plate at a 40 deg angle of attack have been compared with corresponding results from a theory for fully viscous shock layers (FVSLs). Using the 13-moment equations for a Maxwell gas, it is demonstrated that nonequilibrium and equilibrium FVSL flows can be correlated. With the exception of the pressure density, most of the flow properties along a streamline (including heat flux and shear and normal stresses) are correctly predicted to leading order by the Navier-Stokes model.

Cheng, H. K.↗

Hypersonic rarefied flow about plates at incidence

The direct-simulation Monte Carlo method has been used in a numerical study of the transitional flow about two plate configurations at incidence; one of the two plates, both of which are 12 m long, has zero thickness, while the other has a thickness of 0.5 m and a node radius of 0.5 m. The flow conditions simulated are those of the Space Shuttle Orbiter during 7.5 km/hr reentry, in the 200-100 km altitude range encompassing most of the transitional flow for this vehicle. The results obtained clearly demonstrate that transitional effects are significant even at those altitudes where the flow about a typical space vehicle has been considered free-molecular.

Dogra, Virendra K.↗

Assessment of chemical nonequilibrium in rarefied hypersonic flow

Calculations are presented of chemical nonequilibrium in the rarefied hypersonic flowfields surrounding reentering spacecraft in the upper atmosphere. The direct simulation Monte Carlo method is the solution technique employed. Methods by which chemical reactions are included in this simulation scheme are focused upon. In particular, the coupling of vibrational and dissociative nonequilibrium is investigated. The assumption of a fully excited mode has very little effect on the calculated results when the vibrational mode does not contribute to the energy available for reaction. A steric factor for dissociative reactions is developed which includes the vibrational energy of the dissociating molecule. Calculations are presented for the flow along the stagnation streamline of a blunt body for conditions representative of the entry and exit phases of the Aeroassisted Flight Experiment vehicle at an altitude of 90 km. It is found that both the flow properties and surface quantities are significantly affected by the introduction of the steric factor.

Boyd, Iain D.↗

Efficient Monte Carlo simulation of rarefied flow in a small nozzle

The flow of carbon dioxide through a small conical nozzle is computed using two different numerical approaches. The first of these is the direct simulation Monte Carlo method in which the physical phenomena are modeled at the molecular level. The second technique solves the Navier-Stokes equations of fluid flow. The purpose of the paper is to compare the solutions generated by these two methods for the conditions prevailing in a small, low density nozzle. The comparison between the two sets of solutions reveals only small differences. However, the flow properties in the thick viscous layer close to the nozzle at the exit plane show important deviations which have implications for the subsequent determination of the plume backflow.

Boyd, Iain D.↗

Hypersonic flow - Physical issues and computational simulation

Computational results are presented on hypersonic flows, with emphasis on thermal and chemical nonequilibrium effects and the role played by grids. The continuum equations and their boundary conditions are presented, and the direct simulation Monte Carlo method is described. The importance of using the correct computational and physical models in describing hypersonic nonequilibrium flows is demonstrated.

Hassan, H. A.↗

Analysis of rotational nonequilibrium in standing shock waves of nitrogen

The one-dimensional standing shock wave is the simplest flow in which nonequilibrium effects may be considered. The direct-simulation Monte Carlo method (DSMC) has yielded excellent agreement with results reported for the case of nitrogen flow at Mach 1.7. The DSMC technique is presently used in conjunction with a variable energy transfer probability model in which energy is transferred between the translational and rotational modes via the Borgnakke-Larsen (1975) phenomenological model.

Boyd, Iain D.↗

Three-dimensional rarefied flow simulations for the Aeroassist Flight Experiment vehicle

Hypersonic rarefied flow about the Aeroassist Flight Experiment vehicle has been investigated using a three-dimensional direct simulation Monte Carlo method. Calculations are performed for the transitional flows encountered during the vehicle's atmospheric entry for altitudes of 110 and 100 km with an entry velocity of 9.9 km/s. The simulations are performed using a five-species reacting gas model that account for rotational and vibrational internal energies. The solutions indicate that dissociation is important at altitudes of 110 km and below. Results are presented for surface pressures, convective heating, flowfield structure, and aerodynamic coefficient variations with altitude.

Celenligil, M. Cevdet↗

New approach in direct-simulation of gas mixtures

Results are reported for an investigation of a new direct-simulation Monte Carlo method by which energy transfer and chemical reactions are calculated. The new method, which reduces to the variable cross-section hard sphere model as a special case, allows different viscosity-temperature exponents for each species in a gas mixture when combined with a modified Larsen-Borgnakke phenomenological model. This removes the most serious limitation of the usefulness of the model for engineering simulations. The necessary kinetic theory for the application of the new method to mixtures of monatomic or polyatomic gases is presented, including gas mixtures involving chemical reactions. Calculations are made for the relaxation of a diatomic gas mixture, a plane shock wave in a gas mixture, and a chemically reacting gas flow along the stagnation streamline in front of a hypersonic vehicle. Calculated results show that the introduction of different molecular interactions for each species in a gas mixture produces significant differences in comparison with a common molecular interaction for all species in the mixture. This effect should not be neglected for accurate DSMC simulations in an engineering context.

Chung, Chan-Hong↗

Numerical and experimental investigations of rarefied nozzle and plume flows of nitrogen

Numerical and experimental investigations are performed for the rarefied flow of nitrogen through a small nozzle which is expanded into near-vacuum conditions. Two different numerical studies are undertaken: the first employs a continuum approach in solving the Navier-Stokes equations, and the second employs a particle approach through use of the direct simulation Monte Carlo method (DSMC). The experimental investigation concerns the measurement of pressure, using a Pitot tube, in the nozzle exit plane and near-field of the plume. Comparison of the experimental and numerical data at the nozzle exit reveals that the DSMC technique provides the more accurate description of the expanding flow. It is discovered that the DSMC solutions are quite sensitive to the model employed to simulate the interaction between the gas and the nozzle wall surface. It is concluded that the simplistic fully diffuse model is quite satisfactory for the present application.

Boyd, Iain D.↗

Analysis of nozzle lip and backflow expansion of a small hydrazine thruster

Results are presented for computations made with the direct simulation Monte Carlo method for the expansion of the thick boundary layer of a small nozzle. The objective in the investigation is to assess the boundary-layer characteristics that most affect the flow properties of the expanding gas. In the computations, large degrees of species separation and thermal nonequilibrium are observed as the gas expands around the nozzle lip. These aspects of the flow indicate the necessity of treating the problem through a discrete particle approach. Assessment is made of the sensitivity of the calculated results to the form of the boundary layer initially assumed. It is found that the amount of backflow is related to the flow temperature close to the wall. Therefore, the opportunity presents itself for the reduction of the impingement potential of such thrusters through appropriate design considerations. Experimental procedures for verification of the type of computations undertaken in the current study are discussed. It is proposed that the calculations made in the backflow region offer the best opportunity.

Boyd, Iain D.↗