A Second-Order Monolithic Scheme for the Coupled Stokes--Biot Model Using Total Pressure

Technische Samenvatting: Een tweede-orde monolithisch schema voor het gekoppelde Stokes–Biot-model met behulp van totale druk

Probleemstelling
Het artikel behandelt de numerieke benadering van vloeistof–poroelastische structuurinteractie-systemen (FPSI), die de tijdafhankelijke Stokes-vergelijkingen (die een viskeuze incompressibele vloeistof beschrijven) koppelen aan de Biot-consolidatievergelijkingen (die een verzadigd poroelastisch medium beschrijven). Een primaire uitdaging bij het simuleren van deze systemen is het fenomeen van volumetrische vergrendeling (volumetric locking). In klassieke twee-veld verplaatsing–drukformuleringen van de Biot-vergelijkingen vertonen standaard eindige-elementen-discretisaties vaak een gebrek aan robuustheid wanneer de Lamé-parameter groot is (bijna incompressibele poreuze skeletten) of wanneer de permeabiliteit klein is. Dit leidt tot spookachtige drukoscillaties en kunstmatig stijve verplaatsingsresponsen die niet enkel door meshverfijning kunnen worden opgelost. Bovendien vertrouwen bestaande hogere-orde tijdsdiscretisatiemethoden voor gekoppelde Stokes–Biot-systemen vaak op operator-splitting of leggen zij restrictieve stabiliteitsvoorwaarden op die afhankelijk zijn van fysische parameters (bijv. de Beavers–Joseph–Saffman wrijvingsparameter of de Lamé-parameter), wat hun toepasbaarheid in uitdagende regimes beperkt.

Methodologie
De auteurs stellen een volledig impliciet, monolithisch, tweede-orde tijdsdiscretisatie voor met behulp van de Backward Differentiation Formula van orde 2 (BDF2). De kern van de methodologie bestaat uit het herformuleren van het Biot-subsysteem met behulp van een drie-veld totale-drukformulering. In deze aanpak wordt de totale druk ($\xi = \alpha p_p - \lambda_p \nabla \cdot \eta$) geïntroduceerd als een extra primaire onbekende naast de vaste verplaatsing ($\eta$) en de poriedruk ($p_p$).

Belangrijke methodologische kenmerken zijn onder meer:

• Monolithische Koppeling: Alle variabelen (vloeistofsnelheid, vloeistofdruk, vaste verplaatsing, totale druk en poriedruk) en alle interfacecondities (massa-behoud, spanningsbalans en Beavers–Joseph–Saffman slip) worden in elk tijdstap tegelijkertijd opgelost als een enkel lineair systeem. Dit elimineert operator-splitting-fouten en interface-extrapolatie.
• Totale-drukvariabele: Het gebruik van de totale druk creëert een parameter-robuuste saddle-point structuur die volumetrische vergrendeling onderdrukt op zowel continue als discrete niveaus, met name in bijna incompressibele regimes.
• Discrete Energieanalyse: Stabiliteit wordt vastgesteld met behulp van de BDF2 G-stabiliteitsidentiteit, een discrete energie-argument dat geen behoefte heeft aan CFL-type restricties of parameter-afhankelijke tijdstaplimieten.
• Eindige-elementen-discretisatie: De ruimtelijke discretisatie maakt gebruik van een gemengd variationeel kader gebaseerd op discrete inf–sup stabiliteit, waarbij Taylor–Hood-paren worden gebruikt voor de vloeistof- en poroelastische regio's om stabiliteit te garanderen.

Belangrijkste Bijdragen
Het artikel levert de volgende specifieke bijdragen:

1. Onvoorwaardelijke Energie-stabiliteit: De auteurs bewijzen dat het monolithische BDF2-schema onvoorwaardelijk stabiel is voor willekeurige tijdstapgroottes. De stabiliteitsschatting berust op de G-stabiliteitsidentiteit en is geldig zonder restricties die afhankelijk zijn van de Lamé-parameter, de Biot–Willis-coëfficiënt of de opslagcoëfficiënt.
2. Parameter-robuustheid: De analyse toont robuustheid aan ten opzichte van kritische fysische parameters. De schattingen blijven geldig wanneer de Lamé-parameter $\lambda_p \to \infty$ (bijna incompressibel limiet) en de opslagcoëfficiënt $s_0 \to 0$. Dit wordt bereikt door de totale-drukvariabele, die ervoor zorgt dat de energie uniform begrensd blijft.
3. Optimale Convergentiegraden: De auteurs leiden a priori foutschattingen af die tweede-orde convergentie in de tijd en optimale-orde convergentie in de ruimte ($O(\Delta t^2 + h^k)$) laten zien in de bijbehorende energienormen. De convergentieconstanten zijn onafhankelijk van de fysische parameters in de gebruikte gewogen normen.
4. Afhandeling van de Differentiaal-Algebraïsche Structuur: Het schema is geformuleerd om de differentiaal-algebraïsche aard van het gekoppelde systeem te behandelen, waarbij de vloeistofdruk wordt bepaald door incompressibiliteit en de totale druk via een gedifferentieerde opslagterm binnenkomt, zonder dat hiervoor aparte algebraïsche solvers of splitting nodig zijn.

Resultaten
Numerieke experimenten werden uitgevoerd om de theoretische bevindingen te verifiëren:

• Temporele Accuratesse: Tests bevestigden tweede-orde convergentie in de tijd voor alle variabelen (snelheid, verplaatsing, drukken) bij het gebruik van een Backward Euler opstartstap gevolgd door BDF2.
• Ruimtelijke Accuratesse: Ruimtelijke convergentiestudies toonden optimale snelheden die consistent zijn met de polynoomgraden van de eindige-elementenruimtes.
• Parameter-robuustheidstests:
  • Lamé-parameter ($\lambda_p$): De methode behield stabiliteit en nauwkeurigheid naarmate $\lambda_p$ met ordes van grootte toenam (benadering van het incompressibele limiet), wat het vergrendelingsvrije karakter van de formulering bevestigt.
  • Opslagcoëfficiënt ($s_0$): Het schema bleef robuust in het singuliere limiet $s_0 = 0$.
  • Vloeistofviscositeit ($\mu_f$): Variaties in viscositeit beïnvloedden de foutmarges, maar veranderden de asymptotische convergentiesnelheden niet.
  • Permeabiliteit ($\kappa$): Hoewel de methode robuust was voor matige permeabiliteit, observeerden de auteurs een verslechtering in de convergentiesnelheid van de poriedruk in het extreme regime van lage permeabiliteit ($\kappa = 10^{-6}$). Dit wordt toegeschreven aan de gesingularde perturbeerde aard van de Darcy-operator en de aanwezigheid van grenslagen, wat consistent is met de theoretische stabiliteitsschattingen die gewogen zijn door permeabiliteit.

Betekenis en Claims
Het artikel claimt dat het voorgestelde monolithische BDF2-schema met totale-drukformulering een stabiel, accuraat en betrouwbaar kader biedt voor het simuleren van gekoppelde Stokes–Biot-problemen over een breed scala aan fysisch relevante parameterregimes. Door operator-splitting te vermijden, elimineert de methode splitting-fouten, en door de totale-drukvariabele te gebruiken, overwint het de vergrendelingsproblemen die klassieke formuleringen in bijna incompressibele scenario's teisteren. De auteurs benadrukken dat hun stabiliteits- en convergentieanalyse parameter-robuust is, wat betekent dat de theoretische garanties niet verslechteren wanneer sleutelparameters (zoals stijfheid of compressibiliteit) extreem worden, een significante verbetering ten opzichte van eerdere methoden die restrictieve tijdstapcondities vereisten die afhankelijk waren van deze parameters. Het werk legt een rigoureuze basis voor hoogwaardige, volledig gekoppelde simulaties van vloeistof–poroelastische interacties.