# Direct simulation Monte Carlo

Direct Simulation Monte Carlo (DSMC) method uses probabilistic (Monte Carlo) simulation to solve the Boltzmann equation for finite Knudsen number fluid flows.

The DSMC method was proposed by Prof. Graeme Bird,[1][2][3] Emeritus Professor of Aeronautics, University of Sydney. DSMC is a numerical method for modeling rarefied gas flows, in which the mean free path of a molecule is of the same order (or greater) than a representative physical length scale (i.e. the Knudsen number Kn is greater than 1). In supersonic and hypersonic flows rarefaction is characterized by Tsein's parameter, which is equivalent to the product of Knudsen number and Mach number (KnM) or M$^2$/Re, where Re is the Reynolds number.[4][5] In these rarefied flows, the Navier-Stokes equations can be inaccurate. The DSMC method has been extended to model continuum flows (Kn < 1) and the results can be compared with Navier stokes solutions.

The DSMC method models fluid flows using simulation molecules which represent a large number of real molecules in a probabilistic simulation to solve the Boltzmann equation. Molecules are moved through a simulation of physical space in a realistic manner that is directly coupled to physical time such that unsteady flow characteristics can be modeled. Intermolecular collisions and molecule-surface collisions are calculated using probabilistic, phenomenological models. Common collision models include the Hard Sphere model, the Variable Hard Sphere (VHS) model, and the Variable Soft Sphere (VSS) model. The fundamental assumption of the DSMC method is that the molecular movement and collision phases can be decoupled over time periods that are smaller than the mean collision time.

Currently the DSMC method has been applied to the solution of flows ranging from estimation of the Space Shuttle re-entry aerodynamics, to the modeling micro-electro-mechanical systems (MEMS).

## DSMC Software

Multiple implementations of the DSMC method exist:

• DS1V, DS2V and DS3V are the original DSMC programs written by Prof. Bird. These programs have a visual user interface that can be used for configuration and post processing.
• dsmcFoam is a DSMC solver for 2D and 3D flows. dsmcFoam is part of the open source CFD package OpenFOAM. For more details, see Scanlon, T., Roohi, E., White, C., Darbandi, M., Reese, J., An Open Source, Parallel DSMC Code for Rarefied Gas Flows in Arbitrary Geometries , Computers and Fluids , Volume 39, 2010, Pages 2078-2089, http://e.roohi.profcms.um.ac.ir/index.php?option=com_content&view=article&id=8116&Itemid=29226
• MONACO is a DSMC solver devised at Cornell University by Dr. Stefan Dietrich and Prof. Iain Boyd's Nonequilibrium Gas and Plasma Dynamics Laboratory at the University of Michigan.[6]
• PI-DSMC is a commercial DSMC software package for 2D and 3D flows.
• SMILE (Statistical Modeling in Low-density Environment) is a general purpose 2D/3D parallel DSMC software system developed since 1998 by Computational Aerodynamics Laboratory (L7) at the Khristianovich Institute of Theoretical and Applied Mechanics, Siberian Division of the Russian Academy of Sciences. SMILE has been the principal aerodynamic analysis tool for high-altitude stages of reentry of the Mir Space Station as well as many other Russian and European space vehicle projects.
• DAC is a general purpose DSMC code developed by NASA at the Johnson Space and Langley Research Centers. It employs a two level mesh using Cartesian volumes and employs the cut cell algorithm developed by Prof. Tom Schwartzentruber's group at the University of Minnesota. Both scalar and parallel versions exist with the parallel version using the Message Passing Interface (MPI) and domain decomposition. DAC was designed to handle difficult problems such as complex geometries (for example, the International Space Station) and the plume impingement which may occur during the rendezvous of two spacecraft. DAC is classified as ITAR and distribution is restricted to United States users. Requests for DAC should be directed to the Technology Transfer Office at the NASA Johnson Space Center.
• MGDS is a fully 3D DSMC solver incorporating three level adaptive mesh refinement and a cut cell algorithm developed by Prof. Tom Schwartzentruber's group at the University of Minnesota.
• Molflow is a 3D DSMC simulator currently developed at CERN for the simulation of free molecular flow (vacuum) systems.
• SAMADII/SCiV (Statistical Contact in Vacuum) is general purpose 3D DSMC software system based on multi-GPUs.
• HAP (Hypersonic Aerothermodynamics Particle code) is a DSMC code developed at the U.S. Air Force Research Laboratory for high speed flight and space applications.
• SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer), an Open Source 2 & 3D DSMC simulator optimized for parallel computing and developed at Sandia National Laboratories. Written in C++, SPARTA is designed to be easy to modify or extend with new functionality. Code is distributed under GPL, and available from the project website
• PICLas is a parallel, three-dimensional PIC-DSMC solver developed cooperatively by the Institute of Space Systems and Institute of Aerodynamics and Gas Dynamics at the University of Stuttgart.[7] It is a flexible simulation suite for the computation of reactive plasma flows, where the PIC and DSMC modules can be utilized separately. Application areas include the simulation of electric propulsion systems, atmospheric entry manoeuvres and laser ablation.

## References

1. ^ G. A. Bird, 'Approach to translational equilibrium in a rigid sphere gas', Phys. Fluids, 6, p1518 (1963).
2. ^ G. A. Bird, Molecular Gas Dynamics, Clarendon, Oxford (1976)
3. ^ G. A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Claredon, Oxford (1994)
4. ^ H. S. Tsein, 1948 ‘Superaerodynamics, mechanics of rarefied gases’ J. Aerospace Sci, 13, p342, 1946.
5. ^ M. N. Macrossan, 'Scaling Parameters for Hypersonic Flow: Correlation of Sphere Drag Data'. In: M. S. Ivanov and A. K. Rebrov, 25th International Symposium on Rarefied Gas Dynamics, Siberian Division of the Russian Academy of Sciences, p.759 (2007).
6. ^ Dietrich, S.; Boyd, I.: "Scalar and Parallel Optimized Implementation of the Direct Simulation Monte Carlo Method," Journal of Computational Physics, 126:328-42, 1996.
7. ^ Munz, C.-D., Auweter-Kurtz, M., Fasoulas, S. et al.: "Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows," Comptes Rendus Mécanique 342(10-11), 662–670, 2014. http://dx.doi.org/10.1016/j.crme.2014.07.005