[論文レビュー] Curl-free positive definite form of time-harmonic Maxwells equations well-suitable for iterative numerical solving

V.Е. Moiseenko, O. Ågren|arXiv (Cornell University)|Mar 28, 2020

Electromagnetic Simulation and Numerical Methods参考文献 10被引用数 2

ひとこと要約

本稿では、ベクトルポテンシャル（A, Φ）、磁場（H）、電束密度（D）を用いて、回転なしで正定値となる時間調和 Maxwell 方程式の形式を提案する。回転演算子をラプラシアンおよび発散の勾配項に置き換え、空間的に変化する係数を排除することで、離散系における対角優勢性と正定値性を確保し、共役勾配法による収束を効率的に実現する。この手法により、3次元問題において最適な反復回数が N^{1/3} のスケーリングを示す。

ABSTRACT

A new form of time-harmonic Maxwells equations is developed and proposed for numerical modeling. It is written for the magnetic field strength, electric displacement, vector potential and the scalar potential. There are several attractive features of this form. The first one is that the differential operator acting on these quantities is positive. The second is absence of curl operators among the leading order differential operators. The Laplacian stands for the leading order operator in the equations for the magnetic field strength, vector potential and the scalar potential, while the gradient of divergence stands for the electric displacement. The third feature is absence of space varied coefficients in the leading order differential operators that provides diagonal domination of the resulting matrix of the discretized equations. A simple example is given to demonstrate the applicability of this new form of time-harmonic Maxwells equations.

研究の動機と目的

• 標準的な定式化における回転演算子が引き起こす偽モードを回避する、時間調和 Maxwell 方程式の新たな形式の開発。
• 線形系における符号不定性を解消し、共役勾配法のような高速反復解法の利用を可能にする。
• 一次微分演算子における空間的変動係数を除去することで、離散系における対角優勢性を保証する。
• すべての場の量に同一メッシュを用いることで、有限差分離散化を簡素化し、境界条件の実装を容易にする。
• 3次元数値ベンチマークを通じて、反復解法における収束特性の向上を実証する。

提案手法

• 場の大きさを一貫した離散化に適合させるために、正規化されたベクトルポテンシャル A = a / (ik₀) とスカラーポテンシャル Φ = -iφ / (k₀) を用いる。
• 標準的な時間調和 Maxwell 方程式をポテンシャルを用いて導出し、Coulomb の規範条件 ∇·A = 0 を適用する。
• 弱形式を導出するため、積分 Q = ∫[|LA|² + |LH|² + |LD|² + |LΦ|²] dV の最小化により、二次形式を構築する。
• 得られた弱形式は、一様な構造化グリッド上での標準的な3次元有限差分法により離散化され、体積加重のノード寄与を用いる。
• 境界条件は境界上での Φ = 0 の設定と、境界およびエッジ部における H, D, A, Φ の一貫した差分ステンシルを用いて適用する。
• 線形系 Mx = b は、前処理なしの共役勾配法により解かれ、正規化された残差および誤差を用いて収束がモニタリングされる。

実験結果

リサーチクエスチョン

• RQ1 ベクトルポテンシャルとスカラーポテンシャルを用いて、回転なしで正定値な時間調和 Maxwell 方程式の形式を構築できるか？
• RQ2 一次項における回転演算子および空間的に変動する係数の欠如が、離散系における対角優勢性をもたらすか？
• RQ3 この定式化により、符号不定性を持つ系で一般的に見られる停滞を回避する共役勾配法の効率的収束が可能か？
• RQ4 この定式化下で、3次元問題におけるメッシュサイズに応じた反復回数のスケーリングはいかほどか？
• RQ5 すべての場の成分に同一メッシュを用いて、一貫性があり安定的かつ高精度な有限差分離散化を実装できるか？

主な発見

• 提案された定式化により、正定値な離散系が得られ、共役勾配法を用いた収束保証が可能になる。
• 一次微分演算子として、H, A, Φ に対してラプラシアンが、D に対して発散の勾配が用いられ、主要項から回転演算子が排除される。
• 一次項演算子における空間的変動係数の欠如により、得られる線形系行列における対角優勢性が保証される。
• 数値結果から、共役勾配法の収束が滑らかで高速であることが示され、反復回数に応じて残差が一貫して減少する。
• 離散化誤差の半分に達するまでの最適反復回数は、全メッシュノード数 N に対して概ね N^{1/3} のスケーリングを示し、効率的な反復性能が裏付けられる。
• 正弦波成分を持つ3次元電流源の解析解が、数値解と良好に一致し、精度および一貫性の妥当性が検証された。