[論文レビュー] Fast integral equation methods for linear and semilinear heat equations in moving domains

Jun Wang, Leslie Greengard+2 more|arXiv (Cornell University)|Oct 2, 2019

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

ひとこと要約

本稿では、移動域または複雑な静的領域における線形および準線形熱方程式を解くための高速で高次の空間時間適応型積分方程式法を提示する。熱ポテンシャルの高速ガウス変換ベースのスolverと、適応メッシュの細分化および高次の数値積分を組み合わせることで、行列逆行列を回避しつつ安定性と精度を達成し、剛性問題に対しても明示的時間積分が可能となる。

ABSTRACT

We present a family of integral equation-based solvers for the linear or semilinear heat equation in complicated moving (or stationary) geometries. This approach has significant advantages over more standard finite element or finite difference methods in terms of accuracy, stability and space-time adaptivity. In order to be practical, however, a number of technical capabilites are required: fast algorithms for the evaluation of heat potentials, high-order accurate quadratures for singular and weakly integrals over space-time domains, and robust automatic mesh refinement and coarsening capabilities. We describe all of these components and illustrate the performance of the approach with numerical examples in two space dimensions.

研究の動機と目的

• 複雑で時間依存する幾何形状における熱方程式を解く際の従来の有限差分/有限要素法の限界を克服すること。
• 行列逆行列を回避する積分方程式定式化により、明示スキームの安定性制約を克服すること。
• 任意の幾何形状を有する移動域における熱ポテンシャルの高次・適応型スローラーを開発すること。
• 複雑な境界を持つ2次元領域における線形および準線形熱方程式の効率的で安定的かつ高精度なシミュレーションを可能とすること。
• 本手法を拡散方程式系および変動係数を持つ拡散係数へ拡張する基盤を提供すること。

提案手法

• 放物型ポテンシャル理論を用いて、初期値、体積および層の熱ポテンシャルによる熱方程式の解を定式化する。
• 階層的クラスタリングと適応的細分化を用いた高速ガウス変換（FGT）を実装し、熱ポテンシャルの高速評価を実現する。
• 効率的な空間分割と局所的細分化を可能とするレベル制限付きクアッドツリーデータ構造を用いる。
• 空間時間境界上の特異的および弱特異的積分に対して、高次数値積分則（例：16次ラゲール多項式）を適用する。
• 適応的時間ステップと空間時間メッシュの細分化を組み合わせた予測子補正時間積分スキームを採用する。
• 高速FGTを境界積分方程式と統合し、ディリクレ・ノイマン・ロビン境界条件を効率的に計算する。

実験結果

リサーチクエスチョン

• RQ1 高速積分方程式法は、移動域または複雑な幾何形状を有する熱方程式に対して、高次精度と安定性を達成できるか？
• RQ2 大規模N粒子系において、熱ポテンシャルの評価コストをO(N)にまで低減できるか？
• RQ3 剛性反応項を有する準線形項を扱う際、本手法の性能はいかがなものか？
• RQ4 適応的空間時間メッシュ戦略が、変化する領域における精度と効率に与える影響は何か？
• RQ5 2次元領域において、変動係数または区分的定数係数拡散係数を扱えるか？

主な発見

• 16次ラゲール多項式による数値積分を用いることで、時間方向に2次精度、空間方向に高次精度を達成した。
• 73本の閉曲線を有する周期的ボックスにおいて、T = 1000Δtにおける相対L²誤差は2.7 × 10⁻⁶であり、Δt ≈ Δx²であった。
• 1時間ステップあたりのCPU時間は17.6秒であり、そのうち非階層的熱ポテンシャルに10.8秒、層ポテンシャルに6.27秒が費やされた。
• p=2のフジタモデルは周期的境界条件下で有限時間内に爆発的挙動を示し、ODE解1/(1/u_c - (t - t_c))により良好にフィットした。
• 周期的ボックスケースにおける解は、拡散による初期段階の減衰を示し、t_c ≈ 0.099以降はほぼ一定に近づく進化を示し、ODEモデルと整合的であった。
• 本手法は、複数連結領域や静的領域において、複雑な境界幾何形状および高解像度メッシュを効果的に処理する柔軟性と効率性を示した。