Virtual element approximation of eigenvalue problems: is the stabilization of the right hand side necessary?

Titolo

Approssimazione agli Elementi Virtuali (VEM) di Problemi agli Autovalori: È Necessaria la Stabilizzazione del Termine di Destra?

1. Il Problema

Il lavoro affronta l'approssimazione numerica di problemi agli autovalori associati a operatori ellittici auto-aggiunti, utilizzando il metodo degli Elementi Virtuali (VEM).
La formulazione variazionale classica cerca autovalori $\lambda$ e autofunzioni $u$ tali che:
$$a(u, v) = \lambda b(u, v) \quad \forall v \in V$$
dove $a(\cdot, \cdot)$ è la forma bilineare associata all'operatore differenziale (es. Laplaciano) e $b(\cdot, \cdot)$ è la forma bilineare del termine di massa (solitamente il prodotto scalare in $L^2$).

Nella discretizzazione VEM standard, le forme bilineari $a_h$ e $b_h$ sono costruite proiettando le funzioni virtuali su spazi polinomiali. Tuttavia, per garantire la stabilità numerica e la non singolarità delle matrici risultanti, è prassi comune introdurre termini di stabilizzazione parametrici.
Il problema centrale indagato è: è necessario stabilizzare la matrice di massa (il termine di destra, associato a $b$) per ottenere una convergenza corretta degli autovalori?
La letteratura precedente ha evidenziato che una stabilizzazione inadeguata o mal scelta può introdurre "autovalori spurii" che inquinano lo spettro, specialmente per la matrice $B$, spostando gli autovalori spurii verso la parte bassa dello spettro (dove si cercano le prime modalità).

2. Metodologia

Gli autori adottano un approccio teorico basato sulla teoria di Descloux–Nassif–Rappaz per l'analisi di convergenza di operatori non compatti, applicata al contesto VEM.

• Setting Assoluto: Viene definito un operatore soluzione $T$ continuo e un operatore discreto $T_h$. La convergenza degli autovalori è garantita se si dimostrano due proprietà fondamentali:
  1. P1 (Convergenza dell'operatore): $\|(T - T_h)f\|_V \to 0$ per $h \to 0$.
  2. P2 (Approssimabilità): Per ogni funzione soluzione $u$, esiste una funzione discreta $v_h$ tale che $\|u - v_h\|_V \to 0$.
• Costruzione VEM:
  • Viene considerata la discretizzazione dell'equazione di Poisson su mesh poligonali.
  • La forma bilineare $a_h$ (lato sinistro) viene stabilizzata secondo lo standard VEM per garantire la coercività.
  • Punto chiave: La forma bilineare $b_h$ (lato destro, massa) viene definita senza termini di stabilizzazione aggiuntivi, utilizzando solo la proiezione $L^2$ sui polinomi ($\Pi^0_k$).
• Analisi del Nucleo: Poiché $b_h$ non è stabilizzata, può avere un nucleo non banale ($K_{b,h} \neq \{0\}$), ovvero esistono funzioni non nulle $w_h$ tali che $b_h(w_h, v_h) = 0$. Gli autori dimostrano che è possibile restringere l'analisi al sottospazio ortogonale $K_{b,h}^\perp$ rispetto al prodotto scalare indotto da $a_h$, evitando così autovalori infiniti e garantendo la ben posta del problema discreto.

3. Contributi Chiave

Il contributo principale del paper è la dimostrazione teorica e numerica che la stabilizzazione del termine di destra non è necessaria per la convergenza degli autovalori, sotto specifiche condizioni:

1. Dimostrazione Teorica per Ordini Bassi: Viene provato che per spazi VEM di ordine basso ($k=1$ e $k=2$), le proprietà P1 e P2 sono soddisfatte anche senza stabilizzazione di $b_h$.
   • La chiave della dimostrazione risiede nel fatto che, per $k=1, 2$, la proiezione del gradiente $\Pi^\nabla_k$ e la proiezione $L^2$ $\Pi^0_k$ coincidono su certi sottospazi o soddisfano condizioni di momento che permettono di controllare l'errore di consistenza senza bisogno di stabilizzazione aggiuntiva.
2. Analisi del Nucleo: Viene mostrato che anche se la matrice $B$ ha un nucleo non banale (cosa che accade per $k \ge 3$ su certe mesh), la restrizione del problema al sottospazio ortogonale garantisce la convergenza agli autovalori corretti, senza generare autovalori spurii nella finestra di interesse.
3. Risultati Numerici per Ordini Elevati: Sebbene la teoria rigorosa sia limitata a $k=1, 2$, i test numerici estesi mostrano che la convergenza ottimale e l'assenza di autovalori spurii persistono anche per ordini superiori ($k=3, 4$) su diverse sequenze di mesh.

4. Risultati Numerici

Gli autori hanno condotto una vasta serie di test numerici su un dominio quadrato unitario, utilizzando quattro tipi di mesh:

• Triangoli (T) e Quadrati (S): Mesh uniformi standard.
• Voronoi (V): Mesh generiche con poligoni di numero variabile di lati.
• Esagoni (H): Mesh prevalentemente esagonali.
• Dyadic (D): Mesh non standard ottenute suddividendo i quadrati.

Risultati principali:

• Dimensione del Nucleo: Per $k=1$ e $k=2$, il nucleo di $B$ è nullo ($dim(K_{b,h})=0$) su quasi tutte le mesh. Per $k \ge 3$, il nucleo diventa non banale, ma ciò non compromette la soluzione.
• Convergenza: Le tabelle riportate nel paper (da 3 a 25) mostrano che gli errori sugli autovalori convergono con il tasso ottimale atteso ($O(h^{2k})$ per l'errore in norma $L^2$ o $O(h^k)$ in energia, a seconda della metrica considerata) per tutti i tipi di mesh e per tutti gli ordini $k$.
• Assenza di Autovalori Spurii: Non sono stati osservati autovalori spurii che inquinassero lo spettro, nemmeno quando la matrice $B$ non era a rango pieno.
• Oscillazioni per $k=4$: Per ordini molto alti e mesh molto fini, i tassi di convergenza mostrano oscillazioni dovute al raggiungimento della precisione di macchina (machine precision), un fenomeno atteso e non indicativo di instabilità del metodo.

5. Significato e Implicazioni

Questo lavoro ha un impatto significativo sulla pratica computazionale del VEM:

1. Semplificazione dell'Implementazione: Elimina la necessità di calibrare parametri di stabilizzazione per la matrice di massa, che è un processo spesso delicato e fonte di errori (come discusso nell'introduzione, una stabilizzazione eccessiva può spostare gli autovalori spurii nella regione di interesse).
2. Robustezza: Dimostra che il metodo VEM è intrinsecamente stabile per problemi agli autovalori anche con una formulazione "minimalista" per il termine di massa, purché si lavori con ordini bassi o si gestisca correttamente il nucleo tramite proiezione.
3. Validazione Pratica: Fornisce evidenze numeriche solide che estendono la validità del risultato anche a ordini di approssimazione superiori a quelli coperti dalla teoria attuale, suggerendo che la stabilizzazione di $B$ potrebbe essere superflua in un contesto più ampio di quanto ipotizzato.

In conclusione, il paper risponde affermativamente alla domanda del titolo: la stabilizzazione del termine di destra non è necessaria per garantire la convergenza e l'accuratezza nella risoluzione di problemi agli autovalori ellittici con VEM, offrendo una procedura più semplice e robusta per gli utenti.