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

ボルツマン方程式に対する漸近保存半ラグランジュ不連続ガラーキン格式

基本情報

論文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

要旨

本論文は、ボルツマン方程式に対する漸近保存半ラグランジュ不連続ガラーキン格式を提案し、多尺度輸送現象を効果的に扱うことができる。主な課題は、半ラグランジュ枠組み内でペナルティ法に対する適切なモーメント更新機構を設計することである。Dingらの研究に触発され、主要な技術は陰的-陽的Runge-Kutta (IMEX-RK)設定においてShu-Osher形式を利用し、適切なモーメント更新過程を構築することで正しい極限系を捕捉することである。理論解析は、IMEX-RK時間積分と新しいモーメント更新ステップの精度条件を確立し、超強制性技術を用いて線形化モデルの安定性を確立した。数値実験は、提案格式が様々なテスト問題において精度、漸近保存性質、および堅牢性を有することを検証している。

研究背景と動機

問題背景

ボルツマン方程式の重要性: ボルツマン方程式は希薄気体動力学を記述する基礎的な動理学方程式であり、連続体モデル（例えばNavier-Stokes方程式）が失効する場合に気体流動の物理を正確に記述することができる。
計算上の課題:
- 高次元位相空間がもたらす計算複雑性
- 衝突演算子の非線形性および非局所性
- 多尺度特性に起因する硬性問題
既存手法の限界:
- 従来のオイラー法はCFL条件に制限され、時間ステップが制限される
- 小Knudsen数εの場合、陽的方法は∆t = O(ε)の厳しい時間ステップ制限を必要とする
- 半ラグランジュ法はCFL制限を回避するが、漸近保存格式の設計において新たな課題に直面する

研究動機

すべての物理尺度で精度と効率を保持する数値格式を開発すること、特に流体極限(ε→0)において正しいオイラー方程式に自動的に退化する格式を開発することが目標である。

核心的貢献

革新的なモーメント更新機構: 半ラグランジュ枠組み内で新しいモーメント更新過程を設計し、従来のオイラー法におけるモーメント更新戦略がSL枠組みで失効する問題を解決した
Shu-Osher形式の応用: Shu-Osher形式をIMEX-RK設定に適応させ、極限系を正確に識別し、対応するモーメント更新を構築することができるようにした
理論解析:
- IMEX-RK格式の精度条件を確立
- 漸近保存性質を証明
- 超強制性技術を用いて線形化モデルの安定性を分析
- 保正性の十分条件を提示
数値検証: 複数のテスト問題を通じて格式の精度、漸近保存性質、および堅牢性を検証

方法の詳細

問題定義

ボルツマン方程式を求解する： $\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に関連するMaxwell分布である。

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$