Finite element methods for isometric embedding of Riemannian manifolds

Titel: Eindige-elementenmethoden voor isometrische inbedding van Riemanniaanse variëteiten

Auteurs: Guangwei Gao, Kaibo Hu, Buyang Li, en Ganghui Zhang.

---

1. Probleemstelling

Het artikel adresseert het fundamentele probleem in de differentiaalmeetkunde van de isometrische inbedding van Riemanniaanse variëteiten.

• Definitie: Gegeven een Riemanniaanse variëteit $M$ met een intrinsieke metriek $g$, is het doel een inbeddingsafbeelding $r: M \to \mathbb{R}^N$ te vinden zodat de pullback van de Euclidische metriek onder $r$ exact overeenkomt met $g$ (d.w.z. $dr \odot dr = g$).
• Specifiek geval: De auteurs focussen op Weyls probleem: het vinden van een globale isometrische inbedding van een tweedimensionale variëteit met positieve Gaussische kromming in de Euclidische ruimte $\mathbb{R}^3$.
• Uitdagingen:
  • Het probleem wordt beschreven door een niet-lineair en degenererend stelsel van partiële differentiaalvergelijkingen (PDE's).
  • Bestaande existentiebewijzen (zoals die van Nash, Lewy en Nirenberg) zijn niet-constructief; ze garanderen dat een oplossing bestaat, maar bieden geen algoritme om deze numeriek te berekenen.
  • De degeneratie komt voort uit de rigiditeit: elke oplossing blijft geldig na compositie met een stijve beweging (rotatie en translatie), wat leidt tot een niet-triviale kern in de lineaire systemen.
  • Er ontbreekt tot nu toe een rigoureuze numerieke analyse met convergentie- en foutenschattingen voor dit probleem.

2. Methodologie

De auteurs stellen een nieuwe numerieke benadering voor die de kloof tussen theoretische existentie en computationele realisatie overbrugt.

A. Dynamische Formulering (Continuïteitsmethode)

In plaats van het statische probleem direct op te lossen, transformeren de auteurs het probleem naar een dynamisch proces, geïnspireerd door de continuïteitsmethode:

1. Ze construeren een continue familie van metrieken $g(t)$ die de beginmetriek (bijv. van een eenheidssfeer) verbindt met de doelmetriek $g$.
2. In plaats van de inbedding $r(t)$ direct te zoeken, differentiëren ze de isometrische voorwaarde $dr \odot dr = g(t)$ naar de tijd.
3. Dit leidt tot een lineair PDE-stelsel voor de snelheid $v = \partial_t r$:
   $$2 dr \odot dv = \partial_t g$$
   Dit lineariseert het probleem aanzienlijk, maar introduceert een degeneratie: de oplossing voor $v$ is uniek alleen tot op infinitesimale stijve bewegingen.

B. Variatieformulering en Saddle-Point Systeem

Om de degeneratie op te lossen en een unieke oplossing te garanderen, stellen de auteurs een nieuwe variatieformulering voor:

• Ze zoeken de snelheid $v$ die orthogonaal is ten opzichte van de ruimte van infinitesimale stijve bewegingen ($\mathcal{R}_M[r]$).
• Dit leidt tot een zadelpunt-systeem (saddle-point system) met een Lagrange-multiplicator $\lambda$ om de orthogonaliteitsvoorwaarde af te dwingen:
  $$\begin{cases} 2 (Dr v, Dr q)_M + (\lambda, q)_M = (\partial_t g, Dr q)_M \\ (v, \mu)_M = 0 \end{cases}$$
  waarbij $q$ testfuncties zijn en $\mu$ elementen uit de ruimte van stijve bewegingen.

C. Eindige Elementen Discretisatie

Het artikel introduceert een specifieke combinatie van eindige elementen die goed past binnen het kader van differentiaalcomplexen:

• Snelheid en Inbedding ($v, r$): Benaderd met hoge-orde Lagrange-elementen (polynoomgraad $k \geq 5$).
• Metriek Tensor: Benaderd met Regge-eindige elementen. Dit zijn speciaal ontworpen elementen voor het discretiseren van metriek-tensoren op simpliciale triangulaties, waarbij de tangentiële-tangentiële componenten continu zijn over de randen.
• Kern: De keuze van Lagrange- en Regge-elementen respecteert de discrete elasticiteitscomplexen, wat essentieel is voor de stabiliteit.

3. Belangrijkste Bijdragen en Theoretische Resultaten

1. Nieuwe Variatieformulering: De eerste formulering die het isometrische inbeddingsprobleem omzet in een goed gesteld lineair systeem voor de snelheid, geschikt voor hoge-orde discretisatie.
2. Korn-ongelijkheid op Variëteiten: De auteurs bewijzen een nieuwe Korn-ongelijkheid voor gesloten tweedimensionale Riemanniaanse variëteiten met positieve Gaussische kromming. Dit is cruciaal om de coerciviteit van het lineaire systeem te garanderen.
3. Discrete Korn-ongelijkheid: Ze bewijzen een discrete versie van deze ongelijkheid voor de eindige elementenruimte. Dit garandeert de welgesteldheid (bestaan en uniciteit) van de numerieke oplossing en karakteriseert de discrete infinitesimale rigiditeit.
4. Convergentie en Foutenschattingen:
   • Ze bewijzen dat het numerieke schema convergeert met een foutenschatting van de orde $O(h^k)$, waarbij $h$ de meshgrootte is en $k$ de polynoomgraad.
   • Een belangrijke technische uitdaging was het ontbreken van volledige $H^1$-coerciviteit (de coerciviteit degradeert in de normale richting). De auteurs overwinnen dit door inverse-ongelijkheden te gebruiken om $W^{1,\infty}$- en $H^1$-normen van de fout te reduceren tot $L^2$-normen, wat mogelijk is dankzij het gebruik van hoge-orde elementen ($k \geq 5$).
5. Toepassing op Ricci-flow: Het raamwerk wordt uitgebreid naar de isometrische inbedding van evoluerende metrieken, specifiek de Ricci-flow. Dit maakt het mogelijk om de evolutie van intrinsieke kromming visueel weer te geven in $\mathbb{R}^3$.

4. Numerieke Experimenten

De auteurs valideren hun theorie met drie voorbeelden:

• Convergentietesten: Voor een evoluerende metriek op een ellipsoïde wordt de theoretische convergentieorde $O(h^k)$ bevestigd voor $k=5$ en $k=6$.
• Inbedding van Revolutiemetrieken: Simulatie van de inbedding van een oppervlak van revolutie dat evolueert van een eenheidssfeer naar een complexer oppervlak, waarbij de geometrische vervorming nauwkeurig wordt weergegeven.
• Visualisatie van de Genormaliseerde Ricci-flow: De methode wordt gebruikt om de evolutie van een oppervlak onder de Ricci-flow te visualiseren. Het resultaat toont hoe het oppervlak evolueert naar een ronde sfeer (constante kromming), wat de effectiviteit van de methode voor dynamische meetkundige stromen demonstreert.

5. Betekenis en Conclusie

Dit werk vormt een fundamentele doorbraak in de computationele differentiaalmeetkunde:

• Het biedt het eerste rigoureuze numerieke raamwerk met bewezen convergentie voor isometrische inbeddingen.
• Het lost het probleem van de degeneratie op door een slimme variatieformulering en het gebruik van specifieke eindige elementen (Regge).
• Het opent de deur voor numerieke toepassingen in diverse gebieden, waaronder:
  • Algemene relativiteit: Berekening van quasi-lokale energie en geometrie van zwarte gaten.
  • Elasticiteit: Modellering van dunne schalen en niet-lineaire buiging.
  • Computergraphics: Vormontwerp en geometrische verwerking.
  • Data-analyse: Inbedding van hoog-dimensionale data in laag-dimensionale variëteiten.

Samenvattend, transformeren de auteurs een klassiek, niet-constructief wiskundig probleem in een robuust, convergent numeriek algoritme, waarbij ze de theoretische diepgang van de differentiaalmeetkunde koppelen aan praktische computationele methoden.