[論文レビュー] Half-space Kinetic Equations with General Boundary Conditions

Qin Li, Jianfeng Lu+1 more|arXiv (Cornell University)|Sep 10, 2015

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

ひとこと要約

本稿は、一般の境界条件（入射データおよびさまざまな反射型、例えば拡散・鏡面・Maxwell型）を含む線形半無限領域運動方程式に対する数値解法を開発する。減衰項の付加・除去手順を用いて、解の適切な定義（well-posedness）を確立し、スキームの準最適性（quasi-optimality）を証明した。この手法は、多成分系・多周波数・線形化BGKモデルにおいて高い精度で検証された。

ABSTRACT

We study half-space linear kinetic equations with general boundary conditions that consist of both given incoming data and various type of reflections, extending our previous work [LLS14] on half-space equations with incoming boundary conditions. As in [LLS14], the main technique is a damping adding-removing procedure. We establish the well-posedness of linear (or linearized) half-space equations with general boundary conditions and quasi-optimality of the numerical scheme. The numerical method is validated by examples including a two-species transport equation, a multi-frequency transport equation, and the linearized BGK equation in 2D velocity space.

研究の動機と目的

• 入射データと複数の反射型（拡散・鏡面・Maxwell型など）を組み合わせた一般の境界条件を有する半無限領域運動方程式の解析および数値解法を拡張すること。
• Maxwell型およびCercignani-Lampis型の反射を含む、このような一般の境界条件下での線形化半無限領域方程式の適切な定義（well-posedness）を確立すること。
• 減衰項の付加・除去手順に基づく、安定で準最適な数値スキームを構築すること。
• 物理的に重要なモデル、特に二成分輸送・多周波数輸送・2次元線形化BGK方程式に対して、この手法を検証すること。

提案手法

• 複雑な境界条件を有する半無限領域運動方程式の適切な定義と数値解法を扱うために、減衰項の付加・除去手順を導入する。
• 解の安定化と数値計算の可能性を高めるために、減衰項を追加した運動方程式のダミー版を適用する。
• 同次方程式と境界条件を満たす特別な解 $ g_0 $ を構築し、正確な解を $ f_h = f - c_h g_0 + c_h X_0 $ と表現可能にする。ここで $ f $ は減衰系の解である。
• $ ilde{ ho} $ の核空間への $ L^2 $-射影作用素 $ ilde{ ho} $ と $ ilde{ ho}^ot $ を用い、$ ilde{ ho}_1 $ を定義して核空間を正固有値・負固有値・零固有値部分空間に分解する。
• 周波数変数 $ ilde{ ho} $ における台形則による離散化を適用し、多成分系に還元することで基底関数の構成が可能になる。
• 拡散反射および鏡面反射をモデル化するため、境界作用素 $ ilde{ ho} $ と $ ilde{ ho}_d $ を取り扱い、$ (I + ar{ ilde{ ho}})^{-1} $ および $ (I + ar{ ilde{ ho}})^{-1}(I - ar{ ilde{ ho}}) $ の可逆性を導出する。

実験結果

リサーチクエスチョン

• RQ1 入射データと多様な反射機構を含む一般の境界条件を有する線形半無限領域運動方程式に対して、適切な定義（well-posedness）をどのように確立できるか？
• RQ2 核空間の分解 $ H^+ \bigoplus H^- \bigoplus H^0 $ は、無限遠における境界条件の適切な特定にどのように寄与するか？
• RQ3 減衰項の付加・除去手順を体系的に応用することで、このような方程式に対する準最適な数値スキームを構築できるか？
• RQ4 多成分系・多周波数・多次元速度構造を有するモデルに対して、この手法はどのように性能を示すか？
• RQ5 正確な解 $ f_h = X_0 $ が既知の場合、この手法はどの程度正確に解を回復できるか？

主な発見

• 数値スキームにおいて準最適性が達成され、対数因子を除いて最適な収束率が保証される。
• 二成分輸送方程式において、$ h = X_0 $ の場合、数値解は解析解と一致し、$ \tilde{\rho} $ 方向の誤差は $ 10^{-10} $ 未満であった。
• 境界条件がMaxwell型で $ (\tilde{\rho}_d, \tilde{\rho}_s) = (0.3, 0.4) $ の場合、正確な解 $ f_h = X_0 $ が回復され、差 $ f_h - X_0 $ は機械精度内に数値的にゼロであった。
• 2次元速度空間における線形化BGK方程式に対しても、本手法は異なる物理的状態においても安定に動作することが示された。
• 特別な解 $ g_0 $ の構築と係数 $ c_h = \frac{\braket{\tilde{\rho} X_0, f}}{\braket{\tilde{\rho} X_0, g_0}} $ の使用により、減衰系からの正確な解の回復が可能になった。
• 境界作用素 $ \tilde{\rho} $ が必要な強制性条件 $ (\tilde{\rho} f, f) \neq 0 $ を満たすことが示され、$ (I + \bar{\tilde{\rho}}) $ の可逆性が保証された。