[論文レビュー] Multilevel methods for nonuniformly elliptic operators

Long Chen, Ricardo H. Nochetto+2 more|arXiv (Cornell University)|Mar 17, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 36被引用数 6

ひとこと要約

本稿では、非一様楕円型作用素、特に分数ラプラシアンをα調和拡張を用いて分析するための多重グリッド法を開発・分析する。アスペクト比を制御する修正されたメッシュグレーディングとVサイクル多重グリッドフレームワーク内でのライン smoother を導入することにより、非等方的メッシュ上でもほぼ一様収束を達成し、数値結果では安定した反復回数と計算コストの低減が確認された。

ABSTRACT

We develop and analyze multilevel methods for nonuniformly elliptic operators whose ellipticity holds in a weighted Sobolev space with an $A_2$--Muckenhoupt weight. Using the so-called Xu-Zikatanov (XZ) identity, we derive a nearly uniform convergence result, under the assumption that the underlying mesh is quasi-uniform. We also consider the so-called $α$-harmonic extension to localize fractional powers of elliptic operators. Motivated by the scheme proposed in [R.H. Nochetto, E. Otarola, and A.J. Salgado. A PDE approach to fractional diffusion in general domains: a priori error analysis. arXiv:1302.0698, 2013] we present a multilevel method with line smoothers and obtain a nearly uniform convergence result on anisotropic meshes. Numerical experiments reveal a competitive performance of our method.

研究の動機と目的

• 分数ラプラシアン $(-\Delta)^s$ に対して $s \in (0,1)$ である離散問題のための効率的な多重グリッドソルバーを開発すること。
• 非局所的作用素を解く挑戦に直面する際、$\alpha$-調和拡張を用いて高次元の円筒領域 $\mathcal{C} = \Omega \times (0,\infty)$ に局所化すること。
• 拡張解が $y=0$ 近傍で特異的挙動を示すことを捉えるために、非等方的でグレーディングされたメッシュ上で、多重グリッド法の収束性を強固に保つこと。
• アスペクト比を制御し、事前漸近的反復増加を排除するために、メッシュグレーディングを修正することで、ライン smoother を用いたVサイクル多重グリッドソルバーの性能を向上させること。

提案手法

• 準一様メッシュ仮定の下で、近似的に一様収束を導出するためのXu–Zikatanov (XZ) 恒等式を用いる。
• Caffarelli–Silvestre の拡張を適用して、非局所的分数ラプラシアン問題を、重み付き微分作用素 $\mathrm{div}(y^\alpha \nabla \mathscr{U}) = 0$ に従う、拡張領域 $\mathcal{C}$ 上の局所的重み付きPDE に変換する。ここで $\alpha = 1 - 2s$ である。
• 拡張次元におけるアスペクト比を制御するため、パラメータ $\gamma$ を用いて導関数 $F'(\xi)$ をバランスさせることで、特に $y=1$ 近傍で有界なアスペクト比を保証する修正されたメッシュグレーディング関数 $F(\xi)$ を導入する。
• 垂直ライン smoother を $\mathcal{V}$-サイクル多重グリッドフレームワークに組み込み、修正されたメッシュによってメッシュサイズに依存しない安定的収束を実現する。
• 元のグレーディングメッシュに近い挙動を $\xi = 0$ 近傍で再現しつつ、$\xi = 1$ 近傍での挙動を改善することで非等方性を制御する、区分的線形かつ $C^1$-滑らかな写像 $F(\xi)$ を構築する。
• 1次元および2次元領域における数値実験を通じて、元のメッシュと修正メッシュを比較し、反復回数、エネルギー誤差、CPU時間の測定を行う。

実験結果

リサーチクエスチョン

• RQ1 重み付きソボレフ空間 $A_2$ 重みを伴う非一様楕円型作用素に対して、多重グリッド法がほぼ一様収束を達成できるか？
• RQ2 非等方的メッシュグレーディングは、$\alpha$-調和拡張問題に対するライン smoother を用いたVサイクル多重グリッドソルバーの性能にどのように影響するか？
• RQ3 $y=0$ 近傍の特異層の解像度を損なわずに、アスペクト比を制御する修正メッシュグレーディング関数を構築可能か？
• RQ4 修正メッシュは、特に事前漸近的領域において、メッシュの細分化に依存しない安定した反復回数をもたらすか？
• RQ5 元のメッシュグレーディング戦略と比較して、計算コストと精度のトレードオフは何か？

主な発見

• $\xi_{\star} = 0.75$ で定義される修正メッシュグレーディング $F(\xi)$ により、$\mathcal{F} = \max |F'(\xi)| \leq 4$ が達成され、アスペクト比が制御され、多重グリッド収束が安定化される。
• 元のグレーディングメッシュ (3.16) では、ライン smoother を用いたVサイクルにおいて、事前漸近的領域で反復回数が $\log J$ よりも速く増加する傾向が強く、特に細かいメッシュで顕著である。
• 修正メッシュでは、すべてのメッシュサイズで反復回数がほぼ一定（例：2次元、$s=0.15$ では10〜11回）を保ち、一様収束が示された。
• 2次元ケースで $h_{\mathscr{T}_\Omega} = 1/512$ の場合、修正メッシュは元のメッシュと比較してCPU時間を7.34秒から1.71秒に削減したが、精度は同等を維持した。
• すべてのテストケースでエネルギー誤差は0.01未満を保ち、修正メッシュはわずかに良好な条件数と反復回数の低減を達成した。
• $s=0.15$ の場合、元のメッシュグレーディングでは高い非等方性が生じるが、本手法はその場合でもロバスト性を維持する。一方、$s=0.3, 0.6, 0.8$ のような大きな $s$ の場合、元のメッシュは十分であり、修正は必要ない。