An Approximate Bounded Cochain Projection

Sintesi Tecnica: Una Proiezione di Cochaine Limitata Approssimata

Enunciato del Problema
Il documento affronta la costruzione di un proiettore, $\Pi^{(k)}$, che mappa un complesso di Hilbert di $k$-forme differenziali a dimensione infinita, $H\Lambda^{(k)}(\Omega)$, su un sottocomplesso polinomiale a pezzi a dimensione finita, $H\Lambda^{(k)}_h(\Omega)$. L'obiettivo è soddisfare tre proprietà specifiche caratteristiche di una "proiezione di cochaine limitata" (bounded cochain projector):

1. Idempotenza: $\Pi^{(k)}\Pi^{(k)} = \Pi^{(k)}$.
2. Limitatezza Uniforme: Il proiettore è limitato nella norma Sobolev ($H\Lambda^{(k)}$) con una costante indipendente dal parametro di mesh $h$.
3. Proprietà di Commutatività: Il proiettore commuta con la derivata esterna ($d$), tale che $d\Pi^{(k)} = \Pi^{(k+1)}d$.

Sebbene le proiezioni di cochaine limitate siano ben stabilite per domini contrattili (dove il complesso de Rham è esatto), il documento evidenzia una difficoltà su domini non contrattili. Su tali domini, il complesso possiede una coomologia non banale (forme armoniche). Le costruzioni standard spesso falliscono nel mantenere la limitatezza uniforme o l'esatta commutatività quando lo spazio a dimensione finita non cattura perfettamente la topologia del dominio (ovvero, quando il divario tra gli spazi armonici a dimensione infinita e a dimensione finita è non nullo).

Metodologia
Gli autori propongono una costruzione basata su una formulazione mista aumentata e un complesso discreto a due livelli.

1. Formulazione Mista Aumentata:
   Inveve di fare affidamento su operatori di smoothing seguiti da interpolazione (che possono mancare di limitatezza uniforme nella norma $H\Lambda$), gli autori definiscono il proiettore $\Pi^{(k)}$ tramite un problema variazionale. Dato $\phi^{(k)} \in H\Lambda^{(k)}(\Omega)$, il proiettore cerca $\lambda^{(i)} \in H\Lambda^{(i)}_h(\Omega)$ e una variabile armonica $h^{(k+1)} \in H^{(k+1)}_h$ soddisfacendo un sistema di equazioni simmetrico.

   • Questo sistema introduce un moltiplicatore di Lagrange, $h^{(k+1)}$, che impone il vincolo che la derivata proiettata appartenga allo spazio delle forme esatte.
   • Gli autori dimostrano la ben-posedness di questo sistema utilizzando condizioni di stabilità, inclusa la coercitività della forma bilineare sulle forme chiuse e le condizioni di inf-sup per la derivata esterna e lo spazio armonico.

2. Proprietà del Sistema Aumentato:

   • Idempotenza: Dimostrata mostrando che se l'input è già nello spazio a dimensione finita, la soluzione del sistema misto è l'input stesso.
   • Limitatezza: Dimostrata decomponendo la soluzione in componenti chiuse e ortogonali, utilizzando la disuguaglianza di Poincaré per limitare le norme $L^2$ e di gradiente (graph norm).
   • Proprietà di Commutatività: Il documento dimostra che il proiettore commuta con $d$ se il dominio è contrattile o se lo spazio a dimensione finita è "topologicamente conforme" (ovvero, il divario tra lo spazio armonico discreto e quello continuo è zero). Se lo spazio non è topologicamente conforme, il moltiplicatore armonico $h^{(k+1)}$ è non nullo, rappresentando l'incoerenza topologica, e la proprietà di commutatività è soddisfatta solo fino all'ordine di grandezza di questo errore.

3. Proiezione di Cochaine Limitata Approssimata (Costruzione a Due Livelli):
   Per affrontare i domini non contrattili dove lo spazio discreto potrebbe non essere topologicamente conforme, gli autori propongono la composizione di due operatori:

   • Passaggio 1: Un proiettore $\Pi^{(k)}$ che mappa lo spazio a dimensione infinita verso uno spazio a dimensione finita arricchito $\tilde{H}\Lambda^{(k)}_h(\Omega)$. Questo spazio arricchito è costruito per essere topologicamente conforme (gap = 0).
   • Passaggio 2: Un proiettore canonico $\Pi^{(k)}_C$ che mappa lo spazio arricchito $\tilde{H}\Lambda^{(k)}_h(\Omega)$ verso lo spazio computazionale target $H\Lambda^{(k)}_h(\Omega)$.
     La composizione $\Pi^{(k)}_C \Pi^{(k)}$ funge da proiettore di cochaine limitata approssimata finale.

Contributi Chiave e Risultati

• Costruzione di un Proiettore Limitato: Il documento fornisce una costruzione rigorosa di un proiettore che è idempotente e uniformemente limitato nella norma Sobolev sia per domini contrattili che non contrattili.
• Gestione dei Domini Non Contrattili: Gli autori mostrano che su domini non contrattili, la proprietà di commutatività è soddisfatta esattamente se lo spazio armonico discreto è topologicamente conforme. Se non lo è, l'errore nella commutatività è quantificato dal moltiplicatore armonico $h^{(k+1)}$, che misura la discrepanza topologica.
• Stabilità Discreta: A differenza di precedenti costruzioni che si affidano a regolarizzazioni continue, questo approccio utilizza una formulazione mista che evita le difficoltà analitiche associate ai destri continui della derivata esterna.
• Commutatività Approssimata: Il documento dimostra che, utilizzando un approccio a due livelli (arricchendo lo spazio per risolvere la topologia prima di proiettare sulla mesh finale), si può ottenere un proiettore che è limitato e commuta con la derivata esterna fino a un'accuratezza arbitraria.

Significatività e Rivendicazioni
Il documento afferma di presentare la costruzione di un proiettore che soddisfa le tre proprietà definenti di una proiezione di cochaine limitata. Nello specifico, asserisce che:

• Su domini contrattili, il proiettore proposto è idempotente, uniformemente limitato e commuta esattamente con la derivata esterna.
• Su domini non contrattili, il proiettore rimane idempotente e uniformemente limitato. La proprietà di commutatività è soddisfatta "fino a un'accuratezza arbitraria", il che significa che l'errore può essere reso arbitrariamente piccolo raffinando la risoluzione topologica dello spazio intermedio (ad esempio, tramite $p$-raffinamento fino a quando il moltiplicatore armonico svanisce alla precisione di macchina).

Il lavoro si inserisce nel framework della Finite Element Exterior Calculus (FEEC), perfezionando le idee esistenti riguardanti le costruzioni locali e le condizioni al contorno. Si distingue offrendo una costruzione che è "completamente locale", nel senso di evitare destri continui, affidandosi invece a problemi discreti posti su un raffinamento ausiliario per gestire i vincoli topologici. Gli autori posizionano questo lavoro come un metodo per garantire stabilità e consistenza nei metodi numerici per forme differenziali su geometrie complesse.