Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Technische Samenvatting: Stroomfunctie – druk virtuele elementmethoden voor het Stokes–Darcy-interfaceprobleem

Probleemstelling
Het artikel behandelt de numerieke simulatie van het gekoppelde Stokes–Darcy-interfaceprobleem, een fundamentele uitdaging in de stromingsleer die de interactie tussen vrije stromingsgebieden (Stokes) en poreuze media (Darcy) beheerst. De koppeling wordt gerealiseerd door continuïteit van normale snelheden (massa-behoud), evenwicht van normale spanningen en de Beavers–Joseph–Saffman tangentiële slipconditie. Hoewel diverse Eindige Element Methode (FE), Discontinuous Galerkin (dG) en Hybrid High–Order (HHO) methoden bestaan, die vaak vertrouwen op snelheid–drukformuleringen die inf–sup stabiele paren vereisen, richt dit werk zich op een specifieke herformulering in tweedimensionale enkelvoudig samenhangende domeinen. Hier wordt de onsamendrukbaarheidsvoorwaarde identiek voldaan door de introductie van een scalaire stroomfunctie, waardoor de drukvariabele uit de Stokes-vergelijkingen wordt geëlimineerd en het systeem wordt gereduceerd tot een vierde-orde biharmonisch probleem gekoppeld aan een tweede-orde Darcy-drukvergelijking.

Methodologie
De auteurs stellen een nieuwe Virtual Element Method (VEM) discretisatie voor van deze stroomfunctie–drukformulering. De methodologie maakt gebruik van het vermogen van VEM om algemene polygonale roosters te verwerken, wat bijzonder voordelig is voor complexe interfacegeometrieën en niet-aangesloten roosters (non-matching meshes).

1. Formulering: Het continue probleem wordt gegoten in een primaal-primaal vorm. Het Stokes-domein gebruikt een stroomfunctie $\chi \in H^2$, terwijl het Darcy-domein een drukpotentiaal $\phi \in H^1$ gebruikt. De zwakke formulering resulteert in een verstoord saddle-point probleem.
2. Discretisatie:
   • Stokes-domein: Een $C^1$-conforme laagste-orde VEM-ruimte wordt gehanteerd voor de stroomfunctie. Dit omvat kwadratische polynomen ($P_2$) en vereist $H^2$-conformiteit, bereikt door verbeterde virtuele ruimten met vrijheidsgraden (DoF's) gedefinieerd door knoopwaarden en knoopgradiënten.
   • Darcy-domein: Een $C^0$-conforme laagste-orde VEM-ruimte wordt gebruikt voor de druk, gebaseerd op lineaire polynomen ($P_1$).
   • Koppeling: De interfacecondities worden natuurlijk afgedwongen binnen de zwakke formulering. De methode faciliteert niet-aangesloten roosters bij de interface zonder herverdeling van het rooster (remeshing).
3. Implementatie: Het schema is geïmplementeerd met behulp van de open-source objectgeoriënteerde bibliotheek vem++. Een specifiek kenmerk, de vem::vemMesh2dDofFilter klasse, wordt gebruikt om DoF's over verschillende subdomeinen te beheren, waardoor het simultaan definiëren van verschillende VEM-ruimten op een enkel rooster mogelijk is en het filteren op basis van elementmarkers mogelijk wordt.
4. Stabilisatie: De discrete bilineaire vormen bevatten stabilisatietermen ($S_S^K$ en $S_D^K$) om coerciviteit te waarborgen op de kern van de projectie-operatoren, conform standaard VEM-praktijken (bijv. diagonaal recept en DOFI-DOFI).

Belangrijkste Bijdragen
Het artikel zet vijf primaire bijdragen uiteen:

• Continue Analyse: Ontwikkeling van een continue analyse voor het Stokes–Darcy-interfaceprobleem in stroomfunctie–drukformulering, die inherent de noodzaak voor inf–sup stabiliteitsvoorwaarden vermijdt.
• Discreet Schema: Ontwerp van een $C^1$–$C^0$ laagste-orde VEM-discretisatie. De auteurs benadrukken de computationele efficiëntie van dit schema, waarbij zij opmerken dat het slechts $3N_S + N_D + 4N_\Sigma$ vrijheidsgraden omvat (waarbij $N_S$, $N_D$ en $N_\Sigma$ respectievelijk de aantallen interne knooppunten in het Stokes-domein, het Darcy-domein en op de interface zijn).
• Foutanalyse: Afleiding van a priori foutschattingen voor de voorgestelde discretisatie, waarbij optimale convergentiesnelheden worden bewezen onder standaard regulariteitsaannames.
• Software Implementatie: Integratie van nieuwe functionaliteit in de vem++ bibliotheek om de eenvoudige implementatie van interfaceproblemen via mesh filtering te faciliteren.
• Numerieke Validatie: Uitvoering van numerieke experimenten om de theoretische bevindingen te bevestigen.

Resultaten
De numerieke experimenten valideren de theoretische convergentiesnelheden en de robuustheid van de methode:

• Gemanifesteerde Oplossingen: Twee experimenten met gemanifesteerde oplossingen (één met nul interfacewaarden, één met niet-nul) op diverse typen roosters (kwadrilateraal, niet-convex, Voronoi, en geperturbeerde Voronoi) tonen een asymptotische lineaire foutafname ($O(h)$), consistent met de theoretische voorspellingen voor het laagste-orde schema.
• Dead-End Filtratie: Een simulatie van een kwart-annulus dead-end filter slaagt erin resultaten uit eerdere studies te reproduceren, wat de capaciteit van de methode aantoont om gebogen interfaces en gemengde randvoorwaarden te behandelen.
• Bioartificiële Organen: Een simulatie van bloedstroom in een bioartificieel orgaan (encapsulatie-gebaseerd) met variërende permeabiliteitsconstanten illustreert het vermogen van de methode om complexe fysiologische stromingen en de invloed van permeabiliteit op de flux over de interface te modelleren.

Betekenis en Claims
Het artikel positioneert het werk als een brug in de literatuur met betrekking tot de uitbreiding van stroomfunctie-formuleringen naar gekoppelde Stokes–Darcy problemen met behulp van Virtual Elements. De auteurs beweren dat de voorgestelde structuur een "volledig inf–sup-vrije strategie" biedt, wat het ontwerp van discrete ruimten vereenvoudigt vergeleken met traditionele snelheid–druk benaderingen. Het vermogen om irreguliere interfacegeometrieën en niet-aangesloten roosters te hanteren zonder herverdeling van het rooster wordt benadrukt als een belangrijk voordeel van het VEM-framework.

De auteurs blijven bescheiden over toekomstige uitdagingen en merken expliciet op dat het uitbreiden van deze stroomfunctie-formulering naar meervoudig samenhangende domeinen (waar de stroomfunctie niet uniek bepaald wordt door standaard randvoorwaarden) en het schalen van het framework naar drie dimensies (wat vectorpotentialen en hogere-orde $C^1$ ruimten voor vector-georiënteerde functies vereist) openstaande uitdagingen blijven. Daarnaast worden de integratie van stabilisatie-vrije technieken en de afleiding van robuuste a posteriori foutschatters geïdentificeerd als richtingen voor toekomstig onderzoek.