Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct simulation Monte Carlo (DSMC)”

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 NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗
NASA NTRS1995

Effects of Chemistry on Blunt-Body Wake Structure

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using direct simulation Monte Carlo (DSMC) and Navier-Stokes calculations. Particular emphasis is given to the effects of chemistry on the near-wake structure and on the surface quantities and the comparison of the DSMC results with the Navier-Stokes calculations. The flow conditions simulated are those experienced by a space vehicle at an altitude of 85 km and a velocity of 7 km/s during Earth entry. A steady vortex forms in the near wake for these freestream conditions for both chemically reactive and nonreactive air gas models. The size (axial length) of the vortex for the reactive air calculations is 25% larger than that of the nonreactive air calculations. The forebody surface quantities are less sensitive to the chemistry than the base surface quantities. The presence of the afterbody has no effect on the forebody flow structure or the surface quantities. The comparisons of DSMC and Navier-Stokes calculations show good agreement for the wake structure and the forebody surface quantities.

Dogra, Virendra K.↗
NASA NTRS1997

Rarefied Flows of Planetary Entry Capsules

The impact of rarefaction on entry capsules and spacecraft aerothermodynamics is emphasized for various aeroassist missions. The capability of the direct simulation Monte Carlo (DSMC) method to simulate such flows is demonstrated through examples of validation studies and applications. Several space flight projects and ground-based experiments are reviewed for which rarefaction effects have significant effect on spacecraft performance or experimental measurements. This review clearly demonstrates the significant role that the DSMC method plays in characterizing such flows.

James N Moss↗
NASA NTRS1998

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗
NASA NTRS1999

Study of the Transition Flow Regime using Monte Carlo Methods

This NASA Cooperative Agreement presents a study of the Transition Flow Regime Using Monte Carlo Methods. The topics included in this final report are: 1) New Direct Simulation Monte Carlo (DSMC) procedures; 2) The DS3W and DS2A Programs; 3) Papers presented; 4) Miscellaneous Applications and Program Modifications; 5) Solution of Transitional Wake Flows at Mach 10; and 6) Turbulence Modeling of Shock-Dominated Fows with a k-Enstrophy Formulation.

Hassan, H. A.↗
NASA NTRS1994

Aerothermodynamics of Pyrolizing Surfaces in Hypersonic Rarefied Flows

Direct simulation Monte Carlo (DSMC) calculations of rarefied flows about entry bodies typically employ a fixed surface temperature or a radiative-equilibrium energy balance to compute that temperature. Such boundary conditions neglect any effects of heat capacitance and heat conduction in the spacecraft heat shield and, therefore, provide an upper bound for the surface temperature. Such calculations also neglect pyrolysis from the heat shield which can be significant for a high-energy incident flow at very low densities. Accurate prediction of both heating and aerodynamic forces requires including pyrolysis and surface heat transfer in the models for gas-surface interaction employed in DSMC methods. Although these physical models have long appeared in various continuum flow calculation codes, they have only recently appeared in DSMC codes which are required to simulate rarefied flows during entry at high altitudes. In the current implementation, routines from the widely distributed Charring Material Thermal Response and Ablation (CMA) program are coupled into a DSMC code to calculate the one-dimensional heat transfer into the carbon phenolic heat shield at each point on a vehicle surface. Temperature-dependent material properties, surface re-radiation, and in-depth pyrolysis were included in the calculation, but surface ablation was neglected. Sample calculations for entry of the Galileo probe into the atmosphere of Jupiter demonstrate that including pyrolysis in the model leads to significant differences in predicted aerodynamics. Granted, the drag coefficient does not depend strongly on the surface temperature which can itself be significantly below the radiative equilibrium value during entry. However, the surface mass flux due to pyrolysis of the material is significant once the probe drops to altitudes characterized by transition flow. This leads to a noticeable increase in drag and a decrease in heating compared to a body without pyrolysis.

Haas, Brian L.↗
NASA NTRS2001

Hypersonic Flows About a 25 degree Sharp Cone

This paper presents the results of a numerical study that examines the surface heating discrepancies observed between computed and measured values along a sharp cone. With Mach numbers of an order of 10 and the freestream length Reynolds number of an order of 10 000, the present computations have been made with the direct simulation Monte Carlo (DSMC) method by using the G2 code of Bird. The flow conditions are those specified for two experiments conducted in the Veridian 48-inch Hypersonic Shock Tunnel. Axisymmetric simulations are made since the test model was assumed to be at zero incidence. Details of the current calculations are presented, along with comparisons between the experimental data, for surface heating and pressure distributions. Results of the comparisons show major differences in measured and calculated results for heating distributions, with differences in excess of 25 percent for the two cases examined.

Moss, James N.↗
NASA NTRS2002

Hypersonic Shock Interactions About a 25 deg/65 deg Sharp Double Cone

This paper presents the results of a numerical study of shock interactions resulting from Mach 10 air flow about a sharp double cone. Computations are made with the direct simulation Monte Carlo (DSMC) method by using two different codes: the G2 code of Bird and the DAC (DSMC Analysis Code) code of LeBeau. The flow conditions are the pretest nominal free-stream conditions specified for the ONERA R5Ch low-density wind tunnel. The focus is on the sensitivity of the interactions to grid resolution while providing information concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗
NASA NTRS2004

Radiolytic Model for Chemical Composition of Europa's Atmosphere and Surface

The overall objective of the present effort is to produce models for major and selected minor components of Europa s neutral atmosphere in 1-D versus altitude and in 2-D versus altitude and longitude or latitude. A 3-D model versus all three coordinates (alt, long, lat) will be studied but development on this is at present limited by computing facilities available to the investigation team. In this first year we have focused on 1-D modeling with Co-I Valery Shematovich s Direct Simulation Monte Carlo (DSMC) code for water group species (H2O, O2, O, OH) and on 2-D with Co-I Mau Wong's version of a similar code for O2, O, CO, CO2, and Na. Surface source rates of H2O and O2 from sputtering and radiolysis are used in the 1-D model, while observations for CO2 at the Europa surface and Na detected in a neutral cloud ejected from Europa are used, along with the O2 sputtering rate, to constrain source rates in the 2-D version. With these separate approaches we are investigating a range of processes important to eventual implementation of a comprehensive 3-D atmospheric model which could be used to understand present observations and develop science requirements for future observations, e.g. from Earth and in Europa orbit. Within the second year we expect to merge the full water group calculations into the 2-D version of the DSMC code which can then be extended to 3-D, pending availability of computing resources. Another important goal in the second year would be the inclusion of sulk and its more volatile oxides (SO, SO2).

Cooper, John F.↗
NASA NTRS2004

Simulations of Hyperthermal Oxygen Beam Exposures

Pulsed sources of hyper-thermal O-atoms are now being extensively used to simulate low-earth orbit (LEO) surface exposure environments. The peak flux of these sources is many orders of magnitude larger than the corresponding LEO flux. Although it is desirable to accelerate the test by using higher fluxes than found in LEO, even commonly used fluxes are large enough to produce multi-collision effects by causing a build-up of gas at the sample surface. In this paper we characterize the physical consequences to the experiment using the direct simulation Monte Carlo (DSMC) method, DSMC allows us to extract the distributions of energy and impact angle for the O-atoms that reach the surface, and to record how strongly the gas build-up at the target assembly deflects flux from downstream instrumentation. By considering a range of source fluxes, we determine the onset conditions and severity of these multi-collision effects. We find that even at common experimental fluxes with a normally incident beam string a flat surface sample, the energy distribution of incident O-atoms broadens and develops a significant low-energy tail. The angular distributions also broaden significantly. The number of O-atoms that reach downstream instrumentation is decreased by approximately 50%. These simulations will aid in the calibration of ground-based O-atom measurements,a nd provide a model for the energy and angular distributions of O-atoms that actually impinge on surface samples.

Finchum, A.↗
NASA NTRS2004

Shock-Wave/Boundary-Layer Interactions in Hypersonic Low Density Flows

Results of numerical simulations of Mach 10 air flow over a hollow cylinder-flare and a double-cone are presented where viscous effects are significant. The flow phenomena include shock-shock and shock- boundary-layer interactions with accompanying flow separation, recirculation, and reattachment. The purpose of this study is to promote an understanding of the fundamental gas dynamics resulting from such complex interactions and to clarify the requirements for meaningful simulations of such flows when using the direct simulation Monte Carlo (DSMC) method. Particular emphasis is placed on the sensitivity of computed results to grid resolution. Comparisons of the DSMC results for the hollow cylinder-flare (30 deg.) configuration are made with the results of experimental measurements conducted in the ONERA RSCh wind tunnel for heating, pressure, and the extent of separation. Agreement between computations and measurements for various quantities is good except that for pressure. For the same flow conditions, the double- cone geometry (25 deg.- 65 deg.) produces much stronger interactions, and these interactions are investigated numerically using both DSMC and Navier-Stokes codes. For the double-cone computations, a two orders of magnitude variation in free-stream density (with Reynolds numbers from 247 to 24,7 19) is investigated using both computational methods. For this range of flow conditions, the computational results are in qualitative agreement for the extent of separation with the DSMC method always predicting a smaller separation region. Results from the Navier-Stokes calculations suggest that the flow for the highest density double-cone case may be unsteady; however, the DSMC solution does not show evidence of unsteadiness.

Moss, James N.↗
NASA NTRS1998

Rarefied Aerothermodynamic Predictions for Mars Global Surveyor

Mars Global Surveyor (MGS) successfully completed its first phase of aerobraking in early 1998 and is the rst planetary mission to use aerobraking as a primary means of customizing its orbit to achieve its mission objectives. The aerobraking requirements together with post-launch anomalies presented a unique challenge to provide accurate predictions of the aerothermodynamic environment of the spacecraft in the rare transitional flow regime. Direct Simulation Monte Carlo (DSMC) and free molecular techniques were used to provide heating and aerodynamic predictions and to investigate a variety of rarefied flow phenomena across the regime; MGS is the first major planetary mission in which rare flow predictions have played such a critical role all the way through design, mission planning, and operational phases. This paper summarizes these studies with emphasis on transitional-flow and gas-surface interaction phenomena.

Wilmoth, R. G.↗
NASA NTRS1997

Numerical Study of a Continuum Sonic Jet Interacting with a Rarefield Flow

The results of a numerical study with flow and boundary conditions based on an experiment of a continuum sonic jet interacting with rarefied flow about a sharp leading edge flat plate at zero incidence are presented. Comparisons are made between computational fluid dynamics (CFD) and direct simulation Monte Carlo (DSMC) solutions which provide an assessment of applying each technique to the flow conditions of the experiment. An analysis of the CFD results revealed a correlation between the interaction interface of the jet continuum breakdown surface and a non-dimensional parameter derived from jet and free stream flow conditions. Using the breakdown surface from the correlation, the continuum jet was uncoupled from the interaction, thus allowing an uncoupled CFD-DSMC solution to be obtained. Also, a nearest neighbor collision algorithm, similar to the subcell technique, was implemented in the DSMC solution technique. The comparison between CFD and DSMC results shows good qualitative agreement in the interaction region and good quantitative agreement elsewhere.

Glass, Christoper E.↗
NASA NTRS1989

Direct Simulation of Reentry Flows with Ionization

The Direct Simulation Monte Carlo (DSMC) method is applied in this paper to the study of rarefied, hypersonic, reentry flows. The assumptions and simplifications involved with the treatment of ionization, free electrons and the electric field are investigated. A new method is presented for the calculation of the electric field and handling of charged particles with DSMC. In addition, a two-step model for electron impact ionization is implemented. The flow field representing a 10 km/sec shock at an altitude of 65 km is calculated. The effects of the new modeling techniques on the calculation results are presented and discussed.

Carlson, Ann B.↗
NASA NTRS1991

Radiation Modeling with Direct Simulation Monte Carlo

Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available rate data. The new method is compared with the current modeling technique and both techniques are compared with available experimental data. The differences in the results are evaluated. The test case is based on experimental measurements from the AVCO-Everett Research Laboratory electric arc-driven shock tube of a normal shock wave in air at 10 km/s and .1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to di erent ow conditions.

Carlson, Ann B.↗
NASA NTRS2006

Blunt Body Aerodynamics for Hypersonic Low Density Flows

Numerical simulations are performed for the Apollo capsule from the hypersonic rarefied to the continuum regimes. The focus is on flow conditions similar to those experienced by the Apollo 6 Command Module during the high altitude portion of its reentry. The present focus is to highlight some of the current activities that serve as a precursor for computational tool assessments that will be used to support the development of aerodynamic data bases for future capsule flight environments, particularly those for the Crew Exploration Vehicle (CEV). Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction; that is, free molecular to continuum conditions. Also, aerodynamic data are presented that shows their sensitivity to a range of reentry velocities, encompassing conditions that include reentry from low Earth orbit, lunar return, and Mars return velocities (7.7 to 15 km/s). The rarefied results obtained with direct simulation Monte Carlo (DSMC) codes are anchored in the continuum regime with data from Navier-Stokes simulations.

Moss, James N.↗
NASA NTRS2006

Cassini-Huygens Aerodynamics with Comparison to Flight

An analysis has been performed on the aerodynamics of the Cassini spacecraft as it passed through the atmosphere of Titan with and without the Huygens probe attached. The free stream density of Titan s atmosphere was measured by two methods. However, these methods resulted in very different values of density, one result being 3-5 times higher than the other. In an attempt to understand the source of this discrepancy and verify the assumptions made by and the original engineering analysis performed by the Jet Propulsion Laboratory (JPL), free molecular and direct simulation Monte Carlo (DSMC) analyses were performed for two atmospheric passes. The drag coefficient was calculated using an area based on a Monte Carlo area determination program. Although the source of the discrepancy has not been determined, it has been confirmed that the original, simplified force-and-moment analysis performed by JPL produced results that were comparable to the high-fidelity DSMC analysis and that the source of the discrepancy lies elsewhere.

Liechty, Derek S.↗
NASA NTRS2006

Orion Aerodynamics for Hypersonic Free Molecular to Continuum Conditions

Numerical simulations are performed for the Orion Crew Module, previously known as the Crew Exploration Vehicle (CEV) Command Module, to characterize its aerodynamics during the high altitude portion of its reentry into the Earth's atmosphere, that is, from free molecular to continuum hypersonic conditions. The focus is on flow conditions similar to those that the Orion Crew Module would experience during a return from the International Space Station. The bulk of the calculations are performed with two direct simulation Monte Carlo (DSMC) codes, and these data are anchored with results from both free molecular and Navier-Stokes calculations. Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction, that is, for free molecular to continuum conditions (Knudsen numbers of 111 to 0.0003). Also included are aerodynamic data as a function of angle of attack for different levels of rarefaction and results that demonstrate the aerodynamic sensitivity of the Orion CM to a range of reentry velocities (7.6 to 15 km/s).

Moss, James N.↗