On Runge-Kutta convolution quadrature based fractional variational integrators

Resumen Técnico: Integradores Variacionales Fraccionarios Basados en la Cuadratura de Convolución de Runge-Kutta

Planteamiento del Problema
Las ecuaciones diferenciales fraccionarias son esenciales para modelar sistemas con memoria a largo plazo, como materiales viscoelásticos y la difusión anómala. Si bien el cálculo fraccionario extiende los principios variacionales para incluir sistemas disipativos (amortiguamiento fraccionario), la construcción de integradores numéricos de alto orden para estos sistemas sigue siendo un desafío. Un trabajo previo por [23] combinó la cuadratura de convolución de diferenciación hacia atrás (BDFCQ) con métodos variacionales de Galerkin para la parte conservativa del Lagrangiano. Sin embargo, este enfoque sufre de un efecto de "saturación": a pesar del potencial de BDRCQ para alcanzar una precisión de orden superior, los integradores variacionales fraccionarios (FVI) resultantes están limitados a una precisión de segundo orden. Esta limitación surge porque la estructura de pasos múltiples de BDFCQ opera únicamente sobre los valores nodales principales, fallando al no utilizar la información de las etapas internas inherente a las discretizaciones de Galerkin de orden superior de la acción conservativa. En consecuencia, existe la necesidad de integradores variacionales fracionarios donde la discretización de la parte fraccionaria sea estructuralmente compatible con la estructura de etapas internas de la discretización conservativa.

Metodología
Los autores proponen una nueva clase de Integradores Variacionales Fraccionarios (FVI) que combinan la Cuadratura de Convolución de Runge-Kutta (RKCQ) con métodos variacionales de Galerkin de orden superior. La metodología se basa en tres pilares:

1. Marco Variacional Aumentado: Siguiendo el principio de Hamilton restringido [22], el sistema se modela duplicando las variables de estado ($x$ e $y$) para incorporar el amortiguamiento fraccionario. El funcional de acción se divide en una parte conservativa ($L_{con}$) y una parte fraccionaria ($L_{fra}$).
2. Estrategia de Discretización:
   • Parte Conservativa: Aproximada mediante métodos variacionales de Galerkin de orden superior con interpolación polinómica de grado $s$ y cuadratura de Lobatto.
   • Parte Fraccionaria: Aproximada mediante RKCQ basada en métodos Lobatto IIIC. A diferencia de BDFCQ, la RKCQ evalúa la integral de convolución en las etapas internas de Runge-Kutta ($t_k + c_i h$). Esto permite que la aproximación de la derivada fraccionaria se alinee perfectamente con los nodos internos utilizados en la discretización de Galerkin de la acción conservativa.
3. Derivación de las Ecuaciones: Los autores derivan las ecuaciones de Euler-Lagrange Discretas Fraccionarias (FDEL) aplicando el principio de Hamilton restringido discreto. Clave en esta derivación son las nuevas identidades discretas establecidas para RKCQ:
   • Una propiedad de semigrupo discreto para los pesos de convolución (Lema 4.1).
   • Fórmulas de integración por partes discretas asimétricas (Lemmas 4.2 y 4.3), que son necesarias para manejar la naturaleza no local de las derivadas fraccionarias dentro del marco variacional.

Contribuciones Clave
El artículo realiza cuatro contribuciones primarias:

1. Construcción de Integradores Compatibles: Los autores construyen FVIs donde la estructura de etapas de la RKCQ (específicamente Lobatto IIIC) coincide con los nodos internos de la discretización de Galerkin. Esta compatibilidad estructural evita la saturación de orden observada en los enfoques previos basados en BDFCQ.
2. Fundamentos Teóricos: El artículo establece las identidades discretas necesarias (propiedad de semigrupo e integración por partes) requeridas para derivar las ecuaciones FDEL dentro del entorno RKCQ.
3. Análisis de Convergencia:
   • Teórica: Se proporciona una prueba rigurosa de convergencia de segundo orden para el esquema Lobatto IIIC de 2 etapas aplicado a un Lagrangiano cuadrático con $\beta = 2\alpha \in (0, 1]$.
   • Numérica: Experimentos numéricos demuestran que el método alcanza una precisión de orden superior (4.º y 6.º orden) para esquemas Lobatto IIIC de 3 y 4 etapas, respectivamente, cuando se aplica a problemas suaves.
4. Esquema de Punto Medio: Se incluye una breve discusión e implementación de un integrador fraccionario basado en el punto medio (MIDCQ), que también logra una precisión de segundo orden.

Resultados
Los esquemas propuestos fueron probados en dos modelos: un oscilador armónico forzado bidimensional y la ecuación de Bagley-Torvik fraccionaria.

• Precisión: Para el esquema Lobatto IIIC de 2 etapas, los resultados numéricos confirmaron la convergencia de segundo orden teórica. Para los esquemas de 3 y 4 etapas, los órdenes de convergencia observados fueron 4 y 6, respectivamente, coincidiendo con las expectativas teóricas para los métodos de Galerkin subyacentes.
• Decaimiento de Energía: Los integradores reprodujeron con éxito el comportamiento esperado de decaimiento de energía de los sistemas con amortiguamiento fraccionario. En simulaciones de largo plazo, el FVI propuesto mantuvo mejores características de decaimiento de energía en comparación con los métodos RK4 clásicos.
• Dependencia de la Regularidad: El estudio destacó que la convergencia de alto orden está supeditada a la regularidad de la solución. Para problemas con forzamientos discontinuos (por ejemplo, funciones escalón), el orden de convergencia se saturó en primer orden, lo cual es consistente con los requisitos de regularidad de la teoría de RKCQ subyacente. Sin embargo, para forzamientos suaves, se alcanzó la plena precisión de orden superior.
• Comparación con BDFCQ: Los resultados confirman que el enfoque RKCQ supera la limitación de orden del enfoque BDFCQ al utilizar plenamente la información de las etapas internas.

Significancia y Reivindicaciones
El artículo sostiene que la principal significancia de este trabajo reside en la compatibilidad estructural entre la aproximación de la derivada fraccionaria y la discretización conservativa. Al utilizar RKCQ basada en Lobatto IIIC, el método incorpora naturalmente los mismos nodos de etapa utilizados por el Lagrangiano discreto de Galerkin. Esto permite que tanto los términos conservativos como los fraccionarios se evalúen en una malla común, evitando así la saturación de orden que afecta a los enfoques basados en pasos múltiples como BDFCQ.

Los autores señalan modestamente que, si bien la evidencia numérica respalda la convergencia de 4.º y 6.º orden para esquemas de más etapas, una prueba teórica rigurosa para estos órdenes superiores sigue siendo un problema abierto. Además, el costo computacional de la implementación actual es de $O(N^2)$ debido a la dependencia histórica de los pesos de convolución, aunque los autores reconocen que técnicas de convolución rápida podrían reducir esto a $O(N \log^2 N)$ en trabajos futuros. El artículo concluye que este enfoque proporciona un marco robusto para la construcción de integradores de alto orden y preservadores de estructura para sistemas lagrangianos con amortiguamiento fraccionario.