Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshes

Sintesi Tecnica: Metodo DG di Trefftz Incorporato per Problemi di Reazione–Diffusione su Mesh Anisotropi

Enunciato del Problema
Il documento affronta la soluzione numerica di problemi di reazione–diffusione singolarmente perturbati della forma $-\varepsilon^2 \Delta u + u = f$ in un dominio limitato $\Omega \subset \mathbb{R}^2$, con condizioni al contorno di Dirichlet. Il focus è sul regime in cui $0 < \varepsilon \ll 1$, uno scenario in cui le soluzioni presentano strati limite (boundary layers) netti e caratteristiche fortemente anisotrope. Le mesh isotrope standard sono spesso inefficienti per tali problemi, rendendo necessaria l'impiego di un raffinamento di mesh anisotropo (elementi quadrilateri allungati in direzioni specifiche) per risolvere le caratteristiche localizzate con meno gradi di libertà. Sebbene i metodi Discontinuous Galerkin (DG) siano flessibili per discretizzazioni ad alto ordine su mesh generali, le formulazioni DG a prodotto tensoriale standard risultano in sistemi algebrici globali di grandi dimensioni. I metodi di Trefftz, che utilizzano spazi di approssimazione locali contenuti nel nucleo dell'operatore differenziale, offrono una riduzione dei gradi di libertà (DoF); tuttavia, la costruzione di spazi di Trefftz idonei per PDE ellittiche generali è non banale, particolarmente su elementi quadrilateri anisotropi e curvi.

Metodologia
Gli autori propongono un metodo Discontinuous Galerkin (DG) di Trefftz Incorporato su misura per mesh quadrilateri anisotropi. Il metodo è costruito attraverso le seguenti fasi:

1. Mesh e Spazi: Il dominio è partizionato in elementi quadrilateri curvilinei. Ogni elemento $K$ è mappato da un quadrato unitario tramite una trasformazione affine (per gestire l'anisotropia) e una diffeomorfismo $C^1$ (per gestire la curvatura). Lo spazio discreto consiste in polinomi a prodotto tensoriale nelle coordinate ausiliarie.
2. Costruzione di Trefftz Incorporata: Invece di costruire esplicitamente lo spazio di Trefftz globale, il metodo impiega una strategia di embedding. Uno spazio di test locale $L_h(K)$ è definito su ogni elemento. Questo spazio consiste in polinomi a prodotto tensoriale che si annullano sui due bordi più lunghi dell'elemento di pullback.
3. Operatore Locale e Vincoli: Un operatore locale di reazione–diffusione $A_K$ è definito. Lo spazio di Trefftz incorporato $T_h(K)$ è definito come il sottospazio dello spazio DG standard dove il residuo della PDE, testato contro funzioni in $L_h(K)$, svanisce elemento per elemento: $(A_K v_h, q_h)_K = 0$ per ogni $q_h \in L_h(K)$.
4. Formulazione Globale: Il metodo globale cerca una soluzione $u_h$ soddisfacendo due condizioni:
   • Vincoli Locali: $(A_K u_h, v_L)_K = \ell_K(v_L)$ per ogni funzione di test locale $v_L \in L_h(K)$.
   • Consistenza Globale: La forma bilineare DG a penalità interiore standard $a_h(u_h, v_T) = \ell_h(v_T)$ è soddisfatta per tutte le funzioni di test globali $v_T$ nello spazio di Trefftz.
5. Implementazione: In pratica, lo spazio di Trefftz non viene assemblato esplicitamente. Invece, il metodo utilizza la Decomposizione ai Valori Singolari (SVD) delle matrici dell'operatore locale per costruire un embedding, mantenendo solo i DoF globali associati allo spazio di Trefftz.

Contributi Chiave
Il documento presenta diversi contributi tecnici specifici:

• Design degli Spazi di Test Locali: La principale novità è il design di uno spazio di test locale compatibile con basi polinomiali a prodotto tensoriale. A differenza dei lavori precedenti che utilizzano basi simpliciali, questo approccio consente una costruzione sistematica su mesh quadrilateri ed esaedri.
• Stabilità su Elementi Anisotropi: Gli autori dimostrano stabilità e quasi-ottimalità su elementi quadrilateri anisotropi, potenzialmente curvi. Un aspetto crucialo è la scelta specifica dello spazio di test (che si annulla sui bordi più lunghi), che permette la derivazione di disuguaglianze inverse anisotrope e disuguaglianze di Poincaré–Steklov essenziali per la stabilità.
• Analisi dell'Errore: Il documento deriva stime di errore a priori anisotrope nella norma dell'energia. Queste stime catturano il comportamento della soluzione nelle reali direzioni di anisotropia, fornendo limiti più stretti rispetto alle stime isotrope.
• Riduzione della Dimensione: Il metodo riduce la dimensione dello spazio di approssimazione locale da $(p+1)^2$ (prodotto tensoriale standard) a $2p+2$, riducendo significativamente il numero di gradi di libertà globalmente accoppiati rispetto ai metodi DG e HDG standard, particolarmente in 3D.

Risultati
Esperimenti numerici validano i risultati teorici:

• Convergenza $h$: Su domini quadrati e circolari con strati limite, il metodo di Trefftz incorporato raggiunge tassi di convergenza ottimali nella norma dell'energia. Richiede significativamente meno gradi di libertà per raggiungere la stessa accuratezza dei metodi DG standard.
• Domini Curvi: Il metodo opera in modo robusto su mesh curve, inclusi casi con elementi misti triangolari e quadrilateri (dove gli elementi triangolari sono trattati tramite tecniche di embedding standard).
• Convergenza $hp$: Su una mesh a nove elementi adattata agli strati limite, il metodo mostra una convergenza esponenziale. Tuttavia, gli autori notano che per gradi polinomiali elevati ($p \gtrsim 5$) e grandi dimensioni degli strati limite, l'embedding standard può diventare instabile. Ciò viene mitigato modificando lo spazio di test negli elementi d'angolo per rilassare la condizione di Trefftz, ripristinando la stabilità e la convergenza esponenziale.
• Diffusione Anisotropa: Il metodo è testato su un problema di pura diffusione anisotropa (senza termine di reazione), che si colloca al di fuori del rigoroso framework teorico. I risultati mostrano che il metodo di Trefftz incorporato produce costantemente errori inferiori rispetto al DG standard sia sotto strategie di raffinamento $h$ che $p$.

Significatività e Rivendicazioni
Il documento rivendica di fornire un framework generale per i metodi DG di tipo Trefftz che rimane efficace anche quando il nucleo dell'operatore della PDE non è esplicitamente noto o è banale. Estendendo con successo l'approccio di Trefftz incorporato alle mesh quadrilateri anisotrope, il lavoro apre la strada a discretizzazioni ad alto ordine su mesh esaedriche per problemi singolarmente perturbati. Gli autori sottolineano che il metodo offre un'alternativa competitiva all'HDG, producendo i sistemi globali più piccoli in tre dimensioni pur mantenendo le proprietà di approssimazione delle discretizzazioni ad alto ordine. L'analisi è modesta riguardo alla robustezza $p$, riconoscendo che la stabilità per gradi polinomiali molto elevati su mesh grossolane con grandi strati limite può richiedere modifiche specifiche, come dimostrato dagli esperimenti numerici.