A simple distribution function-based gas-kinetic scheme for simulation of viscous incompressible and compressible flows

• Published in: Journal of computational physics 2014-10, Vol.274, p.611-632
• Main Authors: Yang, L.M., Shu, C., Wu, J.
• Format: Article
• Language: English
• Subjects: Circular function; Circularity; Computational fluid dynamics; Computer simulation; Fluid flow; Gas-kinetic scheme; Mathematical analysis; Mathematical models; Maxwellian distribution function; Navier-Stokes equations; Viscous flow; Viscous flows
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2014.06.033

In this work, a simple distribution function-based gas-kinetic scheme for simulation of viscous flows is presented. The work applies the finite volume method to discretize the governing differential equations, and inviscid and viscous fluxes at the cell interface are evaluated simultaneously by local reconstruction of solution for the continuous Boltzmann equation. Differently from the conventional gas-kinetic scheme [13–15], in the present work, the Maxwellian distribution function is simplified by a simple distribution function, and integrals in the infinity domain of phase space are reduced to integrals around a circle. As a consequence, the computational efficiency is greatly improved. Since the simple distribution function is defined on the circle, for simplicity, it is termed as circular function hereafter. The present work is the extension of our previous work [20], where the circular function-based gas-kinetic scheme is presented to simulate inviscid flows. Only the equilibrium distribution function is considered in [20]. To solve viscous flows, the non-equilibrium part of density distribution function has to be considered. One of major contributions in this work is to present a simple way to compute the non-equilibrium part of the distribution function. It can be calculated by the difference of equilibrium distribution functions at the cell interface and its surrounding point. As a result, the formulations for computing the conservative flow variables and fluxes at the cell interface can be given explicitly. The present solver can simulate both incompressible and compressible viscous flows. To validate the proposed new gas-kinetic scheme, several incompressible and compressible viscous flows are simulated. Numerical results showed that the circular function-based gas-kinetic scheme can provide accurate numerical results with the same computational cost as that needed by conventional Navier–Stokes solver.