Multiscale kinetic inviscid flux extracted from a gas-kinetic scheme for simulating incompressible and compressible flows

作     者：Sha Liu Junzhe Cao Chengwen Zhong 

作者机构：National Key Laboratory of Science and Technology on Aerodynamic Design and Research Northwestern Polytechnical University Xi'an Shaanxi 710072 China School of Aeronautics Northwestern Polytechnical University Xi'an Shaanxi 710072 China 

出 版 物：《Physical Review E》 (物理学评论E辑：统计、非线性和软体物理学)

年 卷 期：2020年第102卷第3期

页      面：033310-033310页

核心收录：

学科分类：07[理学] 070203[理学-原子与分子物理] 0702[理学-物理学] 

基　　金：111 Project of China National Numerical Wind-Tunnel Project of China National Natural Science Foundation of China, NSFC, (11702223, 11902264, 11902266) Higher Education Discipline Innovation Project, (B17037) National Numerical Wind Tunnel Project of China 

主　　题：Aerodynamics Compressible flows High-speed flow Kinetic theory Shock waves 

摘      要：A kinetic inviscid flux (KIF) is proposed for simulating incompressible and compressible flows. It is constructed based on the direct modeling of multiscale flow behaviors, which is used in the gas-kinetic scheme (GKS), the unified gas-kinetic scheme (UGKS), the discrete unified gas-kinetic scheme (DUGKS), etc. In KIF, the discontinuities (such as the shock wave) that cannot be well resolved by mesh cells are mainly solved by the kinetic flux vector splitting (KFVS) method representing the free transport mechanism (or microscale mechanism), while in other flow regions that are smooth, the flow behavior is solved mainly by the central-scheme-like totally thermalized transport (TTT). The weights of KFVS and TTT in KIF is automatically determined by those in the theory of direct modeling. Two ways of choosing the weights in KIF are proposed, which are actually the weights adopted in the UGKS and the DUGKS, respectively. By using the test cases of the Sod shock tube, the rarefaction wave, the boundary layer of a flat plate, the cavity flow, hypersonic flow over a circular cylinder, the shock and turbulent boundary iteration, and transonic flow over a three-dimensional M6 wing, the validity and accuracy of the present method are examined. The KIF does not suffer from the carbuncle phenomenon, and does not introduce extra numerical viscosity in smooth regions. Especially in the case of hypersonic cylinder, it gives quite sharp and clear density and temperature contours. The KIF can be viewed as an inviscid-viscous splitting version of the GKS. By doing this splitting, it is easy to be used in traditional computational fluid dynamics frameworks. It can also be classified as a type in the numerical schemes based on the kinetic theory that are represented by the works of Sun et al. [Adv. Appl. Math. Mech. 8, 703 (2016)] and Ohwada et al. [J. Comput. Phys. 362, 131 (2018)], except the weights are determined by the weights of direct modeling.