[論文レビュー] A Numerical Study on the Weak Galerkin Method for the Helmholtz Equation with Large Wave Numbers

Lin Mu, Junping Wang+2 more|arXiv (Cornell University)|Nov 2, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 40被引用数 15

ひとこと要約

本稿は、大きな波数を有するヘルムホルツ方程式を解くための弱ガレルキン（WG）有限要素法を調査し、そのロバスト性、精度、および実装の簡便さを示している。この手法は、勾配を弱く定義した不連続な多項式近似を用い、境界を越えて連続性を要求しないため、高波数においても安定した収束を達成する。相対的 $ H^1 $ 誤差は $ kh = 0.25 $ のときほとんど増加せず、汚染誤差の効果的な制御を示している。

ABSTRACT

Weak Galerkin (WG) refers to general finite element methods for partial differential equations in which differential operators are approximated by weak forms through the usual integration by parts. In particular, WG methods allow the use of discontinuous finite element functions in the algorithm design. One of such examples was recently introduced by Wang and Ye for solving second order elliptic problems. The goal of this paper is to apply the WG method of Wang and Ye to the Helmholtz equation with high wave numbers. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Our numerical experiments indicate that weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.

研究の動機と目的

• 高波数を有するヘルムホルツ方程式を解くための弱ガレルキン（WG）有限要素法の性能を評価すること。
• 凸領域および非凸領域における均質および非均質媒質において、WG法の精度、収束性、ロバスト性を評価すること。
• 特にパrameter感度および汚染誤差の観点から、既存の不連続ガレルキン法と比較してWG法の性能を評価すること。
• メッシュの細分化および多項式次数が、高波数問題における解の精度に与える影響を調査すること。

提案手法

• WG法は、解に不連続な区分的多項式関数を用い、部分積分を用いて勾配作用素の弱い形を定式化する。
• 各要素に対して離散的弱勾配が構築され、要素境界を越えて連続性を要求しない局所的微分作用素の近似が可能になる。
• 変分式は、各要素上でヘルムホルツ方程式を残差の意味で強制し、ソモフェルト放射条件を近似するロビン型吸収境界条件を用いる。
• 標準的な不連続ガレルキンスキームとは異なり、罰則パrameterを必要としないため、実装が簡素化され、パrameterチューニングが軽減される。
• 数値実験では、区分的定数および区分的線形要素を用い、波数およびメッシュサイズを変化させた収束性をテストする。
• 相対的 $ H^1 $ 誤差を計算して収束率および汚染効果を評価し、波数のロバスト性を評価するために $ kh $ を一定に保つ。

実験結果

リサーチクエスチョン

• RQ1 弱ガレルキン法は、大きな波数においてもヘルムホルツ方程式の精度と収束性を維持するか？
• RQ2 罰則を伴う不連続ガレルキン法と比較して、WG法は誤差およびパrameter依存性の観点でどのように性能を発揮するか？
• RQ3 WG法は、高波数問題における汚染効果をどの程度軽減するか？
• RQ4 メッシュの細分化および多項式次数は、振動的ヘルムホルツ解に対するWG解の精度にどのように影響するか？

主な発見

• WG法は波数 $ k = 100 $ まで安定した収束を達成し、$ kh = 0.25 $ のとき相対的 $ H^1 $ 誤差が低く保たれるため、汚染効果が最小限であることが示された。
• $ kh = 0.25 $ のとき、収束曲線はほぼ平坦でわずかな勾配を示し、波数の増加に対してもロバストであることが示された。
• $ k = 100 $ および $ h = 1/200 $ のとき、WG法は正確な解と非常に良好に一致し、メッシュを細かくするに従い視覚的な収束が向上した。
• $ k = 100 $ および $ h = 1/60 $ のとき、区分的線形要素を用いることで、WG法は解の振幅および高周波数の振動を正確に捉えた。
• 罰則パrameterを必要としないため、標準的なDG法よりも実装が簡素化され、同等の精度を達成した。
• $ k = 240 $ でさえ、$ kh $ を一定に保つ限り誤差は $ k $ とともに単調に増加するが、その勾配は小さく保たれ、汚染誤差が制御されていることが示された。