[論文レビュー] Kernel Based High Order "Explicit" Unconditionally-Stable Scheme for Nonlinear Degenerate Advection-Diffusion Equations

Andrew Christlieb, Wei Guo+1 more|arXiv (Cornell University)|Jul 28, 2017

Differential Equations and Numerical Methods参考文献 25被引用数 8

ひとこと要約

本稿では、線形でない非線形劣化拡散方程式に対して、行列の方法の転置（MOL T）におけるカーネルに基づく定式化を用いて、高次かつ無条件安定な数値スキームを提案する。連続的な畳み込み、高次WENO再構成、および明示的で強安定性を保つルンゲ＝クッタ時間積分を組み合わせることで、3次までの精度を達成し、標準的な明示スキームよりも著しく大きな時間ステップを許容する。これにより、安定性を維持しながら不連続性付近の不自然な振動を回避することができ、計算コストを大幅に削減できる。

ABSTRACT

In this paper, we present a novel numerical scheme for solving a class of nonlinear degenerate parabolic equations with non-smooth solutions. The proposed method relies on a special kernel based formulation of the solutions found in our early work on the method of lines transpose and successive convolution. In such a framework, a high order weighted essentially non-oscillatory (WENO) methodology and a nonlinear filter are further employed to avoid spurious oscillations. High order accuracy in time is realized by using the high order explicit strong-stability-preserving (SSP) Runge-Kutta method. Moreover, theoretical investigations of the kernel based formulation combined with an explicit SSP method indicates that the combined scheme is unconditionally stable and up to third order accuracy. Evaluation of the kernel based approach is done with a fast $\mathcal{O}(N)$ summation algorithm. The new method allows for much larger time step evolution compared with other explicit schemes with the same order accuracy, leading to remarkable computational savings.

研究の動機と目的

• 滑らかでない解（不連続性や鋭いフロントを伴うものなど）を有する非線形劣化放物型方程式に対して、高次数の数値スキームを開発すること。
• 特に拡散支配の問題においてΔt ∝ Δx²とされる標準的な明示スキームの深刻な時間ステップ制限を克服すること。
• 空間および時間の両方で高次精度を維持しながら、無条件安定な時間積分を達成すること。
• 高次WENO再構成と非線形フィルタを用いて、不連続性付近の振動を回避すること。
• カーネルに基づく定式化のための高速O(N)和分アルゴリズムを用いて、計算の効率化を実現すること。

提案手法

• 本手法は、行列の方法の転置（MOL T）から導かれたカーネルに基づく積分的定式化を採用し、各時間ステップで偏微分方程式を境界値問題に変換する。
• 連続的な畳み込み（解の展開）を用いて、解を高次空間微分項を含む級数として表現することで、高次精度を実現する。
• 不連続性付近での解像度を保つために、高次重み付き本質的に非振動（WENO）再構成を空間微分項に適用する。
• 急勾配領域や不連続性領域における不自然な振動を抑制するために、非線形フィルタを適用する。
• 時間積分には、安定性と精度を保つ高次強安定性保持（SSP）ルンゲ＝クッタスキームを用いる。
• カーネルに基づく積分を評価するための高速O(N)和分アルゴリズムを用い、計算コストをO(N²)からO(N)に削減する。

実験結果

リサーチクエスチョン

• RQ1 MOL Tフレームワークにおけるカーネルベース定式化は、非線形劣化拡散方程式に対して、無条件安定な時間積分を達成できるか？
• RQ2 明示的時間積分スキームを用いても、時間ステップ制限なしに空間および時間で3次までの高次精度を達成できるか？
• RQ3 WENO再構成と非線形フィルタは、不連続性付近の振動を効果的に抑制しつつ、高次精度を維持できるか？
• RQ4 提案スキームは、同等の次数の精度を持つ標準的な明示スキームよりも著しく大きな時間ステップを許容できるか？
• RQ5 高速O(N)和分アルゴリズムは、カーネルベース定式化の精度と安定性を保持しながら、高い計算効率を実現できるか？

主な発見

• 提案スキームは理論的に無条件安定であり、空間および時間で最高3次までの精度を有することが保証されている。
• 標準的な明示スキームと比較して、Δt ∝ Δx²よりもはるかに大きな時間ステップを許容でき、計算コストの大幅な削減が可能である。
• 数値実験により、設計された高次精度および不連続性や鋭いフロント付近での非振動的挙動が確認された。
• カーネルベース定式化と高速O(N)和分アルゴリズムの組み合わせにより、計算の効率が向上し、計算量の複雑性がO(N²)からO(N)に低下した。
• WENOと非線形フィルタの適用により、急勾配または不連続性を有する解において、不自然な振動が効果的に抑制された。
• 本スキームは、多孔質媒体の流れ、プラズマ輸送、放射線輸送をモデル化する方程式を含め、広範な方程式に適用可能である。これは、その頑健さと安定性のおかげである。