Modeling and computation for non-equilibrium gas dynamics: beyond kinetic relaxation model
The non-equilibrium gas dynamics is described by the Boltzmann equation, which can be solved numerically through the deterministic and stochastic methods. Due to the complicated collision term of the Boltzmann equation, many kinetic relaxation models have been proposed and used in the past seventy years for the study of rarefied flow. In order to develop a multiscale method for the rarefied and continuum flow simulation, by adopting the integral solution of the kinetic model equation a DVM-type unified gas-kinetic scheme (UGKS) has been constructed. The UGKS models the gas dynamics on the cell size and time step scales while the accumulating effect from particle transport and collision has been taken into account within a time step. Under the UGKS framework, a unified gas-kinetic wave-particle (UGKWP) method has been further developed for non-equilibrium flow simulation, where the time evolution of gas distribution function is composed of analytical wave and individual particle. In the highly rarefied regime, particle transport and collision will play a dominant role. Due to the single relaxation time model for particle collision, there is a noticeable discrepancy between the UGKWP solution and the full Boltzmann or DSMC result, especially in the high Mach and Knudsen number cases. In this paper, besides the kinetic relaxation model, a modification of particle collision time according to the particle velocity will be implemented in UGKWP. As a result, the new model greatly improves the performance of UGKWP in the capturing of non-equilibrium flow. There is a perfect match between UGKWP and DSMC or Boltzmann solution in the highly rarefied regime. In the near continuum and continuum flow regime, the UGKWP will gradually get back to the macroscopic variables based Navier-Stokes flow solver at small cell Knudsen number.