[论文解读] 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数据，以及在另一子集 S_out 上测量的法向迹，定义部分Dirichlet-to-Neumann（DN）映射，时间固定为 T_0。
• 在适当的函数空间 H_in,α,T_0 中证明IBVP的适定性，确保DN映射有界地从该空间映射到 L^2(S_out)。
• 分析与DN映射相关的边界谱数据（BSD），包括椭圆算子 A = -Δ_{g,μ} + q 的特征值及特征函数的边界迹。
• 利用边界谱数据的规范等价性，推导出底层流形的等距同构关系，以及在群变换下系数的等价性。
• 应用谱不等式与几何控制条件，将问题约化为已知的双曲反问题结果，特别利用了 α = 2 时的双曲DN映射。

实验结果

研究问题

• RQ1 能否从时间分数阶扩散方程相关的部分Dirichlet-to-Neumann映射唯一确定黎曼流形 (M,g)（至等距同构）？
• RQ2 在部分边界测量下，光滑系数 μ（权重）、a（导热率）和 q（势能）在多大程度上可同时被恢复？
• RQ3 在最小几何假设下，由DN映射导出的边界谱数据是否能唯一确定流形与系数？
• RQ4 流形的几何条件（如谱不等式或非捕获性）在何种条件下能确保 α ∈ (0,1) ∪ (1,2) 情形下反问题的唯一性？
• RQ5 系数恢复过程如何依赖于边界上输入（S_in）与输出（S_out）子集的相对位置？

主要发现

• 当系数 μ 与 q 光滑且为正时，黎曼流形 (M,g) 可通过部分Dirichlet-to-Neumann映射唯一确定（至等距同构）。
• 在平坦情形（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)，结果依然成立，该条件在非捕获流形及几何控制条件下可验证。