[論文レビュー] Perfectly matched layer for second-order time-domain elastic wave equation: formulation and stability

Hisham Assi, R.S.C. Cobbold|arXiv (Cornell University)|Dec 13, 2013

Electromagnetic Simulation and Numerical Methods参考文献 32被引用数 3

ひとこと要約

本稿では、2次元異方性固体の弾性波シミュレーションにおける完全一致層（PML）のコンパクトな2次時間領域定式化を提示する。複素座標スケーリングと2パラメータのスケーリング関数を用いる。不安定な異方性材料に対してはスケーリングパラメータを増加させることで優れた安定性を達成し、6本の式（2本の2次偏微分方程式と4本の補助方程式）でのみ構成されるため計算コストを低減。長期的な安定シミュレーションが可能で、数値的反射を最小限に抑える。

ABSTRACT

A time domain system of equations is proposed to model elastic wave propagation in an unbounded two-dimensional anisotropic solid using perfectly matched layer (PML). Starting from a system of first-order frequency domain stress-velocity equations and using complex coordinate stretching approach with a two-parameter stretch function, a second-order formulation is obtained. The final system, which consists of just two second order equations along with four auxiliary equations, is smaller than existing formulations, thereby simplifying the problem and reducing the computational cost. The discrete stability of the solutions for a given mesh size is examined with the help of a plane-wave analysis of the corresponding continuous problem. It is shown that increasing the scaling parameter of the stretch function leads to significant stability improvements for certain anisotropic media that have known issues. Numerical computations for different isotropic and anisotropic media are used to illustrate the results.

研究の動機と目的

• 無限大領域における弾性波シミュレーションにおける不自然な反射を抑える安定なPML手法を提供すること。
• 時間領域PML定式化における方程式数を最小限に抑えることで、計算コストを低減し、実装を簡素化すること。
• 古典的PML定式化で不安定とされる異方性媒体における離散的安定性を向上させること。
• 複雑な多パラメータチューニングを必要とせず、スケーリングパラメータの調整により安定性を向上させる単純で効果的な手法を提供すること。

提案手法

• 周波数領域の応力-速度1次方程式から、2パラメータのスケーリング関数を用いた複素座標スケーリングを適用し、2次時間領域偏微分方程式系を導出する。
• 2次構造を維持し、時間領域での実装を可能にするために4本の補助方程式を導入する。
• 平面波解析を用いて離散的安定性を評価し、PML領域内の不安定モードを同定する。
• PML全域で1から高い値（例：20、90）に増加する可変スケーリングパラメータ（α̃j）を導入し、不安定な波モードを抑制する。
• COMSOL Multiphysicsを用いた有限要素法（FEM）を用いて数値的妥当性を検証し、等方性および異方性媒体における波動伝播をシミュレートする。
• PML係数に多項式変換を適用し、滑らかな変化を確保し、物理領域境界での数値的反射を低減する。

実験結果

リサーチクエスチョン

• RQ1 既存手法と比較して方程式数を削減した2次時間領域PML定式化を、2次元異方性固体における弾性波伝播に適用可能か？
• RQ2 スケーリング関数のパラメータの選定が、異方性媒体におけるPMLの離散的安定性にどのように影響するか？
• RQ3 スケーリングパラメータα̃jを増加させることで、もともと不安定な異方性媒体をPML内で安定化できるか？多パラメータの複雑な定式化を必要としないか？
• RQ4 長期間の安定性と計算効率の観点から、本定式化は古典的PMLをどの程度上回るか？
• RQ5 メッシュサイズと波の入射角が、2次定式化におけるPMLの安定性に及ぼす影響はいかほどか？

主な発見

• 提案された定式化は、2次偏微分方程式2本と補助方程式4本の合計6本の式のみを用いるため、弾性波に対する時間領域PML定式化としては報告された中で最小である。
• 材料IVでは、最適化されたスケーリングパラメータを用いることで、20 msのシミュレーション後でもPMLが安定していることが確認され、長期的安定性を示している。
• 材料IIIでは、スケーリングパラメータα̃jを1から20および90に増加させることで、特にx₁方向における不安定性が顕著に抑制されたが、x₂方向には残留不安定性が依然として存在した。
• 材料Vでは、α̃jを10および1に増加させることで、あらゆる方向で安定性が向上し、低値での観察された不安定性が解消された。
• 平面波解析により、α̃jの値が高くなると、不安定な固有値が離散的に解像可能な空間周波数範囲から外れることを確認し、安定性の向上を説明できた。
• エネルギー変化プロット（図11）では、α̃jの値が高いほど物理領域内での反射エネルギーが明確に減少しており、改善された吸収特性と安定性を裏付けた。