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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

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

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗
NASA NTRS1992

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗
NASA NTRS1991

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗
NASA NTRS1992

Hypersonic rarefied flow past spheres including wake structure

Results of a numerical study using the direct simulation Monte Carlo method are presented for hypersonic rarefied flow past spheres. The flow conditions considered are those corresponding to low density wind tunnel test conditions. The set of the experimental conditions for the calculations encompasses the transitional to near-continuum flow regimes. Comparison of the calculated drag with experimental results shows good agreement to well within the experimental error. Particular attention is focused on the wake structure. Calculations show that the wake is very rarefied with considerable thermal nonequilibrium for all the cases considered. No flow separation is observed in the wake for the near-continuum case where a vortex has been predicted by Navier-Stokes type calculations.

Dogra, Virendra K.↗
NASA NTRS1992

Internal structure of shock waves in disparate mass mixtures

The detailed flow structure of a normal shock wave for a gas mixture is investigated using the direct-simulation Monte Carlo method. A variable diameter hard-sphere (VDHS) model is employed to investigate the effect of different viscosity temperature exponents (VTE) for each species in a gas mixture. Special attention is paid to the irregular behavior in the density profiles which was previously observed in a helium-xenon experiment. It is shown that the VTE can have substantial effects in the prediction of the structure of shock waves. The variable hard-sphere model of Bird shows good agreement, but with some limitations, with the experimental data if a common VTE is chosen properly for each case. The VDHS model shows better agreement with the experimental data without adjusting the VTE. The irregular behavior of the light-gas component in shock waves of disparate mass mixtures is observed not only in the density profile, but also in the parallel temperature profile. The strength of the shock wave, the type of molecular interactions, and the mole fraction of heavy species have substantial effects on the existence and structure of the irregularities.

Chung, Chan-Hong↗
NASA NTRS1992

Decoupled predictions of radiative heating in air using a particle simulation method

A particle simulation technique with decoupled radiation is used to estimate the radiative emission along the stagnation streamline and the radiative heating at the stagnation point of a blunt-nosed vehicle entering the earth's atmosphere at hypersonic speed. The direct simulation Monte Carlo (DSMC) method is used to compute the fluid mechanics of the weakly ionized flow. The radiative emission is computed using the NEQAIR computer code with the radiation decoupled from the flowfield solution. It is concluded that the new models make it possible to decrease the predicted total radiative heating at the stagnation point of the vehicle by a factor of 15. A comparison of the DSMC approach with a continuum flow model shows that the total predicted radiative heating estimates agree within a factor of 2.

Boyd, Iain D.↗
NASA NTRS1992

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗
NASA NTRS1992

DSMC analysis of species separation in rarefied nozzle flows

The direct-simulation Monte Carlo method has been used to investigate the behavior of a small amount of a harmful species in the plume and the backflow region of nuclear thermal propulsion rockets. Species separation due to pressure diffusion and nonequilibrium effects due to rapid expansion into a surrounding low-density environment are the most important factors in this type of flow. It is shown that a relatively large amount of the lighter species is scattered into the backflow region and the heavier species becomes negligible in this region due to the extreme separation between species. It is also shown that the type of molecular interaction between the species can have a substantial effect on separation of the species.

Chung, Chan-Hong↗
NASA NTRS1992

Direct simulation with vibration-dissociation coupling

The majority of implementations of the Direct Simulation Monte Carlo (DSMC) method of Bird do not account for vibration-dissociation coupling. Haas and Boyd have proposed the vibrationally-favored dissociation model to accomplish this task. This model requires measurements of induction distance to determine model constants. A more general expression has been derived that does not require any experimental input. The model is used to calculate one-dimensional shock waves in nitrogen and the flow past a lunar transfer vehicle (LTV). For the conditions considered in the simulation, the influence of vibration-dissociation coupling on heat transfer in the stagnation region of the LTV can be significant.

Hash, David B.↗
NASA NTRS1992

A new two-temperature dissociation model for reacting flows

A new two-temperature dissociation model for flows undergoing compression is derived from kinetic theory. The model minimizes uncertainties associated with the two-temperature model of Park. The effects of the model on AOTV type flowfields are examined and compared with the Park model. Calculations are carried out for flows with and without ionization. When considering flows with ionization, a four temperature model is employed. For Fire II conditions, the assumption of equilibrium between the vibrational and electron-electronic temperatures is somewhat poor. A similar statement holds for the translational and rotational temperatures. These trends are consistent with results obtained using the direct simulation Monte Carlo method.

Olynick, David R.↗
NASA NTRS1992

Evaluation of thermochemical models for particle and continuum simulations of hypersonic flow

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered that are dominated in turn by vibrational relaxation, weak dissociation, strong dissociation and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Boyd, Iain D.↗
NASA NTRS1992

Measurement and analysis of a small nozzle plume in vacuum

Measurements of Pitot pressure and flow angle were made in the plume of a nozzle flowing nitrogen and exhausting to a vacuum. The measurements were compared to results from a numerical simulation of the flow that was based on kinetic theory and used the direct-simulation Monte Carlo (DSMC) method. Numerical results were compared with measurements made in the plume at various axial and radial stations. Total pressure measurements were made with Pitot tubes sized for specific regions of the plume. Flow angle measurements were made with a conical probe. The measurement area for flow angle extended to 160 mm (5 exit diameters) downstream of the nozzle exit plane and radially to 60 mm (1.9 exit diameters) from the plume axis. The total pressure measurements extended 480 mm (16 exit diameters) downstream and radially to 60 mm. Comparisons of computed results from the DSMC method with measurements of flow angle displayed improved agreement with increasing distance from the exit plane. Pitot pressures computed from the DSMC method were in reasonably good agreement with experimental results over the entire measurement area.

Penko, Paul F.↗
NASA NTRS1992

Investigation of gas/particle heat transfer rates in solid rocket motors

The ability of the current Nusselt number prediction technique developed by Kavanau (1955) to accurately predict alumina particle heat transfer rates in solid rocket nozzles and plumes is investigated. For the solid rocket motors SRMS) analyzed, the transitional regime is the dominant regime for the majority of particles in the flowfield. The analytical approach to determine accuracy of the Kavanau correlation utilized the G2R Direct Simulation Monte Carlo code. With this method, both sphere drag, and heat transfer rates were predicted. The sphere drag prediction were compared to the Hermsen, and Henderson drag correlations, while the heat transfer results were compared to the current theory. Results have indicated, that the predicted drag coefficient is bounded by the drag correlations considered. However, the Nusselt number varies significantly from the extrapolated profile through all flowfield regimes.

Moylan, B.↗
NASA NTRS1991

A review of the kinetic detail required for accurate predictions of normal shock waves

Several aspects of the kinetic models used in the collision phase of Monte Carlo direct simulations have been studied. Accurate molecular velocity distribution function predictions require a significantly increased number of computational cells in one maximum slope shock thickness, compared to predictions of macroscopic properties. The shape of the highly repulsive portion of the interatomic potential for argon is not well modeled by conventional interatomic potentials; this portion of the potential controls high Mach number shock thickness predictions, indicating that the specification of the energetic repulsive portion of interatomic or intermolecular potentials must be chosen with care for correct modeling of nonequilibrium flows at high temperatures. It has been shown for inverse power potentials that the assumption of variable hard sphere scattering provides accurate predictions of the macroscopic properties in shock waves, by comparison with simulations in which differential scattering is employed in the collision phase. On the other hand, velocity distribution functions are not well predicted by the variable hard sphere scattering model for softer potentials at higher Mach numbers.

Muntz, E. P.↗
NASA NTRS1991

Enhancements and applications of DSMC for hypersonic rarefied flows

A review of the direct simulation Monte Carlo (DSMC) method of Bird is presented. The DSMC method provides the capability of simulating real gas flows in the rarefied flow regime. Recent developments and applications of the method for hypersonic flows are reported for both ground-based tests and during entry. The flows considered were three-dimensional, and both an existing general code and a new code under development were used.

Moss, James N.↗
NASA NTRS1991

Application of the DSMC method to hypersonic flow about a delta wing

The direct simulation Monte Carlo method was applied to a hypersonic rarefied flow about a delta wing at angle of incidence. Three-dimensional calculations were performed using a single-species gas model that accounts for energy exchange with the internal modes. Computed results were compared with the available experimental data obtained in a nitrogen wind-tunnel experiment. Discrepancies were shown to exist between the experimental and calculated results. These differences may be caused by the inadequately defined freestream conditions.

Celenligil, M. C.↗
NASA NTRS1992

Experimental and numerical investigations of low-density nozzle and plume flows of nitrogen

New experimental data are used to show that the direct simulation Monte Carlo (DSMC) method provides an accurate description of a low-density flow in the nozzle and the near-field expansion of a small rocket for two slightly different experimental configurations. These results verify the DSMC method in an expansion flow for the first time. From a number of different gas/surface interaction models, it is found that fully diffuse reflection gives the best agreement with experiment.

Boyd, Iain D.↗