A penalized ϕ-FEM scheme for the Poisson Dirichlet problem

Technische Zusammenfassung: Ein penalisiertes $\phi$-FEM-Schema für das Poisson-Dirichlet-Problem

Problemstellung
Die Arbeit befasst sich mit der numerischen Lösung der Poisson-Gleichung mit homogenen Dirichlet-Randbedingungen auf komplexen oder sich entwickelnden Geometrien. Die primäre Herausforderung bei solchen Problemen ist die Erzeugung randangepasster Gitter, was einschränkend wirken kann. Die Autoren konzentrieren sich auf das Poisson-Problem:
$$\begin{cases} -\Delta u = f & \text{in } \Omega, \\ u = 0 & \text{auf } \Gamma, \end{cases}$$
wobei die Domäne $\Omega$ und ihr Rand $\Gamma$ implizit durch eine Level-Set-Funktion $\phi$ beschrieben werden, sodass $\Omega = \{\phi < 0\}$ und $\Gamma = \{\phi = 0\}$.

Methodik
Die Autoren schlagen eine penalisierte Variante der $\phi$-FEM ($\phi$-Finite-Elemente-Methode) vor, einer unangepassten Finite-Elemente-Methode, die randangepasste Gitter vermeidet. Im Gegensatz zur ursprünglichen $\phi$-FEM-Formulierung, welche die Unbekannte über der gesamten Domäne über $u = \phi w$ neu definiert, erzwingt das vorgeschlagene Schema Randbedingungen schwach durch einen Penalisierungsterm.

Zentrale methodische Merkmale sind:

1. Lokale Anforderung an die Level-Set-Funktion: Die Methode erfordert die Level-Set-Funktion $\phi$ nur auf den Zellen, die dem Rand benachbart sind ($\Omega_h^\Gamma$), anstatt global. Dies wird durch die Einführung einer Hilfsvariablen $p$ erreicht, die auf den randbenachbarten Zellen definiert ist, sodass $u = \phi p$ mittels Penalierung erzwungen wird.
2. Variationsformulierung: Das Schema sucht nach $(u_h, p_h) \in V_h^{(k)} \times Q_h^{(k)}$, die eine Variationsgleichung erfüllen, welche Folgendes umfasst:
   • Standard-Gradiententerme über dem Hintergrundgitter $\Omega_h$.
   • Einen Randintegralterm.
   • Einen Penalisierungsterm $\frac{\gamma}{h^2} \int_{\Omega_h^\Gamma} (u_h - \frac{1}{h}\phi_h p_h)(v_h - \frac{1}{h}\phi_h q_h)$, um die Beziehung $u \approx \phi p$ zu erzwingen.
   • Einen „Ghost Penalty“-Stabilisierungsterm ($G_h^{lhs}$), der Sprünge von Normalenableitungen und Laplacen auf den randbenachbarten Facetten und Zellen beinhaltet, um Koerzitivität zu gewährleisten.
3. Stabilisierung: Das Schema nutzt eine Ghost-Penalty-Technik mit den Parametern $\gamma$ und $\sigma_D$, um das System zu stabilisieren, insbesondere für Elemente, die vom Rand geschnitten werden.

Wesentliche Beiträge und theoretische Ergebnisse
Die Arbeit liefert eine rigorose numerische Analyse dieses penalisierten Schemas, das zuvor zwar eingeführt, aber nicht analysiert wurde. Die wichtigsten theoretischen Beiträge sind:

• Koerzitivität: Die Autoren beweisen, dass die zugehörige bilineare Form unter geeigneten Annahmen über das Gitter und die Level-Set-Funktion koerzitiv ist, sofern die Stabilisierungsparameter $\gamma$ und $\sigma_D$ ausreichend groß gewählt werden.
• A-priori-Feşabschätzungen:
  • $H^1$-Halbnorm: Das Schema erreicht eine optimale Konvergenzordnung von $O(h^k)$, wobei $k$ der Polynomgrad des Finite-Elemente-Raums ist.
  • $L^2$-Norm: Die Analyse ergibt eine quasi-optimale Konvergenzrate von $O(h^{k+1/2})$.
  • Regularität: Die Analyse erfordert eine geringere Regularität der Lösung $u$ als die ursprüngliche $\phi$-FEM-Analyse in der bisherigen Literatur.
• Konditionierung: Theoretische Ergebnisse und numerische Experimente legen nahe, dass sich die Konditionszahl der zugehörigen Matrix als $O(h^{-2})$ skaliert.

Numerische Ergebnisse
Die theoretischen Befunde werden durch numerische Experimente validiert, die in Python unter Verwendung der Dolfinx-Bibliothek implementiert wurden:

1. 2D-Komplexe Geometrie: Eine leberähnliche Geometrie wurde getestet. Das penalierte Schema („Dual $\phi$-FEM“) zeigte optimale Konvergenzraten sowohl in der $L^2$- als auch in der $H^1$-Norm, vergleichbar mit der Standard-FEM (fitted FEM) und der ursprünglichen „Direkten $\phi$-FEM“.
2. 3D-Geometrie: Eine sphärische Domäne wurde mit einer künstlich erzeugten Lösung getestet. Auch hier wurde eine optimale Konvergenz beobachtet.
3. Leistungsvergleich: Während beide $\phi$-FEM-Varianten schneller als die Standard-FEM (fitted FEM) waren (aufgrund der Vermeidung komplexer Netzgenerierung), verursachte das penalierte Schema aufgrund der vergrößerten Systemgröße durch die zusätzliche Variable $p$ höhere Rechenkosten als die direkte $\phi$-FEM. Das penalierte Schema bietet jedoch eine höhere Robustheit gegenüber der Level-Set-Funktion, da es keine globalen Gradienten oder zweiten Ableitungen von $\phi$ erfordert.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass dieser penalierte Ansatz erfolgreich die Robustheit gemischter Formulierungen mit der Flexibilität von Immersed-Finite-Elemente-Methoden kombiniert. Ihre primäre Bedeutung liegt in:

• Reduzierte Abhängigkeit von der Level-Set-Funktion: Durch die Beschränkung der Anforderungen an die Level-Set-Funktion (und deren Ableitungen) auf die randbenachbarten Zellen ist die Methode robuster gegenüber Singularitäten der Level-Set-Funktion (z. B. bei Signed Distance Functions) und erfordert keine globale Rekonstruktion der unbekannten Lösung.
• Kompatibilität: Die lokale Natur der Level-Set-Anforderung macht das Schema kompatibel mit Neumann-Randbedingungen, wie in verwandten Arbeiten dargelegt.
• Validierung: Die Arbeit liefert die erste rigorose Fehleranalyse für diese spezifische penalierte Variante und bestätigt deren theoretische Konvergenzeigenschaften sowie die praktische Anwendbarkeit für komplexe Geometrien.

Die Autoren merken an, dass während die theoretische $L^2$-Abschätzung quasi-optimal ist, numerische Experimente im letzten Abschnitt darauf hindeuten, dass in der Praxis eine optimale Konvergenzordnung von $O(h^{k+1})$ erreicht wird.