[論文レビュー] Global uniqueness in an inverse problem for time fractional diffusion equations

Yavar Kian, Lauri Oksanen+2 more|arXiv (Cornell University)|Jan 5, 2016

Numerical methods in inverse problems参考文献 29被引用数 3

ひとこと要約

本稿は、リーマン多様体上での時間分数拡散方程式の逆問題において、全境界測定ではなく部分境界測定に基づいて、多様体構造（等長を除く）および滑らかな係数（拡散率、ポテンシャル、重み）が一意に特定可能であることを確立している。主な結果は、最小限の幾何的仮定のもとで、Dirichlet-to-Neumann写像と境界スペクトルデータを分析することに依拠している。

ABSTRACT

Given $(M,g)$, a compact connected Riemannian manifold of dimension $d \geq 2$, with boundary $\partial M$, we consider an initial boundary value problem for a fractional diffusion equation on $(0,T) imes M$, $T>0$, with time-fractional Caputo derivative of order $α\in (0,1) \cup (1,2)$. We prove uniqueness in the inverse problem of determining the smooth manifold $(M,g)$ (up to an isometry), and various time-independent smooth coefficients appearing in this equation, from measurements of the solution on a subset of $\partial M$ at fixed time. In the "flat" case where $M$ is a compact subset of $\mathbb R^d$, two out the three coefficients $ρ$ (weight), $a$ (conductivity) and $q$ (potential) appearing in the equation $ρ\partial_t^αu- extrm{div}(a abla u)+ q u=0$ on $(0,T) imes Ω$ are recovered simultaneously.

研究の動機と目的

• 時間分数拡散方程式におけるリーマン多様体の構造と時間に依存しない係数が、固定時刻における部分境界測定から一意に回復可能かどうかを特定すること。
• 固定時刻における境界の部分集合上で定義されるDirichlet-to-Neumann写像が、元の幾何的および物理的パラメータをどの程度一意に特定できるかを調査すること。
• 古典的結果（放物型・双曲型の場合）を一般化し、分数階数 α ∈ (0,1) ∪ (1,2) の場合に逆問題の同定結果を確立すること。
• 制限されたデータから、同時に拡散率、ポテンシャル、重みの複数の係数を回復することの可能性を分析すること、特に M が R^d のコンパクト部分集合である平坦な場合を対象とすること。
• DN写像から導かれる境界スペクトルデータが、多様体の等長性と係数のゲージ同値性を示す条件を探索すること。

提案手法

• コンパクトなリーマン多様体 (M,g) の境界を持つ時間分数拡散方程式（Caputo微分の階数 α ∈ (0,1) ∪ (1,2)）の初期境界値問題を定式化する。
• 境界の開部分集合 S_in に台を持つDirichletデータと、固定時刻 T_0 における別の部分集合 S_out でのノイマントレースを用いて、部分Dirichlet-to-Neumann（DN）写像を定義する。
• 適切な関数空間 H_in,α,T_0 におけるIBVPの適切性を確立し、DN写像がこの空間から L^2(S_out) への有界写像であることを保証する。
• 楕円型作用素 A = -Δ_{g,μ} + q に関連する境界スペクトルデータ（BSD）を分析し、固有値および固有関数の境界へのトレースを含む。
• 境界スペクトルデータのゲージ同値性を用いて、元の多様体が等長であること、および係数が群変換のもとで同値であることを導出する。
• スペクトル不等式と幾何的制御条件を適用し、特に双曲型逆問題における既知の結果に帰着することで、問題を簡略化する。特に、α = 2 の場合の双曲型DN写像を活用する。

実験結果

リサーチクエスチョン

• RQ1 時間分数拡散方程式に関連する部分Dirichlet-to-Neumann写像から、リーマン多様体 (M,g) が等長を除いて一意に特定可能かどうか。
• RQ2 部分境界測定から、滑らかで正の係数 μ（重み）、a（拡散率）、q（ポテンシャル）をどの程度同時に回復可能か。
• RQ3 最小限の幾何的仮定のもとで、DN写像から導かれる境界スペクトルデータが、多様体および係数を一意に特定可能かどうか。
• RQ4 多様体の幾何的性質（例：スペクトル不等式、非捕獲性）が、α ∈ (0,1) ∪ (1,2) の場合の逆問題における一意性を保証する条件は何か。
• RQ5 境界上の入力（S_in）と出力（S_out）の部分集合の相対的位置関係が、係数回復にどのように影響するか。

主な発見

• 滑らかで正の係数 μ と q を仮定した場合、部分Dirichlet-to-Neumann写像からリーマン多様体 (M,g) は等長を除いて一意に特定可能である。
• 平坦な場合（M ⊂ R^d）では、多様体がコンパクトで係数が滑らかであれば、μ、a、q はDN写像から同時に回復可能である。
• 境界スペクトルデータ（BSD）がゲージ同値である場合、対応する多様体 (M_k,g_k) は等長であり、係数 (μ_k,q_k) はゲージ変換によって関連づけられる。
• スペクトル不等式 (5.5) を満たす場合、双曲型（α = 2）のDN写像は多様体構造を一意に決定する。この結果は、還元によって分数階数の場合に拡張可能である。
• すべての n および (x,y) ∈ S_out × S_in に対して λ_{1,n} = λ_{2,n} かつ Θ_{1,n}(x,y) = Θ_{2,n}(x,y) が成り立つならば、多様体は等長であり、係数はゲージ同値である。
• S_in と S_out が互いに交わらない場合でも、スペクトル不等式 (5.5) が満たされていれば、結果は成り立つ。これは、非捕獲多様体および幾何的制御条件のもとで確認可能である。