[论文解读] 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）映射重构 α 的 Lipschitz 稳定性。

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 映射之间的联系。
• 将经典的基于单调性的反问题框架扩展至具有各向异性电导率的无界几何形变域。

提出的方法

• 通过微分同胚 φ_θ(x) = T_θ(x₃)(x′, x₃) 将扭曲波导 Ω_θ 上的边值问题转化为直波导 Ω 上的等价问题，其中 T_θ 为旋转矩阵。
• 通过度量变换推导出各向异性电导率张量 A(x′, α(x₃)) = A₀(x′, α(x₃))，其中 A₀ 显式表示为 x′ 和 t = α(x₃) 的函数。
• 在直波导 Ω 的边界上定义 DN 映射 Λ_α: f ↦ A∇u·ν，其中 u 是满足 div(A∇u) = 0 在 Ω 内且 u|∂Ω = f 的解。
• 利用 Lax-Milgram 引理由 f ∈ Ḣ¹/²(∂Ω) 保证 H¹(Ω) 中解的存在性与唯一性，从而确保 DN 映射定义良好。
• 通过算子 I_θ 和 J_θ 建立 ∂Ω 与 ∂Ω_θ 上迹空间之间的同构关系，实现扭曲域与直波导之间边界数据的映射。
• 证明扭曲域上的 DN 映射与 Λ_α 之间通过关系式 Λ_α = J_θ* Λ̃_θ J_θ 实现酉等价，从而可在直波导上进行分析。

实验结果

研究问题

• RQ1 能否从扭曲圆柱形波导中与各向异性电导率相关的 DN 映射唯一确定未知的扭曲函数 α(x₃)？
• RQ2 所推导的电导率张量 A₀(x′, t) 是否满足标准单调性条件 ∂_t A₀(x,t) ≥ cI？若不满足，是否仍可建立稳定性？
• RQ3 当 α 为常数或接近已知常数时，该反问题是否可实现 Lipschitz 稳定性？
• RQ4 电导率张量中由几何引起的各向异性如何影响反问题的单射性与稳定性？
• RQ5 能否通过变量变换将扭曲域上的 DN 映射等价地分析为直波导上的问题？这对反问题理论有何启示？

主要发现

• 扭曲波导中的各向异性电导率反问题等价于在直波导上具有变系数各向异性电导率张量 A₀(x′, α(x₃)) 的问题。
• DN 映射 Λ_α 在 Ḣ¹/²(∂Ω) 到其对偶空间 Ḣ⁻¹/²(∂Ω) 上有界且定义良好，确保了数学一致性。
• 当 α 为常数时，本文建立了从 DN 映射重构 α 的 Lipschitz 稳定性。
• 当 DN 映射被适当近似替代时，只要 α 足够接近已知常数，稳定性结果依然成立。
• 扭曲域与直波导之间的变换保持了反问题的本质结构，使得可在直波导上使用标准分析工具。
• 当前设定下，单调性条件 ∂_t A₀ ≥ cI 不成立，因此需要新方法；作者通过几何变换与小扰动下的稳定性分析成功实现了这一目标。