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

Sintesi Tecnica: Un Precondizionatore Efficiente per la Poromeccanica di Contatto a Dimensioni Miste Basato sullo Schema di Fixed Stress Splitting

Definizione del Problema
La simulazione numerica della poromeccanica di contatto delle fratture è critica per applicazioni quali il sequestro di CO2, la produzione di energia geotermica e lo stoccaggio sotterraneo di gas. Questi problemi coinvolgono una fisica fortemente accoppiata: flusso di fluidi, deformazione della matrice (poromeccanica di Biot) e meccanica di contatto con attrito all'interno delle fratture. La modellazione accurata è resa difficile dalla geometria complessa delle fratture (spesso rappresentate come oggetti a dimensione inferiore immersi in un dominio a dimensione superiore) e dalla natura non regolare delle condizioni di contatto (stati di adesione, scorrimento e apertura governati dall'attrito di Coulomb e dai vincoli di non penetrazione).

Quando risolti utilizzando l'avanzamento temporale completamente implicito e il metodo di Newton, questi problemi generano grandi sistemi lineari completamente accoppiati. L'efficienza e la robustezza del solutore dipendono fortemente dal precondizionatore. I metodi esistenti affrontano ostacoli specifici:

1. Struttura a Punto di Sella: La formulazione dell'Augmented Lagrange della meccanica di contatto introduce una struttura a punto di sella con una sottomatrice singolare, il che complica l'applicazione delle tecniche di precondizionamento standard.
2. Complessità di Accoppiamento: Lo schema di fixed stress splitting, un metodo consolidato per il disaccoppiamento della poromeccanica, è difficile da estendere ai mezzi fratturati con meccanica di contatto. I solutori a punto di sella standard assumono spesso strutture matriciali specifiche (ad esempio, un operatore di tipo Laplaciano) che non si mantengono quando l'operatore incapsula il problema poromeccanico accoppiato.
3. Scalabilità: I precondizionatori devono rimanere efficienti man mano che le dimensioni del problema crescono per includere ampie reti di fratture tipiche delle simulazioni reali di serbatoi.

Metodologia
Gli autori propongono un precondizionatore robusto e scalabile che combina approssimazioni del complemento di Schur annidato con una trasformazione lineare. L'approccio estende lo schema di fixed stress sia ai subdomini della matrice che alle fratture, affrontando al contempo la natura singolare della meccanica di contatto.

La metodologia procede attraverso le seguenti fasi:

1. Trasformazione Lineare (Pre-processing): Per gestire la singolarità nella sottomatrice della meccanica di contatto (specificamente per gli stati di adesione della frattura), viene applicata una matrice di trasformazione lineare a destra ($Q_r$) al Jacobiano. Questa trasformazione rende invertibile la sottomatrice di contatto accoppiandola con le equazioni di bilancio delle forze di interfaccia, ispirandosi a tecniche utilizzate in altri problemi a punto di sella.
2. Decomposizione del Complemento di Schur Annidato: Il precondizionatore approssima il fattore triangolare superiore della decomposizione LDU del Jacobiano trasformato attraverso tre livelli di eliminazione del complemento di Schur:
   • Livello 1 (Meccanica di Contatto): Elimina i gradi di libertà della meccanica di contatto. La sottomatrice di contatto trasformata viene invertita esattamente (o tramite approssimazione a blocchi diagonali) per disaccoppiare le condizioni di contatto dalle forze di interfaccia.
   • Livello 2 (Bilancio del Momento): Applica un'approssimazione basata su fixed stress per disaccoppiare le equazioni del bilancio del momento. Ciò comporta la costruzione di una matrice di stabilizzazione ($D_{55}$) con coefficienti specifici per la matrice ambiente ($l_{mat}$) e i subdomini di frattura ($l_{frac}$). Le equazioni accoppiate del momento e del bilancio delle forze di interfaccia sono risolte monolithicamente utilizzando l'Algebraic Multigrid (AMG).
   • Livello 3 (Flusso di Interfaccia): Elimina il flusso di fluido di interfaccia approssimando la sottomatrice del flusso di interfaccia con la sua diagonale (viene utilizzato ILU(0) per il sistema di flusso di interfaccia). Il complemento di Schur risultante, che rappresenta il bilancio della massa del fluido, viene risolto monolithicamente con AMG.
3. Analisi Teorica: Gli autori forniscono una giustificazione teorica per lo schema di fixed stress splitting in questo contesto. Analizzando un modello linearizzato semplificato, dimostrano che lo schema iterativo sequenziale è una contrazione sotto specifiche assunzioni riguardanti le proprietà dei materiali, i coefficienti di attrito e la stabilità dello stato di contatto. Questa analisi valida la scelta dei coefficienti di stabilizzazione utilizzati nel precondizionatore.

Contributi Chiave

• Design di un Nuovo Precondizionatore: Il documento introduce un precondizionatore che integra con successo lo schema di fixed stress splitting con la meccanica di contatto nella poromeccanica a dimensioni miste. Esso supera la singolarità della sottomatrice di contatto tramite una trasformazione lineare e gestisce l'accoppiamento complesso attraverso complementi di Schur annidati.
• Prova di Convergenza Teorica: Gli autori derivano garanzie di convergenza per lo schema di fixed stress splitting applicato alla poromeccanica di contatto fratturata. Dimostrano che lo schema è una contrazione sotto ipotesi fisicamente ragionevoli, collegando il tasso di convergenza ai coefficienti di attrito e ai parametri di elasticità di contatto.
• Scalabilità e Robustezza: Il metodo è progettato per essere scalabile, utilizzando AMG per i subproblemi dominanti (elasticità e bilancio della massa del fluido) e ILU per i flussi di interfaccia, rendendolo adatto a problemi su larga scala con complesse reti di fratture.

Risultati
Il precondizionatore è stato validato attraverso una serie di esperimenti numerici utilizzando il framework PorePy:

• Sensibilità ai Parametri: I test su un modello a singola frattura hanno mostrato che il numero di iterazioni dipende principalmente dal coefficiente di attrito ($F$), dove un attrito più elevato richiede più iterazioni, in linea con l'analisi teorica. La rigidezza normale della frattura ($K_n$) ha avuto scarso effetto sulla convergenza.
• Raffinamento della Griglia (2D): In un modello 2D con una singola frattura, il numero di iterazioni GMRES è rimasto stabile utilizzando solutori diretti. Sostituendo i solutori diretti con AMG, il numero di iterazioni è aumentato solo leggermente (di un fattore di ~2) mentre i gradi di libertà aumentavano di tre ordini di grandezza, dimostrando una buona scalabilità.
• Complessità Geometrica (2D): Esperimenti con molteplici fratture intersecanti hanno mostrato che il precondizionatore mantiene l'efficienza, con solo un piccolo aumento delle iterazioni GMRES per i solutori AMG al raffinarsi della griglia.
• Scalabilità 3D: Una simulazione 3D con molteplici fratture e quasi $10^6$ gradi di libertà ha dimostrato che il conteggio delle iterazioni GMRES è rimasto quasi costante con il raffinamento della griglia, confermando la scalabilità del metodo per problemi su larga scala. Non sono stati osservati fallimenti del solutore lineare nel caso 3D.

Significatività
Il documento sostiene che il precondizionatore proposto offre una soluzione robusta e scalabile per la problematica complessa della poromeccanica di contatto con attrito in mezzi fratturati. Estendendo lo schema di fixed stress per gestire la struttura a punto di sella della meccanica di contatto tramite trasformazione lineare e complementi di Schur annidati, il metodo consente la simulazione efficiente di complessi processi sotterranei. Gli autori sottolineano che l'approccio è applicabile a qualsiasi discretizzazione della (poro)meccanica di contatto che utilizzi approcci basati su moltiplicatori di Lagrange, fornendo una via per la modellazione accurata e computazionalmente sostenibile della sismicità indotta, del trasporto di fluidi e della deformazione in serbatoi fratturati.