In this paper, we propose an efficient, high order accurate and asymptotic-preserving (AP) semi-Lagrangian (SL) method for the BGK model with constant or spatially dependent Knudsen number. The spatial discretization is performed by a mass conservative nodal discontinuous Galerkin (NDG) method, while the temporal discretization of the stiff relaxation term is realized by stiffly accurate diagonally implicit Runge-Kutta (DIRK) methods along characteristics. Extra order conditions are enforced for asymptotic accuracy (AA) property of DIRK methods when they are coupled with a semi-Lagrangian algorithm in solving the BGK model. A local maximum principle preserving (LMPP) limiter is added to control numerical oscillations in the transport step. Thanks to the SL and implicit nature of time discretization, the time stepping constraint is relaxed and it is much larger than that from an Eulerian framework with explicit treatment of the source term. Extensive numerical tests are presented to verify the high order AA, efficiency and shock capturing properties of the proposed schemes.
翻译:在本文中,我们建议对BGK模型采用高效、高顺序、准确和无症状保存半Lagrangian(SL)方法,该模型具有恒定或空间依赖的Knudsen编号。空间分解方法由大规模保守保守的节点不连续的Galerkin(NDG)方法进行,而僵硬的放松期的暂时分解则通过精确准确的对角隐含的龙格-库塔(DIRK)方法及其特点来实现。对于DIRK方法的无症状精确性(AAA)实施额外命令条件,如果这些方法与解决BGK模型的半Lagrangian算法相结合。在控制运输步骤中的数字分解时增加了一个本地最大原则保存限制(LMPP) 。由于SLL和时间分解的隐含性质,时间缓冲限制松了,从明确处理源术语的Eulerian框架中要大得多。为了核实高顺序的AAAA、效率和冲击性,进行了广泛的数字测试。