  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.
