$hp$-a posteriori error estimates for hybrid high-order methods applied to biharmonic problems

Problematica
Questo articolo affronta la stima dell'errore a posteriori per i metodi Hybrid High-Order (HHO) applicati al problema biarmonico ($\Delta^2 u = f$) con condizioni al contorno incastrate ($u = \partial_n u = 0$ su $\partial\Omega$). Il problema è posto su domini politopali Lipschitz in due e tre dimensioni con topologia generale (potenzialmente contenenti fori). L'obiettivo è derivare stimatori dell'errore $hp$-basati sui residui che siano robusti rispetto alla dimensione della mesh ($h$), al grado polinomiale ($p$ o $k$) e alla topologia del dominio.

Metodologia
Gli autori sviluppano un framework per limitare l'errore totale decomponendolo in componenti conformi e non conformi.

1. Decomposizione dell'errore: L'errore è diviso in una parte conforme (misurata nella norma energetica di una funzione ricostruita in $H^2_0$) e una parte non conforme (la distanza tra la soluzione discreta e lo spazio conforme).
2. Controllo dell'errore non conforme:
   • Per limitare l'errore non conforme, gli autori costruiscono una ricostruzione conforme $H^2_0$ della soluzione discreta utilizzando una partizione dell'unità $C^1$.
   • Questa partizione è costruita utilizzando spazi di elementi finiti conformi $C^1$ (specificamente, la costruzione tratta da [40] basata sugli Alfeld splitting della mesh).
   • Fondamentalmente, gli autori impiegano decomposizioni di Helmholtz locali su patch di vertici (stelle). A differenza di approcci precedenti che si basavano su decomposizioni di Helmholtz globali (che introducono costanti dipendenti dal numero di fori nel dominio), questo approccio locale garantisce che la costante di stabilità dipenda solo dalla regolarità della forma della mesh (shape-regularity).
3. Controllo dell'errore conforme:
   • Vengono derivati due distinti limiti superiori per l'errore conforme, ciascuno dei quali utilizza un operatore di interpolazione specifico:
     • Setting 1: Utilizza l'operatore di interpolazione canonico degli elementi finiti ibridi ($C^{k+2}_T$). Questo approccio si basa sulla proprietà per cui la composizione di questo operatore con la ricostruzione HHO produce una proiezione $H^2$-ellittica. Il limite risultante coinvolge solo il termine di stabilizzazione e l'oscillazione dei dati, ma si basa su un'assunzione numerica riguardante l'errore di interpolazione (Assunzione 3.4).
     • Setting 2: Utilizza l'operatore di interpolazione $hp$ modificato di Babuška–Suri ($I^{k+2}_{mBS}$). Questo approccio non richiede assunzioni aggiuntive ma introduce termini di indicatore di errore supplementari, specificamente residui di volume e salti di flusso normali/tangenziali.
4. Stimatori del residuo: Il documento definisce indicatori di errore locali comprendenti termini di stabilizzazione, residui di volume, salti di flusso normali e salti tangenziali, scalati da appropriati fattori $hp$.

Contributi Chiave

• Primo stimatore $hp$ per HHO su problemi biarmonici: Questo lavoro fornisce il primo stimatore dell'errore a posteriori basato sui residui per metodi HHO applicati al problema biarmonico su domini generici in 2D e 3D.
• Decomposizione di Helmholtz locale con robustezza topologica: Il documento stabilisce una nuova stima di stabilità per le decomposizioni di Helmholtz locali su stelle di vertici. Dimostra che la costante di stabilità dipende esclusivamente dalla regolarità della forma della mesh, evitando così il deterioramento delle costanti associato al numero di fori nel dominio, un problema noto negli approcci di decomposizione globale.
• Limiti superiori duali: La derivazione di due distinti limiti superiori per l'errore conforme offre flessibilità. Un limite è "pulito" (basato solo su stabilizzazione e oscillazione) ma condizionato a un'assunzione verificata numericamente; l'altro è incondizionato ma include termini di residuo più complessi.
• Estensione a topologie generali: La metodologia è esplicitamente progettata per gestire domini con topologia arbitraria (molteplici fori), una caratteristica spesso difficile da affrontare nella stima $hp$ per problemi del quarto ordine.

Risultati

• Limiti teorici: Il teorema principale (Teorema 5.6) stabilisce che l'errore totale (norma dell'Hessiano rotto più stabilizzazione) è limitato dalla somma degli indicatori di errore proposti. Si dimostra che i limiti sono $h$-ottimali e $p$-subottimali (di al massimo un ordine, o potenzialmente $1/2$ ordine sotto l'assunzione verificata).
• Verifica numerica:
  • Verifica dell'assunzione: Esperimenti numerici in 2D e 3D confermano la validità dell'Assunzione 3.4 riguardante l'errore di interpolazione dell'operatore ibrido canonico.
  • Tassi di convergenza: Test su soluzioni regolari in 3D dimostrano che sia l'errore energetico che lo stimatore convergono al tasso ottimale $O(\text{DoFs}^{-(k+1)/d})$.
  • Indici di effettività: Gli stimatori forniscono indici di effettività moderati (che variano approssimativamente tra 1.2 e 2.0) attraverso vari gradi polinomiali e raffinamenti della mesh.
  • Robustezza topologica: Esperimenti su domini 2D con un numero crescente di fori (0, 1, 4, 8) mostrano che gli indici di effettività rimangono stabili, confermando la tesi teorica secondo cui lo stimatore è robusto rispetto alla topologia del dominio.
  • Performance adattiva: Algoritmi $h$-adattivi guidati dallo stimatore recuperano con successo i tassi di convergenza ottimali per soluzioni singolari sia in 2D (dominio a forma di L) che in 3D (singolarità puntuale), superando il raffinamento uniforme.

Significatività
Gli autori pongono questo lavoro come il primo a impiegare una partizione dell'unità $C^1$ combinata con decomposizioni di Helmholtz locali su stelle di vertici per limitare gli errori non conformi nell'ambito $hp$ per problemi biarmonici. La significatività primaria risiede nella capacità di controllare l'errore non conforme senza incorrere in costanti di stabilità che si deteriorano con la complessità topologica del dominio (numero di fori). Inoltre, gli autori osservano che, sebbene l'approccio sia sviluppato per i metodi HHO, la tecnica per limitare gli errori non conformi tramite decomposizioni di Helmholtz locali è potenzialmente applicabile ad altri metodi non conformi, come i metodi Interior Penalty Discontinuous Galerkin (IPDG). Il documento non pretende di fornire limiti inferiori dell'errore, notando che tali limiti seguirebbero adattamenti di argomentazioni esistenti ma presenterebbero probabilmente una subottimalità $p$ pessimistica a causa delle limitazioni delle attuali tecniche $hp$ per la costruzione di specifiche funzioni bubble ed operatori di estensione.