On the Stability of Inverse Conductivity Problem for Polyhedral Inclusions under a Single Measurement

Resumen Técnico: Estabilidad del Problema de Conductividad Inversa para Inclusiones Poliedrales bajo una Única Medición

Formulación del Problema
Este artículo aborda el problema de la conductividad inversa, un caso particular del problema clásico de Calderón, dentro de un dominio acotado Lipschitz $\Omega \subset \mathbb{R}^3$. El objetivo es determinar la forma y la ubicación de una inclusión desconocida $D$ incrustada en un medio homogéneo e isótropo. La distribución de conductividad $\sigma$ se modela como una función constante a trozos con un salto a través de la interfaz $\partial D$:
$$\sigma = 1 + (k-1)\chi_D$$
donde $k > 1$ representa el contraste de conductividad y $\chi_D$ es la función indicadora de la inclusión. Se asume que la inclusión $D$ es un poliedro convexo. El problema inverso se formula utilizando una única medición parcial en la frontera: los datos de Neumann $\partial_\nu u|_{\Gamma_0}$ en un subconjunto abierto $\Gamma_0 \subset \partial\Omega$, dada una condición de frontera de Dirichlet no trivial $f \in H^{1/2}(\partial\Omega)$. El problema directo está gobernado por la ecuación elíptica $\text{div}(\sigma \nabla u) = 0$.

Metodología y Ingredientes Técnicos
El análisis se basa en tres herramientas matemáticas principales para cerrar la brecha entre el error de medición en la frontera y el error de reconstrucción geométrica:

1. Descomposición de Singularidades: Los autores utilizan la teoría de ecuaciones elípticas en dominios no suaves para descomponer la solución $u$ cerca de la interfaz poliedral en una parte regular y partes singulares asociadas a vértices y aristas.

   • Cerca de un vértice $v$, la solución se comporta como una combinación lineal de funciones singulares de la forma $\rho^\lambda \Phi(\vartheta)$, donde $\lambda$ son exponentes singulares determinados por la geometría local y el contraste de conductividad.
   • Cerca de una arista $E$, la solución involucra términos como $r^\mu \phi(\theta)$, donde $\mu$ son exponentes singulares de arista.
   • El artículo impone suposiciones genéricas para asegurar que los coeficientes singulares principales sean no degenerados y que los términos logarítmicos (que surgen de exponentes enteros) estén ausentes.

2. Propagación de la Pequeñez: Esta técnica se utiliza para estimar la diferencia entre soluciones correspondientes a dos inclusiones diferentes ($u$ y $u'$) en el interior del dominio basándose en la discrepancia en sus datos de Cauchy en la frontera accesible $\Gamma_0$. Mediante la aplicación iterativa de la desigualdad de las tres esferas a lo largo de una cadena de bolas que conecta la frontera con el interior, los autores derivan estimaciones de estabilidad logarítmica para la diferencia de soluciones $u - u'$ y su gradiente en regiones alejadas de la inclusión.

3. Análisis Microlocal y Óptica Geométrica Compleja (CGO): Para derivar la estimación final de estabilidad, los autores construyen una identidad integral específica que involucra una solución CGO $u_0(x) = e^{\zeta \cdot (x-x_O)}$ con un vector complejo $\zeta$. Esta función de prueba se inserta en una formulación variacional sobre una región cilíndrica local cuidadosamente seleccionada adyacente a una arista del poliedro. El análisis implica estimar la identidad integral separando las contribuciones de las singularidades de vértice, las singularidades de arista y la parte regular, aprovechando las propiedades específicas de decaimiento de la solución CGO.

Resultados Clave
El resultado principal del artículo es una estimación cuantitativa de estabilidad logarítmica para la distancia de Hausdorff entre dos inclusiones poliedrales convexas $D$ y $D'$. Sea $\varepsilon$ el error de medición en la norma $H^{-1/2}(\Gamma_0)$ de la diferencia en los datos de Neumann. Bajo suposiciones geométricas y analíticas específicas (incluyendo cotas a priori sobre la geometría del poliedro, la distancia de la inclusión desde la frontera y el contraste de conductividad), los autores demuestran:
$$d_H(D, D') \leq C (\ln |\ln \varepsilon|)^{-\kappa}$$
donde $C$ y $\kappa$ son constantes positivas que dependen únicamente de los datos a priori (tales como la geometría del dominio, la longitud mínima de la arista, los ángulos de apertura de vértice/arista y el contraste de conductividad $k$).

La demostración procede mediante:

• Establecer que la distancia de Hausdorff se alcanza en un vértice.
• Utilizar la identidad integral para relacionar el error de medición en la frontera con el comportamiento local de la solución cerca de la arista.
• Equilibrar el decaimiento de la solución CGO contra el crecimiento de los términos singulares para optimizar el parámetro $\tau$ en la función de prueba.
• Combinar las estimaciones de propagación de la pequeñez con las estimaciones integrales locales para resolver la distancia geométrica $h$ (relacionada con $d_H$) en términos de $\varepsilon$.

Significado y Contribución
El artículo contribuye a la literatura sobre problemas inversos abordando la estabilidad de la determinación de inclusiones poliedrales a partir de una única medición en la frontera. Si bien los resultados de unicidad para inclusiones poliedrales en tres dimensiones han sido establecidos (por ejemplo, en [4]), los resultados de estabilidad para esta configuración específica bajo una única medición permanecían abiertos previamente.

Los autores sitúan su trabajo dentro del contexto del problema de Calderón, señalando que, si bien las conductividades suaves pueden recuperarse de manera estable con infinitas mediciones (operador de Dirichlet-a-Neumann), las conductividades constantes a trozos con saltos son más realistas para aplicaciones como la tomografía de impedancia eléctrica (EIT). Los resultados de estabilidad previos para inclusiones a menudo dependían de múltiples mediciones o suposiciones de estabilidad local (por ejemplo, perturbaciones $C^1$). Este trabajo demuestra que incluso con una única medición, se puede lograr una estimación de estabilidad logarítmica para formas poliedrales explotando sus singularidades geométricas específicas. La metodología combina el análisis de singularidades, la continuación única y técnicas microlocales para manejar la naturaleza no suave de la interfaz poliedral, proporcionando una base teórica rigurosa para la reconstrucción de dichas inclusiones.