A numerical approach for solving the time-Fractional Mobile-Immobile Transport Equation

Resumen Técnico: Un enfoque numérico para resolver la ecuación de transporte móvil-inmóvil de tiempo fraccionario

Planteamiento del problema
El artículo aborda la solución numérica de la ecuación de transporte móvil-inmóvil (MIM) de tiempo fraccionario unidimensional, un modelo utilizado para describir la difusión anómala en medios porosos heterogéneos. A diferencia de los modelos clásicos de Fick, el modelo MIM tiene en cuenta el transporte de solutos retrasado causado por la transferencia de masa entre regiones móviles e inmóviles, gobernada por una derivada de Caputo de orden $\alpha \in (0, 1)$. El problema de valor inicial-frontera (IBVP) específico considerado es:
$$\lambda_1 \frac{\partial u}{\partial t} + \lambda_2 {}^C D^\alpha_{0,t} u + Lu = f(x,t),$$
donde $L$ es un operador elíptico espacial, y el sistema está sujeto a condiciones de contorno de Dirichlet cero. El desafío principal radica en la naturaleza no local de la derivada fraccionaria y la necesidad de una alta precisión tanto en el espacio como en el tiempo para capturar eficientemente los efectos de memoria y los comportamientos de transporte anómalo.

Metodología
El autor desarrolla dos esquemas numéricos totalmente discretos que combinan las discretizaciones espacial y temporal dentro de un marco unificado:

1. Discretización espacial: Se emplea el método de Galerkin discontinuo de penalización interior no simétrica (NIPG). Este enfoque utiliza un espacio de Sobolev roto y permite discontinuidades a través de los límites de los elementos. El operador espacial se aproxima mediante una forma bilineal $B_h(\cdot, \cdot)$ que incluye un término de penalización para asegurar la estabilidad y la coercitividad.
2. Discretización temporal: Se proponen dos estrategias distintas para aproximar las derivadas temporales (tanto la de primer orden como la de Caputo):
   • Esquema 1 (CN-L1-NIPG): Combina el método NIPG con un marco de Crank-Nicolson (CN) aplicado a la fórmula L1 estándar para la derivada de Caputo. Este esquema apunta a una precisión temporal de orden $(2-\alpha)$.
   • Esquema 2 (L2-1$\sigma$-NIPG): Combina el método NIPG con la fórmula L2-1$\sigma$ (introducida por Alikhanov) para la derivada de Caputo, también dentro de un marco CN. Este esquema busca una precisión temporal de segundo orden.

El análisis se basa en pruebas de estabilidad rigurosas utilizando normas de energía y estimaciones de error derivadas mediante proyecciones elípticas e desigualdades de tipo Grönwall discretas.

Contribuciones clave

• Esquemas novedosos: El artículo introduce dos formulaciones totalmente discretas distintas (CN-L1-NIPG y L2-1$\sigma$-NIPG) específicamente para el modelo MIM de tiempo fraccionario, una combinación que ha recibido poca atención en la literatura existente en comparación con las ecuaciones de difusión fraccionaria más simples.
• Análisis teórico:
  • Estabilidad: Se establece la estabilidad incondicional para ambos esquemas. El análisis demuestra que las soluciones numéricas permanecen acotadas por los datos iniciales y los términos fuente, independientemente del tamaño de la malla.
  • Convergencia: Se derivan estimaciones de error óptimas. El esquema CN-L1-NIPG logra un límite de error de $O(\tau^{2-\alpha} + h^{k+1})$, mientras que el esquema L2-1$\sigma$-NIPG logra $O(\tau^2 + h^{k+1})$, donde $\tau$ es el paso de tiempo, $h$ es el tamaño de la malla espacial y $k$ es el grado polinómico. Estas tasas dependen de la regularidad de la solución exacta.
• Límites de error exhaustivos: El análisis de error contabiliza explícitamente la regularidad de la solución, proporcionando una base teórica para las tasas de convergencia observadas.

Resultados
Se realizaron experimentos numéricos en un problema de prueba con una solución exacta conocida para validar las predicciones teóricas.

• CN-L1-NIPG: Los resultados numéricos confirmaron el orden de convergencia teórico de $(2-\alpha)$ en la dirección temporal para varios valores de $\alpha \in (0, 1)$.
• L2-1$\sigma$-NIPG: Los experimentos demostraron que este esquema logra una precisión temporal de segundo orden, superando al esquema basado en L1 en términos de tasa de convergencia para la derivada temporal.
• Robustez: Ambos esquemas demostraron eficiencia y robustez a través de diferentes órdenes fraccionarios y refinamientos de malla, mostrando el método L2-1$\sigma$ una precisión temporal superior.

Significancia y afirmaciones
El artículo afirma proporcionar un marco computacional fiable y preciso para simular procesos de transporte fraccionario complejos en medios porosos. Al combinar la flexibilidad del método NIPG con discretizaciones temporales de alto orden, los esquemas propuestos ofrecen un enfoque unificado para modelar fenómenos de transporte no Fickianos. El autor sostiene que su trabajo llena un vacío en la literatura con respecto a los métodos numéricos de alto orden para el modelo MIM, ofreciendo una herramienta que es tanto teóricamente sólida como prácticamente eficiente para aplicaciones en hidrología, transporte de contaminantes y dispersión de sedimentos. El estudio concluye que estos métodos allanan el camino para futuras extensiones a modelos de orden variable y de mayor dimensión.