[論文レビュー] Stability and error analysis for the Helmholtz equation with variable coefficients

Ivan G. Graham, Stefan Sauter|arXiv (Cornell University)|Mar 2, 2018

Seismic Imaging and Inversion Techniques参考文献 44被引用数 8

ひとこと要約

本稿は、一意的拡張原理とフレドホルム補題を用いて、変数係数を伴うヘルムホルツ方程式の適切性と有限要素誤差解析を確立し、解作用素に対する事前評価を導出する。安定性定数は係数の分散に指数関数的に増加し、その上限は鋭いものであり、高周波数領域における変動する波動速度の下で解析的例と数値実験によって検証されている。

ABSTRACT

We discuss the stability theory and numerical analysis of the Helmholtz equation with variable and possibly non-smooth or oscillatory coefficients. Using the unique continuation principle and the Fredholm alternative, we first give an existence-uniqueness result for this problem, which holds under rather general conditions on the coefficients and on the domain. Under additional assumptions, we derive estimates for the stability constant (i.e., the norm of the solution operator) in terms of the data (i.e. PDE coefficients and frequency), and we apply these estimates to obtain a new finite element error analysis for the Helmholtz equation which is valid at high frequency and with variable wave speed. The central role played by the stability constant in this theory leads us to investigate its behaviour with respect to coefficient variation in detail. We give, via a 1D analysis, an a priori bound with stability constant growing exponentially in the variance of the coefficients (wave speed and/or diffusion coefficient). Then, by means a family of analytic examples (supplemented by numerical experiments), we show that this estimate is sharp

研究の動機と目的

• 一意的拡張原理とフレドホルム補題を用いて、事前にエネルギー境界を要件としない、変数係数（非滑らかまたは振動的である可能性あり）を伴うヘルムホルツ方程式の解の存在と一意性を確立する。
• 特に高周波数領域において、係数の変動と周波数の関数としての安定性定数の事前評価を導出する。
• 安定性定数が収束に与える影響を考慮し、変数係数を伴うヘルムホルツ問題への有限要素誤差解析を拡張する。
• 係数の分散に伴う安定性定数の指数関数的増加推定の鋭さを調査する。
• 解析的例と数値実験（特に係数構造への摂動に対する感度を含む）により理論的境界を検証する。

提案手法

• 事前にエネルギー境界を要件としない適切性の証明に、一意的拡張原理とフレドホルム補題を用いる。
• 係数の境界と周波数の関数としての解作用素ノルム（安定性定数）の事前評価を導出する。特に非振動的波動速度に対して有効である。
• 安定性評価を適合ガラーキン法における有限要素誤差解析に適用し、修正されたメッシュ条件の下での収束を示す。
• 1次元解析を実施し、安定性定数が係数の分散に指数関数的に増加することを示す事前評価を導出する。
• 区分的定数係数を伴う解析的例を構築し、指数関数的境界の鋭さを示す。
• 数値実験を実施し、理論的推定の妥当性を検証する。特に係数位置の微小摂動に対する感度を含む。

実験結果

リサーチクエスチョン

• RQ1 標準的な事前境界が成立しない場合、変数係数を伴うヘルムホルツ方程式はどのような条件下で適切性を有するか？
• RQ2 解作用素の安定性定数は、係数関数の分散（例えば波動速度または拡散係数）にどのように依存するか？
• RQ3 高周波数挙動を示す変数係数を伴うヘルムホルツ方程式に対し、有限要素誤差解析を拡張できるか？
• RQ4 安定性定数の係数分散に伴う指数関数的依存性は鋭いか？また、数値的に検証可能か？
• RQ5 解は、特に不安定性閾値付近において、係数構造の微小摂動に対してどの程度感度を示すか？

主な発見

• 1次元解析により、安定性定数が係数の分散に指数関数的に増加することが示され、鋭い事前評価が得られた。
• 係数分散の増加に伴う安定性定数の指数関数的増加推定は、漸増する係数分散を伴う解析的例群によって鋭さが確認された。
• 数値実験により、解は係数構造の微小摂動に対して極めて感度を示すことが明らかになった。特に不安定性閾値付近で顕著であった。
• 不安定性は境界データに依存し、一部の設定では安定である一方、他の設定では解の勾配に指数関数的増加が観察された。
• 係数ジャンク位置における$10^{-6}$オーダーの摂動が、不安定な問題を安定な問題に変換する可能性があることが示され、不安定性の繊細さが浮き彫りになった。
• 変数係数を伴うヘルムホルツ問題への有限要素誤差解析が拡張され、安定性定数の指数関数的増加を含む安定性条件の下で収束が確立された。