Variant of gas kinetic flux solver for flows beyond Navier-Stokes level

In this paper, a variant of gas kinetic flux solver (GKFS) is presented for simulation of flows beyond the Navier-Stokes (NS) level. The method retains the framework of GKFS and reconstructs the numerical fluxes by the moments of distribution function at the cell interface, which is given from the local solution of the Boltzmann equation. In the conventional GKFS, the first-order Chapman-Enskog (CE) expansion is utilized to approximate the initial distribution function. By using the differential chain rule, it was found that the CE expansion form could be linked to the stress tensor and the heat flux. For flows in the NS level, the stress tensor and heat flux can be simply calculated from the linearized constitutive relationship and Fourier's law, respectively. However, for flows beyond the NS level, due to the strong nonequilibrium effect, the linearized constitutive relationship and Fourier's law are insufficient to predict the stress tensor and the heat flux. To overcome this difficulty, this paper introduces correction terms to the stress tensor and heat flux in the initial distribution function. These correction terms will take effect in the strong nonequilibrium region for flows beyond the NS level. To avoid finding complex expressions or solving complicated partial differential equations for the correction terms, a simple and iterative procedure is proposed to update the correction terms based on the framework of GKFS. The proposed method is validated by three benchmark cases which cover the flow from the continuum regime to the transition regime. Numerical results show that the present solver can provide accurate solution in the continuum regime. It is indeed the correction terms that take effect in the strong nonequilibrium region for flows beyond the NS level, which enables the present solver to capture the nonequilibrium phenomenon with reasonable accuracy for rarefied flows at moderate Knudsen number.