Evolving finite elements for advection diffusion with an evolving interface

Resumen Técnico: Elementos Finitos Evolutivos para Advección-Difusión con una Interfaz Evolutiva

Planteamiento del Problema
Este artículo aborda la aproximación numérica de ecuaciones parabólicas definidas en un dominio estacionario $\Omega$ que contiene una interfaz móvil $\Gamma(t)$. La interfaz divide el dominio en dos subdominios, $\Omega_1(t)$ y $\Omega_2(t)$. Las ecuaciones gobernantes modelan el transporte y la difusión de campos escalares $u_1$ y $u_2$ dentro de estos subdominios, sujetos a:

• Dinámicas de advección-difusión-reacción dentro de cada subdominio con coeficientes dependientes del tiempo ($A_i, B_i, C_i$).
• Una condición de frontera de Dirichlet en la frontera exterior $\partial\Omega$.
• Continuidad de la solución a través de la interfaz ($u_1 = u_2$).
• Un salto de flujo prescrito a través de la interfaz ($A_1 \partial_{\nu} u_1 - A_2 \partial_{\nu} u_2 = g$).

Los desafíos principales identificados son la evolución de los subdominios a lo largo del tiempo y la presencia del salto de flujo a través de la interfaz móvil. Si bien existen aplicaciones en dinámica de fluidos, ciencia de materiales y biología celular, el artículo se centra en la formulación matemática y numérica del problema.

Metodología
Los autores proponen un enfoque Lagrangiano-Euleriano Arbitrario (ALE) que utiliza Espacios de Elementos Finitos Evolutivos (ESFEM) basados en el marco abstracto establecido en [21].

1. Marco Analítico Funcional: El problema continuo se formula dentro del contexto de espacios de Hilbert evolutivos. Los autores definen pares compatibles de espacios $(X(t), \phi_t)$ donde $\phi_t$ es un mapa de flujo que transporta funciones desde un dominio de referencia hacia el dominio evolutivo. Esto permite la definición de derivadas materiales débiles y la derivación de una formulación débil bien planteada (Teorema 2.16).
2. Construcción de la Malla:
   • Malla Inicial: El dominio se triángula inicialmente utilizando símplices lineales por partes. Crucialmente, los vértices de las facetas de la interfaz están restringidos a caer exactamente sobre la interfaz inicial $\Gamma(0)$ (Suposición M5).
   • Mapeo Isoparamétrico: Para lograr una precisión de orden superior, la malla lineal se mapea a una malla isoparamétrica utilizando un difeomorfismo $\Psi_h$. Este mapa empuja las facetas de la interfaz hacia la interfaz curva exacta $\Gamma(0)$ preservando la continuidad de la malla.
   • Evolución Temporal: La malla evoluciona moviendo los nodos lagrangianos de acuerdo con un campo de velocidad global conocido y suave $w$. El flujo discreto $\Phi^h$ se construye como una interpolación polinómica del flujo continuo, asegurando que la malla permanezca como una subdivisión admisible conforme en todo momento.
3. Mapeo de Elevación: Se introduce un operador de elevación $\Lambda_h$ para mapear funciones desde el dominio computacional discreto (que aproxima la interfaz curva) de vuelta al dominio continuo exacto. Esto es esencial para comparar la solución discreta con la solución exacta en el análisis de error.
4. Esquema Discreto: El método emplea elementos finitos isoparamétricos de orden $k$. El problema variacional discreto se formula utilizando la derivada material discreta y los coeficientes elevados. El esquema se resuelve mediante un método de pasos temporales (específicamente Fórmulas de Diferencia hacia Atrás, BDF, en los experimentos numéricos).

Contribuciones Clave

• Buen Planteamiento: El artículo proporciona un marco analítico funcional riguroso que prueba la existencia y unicidad de la solución del problema continuo (Teorema 2.16), incluyendo resultados de regularidad adicionales bajo suposiciones específicas sobre los datos.
• Marco ESFEM Isoparamétrico: Los autores extienden el marco abstracto de elementos finitos evolutivos para manejar problemas de interfaz de dos lados con saltos de flujo. Detallan la construcción de elementos isoparamétricos que aproximan con precisión la frontera curva, permitiendo una convergencia de orden superior.
• Cotas de Error Óptimas: El resultado teórico principal (Teorema 5.1) establece cotas de error de orden óptimo. Específicamente, el error en la norma $L^2$ está acotado por $O(h^{k+1})$, y el error en la semi-norma $H^1$ está acotado por $O(h^k)$, donde $h$ es el tamaño de la malla y $k$ es el grado polinómico. Estas cotas se mantienen hasta una constante y coinciden con el error de interpolación de una solución suave conocida.
• Implementación y Verificación: El artículo incluye una implementación robusta utilizando el paquete Firedrake y GMSH. Los experimentos numéricos en 2D y 3D verifican las tasas de convergencia teóricas para varios órdenes polinómicos.

Resultados
Los experimentos numéricos demuestran la eficacia del método en un problema modelo con una solución exacta conocida que es continua pero no diferenciable a través de la interfaz.

• Convergencia: Para una discretización espacial de orden $k$, los resultados numéricos confirman una convergencia de orden $k+1$ en la norma $L^2$ y de orden $k$ en la semi-norma $H^1$.
• Robustez: Los resultados se mantienen tanto para dominios 2D como 3D, incluso cuando las mallas se generan de forma independiente (en lugar de mediante refinamiento uniforme), lo que indica la estabilidad del método.
• Disponibilidad del Código: Los autores proporcionan el código de simulación para permitir la reproducibilidad y un uso adicional por parte de la comunidad.

Significado y Alcance
El artículo afirma proporcionar un marco teórico y práctico completo para resolver problemas de advección-difusión con interfaces móviles utilizando elementos finitos isoparamétricos de alto orden. El significado radica en cerrar la brecha entre la teoría abstracta de espacios evolutivos y la implementación práctica para problemas de interfaz con saltos de flujo.

Los autores mantienen un alcance modesto respecto al campo de velocidad: asumen que un campo de velocidad global y suave se da a priori. Declaran explícitamente que el trabajo no aborda cómo computar dicho campo de velocidad (por ejemplo, mediante re-parametrización para evitar el enredo de la malla), señalando que tales métodos existen en la literatura [8, 18, 19]. Además, aunque el método podría aplicarse teóricamente a sistemas de EDP acoplados donde la evolución del dominio es parte de la solución, el análisis actual no cubre tales problemas acoplados. El enfoque permanece estrictamente en la aproximación del problema de transmisión parabólico dado una evolución de dominio conocida.