High Order Conservative Semi-Lagrangian Scheme for the BGK Model of the Boltzmann Equation
作者机构：Department of Mathematics and Computer ScienceUniversity of Catania95125 CataniaItaly Department of MathematicsSungkyunkwan UniversitySuwon 440-746Republic of Korea
出 版 物：《Communications in Computational Physics》 (计算物理通讯（英文）)
年 卷 期：2021年第29卷第1期
页 面：1-56页
核心收录：
学科分类：07[理学] 0704[理学-天文学] 0701[理学-数学] 0702[理学-物理学] 070101[理学-基础数学]
基 金：supported by Samsung Science and Technology Foundation under Project Number SSTF-BA1801-02
主 题：BGK model Boltzmann equation semi-Lagrangian scheme conservative correction discrete Maxwellian
摘 要：In this paper,we present a conservative semi-Lagrangian finite-difference scheme for the BGK *** semi-Lagrangian finite difference schemes,coupled with an L-stable treatment of the collision term,allow large time steps,for all the range of Knudsen number[17,27,30].Unfortunately,however,such schemes are not *** of conservation is analyzed in detail,and two main sources are identified as its ***,when using classical continuous Maxwellian,conservation error is negligible only if velocity space is resolved with sufficiently large number of grid ***,for a small number of grid points in velocity space such error is not negligible,because the parameters of the Maxwellian do not coincide with the discrete ***,the non-linear reconstruction used to prevent oscillations destroys the translation invariance which is at the basis of the conservation properties of the *** a consequence the schemes show a wrong shock speed in the limit of small Knudsen *** treat the first problem and ensure machine precision conservation of mass,momentum and energy with a relatively small number of velocity grid points,we replace the continuous Maxwellian with the discrete Maxwellian introduced in[22].The second problem is treated by implementing a conservative correction procedure based on the flux difference form as in[26].In this way we can construct conservative semi-Lagrangian schemes which are Asymptotic Preserving(AP)for the underlying Euler limit,as the Knudsen number *** effectiveness of the proposed scheme is demonstrated by extensive numerical tests.
维普期刊数据库