Error estimates for the interpolation and approximation of gradients and vector fields on protected Delaunay meshes in $\mathbb{R}^d$

Technische Zusammenfassung: Fehlerschätzungen für Interpolation und Approximation auf geschützten Delaunay-Gittern

Problemstellung
Die Arbeit adressiert eine kritische Lücke im mathematischen Verständnis der Interpolation und Approximation auf simplizialen Gittern in Dimensionen $d > 2$. Während Standard-Delaunay-Gitter in zwei Dimensionen ($\mathbb{R}^2$) weit verbreitet und gut verstanden sind, ist ihr Verhalten in höheren Dimensionen weniger günstig. Speziell enthalten Standard-Delaunay-Triangulationen in $\mathbb{R}^d$ ($d>2$) häufig „Sliver“-Elemente (Sliver-Simplizes) – Simplizes mit einem sehr kleinen Verhältnis von minimaler Höhe zu maximaler Kantenlänge. Diese Slivers können eine Dicke aufweisen, die beliebig nah an Null liegt, was die Qualität der Gradienten- und Vektorfeld-Interpolation und -Approximation verschlechtert.

Die bestehende Literatur liefert robuste Fehlerschätzungen für stückweise lineare Interpolation auf Standard-Delaunay-Gittern (z. B. Minimierung von Rauheit oder spezifischer Funktionale), jedoch sind Ergebnisse für die hochgradige, stückweise polynomielle Interpolation spärlich, beruhen oft auf impliziten gitterabhängigen Konstanten oder sind nicht explizit auf das Delaunay-Setting zugeschnitten. Zudem minimieren Standard-Delaunay-Gitter zwar bestimmte Funktionale, garantieren jedoch in höheren Dimensionen keine untere Schranke für die Elementdicke, was zu einer Instabilität der Fehlerschätzungen führen kann.

Methodik
Die Autoren entwickeln explizite Fehlerschätzungen für die Interpolation und Approximation von Gradienten und Vektorfeldern auf geschützten Delaunay-Gittern. Ein geschütztes Delaunay-Gitter erfüllt ein modifiziertes leere-Kugel-Kriterium, bei dem die Umkugel jedes Simplizes durch eine sphärische Pufferregion mit dem Radius $r \ge \delta$ (wobei $\delta$ der „Schutz“-Parameter ist) ergänzt wird. Diese Protektion verhindert die Entstehung pathologischer Slivers.

Die Methodik verläuft wie folgt:

1. Rauheit und Norm-Äquivalenz: Die Autoren erweitern das Konzept der „Rauheit“ (ursprünglich von Rippa für $\mathbb{R}^2$ definiert) auf $\mathbb{R}^d$. Sie definieren ein „Rauheitsfunktional“ $\Psi(\nabla v)$ basierend auf einer gewichteten Summe der quadrierten Richtungsableitungen entlang der Kanten der Simplizes. Sie beweisen, dass dieses Funktional äquivalent zur Standard-$L^2$-Norm des Gradienten ist, wobei die Äquivalenzkonstanten von der Gittergeometrie (speziell der minimalen Dicke $C_\Xi$ und einem maximalen Kantenlängen-Funktional $C_\Upsilon$) abhängen.
2. Ableitung des Interpolationsfehlers: Unter Verwendung klassischer Interpolationstheorie (Lagrange-Interpolation) und der etablierten Äquivalenz leiten die Autoren obere Schranken für den Gradienten-Interpolationsfehler ab. Sie setzen diesen Fehler in Beziehung zur Hesse-Matrix der Funktion, der Gitterweite $h$ und spezifischen gitterabhängigen Konstanten.
3. Approximationsfehler für elliptische Probleme: Für elliptische Probleme (speziell $-\nabla^2 u = f$) nutzen die Autoren das Cea-Lemma, um den Gradienten-Approximationsfehler ($\nabla(u - u_h)$) zu begrenzen. Sie verknüpfen diesen Fehler mit dem Gitterregularitätsparameter $\sigma(K)$ und dem maximalen Radius der minimalen Einbettungskugel $R_{\max}$.
4. Generalisierung auf Vektorfelder: Die Analyse wird von skalaren Gradienten auf allgemeine Vektorfelder generalisiert, wobei gezeigt wird, dass ähnliche Fehlergrenzen auch hier gelten.

Wesentliche Beiträge und Ergebnisse
Die Arbeit liefert explizite, geometrisch abhängige Fehlerschätzungen für die hochgradige polynomielle Interpolation und Approximation. Die Hauptergebnisse sind:

• Explizite Fehlergrenzen: Die Autoren leiten Ungleichungen der Form ab:
  $$\|\nabla v - I_h(\nabla v)\|_{L^2(\Omega)} \le C \tilde{C}(\mathcal{T}_h) h^r \|\nabla(\nabla v)\|_{L^{2\varrho}(\Omega)}$$
  wobei $\tilde{C}(\mathcal{T}_h)$ eine explizit durch Gittereigenschaften definierte Konstante ist (Rajans Funktional $\Theta$, minimale Dicke $C_\Xi$, Durchmesserverhältnis $C_\Delta$ und maximaler Radius der minimalen Einbettungskugel $R_{\max}$).
• Rolle der Protektion: Die Analyse zeigt, dass die Qualität der Interpolation stark von der minimalen Dicke ($C_\Xi$) der Simplizes abhängt.
  • In Standard-Delaunay-Gittern ($d>2$) kann $C_\Xi$ gegen Null gehen, wodurch die Fehlerkonstanten explodieren.
  • In geschützten Delaunay-Gittern ist $C_\Xi$ durch den Schutzparameter $\delta$ nach unten beschränkt (speziell $C_\Xi \ge \delta^2 / (8d\varepsilon^2)$).
  • Folglich ist der Gitterregularitätsparameter $C_\sigma$ (der umgekehrt proportional zur Dicke ist) auf geschützten Gittern nach oben beschränkt.
• Optimierung der Gitterparameter: Die Arbeit bestätigt, dass Standard-Delaunay-Gitter das Rajan-Funktional $\Theta$ (gewichtete Summe der quadrierten Kantenlängen) und den maximalen Radius der minimalen Einbettungskugel $R_{\max}$ minimieren. Geschützte Delaunay-Gitter sind jedoch essenziell, da sie über den Schutzparameter $\delta$ eine untere Schranke für die minimale Dicke $C_\Xi$ garantieren, welche Standard-Gitter in höheren Dimensionen nicht bieten können. Diese Kombination stellt sicher, dass die kritischen Konstanten in den Fehlerschätzungen kontrolliert bleiben.

Bedeutung und Ansprüche
Die Arbeit behauptet, dass ihr primärer Beitrag in der Bereitstellung expliziter, interpretierbarer Fehlerschätzungen für die hochgradige Interpolation und Approximation auf geschützten Delaunay-Gittern in beliebigen Dimensionen besteht.

• Schließung der Lücke: Im Gegensatz zu früheren Ergebnissen, die sich auf stückweise lineare Interpolation konzentrierten oder implizite Konstanten verwendeten, setzt dieses Werk die hochgradigen Fehlergrenzen explizit mit geometrischen Gittereigenschaften (Dicke und Regularität) in Verbindung.
• Rechtfertigung der Protektion: Die Ergebnisse rechtfertigen mathematisch den Einsatz von geschützten Delaunay-Gittern in höheren Dimensionen. Durch Sicherstellung eines nicht-verschwindenden Schutzes $\delta$ wird gewährleistet, dass die minimale Dicke von Null verschieden bleibt, wodurch verhindert wird, dass die Fehlerkonstanten beliebig groß werden.
• Praktische Implikation: Die Autoren geben an, dass diese Schätzungen „mathematisch fundierte Ansätze zur Reduzierung numerischer Fehler“ in PDE-Simulationen liefern. Sie argumentieren, dass Praktiker, die Finite-Elemente-Methoden für die Gradienteninterpolation oder Vektorfeld-Approximation einsetzen, geschützte Delaunay-Gitter ernsthaft in Betracht ziehen sollten, um Stabilität und Genauigkeit zu gewährleisten, insbesondere wenn hochgradige Polynome verwendet werden.

Die Arbeit bleibt in ihrem Umfang bescheiden und merkt an, dass diese Ergebnisse für hinreichend glatte Funktionen (stückweise $W^{2,p}$ oder $W^{1,m}$) gelten und dass für weniger reguläre Lösungen andere Strategien erforderlich sein könnten. Die Arbeit schlägt keine neuen Algorithmen zur Generierung von Gittern vor, sondern liefert vielmehr die theoretische Fehleranalyse, die den Einsatz bestehender Techniken zur Generierung geschützter Gitter validiert.