[論文レビュー] Continuous Interior Penalty Finite Element Method for Helmholtz Equation with High Wave Number: One Dimensional Analysis

Lingxue Zhu, Erik Burman+1 more|arXiv (Cornell University)|Nov 7, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 16被引用数 7

ひとこと要約

本稿は、高波数における1次元ヘルムホルツ方程式に対する連続内部ペナルティ（CIP）有限要素法を分析する。適切に選ばれた実数値のペナルティパラメータを用いることで、通常のFEMでは見られるポリューション誤差（高kにおける誤差）が解消され、H¹ノルムにおける最適収束が得られ、誤差は $ kh\|f\| $ で有界となる。この結果は、厳密な誤差推定と数値実験によって検証されている。

ABSTRACT

This paper addresses the properties of Continuous Interior Penalty (CIP) finite element solutions for the Helmholtz equation. The $h$-version of the CIP finite element method with piecewise linear approximation is applied to a one-dimensional model problem. We first show discrete well posedness and convergence results, using the imaginary part of the stabilization operator, for the complex Helmholtz equation. Then we consider a method with real valued penalty parameter and prove an error estimate of the discrete solution in the $H^1$-norm, as the sum of best approximation plus a pollution term that is the order of the phase difference. It is proved that the pollution can be eliminated by selecting the penalty parameter appropriately. As a result of this analysis, thorough and rigorous understanding of the error behavior throughout the range of convergence is gained. Numerical results are presented that show sharpness of the error estimates and highlight some phenomena of the discrete solution behavior.

研究の動機と目的

• 高波数における1次元ヘルムホルツ方程式に対するCIP有限要素法の安定性と収束性を分析すること。
• 安定化パラメータが振幅誤差と位相誤差、特にポリューション効果に与える影響を調査すること。
• 波数とメッシュサイズの依存性を明示的に取り入れた事前漸近的誤差推定を導出すること。
• ポリューション誤差が最適ペナルティパラメータ選択によってどのように消去可能かの条件を同定すること。
• 理論的結果を数値実験により検証し、鋭い誤差推定と離散解の挙動を示すこと。

提案手法

• 1次元ヘルムホルツモデル問題に、複素吸収境界条件を伴うhバージョンのCIP-FEM（一次元有限要素）を適用する。
• 離散的well-posednessを保証し、誤差推定を導出するためにペナルティパラメータに虚数部を導入し、誤差のオーダーを $ kh + \min(1, k^3 h^2) $ とする。
• 実数値ペナルティパラメータの場合、離散グリーン関数を構築し、誤差を最良近似誤差と位相誤差の成分に分解する。
• 位相誤差推定を導出し、$ |k_h^- - k| $ に比例することを示し、ポリューション項が離散波数と正確な波数のずれに依存することを明らかにする。
• ペナルティパラメータ $ \gamma $ が最適値 $ \gamma_o $ の $ O(h) $-近傍に選ばれるとき、ポリューション誤差が消滅し、最適誤差推定 $ \|(u - u_h)'’\lesssim kh\|f\| $ が得られることを証明する。
• 基本解系理論を用いてグローバル離散系の明示的解法を導出し、安定性および誤差解析を支援する。

実験結果

リサーチクエスチョン

• RQ1 CIP-FEMは高波数における1次元ヘルムホルツ方程式に対して、安定性と収束性の観点でどのように動作するか？
• RQ2 安定化パラメータの実部または虚部が振幅精度と位相精度に与える影響は何か？
• RQ3 CIP-FEMフレームワークにおいてポリューション誤差を消去可能か？ もし可能であれば、ペナルティパラメータにどのような条件下で可能か？
• RQ4 離散波数と正確な波数の関係は何か？ そして、これは誤差にどのように影響するか？
• RQ5 収束のための臨界な自由度数 $ N_c $ は波数 $ k $ に対してどのようにスケーリングされるか？

主な発見

• 虚数ペナルティパラメータを用いたCIP-FEMでは、H¹セミノルムにおける誤差は $ \|(u - u_h)'’\lesssim kh + \min(1, k^3 h^2)\|f\| $ で有界であり、well-posednessとロバスト性を保証する。
• ペナルティパラメータが実数値の場合、誤差は最良近似誤差とポリューション項（$ |k_h^- - k| $ に比例）に分解される。
• ペナルティパラメータ $ \gamma $ が最適値 $ \gamma_o $ の $ O(h) $-近傍に選ばれるとき、ポリューション誤差は消滅し、最適誤差推定 $ \|(u - u_h)'’\lesssim kh\|f\| $ が得られる。
• 収束のための臨界自由度数 $ N_c $ は、$ N_c \approx k/\pi $ で与えられ、最良近似と $ \gamma = \gamma_o $ のCIP-FEMの理論的予測と一致する。
• $ \gamma = -1/12 $ の場合、予測される $ N_c \sim (k^5/720)^{1/4} $ であり、数値結果は大波数において高い精度を確認している。
• 数値結果は誤差推定の鋭さを検証し、最適 $ \gamma $ の場合、相対誤差比 $ e_c/e_{ba} $ が有界に保たれることを示している。一方、非最適選択ではこの比が増大し、ポリューションの不在が確認される。