Projection-free approximation of flows of harmonic maps with quadratic constraint accuracy and variable step sizes

Titel

Projectievrije benadering van stromen van harmonische kaarten met kwadratische constraint-accuraatheid en variabele tijdstappen.

1. Probleemstelling

Het artikel richt zich op de numerieke approximatie van stromen van harmonische kaarten (harmonic map flows) naar een eenheidsbol. Dit is een fundamenteel probleem in de meetkundige analyse en verschijnt in diverse toepassingen, zoals magnetisatiedynamica (Landau-Lifshitz-Gilbert vergelijking), modellen voor nematische vloeibare kristallen en meetkundige evolutievergelijkingen.

De kern van het probleem is het vinden van een vectorveld $u$ dat de Dirichlet-energie minimaliseert onder de unit-length constraint (eenheidslengte-beperking):
$$|u(x)| = 1 \quad \text{bijna overal in } \Omega.$$
De bijbehorende Euler-Lagrange-vergelijkingen en de gradient flow (warmtestroom) vereisen dat deze constraint exact wordt gehandhaafd. Traditionele numerieke methoden hebben vaak te kampen met:

• Projectie-stappen: Methoden die de constraint na elke stap projecteren, vereisen vaak restricties op de tijdstapgrootte of het type mesh.
• Linearisatie: Newton-methoden convergeren vaak lokaal en vereisen de oplossing van niet-lineaire systemen.
• Constraint schending: Methoden die de constraint niet exact handhaven, kunnen leiden tot fouten die lineair groeien met de tijdstap ($O(\tau)$), wat de nauwkeurigheid beperkt.

2. Methodologie

De auteurs ontwikkelen en analyseren een projectievrije, lineair impliciete methode gebaseerd op de zwakke formulering van Alouges.

De $(\theta, \mu)$-methode

In plaats van een projectie uit te voeren, wordt de constraint geïmplementeerd via een linearisatie van de orthogonaliteitsvoorwaarde. De methode gebruikt twee parameters:

• $\theta$: Regelt de energie-dissipatie (tijdinterpolatie).
• $\mu$: Regelt hoe de constraint wordt benaderd (extrapolatie).

De iteratieve stap voor $n \geq 2$ zoekt een differentiequotiënt $dt u_n$ die voldoet aan:

1. Linearized unit-length condition: $dt u_n \cdot \hat{u}_{n-1+\mu} = 0$, waarbij $\hat{u}_{n-1+\mu}$ een geëxtrapoleerde waarde is.
2. Variationale vergelijking: $(dt u_n, v)_* + (\nabla u_{n-1+\theta}, \nabla v) = 0$ voor alle testfuncties $v$ die orthogonaal zijn op $\hat{u}_{n-1+\mu}$.

Hierbij is $u_{n-1+\theta}$ een geïnterpoleerde waarde en $\hat{u}_{n-1+\mu}$ een geëxtrapoleerde waarde gebaseerd op vorige stappen.

Speciale gevallen

• Implicit Euler ($\theta=1, \mu=0$): Eerste-orde nauwkeurigheid in de constraint.
• Midpoint-methode ($\theta=1/2, \mu=1/2$): De voorgestelde methode die tweede-orde nauwkeurigheid belooft.
• Gewijzigde Implicit Euler ($\theta=1, \mu=1/2$): Behoudt de energie-afname van Euler maar bereikt tweede-orde constraint-accuraatheid.

Variabele tijdstappen

Een belangrijk aspect van het artikel is de uitbreiding naar variabele tijdstappen ($\tau_n$). Dit is cruciaal voor:

• Versnelling van de convergentie naar stationaire punten (door de stapgrootte te vergroten).
• Verbetering van de nauwkeurigheid bij singulariteiten (door de stapgrootte te verkleinen).
  De auteurs bewijzen dat de methode onvoorwaardelijk energiestabiel blijft en de constraint-eigenschappen behoudt, zelfs met variabele stappen. Dit is een significant voordeel ten opzichte van BDF-methoden (zoals BDF2), waar stabiliteit bij variabele stappen een delicaat probleem is.

3. Belangrijkste Bijdragen

1. Kwadratische Constraint-Accuraatheid:
   De methode (specifiek de midpoint-variant met $\theta=\mu=1/2$) bereikt een fout in de constraint van orde $O(\tau^2)$, mits een scherpe discrete regulariteitsvoorwaarde geldt. Dit is een verbetering ten opzichte van de standaard lineair impliciete Euler-methode ($O(\tau)$).

2. Onvoorwaardelijke Energie-Stabiliteit:
   De methode garandeert dat de discrete Dirichlet-energie monotoon afneemt, ongeacht de grootte van de tijdstap. Dit is essentieel voor de convergentie naar stationaire oplossingen.

3. Flexibiliteit voor Variabele Stappen:
   De auteurs presenteren een analyse en algoritmen voor variabele tijdstappen die de stabiliteit en constraint-eigenschappen behouden. Ze introduceren een update-regel ($\tau_{n+1} = \tau_n \sqrt{1 + c\tau_n}$) om de stapgrootte te optimaliseren tijdens de iteratie.

4. Vergelijking met Bestaande Methodes:
   De methode wordt vergeleken met de lineair impliciete Euler-methode en de twee-staps BDF2-methode. De resultaten tonen aan dat de nieuwe methode vergelijkbare stabiliteit biedt als BDF2, maar met een robuustere prestatie bij variabele stappen en een lagere constraint-fout.

4. Numerieke Resultaten

De auteurs voeren uitgebreide numerieke experimenten uit in 2D met lineaire eindige elementen.

• Constante Stappen:

  • Voor de $H^1$-gradient flow vertoont de midpoint-methode tweede-orde convergentie in zowel de $L^1$- als $L^\infty$-fout voor de constraint, terwijl de Euler-methode slechts eerste orde toont.
  • De constraint-fout van de midpoint-methode is ongeveer 3 keer kleiner dan die van de BDF2-methode bij dezelfde stapgrootte.
  • Voor de $L^2$-gradient flow is de convergentie iets trager (pre-asymptotisch gedrag), maar de tweede orde wordt uiteindelijk bereikt.

• Singulariteiten:
  Bij testen met een initiële waarde die leidt tot singulariteiten (blow-up), gedraagt de methode zich robuust. De constraint-fout convergeert superlineair, hoewel de tweede orde wordt beperkt door de onbegrensde regulariteit van de oplossing.

• Variabele Stappen:

  • Versnelling: Bij het minimaliseren van energie (H1-flow) leidt het gebruik van variabele stappen (waarbij de stapgrootte toeneemt naarmate de iteratie vordert) tot een aanzienlijke verbetering. Voor een vergelijkbaar aantal iteraties is de constraint-fout met variabele stappen veel kleiner (orde $10^{-8}$ versus $10^{-6}$ bij constante stappen).
  • Adaptiviteit: Een adaptieve strategie die de stapgrootte verkleint bij snelle dynamiek (bij singulariteiten) en vergroot bij rustige fasen, reduceert de totale rekentijd met ongeveer 20% zonder verlies aan nauwkeurigheid.

5. Betekenis en Conclusie

Dit werk biedt een krachtig alternatief voor bestaande methoden om harmonische kaartstromen te simuleren. De belangrijkste doorbraken zijn:

• Het elimineren van de noodzaak voor dure projectiestappen of restrictieve mesh-voorwaarden.
• Het bereiken van tweede-orde nauwkeurigheid in de unit-length constraint, wat essentieel is voor lange-termijn simulaties.
• Het mogelijk maken van variabele tijdstappen zonder in te leveren op stabiliteit, wat de efficiëntie van de simulatie aanzienlijk verhoogt, vooral bij problemen met singulariteiten of langzame convergentie naar stationaire toestanden.

De methode is niet alleen van toepassing op harmonische kaarten, maar kan worden uitgebreid naar een breed scala aan geometrisch beperkte partiële differentiaalvergelijkingen, zoals in micromagnetisme en vloeibaar kristaltheorie.