Error analysis for a Finite Element Discretization of a corotational harmonic map heat flow problem

Titel: Fehleranalyse für eine Finite-Elemente-Diskretisierung eines korotationalen harmonischen Abbildungs-Heat-Flow-Problems

Autoren: Nam Anh Nguyen und Arnold Reusken

1. Problemstellung

Das Papier untersucht das harmonische Abbildungs-Heat-Flow-Problem (HMHF) für einen speziellen korotationalen Fall.

• Allgemeines HMHF: Gesucht ist eine Abbildung $u: \Omega \to S^2$ (Einheitskugel), die als Gradientenfluss der Dirichlet-Energie $E(v) = \frac{1}{2}\int_\Omega |\nabla v|^2 dx$ unter der Nebenbedingung $|u|=1$ definiert ist. Die Gleichung lautet $\partial_t u = \Delta u + |\nabla u|^2 u$.
• Spezifischer Fall: Der Fokus liegt auf dem Fall, wo $\Omega$ die Einheitskreisscheibe in $\mathbb{R}^2$ ist und die Lösung eine korotationale Symmetrie aufweist. Durch Verwendung von Polarkoordinaten $(r, \psi)$ lässt sich das Vektorfeld $u$ auf eine skalare Funktion $u(t, r)$ reduzieren.
• Reduzierte Gleichung: Das Problem wird auf eine skalare parabolische Gleichung auf dem Intervall $I=[0,1]$ reduziert:
  $$\partial_t u = \partial_{rr} u + \frac{1}{r} \partial_r u - \frac{\sin(2u)}{2r^2}$$
  mit Randbedingungen $u(t,0)=u(t,1)=0$.
• Bedeutung: Dieser Fall ist fundamental für das Verständnis von Singularitäten (Blow-up) in der HMHF-Theorie. Für bestimmte Anfangsdaten existiert eine globale glatte Lösung, während bei anderen Daten die Lösung in endlicher Zeit singulär wird (die Ableitung an $r=0$ geht gegen unendlich).

2. Methodik

Die Autoren entwickeln und analysieren ein numerisches Verfahren zur Diskretisierung des reduzierten Problems (1.4).

• Raumdiskretisierung: Es wird eine $H^1$-konforme Finite-Elemente-Methode verwendet. Der Raum $S^k_{h,0}$ besteht aus stückweisen Polynomen vom Grad $k \ge 1$, die die homogenen Randbedingungen erfüllen.
• Zeitdiskretisierung: Es wird ein semi-implizites Euler-Verfahren (Backward Euler) angewendet.
  • Der lineare Teil (Diffusion) wird implizit behandelt.
  • Der nichtlineare Term $f(u)/r^2$ (mit $f(u) = u - \frac{1}{2}\sin(2u)$) wird explizit aus dem vorherigen Zeitschritt berechnet (semi-implizit).
  • Vorteil: Dies führt in jedem Zeitschritt zu einem linearen Gleichungssystem, was den Rechenaufwand im Vergleich zu voll-impliziten Newton-Verfahren reduziert.
• Variationsformulierung: Das Problem wird in einem gewichteten Raum $L^2_r$ und $H^1_r$ formuliert, wobei die Gewichtung durch den Faktor $r$ (aus den Polarkoordinaten) erfolgt.

3. Schlüsselbeiträge und theoretische Analyse

Der Hauptbeitrag des Papers ist eine rigorose Fehleranalyse, die optimale Konvergenzordnungen unter der Annahme glatter Lösungen liefert.

• Energiezerlegung und Konvexität:
  • Die Energie $E(u)$ ist nicht konvex. Die Autoren nutzen jedoch eine Energie-Splitting-Strategie:
    $$E(v) = \frac{1}{2}\|v\|_{H^1_r}^2 - \int_0^1 F(v(r)) \frac{dr}{r}$$
    wobei $F(v) = \frac{1}{2}(v^2 - \sin^2 v)$.
  • Ein entscheidender Befund ist, dass die Funktion $F$ konvex ist ($F''(z) = 1-\cos(2z) \ge 0$). Diese verborgene Konvexität ist essenziell für die Stabilitätsanalyse.
• Diskrete Energiedissipation:
  • Es wird bewiesen, dass das diskrete Verfahren die Energiedissipationseigenschaft des kontinuierlichen Problems nachahmt (Theorem 3.1).
  • Dies führt zu einer diskreten Stabilität: Die diskrete Lösung bleibt im Bereich, in dem die Energie beschränkt ist, was eine $L^\infty$-Schranke ($|u| \le \pi$) garantiert.
• Fehleranalyse:
  • Die Analyse nutzt die klassische Aufspaltung des Fehlers mittels der Galerkin-Projektion.
  • Eine große Herausforderung ist der Term $1/r^2$, der bei $r \to 0$ singulär wird. Die Autoren beweisen spezielle Abschätzungen für die Galerkin-Projektion von $u/r$ in der $L^\infty$-Norm (Lemma 4.3).
  • Durch Kombination der Konvexitätseigenschaft von $F$ mit diesen Abschätzungen gelingt es, den Linearisierungsfehler zu kontrollieren.
• Ergebnis: Es werden optimale Fehlerabschätzungen hergeleitet:
  $$\|u_h - u\|_{H^1_r} \le C(\tau + h^k)$$
  wobei $\tau$ die Zeitschrittweite und $h$ die Gitterweite ist.

4. Numerische Ergebnisse

Die theoretischen Vorhersagen werden durch numerische Experimente validiert:

• Konvergenzordnung: Für lineare ($k=1$) und quadratische ($k=2$) Elemente wurden die Konvergenzraten in den Normen $\|\cdot\|_{H^1_r}$ und $\|\cdot\|_{L^2_r}$ gemessen.
  • Die Ergebnisse zeigen die erwartete optimale Konvergenzordnung ($O(h^k)$ im $H^1$-Fehler und $O(h^{k+1})$ im $L^2$-Fehler).
  • Die zeitliche Konvergenz ist von erster Ordnung ($O(\tau)$), wie vom Euler-Verfahren erwartet.
• Energiedissipation: Die Simulationen bestätigen die theoretische Vorhersage, dass die diskrete Energie monoton abnimmt.

5. Signifikanz und Einschränkungen

• Beitrag zur Literatur: Während es bereits Arbeiten zu Finite-Differenzen-Methoden für dieses Problem gibt, ist dies eine der ersten detaillierten Fehleranalysen für eine Finite-Elemente-Diskretisierung im korotationalen Fall.
• Einschränkungen:
  • Die Analyse setzt glatte Lösungen voraus (keine Blow-up-Szenarien).
  • Die Methode und die Analyse lassen sich auf höhergradige korotationale Fälle erweitern, scheitern jedoch an der äquivarianten Version (m-equivariant) oder der Landau-Lifshitz-Gilbert (LLG) Gleichung. Der Grund liegt darin, dass diese Probleme die spezifische konvexe Struktur und die entsprechende Energie-Splitting-Eigenschaft nicht besitzen, die für die Stabilitätsbeweise in diesem Papier entscheidend sind.

Fazit: Das Papier liefert einen robusten theoretischen Rahmen für die numerische Simulation des korotationalen HMHF-Problems mittels Finite-Elemente-Methoden und beweist, dass die gewählte semi-implizite Diskretisierung sowohl stabil als auch optimal konvergent ist, solange die Lösung glatt bleibt.