[論文レビュー] An inverse anisotropic conductivity problem induced bytwisting a homogeneous cylindrical domain

Mourad Choulli, Éric Soccorsi|arXiv (Cornell University)|Sep 25, 2012

Numerical methods in inverse problems参考文献 27被引用数 5

ひとこと要約

この論文は、関数 α(x₃) によってパラメトリックに表される幾何的ねじれによって生じる伝導度テンソルを有する、ねじれのある無限円筒波ガイドにおける逆異方性伝導度問題を研究している。変数変換により問題を直線的波ガイドに変換することで、α が定数または既知の定数に十分近いと仮定したもとで、Dirichlet-to-Neumann（DN）写像からの α の再構成におけるリプシッツ安定性を証明している。

ABSTRACT

We consider the inverse problem of determining the unknown function $α: \mathbb{R} ightarrow \mathbb{R}$ from the DN map associated to the operator $\mbox{div}(A(x',α(x\_3)) abla \cdot)$ acting in the infinite straight cylindrical waveguide $Ω=ω imes \mathbb{R}$, where $ω$ is a bounded domain of $\mathbb{R}^2$. Here $A=(A\_{ij}(x))$, $x=(x',x\_3) \in Ω$, is a matrix-valued metric on $Ω$ obtained by straightening a twisted waveguide. This inverse anisotropic conductivity problem remains generally open, unless the unknown function $α$ is assumed to be constant. In this case we prove Lipschitz stability in the determination of $α$ from the corresponding DN map. The same result remains valid upon substituting a suitable approximation of the DN map, provided the function $α$ is sufficiently close to some {\it a priori} fixed constant.

研究の動機と目的

• 幾何的ねじれに起因する異方性伝導度を有する無限円筒波ガイドにおける逆問題を調査すること。
• ねじれ速度を表す未知関数 α(x₃) が、境界上のDirichlet-to-Neumann（DN）写像から一意に再構成可能かどうかを特定すること。
• α が定数または既知の定数に十分近い場合に、α の再構成における安定性推定を確立すること。
• 変数変換を用いてねじれ領域におけるDN写像と、直線的波ガイドにおける変換されたDN写像との関係を分析すること。
• 古典的な単調性に基づく逆問題枠組みを、非有界かつ幾何的に変形された異方性伝導度を有する領域へと拡張すること。

提案手法

• 波ガイド Ω_θ における境界値問題を、回転行列 T_θ を用いた微分同相写像 φ_θ(x) = T_θ(x₃)(x′, x₃) を用いて直線的波ガイド Ω に変換する。
• 計量変換から、A(x′, α(x₃)) = A₀(x′, α(x₃)) として異方性伝導度テンソルを導出する。ここで A₀ は x′ と t = α(x₃) の関数として明示的に与えられる。
• u が Ω において div(A∇u) = 0 を満たし、∂Ω で u|∂Ω = f を満たすとき、境界上の DN 写像 Λ_α: f ↦ A∇u·ν を定義する。
• Lax-Milgram補題を用いて、f ∈ Ḣ¹/²(∂Ω) のもとで H¹(Ω) 内に解の存在と一意性を確立し、DN 写像が適切に定義されることを保証する。
• 境界データをねじれ領域 ∂Ω_θ と直線的領域 ∂Ω の間で写像する演算子 I_θ および J_θ を用いて、境界上の追跡空間の同型性を確立する。
• ねじれ領域における DN 写像が、関係式 Λ_α = J_θ* Λ̃_θ J_θ によりユニタリ同値であることを証明し、直線的領域における解析が可能になる。

実験結果

リサーチクエスチョン

• RQ1 ねじれのある円筒波ガイドにおける異方性伝導度に伴うDN写像から、未知のねじれ関数 α(x₃) が一意に特定可能か。
• RQ2 導出された伝導度テンソル A₀(x′, t) に対して、標準的な単調性条件 ∂_t A₀(x,t) ≥ cI が成り立つか。もし成り立たない場合、安定性は依然として確立可能か。
• RQ3 α が定数または既知の定数に十分近い場合に、逆問題におけるリプシッツ安定性が達成可能か。
• RQ4 伝導度テンソルにおける幾何的誘導異方性が、逆問題の単射性および安定性に与える影響は何か。
• RQ5 変数変換を用いてねじれ領域におけるDN写像を、直線的領域で同等に解析可能か。その結果、逆問題理論にどのような意味があるか。

主な発見

• ねじれ波ガイドにおける逆異方性伝導度問題は、変数異方性伝導度テンソル A₀(x′, α(x₃)) を有する直線的波ガイドにおける問題と同値である。
• DN 写像 Λ_α は Ḣ¹/²(∂Ω) からその双対空間 Ḣ⁻¹/²(∂Ω) への有界線形写像として適切に定義されており、数学的整合性を保証する。
• α が定数の場合、論文は DN 写像からの α の再構成におけるリプシッツ安定性を確立している。
• DN 写像が適切な近似に置き換えられても、α が既知の定数に十分近い限り、安定性結果は成立する。
• ねじれ領域と直線的領域との間の変換は、逆問題の本質的構造を保ち、標準的手法を直線的領域に適用可能にする。
• 本設定では単調性条件 ∂_t A₀ ≥ cI が成立しないが、幾何的変換と小領域摂動における安定性解析により、著者らは新たな手法を確立した。