An Approximate Bounded Cochain Projection

Technisches Resümee: Eine approximative gebundene Cochain-Projektion

Problemstellung
Die Arbeit befasst sich mit der Konstruktion eines Projektors $\Pi^{(k)}$, der einen unendlichdimensionalen Hilbert-Komplex von Differentialformen $k$-ten Grades, $H\Lambda^{(k)}(\Omega)$, auf einen endlichdimensionalen, stückweise polynomielle Unterkomplex $H\Lambda^{(k)}_h(\Omega)$ abbildet. Das Ziel ist es, drei spezifische Eigenschaften zu erfüllen, die charakteristisch für eine „gebundene Cochain-Projektion“ sind:

1. Idempotenz: $\Pi^{(k)}\Pi^{(k)} = \Pi^{(k)}$.
2. Gleichmäßige Beschränktheit: Der Projektor ist in der Sobolev-Norm ($H\Lambda^{(k)}$) beschränkt, wobei die Konstante unabhängig vom Netzparameter $h$ ist.
3. Kommutativitätseigenschaft: Der Projektor kommutiert mit der äußeren Ableitung ($d$), sodass $d\Pi^{(k)} = \Pi^{(k+1)}d$ gilt.

Während gebundene Cochain-Projektionen für kontraktile Domänen (wo der de Rham-Komplex exakt ist) gut etabliert sind, hebt die Arbeit eine Schwierigkeit auf nicht-kontraktilen Domänen hervor. Auf solchen Domänen besitzt der Komplex eine nicht-triviale Kohomologie (harmonische Formen). Standardkonstruktionen scheitern oft daran, die gleichmäßige Beschränktheit oder die exakte Kommutativität aufrechtzuerhalten, wenn der endlichdimensionale Raum die Topologie der Domäne nicht perfekt erfasst (d. h. wenn die Lücke zwischen den unendlichdimensionalen und den endlichdimensionalen harmonischen Räumen ungleich Null ist).

Methodik
Die Autoren schlagen eine Konstruktion vor, die auf einer augmentierten gemischten Formulierung und einem zweistufigen diskreten Komplex basiert.

1. Augmentierte gemischte Formulierung:
   Anstatt sich auf Glättungsoperatoren zu verlassen, die anschließend interpolieren (was eine gleichmäßige Beschränktheit in der $H\Lambda$-Norm vermissen lassen kann), definieren die Autoren den Projektor $\Pi^{(k)}$ über ein Variationsproblem. Gegeben sei $\phi^{(k)} \in H\Lambda^{(k)}(\Omega)$, sucht der Projektor $\lambda^{(i)} \in H\Lambda^{(i)}_h(\Omega)$ und eine harmonische Variable $h^{(k+1)} \in H^{(k+1)}_h$, die ein symmetrisches Gleichungssystem erfüllen.

• Dieses System führt einen Lagrange-Multiplikator, $h^{(k+1)}$, ein, der die Bedingung erzwingt, dass die projektierte Ableitung im Raum der exakten Formen liegt.
• Die Autoren beweisen die Wohldefiniertheit dieses Systems unter Verwendung von Stabilitätsbedingungen, einschließlich der Koerzitivität der bilinearen Form auf geschlossenen Formen sowie Inf-Sup-Bedingungen für die äußere Ableitung und den harmonischen Raum.

2. Eigenschaften des augmentierten Systems:

• Idempotenz: Bewiesen durch den Nachweis, dass, falls das Eingangsobjekt bereits im endlichdimensionalen Raum liegt, die Lösung des gemischten Systems das Eingangsobjekt selbst ist.
• Beschränktheit: Bewiesen durch die Zerlegung der Lösung in geschlossene und orthogonale Komponenten unter Nutzung der Poincaré-Ungleichung, um die $L^2$- und Graph-Normen zu beschränken.
• Kommutativitätseigenschaft: Die Arbeit zeigt, dass der Projektor mit $d$ kommutiert, wenn die Domäne kontraktil ist oder wenn der endlichdimensionale Raum „topologisch konform“ ist (d. h. die Lücke zwischen dem diskreten und dem kontinuierlichen harmonischen Raum Null ist). Wenn der Raum nicht topologisch konform ist, ist der harmonische Multiplikator $h^{(k+1)}$ ungleich Null, was die topologische Inkonsistenz repräsentiert, und die Kommutativitätseigenschaft gilt nur bis zur Größenordnung dieses Fehlers.

3. Approximative gebundene Cochain-Projektion (Zweistufige Konstruktion):
   Um das Problem auf nicht-kontraktilen Domänen anzugehen, bei denen der diskrete Raum möglicherweise nicht topologisch konform ist, schlagen die Autoren eine Komposition zweier Operatoren vor:

• Schritt 1: Ein Projektor $\Pi^{(k)}$, der den unendlichdimensionalen Raum auf einen angereicherten endlichdimensionalen Raum $\tilde{H}\Lambda^{(k)}_h(\Omega)$ abbildet. Dieser angereicherte Raum ist so konstruiert, dass er topologisch konform ist (Lücke = 0).
• Schritt 2: Ein kanonischer Projektor $\Pi^{(k)}_C$, der den angereicherten Raum $\tilde{H}\Lambda^{(k)}_h(\Omega)$ auf den Ziel-Rechenraum $H\Lambda^{(k)}_h(\Omega)$ abbildet.
  Die Komposition $\Pi^{(k)}_C \Pi^{(k)}$ dient als die finale approximative gebundene Cochain-Projektion.

Zentrale Beiträge und Ergebnisse

• Konstruktion eines beschränkten Projektors: Die Arbeit liefert eine rigorose Konstruktion eines Projektors, der idempotent und in der Sobolev-Norm gleichmäßig beschränkt ist, sowohl für kontraktile als auch für nicht-kontratile Domänen.
• Umgang mit nicht-kontraktilen Domänen: Die Autoren zeigen, dass auf nicht-kontraktilen Domänen die Kommutativitätseigenschaft exakt erfüllt ist, falls der diskrete harmonische Raum topologisch konform ist. Falls nicht, wird der Fehler in der Kommutativität durch den harmonischen Multiplikator $h^{(k+1)}$ quantifiziert, welcher die topologische Diskrepanz misst.
• Diskrete Stabilität: Im Gegensatz zu früheren Konstruktionen, die auf kontinuierlichen Regularisierungen beruhen, verwendet dieser Ansatz eine gemischte Formulierung, die analytische Schwierigkeiten vermeidet, die mit kontinuierlichen rechten Inversen der äußeren Ableitung verbunden sind.
• Approximative Kommutativität: Die Arbeit demonstriert, dass durch die Verwendung eines zweistufigen Ansatzes (die Anreicherung des Raums zur Auflösung der Topologie vor der Projektion auf das endgültige Netz) ein Projektor erreicht werden kann, der beschränkt ist und mit der äußeren Ableitung bis zu einer beliebigen Genauigkeit kommutiert.

Bedeutung und Behauptungen
Die Arbeit beansprucht, die Konstruktion eines Projektors zu präsentieren, der die drei definierenden Eigenschaften einer gebundenen Cochain-Projektion erfüllt. Insbesondere stellt sie fest, dass:

• Auf kontraktilen Domänen der vorgeschlagene Projektor idempotent ist, gleichmäßig beschränkt ist und exakt mit der äußeren Ableitung kommutiert.
• Auf nicht-kontraktilen Domänen der Projektor idempotent und gleichmäßig beschränkt bleibt. Die Kommutativitätseigenschaft wird „bis zu einer beliebigen Genauigkeit“ erfüllt, was bedeutet, dass der Fehler durch Verfeinerung der topologischen Auflösung des Zwischenraums (z. B. mittels $p$-Verfeinerung, bis der harmonische Multiplikator bis zur Maschinengenauigkeit verschwindet) beliebig klein gemacht werden kann.

Die Arbeit baut auf dem Framework der Finite Element Exterior Calculus (FEEC) auf und verfeinert bestehende Ideen bezüglich lokaler Konstruktionen und Randbedingungen. Sie unterscheidet sich dadurch, dass sie eine Konstruktion bietet, die „vollständig lokal“ ist, indem sie kontinuierliche rechte Inverse vermeidet und stattdessen auf diskreten Problemen auf einem Hilfsnetz basiert, um topologische Zwänge zu handhaben. Die Autoren positionieren dies als Methode, um Stabilität und Konsistenz in numerischen Verfahren für Differentialformen auf komplexen Geometrien zu gewährleisten.