Implicit discontinuous Galerkin method for the Boltzmann equation

Wei Su, Peng Wang, Yonghao Zhang, Lei Wu*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

19 Citations (Scopus)
83 Downloads (Pure)

Abstract

An implicit high-order discontinuous Galerkin (DG) method is developed to find the steady-state solution of rarefied gas flow described by the Boltzmann equation with full collision operator. In the physical space, the velocity distribution function is approximated by the piecewise polynomial of degree up to 4, while in the velocity space the fast spectral method is incorporated into the DG discretization to evaluate the collision operator. A specific polynomial approximation for the collision operator is proposed to reduce the computational cost of the fast spectral method by K times, where for two-dimensional problem K is 15 when the DG with 4th-order polynomial is used on triangular mesh. Based on the first-order upwind scheme, a sweeping technique is employed to solve the local linear equations resulting from the DG discretization sequentially over spatial elements. This technique can preserve stability of the scheme and requires no nonlinear limiter in solving hypersonic rarefied gas flows when the shock wave structure is fully resolved by fine spatial grid. Moreover, without assembling large sparse linear system, the computational cost in terms of memory consumption and CPU time is significantly reduced. Five different one/two-dimensional tests including low-speed microscale flows and hypersonic rarefied gas flows are used to assess the accuracy and efficiency of proposed approach. Our results show that, DG scheme of different order of approximating polynomial requires the same number of iterative steps to obtain the steady-state solution with the same order of accuracy; and the higher order the scheme is, the fewer spatial elements are needed, thus leading to less CPU time. Besides, our method can be faster than the finite difference solver by about one order of magnitude. The produced solutions can be used as benchmark data for assessing the accuracy of other gas kinetic solvers for the Boltzmann equation and gas kinetic models that simplify the Boltzmann collision operator.

Original languageEnglish
Article number39
JournalJournal of Scientific Computing
Volume82
Issue number2
DOIs
Publication statusPublished - 3 Feb 2020

Keywords

  • Boltzmann equation
  • discontinuous Galerkin
  • fast spectral method
  • high-order discretization
  • implicit scheme

Fingerprint

Dive into the research topics of 'Implicit discontinuous Galerkin method for the Boltzmann equation'. Together they form a unique fingerprint.

Cite this