[論文レビュー] Efficient PML for the wave equation

Marcus J. Grote, Imbo Sim|arXiv (Cornell University)|Jan 2, 2010

Electromagnetic Simulation and Numerical Methods参考文献 19被引用数 52

ひとこと要約

本稿では、2次元および3次元の2階波動方程式に対する、新規で効率的な完全一致層（PML）定式化を提示する。2次元では2つの補助変数、3次元では4つの補助変数を必要とし、波動方程式を1階系に再定式化する必要がなくなる。標準的な有限差分法や有限要素法を用いた数値シミュレーションにおいて、強い安定性と長時間にわたる精度を確保し、不均質な媒質や複雑な幾何形状に対しても頑健であることが示された。

ABSTRACT

In the last decade, the perfectly matched layer (PML) approach has proved a flexible and accurate method for the simulation of waves in unbounded media. Most PML formulations, however, usually require wave equations stated in their standard second-order form to be reformulated as first-order systems, thereby introducing many additional unknowns. To circumvent this cumbersome and somewhat expensive step, we instead propose a simple PML formulation directly for the wave equation in its second-order form. Inside the absorbing layer, our formulation requires only two auxiliary variables in two space dimensions and four auxiliary variables in three space dimensions; hence it is cheap to implement. Since our formulation requires no higher derivatives, it is also easily coupled with standard finite difference or finite element methods. Strong stability is proved while numerical examples in two and three space dimensions illustrate the accuracy and long time stability of our PML formulation.

研究の動機と目的

• 1階系への再定式化に伴う計算コストの増加を回避する、2階波動方程式のための安定的で効率的なPML定式化の開発。
• 吸収層に必要な補助変数の数を最小限に抑えること——2次元では2つ、3次元では4つ——これにより実装の複雑さとメモリコストを低減。
• ラプラス変換解析を用いて、連続的設定におけるPML定式化の強い安定性と適切に定義された問題の性質を保証。
• 高階数の微分を導入しないことで、標準的な有限差分法や有限要素法との簡単な結合を可能に。
• 不均質な媒質や点源を含む2次元および3次元の設定において、長時間にわたる安定性と精度を示す数値的シミュレーションの実施。

提案手法

• 波動方程式を2階形式で表し、複素座標変換を適用することで、ラプラス変換領域におけるPML定式化を導出。
• 変換領域に補助変数を導入し、2階構造を保ちつつ変数数を最小限に抑えるように慎重に選定。
• 吸収層に垂直な方向にのみ変換を適用し、出射波の指数的減衰を実現するとともに、反射を発生させない。
• 標準的なエネルギー法と半群理論を用いて安定性を証明し、有界でないPMLを伴う連続的初期値問題の適切に定義された問題の性質を確立。
• 空間方向には2次中央差分法、時間方向にはループフロッグスキームを用いて、PML変更済み波動方程式を離散化。
• PMLは、計算領域内ではゼロ、吸収層内では正の値をとる減衰プロファイルζiを用い、最適な減衰パラメータは数値実験により調整。

実験結果

リサーチクエスチョン

• RQ1 2階波動方程式のPML定式化を、1階系への変換を回避できる形で構築することは可能か？
• RQ2 2次元および3次元において、安定性と精度を維持するために必要な補助変数の最小数は何か？
• RQ3 不均質な媒質や複雑な波フロントが存在する状況でも、提案されたPML定式化は長時間にわたって安定性を保つのか？
• RQ4 従来のPML定式化と比較して、反射誤差と長時間安定性の観点で、本手法はどのように性能を発揮するか？
• RQ5 高階数微分を導入せずに、有限要素法や他の弱形式に容易に拡張可能か？

主な発見

• 提案されたPML定式化は、2次元で2つの補助変数、3次元で4つの補助変数を必要とし、従来の2階PML定式化と比較して未知数の数を顕著に削減した。
• 2次元および3次元の数値結果から、時間t=8の間、L²誤差が7桁減少することが確認され、長時間安定性が示された。
• 異なる減衰係数（ζ̄i）に対しても本手法は頑健であり、長時間にわたるシミュレーションでも不安定性は観測されなかった。
• 2次元では点源が円形波を生成し、外部に拡散する際、わずかな不自然な反射しか発生せず、スナップショットと誤差の減少から確認された。
• 波速が変化する不均質媒質においても、斜めに入射する波を効果的に吸収し、反射や不安定性を発生させなかった。
• 点源の3次元シミュレーションでは、球面波が外部に拡散するが、t=1の時点で目に見える反射や数値的爆発は観測されず、安定性が確認された。