Error estimates for completely discrete FEM in energy-type and weaker norms

Sintesi Tecnica di "Error estimates for completely discrete FEM in energy-type and weaker norms"

Enunciato del Problema
Il saggio affronta l'analisi numerica dei Metodi degli Elementi Finiti (FEM) per problemi ai valori al contorno che coinvolgono equazioni differenziali ellittiche lineari del secondo ordine con una struttura di diffusione-convezione-reazione. I problemi sono definiti su domini poliedrici Lipschitziani limitati in $\mathbb{R}^d$ ($d \in \{2, 3\}$) e presentano condizioni al contorno miste di primo, secondo e terzo tipo. Un focus centrale sono le formulazioni "completamente discrete", dove gli spazi discreti degli elementi finiti non sono necessariamente conformi (ovvero non sono immersi nello spazio della soluzione continua della formulazione debole) e i problemi discreti non soddisfano necessariamente le standard proprietà di consistenza. Gli autori mirano a fornire un quadro astratto unificato per derivare stime di errore sia in norme di tipo energia che in norme più deboli (come $L^2$) senza assumere a priori la conformità o la consistenza.

Metodologia
Gli autori sviluppano un quadro variazionale generale basato su due spazi di Banach reali, $U$ e $V$, e i loro corrispondenti discreti $U_h$ e $V_h$. La metodologia procede attraverso i seguenti passaggi:

1. Analisi dell'Errore Astratta:

   • Stabilità e Consistenza: Il saggio generalizza il secondo lemma di Strang per gestire discretizzazioni non conformi e non consistenti. Definisce la consistenza su un sottospazio $W \subset U$ richiedendo che la forma bilineare discreta $a_h$ possa essere estesa in modo che la soluzione continua $u$ soddisfi esattamente l'equazione discreta (fino alla forma lineare discreta $\ell_h$).
   • Stime della Norma Energia: Il Teorema 3 fornisce una stima dell'errore nella norma energia discreta $\|u - u_h\|_{U_h}$. Questa stima dipende dall'errore di approssimazione (migliore approssimazione nello spazio discreto), dalla costante di stabilità (condizione inf-sup) e da un termine di errore di consistenza che misura lo scostamento della soluzione continua dalla formulazione discreta.
   • Stime in Norme Più Deboli (Aubin-Nitsche Generalizzato): Il saggio estende il classico argomento di dualità di Aubin-Nitsche al contesto non conforme/non consistente. Il Teorema 5 deriva una stima dell'errore in una norma più debole $\|u - u_h\|_{Z'}$ (dove $Z$ è uno spazio di Banach riflessivo). Questo limite decompone l'errore in quattro termini che coinvolgono:
     • Il prodotto degli errori di approssimazione primale e duale.
     • Gli errori di consistenza dei problemi primale e duale.
     • Gli errori di approssimazione derivanti dalla differenza tra le forme bilineari/lineari continue e discrete.
   • L'analisi si basa sulla solvibilità e sulla regolarità di un problema aggiunto e del suo corrispondente discreto.

2. Applicazione a Discretizzazioni Specifiche:
   La teoria astratta è applicata a due concrete discretizzazioni stabilizzate per l'equazione di diffusione-convezione-reazione:

   • Esempio I: Discretizzazione Crouzeix-Raviart (CR1): Un metodo non conforme che utilizza polinomi lineari a tratti che sono continui nei baricentri delle facce degli elementi ma discontinui attraverso le facce stesse. Il metodo include termini di stabilizzazione upwind per la parte convettiva. Gli autori dimostrano la coercitività e derivano stime di errore ottimali nella norma di tipo $H^1$ e nella norma $L^2$.
   • Esempio II: Metodi Discontinui a Penalità Interna (IPG): Una famiglia di metodi (che include SIPG, NIPG e IIPG) che utilizzano polinomi a tratti discontinui di grado $k$. La formulazione include termini di penalità per i salti nella soluzione e nei flussi, oltre alla stabilizzazione upwind. Gli autori stabiliscono la stabilità (coercitività) sotto specifiche condizioni sul parametro di penalità e sulla regolarità della mesh, e dimostrano tassi di convergenza.

Contributi Chiave e Risultati

• Quadro Unificato: Il saggio presenta un approccio unificato che non richiede che gli spazi discreti siano sottospazi degli spazi continui (non conformità) né che la formulazione discreta sia consistente con quella continua. Ciò consente l'analisi di una classe più ampia di metodi stabilizzati.
• Tecnica di Aubin-Nitsche Generalizzata: Un contributo significativo è l'estensione dell'argomento di dualità ai contesti completamente discreti. Il Teorema 5 fornisce un limite rigoroso per gli errori in norme più deboli, tenendo esplicitamente conto degli errori di consistenza sia nel problema primale che in quello aggiunto.
• Stime di Errore Ottimali:
  • Per il metodo Crouzeix-Raviart, il saggio dimostra una convergenza del primo ordine ottimale nella norma energia ($O(h)$) e una convergenza del secondo ordine nella norma $L^2$ ($O(h^2)$), assumendo una sufficiente regolarità della soluzione e dei coefficienti.
  • Per i metodi IPG, sono stabiliti tassi di convergenza ottimali. Nella norma energia, l'errore è $O(h^m)$ per il grado polinomiale $k$ (con $m \le k$). Nelle norme più deboli, il saggio dimostra che per la versione simmetrica IPG (SIPG, $\theta = -1$), il problema aggiunto rimane consistente, permettendo stime di errore ottimali in norme negative ($O(h^{k+m+1})$). Per le varianti non simmetriche ($\theta \neq -1$), si dimostra che il problema aggiunto è non consistente, richiedendo una stima aggiuntiva del termine di errore di consistenza.
• Analisi della Stabilità: Gli autori forniscono dimostrazioni dettagliate della coercitività delle forme bilineari discrete, identificando le condizioni necessarie sui parametri di penalità e sulla geometria della mesh (regolarità di forma e di contatto) per garantire la stabilità.

Significatività e Rivendicazioni
Gli autori affermano che, sebbene alcuni aspetti dell'esposizione (come la loro versione del secondo lemma di Strang e le varianti dell'argomento di dualità) non siano del tutto nuovi, il valore del lavoro risiede nel quadro astratto unificato. Questo quadro chiarisce e quantifica l'interazione tra stabilità, errori di approssimazione ed errori di consistenza senza fare affidamento sulle restrittive assunzioni di conformità o consistenza.

Il saggio sostiene che tale approccio:

1. Permette l'analisi di metodi a elementi finiti più generali, inclusi quelli che sono intrinsecamente non conformi o non consistenti.
2. Fornisce un modo sistematico per derivare stime di errore ottimali in norme più deboli generalizzando la tecnica di Aubin-Nitsche a questi contesti più ampi.
3. Offre linee guida per impostare e analizzare nuovi schemi stabili e convergenti per problemi di diffusione-convezione-reazione con condizioni al contorno miste.

Gli autori notano con modestia che i risultati teorici sono applicati a due esempi specifici (CR1 e IPG) per dimostrare l'utilità del quadro, ma i risultati astratti sono destinati ad essere applicabili a una gamma più ampia di modelli e tecniche numeriche. Il lavoro costruisce ed estende la letteratura precedente (ad esempio, [DPD18], [DPD21], [DPE12]) evitando l'inclusione di operatori di interpolazione nella definizione di consistenza e gestendo le condizioni al contorno miste e i termini di convezione in modo più generale.