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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Analysis of a semiclassical model for rotational transition probabilities

A semiclassical model proposed by Pearson and Hansen (1974) for computing collision-induced transition probabilities in diatomic molecules is tested by the direct-simulation Monte Carlo method. Specifically, this model is described by point centers of repulsion for collision dynamics, and the resulting classical trajectories are used in conjunction with the Schroedinger equation for a rigid-rotator harmonic oscillator to compute the rotational energy transition probabilities necessary to evaluate the rotation-translation exchange phenomena. It is assumed that a single, average energy spacing exists between the initial state and possible final states for a given collision.

Deiwert, G. S.↗
NASA NTRS1978

One-dimensional outgassing problem

A predominantly numerical study, using the direct-simulation Monte Carlo method, of the flow from outgassing surfaces with plane, cylindrical, or spherical geometry is presented. The numerical results for the plane flow case are compared with existing theoretical and experimental results and also with a new formulation of the jump conditions across the Knudsen layer between the surface and the continuum outflow. This outflow is sonic for the free expansion case but subsonic in the presence of back-pressure effects. The self-scattered return flux is 19% of the outgassed flux in the sonic outflow case and increases at lower Mach numbers. Thermal accomodation and reflection of this return flux produces an outgassing cooling effect. The cylindrical and spherical geometries introduce an additional Knudsen number defined by the ratio of the mean free path in the effusing gas to the radius. The plane flow results apply when this Knudsen number is less than 0.001 and, when it is of order unity or above, the return flux is inversely proportional to it.

Tran Cong, T.↗
NASA NTRS1978

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases - Translational rotational energy relaxation

Theoretical results pertaining to internally excited translational-rotational energy relaxation in a spatially uniform diatomic gas far removed from solid boundaries are obtained by solving the Boltzmann equation by means of the Monte Carlo direct simulation method. The analysis is based on calculations involving three different types of initial conditions: equilibrium, nonequilibrium-equipartition (i.e., equipartition is satisfied, but the distributions are perturbed), and nonequilibrium-nonequipartition (i.e., both equipartition and the distributions are perturbed). Results of monatomic-gas simulations are also included to facilitate comparisons with the coupled translational-rotational relaxation simulations, and some simulations for a normal shock-wave structure are briefly examined. The results show that: (1) single-step transitions are the significant mechanisms of intermodal energy transfer; (2) translational-rotational transitions are coupled most efficiently for low-lying states of rotationally excited molecules and least efficiently for highly rotationally excited molecules; and (3) relaxation occurs via a successive set of distributions that are not Maxwell-Boltzmann (nonlocal Maxwellian).

Yoshikawa, K. K.↗
NASA NTRS1979

Solution of Boltzmann equation for highly nonequilibrium diatomic gases rotational translational energy relaxation

The direct simulation Monte Carlo method is applied to solve the Boltzmann equation for collisions between internally excited diatomic gases in highly nonequilibrium states. The semiclassical transition probability is incorporated in the simulation for energy exchange between rotational and translational energy. The results provide details on the fundamental mechanisms of gas kinetics where analytical methods are impractical. The validity of the local Maxwellian assumption and relaxation time, rotational-translational energy transition, and a velocity analysis of the inelastic collision are discussed in detail.

Yoshikawa, K. K.↗
NASA NTRS1981

Spacecraft outgas ambient flow interaction

An investigation has been made of the extent to which outgassed or emitted molecules return to a spacecraft as a result of intermolecular collisions with ambient freestream molecules. The governing parameters are discussed and heuristic arguments are used to predict the combination of parameters that is most likely to describe the return flux ratio. This is then tested by an extensive set of numerical calculations using the direct simulation Monte Carlo method. Computations have been made for the flow past spheres and past circular cylinders with their axes normal in the stream. The numerical results lead to empirical expressions that can be used for engineering estimates in spacecraft contamination studies.

Bird, G. A.↗
NASA NTRS1982

Effects of translational nonequilibrium on vacuum plume expansions

The plume flowfield of a helium vent nozzle exhausting into a vacuum is defined by two techniques: the method of characteristics and the direct-simulation Monte Carlo method. The method of characteristics is shown to severely underpredict gas static temperatures due to the assumption of translational equilibrium inherent in the continuum formulation. Mass flow per unit solid angle is shown to be accurately predicted for this flowfield by the method of characteristics until the flow angle approaches within 10 deg of the maximum Prandtl-Meyer expansion angle. Improved treatment of the noncontinuum flow within the nozzle (near the lip) is postulated to have minor effects on the accuracy of this continuum method in comparison with translational nonequilibrium in the external flowfield. Possible treatment of translational nonequilibrium by the method of characteristics is discussed.

Guernsey, C. S.↗
NASA NTRS1984

Direct simulation of transitional flow for hypersonic reentry conditions

This paper presents results of flowfield calculations for typical hypersonic reentry conditions encountered by the nose region of the Space Shuttle Orbiter. Most of the transitional flow regime is covered by the altitude range of 150 to 92 km. Calculations were made with the Direct Simulation Monte Carlo (DSMC) method that accounts for translational, rotational, vibrational, and chemical nonequilibrium effects. Comparison of the DSMC heating results with both Shuttle flight data and continuum predictions showed good agreement at the lowest altitude considered. However, as the altitude increased, the continuum predictions, which did not include slip effects, departed rapidly from the DSMC results by overpredicting both heating and drag. The results demonstrate the effects of rarefaction on the shock and the shock layer, along with the extent of the slip and temperature jump at the surface. Also, the sensitivity of the flow structure to the gas-surface interaction model, thermal accommodation, and surface catalysis are studied.

Moss, J. N.↗
NASA NTRS1986

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

Direct simulation of hypersonic flows over blunt slender bodies

Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.

Moss, J. N.↗
NASA NTRS1986

Hypersonic low-density solutions of the Navier-Stokes equations with chemical nonequilibrium and multicomponent surface slip

Solutions of the Navier-Stokes equations with chemical nonequilibrium and multicomponent surface slip are presented along the stagnation streamline under low-density hypersonic flight conditions. The conditions analyzed are those encountered by the nose region of the Space Shuttle Orbiter during reentry. A detailed comparison of the Navier-Stokes (NS) results is made with the viscous shock-layer (VSL) and Direct Simulation Monte Carlo (DSMC) predictions. With the inclusion of surface-slip boundary conditions in NS calculations, the surface heat transfer and other flow field quantities adjacent to the surface are predicted favorably with the DSMC calculations from 75 km to 115 km in altitude. Therefore, the practical range for the applicability of Navier-Stokes solutions is much wider than previously thought. This is appealing because the continuum (NS and VSL) methods are commonly used to solve the fluid flow problems and are less demanding in terms of computer resource requirements than the noncontinuum (DSMC) methods. The NS solutions agree well with the VSL results for altitudes less than 92 km. An assessment is made of the frozen flow approximation employed in the VSL calculations.

Gupta, R. N.↗
NASA NTRS1986

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

Direct simulation of hypersonic transitional flow

A review of recent calculations obtained with Bird's direct simulation Monte Carlo (DSMC) method are presented for the transitional flowfield encountered at reentry conditions. Consequently, the emphasis is to simulate the real-gas effects resulting from a highly energetic nonequilibrium flow. The DSMC calculations for both wide-angle and slender bodies are compared with experimental and continuum calculations. The wide angle body calculations simulate the windward nose of the Shuttle Orbiter and projected aeroassisted orbital transfer vehicles. For the slender, blunted configurations, both two-dimensional and axisymmetric results are presented. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant.

Moss, J. N.↗
NASA NTRS1984

Investigation of the external flow analysis for density measurements at high altitude

The results of analysis performed on the external flow around the shuttle orbiter nose regions at the Shuttle Upper Atmosphere Mass Spectrometer (SUMS) inlet orifice are presented. The purpose of the analysis is to quantitatively characterize the flow conditions to facilitate SUMS flight data reduction and subsequent determination of orbiter aerodynamic force coefficients in the hypersonic rarefied flow regime. Experimental determination of aerodynamic force coefficients requires accurate simultaneous measurement of forces (or acceleration) and dynamic pressure along with independent knowledge of density and velocity. The SUMS provides independent measurement of dynamic pressure; however, it does so indirectly and requires knowledge of the relationship between measured orifice conditions and the dynamic pressure which can only be determined on the basis of molecule or theory for a winged configuration. Monte Carlo direct simulation computer codes were developed for both the flow field solution at the orifice and for the internal orifice flow. These codes were used to study issues associated with geometric modeling of the orbiter nose geometry and the modeling of intermolecular collisions including rotational energy exchange and a preliminary analysis of vibrational excitation and dissociation effects. Data obtained from preliminary simulation runs are presented.

Bienkowski, G. K.↗
NASA NTRS1987

Probabilistic SSME blades structural response under random pulse loading

The purpose is to develop models of random impacts on a Space Shuttle Main Engine (SSME) turbopump blade and to predict the probabilistic structural response of the blade to these impacts. The random loading is caused by the impact of debris. The probabilistic structural response is characterized by distribution functions for stress and displacements as functions of the loading parameters which determine the random pulse model. These parameters include pulse arrival, amplitude, and location. The analysis can be extended to predict level crossing rates. This requires knowledge of the joint distribution of the response and its derivative. The model of random impacts chosen allows the pulse arrivals, pulse amplitudes, and pulse locations to be random. Specifically, the pulse arrivals are assumed to be governed by a Poisson process, which is characterized by a mean arrival rate. The pulse intensity is modelled as a normally distributed random variable with a zero mean chosen independently at each arrival. The standard deviation of the distribution is a measure of pulse intensity. Several different models were used for the pulse locations. For example, three points near the blade tip were chosen at which pulses were allowed to arrive with equal probability. Again, the locations were chosen independently at each arrival. The structural response was analyzed both by direct Monte Carlo simulation and by a semi-analytical method.

Shiao, Michael↗
NASA NTRS1987

Nonequilibrium effects for hypersonic transitional flows

Presented are the results of numerical simulations of hypersonic flow about blunt cones and hemispherical nose configurations for reentry velocities of 7.5 and 10 km/s. Cone half angles 0, 5, and 10 deg are considered at zero angle of incidence; however, the focus is for the 5 deg cone. The body size and altitude ranges considered (70 to 110 km) are such that the flow is in the transitional regime. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method of Bird. The DSMC results are compared with those obtained with viscous shock-layer and Navier-Stokes methods. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant. The current calculations show that the binary scaling similitude provides a means of correlating the blunt body surface quantities in the hypersonic, transitional regime. Furthermore, for the higher velocity entry conditions, the results highlight some of the concerns in the application of multitemperature continuum formulations, particularly the use of some proposed functional relations for the chemical rate constants under thermodynamic nonequilibrium conditions.

Moss, James N.↗
NASA NTRS1987

Direct simulation of stagnation streamline flow for hypersonic reentry

The direct simulation Monte Carlo (DSMC) method is applied to simulate one-dimensional flow along the stagnation streamline. The freestream conditions considered are those encountered by the nose region of the Space Shuttle Orbiter during the hypersonic reentry. The range of altitudes (75 to 92 km) considered in the present calculations covers continuum to the less rarefied portion of the transition flow regime. The calculations account for thermal as well as the chemical nonequilibrium effects. The attention is focused on the flow structure along the stagnation streamline for different specified shock locations at a given altitude. The effect of shock location on the stagnation-point gas composition and heat transfer is analyzed for specific altitude; consequently, these results are appropriate for hemispherical bodies with different nose radii. Finally, comparison is made between the present DSMC results and the viscous shock-layer calculations, which show good agreement at the lowest altitude.

Dogra, Virendra K.↗
NASA NTRS1986

Numerical simulations of rarefied reentry flows

A review of recent calculations obtained with Bird's direct simulation Monte Carlo (DSMC) method is presented for the transitional flowfield encountered at reentry conditions. Consequently, the emphasis is the real-gas effects resulting from a highly energetic nonequilibrium flow. The DSMC calculations for both wide-angle and slender bodies are compared with continuum calculations. The wide-angle body calculations simulate the windward nose of the Shuttle Orbiter and projected aeroassisted orbital transfer vehicles. For the slender, blunted configurations, both two-dimensional and axisymmetric results are presented. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant.

Moss, James N.↗
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.↗