On Runge-Kutta convolution quadrature based fractional variational integrators

Résumé Technique : Intégrateurs Variationnels Fractionnaires Basés sur la Quadrature de Convolution de Runge-Kutta

Énoncé du Problème
Les équations différentielles fractionnaires sont essentielles pour modéliser des systèmes dotés d'une mémoire à long terme, tels que les matériaux viscoélastiques et la diffusion anomale. Bien que le calcul fractionnaire étende les principes variationnels aux systèmes dissipatifs (amortissement fractionnaire), la construction d'intégrateurs numériques d'ordre élevé pour ces systèmes reste un défi. Des travaux antérieurs [23] ont combiné la quadrature de convolution par différenciation de type BDF (BDFCQ) avec des méthodes variationnelles de Galerkin pour la partie conservatrice du lagrangien. Cependant, cette approche souffre d'un effet de « saturation » : malgré le potentiel de précision d'ordre élevé de la BDFCQ, les intégrateurs variationnels fractionnaires (FVI) qui en résultent sont limités à une précision de second ordre. Cette limitation provient de la structure multi-pas de la BDFCQ qui n'opère que sur les valeurs nodales principales, échouant ainsi à utiliser l'information des étages internes inhérente aux discrétisations de Galerkin d'ordre supérieur de l'action conservatrice. Par conséquent, il existe un besoin pour des intégrateurs variationnels fractionnaires où la discrétisation de la partie fractionnaire est structurellement compatible avec la structure d'étages internes de la discrétisation conservatrice.

Méthodologie
Les auteurs proposent une nouvelle classe d'Intégrateurs Variationnels Fractionnaires (FVI) qui combinent la Quadrature de Convolution de Runge-Kutta (RKCQ) avec des méthodes variationnelles de Galerkin d'ordre supérieur. La méthodologie repose sur trois piliers :

1. Cadre Variationnel Augmenté : Suivant le principe de Hamilton restreint [22], le système est modélisé en doublant les variables d'état ($x$ et $y$) pour incorporer l'amortissement fractionnaire. Le fonctionnel d'action est divisé en une partie conservatrice ($L_{con}$) et une partie fractionnaire ($L_{fra}$).
2. Stratégie de Discrétisation :
   • Partie Conservatrice : Approchée par des méthodes variationnelles de Galerkin d'ordre supérieur avec une interpolation polynomiale de degré $s$ et une quadrature de Lobatto.
   • Partie Fractionnaire : Approchée par RKCQ basée sur des méthodes Lobatto IIIC. Contrairement à la BDFCQ, la RKCQ évalue l'intégrale de convolution aux étages de Runge-Kutta internes ($t_k + c_i h$). Cela permet à l'approximation de la dérivée fractionnaire de s'aligner parfaitement avec les nœuds internes utilisés pour la discrétisation de Galerkin de l'action conservatrice.
3. Dérivation des Équations : Les auteurs dérivent les équations de l'Euler-Lagrange Discrètes Fractionnaires (FDEL) en appliquant le principe de Hamilton restreint discret. La clé de cette dérivation réside dans les nouvelles identités discrètes établies pour la RKCQ :
   • Une propriété de semi-groupe discrète pour les poids de convolution (Lemme 4.1).
   • Des formules d'intégration par parties discrètes asymétriques (Lemmes 4.2 et 4.3), nécessaires pour traiter la nature non locale des dérivées fractionnaires au sein du cadre variationnel.

Contributions Clés
Le document présente quatre contributions principales :

1. Construction d'Intégrateurs Compatibles : Les auteurs construisent des FVI où la structure d'étage de la RKCQ (spécifiquement Lobatto IIIC) correspond aux nœuds internes de la discrétisation de Galerkin. Cette compatibilité structurelle évite la saturation d'ordre observée dans les approches précédentes basées sur la BDFCQ.
2. Fondements Théoriques : Le papier établit les identités discrètes nécessaires (propriété de semi-groupe et intégration par parties) requises pour dériver les équations FDEL dans le cadre de la RKCQ.
3. Analyse de Convergence :
   • Théorique : Une preuve rigoureuse de convergence de second ordre est fournie pour le schéma Lobatto IIIC à 2 étages appliqué à un lagrangien quadratique avec $\beta = 2\alpha \in (0, 1]$.
   • Numérique : Des expériences numériques démontrent que la méthode atteint une précision d'ordre supérieur (4 et 6 pour les schémas Lobatto IIIC à 3 et 4 étages, respectivement) lorsqu'elle est appliquée à des problèmes lisses.
4. Schéma de Point Milieu : Une brève discussion et une implémentation d'un intégrateur fractionnaire basé sur le point milieu (MIDCQ) sont incluses, lequel atteint également une précision de second ordre.

Résultats
Les schémas proposés ont été testés sur deux modèles : un oscillateur harmonique forcé en deux dimensions et l'équation de Bagley-Torvik fractionnaire.

• Précision : Pour le schéma Lobatto IIIC à 2 étages, les résultats numériques ont confirmé la convergence de second ordre théorique. Pour les schémas à 3 et 4 étages, les ordres de convergence observés étaient respectivement de 4 et 6, correspondant aux attentes théoriques des méthodes de Galerkin sous-jacentes.
• Décroissance de l'Énergie : Les intégrateurs ont reproduit avec succès le comportement de décroissance d'énergie attendu des systèmes fractionnairement amortis. Dans les simulations à long terme, le FVI proposé a maintenu de meilleures caractéristiques de décroissance d'énergie par rapport aux méthodes RK4 classiques.
• Dépendance à la Régularité : L'étude a souligné que la convergence d'ordre élevé est contingente sur la régularité de la solution. Pour les problèmes avec des forçages discontinus (ex: fonctions échelon), l'ordre de convergence a saturé à premier ordre, ce qui est cohérent avec les exigences de régularité de la théorie RKCQ sous-jacente. Cependant, pour des forçages lisses, la pleine précision d'ordre élevé a été réalisée.
• Comparaison avec la BDFCQ : Les résultats confirment que l'approche RKCQ surmonte la limitation d'ordre de l'approche BDFCQ en utilisant pleinement l'information des étages internes.

Signification et Revendications
Le papier affirme que la principale signification de ce travail réside dans la compatibilité structurelle entre l'approximation de la dérivée fractionnaire et la discrétisation conservatrice. En utilisant la RKCQ basée sur Lobatto IIIC, la méthode incorpore naturellement les mêmes nœuds d'étage utilisés par le lagrangien discret de Galerkin. Cela permet d'évaluer les termes conservatifs et fractionnaires sur une grille commune, évitant ainsi la saturation d'ordre qui affecte les approches basées sur les multi-pas comme la BDFCQ.

Les auteurs notent modestement que, bien que les preuves numériques soutiennent une convergence d'ordre 4 et 6 pour les schémas à plus d'étages, une preuve théorique rigoureuse pour ces ordres supérieurs demeure un problème ouvert. De plus, le coût computationnel de l'implémentation actuelle est de $O(N^2)$ en raison de la dépendance historique des poids de convolution, bien que les auteurs reconnaissent que des techniques de convolution rapide pourraient potentiellement réduire cela à $O(N \log^2 N)$ dans des travaux futurs. Le papier conclut que cette approche fournit un cadre robuste pour la construction d'intégrateurs d'ordre élevé et préservant la structure pour les systèmes lagrangiens fractionnairement amortis.