[論文レビュー] The probe and enclosure methods for inverse obstacle scattering problems. The past and present

Masaru Ikehata|arXiv (Cornell University)|Feb 1, 2010

Numerical methods in inverse problems参考文献 39被引用数 7

ひとこと要約

本稿では、境界測定値から未知の障害物の形状と位置を再構成することを目的とした、逆障害散乱問題におけるプローブ法とエンクロージャ法を提示する。固定周波数および時間領域データの両方の理論的基盤を確立し、指標関数と漸近解析を用いて、最小限のデータで障害物の凸包を正確に再構成し、条件を満たせば正確な幾何構造を回復できることを示している。

ABSTRACT

More than ten years ago Ikehata discovered two mathematical methods for the purpose of extracting information about the location and shape of unknown discontinuity embedded in a known background medium from observation data. The methods are called the probe and enclosure methods. This paper presents their past and recent applications to inverse obstacle scattering problems of acoustic wave.

研究の動機と目的

• 音響波動伝播における逆障害散乱問題のプローブ法とエンクロージャ法をレビューし、拡張すること。
• 特に1つの入射波または有限時間データのみを用いた場合に、限られた境界測定値から未知の障害物を再構成するという課題に取り組むこと。
• 新規の漸近解析と積分恒等式を用いて、障害物の幾何構造の再構成に理論的根拠を与えること。
• 時間依存およびマルチフィジックス逆問題への応用における未解決問題と今後の研究方向性を特定すること。

提案手法

• プローブ法は、ドリル状の関数（ドリル列）を用いて領域内部をプローブし、ヘルムホルツ方程式の解の $ L^2 $-ノルムから得られる指標関数を用いる。
• エンクロージャ法は、ディリクレ・トゥ・ノイマン写像に基づく双対性の議論と、指数関数的に増大する解の漸近解析を用い、境界積分における正の条件を導出する。
• 時間領域問題では、波動場のラプラス変換を用いて指数的重み付き解 $ w(x;\tau) $ を構成し、$ \tau \to \infty $ における漸近解析を可能にする。
• 重要な方程式には、形の $ \int_{\partial\Omega} (\partial v/\partial\nu \cdot w - \partial w/\partial\nu \cdot v) dS $ の境界積分が含まれ、その符号と対数的減衰率から幾何的情報を得る。
• この方法は、障害物領域と背景領域における解の差の表現式に依存し、エネルギー推定により勾配の振る舞いを制御する。
• 時間領域の場合、解 $ w(x;\tau) $ は波動場 $ u(x,t) $ のラプラス変換を用いて構成され、$ \tau \to \infty $ における漸近的挙動から、観測境界からの障害物までの距離が明らかになる。

実験結果

リサーチクエスチョン

• RQ1 プローブ法は、固定周波数で与えられた波数における境界測定値から、未知の障害物の形状と位置を再構成できるか？
• RQ2 エンクロージャ法は、無限個のデータ点ではなく、1つの平面波入射のみを用いても適応可能か？
• RQ3 有限時間の波動場データから、エンクロージャ法が幾何的情報を抽出できる条件は何か？
• RQ4 1つの初期源のみを用いて、透過性または音硬い障害物の凸包または正確な形状をエンクロージャ法が特定できるか？
• RQ5 時間領域におけるエンクロージャ法の指標関数の漸近的挙動は何か？ また、これと障害物から観測面までの距離とはどのように関係するか？

主な発見

• プローブ法は、$ k^2 $ が $ B_R \setminus \overline{D} $ におけるヘルムホルツ作用素のディリクレ固有値でない限り、境界データ $ \Phi_D(x,y;k) $ から障害物 $ D $ を正確に再構成できる。
• 無限個のデータを用いたエンクロージャ法は、ディリクレ・トゥ・ノイマン写像の差 $ \Lambda_0 - \Lambda_D $ を用い、ヘルムホルツ方程式の解を含む境界積分の正の性質により $ D $ を再構成する。
• 1つの入射源を用いた時間領域問題では、エンクロージャ法により、$ \lim_{\tau \to \infty} \frac{1}{2\tau} \log \left| \int_{\partial\Omega} (\partial v/\partial\nu \cdot w - \partial w/\partial\nu \cdot v) dS \right| = -\text{dist}(D, B) $ を用いて、障害物から境界までの距離が得られる。
• 媒質対比 $ h(x) $ についての条件 (A1) または (A2) を満たす場合、境界積分の符号により障害物の凸包が特定され、漸近的減衰率により境界からの距離が得られる。
• この方法は、音硬い障害物および透過性障害物の両方へ適用可能であり、後者では $ \nabla v \cdot \nabla v $ 上で媒質対比 $ h(x) $ が定符号条件を満たす必要がある。
• ソースの位置が未知であっても、ソースが観測領域の外部に compact に支持され、時間区間 $ T $ が十分に大きい限り、このアプローチはロバストである。