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 145 records · Page 8

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

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

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

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

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

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

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

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↗

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

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

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

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

Direct simulation of aerothermal loads for an aeroassist flight experiment vehicle

Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for the hypersonic flow about an elliptically blunted cone. The flow conditions are those for a proposed Aeroassist Flight Experiment (AFE) vehicle. The altitude range considered is that from 130 to 90 km which encompasses most of the transitional flow regime for the AFE vehicle, that is, the region bounded by free molecular and continuum flow. Freestream velocities of 9.9 to 7.5 km/sec are considered. The numerical simulations show that noncontinuum effects are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of about 130 km. Results presented highlight the thermal and chemical nonequilibrium nature of the flowfield and the impact of these effects on the surface heating and body drag. A calculation which included the additional effects of ionization and thermal radiation demonstrates that the inclusion of such efects would not significantly alter the surface quantities calculated in the present study. The radiative heating is negligible when compared with the convective heating, and the same would be true for the other conditions considered.

Dogra, Virendra K.↗

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

A continuum analysis of chemical nonequilibrium under hypersonic low-density flight conditions

Results of employing the continuum model of Navier-Stokes equations under the low-density flight conditions are presented. These results are obtained with chemical nonequilibrium and multicomponent surface slip boundary 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 new surface-slip boundary conditions in NS calculations, the surface heat transfer and other flowfield quantities adjacent to the surface are predicted favorably with the DSMC calculations from 75 km to 115 km in altitude. This suggests a much wider practical range for the applicability of Navier-Stokes solutions 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.

Gupta, R. N.↗

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

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

Chapman, Dean R.↗

Direct simulation of three-dimensional flow about the AFE vehicle at high altitudes

Three-dimensional hypersonic rarefied flow about the Aeroassist Flight Experiment (AFE) vehicle was studied using the direct simulation Monte Carlo (DSMC) technique. Results are presented for the transitional flow regime encountered between 120 and 200 km altitudes with a reentry velocity of 9.92 km/s. In the simulations, a five-species reacting real-gas model that accounts for internal energies (rotational and vibrational) is used. The results indicate that the transitional effects are significant even at an altitude of 200 km and influence the overall vehicle aerodynamics. For the cases considered, the aerodynamic coefficients, surface pressures, convective heating, and flow field structure variations with rarefaction effects are presented.

Celenligil, M. Cevdet↗