Discrete unified gas kinetic scheme for all Knudsen number flows. II. Thermal compressible case

This paper is a continuation of our work on the development of multiscale numerical scheme from low-speed isothermal flow to compressible flows at high Mach numbers. In our earlier work [Z.L. Guo et al., Phys. Rev. E 88, 033305 (2013)], a discrete unified gas kinetic scheme (DUGKS) was developed for low-speed flows in which the Mach number is small so that the flow is nearly incompressible. In the current work, we extend the scheme to compressible flows with the inclusion of thermal effect and shock discontinuity based on the gas kinetic Shakhov model. This method is an explicit finite-volume scheme with the coupling of particle transport and collision in the flux evaluation at a cell interface. As a result, the time step of the method is not limited by the particle collision time. With the variation of the ratio between the time step and particle collision time, the scheme is an asymptotic preserving (AP) method, where both the Chapman-Enskog expansion for the Navier-Stokes solution in the continuum regime and the free transport mechanism in the rarefied limit can be precisely recovered with a second-order accuracy in both space and time. The DUGKS is an idealized multiscale method for all Knudsen number flow simulations. A number of numerical tests, including the shock structure problem, the Sod tube problem in a whole range of degree of rarefaction, and the two-dimensional Riemann problem in both continuum and rarefied regimes, are performed to validate the scheme. Comparisons with the results of direct simulation Monte Carlo (DSMC) and other benchmark data demonstrate that the DUGKS is a reliable and efficient method for multiscale flow problems.