Lipschitz stability for an inverse source problem of the wave equation with kinetic boundary conditions

Resumen Técnico: Estabilidad de Lipschitz para un Problema de Fuente Inversa de la Ecuación de Onda con Condiciones de Contorno Cinéticas

Planteamiento del Problema
Este artículo aborda un problema de fuente inversa para la ecuación de onda definida en un dominio acotado $\Omega \subset \mathbb{R}^n$ con una frontera suave $\Gamma = \Gamma_0 \cup \Gamma_1$. El sistema incorpora condiciones de contorno mixtas: una condición de Dirichlet en $\Gamma_0$ y una condición de contorno dinámica (cinética) en $\Gamma_1$. La condición dinámica tiene en cuenta el momento de la onda en la frontera y la interacción entre el volumen y la superficie, modelada por la ecuación $\partial^2_t y_\Gamma - \delta \Delta_\Gamma y_\Gamma + d\partial_\nu y + q_\Gamma(x)y_\Gamma = g(x,t)$, donde $y_\Gamma$ es la traza de la solución en la frontera, $\Delta_\Gamma$ es el operador de Laplace-Beltrami y $\partial_\nu$ denota la derivada normal.

El objetivo principal es establecer la estabilidad de Lipschitz global para la determinación simultánea de dos términos de forzamiento desconocidos, $f(x,t)$ en el volumen y $g(x,t)$ en la frontera $\Gamma_1$, a partir de una única medición parcial del flujo de frontera $\partial_\nu y$ en un subconjunto $\gamma \subset \Gamma_0$. Se asume que los términos de la fuente pertenecen a un conjunto admisible $S(C_0)$ donde sus derivadas temporales están acotadas por sus valores iniciales.

Metodología
La metodología central se basa en la derivación de una estimación de Carleman global aguda para el sistema de ondas con condiciones de contorno dinámicas. Los autores emplean un enfoque refinado que consta de los siguientes pasos:

1. Corrección de Funciones de Peso Anteriores: Los autores identifican un error crítico en la función de peso $\zeta(x,t) = |x-x_0|^2 - \beta t^2 + C_0$ utilizada en el trabajo previo de Gal y Tebou [16]. Se demuestra que la identidad $\nabla^2_\Gamma \zeta = 2I$ afirmada en [16] es incorrecta para dominios estándar (por ejemplo, una esfera). Para resolver esto, los autores introducen una nueva función de peso $\psi(x,t) = \psi_0(x) - \beta t^2 + C_1$, donde $\psi_0(x) = \mu^2(x)$ y $\mu$ es la función de Minkowski (función de medida) asociada a un conjunto convexo $\Omega_1$ contenido dentro de $\Omega$. Esta construcción asegura que $\nabla_\Gamma \psi = 0$ y $\nabla^2_\Gamma \psi = 0$ en $\Gamma_1$, satisfaciendo las condiciones de pseudo-convexidad necesarias sin depender de la identidad errónea.

2. Derivación de la Estimación de Carleman: Utilizando la nueva función de peso, los autores derivan una estimación de Carleman (Teorema 3.1) para el sistema. Un logro técnico clave es la capacidad de mantener el término que involucra el gradiente tangencial $|\nabla_\Gamma y_\Gamma|^2$ en el lado izquierdo de la desigualdad, siempre que se cumpla la condición $\delta > d$. Esto se logra seleccionando cuidadosamente el parámetro $\beta$ y utilizando las propiedades de la nueva función de peso, lo que evita la necesidad de funciones de corte (cut-off) en el tiempo que eran requeridas anteriormente.

3. Prueba de Estabilidad: El problema inverso se resuelve extendiendo un argumento reciente (evitando funciones de corte) al entorno de la frontera dinámica. Al diferenciar el sistema con respecto al tiempo y aplicar la estimación de Carleman derivada a la solución diferenciada en el tiempo, los autores relacionan las normas $L^2$ de los términos de la fuente con la observación de la derivada temporal del flujo de frontera $\partial_t \partial_\nu y$.

4. Aplicación de Controlabilidad: Como aplicación secundaria, la estimación de Carleman establecida se utiliza mediante un argumento de dualidad estándar para probar la controlabilidad exacta de frontera para el sistema con un control que actúa solo en la parte de Dirichlet $\gamma \subset \Gamma_0$.

Resultados Clave

• Teorema 3.1 (Estimación de Carleman): El artículo establece una estimación de Carleman global para la ecuación de onda con condiciones de contorno dinámicas. Esta estimación corrige y mejora el resultado en [16] mediante el uso de una función de peso válida que asegura la pseudo-convexidad en la frontera. Crucialmente, bajo la condición $\delta > d$, la estimación absorbe el término del gradiente tangencial $|\nabla_\Gamma y_\Gamma|^2$ en el lado izquierdo, proporcionando una desigualdad más aguda.
• Teorema 4.1 (Estabilidad de Lipschitz): Para un tiempo $T$ mayor que un tiempo mínimo específico $T^*$ (dependiente de la geometría y las velocidades de la onda), y asumiendo $\delta > d$, los autores demuestran una estimación de estabilidad de Lipschitz global:
  $$\|f\|_{L^2} + \|g\|_{L^2} \leq C \|\partial_t \partial_\nu y\|_{L^2(\gamma \times (0,T))}$$
  Este resultado garantiza la recuperación única y estable de los términos de la fuente a partir de una única medición de frontera.
• Teorema 5.5 (Controlabilidad Exacta): Se demuestra que el sistema es exactamente controlable en un tiempo $T > 2T^*$ utilizando una única fuerza de control que actúa en la parte de la frontera de Dirichlet $\gamma$. Esto mejora los resultados anteriores que requerían dos fuerzas de control.

Significación y Reivindicaciones
Los autores afirman que este trabajo proporciona el primer resultado de estabilidad de Lipschitz global para un problema hiperbólico inverso con condiciones de contorno dinámicas/Dirichlet generales. La significación radica en:

1. Corrección Rigurosa: El artículo corrige un error fundamental en el análisis de la función de peso del trabajo seminal de Gal y Tebou [16], el cual se había basado en una identidad de la Hessiana incorrecta para la frontera.
2. Agudeza (Sharpness): La estimación de Carleman derivada es "aguda" en el sentido de que coloca los términos de la derivada tangencial en el lado izquierdo (bajo la condición $\delta > d$), mientras que las estimaciones anteriores los colocaban en el lado derecho o requerían supuestos adicionales.
3. Avance Metodológico: El artículo demuestra que una función de peso diferente, utilizada originalmente para ecuaciones de Schrödinger, puede adaptarse con éxito a las ecuaciones de onda con condiciones de contorno dinámicas, superando las dificultades relacionadas con el acoplamiento de la dinámica de volumen y superficie.
4. Optimalidad: El resultado de controlabilidad se señala como óptimo porque logra la controlabilidad exacta con un solo control en la parte de Dirichlet, contrastando con trabajos anteriores que utilizaban dos controles.

El artículo concluye identificando problemas abiertos, específicamente el caso donde $\delta \leq d$ (donde la actual estimación de Carleman falla al absorber ciertos términos) y la extensión de estos resultados al caso donde las condiciones de contorno dinámicas se imponen en toda la frontera $\Gamma$.