Discrete FEM-BEM coupling with the Generalized Optimized Schwarz Method

技術要約：一般化最適化シュワルツ法を用いた離散FEM-BEM結合

問題設定
本論文は、有界領域および無界領域の両方における、ヘルムホルツ方程式によってモデル化される音響波伝搬問題の数値解法に取り組んでいる。目的は、不均質材料や、クロスポイント（2つ以上のサブドメインが共有する点）を含む複雑な幾何学的形状を扱うことができる、汎用的かつ効率的なソルバーを開発することである。著者らは、不均質で有界な領域に適した有限要素法（FEM）と、区分的に均質で無界な領域に最適な境界要素法（BEM）を結合させることに焦材している。対処すべき特定の課題は、古典的なFEM-BEM定式化において発生し得る「偽共鳴（spurious resonances）」（物理的にはウェルポーズドであるが、変分定式化が非自明な核を持つ波数）における不良設定性の可能性である。

手法
著者らは、一般化最適化シュワルツ法（GOSM）に基づく非重複領域分割（DD）アプローチを提案している。その手法は以下の通りである：

1. 部分構造化（Substructuring）: ドメインを非重複のサブドメインに分割する。反復ソルバーの未知数は、体積ではなく界面（スケルトン）に関連付けられる。
2. 離散フレームワーク: 本解析は、ガラーキン近似（FEMおよびBEM）を用いた完全離散設定において、[20]で示された連続理論を拡張するものである。
3. 演算子定式化: 全体問題を $(Id + \Pi S)q = \text{rhs}$ という形式の構造化方程式へと再定式化する。ここで：
   • $S$ は、最適化インピーダンス境界条件を用いて局所問題を解く散乱演算子である。
   • $\Pi$ は、非局所的な交換演算子（古典的な最適化シュワルツ法で使用される標準的なスワップ演算子の代わりとなるもの）であり、クロスポイントを処理するように設計されている。
   • $q$ は、インピーダンス・トレースを表す。
4. FEM-BEM結合: 本研究では、3つの古典的なFEM-BEM結合（Costabel、Johnson-Nédélec、Bielak-MacCamy）を具体的に分析している。また、これらの結合に関連する境界演算子 $A_\Gamma$ に関する「符号仮定（sign assumption）」（物理的な散逸に関連する）について調査している。
5. 反復ソルバー: 符号仮定が成立する場合（特にCostabel結合の場合）、演算子 $\Pi S$ は縮小写像となる。これにより、幾何級数的な収束が保証されたリチャードソン反復法の使用が可能となる。符号仮定が成立しない結合（Johnson-Nédélec、Bielak-MacCamy）については、GMResソルバーが提案されている。
6. 伝達演算子: 数値実験では、局所的なDesprés演算子、非局所的なYukawa超特異演算子、および正のDirichlet-to-Neumann写像から導出されたSchur補行列ベースの演算子を含む、様々な伝達演算子が利用されている。

主な貢献

• 離散的ウェルポーズド性: 本論文は、古典的なFEM-BEM結合の構造化定式化が、古典的な偽共鳴とは異なる波数においてウェルポーズドであることを確立している。
• 核の次元関係: 初期の変分定式化、局所問題、および構造化定式化の間の核の次元の間に、明示的な関係が確立されている（定理7.3）。
• 偽共鳴における堅牢性: 構造化された定式化（Costabel FEM-BEM結合）が、偽共鳴においてもウェルポーズドであり続けるという重要な理論的結果を示している。これは、構造化演算子の核の次元が局所問題の核の次元と一致し、特定の条件下で構造化システムが自明な核を持つためである。
• 収束理論: Costabel結合（符号仮定を満たす）について、構造化演算子が「恒等写像 + 縮小写像」の形式をとることを証明している。これにより、リチャードソン法の幾何級数的な収束が保証される。
• クロスポイントの処理: 本手法は、クロスポイントごとに追加の未知数を必要としたり、これらの点での局所問題を結合したりすることなく、クロスポイントに対して堅牢であることが示されている。これは他のDD手法における一般的な困難点である。
• h-堅牢性: 理論的な推定と数値実験は、非局所的な伝達演算子（YukawaまたはSchur補行列型）を使用する場合、メッシュサイズ $h$ に対して頑健な（h-robustな）幾何級数的な収束が達成可能であることを示している。

結果
数値実験は、均質および不均質問題、ならびに不透過な障害物や開いた空洞（クロスポイントを導入）を含む問題について、2次元で実施された。

• 反復回数: 収束に必要な反復回数が、メッシュサイズ（$h$）および波数（$\kappa$）に対してどのように変化するかを分析した。
  • 局所的な伝達演算子（Després）を用いた構成では、メッシュの細分化に伴い、反復回数が準二次的に増加した（$O(h^{1.69} \dots O(h^{1.78})$）。
  • 非局所的な伝達演算子（Yukawa、Schur補行列）を用いた構成では、メッシュの細分化に対して反復回数が準定数となり、h-堅牢性が確認された。
  • 波数に関しては、非局所的な構成は一般に準定数的または劣線形な増加を示したが、局所的な構成は劣線形な増加を示した。
• 結合の比較: 理論的な幾何級数的収束の保証は主にCostabel結合に対して導出されたが、数値結果は、テストされたシナリオにおいてJohnson-NédlecおよびBielak-MacCamy結合が同様の反復回数を与えることを示した。
• クロスポイントの堅牢性: クロスポイントを含む矩形空洞テストケースにおいて、本手法はh-堅牢性を維持し、メッシュの細分化に伴い反復回数が安定していることが示された。
• 伝達演算子の選択: 「Y-S」構成（BEMにYukawa演算子、FEMにSchur補行列を使用）が、一般に最も低い反復回数を示した。本研究は、伝達演算子のタイプをサブドメインの性質（疎な領域か密な領域か）に合わせることが、実装効率において有益であることを強調している。

意義と主張
本論文は、保証された収束境界を持つ、初のh-堅牢な反復アルゴリズムであると主張している。古典的なFEM-BEM結合を構造化形式へと再定式化することにより、著者らは、クロスポイントや偽共鳴が存在する場合でも、これらの古典的な手法をドメイン分割技術を用いて効率的に解けることを実証している。本研究は、連続的なドメイン分割理論と離散的なFEM-BEM実装の間の溝を埋めるものであり、理論的に健全かつハイパフォーマンスコンピューティングへの応用において数値的に効率的な、これらの手法を結合するための体系的な戦略を提供している。著者らは、本アプローチがクロスポイントにおける局所問題の結合という複雑さを回避しており、これが他のマルチドメインFEM-BEM定式化に対する明確な利点であると強調している。