A Divergence-Free Scott-Vogelius Finite Element Method for the Surface Stokes Problem

Technische Samenvatting: Een Divergentievrije Scott-Vogelius Eindige Elementenmethode voor het Oppervlakte-Stokesprobleem

Probleemstelling
Het artikel behandelt het stationaire oppervlakte-Stokesprobleem op een glad, gesloten, oriënteerbaar twee-dimensionaal oppervlak $\gamma \subset \mathbb{R}^3$. De beheersende vergelijkingen omvatten een tangentieel snelheidveld $u$ en een druk $p$, onderworpen aan de beperkingen dat $u$ tangentieel aan het oppervlak is ($u \cdot \nu = 0$) en dat de stroming onsamendrukbaar is ($\text{div}_\gamma u = 0$).

Het construeren van eindige elementenmethoden voor dit probleem brengt twee primaire uitdagingen met zich mee:

1. Tangentialiteit: Het sterk afdwingen van de tangentiële beperking op een discrete geometrie leidt vaak tot "locking"-verschijnselen wanneer standaard $H^1$-conforme vectorruimten worden gebruikt, wat de benaderingseigenschappen verslechtert.
2. Onsamendrukbaarheid: Standaard snelheid-druk paren (zoals Taylor-Hood) dwingen de divergentievrije beperking alleen zwak af. Hoewel er methoden bestaan die tangentialiteit sterk afdwingen (bijv. $H(\text{div})$-conforme BDM-ruimten), resulteren deze vaak in grote algebraïsche systemen of ontbreekt een volledige stabiliteitsanalyse op oppervlakken. Daarentegen, methoden die tangentialiteit sterk afdwingen terwijl ze $H^1$-conformiteit behouden (bijv. MINI of Taylor-Hood op oppervlakken), dwingen onsamendrukbaarheid doorgaans alleen zwak af.

De auteurs streven ernaar een methode te construeren die zowel de tangentiële als de onsamendrukbaarheidsbeperkingen exact afdwingt, de dimensionaliteit van standaard 2D Euclidische discretisaties behoudt, en de implementatiecomplexiteit van macro-elementconstructies op gekromde oppervlakken vermijdt.

Methodologie
De voorgestelde aanpak breidt het oppervlakte eindige elementenkader van Kone et al. [10, 11] uit naar het Scott-Vogelius paar. De methodologie omvat de volgende sleutelcomponenten:

• Triangulatie: De methode maakt gebruik van een Clough-Tocher (CT) verfijning van een polyhedrale benadering van het oppervlak. Elk driehoek in het initiële mesh wordt gesplitst door de barycentrische coördinaten te verbinden met de hoekpunten. Deze verfijning verlicht de vereisten voor de polynoomgraad voor stabiliteit (waardoor $r \ge 2$ is toegestaan) en elimineert "bijna-singuliere" hoekpunten.
• Eindige Elementenruimten:
  • Snelheid ($V_h$): Gedefinieerd als continue, stuksgewijze polynomen van graad $r$ op de CT-verfijnde oppervlaktetriangulatie. Om de zwakke continuïteit over de randen van de elementen af te dwingen (nodig voor $H(\text{div})$-conformiteit op het oppervlak), gebruiken de auteurs een hulp-Piola-transformatie bij de Lagrange vrijheidsgraden (DOFs). Dit verbindt de lokale ruimten zonder dat inter-element integratie of stabilisatietermen nodig zijn.
  • Druk ($Q_h$): Gedefinieerd als discontinue stuksgewijze polynomen van graad $r-1$ op hetzelfde CT-verfijnde mesh.
• Implementatiestrategie: In tegenstelling tot eerdere isoparametrische Scott-Vogelius methoden gebaseerd op macro-element definities (waarbij basisfuncties worden gedefinieerd op een gesubdivideerd referentie-element en globaal worden gemapt), definieert dit artikel de eindige elementenruimten direct op de verfijnde oppervlaktetriangulatie. Deze aanpak vermijdt de complexe datastructuren en gespecialiseerde kwadratuurregels die vereist zijn door macro-element implementaties, wat leidt tot een meer praktische regeling die compatibel is met standaard eindige elementensoftware.
• Stabiliteitsanalyse: De auteurs bewijzen de discrete inf-sup conditie (Ladyzhenskaya-Babuška-Brezzi conditie). Het bewijs maakt gebruik van een perturbatie-argument:
  1. Stabiliteit wordt eerst vastgesteld voor een hulp-"macro-element" paar gedefinieerd op een iets andere hoog-orde triangulatie.
  2. Het verschil tussen de voorgestelde directe verfijningsruimte en de hulp-"macro-element" ruimte wordt gekwantificeerd als een $O(h)$ perturbatie.
  3. Stabiliteit wordt van het hulp-paar overgedragen naar het voorgestelde paar via deze perturbatieanalyse.

Belangrijkste Bijdragen en Resultaten

1. Exacte Divergentievrije en Tangentiële Discretisatie: Het geconstrueerde Scott-Vogelius paar op gekromde Clough-Tocher triangulaties dwingt zowel de tangentialiteit ($u \cdot \nu = 0$) als de onsamendrukbaarheid ($\text{div}_\gamma u = 0$) exact af.
2. Vereenvoudigde Implementatie: Door de ruimten direct op de verfijnde triangulatie te definiëren in plaats van via macro-elementen, bieden de auteurs een schema dat aanzienlijk eenvoudiger te implementeren is, terwijl het de theoretische voordelen van het Scott-Vogelius paar behoudt.
3. Stabiliteit en Convergentie:
   • De methode is bewezen inf-sup stabiel met een constante die onafhankelijk is van de mesh-grootte.
   • Optimale convergentieorde wordt afgeleid in het isoparametrische regime ($k=r$). Specifiek convergeert de snelheid met orde $O(h^r)$ in de $H^1$-norm en $O(h^{r+1})$ in de $L^2$-norm, terwijl de druk convergeert met orde $O(h^r)$ in de $L^2$-norm.
   • De analyse bevat een bewijs van stabiliteit voor het bijbehorende macro Scott-Vogelius paar op oppervlakken en het oppervlak $P_r - P_0$ paar als bijproducten.
4. Druk-robuustheidsanalyse: Het artikel onderzoekt druk-robuustheid (waarbij de snelheidsfout onafhankelijk is van de drukbenadering). De auteurs betogen dat voor oppervlakte-Stokesproblemen druk-robuustheid niet automatisch gegarandeerd is, zelfs voor exact divergentievrije methoden, vanwege geometrische benaderingsfouten en de extensie van de bronterm. Ze bieden een kader voor het construeren van bron-extensies die druk-robuustheid behouden en tonen via numerieke experimenten aan dat standaard extensies deze eigenschap kunnen vernietigen.

Numerieke Experimenten
De auteurs valideren de theorie via drie numerieke experimenten op de eenheidssfeer en een torus:

• Test 1 (Divergentievrije stroming): Bevestigt optimale convergentieordes ($O(h^{r+1})$ voor snelheid $L^2$, $O(h^r)$ voor snelheid gradiënt en druk) en toont aan dat de oplossing divergentievrij is tot op de machineprecisie.
• Test 2 (Geen-stroming/Nul-snelheid): Demonstreert dat hoewel de druk optimaal convergeert, de snelheidbenadering niet-nul is door geometrische fouten, wat de afwezigheid van automatische druk-robuustheid benadrukt. Het experiment observeert ook superconvergentie in de snelheid $L^2$-fout ($O(h^{r+2})$), wat suggereert dat de theoretische $L^2$-schatting mogelijk niet scherp is.
• Test 3 (Torus): Bevestigt optimale convergentieordes op een complexere geometrie, hoewel met enige degradatie toegeschreven aan systeemconditionering en afrondingsfouten.

Betekenis
Het artikel claimt een gat in de literatuur over oppervlaktevloeistofdynamica te vullen door een methode te bieden die de voordelen van exact divergentievrije discretisaties (behoudswetten, potentieel voor druk-robuustheid) combineert met de praktische bruikbaarheid van standaard eindige elementenimplementaties. Door af te stappen van de macro-element constructie, bieden de auteurs een "aanzienlijk eenvoudigere en meer praktische implementatie" die de theoretische strengheid van het Scott-Vogelius paar op gekromde oppervlakken behoudt. Het werk breidt de toepasbaarheid van divergentievrije methoden uit naar oppervlakteproblemen, die cruciaal zijn voor het modelleren van fenomenen zoals lipiden-bilayers en dunne films.