Second order explicit splitting scheme for fluid-poroelastic structure interaction problems

Technische Samenvatting: Tweede-orde expliciet splitsingsschema voor vloeistof-poroelastische structuurinteractieproblemen

Probleemstelling
Het artikel behandelt de numerieke simulatie van vloeistof–poroelastische structuurinteractieproblemen (FPSI), specifiek het tijdgebonden Stokes–Biot-systeem op vaste domeinen. Deze problemen komen voor in contexten zoals biologische weefselperfusie en stroming door vervormbare poreuze materialen. Een centrale uitdaging in dit veld is de efficiënte en stabiele behandeling van de interfacekoppelingvoorwaarden, die continuïteit van de normale flux en de balans van normale en tangentiële spanningen afdwingen. Hoewel monolithische methoden robuust zijn, vereisen zij het oplossen van grote gekoppelde systemen bij elke tijdstap. Gepartitioneerde methoden zijn computationeel aantrekkelijk omdat ze het mogelijk maken om gespecialiseerde solvers voor de vloeistof (Stokes) en de poroelastische (Biot) subproblemen onafhankelijk te gebruiken. Echter, expliciete of losjes gekoppelde gepartitioneerde schema's worstelen vaak met stabiliteit, met name wanneer de interfacekoppeling expliciet wordt behandeld.

Methodologie
De auteurs stellen een volledig gediscretiseerd, tweede-orde expliciet splitsingsschema voor het Stokes–Biot-probleem voor. De methodologie steunt op de volgende componenten:

1. Robin-herformulering: De interfacevoorwaarden worden herschreven in een Robin-vorm met penalty-parameters $\gamma > 0$ (tangentieel) en $L > 0$ (normaal). Deze herformulering ontkoppelt de subproblemen terwijl een consistente koppeling wordt behouden die geschikt is voor gepartitioneerde algoritmen.
2. Tijddiscretisatie: Het schema maakt gebruik van de tweede-orde Backward Differentiation Formula (BDF2) voor de tijdsontwikkeling van de subdomeinvergelijkingen.
3. Expliciete extrapolatie: Om een volledig expliciet, paralleliseerbaar schema te bereiken, worden de interfacegegevens die nodig zijn voor de Robin-randvoorwaarden op tijdstip $t_{n+1}$ benaderd met behulp van tweede-orde Adams–Bashforth (AB2) extrapolatie gebaseerd op gegevens van voorgaande tijdstappen ($t_n$ en $t_{n-1}$).
4. Algoritmische structuur: In elke tijdstap worden de Stokes- en Biot-subproblemen onafhankelijk en parallel opgelost. De kinematische relatie tussen de structuurverplaatsing en de snelheid wordt afgehandeld via een BDF2-relatie binnen het Biot-subprobleem.

Belangrijkste Bijdragen
De primaire bijdragen van het werk zijn theoretisch en analytisch van aard:

• Discrete Stabiliteitsanalyse: De auteurs leveren een rigoureus stabiliteitsbewijs voor het voorgestelde schema. Door gebruik te maken van BDF2-energie-identiteiten en een scherpe decompositie van de geëxtrapoleerde interfacetermen, stellen zij een gesloten stabiliteitsgrens vast. Deze grens houdt stand onder een parabolische Courant–Friedrichs–Lewy (CFL) conditie die de tijdstap $\Delta t$ en de roostergrootte $h$ koppelt (specifiek $\Delta t \leq c^* h^2$).
• A Priori Foutschattingen: Het artikel leidt een a priori foutschatting af met behulp van een projectiegebaseerd framework. Dit omvat:
  • Een Fortin-projectie voor de vloeistofsnelheid om de divergentievrije beperking te hanteren.
  • Ritz-type projecties voor de poroelastische verplaatsing en de poriedruk.
  • Een zorgvuldige analyse van consistentiefouten voortvloeiend uit de BDF2-tijddiscretisatie, de AB2-interfaceextrapolatie en de geprojecteerde kinematische relatie.
• Convergentiegraden: De analyse toont aan dat de totale fouten in de vloeistofsnelheid, de structuurverplaatsing, de poriedruk en de elastische verplaatsing begrensd worden door $C(h^k + \Delta t^2)$ in bulk-energienormen voor $1 \leq k \leq 3$. Dit bevestigt een tweede-orde temporele nauwkeurigheid en een optimale ruimtelijke convergentie.

Numerieke Resultaten
De theoretische bevindingen worden gevalideerd door twee numerieke experimenten:

1. Vast-domein Benchmark (Manufactured Solution): Een testgeval met bekende exacte oplossingen werd gebruikt om de convergentiegraden te verifiëren. De resultaten bevestigden de voorspelde tweede-orde convergentie in de tijd voor alle variabelen (vloeistofsnelheid/druk, structuurverplaatsing/snelheid, en poriedruk). Ruimtelijke convergentietests op verfijnde roosters toonden optimale convergentiegraden die consistent zijn met de gekozen eindige elementruimtes (Taylor–Hood $P_2/P_1$ voor Stokes en $P_2$-elementen voor Biot-variabelen).
2. Bewegend-domein Voorbeeld (Navier–Stokes–Biot): Om de toepasbaarheid van de gepartitioneerde strategie buiten de vaste-domein Stokes–Biot setting aan te tonen, hebben de auteurs het schema toegepast op een 2D bloedstroomprobleem in een bewegend domein, beheerst door de Navier–Stokes-vergelijkingen gekoppeld aan het Biot-systeem. Met behulp van een Arbitrary Lagrangian–Eulerian (ALE) formulering heeft de simulatie succesvol drukgolfpropagatie en structurele vervorming gevangen, waarbij de resultaten consistent zijn met eerdere literatuur.

Betekenis en Claims
Het artikel claimt een wiskundig gerechtvaardigde, tweede-orde nauwkeurige en paralleliseerbare expliciete gepartitioneerde methode voor vloeistof–poroelastische interactie te bieden. De significantie ligt in het overwinnen van de stabiliteitsuitdagingen die gewoonlijk geassocieerd worden met expliciete koppeling in FPSI-problemen. Door stabiliteit te bewijzen onder een parabolische CFL-conditie en rigoureuze foutgrenzen af te leiden, stelt het werk vast dat expliciete splitting zowel efficiënt als nauwkeurig kan zijn zonder wiskundige strengheid op te offeren.

De auteurs merken op dat hoewel de stabiliteits- en foutanalyse strikt is voor het vaste-domein Stokes–Biot-systeem, het numerieke voorbeeld suggereert dat de strategie robuust genoeg is om uitgebreid te worden naar meer algemene bewegende-domeinstromingsproblemen, zoals die met de Navier–Stokes-vergelijkingen. Zij concluderen dat het schema een praktische compromis biedt tussen modulariteit, parallel efficiëntie en bewijsbare nauwkeurigheid, waarbij zij bescheiden erkennen dat verdere arbeid nodig is om de CFL-restrictie te versoepelen of het framework uit te breiden naar geometrisch niet-lineaire bewegende-interfaceproblemen.