An inverse problem for semilinear wave equations on metric tree graphs

Resumen Técnico: Problema Inverso para Ecuaciones de Onda Semilineales en Grafos Métricos

1. Planteamiento del Problema

El artículo aborda un problema inverso dinámico en el contexto de ecuaciones diferenciales parciales (EDP) definidas sobre grafos métricos en forma de árbol (metric trees). El objetivo principal es recuperar la estructura geométrica del grafo y los coeficientes desconocidos de una ecuación de onda semilineal a partir de mediciones de frontera.

La Ecuación Directa:
Se considera la siguiente ecuación de onda semilineal en un grafo métrico $\Omega$ hasta un tiempo $T > 0$:
$$(\partial_t^2 - \partial_x^2)u(x, t) + q(x)u(x, t) + a(x, t)u^3(x, t) = 0$$
donde:

• $x \in \Omega \setminus V$ (puntos interiores de las aristas).
• $q(x)$ es un potencial independiente del tiempo.
• $a(x, t)$ es un coeficiente dependiente del tiempo asociado al término no lineal cúbico.
• Se imponen condiciones de Kirchhoff-Neumann en los vértices internos (continuidad de la solución y suma nula de las derivadas normales hacia afuera).
• En un subconjunto de vértices de frontera (hojas) $\Gamma_0$, se impone una condición de Dirichlet (fuente de control $f$).
• En una hoja restante $\gamma_0$, se impone una condición de Dirichlet homogénea ($u=0$).

El Problema Inverso:
Dado el operador Dirichlet-a-Neumann ($\Lambda_T$) definido en el conjunto de hojas $\Gamma_0 = \Gamma \setminus \{\gamma_0\}$, ¿es posible determinar:

1. La estructura del grafo (conectividad y longitudes de las aristas).
2. El potencial $q(x)$.
3. El coeficiente no lineal $a(x, t)$.

El problema se resuelve salvo isometría del grafo, y la recuperación de los coeficientes es válida en un dominio temporal y espacial específico que depende de $T$.

2. Metodología

Los autores emplean una combinación de técnicas de óptica geométrica, linealización y un método recursivo conocido como "leaf-peeling" (pelado de hojas) o "limpieza de hojas".

A. Linealización y Diferenciación Fraccional:
El núcleo de la estrategia es la linealización del operador de Dirichlet-a-Neumann. Dado que la ecuación es semilineal ($u^3$), los autores utilizan derivadas de Fréchet (o diferencias finitas de orden superior) con respecto a la amplitud de la fuente $f$.

• La primera linealización recupera el problema lineal ($\partial_t^2 - \partial_x^2 + q$), permitiendo recuperar la estructura del árbol y $q$ usando resultados existentes (método de leaf-peeling para ecuaciones lineales).
• La tercera derivada (o linealización de tercer orden) es crucial para recuperar el coeficiente no lineal $a(x, t)$. Se demuestra que:
  $$\partial_{\varepsilon_1}\partial_{\varepsilon_2}\partial_{\varepsilon_3} \Lambda_T (\varepsilon_1 f_1 + \varepsilon_2 f_2 + \varepsilon_3 f_3) \big|_{\varepsilon=0} = \partial_w |_{\Gamma_0}$$
  donde $w$ satisface una ecuación de onda lineal forzada por el término $-6a v_1 v_2 v_3$, siendo $v_j$ soluciones de la ecuación lineal homogénea.

B. Identidad Integral y Soluciones de Óptica Geométrica:
Mediante integración por partes, se establece una identidad integral que relaciona los datos de frontera con el coeficiente desconocido:
$$\int_0^T \int_\Omega 6 a v_0 v_1 v_2 v_3 \, dx dt = \int_0^T \int_{\Gamma_0} f_0 \partial_n w \, dt$$
donde $v_0, v_1, v_2, v_3$ son soluciones de la ecuación lineal.
Para extraer el valor puntual de $a(x,t)$, los autores construyen soluciones de óptica geométrica (paquetes de ondas de alta frecuencia) con fuentes oscilantes del tipo $e^{i h^{-1} \theta(t)}$. Estas ondas viajan a lo largo del grafo, se reflejan en los vértices y en los extremos, y sus soportes se intersectan en regiones específicas. Al hacer tender el parámetro de frecuencia $h \to 0$ y el ancho del soporte $b \to 0$, la integral se concentra en el punto de intersección de las cuatro ondas, permitiendo recuperar $a(x,t)$ puntualmente.

C. Estrategia Recursiva (Leaf-Peeling / Sheaf):
Para grafos generales, el problema se resuelve paso a paso:

1. Se identifican subgrafos tipo "sheaf" (un vértice interno conectado a varias hojas, donde todas menos una tienen control/observación).
2. Se demuestra que, en un grafo estrella, es posible recuperar $a$ en las aristas incidentes a las hojas controladas y en la arista incidente a la hoja no controlada, dentro de un dominio temporal específico.
3. Se utiliza un argumento de inducción: una vez recuperado el coeficiente en un "sheaf", este subgrafo se trata como una nueva "hoja" efectiva, permitiendo avanzar recursivamente hacia el interior del árbol hasta cubrir toda la estructura.

3. Contribuciones Clave y Resultados Principales

Teorema Principal (Teorema 1):
Sea $\Omega$ un árbol métrico finito y $\Gamma$ el conjunto de sus hojas. Si se tiene acceso al mapa Dirichlet-a-Neumann en $\Gamma_0 = \Gamma \setminus \{\gamma_0\}$ hasta un tiempo $T > 2D(\gamma_0)$ (donde $D(\gamma_0)$ es la distancia máxima desde $\gamma_0$ a cualquier otra hoja), entonces:

1. La estructura métrica del árbol (longitudes y conectividad, salvo vértices de grado 2) y el potencial $q(x)$ están únicamente determinados.
2. El coeficiente no lineal $a(x, t)$ está únicamente determinado en el dominio:
   $$\{(x, t) \in \Omega \times [0, T] : D(\gamma_0) + d(x, \gamma_0) \le t \le T - D(\gamma_0) + d(x, \gamma_0)\}$$
   Este dominio representa la región del espacio-tiempo accesible a las ondas que viajan desde las hojas controladas, se reflejan y regresan, evitando interferencias complejas que ocultarían información.

Resultados Técnicos Adicionales:

• Existencia y Unicidad: Se establece la existencia y unicidad de soluciones para el problema directo semilineal en espacios de Sobolev adecuados ($E_1$), bajo la suposición de fuentes pequeñas (Lema 6.2).
• Generalización de la No Linealidad: Aunque el papel se centra en $u^3$, los autores notan que el método funciona para no linealidades de la forma $u^m$ con $m \ge 3$.
• Robustez ante Reflexiones: A diferencia de problemas en variedades continuas, en grafos las ondas interactúan de manera compleja en los vértices. El artículo proporciona un procedimiento constructivo para aislar las interacciones en cada arista, demostrando que la información no se pierde irreversiblemente en los nodos internos si se tiene control en todas menos una hoja.

4. Significado e Impacto

Este trabajo representa un avance significativo en la teoría de problemas inversos por varias razones:

1. Primera extensión a no linealidades en grafos: Hasta la fecha, la literatura sobre problemas inversos en grafos métricos se limitaba casi exclusivamente a ecuaciones lineales. Este artículo es pionero en abordar ecuaciones semilineales, abriendo la puerta a aplicaciones en física de la materia condensada, guías de onda no lineales y fibras ópticas modeladas como redes.
2. Método Constructivo: A diferencia de enfoques puramente teóricos de unicidad, el artículo ofrece un procedimiento constructivo (basado en la linealización de alto orden y el método de pelado de hojas) para recuperar los coeficientes.
3. Superación de la "Ceguera" de los Vértices: Demuestra que, incluso con la complejidad de las reflexiones en vértices de grado mayor a 2, es posible recuperar información completa del interior del grafo si se dispone de control en todas las hojas menos una.
4. Aplicabilidad Práctica: Los resultados sugieren que es posible caracterizar redes físicas (como redes de tuberías, circuitos eléctricos no lineales o redes neuronales simplificadas) midiendo solo en la mayoría de sus extremos, lo cual es altamente relevante para aplicaciones de ingeniería y diagnóstico no destructivo.

En resumen, el artículo establece un marco teórico sólido y un algoritmo para resolver problemas inversos no lineales en estructuras de red, demostrando que la no linealidad, lejos de ser un obstáculo, puede ser explotada (a través de la linealización de orden superior) para recuperar coeficientes dependientes del tiempo que serían inaccesibles en el caso lineal.