A parallel Runge–Kutta discontinuous Galerkin solver for rarefied gas flows based on 2D Boltzmann kinetic equations

Abstract The high-order Runge–Kutta Discontinuous Galerkin (RKDG) method is applied to solve the 2D Boltzmann kinetic equations. A conservative DG type discretization of the non-linear collision relaxation term is formulated for both the Bhatnagar–Gross–Krook and the ellipsoidal statistical kinetic models. Verification is carried out for a steady and an unsteady oscillatory 1D Couette flows, a 2D conduction problem as well as for a 2D long microchannel flow by comparison with the DSMC and analytical solutions. The computational performance of the RKDG method is compared with a widely used second-order finite volume method. The RKDG method has up to 3rd-order spatial accuracy and up to 4th-order time accuracy and is more efficient than the finite volume approach. The parallelization by domain decomposition in physical space is implemented and parallel performance is evaluated. It is shown that 2nd order RKDG is over 15 times faster than the 2nd-order FVM method for the Couette flow test case. The high-order RKDG method is especially attractive for solution of low-speed and unsteady rarefied flows.

[1]  John E. Sader,et al.  High accuracy numerical solutions of the Boltzmann Bhatnagar-Gross-Krook equation for steady and oscillatory Couette flows , 2012 .

[2]  Chi-Wang Shu,et al.  A numerical study for the performance of the Runge-Kutta discontinuous Galerkin method based on different numerical fluxes , 2006, J. Comput. Phys..

[3]  Juan-Chen Huang,et al.  Rarefied Flow Computations Using Nonlinear Model Boltzmann Equations , 1995 .

[4]  Frank P. Incropera,et al.  Fundamentals of Heat and Mass Transfer , 1981 .

[5]  Jonathan M. Burt,et al.  A low diffusion particle method for simulating compressible inviscid flows , 2007, J. Comput. Phys..

[6]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[7]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[8]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[9]  N. Hadjiconstantinou.  Oscillatory shear-driven gas flows in the transition and free-molecular-flow regimes , 2005 .

[10]  Felix Sharipov,et al.  Data on Internal Rarefied Gas Flows , 1998 .

[11]  Chi-Wang Shu,et al.  Numerical Comparison of WENO Finite Volume and Runge–Kutta Discontinuous Galerkin Methods , 2001, J. Sci. Comput..

[12]  Matthias K. Gobbert,et al.  A Galerkin Method for the Simulation of the Transient 2-D/2-D and 3-D/3-D Linear Boltzmann Equation , 2007, J. Sci. Comput..

[13]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

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

[15]  Vladimir Titarev,et al.  Efficient Deterministic Modelling of Three-Dimensional Rarefied Gas Flows , 2012 .

[16]  Lorenzo Pareschi,et al.  An Implicit Monte Carlo Method for Rarefied Gas Dynamics , 1999 .

[17]  Sergey F. Gimelshein,et al.  Numerical Modeling of Axisymmetric and Three-Dimensional Flows in Microelectromechanical Systems Nozzles , 2002 .

[18]  V. V. Aristov,et al.  Construction of a Unified Continuum/Kinetic Solver for Aerodynamic Problems , 2005 .

[19]  D. P. Giddens,et al.  The discrete ordinate method for the linearized boundary value problems in kinetic theory of gases. , 1967 .

[20]  Clint Dawson,et al.  Time step restrictions for Runge-Kutta discontinuous Galerkin methods on triangular grids , 2008, J. Comput. Phys..

[21]  Michael N. Macrossan,et al.  Hypersonic Flow over a Wedge with a Particle Flux Method , 2005 .

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

[23]  Juan-Chen Huang.  A conservative discrete ordinate method for model Boltzmann equations , 2011 .

[24]  Alina Alexeenko,et al.  High-Order Discontinuous Galerkin Method for Boltzmann Model Equations , 2008 .

[25]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[26]  Luc Mieussens,et al.  Numerical comparison of Bhatnagar–Gross–Krook models with proper Prandtl number , 2004 .

[27]  Francis J. Alexander,et al.  Cell size dependence of transport coefficients in stochastic particle algorithms , 1998 .

[28]  K. Breuer,et al.  Gaseous slip flow in long microchannels , 1997 .

[29]  Jonathan M. Burt,et al.  A hybrid particle approach for continuum and rarefied flow simulation , 2009, J. Comput. Phys..

[30]  V. V. Aristov,et al.  Unified solver for rarefied and continuum flows with adaptive mesh and algorithm refinement , 2007, J. Comput. Phys..

[31]  Lowell L. Baker,et al.  Variance reduction for Monte Carlo solutions of the Boltzmann equation , 2005 .

[32]  Lowell L. Baker,et al.  Variance‐reduced Monte Carlo solutions of the Boltzmann equation for low‐speed gas flows: A discontinuous Galerkin formulation , 2008 .

[33]  Matthias Laux.  Local Time Stepping with Automatic Adaptation for the DSMC Method , 1998 .

[34]  A. Alekseenko,et al.  Numerical properties of high order discrete velocity solutions to the BGK kinetic equation , 2011 .

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

[36]  Wim Desmet,et al.  CFL Conditions for Runge-Kutta discontinuous Galerkin methods on triangular grids , 2011, J. Comput. Phys..

[37]  L. Mieussens.  Discrete-Velocity Models and Numerical Schemes for the Boltzmann-BGK Equation in Plane and Axisymmetric Geometries , 2000 .

[38]  Nicolas G. Hadjiconstantinou,et al.  Low-noise Monte Carlo simulation of the variable hard sphere gas , 2011 .

[39]  Alina A. Alexeenko,et al.  Numerical and experimental study of gas flows in 2D and 3D microchannels , 2008 .

[40]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[41]  G. R. Cowper,et al.  Gaussian quadrature formulas for triangles , 1973 .

[42]  B. Perthame,et al.  The Gaussian-BGK model of Boltzmann equation with small Prandtl number , 2000 .

[43]  Sruti Chigullapalli,et al.  Entropy considerations in numerical simulations of non-equilibrium rarefied flows , 2010, J. Comput. Phys..

[44]  Koji Morinishi,et al.  Numerical simulation for gas microflows using Boltzmann equation , 2006 .

[45]  Lowell H. Holway,et al.  New Statistical Models for Kinetic Theory: Methods of Construction , 1966 .

[46]  Zhi-Hui Li,et al.  Gas kinetic algorithm using Boltzmann model equation , 2004 .

[47]  Wolfgang Wagner,et al.  Deviational particle Monte Carlo for the Boltzmann equation , 2008, Monte Carlo Methods Appl..

[48]  Hongwei Liu,et al.  A Runge-Kutta discontinuous Galerkin method for viscous flow equations , 2007, J. Comput. Phys..

[49]  S. Gimelshein,et al.  Reconsideration of low Reynolds number flow-through constriction microchannels using the DSMC method , 2005, Journal of Microelectromechanical Systems.