2025-11-10T02:36:53.335691

Asymptotic-preserving semi-Lagrangian discontinuous Galerkin schemes for the Boltzmann equation

Cai, Hao, Liu et al.

In this work, we present an asymptotic-preserving semi-Lagrangian discontinuous Galerkin scheme for the Boltzmann equation that effectively handles multi-scale transport phenomena. The main challenge lies in designing appropriate moments update for penalization within the semi-Lagrangian framework. Inspired by [M. Ding, J. M. Qiu, and R. Shu, Multiscale Model. Simul. 21 (2023), no. 1, 143--167], the key ingredient is utilizing the Shu-Osher form of the scheme in the implicit-explicit Runge-Kutta (IMEX-RK) setting, which enables us to capture the correct limiting system by constructing an appropriate moments update procedure. Our theoretical analysis establishes accuracy order conditions for both the IMEX-RK time integration and the new moments update step. We also employ hypocoercivity techniques to establish stability for the linearized model. Numerical experiments for various test problems validate our proposed scheme's accuracy, asymptotic-preserving property, and robustness in various regimes, which demonstrates its effectiveness for multi-scale kinetic simulations.

academic

Asymptotic-preserving semi-Lagrangian discontinuous Galerkin schemes for the Boltzmann equation

基本信息

论文ID: 2510.14375
标题: Asymptotic-preserving semi-Lagrangian discontinuous Galerkin schemes for the Boltzmann equation
作者: Xiaofeng Cai, Zhen Hao, Liu Liu, Jiayu Wan
分类: math.NA cs.NA
发表时间: 2025年10月16日提交
论文链接: https://arxiv.org/abs/2510.14375

摘要

本文提出了一种针对Boltzmann方程的渐近保持半拉格朗日不连续伽辽金格式，能够有效处理多尺度输运现象。主要挑战在于在半拉格朗日框架内为惩罚方法设计适当的矩更新机制。受Ding等人工作启发，关键技术是在隐式-显式Runge-Kutta (IMEX-RK)设置中利用Shu-Osher形式，通过构造适当的矩更新过程来捕捉正确的极限系统。理论分析建立了IMEX-RK时间积分和新矩更新步骤的精度阶条件，并采用超强制性技术建立了线性化模型的稳定性。数值实验验证了所提格式在各种测试问题上的精度、渐近保持性质和鲁棒性。

研究背景与动机

问题背景

Boltzmann方程的重要性: Boltzmann方程是描述稀薄气体动力学的基础动理学方程，在连续介质模型如Navier-Stokes方程失效的情况下能够正确描述气体流动物理。
计算挑战:
- 高维相空间带来的计算复杂性
- 碰撞算子的非线性、非局部性质
- 多尺度特性导致的刚性问题
现有方法局限性:
- 传统欧拉方法受CFL条件限制，时间步长受限
- 在小Knudsen数ε情况下，显式方法需要∆t = O(ε)的严格时间步限制
- 半拉格朗日方法虽然避免CFL限制，但在设计渐近保持格式时面临新挑战

研究动机

开发能够在所有物理尺度下都保持精度和效率的数值格式，特别是在流体极限(ε→0)下能够自动退化到正确的欧拉方程。

核心贡献

创新的矩更新机制: 在半拉格朗日框架内设计了新的矩更新过程，解决了传统欧拉方法中矩更新策略在SL框架下失效的问题
Shu-Osher形式的应用: 将Shu-Osher形式适配到IMEX-RK设置中，使得能够正确识别极限系统并构造相应的矩更新
理论分析:
- 建立了IMEX-RK方案的精度阶条件
- 证明了渐近保持性质
- 利用超强制性技术分析了线性化模型的稳定性
- 给出了保正性的充分条件
数值验证: 通过多个测试问题验证了格式的精度、渐近保持性质和鲁棒性

方法详解

任务定义

求解Boltzmann方程： $\partial_t f + v \cdot \nabla_x f = \frac{1}{\varepsilon} Q(f)$

其中f(x,v,t)是粒子分布函数，ε是Knudsen数，Q(f)是碰撞算子。目标是设计在所有ε值下都高效准确的数值格式。

模型架构

1. BGK惩罚方法

采用BGK惩罚分解： $Q(f) = Q(f) - \beta Q_{BGK}(f) + \beta Q_{BGK}(f)$ 其中 $Q_{BGK}(f) = M_f - f$ ， $M_f$ 是与f相关的Maxwellian分布。

2. 半拉格朗日离散化

基于特征Galerkin弱形式： $\int_{I} f(t_2,x,v)\Psi(x)dx = \int_{I-v\Delta t} f_0(x,v)\Psi(x+v\Delta t)dx + \text{碰撞项积分}$

3. 节点DG空间离散

使用Lagrange多项式基函数： $f(x,v) = \sum_{j=1}^{N_x} \sum_{p_j=0}^k f(x_{j,p_j},v)\ell_{j,p_j}(x)$

4. IMEX-RK时间离散

s阶IMEX-RK格式： $f^{(i)} = \tilde{S}_{i,0}[f^n] + \Delta t \sum_{j=1}^{i-1} a_{ij}S_{i,j}[\frac{1}{\varepsilon}G_P(f^{(j)})] + \Delta t \sum_{j=1}^i \tilde{a}_{ij}\tilde{S}_{i,j}[\frac{1}{\varepsilon}Q_P(f^{(j)})]$

技术创新点

1. 新的矩更新策略

传统欧拉方法中，矩更新通过积分消除刚性项： $U^{n+1} = U^n - \langle v \cdot \nabla_x f^n, \phi(v) \rangle$

但在SL框架中，由于算子 $S_{v,\Delta t}$ 的存在，刚性项无法消除。本文提出基于极限格式的矩更新： $U^{(i)} = (1-\bar{\tilde{A}}_{i-1}\hat{\tilde{A}}^{-1}_{(i-2)}e_{(i)})\langle S_{i,0}[f^n], \phi(v)\rangle + \bar{\tilde{A}}_{i-1}\hat{\tilde{A}}^{-1}_{(i-2)}\langle \tilde{S}^i[\hat{F}_{i-1}], \phi(v)\rangle$