On Runge-Kutta convolution quadrature based fractional variational integrators

Technische Samenvatting: Runge-Kutta Convolutie Kwadratuur Gebaseerde Fractionele Variationele Integratoren

Probleemstelling
Fractionele differentiaalvergelijkingen zijn essentieel voor het modelleren van systemen met langetermijngeheugen, zoals visco-elastische materialen en anomalieus diffusie. Hoewel fractionele calculus de variatiemethoden uitbreidt naar dissipatieve systemen (fractionele demping), blijft het construeren van numerieke integratoren van hoge orde voor deze systemen een uitdaging. Voorgaand werk door [23] combineerde backward-differentiation convolution quadrature (BDFCQ) met Galerkin-variatiemethoden voor het conservatieve deel van de Lagrangiaan. Echter, deze aanpak lijdt onder een "verzadigingseffect": ondanks het potentieel van BDFCQ voor hogere-orde nauwkeurigheid, blijven de resulterende fractionele variationele integratoren (FVIs) beperkt tot tweede-orde nauwkeurigheid. Deze beperking komt voort uit het feit dat de multistep-structuur van BDFCQ alleen werkt op de hoofd-knoopwaarden, waardoor de interne stage-informatie die inherent is aan hogere-orde Galerkin-discretisaties van de conservatieve actie, niet wordt benut. Daarom is er behoefte aan fractionele variationele integratoren waarbij de discretisatie van het fractionele deel structureel compatibel is met de interne stage-structuur van de conservatieve discretisatie.

Methodologie
De auteurs stellen een nieuwe klasse van Fractionele Variationele Integratoren (FVIs) voor die Runge-Kutta Convolutie Kwadratuur (RKCQ) combineren met hogere-orde Galerkin-variatiemethoden. De methodologie rust op drie pijlers:

1. Geaugmenteerd Variatiekader: Volgens het beperkte Hamilton-principe [22] wordt het systeem gemodelleerd door de toestandsvariabelen ($x$ en $y$) te verdubbelen om fractionele demping te incorporeren. Het actiefunctionaal wordt gesplitst in een conservatief deel ($L_{con}$) en een fractioneel deel ($L_{fra}$).
2. Discretisatiestrategie:
   • Conservatief Deel: Benaderd met hogere-orde Galerkin-variatiemethoden met polynomiale interpolatie van graad $s$ en Lobatto-kwadratuur.
   • Fractioneel Deel: Benaderd met RKCQ gebaseerd op Lobatto IIIC-methoden. In tegenstelling tot BDFCQ evalueert RKCQ de convolutie-integraal op de interne Runge-Kutta stages ($t_k + c_i h$). Dit maakt het mogelijk dat de benadering van de fractionele afgeleide perfect aansluit bij de interne knopen die worden gebruikt voor de Galerkin-discretisatie van de conservatieve actie.
3. Afleiding van Vergelijkingen: De auteurs leiden de Fractionele Discrete Euler-Lagrange (FDEL) vergelijkingen af door het discrete beperkte Hamilton-principe toe te passen. Cruciaal voor deze afleiding zijn de nieuwe discrete identiteiten die zijn vastgesteld voor RKCQ:
   • Een discrete semigroep-eigenschap voor de convolutie-gewichten (Lemma 4.1).
   • Asymmetrische discrete integratie-door-delen-formules (Lemma's 4.2 en 4.3), die noodzakelijk zijn om de niet-lokale aard van de fractionele afgeleiden binnen het variationele kader te behandelen.

Belangrijkste Bijdragen
Het artikel levert vier primaire bijdragen:

1. Constructie van Compatibele Integratoren: De auteurs construeren FVIs waarbij de stage-structuur van de RKCQ (specifiek Lobatto IIIC) overeenkomt met de interne knopen van de Galerkin-discretisatie. Deze structurele compatibiliteit voorkomt de orde-verzadiging die wordt waargenomen bij eerdere BDFCQ-gebaseerde benaderingen.
2. Theoretische Fundamenten: Het artikel stelt de noodzakelijke discrete identiteiten (semigroep-eigenschap en integratie-door-delen) vast die vereist zijn om de FDEL-vergelijkingen binnen de RKCQ-setting af te leiden.
3. Convergentieanalyse:
   • Theoretisch: Een rigoureus bewijs van tweede-orde convergentie wordt geleverd voor het 2-stage Lobatto IIIC-schema toegepast op een kwadratische Lagrangiaan met $\beta = 2\alpha \in (0, 1]$.
   • Numeriek: Numerieke experimenten tonen aan dat de methode hogere-orde nauwkeurigheid (4e en 6e orde) bereikt voor 3- en 4-stage Lobatto IIIC-schema's, respectievelijk, wanneer deze wordt toegepast op gladde problemen.
4. Midpoint-schema: Een korte discussie en implementatie van een midpoint-gebaseerde fractionele integrator (MIDCQ) worden opgenomen, die eveneens tweede-orde nauwkeurigheid bereikt.

Resultaten
De voorgestelde schema's werden getest op twee modellen: een tweedimensionale geforceerde harmonische oscillator en de fractionele Bagley-Torvik vergelijking.

• Nauwkeurigheid: Voor het 2-stage Lobatto IIIC-schema bevestigden de numerieke resultaten de theoretische tweede-orde convergentie. Voor 3- en 4-stage schema's waren de geobserveerde convergentieordes respectievelijk 4 en 6, wat overeenkomt met de theoretische verwachtingen voor de onderliggende Galerkin-methoden.
• Energie-afname: De integratoren reproduceerden succesvol de verwachte energie-afname gedrag van fractioneel gedempte systemen. In langetermijnsimulaties behield de voorgestelde FVI betere energie-afname karakteristieken vergeleken met klassieke RK4-methoden.
• Regulariteitsafhankelijkheid: De studie benadrukte dat hogere-orde convergentie afhankelijk is van de regulariteit van de oplossing. Voor problemen met discontinue forcering (bijv. stapfuncties) verzadigde de convergentieorde bij de eerste orde, wat consistent is met de regulariteitseisen van de onderliggende RKCQ-theorie. Echter, voor gladde forcering werd de volledige hogere-orde nauwkeurigheid gerealiseerd.
• Vergelijking met BDFCQ: De resultaten bevestigen dat de RKCQ-aanpak de ordebeperking van de BDFCQ-aanpak overwint door de interne stage-informatie volledig te benutten.

Betekenis en Claims
Het artikel claimt dat de primaire betekenis van dit werk ligt in de structurele compatibiliteit tussen de benadering van de fractionele afgeleide en de conservatieve discretisatie. Door gebruik te maken van RKCQ gebaseerd op Lobatto IIIC, incorporeert de methode op natuurlijke wijze dezelfde stage-knopen die door de Galerkin discrete Lagrangiaan worden gebruikt. Dit maakt het mogelijk om zowel de conservatieve als de fractionele termen op een gemeenschappelijk rooster te evalueren, waardoor de orde-verzadiging die de multistep-gebaseerde benaderingen zoals BDFCQ teistert, wordt vermeden.

De auteurs merken bescheiden op dat hoewel numeriek bewijs 4e en 6e-orde convergentie ondersteunt voor hogere-orde schema's, een rigoureus theoretisch bewijs voor deze hogere ordes een openstaand probleem blijft. Daarnaast is de computationele kosten van de huidige implementatie $O(N^2)$ vanwege de geschiedenis-afhankelijkheid van de convolutie-gewichten, hoewel de auteurs erkennen dat snelle convolutie-technieken dit in de toekomst potentieel kunnen reduceren naar $O(N \log^2 N)$. Het artikel concludeert dat deze aanpak een robuust kader biedt voor het construeren van hogere-orde, structuur-preserverende integratoren voor fractioneel gedempte Lagrangiaanse systemen.