A Three-Level Cartesian Geometry Based Implementation of the DSMC Method

The data structures and overall algorithms of a newly developed 3-D direct simulation Monte Carlo (DSMC) program are outlined. The code employs an embedded 3-level Cartesian mesh, accompanied by a cut-cell algorithm to incorporate triangulated surface geometry into the adaptively refined Cartesian mesh. Such an approach enables decoupling of the surface mesh from the flow fie ld mesh, which is desirable for near-continuum flows, flows with large density variation, and also for adapti ve mesh refinement (AMR). Two separate data structures are proposed in order to separate geometry data from cell and particle information. The geometry data structure requires little memory so that each partition in a parallel simulation can store the entire mesh, potentially leading to better scalability and efficient AMR for parallel simulations. A simple and efficient AMR algorithm that maintains local cell size and time step consistent with the local mean-free-path and local mean collision time is detailed. The 3-level embedded Cartesian mesh combined with AMR allows increased flexibility for precise control of local mesh size and time-step, both vital for accurate and efficient DSMC simulation. Simulations highlighting the benefits of AMR and variable lo cal time steps will be presented along with DSMC results for 3-D flows with large density variations.

[1]  Ian D. Boyd,et al.  Rotational-Translational Energy Transfer in Rare? ed Nonequilibrium Flows , 1990 .

[2]  Iain D. Boyd,et al.  Strategies for Efficient Particle Resolution in the Direct Simulation Monte Carlo Method , 2000 .

[3]  Stefan Dietrich,et al.  Scalar and Parallel Optimized Implementation of the Direct Simulation Monte Carlo Method , 1996 .

[4]  G. J. LeBeau,et al.  A parallel implementation of the direct simulation Monte Carlo method , 1999 .

[5]  Thomas E. Schwartzentruber,et al.  Parallel Implementation of the Direct Simulation Monte Carlo Method For Shared Memory Architectures , 2010 .

[6]  Michail A. Gallis,et al.  DSMC convergence behavior of the hard-sphere-gas thermal conductivity for Fourier heat flow. , 2005 .

[7]  H. Lugt,et al.  Laminar flow behavior under slip−boundary conditions , 1975 .

[8]  Forrest E. Lumpkin,et al.  Virtual Sub-Cells for the Direct Simulation Monte Carlo Method , 2003 .

[9]  G. Bird Molecular Gas Dynamics and the Direct Simulation of Gas Flows , 1994 .

[10]  Mikhail S. Ivanov,et al.  Statistical simulation of reactive rarefied flows - Numerical approach and applications , 1998 .

[11]  Claus Borgnakke,et al.  Statistical collision model for Monte Carlo simulation of polyatomic gas mixture , 1975 .

[12]  Timothy J. Bartel,et al.  An Investigation of Two-Dimensional CAD Generated Models with Body Decoupled Cartesian Grids for DSMC , 2000 .

[13]  Michail A. Gallis,et al.  Direct simulation Monte Carlo convergence behavior of the hard-sphere-gas thermal conductivity for Fourier heat flow , 2006 .

[14]  M. Berger,et al.  Robust and efficient Cartesian mesh generation for component-based geometry , 1998 .

[15]  Michail A. Gallis,et al.  DSMC Convergence Behavior for Transient Flows. , 2007 .