Linear poroelasticity with solid incompressibility: consistent formulation and scalable numerical solution

Resumo Técnico: Poroelasticidade Linear com Incompressibilidade do Sólido

Enunciado do Problema
O artigo aborda a modelagem numérica de meios poroelásticos onde a fase sólida é incompressível, uma condição relevante para tecidos moles e hidrogéis. Embora as equações de poroelasticidade totalmente não lineares sejam gerais, suas formulações primais frequentemente sofrem de problemas de inconsistência quando a pressão depende da deformação, particularmente ao impor a incompressibilidade do sólido. Abordagens padrão frequentemente dependem de incompressibilidade total ($J=1$) ou quase-incompressibilidade, mas os autores visam derivar rigorosamente um modelo onde apenas a fase sólida é incompressível, evitando as inconsistências encontradas em formulações primais anteriores de poroelasticidade não linear.

Metodologia
Os autores propõem uma linearização consistente das equações de poroelasticidade totalmente não lineares sob a hipótese de incompressibilidade do sólido. A metodologia procede através dos seguintes passos:

1. Derivação do Modelo: Partindo de um modelo não linear com porosidade lagrangiana e um potencial de Helmholtz, os autores impõem a incompressibilidade do sólido via um multiplicador de Lagrange. Para evitar inconsistências decorrentes de termos de terceira ordem e a necessidade de computar o Laplaciano de um multiplicador $L^2$, eles adotam uma formulação de pressão total ($\mu$) e linearizam o sistema em torno de um ponto de equilíbrio.
2. Sistema Linearizado: O modelo resultante é um sistema de quatro campos envolvendo:
   • Deslocamento ($d$)
   • Pressão total do poro ($\mu$)
   • Porosidade ($\phi$)
   • Multiplicador de Lagrange da incompressibilidade ($\lambda$)
     O sistema é formulado como um problema fraco semi-discreto usando uma discretização temporal implícita, levando a uma estrutura de duplo ponto de sela (double saddle-point).
3. Análise Matemática: Os autores estabelecem a bem-postura (well-posedness) dos esquemas semi-discretos e totalmente discretos. Eles verificam a limitação dos operadores e provam a condição inf-sup (condições de Ladyzhenskaya–Babuška–Brezzi). Uma descoberta fundamental é que a condição inf-sup discreta exige apenas que o espaço de teste para o multiplicador de Lagrange ($V_{\lambda, h}$) seja um subespaço do espaço de porosidade ($V_{\phi, h}$). Isso permite o uso de pares de Elementos Finitos (FE) de ordem mais baixa (ex: $P_1$-contínuo para deslocamento e $P_0$-descontínuo para porosidade) sem exigir estabilização complexa ou elementos de ordem superior.
4. Equivalência ao Modelo de Biot: O artigo demonstra que este sistema de quatro campos proposto é matematicamente equivalente a um modelo de Biot padrão de dois campos no limite de incompressibilidade total, com um parâmetro de Lamé volumétrico modificado.
5. Estratégia de Solução Numérica: Para resolver o sistema linearizado de forma eficiente, os autores propõem um precondicionador de blocos escalável baseado no formalismo do complemento de Schur. Eles utilizam o método iterativo de tensão fixa (fixed-stress), interpretado como uma iteração de Uzawa, para aproximar o complemento de Schur. O precondicionador utiliza uma estratégia de partição de blocos:
   • O bloco de deslocamento é tratado via Multigrid Algébrico (AMG).
   • Os blocos restantes (pressão, porosidade, multiplicador) utilizam uma abordagem de Jacobi de blocos com AMG para o operador elíptico de pressão e Jacobi para matrizes de massa.

Principais Contribuições

• Formulação Consistente: Um modelo linearizado que impõe rigorosamente a incompressibilidade do sólido sem os problemas de inconsistência encontrados nas formulações primais padrão de poroelasticidade não linear.
• Estabilidade de Ordem Mais Baixa: A formulação admite uma aproximação FE estável de ordem mais baixa (pares do tipo $P_1-P_0$), que é computacionalmente eficiente e evita a necessidade de multiplicadores de interface ou termos de estabilização complexos.
• Solucionador Escalável: O desenvolvimento de um precondicionador robusto e escalável baseado na aproximação de tensão fixa que mantém contagens de iteração ótimas independentemente do tamanho da malha ou dimensão do problema (testado em 2D e 3D).
• Validação Teórica: Prova de bem-postura e taxas de convergência ($O(h^{k+1})$) para a discretização proposta.

Resultados
Testes numéricos foram conduzidos para validar as afirmações teóricas:

• Teste de Equivalência: Uma comparação do problema de Mandel em estado estacionário mostrou que o modelo de quatro campos proposto produz resultados indistinguíveis do modelo de Biot incompressível padrão (diferenças da ordem da precisão de máquina).
• Testes de Convergência: Testes de soluções manufaturadas em 2D e 3D confirmaram as taxas de convergência teóricas para deslocamento, pressão, porosidade e o multiplicador de Lagrange usando tanto esquemas de primeira ordem ($k=0$) quanto de segunda ordem ($k=1$).
• Testes de Escalabilidade: O desempenho do precondicionador foi avaliado em um benchmark de inchamento (swelling) 3D. O número médio de iterações do GMRES permaneceu constante (aproximadamente 60–70 iterações) conforme os Graus de Liberdade (DoFs) aumentaram de ~1.000 para mais de 33 milhões, demonstrando a escalabilidade do solucionador.

Significância
O artigo afirma que este trabalho fornece um framework consistente e computacionalmente eficiente para modelar fases sólidas incompressíveis em poroelasticidade. Ao derivar o modelo diretamente de equações não lineares e utilizar uma variável de pressão total, os autores unem a física não linear complexa e solucionadores lineares estáveis e práticos. A capacidade de usar elementos de ordem mais baixa reduz o custo computacional, enquanto o precondicionador baseado em tensão fixa garante que a estratégia de solução permaneça escalável para problemas 3D de grande escala. Os autores observam que, embora o modelo seja atualmente linearizado, o framework é projetado para ser consistente com a física não linear subjacente, e trabalhos futuros explorarão a aplicação desta estratégia para sistemas totalmente não lineares e acoplados (ex: termo-poroelasticidade).