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

Resumo Técnico: Um Precondicionador Eficiente para Poromecânica de Contato de Dimensão Mista Baseado no Esquema de Divisão de Tensão Fixa (Fixed Stress Splitting)

Definição do Problema
A simulação numérica do contato poromecânico de fraturas é crítica para aplicações como sequestro de CO2, produção de energia geotérmica e armazenamento subterrâneo de gás. Esses problemas envolvem uma física fortemente acoplada: fluxo de fluidos, deformação da matriz (poromecânica de Biot) e mecânica de contato friccional dentro de fraturas. A modelagem precisa é desafiadora devido à geometria complexa das fraturas (frequentemente representadas como objetos de dimensão inferior inseridos em um domínio de dimensão superior) e à natureza não suave das condições de contato (estados de aderência, deslizamento e abertura, governados pelo atrito de Coulomb e restrições de não penetração).

Quando resolvidos utilizando passos de tempo totalmente implícitos e o método de Newton, esses problemas geram sistemas lineares totalmente acoplados e de grande escala. A eficiência e a robustez do solver dependem fortemente do precondicionador. Métodos existentes enfrentam obstáculos específicos:

1. Estrutura de Ponto de Selado (Saddle-Point): A formulação de Lagrange Ampliado da mecânica de contato introduz uma estrutura de ponto de selado com uma submatriz singular, o que complica a aplicação de técnicas de precondicionamento padrão.
2. Complexidade de Acoplamento: O esquema de divisão de tensão fixa (fixed stress splitting), um método bem estabelecido para o desacoplamento da poromecânica, é difícil de estender para meios fraturados com mecânica de contato. Solvers de ponto de selado padrão geralmente assumem estruturas de matriz específicas (ex: um operador do tipo Laplaciano) que não se sustentam quando o operador encapsula o problema poromecânico acoplado.
3. Escalabilidade: Os precondicionadores devem permanecer eficientes à medida que os tamanhos dos problemas crescem para incluir grandes redes de fraturas, típicas de simulações reais de reservatórios.

Metodologia
Os autores propõem um precondicionador robusto e escalável que combina aproximações de complemento de Schur aninhadas com uma transformação linear. A abordagem estende o esquema de tensão fixa tanto para os subdomínios da matriz quanto para os das fraturas, abordando simultaneamente a natureza singular da mecânica de contato.

A metodologia procede através das seguintes etapas:

1. Transformação Linear (Pré-processamento): Para lidar com a singularidade na submatriz de mecânica de contato (especificamente para estados de aderência da fratura), uma matriz de transformação linear à direita ($Q_r$) é aplicada ao Jacobiano. Essa transformação torna a submatriz de mecânica de contato invertível ao acoplá-la com as equações de balanço de forças de interface, inspirado em técnicas utilizadas em outros problemas de ponto de selado.
2. Decomposição de Complemento de Schur Aninhado: O precondicionador aproxima o fator triangular superior da decomposição LDU em blocos do Jacobiano transformado através de três níveis de eliminação de complemento de Schur:
   • Nível 1 (Mecânica de Contato): Elimina os graus de liberdade da mecânica de contato. A submatriz de contato transformada é invertida exatamente (ou via aproximação de bloco diagonal) para desacoplar as condições de contato das forças de interface.
   • Nível 2 (Balanço de Momento): Aplica uma aproximação baseada em tensão fixa para desacoplar as equações de balanço de momento. Isso envolve a construção de uma matriz de estabilização ($D_{55}$) com coeficientes específicos para a matriz ambiente ($l_{mat}$) e subdomínios de fratura ($l_{frac}$). As equações acopladas de momento e balanço de força de interface são resolvidas monotonicamente usando Multigrid Algébrico (AMG).
   • Nível 3 (Fluxo de Interface): Elimina o fluxo de fluido de interface aproximando a submatriz de fluxo de interface com sua diagonal (utiliza-se ILU(0) para o sistema de fluxo de interface). O complemento de Schur resultante, que representa o balanço de massa de fluido, é resolvido monotonicamente com AMG.
3. Análise Teórica: Os autores fornecem uma justificativa teórica para o esquema de divisão de tensão fixa neste contexto. Ao analisar um modelo linearizado simplificado, eles provam que o esquema iterativo sequencial é uma contração sob suposições específicas relativas às propriedades dos materiais, coeficientes de atrito e estabilidade do estado de contato. Esta análise valida a escolha dos coeficientes de estabilização usados no precondicionador.

Principais Contribuições

• Design de Novo Precondicionador: O artigo introduz um precondicionador que integra com sucesso o esquema de divisão de tensão fixa com a mecânica de contato em poromecânica de dimensão mista. Ele supera a singularidade da submatriz de contato via transformação linear e lida com o complexo acoplamento através de complementos de Schur aninhados.
• Prova de Convergência Teórica: Os autores derivam garantias de convergência para o esquema de divisão de tensão fixa aplicado à poromecânica de contato fraturada. Eles demonstram que o esquema é uma contração sob suposições fisicamente razoáveis, vinculando a taxa de convergência aos coeficientes de atrito e parâmetros de elasticidade de contato.
• Escalabilidade e Robustez: O método foi projetado para ser escalável, utilizando AMG para os subproblemas dominantes (elasticidade e balanço de massa de fluido) e ILU para fluxos de interface, tornando-o adequado para problemas de grande escala com redes de fraturas complexas.

Resultados
O precondicionador foi validado através de uma série de experimentos numéricos utilizando o framework PorePy:

• Sensibilidade de Parâmetros: Testes em um modelo de fratura única mostraram que o número de iterações depende primariamente do coeficiente de atrito ($F$), onde um atrito mais alto exige mais iterações, consistente com a análise teórica. A rigidez normal da fratura ($K_n$) teve pouco efeito na convergência.
• Refinamento de Malha (2D): Em um modelo 2D com uma única fratura, o número de iterações do GMRES permaneceu estável ao utilizar sub-solvers diretos. Ao substituir os solvers diretos por AMG, o número de iterações aumentou ligeiramente (por um fator de ~2) conforme os graus de liberdade aumentavam em três ordens de magnitude, demonstrando boa escalabilidade.
• Complexidade Geométrica (2D): Experimentos com múltiplas fraturas que se interceptam mostraram que o precondicionador manteve a eficiência, com apenas um pequeno aumento nas iterações do GMRES para os subsolvers AMG conforme a malha era refinada.
• Escalabilidade 3D: Uma simulação 3D com múltiplas fraturas e quase $10^6$ graus de liberdade demonstrou que o número de iterações do GMRES permaneceu quase constante com o refinamento da malha, confirmando a escalabilidade do método para problemas de grande escala. Não foram observadas falhas de solver linear no caso 3D.

Significância
O artigo afirma que o precondicionador proposto oferece uma solução robusta e escalável para o desafiador problema da poromecânica de contato friccional em meios fraturados. Ao estender o esquema de tensão fixa para lidar com a estrutura de ponto de selado da mecânica de contato via transformação linear e complementos de Schur aninhados, o método permite a simulação eficiente de processos subsuperficiais complexos. Os autores enfatizam que a abordagem é aplicável a qualquer discretização de (poro)mecânica de contato que utilize abordagens baseadas em multiplicadores de Lagrange, fornecendo um caminho para a modelagem precisa e computacionalmente viável de sismicidade induzida, transporte de fluidos e deformação em reservatórios fraturados.