[論文レビュー] p-refined RBF-FD solution of a Poisson problem

Mitja Jančič, J. Slak+1 more|arXiv (Cornell University)|Jul 8, 2021

Numerical methods in engineering参考文献 24被引用数 8

ひとこと要約

本稿では、多項式スプラインと単項式を組み合わせた空間可変な近似次数を用いたp-リファインド径数基底関数差分法（RBF-FD）を提案し、強い源項を有するポアソン問題を解く。p-リファインド手法は、グローバルに高次の方法と同等の収束速度を達成しながら、低次の方法と同程度の計算コストを維持し、強力な源項を有する領域における効率性と精度を著しく向上させる。

ABSTRACT

Local meshless methods obtain higher convergence rates when RBF approximations are augmented with monomials up to a given order. If the order of the approximation method is spatially variable, the numerical solution is said to be p-refined. In this work, we employ RBF-FD approximation method with polyharmonic splines augmented with monomials and study the numerical properties of p-refined solutions, such as convergence orders and execution time. To fully exploit the refinement advantages, the numerical performance is studied on a Poisson problem with a strong source within the domain.

研究の動機と目的

• 強い内部源項を有するポアソン問題を解く際のp-リファインドRBF-FD法の数値的性能を調査すること。
• 空間的に可変な近似次数（p-リファインド）が、計算コストを著しく増加させることなく収束速度を向上させられるかどうかを評価すること。
• p-リファインド解が、グローバルに高次のRBF-FD法と低次のRBF-FD法と比較して、精度と効率の観点からどのように異なるかを比較すること。
• p-リファインドと空間的に適応的なノード配置を組み合わせることで、解の勾配や源項の強度が大きい領域における精度をさらに向上させられるかを検討すること。
• 局所的な高次近似の利点を実証することで、将来のhp-adaptiveなメッシュレス法の基盤を築くこと。

提案手法

• 多項式スプライン（PHS）に最大次数mの単項式を付加することで収束性を保証し、近似次数を制御するRBF-FDを用いる。
• 空間的に可変な単項式付加：源に近い領域ではm = 6、中間領域ではm = 4、それ以外ではm = 2として、p-リファインドを実現する。
• 次元に依存しないノード生成アルゴリズムを用いて、準均一な散乱ノード配置を生成する。
• 得られた連立一次方程式を、1スレッドで実行するスパース直接解法（Pardiso）で解き、誤差は無限大ノルムで測定する。
• ラグランジュ乗数を用いて、RBFと単項式の両方に正確性を強制する一般化最小二乗法によりスタencil重みを導出する。
• 高次近似領域の半径（c1, c2, c3）を変化させることで、収束性とコストの制御された比較が可能な複数のp-リファインド設定をテストする。

実験結果

リサーチクエスチョン

• RQ1 p-リファインドRBF-FDが、計算コストを著しく増加させることなく、強い内部源項を有するポアソン問題の収束速度を向上させられるか？
• RQ2 高次近似（m = 6, 4, 2）の空間的分布が、全体の収束挙動と誤差低減にどのように影響するか？
• RQ3 p-リファインド解が、グローバルに高次法（m = 6）と同等の収束速度を達成しながら、低次法（m = 2）と同等の計算コストを維持できるか、その程度はどの程度か？
• RQ4 特に高次近似を用いるノードが少数（数パーセント）である場合、p-リファインドの実行時間にどのような影響があるか？
• RQ5 p-リファインドは、解の勾配や源項の強度が大きい領域において、空間的に適応的なノード配置と効果的に組み合わせられるか？

主な発見

• 中央部の小さな領域でm = 6を適用したp-リファインド解（c3）は、収束率k = −3.97を達成し、グローバルにm = 6を適用した解（k = −3.98）とほぼ同一の結果を示した。
• m = 6を適用するノードは約2%、m = 4は約5%にとどまるが、c3設定のp-リファインド解はグローバルにm = 6を適用した解と同等の収束率を達成した。
• c3設定のp-リファインド解は、グローバルにm = 6を適用した解の約2倍の速度で実行され、著しく低い計算コストで同程度の精度を達成した。
• c2設定のp-リファインド解の収束率は、m = 2（k = −3.37 vs. −1.67）よりも著しく優れており、93%のノードがm = 2を用いていながらも、局所的高次近似の有効性を確認した。
• すべてのp-リファインド解の計算時間は、m = 2のケースと同等であり、精度の向上にもかかわらず、ほとんど追加コストがないことを示した。
• 結果から、p-リファインドは、強い源項の周辺に高次スタencilを集中させることで、高次収束挙動を実現しつつ、計算コストを低く抑えることが可能であることが確認された。