[論文レビュー] Domain Decomposition Methods based on quasi-optimal transmission operators for the solution of Helmholtz transmission problems

Yassine Boubendir, Carlos Jerez-Hanckes+2 more|arXiv (Cornell University)|Oct 7, 2017

Electromagnetic Scattering and Analysis参考文献 22被引用数 3

ひとこと要約

本稿では、Dirichlet-to-Neumann (DtN) 演算子の近似から導かれた準最適な伝達演算子を用いた、Helmholtz透過問題に対する非重複領域分割法 (DDM) を提案する。複素化された超特異境界積分演算子または平方根フーリエ多重作用素を用いることで、周波数、対比度、部分領域数にほとんど依存しない、少数の反復で収束する安定した性能を達成する。特に高対比・高周波数領域でも有効である。

ABSTRACT

We present non-overlapping Domain Decomposition Methods (DDM) based on quasi-optimal transmission operators for the solution of Helmholtz transmission problems with piece-wise constant material properties. The quasi-optimal transmission boundary conditions incorporate readily available approximations of Dirichlet-to-Neumann operators. These approximations consist of either complexified hypersingular boundary integral operators for the Helmholtz equation or square root Fourier multipliers with complex wavenumbers. We show that under certain regularity assumptions on the closed interface of material discontinuity, the DDM with quasi-optimal transmission conditions are well-posed. We present a DDM framework based on Robin-to-Robin (RtR) operators that can be computed robustly via boundary integral formulations. More importantly, the use of quasi-optimal transmission operators results in DDM that converge in small numbers of iterations even in the challenging high-contrast, high-frequency regime of Helmholtz transmission problems. Furthermore, the DDM presented in this text require only minor modifications to handle the case of transmission problems in partially coated domains, while still maintaining excellent convergence properties. We also investigate the dependence of the DDM iterative performance on the number of subdomains.

研究の動機と目的

• 高対比・高周波数領域におけるHelmholtz透過問題に対する従来の領域分割法 (DDM) の収束遅延を解消すること。
• 困難なパrameter領域において多くのKrylov反復を要する古典的Robin境界条件や標準的境界積分方程式 (BIE) 形式の限界を克服すること。
• 部分的にコーティングされた領域を含む複雑な材料配置に対しても、収束性が安定したDDMフレームワークを構築すること。
• 伝達演算子における波数の複素化により、部分領域問題の適切な定式化を保証すること。
• 部分領域数や高周波数パrameterの変化に対しても反復回数を最小限に抑え、スケーラブルで効率的な解法戦略を実現すること。

提案手法

• Dirichlet-to-Neumann (DtN) 演算子の近似に基づく準最適伝達演算子を用いた非重複領域分割法 (DDM) を定式化する。
• 複素波数を用いた複素化された超特異境界積分演算子 (BIO) または平方根フーリエ多重作用素を伝達演算子として用いる。
• 伝達演算子における波数の複素化により、 coercivity 条件を満たすようにして、部分領域Helmholtz問題の適切な定式化を保証する。
• 境界積分形式による計算で得られるロビンからロビン (RtR) 表現を用いて、反復の安定性と効率性を確保する。
• 伝達演算子の修正により、完全電導体 (PEC) 界面を持つ部分コーティング領域に対しても適応可能であり、収束性能に影響を与えない。
• Krylov部分空間ソルバを界面問題に適用し、反復回数と遠方場誤差をモニタリングすることで収束を評価する。

実験結果

リサーチクエスチョン

• RQ1 DtN近似に基づく準最適伝達演算子は、Helmholtz透過問題におけるDDMの反復回数を顕著に削減できるか？
• RQ2 準最適伝達演算子を用いたDDMの収束特性は、周波数および材料対比の増加とともにどのように変化するか？
• RQ3 完全電導体 (PEC) 界面を持つ部分コーティング領域に適用した場合、DDMフレームワークは依然として安定的かつ効率的か？
• RQ4 部分領域数の増加が、提案されたDDMにおける反復回数に与える影響は何か？
• RQ5 高対比Helmholtz問題において、標準的CFIESK形式と比較して、提案手法は反復回数と精度の両面で優れているか？

主な発見

• ε₁ = 16、ω = 64 の円形領域において、準最適伝達演算子 $Z_j$ を用いたDDMは23反復で収束したが、CFIESK形式では2271反復を要した。
• 部分コーティング付き円形領域のケースでは、ω = 32 において $Z_j$ を用いたDDMは反復回数を2271（CFIESK）から23に削減し、遠方場誤差は8.5×10⁻³未満であった。
• 領域を4つの四分円部分領域に分割した場合、ω = 32 において $Z_j$ を用いたDDMは173反復を要したが、周波数増加に伴う反復回数の増加は最小限に抑えられた。
• 2つの部分領域に対して $Z_j^{PS}$ を用いたDDMはω = 32で256反復を要したが、$Z_j^a$ では529反復を要した。これは、伝達演算子の選定が収束速度に顕著な影響を与えることを示している。
• 3つの部分領域に分割されたL字型領域において、ω = 32 で $Z_j$ を用いたDDMは447反復を要した。これは、部分領域数の増加に対してもスケーラブルであることを示している。
• すべての形式において遠方場誤差は8.5×10⁻³未満で安定しており、反復回数が少ないにもかかわらず、DDMの精度はCFIESK法と同等またはそれ以上であることが示された。