[論文レビュー] On deterministic approximation of the Boltzmann equation in a bounded domain

Francis Filbet|arXiv (Cornell University)|Jun 6, 2011

Gas Dynamics and Kinetic Theory参考文献 25被引用数 27

ひとこと要約

本稿では、周期的・鏡像的・拡散的境界条件を有する有界領域におけるボルツマン方程式を解く、完全に決定論的で高次の数値法を提示する。空間方向には2次精度の有限体積法を、衝突項にはフーリエスペクトル法を組み合わせることで、速度空間においてスペクトル精度を達成し、Knudsen数の広い範囲（流体力学的極限から希薄な状態まで）における効率的なシミュレーションを可能にする。

ABSTRACT

In this paper we present a fully deterministic method for the numerical solution to the Boltzmann equation of rarefied gas dynamics in a bounded domain for multi-scale problems. Periodic, specular reflection and diffusive boundary conditions are discussed and investigated numerically. The collision operator is treated by a Fourier approximation of the collision integral, which guarantees spectral accuracy in velocity with a computational cost of $N\\,\\log(N)$, where $N$ is the number of degree of freedom in velocity space. This algorithm is coupled with a second order finite volume scheme in space and a time discretization allowing to deal for rarefied regimes as well as their hydrodynamic limit. Finally, several numerical tests illustrate the efficiency and accuracy of the method for unsteady flows (Poiseuille flows, ghost effects, trend to equilibrium).

研究の動機と目的

• 有界領域における時空間依存ボルツマン方程式に対して、マルチスケール挙動を示す完全に決定論的な数値法の開発。
• 周期的境界条件、鏡像反射境界条件、および拡散的（マックスウェル型）境界条件を含む、さまざまな境界条件を正確に扱う。
• 2次精度の有限体積法とスペクトル法を組み合わせることで、空間および速度空間において高次精度を達成する。
• 時間分割法を分離することで、希薄状態から流体力学的極限まで、Knudsen数の全範囲におけるシミュレーションを可能にする。
• 連続体モデルが失敗する、特にゴースト効果および平衡への傾向という既知の漸近的極限に対して、手法の妥当性を検証する。

提案手法

• ボルツマン方程式の輸送項の空間離散化に、2次精度の有限体積法を用いる。
• 衝突積分にフーリエスペクトル法を適用し、$N\log N$ の計算コストで速度空間においてスペクトル精度を確保する。
• 作用素分割を用いて輸送項と衝突項を分離し、時間積分の効率化と並列処理を可能にする。
• 時間離散化を設計し、流体力学的極限を含む複数スケールにわたる安定性と精度を維持する。
• 境界条件はハイブリッドマックスウェル型モデルにより処理：割合 $\alpha$ がマックスウェル分布に再発射され、残りの $1-\alpha$ が鏡像反射される。
• 境界での質量保存則を用いて壁の緩和係数 $\mu(t,x)$ を求め、物理的整合性を保証する。

実験結果

リサーチクエスチョン

• RQ1 完全に決定論的な手法が、有界領域におけるボルツマン方程式に対して、空間および速度空間の両方で高次精度を達成できるか？
• RQ2 小さなKnudsen数を有する周期的温度勾配駆動流れにおいて、ゴースト効果をどれほど正確に捉えることができるか？
• RQ3 流体力学的極限において、数値解が古典的熱伝導方程式ではなく、漸近理論に収束するか？
• RQ4 1つのアルゴリズムで、希薄状態と流体力学的状態の両方を、最小限のグリッド解像度で効率的にシミュレートできるか？
• RQ5 スペクトル法による速度空間の衝突項処理は、従来の数値積分ベースの衝突ソルバーと比較して、精度および計算コストの面で優れているか？

主な発見

• 空間方向で2次精度、速度空間でスペクトル精度を達成しており、比較的少ない速度グリッドポイント数で高解像度の結果を得られる。
• ポアズイユ流れの数値シミュレーションは、理論的予測およびベンチマーク解と良好に一致する。
• ゴースト効果問題において、Knudsen数 $\varepsilon$ が小さくなるにつれて温度場が上昇し、Soneらの漸近理論に収束するが、熱伝導方程式の解には収束しない。
• 極限 $\varepsilon \to 0$ において平均速度場が消えることは、漸近理論と整合しており、誤った流れの発生がないことを確認する。
• 周期的壁温度問題において、非平衡温度プロファイルを正しく捉えており、熱方程式が予測する等温線とは異なる形状を示す。
• 衝突ステップの計算コストは $N\log N$ スケーリングを示し、高解像度の速度グリッドに対しても効率的であり、分離構造のおかげで並列処理が可能となる。