A velocity-space adaptive unified gas kinetic scheme for continuum and rarefied flows

被引:28
作者
Xiao, Tianbai [1 ,2 ]
Liu, Chang [3 ]
Xu, Kun [3 ,4 ,5 ]
Cai, Qingdong [6 ]
机构
[1] Karlsruhe Inst Technol, Dept Math, D-76131 Karlsruhe, Germany
[2] Karlsruhe Inst Technol, Steinbuch Ctr Comp, D-76131 Karlsruhe, Germany
[3] Hong Kong Univ Sci & Technol, Dept Math, Kowloon, Clear Water Bay, Hong Kong, Peoples R China
[4] Hong Kong Univ Sci & Technol, Dept Mech & Aerosp Engn, Kowloon, Clear Water Bay, Hong Kong, Peoples R China
[5] Hong Kong Univ Sci & Technol, Shenzhen Res Inst, Shenzhen 518057, Peoples R China
[6] Peking Univ, Coll Engn, Dept Mech & Engn Sci, Beijing 100871, Peoples R China
基金
美国国家科学基金会;
关键词
Unified gas kinetic scheme; Gas kinetic scheme; Multiscale flow; Non-equilibrium phenomena; Adaptive velocity space; NAVIER-STOKES EQUATIONS; BOLTZMANN; SIMULATIONS; TRANSPORT;
D O I
10.1016/j.jcp.2020.109535
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Compressible flow has intrinsically multiple scale nature due to the large variations of gas density and characteristic scale of flow structure, especially in hypersonic and reentry problems. It is challenging to construct an accurate and efficient numerical algorithm to capture non-equilibrium flow physics across different regimes. In this paper, a unified gas kinetic scheme with adaptive velocity space (AUGKS) for multiscale flow transport will be developed. In near-equilibrium flow region, particle distribution function is close to the Chapman-Enskog expansion and can be formulated analytically, where only macroscopic conservative flow variables are updated. With the emergence of non-equilibrium effect, the AUGKS automatically switches to the original unified gas kinetic scheme (UGKS) with a discrete velocity space to follow the evolution of particle distribution function. A criterion is proposed to quantify the non-equilibrium and is used for the switch between continuous and discrete particle velocity space. Following the scale-dependent local evolving solution, the AUGKS presents the discretized gas dynamic equations directly on the cell size and time step scales, i.e., the so-called direct modeling method. As a result, the scheme is able to capture the cross-scale flow physics from particle transport to hydrodynamic wave propagation, and provides a continuous variation of solutions from the Boltzmann to the Euler equations. Different from conventional DSMC-Fluid hybrid method, the AUGKS does not need a buffer zone to match up kinetic and hydrodynamic solutions. Instead, continuous and discrete particle velocity spaces are automatically and robustly switched, such as translating the continuous Chapman-Enskog distribution function to the discrete grid points in the velocity space. Therefore, the AUGKS is feasible for the numerical simulations with unsteadiness and complex geometries. Compared with the asymptotic preserving (AP) method which solves kinetic equation uniformly over entire computational domain with discretized velocity space, the current velocity-space adaptive unified scheme speeds up the computation, reduces the memory requirement significantly, and maintains the equivalent accuracy for multiscale flow simulations. Many test cases validate the current approach. The AUGKS provides an effective tool for non-equilibrium flow studies. (C) 2020 Elsevier Inc. All rights reserved.
引用
收藏
页数:35
相关论文
共 44 条
[21]   COMPACTNESS IN BOLTZMANNS EQUATION VIA FOURIER INTEGRAL-OPERATORS AND APPLICATIONS .1. [J].
LIONS, PL .
JOURNAL OF MATHEMATICS OF KYOTO UNIVERSITY, 1994, 34 (02) :391-427
[22]   A unified gas-kinetic scheme for continuum and rarefied flows IV: Full Boltzmann and model equations [J].
Liu, Chang ;
Xu, Kun ;
Sun, Quanhua ;
Cai, Qingdong .
JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 314 :305-340
[23]   Fast algorithms for computing the Boltzmann collision operator [J].
Mouhot, Clement ;
Pareschi, Lorenzo .
MATHEMATICS OF COMPUTATION, 2006, 75 (256) :1833-1852
[24]   The kinetic scheme for the full-Burnett equations [J].
Ohwada, T ;
Xu, K .
JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 201 (01) :315-332
[25]   Hybrid Euler/direct simulation Monte Carlo calculation of unsteady slit flow [J].
Roveda, R ;
Goldstein, DB ;
Varghese, PL .
JOURNAL OF SPACECRAFT AND ROCKETS, 2000, 37 (06) :753-760
[26]   Hybrid Euler/particle approach for continuum/rarefied flows [J].
Roveda, R ;
Goldstein, DB ;
Varghese, PL .
JOURNAL OF SPACECRAFT AND ROCKETS, 1998, 35 (03) :258-265
[27]   A modular particle-continuum numerical method for hypersonic non-equilibrium gas flows [J].
Schwartzentruber, T. E. ;
Scalabrin, L. C. ;
Boyd, I. D. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 225 (01) :1159-1174
[28]  
Shakhov E.M., 1972, Fluid Dyn., V3, P95, DOI DOI 10.1007/BF01029546
[29]  
Sun Q., 2003, Information preservation methods for modeling micro-scale gas flows
[30]   An adaptive domain decomposition procedure for Boltzmann and Euler equations [J].
Tiwari, S ;
Klar, A .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 1998, 90 (02) :223-237