A hybrid continuum/particle approach for modeling subsonic, rarefied gas flows

A hybrid approach that combines a continuum approach solving the Navier-Stokes equations and a particle method called the information preservation (IP) method is implemented to simulate subsonic, rarefied gas flows accurately and efficiently. The coupling between the continuum and particle approaches relies on a continuum/particle interface, which is investigated in detail. The IP method preserves information at the macroscopic level that allows the continuum approach to directly use the information from the particle region. In order to correctly generate particles from the continuum region, two strategies are proposed. One strategy adopts a Marshak-type condition, which requires many particles in a cell to control the flux fluctuation due to the microscopic motion of particles. In the second strategy, reservoir cells and buffer cells are used, and this is the approach adopted in our general hybrid scheme. The location of the interface is determined using a continuum breakdown parameter. Studies show that the continuum breakdown parameter suggested by Garcia et al. can be used to determine the location of the interface, but the implementation of the interface also affects the location of the interface. Simulation of a flow over a flat plate shows the performance of the hybrid approach and reveals the effects of the cutoff value to a continuum breakdown parameter. This approach is also applied to study the aerodynamics of a micro-scale airfoil, which is very difficult for a single method. The hybrid approach generally spends less computational time than the IP method for rarefied gas flows, and the numerical performance of the hybrid approach depends on the number ratio of continuum cells to particle cells.

[1]  E. Oran,et al.  A New Hybrid Algorithm: Navier-Stokes as a DSMC Filter , 1997 .

[2]  Alejandro L. Garcia,et al.  Adaptive Mesh and Algorithm Refinement Using Direct Simulation Monte Carlo , 1999 .

[3]  Hassan Hassan,et al.  Assessment of schemes for coupling Monte Carlo and Navier-Stokes solution methods , 1996 .

[4]  Jing Fan,et al.  Statistical simulation of low-speed rarefied gas flows , 2001 .

[5]  S. A. Schaaf,et al.  Flow of rarefied gases , 1961 .

[6]  Iain D. Boyd,et al.  Predicting continuum breakdown in hypersonic viscous flows , 2003 .

[7]  Axel Klar,et al.  An adaptive domain decomposition procedure for Boltzmann and Euler equations , 1998 .

[8]  Quanhua Sun,et al.  A Direct Simulation Method for Subsonic, Microscale Gas Flows , 2002 .

[9]  Jing Fan,et al.  Direct Simulation Methods for Low-Speed Microchannel Flows , 2000 .

[10]  Moulay D. Tidriri,et al.  Coupling Boltzmann and Navier-Stokes Equations by Friction , 1996 .

[11]  I. Boyd,et al.  Computation of Rarefied Gas Flows Around a NACA 0012 Airfoil , 1999 .

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

[13]  David B. Goldstein,et al.  Hybrid Euler/Direct Simulation Monte Carlo Calculation of Unsteady Slit Flow , 2000 .

[14]  S. Tiwari,et al.  Coupling of the Boltzmann and Euler Equations with Automatic Domain Decomposition , 1998 .

[15]  Alejandro L. Garcia,et al.  Generation of the Chapman-Enskog Distribution , 1998 .

[16]  O. Aktas,et al.  A Combined Continuum/DSMC Technique for Multiscale Analysis of Microfluidic Filters , 2002 .

[17]  M. Gad-el-Hak The Fluid Mechanics of Microdevices—The Freeman Scholar Lecture , 1999 .

[18]  Quanhua Sun,et al.  Information preservation methods for modeling micro-scale gas flows , 2003 .

[19]  Graeme A. Bird,et al.  Breakdown of translational and rotational equilibrium in gaseous expansions , 1970 .

[20]  Berni J. Alder Highly discretized dynamics , 1997 .

[21]  T. G. Cowling,et al.  The mathematical theory of non-uniform gases , 1939 .

[22]  Hassan Hassan,et al.  Two-dimensional coupling issues of hybrid DSMC/Navier-Stokes solvers , 1997 .

[23]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[24]  R. Nance,et al.  Role of Boundary Conditions in Monte Carlo Simulation of Microelectromechanical Systems , 1998 .

[25]  M. Gad-el-Hak,et al.  Micro Flows: Fundamentals and Simulation , 2002 .

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

[27]  Graham V. Candler,et al.  Predicting failure of the continuum fluid equations in transitional hypersonic flows , 1994 .

[28]  Quanhua Sun,et al.  Numerical Simulation of Gas Flow over Microscale Airfoils , 2001 .

[29]  I. Boyd,et al.  A new energy flux model in the DSMC-IP Method for Nonequilibrium Flows , 2003 .

[30]  Hassan Hassan,et al.  Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices , 1997 .

[31]  N. Hadjiconstantinou Regular Article: Hybrid Atomistic–Continuum Formulations and the Moving Contact-Line Problem , 1999 .

[32]  Daniel A. Erwin,et al.  Two-dimensional hybrid continuum/particle approach for rarefied flows , 1992 .

[33]  Shigeru Sunada,et al.  Airfoil Section Characteristics at a Low Reynolds Number , 1997 .

[34]  David B. Goldstein,et al.  Hybrid Euler/Particle Approach for Continuum/Rarefied Flows , 1998 .

[35]  Graham V. Candler,et al.  The solution of the Navier-Stokes equations using Gauss-Seidel line relaxation , 1989 .

[36]  S. Senturia Microsystem Design , 2000 .

[37]  Patrick Le Tallec,et al.  Coupling Boltzmann and Navier-Stokes Equations by Half Fluxes , 1997 .

[38]  Hassan Hassan,et al.  A decoupled DSMC/Navier-Stokes analysis of a transitional flow experiment , 1996 .

[39]  D. C. Wadsworth,et al.  One-dimensional hybrid continuum/particle simulation approach for rarefied hypersonic flows , 1990 .

[40]  Eirik Grude Flekkøy,et al.  Hybrid model for combined particle and continuum dynamics , 2000 .

[41]  Charles Hirsch,et al.  Numerical computation of internal & external flows: fundamentals of numerical discretization , 1988 .

[42]  Forrest E. Lumpkin,et al.  Enhanced analyses of plume impingement during Shuttle-Mir docking using a combined CFD and DSMC methodology , 1996 .