[論文レビュー] Optimized Schwarz Methods for Maxwell equations

Victorita Dolean, Martin J. Gander+1 more|ArXiv.org|Oct 17, 2006

Advanced Numerical Methods in Computational Mathematics参考文献 20被引用数 9

ひとこと要約

本稿では、マクスウェル方程式に対する最適化されたシュワーツ法の新しい階層を提案する。古典的シュワーツ法とヘルムホルツ方程式に対する最適化手法との間には驚くべき同等性が存在する。最適化された界面条件の漸近的公式を導出することで、重なりあり・なしの両方の状況で著しく高速な収束が達成された。2次元および3次元の数値実験により検証され、可変係数および非零の導電率を有する現実的なマイクロ波オーブンのシミュレーションも含まれる。

ABSTRACT

Over the last two decades, classical Schwarz methods have been extended to systems of hyperbolic partial differential equations, and it was observed that the classical Schwarz method can be convergent even without overlap in certain cases. This is in strong contrast to the behavior of classical Schwarz methods applied to elliptic problems, for which overlap is essential for convergence. Over the last decade, optimized Schwarz methods have been developed for elliptic partial differential equations. These methods use more effective transmission conditions between subdomains, and are also convergent without overlap for elliptic problems. We show here why the classical Schwarz method applied to the hyperbolic problem converges without overlap for Maxwell's equations. The reason is that the method is equivalent to a simple optimized Schwarz method for an equivalent elliptic problem. Using this link, we show how to develop more efficient Schwarz methods than the classical ones for the Maxwell's equations. We illustrate our findings with numerical results.

研究の動機と目的

• マクスウェル方程式に対する古典的手法を上回る収束速度を達成する最適化シュワーツ法の開発。
• マクスウェル方程式に対する古典的シュワーツ法と、スカラー楕円型問題に対する最適化シュワーツ法との間の隠れた同等性の特定と活用。
• 可変係数および非零導電率問題に適用可能な明示的で漸近的な最適化界面条件の公式の導出。
• 時間調和および時間離散化設定の両方における数値実験を通じた新アルゴリズムの検証、3次元マイクロ波オーブンのシミュレーションも含む。
• 楕円型問題に限らない、マクスウェル方程式のような双曲型系への最適化シュワーツ法の適用範囲の拡張。

提案手法

• 界面での界面条件の解析を可能にするために、特性変数を用いたマクスウェル方程式の再定式化。
• マクスウェル方程式に対する古典的シュワーツ法の収束特性と、低周波数近似を用いたヘルムホルツ方程式に対する最適化シュワーツ法との間の同等性を確立。
• この同等性に基づき、重なりありおよび重なりなしの両方のバリアントを含む、最適化シュワーツ法の階層を構築。
• 可変係数に適応可能で、数値コードに直接実装可能な、最適化された界面パラメータの漸近的公式を導出。
• 有限差分または有限要素離散化を用いて、時間調和および時間離散化マクスウェル方程式にこれらの手法を適用。
• 2次元および3次元の数値実験を通じて性能を検証。非定数の材料特性および鶏肉の非零導電率を有する現実的な3次元マイクロ波オーブン問題も含む。

実験結果

リサーチクエスチョン

• RQ1 なぜ古典的シュワーツ法（特性界面条件を用いる）は、楕円型問題とは異なり、重なりなしでもマクスウェル方程式に対して収束するのか？
• RQ2 マクスウェル方程式に対する古典的シュワーツ法の収束特性は、スカラー楕円型問題に対する最適化シュワーツ法と結びつけることができるか？
• RQ3 マクスウェル方程式に適用する重なりありおよび重なりなしシュワーツ法の最適界面条件は何か？ そして、それらはどのように漸近的に導出可能か？
• RQ4 新しい最適化手法は、可変係数および非零導電率を有する問題において実際の性能をどのように発揮するか？
• RQ5 システム問題とスカラー問題の同等性は、オイラー方程式のような他の双曲型系へ一般化可能か？

主な発見

• 古典的シュワーツ法は、マクスウェル方程式に対して、ヘルムホルツ方程式に対する最適化手法と根本的に同等であるため、重なりなしでも収束する。
• 古典的法の収束速度は、最適界面条件の低周波数近似を用いたヘルムホルツ方程式に対する最適化手法と一致する。
• マクスウェル方程式に対する最適化シュワーツ法の階層が構築され、古典的手法に比べて著しく改善された収束性能が得られた。
• 可変係数問題および非零導電率問題に対しても有効であることが示された、最適化界面条件の漸近的公式が導出された。
• マイクロ波オーブンの3次元シミュレーションでは、0.005 mのメッシュ、8個のサブドメインで、46〜60反復（界面条件に依存）で10⁻⁶の誤差レベルに収束した。
• 非零導電率の場合の完全な理論的解析が欠如しているにもかかわらず、複雑な幾何形状および不均質な材料特性を持つ現実応用においても、本手法は有効であることが確認された。