On Runge-Kutta convolution quadrature based fractional variational integrators

Sintesi Tecnica: Integratori Variazionali Frazionari basati su Integrazione di Convoluzione Runge-Kutta

Problematica
Le equazioni differenziali frazionarie sono essenziali per modellare sistemi con memoria a lungo termine, come i materiali viscoelastici e la diffusione anomala. Sebbene il calcolo frazionario estenda i principi variazionali ai sistemi dissipativi (smorzamento frazionario), la costruzione di integratori numerici ad alto ordine per questi sistemi rimane una sfida. Il lavoro precedente di [23] ha combinato la quadratura di convoluzione basata su differenze retrostanti (BDFCQ) con metodi variazionali di tipo Galerkin per la parte conservativa del Lagrangiano. Tuttavia, questo approccio soffre di un effetto di "saturazione": nonostante il potenziale di accuratezza di ordine superiore della BDFCQ, gli integratori variazionali frazionari (FVI) risultanti sono limitati a un'accuratezza del secondo ordine. Questa limitazione deriva dal fatto che la struttura multistep della BDFCQ opera solo sui valori nodali principali, fallendo nel utilizzare l'informazione degli stadi interni inerente alle discretizzazioni di ordine superiore dei metodi di Galerkin per l'azione conservativa. Di conseguenza, vi è la necessità di integratori variazionali frazionari in cui la discretizzazione della parte frazionaria sia strutturalmente compatibile con la struttura degli stadi interni della discretizzazione conservativa.

Metodologia
Gli autori propongono una nuova classe di Integratori Variazionali Frazionari (FVI) che combinano la Quadratura di Convoluzione Runge-Kutta (RKCQ) con metodi variazionali di Galerkin ad ordine superiore. La metodologia si basa su tre pilastri:

1. Framework Variazionale Aumentato: Seguendo il principio di Hamilton ristretto [22], il sistema è modellato raddoppiando le variabili di stato ($x$ e $y$) per incorporare lo smorzamento frazionario. Il funzionale dell'azione è diviso in una parte conservativa ($L_{con}$) e una parte frazionaria ($L_{fra}$).
2. Strategia di Discretizzazione:
   • Parte Conservativa: Approssimata utilizzando metodi variazionali di Galerkin ad ordine superiore con interpolazione polinomiale di grado $s$ e quadratura di Lobatto.
   • Parte Frazionaria: Approssimata mediante RKCQ basata su metodi Lobatto IIIC. A differenza della BDFCQ, la RKCQ valuta l'integrale di convoluzione ai nodi interni di Runge-Kutta ($t_k + c_i h$). Ciò consente all'approssimazione della derivata frazionaria di allinearsi perfettamente con i nodi interni utilizzati nella discretizzazione di Galerkin dell'azione conservativa.
3. Derivazione delle Equazioni: Gli autori derivano le equazioni di Euler-Lagrange discrete frazionarie (FDEL) applicando il principio di Hamilton ristretto discreto. Fondamentali per questa derivazione sono le nuove identità discrete stabilite per la RKCQ:
   • Una proprietà semigruppo discreta per i pesi di convoluzione (Lemma 4.1).
   • Formule di integrazione per parti discrete asimmetriche (Lemmi 4.2 e 4.3), necessarie per gestire la natura non locale delle derivate frazionarie all'interno del framework variazionale.

Contributi Chiave
Il documento presenta quattro contributi primari:

1. Costruzione di Integratori Compatibili: Gli autori costruiscono FVI in cui la struttura di stadio della RKCQ (specificamente Lobatto IIIC) corrisponde ai nodi interni della discretizzazione di Galerkin. Questa compatibilità strutturale evita la saturazione dell'ordine osservata nei precedenti approcci basati su BDFCQ.
2. Fondamenti Teorici: Il lavoro stabilisce le identità discrete necessarie (proprietà semigruppo e integrazione per parti) richieste per derivare le equazioni FDEL nel contesto della RKCQ.
3. Analisi della Convergenza:
   • Teorica: Viene fornita una prova rigorosa di convergenza del secondo ordine per lo schema Lobatto IIIC a 2 stadi applicato a un Lagrangiano quadratico con $\beta = 2\alpha \in (0, 1]$.
   • Numerica: Esperimenti numerici dimostrano che il metodo raggiunge un'accuratezza di ordine superiore (4° e 6° ordine) per schemi Lobatto IIIC a 3 e 4 stadi, rispettivamente, quando applicato a problemi regolari.
4. Schema Midpoint: Viene inclusa una breve discussione e implementazione di un integratore frazionario basato sul punto medio (MIDCQ), che raggiunge anch'esso un'accuratezza del secondo ordine.

Risultati
Gli schemi proposti sono stati testati su due modelli: un oscillatore armonico forzato bidimensionale e l'equazione di Bagley-Torvik frazionaria.

• Accuratezza: Per lo schema Lobatto IIIC a 2 stadi, i risultati numerici hanno confermato la convergenza del secondo ordine prevista. Per gli schemi a 3 e 4 stadi, gli ordini di convergenza osservati sono stati rispettivamente 4 e 6, in linea con le aspettative teoriche per i sottostanti metodi di Galerkin.
• Decadimento dell'Energia: Gli integratori hanno riprodotto con successo il previsto comportamento di decadimento dell'energia dei sistemi con smorzamento frazionario. Nelle simulazioni a lungo termine, il proposto FVI ha mantenuto migliori caratteristiche di decadimento dell'energia rispetto ai classici metodi RK4.
• Dipendenza dalla Regolarità: Lo studio ha evidenziato che l'alta convergenza di ordine è contingente sulla regolarità della soluzione. Per problemi con forzanti discontinue (es. funzioni a gradino), l'ordine di convergenza si è saturato al primo ordine, in linea con i requisiti di regolarità della teoria RKCQ sottostante. Tuttavia, per forzanti regolari, è stata realizzata la piena accuratezza di ordine superiore.
• Confronto con BDFCQ: I risultati confermano che l'approccio RKCQ supera il limite di ordine della BDFCQ, utilizzando appieno l'informazione degli stadi interni.

Significatività e Rivendicazioni
Il lavoro sostiene che la principale significatività di questa ricerca risieda nella compatibilità strutturale tra l'approssimazione della derivata frazionaria e la discretizzazione conservativa. Utilizzando la RKCQ basata su Lobatto IIIC, il metodo incorpora naturalmente gli stessi nodi di stadio usati dal Lagrangiano discreto di Galerkin. Ciò permette sia ai termini conservativi che a quelli frazionari di essere valutati su una griglia comune, evitando così la saturazione dell'ordine che affligge gli approcci basati su multistep come la BDFCQ.

Gli autori notano con modestia che, sebbene l'evidenza numerica supporti la convergenza del 4° e 6° ordine per schemi a stadi superiori, una prova teorica rigorosa per questi ordini più elevati rimane un problema aperto. Inoltre, il costo computazionale dell'attuale implementazione è $O(N^2)$ a causa della dipendenza storica dei pesi di convoluzione, sebbene gli autori riconoscano che tecniche di convoluzione veloce potrebbero potenzialmente ridurre tale costo a $O(N \log^2 N)$ in lavori futuri. Il documento conclude che questo approccio fornisce un framework robusto per la costruzione di integratori ad alto ordine e preservanti la struttura per sistemi lagrangiani con smorzamento frazionario.