[論文レビュー] A novel fast solver for Poisson equation with the Neumann boundary condition

Zu‐Hui Ma, Weng Cho Chew+1 more|arXiv (Cornell University)|Jul 18, 2012

Electromagnetic Scattering and Analysis被引用数 3

ひとこと要約

本稿では、2次元ポアソン方程式のノイマン境界条件に対する高速な新規スケーラブルな数値解法を提示する。この手法は、解を2段階に分離する：まず、ループ-ツリー分解に基づく高速ツリーソルバーを用いて電束密度場 D を解き、次に別の高速ソルバーでポテンシャル φ を回復する。本手法は O(N) の計算量を達成し、反復ソルバーを用いた伝統的なFEMと比較して、大幅に解法時間を短縮しており、数値例において高い効率性と精度を示している。

ABSTRACT

In this paper we present a novel fast method to solve Poisson equation in an arbitrary two dimensional region with Neumann boundary condition. The basic idea is to solve the original Poisson problem by a two-step procedure: the first one finds the electric displacement field $\mathbf{D}$ and the second one involves the solution of potential $ϕ$. The first step exploits loop-tree decomposition technique that has been applied widely in integral equations within the computational electromagnetics (CEM) community. We expand the electric displacement field in terms of a tree basis. Then, coefficients of the tree basis can be found by the fast tree solver in O(N) operations. Such obtained solution, however, fails to expand the exact field because the tree basis is not completely curl free. Despite of this, the accurate field could be retrieved by carrying out a procedure of divergence free field removal. Subsequently, potential distribution $ϕ$ can be found rapidly at the second stage with another fast approach of O(N) complexity. As a result, the method dramatically reduces solution time comparing with traditional FEM with iterative method. Numerical examples including electrostatic simulations are presented to demonstrate the efficiency of the proposed method.

研究の動機と目的

• 任意の形状の領域における2次元ポアソン方程式のノイマン境界条件を解くための、高速で頑健かつ効率的な数値解法の開発。
• 特に不規則な領域において、計算コストが高く収束が遅いという伝統的なFEMの限界を克服すること。
• 計算電磁気学におけるループ-ツリー分解と高速ツリーソルバーを活用し、静電界におけるポアソン問題に対して近似的な線形計算量 O(N) を達成すること。
• 解析的ベンチマークと複雑な源分布を含む数値例を通じて、手法の精度と効率性を検証すること。

提案手法

• 手法は解を2段階の逐次的ステップに分解する：まず電束密度場 D に対して ∇·D = ρ を解き、次にポテンシャル φ に対して ∇φ = -D/ε を解く。
• 電束密度場 D はツリーベースに展開され、計算電磁気学におけるループ-ツリー分解を活用した O(N) 計算量の高速ツリーソルバーにより係数が計算される。
• ツリーベースは完全に回転成分を含まないが、この欠陥を補うために反復的手順により発散成分を除去する。
• ポテンシャル φ は、修正済み D フィールドに再び高速ツリーソルバーを適用することで再構築され、全体の計算量は O(N) を維持する。
• 発散成分の除去ステップでは、ループ基底のグラム行列からなる対称かつ正定値な連立一次方程式を解き、CG や GMRES を用いて反復的に解く。
• ノイマン問題の要件に従い、一意性を保証するため、基準ポテンシャルをある点に固定する。

実験結果

リサーチクエスチョン

• RQ1 2次元ポアソン方程式のノイマン境界条件を解く高速かつ直接的なソルバーを、問題サイズに比例して線形にスケーリングできるように開発可能か？
• RQ2 計算電磁気学（CEM）で用いられるループ-ツリー分解技術を、不規則な領域におけるポアソン方程式の効率的解法にどのように応用できるか？
• RQ3 本手法は、ノイマン問題に対して、反復ソルバーを用いた伝統的なFEMと比較して、どの程度の計算性能向上を達成するか？
• RQ4 ツリーベースが完全に回転成分を含まない場合、その影響を受けても、手法が精度をどの程度保持できるか？

主な発見

• 複雑度比較図から、CG および GMRES を用いた伝統的なFEMと比較して、本手法は著しく短い解法時間を達成していることが示された。
• 本手法の計算複雑度は O(N) に近づき、主なコストは発散成分の除去ステップであり、これは反復ソルバーの収束に依存する。
• 解析解が存在する単純なノイマン問題において、シミュレートされたポテンシャル分布は、y=0.1 沿いの線分比較により正確な解と一致しており、高い精度を示した。
• 複雑な幾何形状を有する線源問題においても、本手法はFEMの結果と整合性のある妥当なポテンシャル分布を生成し、非長方形領域における堅牢性を裏付けた。
• 反復段階において、GMRES は CG よりも優れた性能を示し、同じ停止基準（0.01）のもとで収束に必要な反復回数が少なかった。
• ある点における基準ポテンシャルを固定することで、ノイマン問題の非一意性を効果的に処理し、一意な解を保証できた。