On a Boundary-Initial Value Problem for Fractional Differential Equation with Sequential Caputo derivatives

Resumen Técnico: Problema de Valor en la Frontera e Inicial para Ecuaciones Diferenciales Fraccionarias con Derivadas Caputo Secuenciales

1. Planteamiento del Problema

El artículo aborda un problema de valor en la frontera e inicial (BIVP) para una ecuación diferencial fraccionaria que involucra derivadas Caputo secuenciales. El modelo matemático estudiado es:

$$D_t^\beta (D_t^\alpha u(x, t)) + D_t^\beta u(x, t) - u_{xx}(x, t) = f(x, t)$$

definido en el dominio $\Omega = \{(x, t) : 0 < x < 1, 0 < t < T\}$, sujeto a:

• Condiciones iniciales: $u(x, 0) = \phi(x)$ y $D_t^\alpha u(x, 0) = \psi(x)$.
• Condiciones de frontera: $u(0, t) = u(1, t) = 0$.

Donde $D_t^\alpha$ y $D_t^\beta$ son derivadas fraccionarias de Caputo de orden $0 < \alpha, \beta < 1$. El término secuencial $D_t^\beta (D_t^\alpha u)$ es crucial, ya que, a diferencia de los operadores fraccionarios estándar, la composición de estos operadores no es conmutativa ni satisface la ley de índices clásica ($D^\beta D^\alpha \neq D^{\alpha+\beta}$), lo que introduce efectos de memoria complejos y múltiples escalas de tiempo. El estudio se centra en el régimen donde $\alpha + \beta > 1$ y $\alpha > \beta$.

2. Metodología

Los autores emplean un enfoque híbrido que combina análisis funcional, teoría de funciones especiales y métodos numéricos discretos:

• Análisis Teórico (Método de Fourier):

  • Se utiliza el método de separación de variables (series de Fourier en senos) para transformar la EDP en una familia de problemas de valor inicial para ecuaciones diferenciales ordinarias (EDO) fraccionarias.
  • Se resuelve la EDO auxiliar utilizando la Transformada de Laplace, derivando la solución exacta en términos de la función de Mittag-Leffler bivariada ($E_2$).
  • Se establecen lemas sobre las propiedades de convergencia y cotas de la función $E_2$ (incluyendo representaciones integrales tipo Hankel) para probar la convergencia uniforme de las series de la solución y sus derivadas.
  • Se demuestra la existencia y unicidad de la solución regular bajo condiciones de regularidad de Hölder para los datos iniciales y la fuente.

• Implementación Numérica (Esquema L1-FEM):

  • Reformulación Secuencial: Se introduce una variable auxiliar $v(x, t) = D_t^\alpha u(x, t)$ para convertir el operador secuencial anidado en un sistema acoplado de primer orden, evitando la discretización directa de términos históricos complejos.
  • Discretización Espacial: Se utiliza el Método de Elementos Finitos (FEM) con elementos lineales ($P_1$) en una malla uniforme espacial.
  • Discretización Temporal: Se emplea un esquema L1 sobre una malla temporal graduada (graded mesh) definida por $t_n = T(n/N)^r$. El parámetro de graduación $r$ se elige para agrupar nodos cerca de $t=0$, compensando la singularidad inicial típica de las soluciones fraccionarias.
  • Se desarrolla un algoritmo paso a paso (Algoritmo 1) para resolver el sistema lineal en cada paso de tiempo.

3. Contribuciones Clave

1. Solución Analítica Exacta: Derivación de la solución explícita del problema en términos de la función de Mittag-Leffler bivariada, una herramienta menos común pero esencial para manejar la estructura secuencial de los operadores.
2. Propiedades de la Función $E_2$: Formulación y demostración de propiedades críticas (convergencia, cotas asintóticas y reglas de diferenciación/integración) de la función de Mittag-Leffler bivariada, que son fundamentales para la prueba de existencia y unicidad.
3. Esquema Numérico Robusto: Desarrollo de un esquema totalmente discreto que combina la reformulación secuencial con el método L1-FEM en mallas graduadas, logrando alta precisión incluso cerca de la singularidad inicial.
4. Análisis de Dinámica del Sistema: Identificación de un comportamiento transicional en la solución: a medida que aumentan los órdenes fraccionarios, el sistema pasa de un decaimiento puramente difusivo a un régimen oscilatorio amortiguado con un "rebote" de energía, indicando una inercia fraccionaria.

4. Resultados Principales

• Teorema de Existencia y Unicidad: Se demuestra que, bajo las condiciones $\alpha + \beta > 1$, $\alpha > \beta$ y regularidad adecuada de los datos ($\phi, \psi \in \mathring{C}^a_2$, $f \in \mathring{C}^a_2(\bar{\Omega})$), existe una única solución regular dada por la serie de Fourier convergente.
• Validación Numérica:
  • Se utilizó una solución manufacturada para verificar la precisión del método. Con una malla graduada ($r \approx 1.6$), el error máximo en la norma $L^2$ alcanzó $\approx 1.73 \times 10^{-6}$.
  • Comparación de Mallas: Se demostró que las mallas graduadas superan significativamente a las uniformes. Mientras que una malla uniforme ofrece una tasa de convergencia subóptima ($\approx 0.75$), la malla graduada restaura la tasa de convergencia a valores superiores a $1.10$, reduciendo el error en un factor de ~7.75 para un número fijo de nodos temporales.
• Estudio Paramétrico:
  • Los parámetros $\alpha$ y $\beta$ no solo afectan la velocidad de decaimiento, sino que reconfiguran fundamentalmente la forma de la evolución de la solución (cruce de curvas en el tiempo).
  • Se observó un efecto de "rebote" en la energía del sistema para valores altos de $\alpha$ y $\beta$, donde la energía cruza el equilibrio y oscila ligeramente antes de estabilizarse.
  • Se identificó una compensación ("trade-off") entre $\alpha$ y $\beta$: diferentes combinaciones de estos parámetros pueden producir el mismo resultado en un punto específico, lo cual es vital para problemas de identificación inversa.

5. Significado e Impacto

Este trabajo es significativo porque:

• Avanza la teoría de operadores secuenciales: Proporciona un marco riguroso para resolver problemas con derivadas secuenciales, que son comunes en modelos de física con múltiples efectos de memoria (como ecuaciones de difusión-onda o telegrafía fraccionaria).
• Resuelve desafíos de regularidad: La propuesta de usar mallas graduadas en combinación con la reformulación secuencial ofrece una solución práctica y eficiente a la baja regularidad temporal en $t=0$, un obstáculo común en la simulación numérica de EDPs fraccionarias.
• Aplicabilidad en Modelado: Los hallazgos sobre la dinámica de "rebote" y la sensibilidad a los parámetros $\alpha$ y $\beta$ ofrecen nuevas perspectivas para modelar sistemas físicos complejos donde la inercia y la memoria juegan roles no triviales, facilitando la calibración de modelos en ingeniería y ciencias naturales.

En resumen, el artículo establece un puente sólido entre la teoría analítica de funciones especiales avanzadas y la implementación numérica eficiente, ofreciendo herramientas tanto para la comprensión teórica como para la simulación práctica de sistemas dinámicos gobernados por derivadas fraccionarias secuenciales.