A Stabilized Mortar Method for Discontinuities in Geological Media with Non-Conforming Grids

1. Problema e Contexto

A simulação numérica precisa da mecânica de falhas e fraturas é crucial para a avaliação de segurança e desempenho de sistemas subsuperficiais (ex.: armazenamento de gás, energia geotérmica, sequestro de carbono). O principal desafio reside na representação discretizada de superfícies de descontinuidade e na simulação robusta do seu comportamento de contato friccional, especialmente quando se utilizam malhas não conformes (onde os elementos dos dois lados da falha não coincidem geometricamente).

Embora métodos como o de Penalidade sejam simples, eles impõem as restrições apenas de forma aproximada. Métodos de Multiplicadores de Lagrange (LMM) oferecem maior precisão, mas introduzem desafios de estabilidade matemática. Especificamente, ao utilizar o Método de Argamassa (Mortar Method) para permitir malhas não conformes, a escolha de espaços funcionais para os deslocamentos e os multiplicadores de Lagrange deve satisfazer a condição de estabilidade inf-sup (ou condição de Babuška-Brezzi).

O artigo foca em uma configuração específica: deslocamentos lineares por partes (Q1) combinados com multiplicadores de Lagrange constantes por partes (P0) definidos apenas em um lado da interface (lado "não-argamassa"). Embora essa escolha de multiplicadores P0 traga vantagens práticas (como tratamento simplificado de condições de contorno e pontos de cruzamento), ela viola a condição de inf-sup, levando a oscilações espúrias nos perfis de tração, especialmente quando o lado não-argamassa é mais refinado que o lado argamassa.

2. Metodologia

Os autores propõem uma formulação de argamassa estabilizada para resolver o problema de contato com atrito (Lei de Coulomb) em meios geológicos.

• Formulação Variacional: O problema é formulado usando multiplicadores de Lagrange para impor as restrições de contato (não penetração e lei de atrito) de forma fraca.
• Espaços Discretos:
  • Deslocamentos: Elementos finitos trilineares (Q1) em hexaedros.
  • Multiplicadores de Lagrange: Funções constantes por partes (P0) definidas no lado não-argamassa.
• Estratégia de Estabilização: Para restaurar a estabilidade inf-sup sem alterar os espaços funcionais, os autores introduzem um termo de estabilização por salto de tração (traction-jump stabilization) nas equações de restrição.
  • O termo adicionado é uma forma bilinear que penaliza o salto da tração através das arestas internas da malha não-argamassa.
  • A formulação é: $R^*_t = R_t + s(\mu, t_h)$, onde $s$ é o termo de estabilização.
• Derivação do Fator de Escala (Sem Parâmetros de Usuário):
  • Um dos maiores desafios em métodos estabilizados é a escolha do coeficiente de estabilização. Os autores desenvolveram uma estratégia algorítmica baseada na teoria de macroelementos.
  • Eles calculam uma aproximação do complemento de Schur local para cada macroelemento (um conjunto de células ao redor de um nó interno).
  • O tensor de escala local $S_{EKL}$ é derivado automaticamente a partir da inversa da matriz de rigidez diagonalizada do macroelemento. Isso elimina a necessidade de parâmetros de ajuste selecionados pelo usuário, garantindo que a estabilização seja forte o suficiente para remover modos espúrios, mas não tão forte a ponto de degradar a precisão.
• Resolução Numérica: O sistema não linear resultante é resolvido utilizando uma estratégia de conjunto ativo (active-set) para determinar as regiões de aderência (stick), deslizamento (slip) e abertura (open), acoplado ao método de Newton-Raphson.

3. Contribuições Principais

1. Estabilização Inf-Sup para P0: Demonstração de que o uso de multiplicadores P0 (constant per element), que são computacionalmente eficientes e robustos para geometrias complexas, pode ser estabilizado para satisfazer a condição inf-sup através de um termo de salto de tração.
2. Algoritmo Automático de Escalonamento: Desenvolvimento de um método baseado em macroelementos para calcular automaticamente a magnitude da estabilização, removendo a subjetividade e a necessidade de calibração manual de parâmetros.
3. Robustez em Malhas Não Conformes: A metodologia demonstra capacidade de lidar com malhas industriais padrão, especificamente malhas de pontos de canto (corner-point grids), que são comuns em simulações de reservatórios e frequentemente geram interfaces não conformes complexas.
4. Superação de Oscilações: A abordagem supera as oscilações de tração que ocorrem em formulações estáveis teoricamente (como multiplicadores duais) quando o lado não-argamassa é mais refinado que o lado argamassa.

4. Resultados Numéricos

O método foi validado através de cinco testes numéricos:

• Teste de Patch (Test 1): Verificou a capacidade de transferir um campo de tração constante através de interfaces não conformes. A formulação estabilizada reproduziu a solução exata, enquanto multiplicadores nodais padrão apresentaram oscilações espúrias.
• Deslizamento Constante (Test 2): Reproduziu com exatidão a solução analítica para uma interface deslizando completamente sob atrito de Coulomb.
• Fratura Única sob Compressão (Test 3): Comparação com solução analítica para uma fratura em domínio infinito. Os resultados mostraram boa concordância, com o termo de estabilização ajudando a suavizar oscilações próximas à ponta da fratura.
• Dois Blocos sob Compressão e Cisalhamento (Test 4): Este foi o teste crítico. Em uma configuração onde o lado não-argamassa era duas vezes mais refinado que o lado argamassa ($r_h = 1/2$), formulações LMM teoricamente estáveis (multiplicadores duais e enriquecimento por bolhas) falharam, exibindo oscilações severas na tração normal. A formulação proposta com estabilização de salto de tração produziu um perfil de tração suave e preciso, correspondendo à solução de referência em malha fina.
• Retirada de Aquífero em Domínio Falhado (Test 5): Aplicação realista em um modelo geológico com malha de pontos de canto (corner-point). O modelo simulou a ativação de uma falha devido à extração de fluidos. O método demonstrou estabilidade e capacidade de simular a evolução suave do estado de contato (de aderência para deslizamento e abertura) em geometrias complexas típicas da indústria de petróleo e gás.

5. Significado e Conclusão

Este trabalho oferece uma solução robusta e prática para um problema clássico na mecânica de contato geológico: a incompatibilidade entre a necessidade de malhas não conformes (exigidas por geometrias complexas e fluxos de dados industriais) e a estabilidade matemática dos métodos de multiplicadores de Lagrange.

A principal inovação é a capacidade de utilizar multiplicadores constantes (P0), que simplificam a implementação de condições de contorno e tratam naturalmente pontos de interseção de falhas, sem sacrificar a estabilidade do sistema. O algoritmo de auto-escalonamento da estabilização torna o método "plug-and-play" para aplicações industriais, eliminando a dependência de parâmetros empíricos.

O método é particularmente significativo para a indústria de reservatórios, pois permite a integração direta de simulações geomecânicas complexas com malhas de fluxo padrão (corner-point), facilitando a avaliação de segurança em projetos de armazenamento de carbono, fraturamento hidráulico e extração de recursos.