A Discontinuous Galerkin Fast Spectral Method for the Full Boltzmann Equation with General Collision Kernels

Abstract The Boltzmann equation, an integro-differential equation for the molecular distribution function in the physical and velocity phase space, governs the fluid flow behavior at a wide range of physical conditions, including compressible, turbulent, as well as flows involving further physics such as non-equilibrium internal energy exchange and chemical reactions. Despite its wide applicability, deterministic solution of the Boltzmann equation presents a huge computational challenge, and often the collision operator is simplified for practical reasons. In this work, we introduce a highly accurate deterministic method for the full Boltzmann equation which couples the Runge–Kutta discontinuous Galerkin (RKDG) discretization in time and physical space (Su et al. (2015) [1] ) and the recently developed fast Fourier spectral method in velocity space (Gamba et al. (2017) [2] ). The novelty of this approach encompasses three aspects: first, the fast spectral method for the collision operator applies to general collision kernels with little or no practical limitations, and in order to adapt to the spatial discretization, we propose here a singular-value-decomposition based algorithm to further reduce the cost in evaluating the collision term; second, the DG formulation employed has high order of accuracy at element-level, and has shown to be more efficient than the finite volume method; thirdly, the element-local compact nature of DG as well as our collision algorithm is amenable to effective parallelization on massively parallel architectures. The solver has been verified against analytical Bobylev–Krook–Wu solution. Further, the standard benchmark test cases of rarefied Fourier heat transfer, Couette flow, oscillatory Couette flow, normal shock wave, lid-driven cavity flow, and thermally driven cavity flow have been studied and their results are compared against direct simulation Monte Carlo (DSMC) solutions with equivalent molecular collision models or published deterministic solutions.

[1]  Alina A. Alexeenko,et al.  Binary scattering model for Lennard-Jones potential: Transport coefficients and collision integrals for non-equilibrium gas flow simulations , 2012 .

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

[3]  Timothy J. Barth,et al.  On Discontinuous Galerkin Approximations of Boltzmann Moment Systems With Levermore Closure , 2006 .

[4]  G. Karniadakis,et al.  Microflows and Nanoflows: Fundamentals and Simulation , 2001 .

[5]  K. Stephani,et al.  Recommended direct simulation Monte Carlo collision model parameters for modeling ionized air transport processes , 2016 .

[6]  J. Torczynski,et al.  Direct simulation Monte Carlo investigation of the Rayleigh-Taylor instability , 2016 .

[7]  Graeme A. Bird,et al.  The DSMC Method , 2013 .

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

[9]  Lei Wu,et al.  Influence of intermolecular potentials on rarefied gas flows: Fast spectral solutions of the Boltzmann equation , 2015, Physics of Fluids.

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

[11]  Alejandro L. Garcia,et al.  Direct Simulation Monte Carlo , 1997 .

[12]  J. R. Torczynski,et al.  Normal solutions of the Boltzmann equation for highly nonequilibrium Fourier flow and Couette flow , 2006 .

[13]  Lei Wu,et al.  Solving the Boltzmann equation deterministically by the fast spectral method: application to gas microflows , 2014, Journal of Fluid Mechanics.

[14]  E. H. van Brummelen,et al.  An entropy stable discontinuous Galerkin finite-element moment method for the Boltzmann equation , 2016, Comput. Math. Appl..

[15]  E. Atlee Jackson,et al.  Kinetic Models and the Linearized Boltzmann Equation , 1959 .

[16]  Luc Mieussens,et al.  A survey of deterministic solvers for rarefied flows (Invited) , 2014 .

[17]  Graeme A. Bird,et al.  Approach to Translational Equilibrium in a Rigid Sphere Gas , 1963 .

[18]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[19]  Lorenzo Pareschi,et al.  Numerical Solution of the Boltzmann Equation I: Spectrally Accurate Approximation of the Collision Operator , 2000, SIAM J. Numer. Anal..

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

[21]  Giacomo Dimarco,et al.  Hybrid Multiscale Methods II. Kinetic Equations , 2008, Multiscale Model. Simul..

[22]  Zhenyuan Tang,et al.  Stable Runge-Kutta discontinuous Galerkin solver for hypersonic rarefied gaseous flow based on 2D Boltzmann kinetic model equations , 2017 .

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

[24]  Yoshio Sone,et al.  Kinetic Theory and Fluid Dynamics , 2002 .

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

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

[27]  Lei Wu,et al.  Assessment and development of the gas kinetic boundary condition for the Boltzmann equation , 2017, Journal of Fluid Mechanics.

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

[29]  Alina Alexeenko,et al.  Revised Variable Soft Sphere and Lennard-Jones Model Parameters for Eight Common Gases up to 2200 K , 2015 .

[30]  Daniel J. Rader,et al.  Direct simulation Monte Carlo: The quest for speed , 2014 .

[31]  A. Bobylev.  Exact solutions of the Boltzmann equation , 1975 .

[32]  Haihu Liu,et al.  A fast iterative scheme for the linearized Boltzmann equation , 2016, J. Comput. Phys..

[33]  Yonghao Zhang,et al.  Deterministic numerical solutions of the Boltzmann equation using the fast spectral method , 2013, J. Comput. Phys..

[34]  Stefan K. Stefanov,et al.  On DSMC Calculations of Rarefied Gas Flows with Small Number of Particles in Cells , 2011, SIAM J. Sci. Comput..

[35]  Giacomo Dimarco,et al.  Numerical methods for kinetic equations* , 2014, Acta Numerica.

[36]  J. Torczynski,et al.  Molecular-Level Simulations of Turbulence and Its Decay. , 2017, Physical review letters.

[37]  Alina A. Alexeenko,et al.  Rarefaction effects for the transonic airfoils in low Reynolds number regime , 2016 .

[38]  Lorenzo Pareschi,et al.  Fast algorithms for computing the Boltzmann collision operator , 2006, Math. Comput..

[39]  C. Cercignani.  The Boltzmann equation and its applications , 1988 .

[40]  H. Grad.  On the kinetic theory of rarefied gases , 1949 .

[41]  Quanhua Sun,et al.  A Hybrid Continuum / Particle Approach for Micro‐Scale Gas Flows , 2003 .

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

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

[44]  David R. Emerson,et al.  Investigation of Heat and Mass Transfer in a Lid-Driven Cavity Under Nonequilibrium Flow Conditions , 2010 .

[45]  Ralf Dresner,et al.  Theoretical Mechanics Of Particles And Continua , 2016 .

[46]  David I. Ketcheson,et al.  Strong stability preserving runge-kutta and multistep time discretizations , 2011 .

[47]  Taku Ohwada,et al.  Structure of normal shock waves: Direct numerical analysis of the Boltzmann equation for hard-sphere molecules , 1993 .

[48]  J. R. Torczynski,et al.  Convergence behavior of a new DSMC algorithm , 2009, J. Comput. Phys..

[49]  Alina A. Alexeenko,et al.  Direct simulation Monte Carlo modeling of H2–O2 deflagration waves , 2018, Combustion and Flame.

[50]  C. Beentjes,et al.  QUADRATURE ON A SPHERICAL SURFACE , 2016 .

[52]  David Emerson,et al.  A high-order moment approach for capturing non-equilibrium phenomena in the transition regime , 2009, Journal of Fluid Mechanics.

[53]  Bernardo Cockburn,et al.  The Runge-Kutta local projection discontinous Galerkin finite element method for conservation laws , 1990 .

[54]  J. Torczynski,et al.  Calculations of the near-wall thermophoretic force in rarefied gas flow , 2002 .

[55]  Alina Alexeenko,et al.  A parallel Runge–Kutta discontinuous Galerkin solver for rarefied gas flows based on 2D Boltzmann kinetic equations , 2015 .

[56]  M. N. Kogan,et al.  Stresses produced in gasses by temperature and concentration inhomogeneities. New types of free convection , 1976 .

[57]  John E. Sader,et al.  Velocity gradient singularity and structure of the velocity profile in the Knudsen layer according to the Boltzmann equation. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[58]  H. Struchtrup,et al.  Regularization of Grad’s 13 moment equations: Derivation and linear analysis , 2003 .

[59]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[60]  Ruo Li,et al.  Globally Hyperbolic Regularization of Grad's Moment System , 2012 .

[61]  W. Wagner.  A convergence proof for Bird's direct simulation Monte Carlo method for the Boltzmann equation , 1992 .

[62]  C. D. Levermore,et al.  Moment closure hierarchies for kinetic theories , 1996 .

[63]  Jingwei Hu,et al.  A Fast Spectral Method for the Boltzmann Collision Operator with General Collision Kernels , 2016, SIAM J. Sci. Comput..

[64]  Hiroaki Matsumoto,et al.  Variable soft sphere molecular model for inverse-power-law or Lennard-Jones potential , 1991 .

[65]  Irene M. Gamba,et al.  Spectral-Lagrangian methods for collisional models of non-equilibrium statistical states , 2009, J. Comput. Phys..