[論文レビュー] On superconvergence of Runge-Kutta convolution quadrature for the wave equation

Jens Markus Melenk, Alexander Rieder|arXiv (Cornell University)|Apr 1, 2019

Electromagnetic Scattering and Analysis参考文献 45被引用数 6

ひとこと要約

この論文は、波動方程式に対するルンゲ＝クッタ型畳み込みクワズィ・クワズィ（RK-CQ）におけるスーパー収束現象を、収束速度を向上させる時間微分形式の導入によって説明する。時間微分を入力とすることで、ディリクレからインピーダンスへの写像に対する新たな推定を用いることで、周波数依存性を |s| の1乗分だけ低減（多角形領域では対数項を含む）し、古典的境界を超える収束性を実現する。

ABSTRACT

The semidiscretization of a sound soft scattering problem modelled by the wave equation is analyzed. The spatial treatment is done by integral equation methods. Two temporal discretizations based on Runge-Kutta convolution quadrature are compared: one relies on the incoming wave as input data and the other one is based on its temporal derivative. The convergence rate of the latter is shown to be higher than previously established in the literature. Numerical results indicate sharpness of the analysis. The proof hinges on a novel estimate on the Dirichlet-to-Impedance map for certain Helmholtz problems. Namely, the frequency dependence can be lowered by one power of $\abs{s}$(up to a logarithmic term for polygonal domains) compared to the Dirichlet-to-Neumann map.

研究の動機と目的

• 波動方程式に対するRK-CQで観測されるスーパー収束現象を説明すること。これは、収束速度が古典的予測を上回る現象である。
• 入射波とその時間微分を用いた2つの時間離散化法を分析すること。
• ディリクレからインピーダンスへの写像に対する鋭い境界を導出することで、改善された収束の理論的基盤を確立すること。
• 球面または半空間幾何における既知の境界を、凸性の仮定なしに一般な滑らかまたは多角形領域へと拡張すること。
• 高次Radau IIAスキームを用いた非凸（L字型）領域における数値実験を通じて、理論的予測の妥当性を検証すること。

提案手法

• 境界積分方程式を用いて波動方程式の散乱問題を定式化し、時間方向の半離散化に畳み込みクワズィを適用する。
• 入力データを入射波の時間微分に変更する修正形式を導入し、これにより変更された畳み込み記号 s⁻¹K(s) を得る。
• ディリクレからインピーダンスへの写像に対する新たな推定を導出し、周波数依存性を |s| の1乗分だけ低減（多角形領域では対数因子を含む）することを示す。
• 改善された境界を用いて、修正されたスキームの収束速度が標準的定式化と比較して1次上昇することを証明する。
• 半離散設定に理論を適用し、空間方向の半離散化を境界要素法で行うものと仮定した上で、時間離散化に焦点を当てる。
• 空間離散化にガラーキン境界要素法を用い、H⁻¹/²(Γ)ノルムを用いて誤差を評価する。この際、ノルム同値性のための演算子 V(1) を使用する。

実験結果

リサーチクエスチョン

• RQ1 なぜ波動方程式に対するRK-CQの収束速度が、特定の定式化において古典的予測を上回るのか？
• RQ2 境界積分作用素に対する改善された境界によって、スーパー収束現象を厳密に説明できるか？
• RQ3 非凸または多角形領域に対しても改善された収束性は成立するのか？また、幾何的形状が収束速度に与える影響は何か？
• RQ4 入射波の時間微分を入力データとして用いることで、標準的定式化と比較して収束次数にどのような影響を与えるか？
• RQ5 理論的推定は鋭いか？高次時間積分スキームに対して、数値的に妥当性を検証できるか？

主な発見

• 入射波の時間微分を用いたRK-CQスキームでは、標準的定式化と比較して収束速度が1次上昇し、記号境界において q+1−μ+1 の次数を達成する。
• ディリクレからインピーダンスへの写像において、周波数依存性はディリクレからノイマンへの写像と比較して |s| の1乗分だけ低減され、多角形領域では対数因子が付加される。
• 数値実験により、予測された収束速度が確認された：3段Radau IIA法では3次、5段法では5次であり、微分を用いた方法では対数項を除き、古典的次数に達する。
• 非凸なL字型領域に対しても、改善された収束性が観測され、理論的分析が凸幾何に限らない強靭性を示している。
• 理論的分析は鋭く、数値結果が予測された収束次数と一致しており、ディリクレからインピーダンスへの写像に対する改善された境界の妥当性が裏付けられる。
• スーパー収束効果は、s⁻¹重み付きディリクレからノイマンへの写像を単位演算子とディリクレからインピーダンスへの写像に分解することで説明され、周波数依存性がより良い境界で与えられる。