[論文レビュー] Method of Lines Transpose: A Fast Implicit Wave Propagator

Matthew F. Causley, Andrew Christlieb+2 more|arXiv (Cornell University)|Jun 28, 2013

Electromagnetic Scattering and Analysis参考文献 5被引用数 15

ひとこと要約

本稿では、非構造格子上で波動方程式を効率的に解くために、Method of Lines Transpose (MOL T) と高速畳み込みを組み合わせた $O(N)$、A安定、2次精度の陰的波動伝搬法を提示する。Lax型補正の代わりに収束性が保証された空間積分法を採用し、ADI分割を活用することで、CFL制限を超えた大きな時間ステップを可能にし、周期的、ディリクレ、ノイマン、流出境界条件を含む多様な境界条件を扱える。複雑な幾何形状でも精度を維持する。

ABSTRACT

As a follow up to \cite{Causley2013}, we provide a detailed description of the numerical implementation of an O(N), A-stable, second order accurate solution of the wave equation, constructed from semi-discrete boundary value problems. We improve on the previous algorithm by replacing the Lax-type correction used in \cite{Causley2013}, which was necessary for convergence when $Δt < Δx/c$, with a more accurate spatial quadrature, which we prove is convergent. We also demonstrate that the resulting solver remains fast even in the case of unstructured meshes, can incorporate domain decomposition, and allows for the implementation of Dirichlet, Neumann, periodic and outflow boundary conditions. Building upon results for the 1d formulation, we utilize alternate direction implicit (ADI) splitting to achieve a fast O(N) solver in higher spatial dimensions. Our solver is built upon line objects and, combined with the flexibility of the integral solver, allows us to solve problems on arbitrary spatial domains, by embedding the boundary in a regular Cartesian mesh. Our solver is designed to couple with particle codes, where scale separation is an issue. We therefore demonstrate the ability of our solver to take time steps well beyond that of the Courant-Friedrichs-Lewy (CFL) stability limit of explicit codes.

研究の動機と目的

• 明示スキームの時間ステップ制限を克服する、高速でA安定かつ2次精度の陰的波動方程式ソルバの開発。
• Method of Lines Transpose (MOL T) を基盤とする陰的定式化により、CFL制限を超えた大きな時間ステップを実現。
• 埋め込み境界処理を用いて、任意の幾何形状におけるディリクレ、ノイマン、周期的、流出境界条件をサポート。
• 1次元のMOL T定式化をADインパルシブ（ADI）分割を用いて多次元に拡張し、$O(N)$の計算量を維持。
• スケール分離が重要なプラズマシミュレーションにおける粒子コードと結合するための点源（デルタ関数）を組み込む。

提案手法

• 時間微分を空間演算子に変換するため、Method of Lines Transpose (MOL T) を用いて半離散波動方程式を境界値問題として定式化。
• 1次元グリーン関数（減衰指数関数）の解析的性質に基づく高速畳み込みアルゴリズムを採用し、$O(N^2)$ から $O(N)$ の計算量に低減。
• 先行研究におけるLax型補正の代わりに、収束性が保証された高次精度空間積分法を採用し、安定性を損なわず精度を向上。
• 1次元BVPを各方向で逐次解くADインパルシブ（ADI）分割を用いて、1次元ソルバを多次元に拡張。
• 周期的・均一でない領域を扱う柔軟性を確保するため、フーリエ基底の代わりに多項式基底関数を採用。
• 波が領域外に反射せずに通過するように数値的に処理する流出境界条件を実装し、2次元シミュレーションで妥当性を検証。

実験結果

リサーチクエスチョン

• RQ1 MOL Tを用いた高速でA安定かつ2次精度の波動ソルバを構築し、明示法のCFL時間ステップ制限を回避できるか？
• RQ2 一般の幾何形状と境界条件に対して、境界積分法の計算コストを$O(N)$に低減しつつ精度を維持できるか？
• RQ3 ADインパルシブ（ADI）分割を用いてMOL T手法を多次元に拡張しても、安定性や収束速度が損なわれないか？
• RQ4 非構造格子または埋め込み格子上での複雑な境界条件（ディリクレ、ノイマン、周期的、流出）に対して、提案手法の性能はいかが？
• RQ5 粒子-場の自己一致（PIC）シミュレーションにおけるプラズマダイナミクスの結合に必要な点源（デルタ関数）を、本手法が正確に扱えるか？

主な発見

• すべてのテストされた幾何形状（1次元および2次元）において、ディリクレ、ノイマン、周期的、流出境界条件を含め、空間的・時間的両方で2次収束を達成。
• グリーン関数の指数的再帰関係に基づく高速畳み込みアルゴリズムにより、$O(N)$の計算量を維持。
• ディリクレ境界条件を課したダブルサークルキャビティにおける数値実験では、$N=1024$ 時に相対$L^2$誤差 $0.000105$ で2次収束を示した。
• 軸上でノイマン条件を課した四分円領域におけるベッセルモードのシミュレーションは、完全な円領域の解と高い精度で一致した。
• 周期的スリット回折格子の例では、流出境界条件に起因する可視な反射が認められず、波の領域外への透過が効果的に実現された。
• CFL = 2 の例（格子例）ではCFL制限を大幅に超える時間ステップが可能であり、時間スケールが異なる問題の効率的シミュレーションが可能となった。