An efficient preconditioner for mixed-dimensional contact poromechanics based on the fixed stress splitting scheme

Technische Samenvatting: Een Efficiënte Preconditioner voor Mixed-Dimensionale Contact-Poromechanica Gebaseerd op het Fixed Stress Splitting-Schema

Probleemstelling
De numerieke simulatie van breukcontact-poromechanica is cruciaal voor toepassingen zoals CO2-sequestratie, geothermische energieproductie en ondergrondse gasopslag. Deze problemen omvatten sterk gekoppelde fysica: vloeistofstroming, matrixdeformatie (Biot-poromechanica) en wrijvingscontactmechanica binnen breuken. Nauwkeurige modellering is uitdagend vanwege de complexe geometrie van breuken (vaak gerepresenteerd als lager-dimensionale objecten ingebed in een hoger-dimensionaal domein) en de niet-gladde aard van contactcondities (stick, slide en open staten beheerst door Coulomb-wrijving en niet-penetratiebeperkingen).

Wanneer deze problemen worden opgelost met volledig impliciete tijdstapmethoden en de Newton-methode, genereren ze grote, volledig gekoppelde lineaire systemen. De efficiëntie en robuustheid van de solver hangen zwaar af van de preconditioner. Bestaande methoden worden geconfronteerd met specifieke hindernissen:

1. Saddle-Point Structuur: De Augmented Lagrange-formulering van contactmechanica introduceert een saddle-point structuur met een singuliere submatrix, wat de toepassing van standaard preconditioning-technieken bemoeilijkt.
2. Koppelingscomplexiteit: Het fixed stress splitting-schema, een goed gevestigde methode voor het ontkoppelen van poromechanica, is moeilijk uit te breiden naar gefractureerde media met contactmechanica. Standaard saddle-point solvers gaan vaak uit van specifieke matrixstructuren (bijv. een Laplacian-achtige operator) die niet standhouden wanneer de operator de gekoppelde poromechanische problematiek omvat.
3. Schaalbaarheid: Preconditioners moeten efficiënt blijven naarmate de probleemomvang groeit om grote breuknetwerken te bevatten die typerend zijn voor real-world reservoir-simulaties.

Methodologie
De auteurs stellen een robuuste en schaalbare preconditioner voor die geneste Schur-complement benaderingen combineert met een lineaire transformatie. De aanpak breidt het fixed stress schema uit naar zowel matrix- als breuksubdomeinen en adresseert tegelijkertijd de singuliere aard van contactmechanica.

De methodologie verloopt via de volgende stadia:

1. Lineaire Transformatie (Preprocessing): Om de singulariteit in de contactmechanica-submatrix (specifiek voor klevende breuktoestanden) aan te pakken, wordt een recht lineaire transformatiematrix ($Q_r$) toegepast op de Jacobiaan. Deze transformatie maakt de contactmechanica-submatrix inverseerbaar door deze te koppelen aan de interface krachtbalansvergelijkingen, geïnspireerd door technieken gebruikt in andere saddle-point problemen.
2. Geneste Schur-complement Decompositie: De preconditioner benadert de bovenste driehoekige factor van de blok LDU-decompositie van de getransformeerde Jacobiaan via drie niveaus van Schur-complement eliminatie:
   • Niveau 1 (Contactmechanica): Elimineert de vrijheidsgraden van de contactmechanica. De getransformeerde contactsubmatrix wordt exact geïnverteerd (of via een blok-diagonale benadering) om contactcondities te ontkoppelen van interfacekrachten.
   • Niveau 2 (Momentum Balans): Past een fixed stress-gebaseerde benadering toe om de momentumbalansvergelijkingen te ontkoppelen. Dit omvat het construeren van een stabilisatiematrix ($D_{55}$) met specifieke coëfficiënten voor de ambient matrix ($l_{mat}$) en de breuksubdomeinen ($l_{frac}$). De gekoppelde momentum- en interface krachtbalansvergelijkingen worden monolithisch opgelost met behulp van Algebraic Multigrid (AMG).
   3. Niveau 3 (Interface Flow): Elimineert de interface vloeistofstroom door de interface flow submatrix te benaderen met de diagonaal (ILU(0) wordt gebruikt voor het interface flow systeem). Het resulterende Schur-complement, dat de interface massa-balans representeert, wordt monolithisch opgelost met AMG.
3. Theoretische Analyse: De auteurs bieden een theoretische rechtvaardiging voor het fixed stress splitting schema in deze context. Door een vereenvoudigd lineair model te analyseren, bewijzen zij dat het sequentiële iteratieve schema een contractie is onder specifieke aannames met betrekking tot materiaaleigenschappen, wrijvingscoëfficiënten en de stabiliteit van de contacttoestand. Deze analyse valideert de keuze van de stabilisatiecoëfficiënten die in de preconditioner worden gebruikt.

Belangrijkste Bijdragen

• Nieuwe Preconditioner Ontwerp: Het artikel introduceert een preconditioner die er succesvol in slaagt de fixed stress splitting methode te integreren met contactmechanica in mixed-dimensionale poromechanica. Het overwint de singulariteit van de contactsubmatrix via een lineaire transformatie en behandelt de complexe koppeling via geneste Schur-complementen.
• Theoretisch Convergentiebewijs: De auteurs leiden convergentiegaranties af voor de toepassing van het fixed stress splitting schema op gefractureerde contact-poromechanica. Zij demonstreren dat het schema een contractie is onder fysiek redelijke aannames, waarbij de convergentiesnelheid wordt gekoppeld aan wrijvingscoëfficiënten en contactelasticiteitsparameters.
• Schaalbaarheid en Robuustheid: De methode is ontworpen om schaalbaar te zijn, waarbij gebruik wordt gemaakt van AMG voor de dominante subproblemen (elasticiteit en vloeistofmassa-balans) en ILU voor interface flows, wat het geschikt maakt voor grootschalige problemen met complexe breuknetwerken.

Resultaten
De preconditioner werd gevalideerd via een reeks numerieke experimenten met het PorePy framework:

• Parametergevoeligheid: Tests op een enkelvoudig breukmodel toonden aan dat het aantal iteraties primair afhangt van de wrijvingscoëfficiënt ($F$), waarbij hogere wrijving meer iteraties vereist, wat consistent is met de theoretische analyse. De normale stijfheid van de breuk ($K_n$) had weinig effect op de convergentie.
• Grid Verfijning (2D): In een 2D-model met een enkele breuk bleef het aantal GMRES-iteraties stabiel bij het gebruik van directe subsolvers. Bij het vervangen van directe solvers door AMG, nam het aantal iteraties slechts licht toe (met een factor ~2) terwijl het aantal vrijheidsgraden met drie ordes van grootte toenam, wat een goede schaalbaarheid aantoont.
• Geometrische Complexiteit (2D): Experimenten met meerdere snijdende breuken toonden aan dat de preconditioner efficiënt bleef, met slechts een kleine toename in GMRES-iteraties voor de AMG subsolvers naarmate het grid werd verfijnd.
• 3D Schaalbaarheid: Een 3D-simulatie met meerdere breuken en bijna $10^6$ vrijheidsgraden toonde aan dat het aantal GMRES-iteraties nagenoeg constant bleef bij gridverfijning, wat de schaalbaarheid van de methode voor grootschalige problemen bevestigt. Er werden geen breakdowns van de lineaire solver waargenomen in de 3D-casus.

Betekenis
Het artikel stelt dat de voorgestelde preconditioner een robuuste en schaalbare oplossing biedt voor het uitdagende probleem van wrijvingscontact-poromechanica in gefractureerde media. Door het fixed stress schema uit te breiden om de saddle-point structuur van contactmechanica te behandelen via lineaire transformatie en geneste Schur-complementen, maakt de methode efficiënte simulatie van complexe ondergrondse processen mogelijk. De auteurs benadrukken dat de aanpak toepasbaar is op elke discretisatie van contact-(poro)mechanica die gebruikmaakt van Lagrange-multiplier-gebaseerde benaderingen, wat een pad biedt voor nauwkeurige en computationeel haalbare modellering van geïnduceerde seismiciteit, vloeistoftransport en deformatie in gefractureerde reservoirs.