Stable $C^1$-conforming finite element methods for a class of nonlinear fourth-order evolution equations

Technische Samenvatting: Stabiele $C^1$-conforme Eindelementmethoden voor een Klasse van Nietlineaire Vierde-orde Evolutievergelijkingen

Probleemformulering
Het artikel behandelt de numerieke analyse van een algemene klasse van nietlineaire vierde-orde evolutie partiële differentiaalvergelijkingen (PDE's) gedefinieerd op een begrensd domein $\Omega \subset \mathbb{R}^d$ ($d \in \{1, 2, 3\}$). De leidende vergelijking wordt gegeven door:
$$\frac{\partial u}{\partial t} - \beta_1 \Delta u + \beta_2 \Delta^2 u = \beta_3(1 - |u|^2)u - \beta_4 u \times \Delta u + \beta_5 \Delta(|u|^2 u) + \beta_6 (j \cdot \nabla)u$$
onderworpen aan Neumann-randvoorwaarden ($\frac{\partial u}{\partial n} = 0, \frac{\partial (\Delta u)}{\partial n} = 0$) en een beginvoorwaarde $u(0) = u_0$. Hierbij kan $u$ scalair-waardig zijn ($m=1$) of vector-waardig ($m=3$). De vergelijking omvat diverse fysische modellen afhankelijk van de coëfficiënten $\beta_i$ en de dimensie $d$, waaronder:

• Micromagnetica: De Landau–Lifshitz–Baryakhtar (LLB) vergelijking (met of zonder spin-torque) en regularisaties van de Landau–Lifshitz–Bloch vergelijking.
• Patroonvorming: De Swift–Hohenberg vergelijking.
• Fasescheiding en Diffusie: Convectieve Cahn–Hilliard vergelijkingen met brontermen, anomalieën in bi-flux diffusiemodellen en gegeneraliseerde populatieverspreidingsmodellen.

De auteurs merken op dat hoewel er numerieke methoden bestaan voor specifieke instanties van deze vergelijkingen (bijv. de scalaire Swift–Hohenberg of de standaard Cahn–Hilliard), een verenigd, algemeen numeriek schema voor het volledige systeem (1.1) waarbij alle coëfficiënten potentieel niet nul zijn, ontbrak in de literatuur.

Methodologie
De auteurs stellen volledig gediscretiseerde eindelementmethoden voor en analyseren deze met behulp van $C^1$-conforme elementen. Deze keuze is noodzakelijk vanwege de vierde-orde ruimtelijke afgeleide ($\Delta^2 u$), die vereist dat de eindelementruimte een deelverzameling is van $H^2(\Omega)$.

• Ruimtelijke Discretisatie: De methode maakt gebruik van stuksgewijze polynomen van graad $r$ (waarbij $r=3$ voor $d=1,2$ en $r=9$ voor $d=3$ om $C^1$-continuïteit te behouden). Specifieke genoemde elementen zijn $C^1$-P3 Hermite-elementen (1D), Hsieh–Clough–Tocher (HCT) macro-elementen (2D) en Zenísek-elementen (3D).
• Temporele Discretisatie: Er worden twee gelineraliseerde tijdstap-schema's voorgesteld:
  1. Een semi-impliciete Euler-methode (eerste orde in de tijd).
  2. Een semi-impliciete BDF2-methode (tweede orde in de tijd).
     In beide schema's worden de nietlineaire termen gelineraliseerd door ze te evalueren op het vorige tijdstap, wat resulteert in een stelsel van lineaire vergelijkingen bij elke stap. Dit vermijdt de computationele kosten van het oplossen van nietlineaire stelsels via de Newton-methode.

Belangrijkste Bijdragen en Theoretische Resultaten
Het artikel biedt een rigoureuze foutanalyse en een stabiliteitsbewijs voor de voorgestelde schema's:

1. Stabiliteit:

   • De schema's zijn bewezen onvoorwaardelijk stabiel in de $H^2$-norm voor dimensies $d=1$ en $d=2$. Dit betekent dat stabiliteit standhoudt zonder restricties op de tijdstapgrootte $k$ ten opzichte van de roostermaat $h$, en zonder de aanname van quasi-uniforme triangulaties.
   • Voor $d=3$ zijn de schema's voorwaardelijk stabiel. Stabiliteit is gegarandeerd indien de begingegevens voldoende klein zijn of indien de tijdstap voldoet aan een specifieke restrictie ten opzichte van de initiële energie. De auteurs erkennen dat $C^1$-conforme elementen in 3D (die een graad van 9 polynomen vereisen) computationeel prohibitief zijn, en suggereren dat niet-conforme methoden praktischer zijn voor $d=3$, hoewel de theoretische analyse voor de volledigheid is opgenomen.

2. Foutschattingen en Convergentie:

   • De auteurs vestigen optimale convergentiesnelheden in de $L^\infty(0, T; H^\beta)$ normen voor $\beta = 0, 1, 2$.
   • Voor het semi-discrete schema is de fout begrensd door $O(h^{r+1-\beta})$.
   • Voor de volledig gediscretiseerde schema's:
     • De semi-impliciete Euler-methode bereikt een fout van $O(h^{r+1-\beta} + k)$.
     • De semi-impliciete BDF2-methode bereikt een fout van $O(h^{r+1-\beta} + k^2)$.
   • De analyse steunt op de ontwikkeling van een gespecialiseerde elliptische projectie, aangepast aan het nietlineaire probleem, en zorgvuldige schatting van nietlineaire termen met behulp van Sobolev-inbeddingen en inverse ongelijkheden.

3. Numerieke Experimenten:

   • De theoretische convergentiesnelheden werden geverifieerd met behulp van de open-source package FreeFEM met HCT-elementen in 2D.
   • Experimenten werden uitgevoerd op de Landau–Lifshitz–Baryakhtar vergelijking (vector-waardig) en de Cahn–Hilliard vergelijking met een bronterm (scalar-waardig).
   • De resultaten bevestigden de voorspelde convergentieordes: 4de orde in $L^2$, 3de orde in $H^1$, en 2de orde in $H^2$ voor de ruimtelijke discretisatie, consistent met het gebruik van kubische ($r=3$) elementen.

Betekenis en Claims
Het artikel claimt het eerste numerieke schema te bieden dat in staat is de oplossing van de algemene nietlineaire vierde-orde systeem (1.1) in al zijn algemeenheid te benaderen. De belangrijkste voordelen zijn:

• Lineariteit: Door de nietlineaire termen te lineariseren, vermijdt de methode de complexiteit en potentiële instabiliteit van het oplossen van nietlineaire stelsels bij elke tijdstap, wat voordelig is voor toepassingen zoals micromagnetische simulaties.
• Stabiliteit: De onvoorwaardelijke stabiliteit in $H^2$ voor $d=1,2$ zorgt voor robuuste prestaties zonder strikte tijdstaprestricties.
• Nauwkeurigheid: Het gebruik van $C^1$-conforme elementen maakt een hogere orde nauwkeurigheid mogelijk bij het vastleggen van oplossingsgradiënten vergeleken met niet-conforme methoden, ondanks de verhoogde computationele kosten van $C^1$-continuïteit.

De auteurs hanteren een bescheiden toon met betrekking tot de 3D-gevallen, waarbij zij expliciet stellen dat hoewel de theoretische analyse volledig is, de praktische toepassing van $C^1$-conforme elementen in 3D momenteel computationeel onhaalbaar is vanwege de vereiste hoge polynoomgraad. Het werk dient als een fundamenteel theoretisch kader voor het oplossen van complexe nietlineaire vierde-orde evolutievergelijkingen voortvloeiend uit de fysica, materiaalkunde en biologie.