[論文レビュー] A fast spectral method for the Boltzmann collision operator with general collision kernels

Irene M. Gamba, J. Haack+2 more|arXiv (Cornell University)|Oct 3, 2016

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

ひとこと要約

本稿では、一般の衝突核を持つボルツマン衝突演算子の高速スペクトル法を提示する。計算複雑性を $O(N^6)$ から $O(MN^4\log N)$ に、メモリ使用量を $O(N^6)$ から $O(MN^4)$ に削減し、硬い球体に限定された従来の高速法とは異なり、任意の衝突核の効率的シミュレーションを可能にする。この手法は、直接スペクトル法に比べて高精度かつ数個のオーダーの高速化を達成し、可変硬球および角度依存核に対応している。

ABSTRACT

We propose a simple fast spectral method for the Boltzmann collision operator with general collision kernels. In contrast to the direct spectral method \cite{PR00, GT09} which requires $O(N^6)$ memory to store precomputed weights and has $O(N^6)$ numerical complexity, the new method has complexity $O(MN^4\log N)$, where $N$ is the number of discretization points in each of the three velocity dimensions and $M$ is the total number of discretization points on the sphere and $M\ll N^2$. Furthermore, it requires no precomputation for the variable hard sphere (VHS) model and only $O(MN^4)$ memory to store precomputed functions for more general collision kernels. Although a faster spectral method is available \cite{MP06} (with complexity $O(MN^3\log N)$), it works only for hard sphere molecules, thus limiting its use for practical problems. Our new method, on the other hand, can apply to arbitrary collision kernels. A series of numerical tests is performed to illustrate the efficiency and accuracy of the proposed method.

研究の動機と目的

• ボルツマン衝突演算子の直接スペクトル法が計算複雑性とメモリ使用量の両方で $O(N^6)$ にスケーリングするという高い計算コストとメモリボトルネックを解消すること。
• 硬い球体に限定された既存の高速法とは異なり、可変硬球および角度依存モデルを含む一般の衝突核に適用可能な高速スペクトル法を開発すること。
• 前処理の要件を低減し、任意の衝突核を有するボルツマン方程式の大規模な実用的シミュレーションを可能にすること。
• マックスウェル、硬い球体、可変軟球モデルにおける数値的検証を通じて、本手法の効率性、精度、スケーラビリティを示すこと。

提案手法

• 本手法は、速度空間における畳み込みとしての衝突演算子の再定式化を基盤とし、フーリエ変換と単位球面上のレーベデフ求積法を用いて角度統合を処理するスペクトル離散化を採用する。
• 衝突演算子を速度空間における畳み込みとして再定式化し、高速フーリエ変換（FFT）を用いた高速計算を可能にする。
• 本手法は径方向と角度方向の成分を分離する：各速度次元に $N$ ポイント、球面上に $M \ll N^2$ ポイントを配置し、$M$ は球面上の求積法の精度を満たすように選定される。
• 事前計算された関数は $O(MN^4)$ のメモリに格納され、直接法の $O(N^6)$ のストレージを回避し、時間ステップ間で再利用可能である。
• 衝突核の構造を活用することで、特に可変硬球モデルにおいて重みの完全な事前計算を回避する。
• 速度差と相対角度を用いた変数変換により、衝突積分を速度空間の畳み込みとして表現し、FFTに基づく効率的評価を可能にする。

実験結果

リサーチクエスチョン

• RQ1 一般の衝突核に対して、$O(N^6)$ より低い計算複雑性を維持しながら高精度を保つことができる高速スペクトル法を設計できるか？
• RQ2 直接スペクトル法の一般核に対する $O(N^6)$ のメモリ要件を排除または著しく低減することは可能か？
• RQ3 本手法を硬い球体分子に限定されず、可変硬球および角度依存核（例：VSSモデル）にまで拡張可能か？
• RQ4 異なる衝突モデルにおいて、本手法の精度と効率は直接スペクトル法と比べてどの程度か？
• RQ5 球面求積法の解像度 $M$ が、高次モーメントの解の精度に及ぼす影響は何か？

主な発見

• 提案手法は $O(MN^4\log N)$ の計算複雑性と $O(MN^4)$ のメモリ使用量を達成し、直接スペクトル法の $O(N^6)$ の複雑性とストレージと比べて顕著な低減を実現した。
• 可変硬球（VHS）モデルでは、前処理が一切不要であり、セットアップ時間とメモリオーバーヘッドが著しく削減された。
• 数値的検証では、マックスウェル、硬い球体、可変軟球（VSS）分子のすべてのモデルで、高次モーメントにおいて三位数の精度を維持した。
• 分布関数の解における誤差は $M$ の増加に従い急速に減少し、$M=74$ の場合に直接法と同等の高精度な結果が得られた。
• 計算複雑性とメモリフットプリントの低減により、$N$ が大きくなると直接スペクトル法に比べて数個のオーダーの高速化が達成された。
• 本手法は、速度依存性と角度依存性を併せ持つ（$B \propto |v-v_*|^\gamma (1+\cos\theta)^\eta$）VSSモデルを効果的に処理でき、硬い球体仮定を越えた広範な適用可能性を示した。