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

Titolo: Su un Problema ai Limiti-Iniziale per Equazioni Differenziali frazionarie con Derivate di Caputo Sequenziali

1. Il Problema Studiato

Il lavoro si concentra sullo studio di un problema ai limiti-iniziale (BIVP) per un'equazione differenziale frazionaria che coinvolge derivate di Caputo sequenziali. L'equazione principale, definita nel dominio $\Omega = \{(x, t) : 0 < x < 1, 0 < t < T\}$, è data da:

$$D^\beta_t (D^\alpha_t u(x, t)) + D^\beta_t u(x, t) - u_{xx}(x, t) = f(x, t)$$

dove:

• $D^\alpha_t$ e $D^\beta_t$ sono derivate di Caputo di ordine $0 < \alpha < 1$ e $0 < \beta < 1$.
• Il termine $D^\beta_t (D^\alpha_t u)$ rappresenta la derivata frazionaria sequenziale.
• $f(x, t)$ è una funzione sorgente assegnata.

Il problema è soggetto a:

• Condizioni iniziali: $u(x, 0) = \phi(x)$ e $D^\alpha_t u(x, 0) = \psi(x)$.
• Condizioni al bordo: $u(0, t) = u(1, t) = 0$.

L'obiettivo è dimostrare l'esistenza e l'unicità della soluzione regolare, derivarne la forma analitica esatta e sviluppare uno schema numerico efficiente per la sua risoluzione.

2. Metodologia

Gli autori adottano un approccio ibrido che combina analisi teorica rigorosa e implementazione numerica avanzata.

A. Analisi Teorica:

• Metodo di Fourier: La soluzione viene cercata sotto forma di serie di seni (sviluppo in autofunzioni), trasformando l'equazione alle derivate parziali in una successione di problemi di Cauchy per equazioni differenziali ordinarie (ODE) frazionarie.
• Funzioni di Mittag-Leffler Bivariate: Per risolvere le ODE risultanti, gli autori utilizzano la funzione di Mittag-Leffler bivariate ($E_2$). Viene dimostrata l'esistenza di una soluzione esplicita in termini di questa funzione, che generalizza le funzioni di Mittag-Leffler classiche e Prabhakar.
• Stime di Convergenza: Viene stabilita la convergenza uniforme delle serie ottenute per la soluzione e le sue derivate (temporali e spaziali) utilizzando le proprietà asintotiche della funzione $E_2$ e disuguaglianze di tipo Hölder.

B. Implementazione Numerica:

• Riformulazione Sequenziale: Per gestire la complessità della derivata sequenziale $D^\beta_t (D^\alpha_t u)$, viene introdotta una variabile ausiliaria $v(x, t) = D^\alpha_t u(x, t)$. Questo trasforma il problema in un sistema accoppiato di equazioni, evitando la discretizzazione diretta di termini storici nidificati.
• Metodo agli Elementi Finiti (FEM): Viene utilizzata una discretizzazione spaziale con elementi lineari (P1) conformi.
• Schema L1 su Griglia Graduata: Per la discretizzazione temporale, viene impiegato lo schema L1 per le derivate di Caputo. Poiché le soluzioni di equazioni frazionarie presentano spesso singolarità in $t=0$, viene adottata una griglia temporale graduata ($t_n = T(n/N)^r$) con un parametro di gradazione $r \ge 1$ (specificamente $r = 2 - \min\{\alpha, \beta\}$) per migliorare l'accuratezza vicino all'origine temporale.

3. Contributi Chiave

1. Soluzione Analitica Esatta: Derivazione della soluzione esatta del problema BIVP in termini di una serie di autofunzioni moltiplicate per funzioni di Mittag-Leffler bivariate, fornendo una rappresentazione chiusa precedentemente non disponibile per questo tipo specifico di operatori sequenziali.
2. Dimostrazione di Esistenza e Unicità: Prova rigorosa dell'esistenza e unicità della soluzione regolare sotto le condizioni $\alpha + \beta > 1$ e $\alpha > \beta$, utilizzando tecniche di analisi funzionale e convergenza di serie.
3. Nuovo Schema Numerico Robusto: Sviluppo di un algoritmo completamente discreto che combina FEM spaziale, schema L1 temporale e griglie graduate, specificamente ottimizzato per gestire la natura non commutativa e le singolarità delle derivate sequenziali.
4. Proprietà della Funzione $E_2$: Formulazione e utilizzo di nuove proprietà e stime per la funzione di Mittag-Leffler bivariate, essenziali sia per la costruzione della soluzione analitica che per la validazione numerica.

4. Risultati

• Validazione Numerica: Lo schema numerico è stato validato utilizzando una soluzione "manufactured" (costruita artificialmente). I risultati mostrano un errore massimo nell'ordine di $10^{-6}$, confermando l'alta precisione del metodo.
• Confronto Griglie Uniformi vs Graduata: Gli esperimenti numerici dimostrano che l'uso di una griglia temporale uniforme ($r=1$) porta a tassi di convergenza subottimali (circa 0.75). Al contrario, l'uso di una griglia graduata ($r \approx 1.6$) ripristina il tasso di convergenza teorico (superiore a 1.1), riducendo l'errore di un fattore di circa 7.75 rispetto alla griglia uniforme a parità di nodi temporali.
• Comportamento Fisico: L'analisi parametrica rivela che gli ordini frazionari $\alpha$ e $\beta$ non influenzano solo la velocità di decadimento, ma modificano qualitativamente la dinamica del sistema. Per valori elevati di $\alpha$ e $\beta$, si osserva una transizione da un decadimento puramente diffusivo a un regime oscillatorio smorzato, con un fenomeno di "rimbalzo" dell'energia (il sistema attraversa lo stato di equilibrio), suggerendo un comportamento inerziale.
• Trade-off Parametrico: È stato identificato che diverse combinazioni di $\alpha$ e $\beta$ possono produrre lo stesso risultato finale, offrendo flessibilità nella modellazione di fenomeni reali.

5. Significato e Impatto

Questo lavoro colma un divario significativo nella letteratura sulle equazioni differenziali frazionarie con operatori sequenziali, che sono cruciali per modellare sistemi con effetti multipli di memoria o scale temporali multiple.

• Teorico: Fornisce un fondamento matematico solido per l'analisi di problemi di valore iniziale-bordo con derivate sequenziali, superando le limitazioni delle leggi degli indici classiche.
• Pratico: Lo schema numerico proposto offre uno strumento affidabile per simulare sistemi fisici complessi (come equazioni del telegrafo frazionarie o processi di diffusione anomala) dove la memoria sequenziale è dominante. La capacità di gestire le singolarità iniziali attraverso griglie graduate rende il metodo superiore rispetto agli approcci standard per una vasta gamma di applicazioni ingegneristiche e fisiche.