Achieving Material Robustness via Symmetric Stress Finite Element Discretizations

Technische Zusammenfassung: Erreichung materieller Robustheit durch symmetrische Spannungs-Finite-Elemente-Diskretisierungen

Problemstellung
In der Kontinuumsmechanik beinhalten Variationsformulierungen wie die Elastizität nach Hellinger–Reissner und Geschwindigkeit–Spannungs-Formulierungen für inkompressible Strömungen symmetrische Spannungstensoren. Die Symmetrie des Cauchy-Spannungstensors ist physikalisch an die Erhaltung des Drehimpulses gebunden. Bei der Diskretisierung dieser Probleme mit $H(\text{div})-konformen$ Finite-Elementen müssen Praktiker wählen, ob sie die Symmetrie stark (indem sie fordern, dass die diskreten Tensoren punktweise symmetrisch sind) oder schwach (indem sie die Symmetrie über einen Lagrange-Multiplikator erzwingen) durchsetzen. Während die schwache Symmetrie in der Literatur aufgrund der wahrgenommenen Komplexität beim Aufbau stark symmetrischer Elemente oft bevorzugt wird, untersucht dieser Beitrag die Auswirkungen dieser Wahl auf die numerische Genauigkeit, insbesondere für Materialien mit anisotropen konstitutiven Gesetzen.

Methodik
Die Autoren entwickeln einen theoretischen Rahmen und eine Reihe numerischer Benchmarks, um stark und schwach symmetrische Schemata zu vergleichen.

1. Mathematische Formulierung: Die Studie betrachtet ein allgemeines linearisiertes statisches Problem sowohl für Fluide als auch für Festkörper, wobei die Spannung $\sigma$ mit der Verzerrung $\varepsilon(u)$ und einem externen Feldterm $\tilde{F}$ (der Anisotropie repräsentiert) in Beziehung steht. Das Problem wird als Sattelpunktssystem formuliert.
2. Diskretisierungsschemata:
   • Stark symmetrisch: Das Johnson-Mercier-Křížek (JMK)-Schema und das Hu-Zhang (HZk)-Schema werden verwendet. Diese setzen $\sigma_h \in \Sigma^{\text{sym}}_h$ punktweise durch.
   • Schwach symmetrisch: Die PEERS- und Arnold-Falk-Winther (AFWk)-Schemata werden verwendet. Diese führen einen Lagrange-Multiplikator $\omega_h$ (der den antisymmetrischen Teil des Gradienten repräsentiert) ein, um die Symmetrie in einem schwachen Sinne durchzusetzen.
3. Benchmarks: Die Autoren konstruieren „hergestellte Lösungen", bei denen die exakte Spannung identisch null ist ($\sigma \equiv 0$), aber die Verschiebungs- und Rotationsfelder nicht-trivial sind. Diese Lösungen werden aus konstitutiven Gesetzen abgeleitet, die inspiriert sind von:
   • Linearen isotropen Festkörpern.
   • Transversal isotropen Festkörpern (faserverstärkte Materialien).
   • Polaren Fluiden (Flüssigkristall-Polymernetzwerke).
     Die Größe des anisotropen Kraftterms wird durch einen Skalierungsparameter $\delta$ gesteuert.
4. Theoretische Analyse: Eine abstrakte Hilbertraum-Theorie wird für allgemeine Sattelpunktprobleme entwickelt. Die Autoren definieren „strukturerhaltende" Diskretisierungen und leiten a-priori-Fehlerabschätzungen ab, um das beobachtete numerische Verhalten zu erklären.

Hauptbeiträge

• Definition der materiellen Robustheit: Der Beitrag führt das Konzept der „materiellen Robustheit" ein, definiert als die Eigenschaft, dass ein numerisches Schema für jede exakte Konfiguration mit null Spannung eine Null-Spannungs-Näherung liefert, unabhängig vom konstitutiven Gesetz oder der Größe der anisotropen Kraft.
• Identifikation von Ausfallmodi: Die Autoren zeigen, dass schwach symmetrische Schemata beliebig schlechte Spannungsapproximationen liefern können (Fehler, die mit dem Kraftparameter $\delta$ skalieren), selbst wenn die exakte Spannung null ist. Dies tritt spezifisch auf, wenn der exakte Rotationstensor $\omega = \text{skw}(\nabla u)$ nicht zum diskreten Lagrange-Multiplikator-Raum $\Xi_h$ gehört.
• Vereinheitlichte Theorie: Eine rigorose Theorie wird präsentiert, die zeigt, dass stark symmetrische Schemata, die bestimmte Kern-Inklusionseigenschaften erfüllen ($\ker B_h \subseteq \ker B$), inhärent materiell robust sind. Im Gegensatz dazu fehlt schwach symmetrischen Schemata diese Invarianz, es sei denn, der diskrete Rotationsraum enthält den exakten Rotationstensor.
• Fehlerabschätzungen: Neue a-priori-Fehlerabschätzungen werden für teilweise strukturerhaltende Diskretisierungen hergeleitet, die explizit zeigen, wie der Spannungsfehler im Fall der schwachen Symmetrie von der Approximationsgüte des Rotationstensors abhängt.

Ergebnisse
Numerische Experimente auf unstrukturierten Gittern bestätigen die theoretischen Vorhersagen:

• Isotroper Fall: Für lineare isotrope Festkörper, bei denen die exakte Rotation konstant ist (und somit leicht darstellbar), performen sowohl starke als auch schwache Schemata ähnlich, wobei die Fehler durch Rundungsfehler dominiert werden.
• Anisotrope Fälle: Für transversal isotrope Festkörper und polare Fluide, bei denen der Rotationstensor komplex ist:
  • Stark symmetrische Schemata (JMK, HZ3): Erzeugen konsistent Spannungsfehler auf dem Niveau der Maschinengenauigkeit (Null-Spannungszustand), unabhängig vom Skalierungsparameter $\delta$.
  • Schwach symmetrische Schemata (PEERS, AFWk): Erzeugen Spannungsfehler, die linear mit $\delta$ skalieren. Selbst bei sehr hochordentlichen Schemata (z. B. AFW3) bleibt der Spannungsfehler groß, wenn der exakte Rotationstensor nicht im diskreten Raum $\Xi_h$ enthalten ist.
• Transientes Verhalten: Ein transientes Beispiel für polare Fluide zeigt, dass zwar beide Schemata anfangs übereinstimmen, das nicht-robuste schwach symmetrische Schema jedoch versagt, sich in den korrekten Null-Spannungs-Gleichgewichtszustand zu entspannen, während das robuste stark symmetrische Schema dies tut.

Bedeutung und Behauptungen
Der Beitrag behauptet, dass die Wahl zwischen starker und schwacher Symmetrieerzwingung nicht lediglich eine Frage der Implementierungskomplexität ist, sondern die physikalische Genauigkeit der numerischen Lösung fundamental beeinflusst. Die Autoren argumentieren, dass die starke Durchsetzung der Symmetrie für die materielle Robustheit entscheidend ist, um sicherzustellen, dass numerische Schemata in spannungsfreien Konfigurationen keine spuriosen Spannungen erzeugen, insbesondere in Anwendungen, die Rotationsinvarianten oder anisotrope Materialien betreffen.

Die Arbeit plädiert für die breitere Anwendung stark symmetrischer Finite-Elemente (wie JMK und Hu-Zhang) und stellt fest, dass diese Elemente zunehmend in Open-Source-Software wie Firedrake und NGSolve verfügbar sind. Die Autoren schließen, dass für Probleme, die empfindlich auf die Erhaltung des Drehimpulses reagieren, stark symmetrische Diskretisierungen eine notwendige Genauigkeitsgarantie bieten, die schwach symmetrische Schemata ohne restriktive Bedingungen an die diskreten Räume nicht bieten können.