[論文レビュー] Continuous Interior Penalty Finite Element Methods for the Helmholtz Equation with Large Wave Number

Haijun Wu|arXiv (Cornell University)|Jun 21, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 32被引用数 6

ひとこと要約

本稿では、虚部が正である複素数のペナルティパラメータを用いた、波数が大きい場合のヘルムホルツ方程式に対する、連続内部ペナルティ有限要素法（CIP-FEM）を提案する。この手法は、非条件安定性および最適な前漸近的誤差推定を達成し、標準的なFEMと比較して汚染誤差を顕著に低減する。誤差境界は $ C_1kh + C_2k^3h^2 $（$ k^3h^2 \leq C_0 $ のとき）および $ C_1kh + \frac{C_2}{\gamma} $（それ以外のとき）である。

ABSTRACT

This paper develops and analyzes some continuous interior penalty finite element methods (CIP-FEMs) using piecewise linear polynomials for the Helmholtz equation with the first order absorbing boundary condition in two and three dimensions. The novelty of the proposed methods is to use complex penalty parameters with positive imaginary parts. It is proved that, if the penalty parameter is a pure imaginary number $ı\ga$ with $0 C_0$ and $kh$ is bounded, where $k$ is the wave number, $h$ is the mesh size, and the $C$'s are positive constants independent of $k$, $h$, and $\ga$. Optimal order $L^2$ error estimates are also derived. The analysis is also applied if the penalty parameter is a complex number with positive imaginary part. By taking $\ga o 0+$, the above estimates are extended to the linear finite element method under the condition $k^3h^2\le C_0$. Numerical results are provided to verify the theoretical findings. It is shown that the penalty parameters may be tuned to greatly reduce the pollution errors.

研究の動機と目的

• 波数が大きい場合のヘルムホルツ問題における有限要素法における汚染誤差という、長年の課題に取り組む。
• メッシュ制約なしに安定かつ適切に定義された有限要素法を、高波数問題に適用可能に開発する。
• 2次元および3次元の線形有限要素法における前漸近的誤差推定を、複素数ペナルティパラメータの導入によって改善する。
• ペナルティパラメータの調整が汚染誤差を顕著に低減できることを示し、数値的精度を向上させる。
• ペナルティパラメータの虚部をゼロに近づける極限において、標準的なFEMの理論的結果を回復する。

提案手法

• ヘルムホルツ方程式と一次吸収境界条件を伴う、区分的線形要素を用いた連続内部ペナルティ有限要素法（CIP-FEM）を導入する。
• 虚部が正である $ \gamma > 0 $ の複素数ペナルティパラメータ $ \mathbf{i}\gamma $ を使用し、非条件安定性および適切に定義された問題を保証する。
• $ H^1 $-ノルム誤差推定を導出：$ k^3h^2 \leq C_0 $ のとき $ C_1kh + C_2k^3h^2 $、$ k^3h^2 > C_0 $ のとき $ C_1kh + \frac{C_2}{\gamma} $ であり、定数は $ k $, $ h $, $ \gamma $ に依存しない。
• 同じ条件下で最適な $ L^2 $-ノルム誤差推定を確立する。
• 一般の虚部が正である複素数ペナルティパラメータを適用し、誤差低減のためのパラメータ調整を可能にする。
• 極限 $ \gamma \to 0^+ $ において、標準的なFEMの誤差推定を回復する。これは $ k^3h^2 \leq C_0 $ の条件下で有効である。

実験結果

リサーチクエスチョン

• RQ1 虚部が正である複素数ペナルティパラメータを用いた連続内部ペナルティFEMは、波数が大きい場合のヘルムホルツ方程式に対して非条件安定性を達成できるか？
• RQ2 CIP-FEMの$ H^1 $-ノルム誤差推定は何か？ また、$ k $, $ h $, $ \gamma $ にどのように依存するか？
• RQ3 線形FEMにおける高波数ヘルムホルツ問題の汚染誤差は、複素数ペナルティパラメータを用いることで顕著に低減可能か？
• RQ4 ペナルティパラメータの選択が、実際の収束特性および相対誤差にどのように影響するか？
• RQ5 CIP-FEMは、より高価な方法（例：IPDG）と同等の精度を、自由度を少なくして達成できるか？

主な発見

• 虚部が正である複素数ペナルティパラメータ $ \mathbf{i}\gamma $ を用いたCIP-FEMは、すべての $ k > 0 $ に対して非条件安定かつ適切に定義されており、メッシュ制約なしに適用可能である。
• $ H^1 $-ノルム誤差推定は、$ k^3h^2 \leq C_0 $ のとき $ C_1kh + C_2k^3h^2 $、$ k^3h^2 > C_0 $ のとき $ C_1kh + \frac{C_2}{\gamma} $ であり、定数は $ k $, $ h $, $ \gamma $ に依存しない。
• 最適な $ L^2 $-ノルム誤差推定が導出され、手法の頑健性が確認された。
• 数値結果から、ペナルティパラメータ $ \gamma $ の調整が汚染誤差を顕著に低減することが示され、CIP-FEMは低 $ kh $ 値でも標準FEMを上回る性能を示した。
• $ k=100 $ の場合、CIP-FEMはわずか 35,971 個の自由度で $ H^1 $-半ノルムにおいて 30% の相対誤差を達成したが、標準FEMは 229,357 個の自由度を必要とした。
• CIP-FEMでは $ h(k,0.5) \geq 0.001 $ を満たす最大波数は $ k_{\text{max}} = 622 $ であり、標準FEMを上回る2倍以上であり、汚染効果の顕著な低減が示された。