Exact solution to a class of problems for the Burgers' equation on bounded intervals

Sintesi Tecnica: Soluzione Esatta per una Classe di Problemi dell'Equazione di Burgers su Intervalli Limitati

Enunciato del Problema
Il lavoro affronta l'equazione di Burgers viscosa, un'equazione differenziale alle derivate parziali non lineare che modella l'equilibrio tra convezione e diffusione, definita su intervalli limitati generici $[l_1, l_2]$. Il problema specifico consiste nel determinare la funzione $w(x, t)$ soggetta a:

1. Equazione Governativa: $w_t - a^2 w_{xx} + w w_x = 0$, dove $a^2$ rappresenta il coefficiente di viscosità cinematica ($1/Re$).
2. Condizione Iniziale: Un profilo iniziale non omogeneo $w(x, 0) = w_0(x)$.
3. Condizioni al Contorno: Condizioni al contorno di Dirichlet fisse $w(l_1, t) = \alpha_1$ e $w(l_2, t) = \alpha_2$.

Gli autori osservano che, sebbene i metodi numerici (differenze finite, elementi finiti, metodi spettrali) e le tecniche semi-analitiche (decomposizione di Adomian, perturbazione dell'omotopia) siano prevalenti, spesso mancano di soluzioni in forma chiusa o richiedono risorse computazionali significative. L'obiettivo è derivare una soluzione operativa esatta che offra un'espressione in forma chiusa nel dominio di Laplace e una soluzione precisa nel dominio del tempo tramite trasformate inverse.

Metodologia
La strategia di soluzione proposta si basa su un processo di trasformazione in due fasi:

1. Trasformazione di Hopf-Cole: L'equazione di Burgers non lineare viene linearizzata utilizzando la classica trasformazione di Hopf-Cole:
   $$w(x, t) = -2a^2 \frac{u_x(x, t)}{u(x, t)}$$
   Ciò converte il problema non lineare in un problema ai valori al contorno lineare per l'equazione di diffusione del calore unidimensionale (reazione-diffusione con termine di reazione nullo):
   $$u_t - a^2 u_{xx} = 0$$
   La condizione iniziale per $u$ è derivata dall'integrale della condizione iniziale $w_0$, e le condizioni di Dirichlet per $w$ si trasformano in condizioni al contorno di Robin omogenee per $u$.

2. Trasformata di Laplace e Soluzione Operativa: L'equazione del calore lineare è risolta nel dominio di Laplace (dominio $p$). Gli autori utilizzano una soluzione operativa esatta precedentemente stabilita per equazioni lineari di reazione-diffusione con coefficienti costanti su intervalli finiti.

   • Il problema è formulato come un problema ai valori al contorno a due punti nel dominio di Laplace.
   • Viene derivata un'espressione esatta per la trasformata di Laplace $U(x, p)$, che coinvolge termini esponenziali dipendenti da $\sqrt{p}$ e un sistema lineare di quattro equazioni per determinare i valori al contorno ($U(l_1, p), U_x(l_1, p), U(l_2, p), U_x(l_2, p)$).
   • Il termine sorgente nell'equazione linearizzata è nullo, semplificando la soluzione operativa in modo che dipenda dalla trasformata della condizione iniziale $R(x, p)$.

3. Recupero nel Dominio del Tempo: La soluzione esatta per la variabile originale $w(x, t)$ è ottenuta implicitamente tramite la trasformata inversa di Laplace:
   $$w(x, t) = -2a^2 \frac{\mathcal{L}^{-1}\{U_x(x, p)\}}{\mathcal{L}^{-1}\{U(x, p)\}}$$
   Il lavoro nota che, sebbene possano esistere inversi analitici in forma chiusa (potenzialmente tramite trasformate di Mellin), vengono impiegati algoritmi numerici altamente efficienti (in particolare l'algoritmo di de Hoog et al. implementato come invlap.m in Matlab) per calcolare la soluzione nel dominio del tempo indipendentemente dalla complessità delle espressioni di Laplace.

Contributi e Risultati Chiave
Il lavoro presenta due esempi illustrativi per validare la metodologia:

• Esempio 1: Un problema con una condizione iniziale $w_0(x) = \frac{2a^2\pi \sin(\pi x)}{\sigma + \cos(\pi x)}$ e condizioni al contorno nulle. Gli autori hanno derivato la soluzione nel dominio di Laplace e applicato la trasformata inversa. La soluzione risultante nel dominio del tempo corrispondeva perfettamente alla nota soluzione esatta classica in forma chiusa (Wood, 2014).
• Esempio 2: Un problema con una condizione iniziale $w_0(x) = \sin(\pi x)$. La soluzione esatta classica per questo caso è nota per essere una serie infinita. Il metodo degli autori ha prodotto una soluzione esatta implicita nel dominio di Laplace. Quando invertita numericamente, i risultati si allineavano strettamente con la soluzione della serie infinita troncata ("Soluzione teorica").

Prestazioni Comparative
Lo studio evidenzia due vantaggi primari del metodo proposto rispetto agli approcci esistenti:

1. Efficienza Computazionale: Nell'Esempio 2, la procedura numerica di inversione di Laplace è risultata essere più di 100 volte più veloce in termini di tempo di esecuzione rispetto alla valutazione della soluzione della serie infinita troncata.
2. Accuratezza: Le soluzioni numeriche di inversione di Laplace hanno mostrato errori trascurabili rispetto alle soluzioni teoriche in serie. Al contrario, le soluzioni a differenze finite (calcolate tramite pdepe() di Matlab) hanno mostrato una discrepanza crescente al diminuire del coefficiente di viscosità, in particolare nei primi passi temporali.

Significato e Affermazioni
Gli autori affermano che il loro lavoro fornisce una soluzione operativa in forma chiusa nel dominio di Laplace e una soluzione esatta implicita nel dominio del tempo. Affermano che questo approccio offre un'"alternativa nuova" rispetto alle soluzioni basate su serie infinite o sulla discretizzazione puramente numerica.

Il significato del lavoro è inquadrato come segue:

• Punto di Riferimento: La soluzione esatta derivata funge da benchmark robusto per convalidare e testare altri metodi numerici e semi-analitici.
• Efficienza: Il metodo è computazionalmente superiore alle espansioni in serie, che possono essere onerose in termini di risorse.
• Generalità: La metodologia gestisce soluzioni regolari così come soluzioni generalizzate con discontinuità nelle derivate.

Il lavoro conclude con modestia, affermando che la metodologia di soluzione sarà un riferimento per lavori futuri, inclusi potenziali estensioni alle equazioni di Burgers forzate in dimensioni superiori e confronti con altre tecniche per equazioni di evoluzione non lineari come l'analisi di Painlevé e la tecnica bilineare di Hirota.