Implicit Discontinuous Galerkin Method for the Boltzmann Equation (1901.01865v1)

Abstract: An implicit high-order discontinuous Galerkin (DG) method is developed to find steady-state solution of rarefied gas flow described by the Boltzmann equation with full collision operator. In the physical space, velocity distribution function is approximated by the piecewise polynomials of degree up to 4, while in the velocity space the fast spectral method is incorporated into the DG to evaluate the collision operator. A specific polynomial approximation for the collision operator is proposed to reduce the computational complexity of the fast spectral method by $K$ times, where for two-dimensional problems $K$ is 15 when DG with 4th-order polynomials are 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 flow when the flow is fully resolved. Moreover, without assembling large sparse linear system, the computational cost in terms of memory and CPU time can be significantly reduced. Five different one/two-dimensional tests including low-speed microscale flows and hypersonic rarefied gas flows are used to verify the proposed approach. Our results show that, DG schemes of different order of approximating polynomial require the same number of iterative steps to obtain the steady-state solution with the same order of accuracy; and the higher order the scheme, fewer spatial elements thus less CPU time is needed. 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.