A Nonlocal $p$-Laplacian Interface Model with Sharp Interface

Resumen Técnico: Un Modelo de Interfaz de p-Laplaciano No Local con Interfaz Aguda

Planteamiento del Problema
El artículo aborda el desafío de formular problemas de interfaz dentro de un marco no local. Los problemas de interfaz clásicos involucran ecuaciones diferenciales parciales (EDP) definidas en un dominio $\Omega$ dividido por una interfaz $\Gamma$, donde la solución o sus derivadas exhiben discontinuidades. Específicamente, los autores consideran la extensión del problema de interfaz de Poisson y del problema de la interfaz del $p$-Laplaciano con condiciones de membrana hacia entornos no locales. Una dificultad principal en las formulaciones no locales es que las condiciones de interfaz locales clásicas (como la continuidad de la solución y los saltos prescritos en el flujo normal) no pueden imponerse directamente. Los enfoques existentes suelen aproximar la interfaz aguda $\Gamma$ con una capa volumétrica de espesor $\delta$, lo que puede complicar el modelado de la difusión independiente en los subdominios y el manejo de los datos de contorno.

Metodología
Los autores proponen una novedosa aproximación no local basada en la energía que conserva una interfaz aguda $\Gamma$. En lugar de modificar la geometría del dominio, el modelo utiliza funcionales de energía específicos donde las condiciones de interfaz se codifican implícitamente.

1. Formulación de la Energía:

   • Para el caso lineal ($p=2$), el funcional de energía local $E(u)$ es aproximado por un funcional no local $E_\delta(u)$. El término de difusión no local se construye utilizando un núcleo $\Lambda(x,y)$ que restringe las interacciones dentro de los subdominios $\Omega_1$ y $\Omega_2$, asegurando una difusión no local independiente en cada región.
   • La condición de interfaz de Dirichlet (salto en el flujo) se gestiona mediante un término de penalización que involucra promedios no locales ($K_{\delta,i} * u$) a lo largo de la interfaz $\Gamma$.
   • La condición de contorno de Neumann se incorpora naturalmente a través de la estructura variacional.
   • El marco se extiende al caso no lineal del $p$-Laplaciano con condiciones de membrana, donde el salto de la solución a través de la interfaz es penalizado. El funcional no local resultante $F_\delta(u)$ comparte la misma forma estructural que $E_\delta(u)$.

2. Análisis Teórico:

   • Existencia y Unicidad: Los autores establecen la existencia y unicidad de los minimizadores para los funcionales no locales $E_\delta$ y $F_\delta$ utilizando el método directo del cálculo de variaciones. Esto se basa en una desigualdad de Poincaré no local recién derivada que tiene en cuenta la penalización del salto en la interfaz.
   • $\Gamma$-Convergencia: La contribución teórica central es la demostración de la $\Gamma$-convergencia. Los autores demuestran que, a medida que el horizonte no local $\delta$ desaparece ($\delta \to 0$), los funcionales no locales $\Gamma$-convergen a sus contrapartes locales ($E$ y $F$).
   • Compacidad: Combinada con los resultados de compacidad para los operadores de $p$-Laplaciano no locales, la $\Gamma$-convergencia implica que los minimizadores de los problemas no locales convergen a los minimizadores de los problemas de interfaz locales en la norma $L^p$.

3. Implementación Numérica:

   • Los autores emplean un método de elementos finitos eficiente basado en la estructura de producto tensorial de los núcleos gaussianos y bases polinómicas multicúbicas.
   • Este enfoque descompone las integrales de alta dimensión en productos de integrales unidimensionales, reduciendo significamente el costo computacional.
   • Se emplean reglas de cuadratura especiales (Gauss-Legendre) diseñadas para manejar las capas de contorno que surgen de las funciones de base cúbica cerca de las interfaces de las celdas.

Resultados Clave

• Preservación de la Interfaz Aguda: A diferencia de los enfoques previos de restricción volumétrica, el modelo propuesto mantiene una interfaz aguda, permitiendo la difusión no local independiente en subdominios disjuntos sin extender los datos de la interfaz hacia el interior (bulk).
• Prueba de Convergencia: Pruebas rigurosas confirman que los minimizadores de los funcionales no locales convergen a los minimizadores locales cuando $\delta \to 0$. Se observa que la tasa de convergencia es de primer orden.
• Validación Numérica: Experimentos numéricos en una y dos dimensiones verifican la convergencia teórica. El error $L^2$ entre la solución no local y la solución local de referencia decae linealmente con $\delta$, y los resultados son robustos ante diferentes resoluciones de malla.

Significado y Reivindicaciones
El artículo sostiene que el modelo de interfaz aguda propuesto ofrece un marco más natural y flexible para los problemas de interfaz en comparación con las aproximaciones volumétricas existentes. Al formular el problema a través de funcionales de energía, las difíciles condiciones de interfaz de flujo normal se satisfacen implícitamente. Los autores destacan que este enfoque:

1. Facilita la extensión a problemas de interfaz más complejos, como se demuestra con la aplicación exitosa al $p$-Laplaciano no lineal con condiciones de membrana.
2. Evita la necesidad de extender los datos de la interfaz (por ejemplo, $g$) de $L^2(\Gamma)$ a $H^1(\Omega)$.
3. Proporciona un vínculo matemáticamente riguroso entre los modelos de interfaz no locales y locales a través de la $\Gamma$-convergencia, asegurando que el modelo no local sea una aproximación válida del problema físico clásico a medida que el horizonte desaparece.

El trabajo no propone nuevas aplicaciones físicas, sino que proporciona una herramienta matemática y numérica robusta para modelar fenómenos de interfaz existentes en campos como la dinámica de fluidos y la ciencia de materiales donde los efectos no locales son relevantes.