A Jacobian-free Newton-Krylov method for cell-centred finite volume solid mechanics

Resumo Técnico: Um Método Newton-Krylov Livre de Jacobiano para Mecânica dos Sólidos de Volumes Finitos Centrados em Células

Definição do Problema
Os métodos tradicionais de volumes finitos (FV) para mecânica dos sólidos baseiam-se predominantemente em abordagens implícitas segregadas (iterações de Picard) ou esquemas totalmente explícitos. Embora os métodos segregados ofereçam simplicidade de implementação e eficiência de memória, eles frequentemente sofrem de convergência lenta, particularmente em problemas envolvendo não linearidade geométrica, grandes rotações ou acoplamentos complexos. Por outro lado, os métodos do tipo Newton-Raphson, comuns na comunidade de elementos finitos, oferecem convergência quadrática, mas exigem a montagem, armazenamento e inversão explícitos da matriz Jacobiana completa. Esse requisito é computacionalmente caro e intensivo em termos de memória, e adaptar as estruturas de FV segregadas existentes para métodos Newton de blocos acoplados não é trivial devido à complexidade de derivar e armazenar o verdadeiro Jacobiano. Consequentemente, há uma necessidade de um solver de blocos acoplados robusto e eficiente que possa ser integrado em estruturas de FV segregadas existentes sem uma refatoração de código significativa.

Metodologia
Os autores propõem uma abordagem Newton-Krylov Livre de Jacobiano (JFNK) adaptada para a mecânica dos sólidos de volumes finitos centrados em células. A metodologia central evita a formação explícita da matriz Jacobiana ($J$) ao aproximar sua ação sobre um vetor ($v$) usando diferenças finitas:
$$Jv \approx \frac{R(u + \epsilon v) - R(u)}{\epsilon}$$
onde $R$ é a função resíduo, $u$ é o vetor solução e $\epsilon$ é um escalar de perturbação. Esta aproximação é combinada com o solver de subespaço de Krylov de Resíduos Mínimos Generalizados (GMRES) para resolver os sistemas lineares que surgem dentro das iterações de Newton.

Para garantir eficiência e robustez, o método utiliza uma estratégia de pré-condicionamento. Em vez de usar o Jacobiano exato (que não está disponível), os autores utilizam um Jacobiano aproximado de estêncil compacto derivado do termo de inércia e de um termo de difusão (baseado na estabilização de Rhie-Chow usada em solvers segregados). Este pré-condicionador ($P$) é conceitualmente semelhante à matriz usada em abordagens segregadas tradicionais, permitindo que o novo método aproveite as rotinas de montagem e armazenamento de matriz existentes em frameworks como o toolbox sololids4foam do OpenFOAM.

A implementação utiliza a biblioteca PETSc para o solver não linear (SNES) e álgebra linear, com interface para o solver de FV solids4foam. Características algorítmicas principais incluem:

• Pré-condicionadores: Comparação de solvers diretos LU (MUMPS), ILU (Incomplete LU) e Multigrid Algébrico (HYPARE BoomerAMG).
• Globalização: Uso de busca de linha (line search) e continuação transiente (preditor-corretor) para direcionar a solução para o domínio de convergência quadrática.
• Estabilização: Um termo do tipo Rhie-Chow é adicionado ao resíduo para suprimir modos de energia zero (efeito checkerboard), com um fator de escala global ($\alpha$) ajustado para a velocidade de convergência.

Principais Contribuições

1. Primeira Aplicação à Mecânica dos Sólidos em FV: Este trabalho representa a primeira proposta de um método JFNK especificamente para a mecânica dos sólidos de volumes finitos, preenchendo uma lacuna entre aplicações de dinâmica de fluidos e mecânica dos sólidos.
2. Integração Fluida: Ao utilizar um Jacobiano aproximado de estêncil compacto como pré-condicionador, o método permite a adoção de convergência do tipo Newton em estruturas de FV segregadas existentes com modificações mínimas de código.
3. Benchmarking Abrangente: O método é rigorosamente testado contra uma abordagem segregada padrão através de uma série de casos variando em dimensionalidade (2-D/3-D), comportamento de material (elasticidade linear, hiperelasticidade, elastoplasticidade) e física (estática, dinâmica, interação fluido-estrutura).

Resultados
O estudo compara o método JFNK com a abordagem segregada utilizando métricas de tempo computacional, uso de memória e robustez:

• Eficiência: O método JFNK supera a abordagem segregada em todos os casos elásticos lineares e não lineares. Acelerações de uma a duas ordens de magnitude foram observadas em muitos casos, com uma aceleração máxima de 341x no caso de membrana de Cook elástica linear em 2-D.
• Memória: Para problemas com menos de $10^4$ graus de liberdade, o uso de memória é comparável ao do método segregado. Para problemas maiores, o uso de memória do JFNK aumenta, atingindo um máximo de 3,6 vezes o do método segregado (no caso da cavidade esférica), principalmente devido aos requisitos de armazenamento do solver GMRES e do pré-condicionador de multigrid.
• Robustez e Limitações:
  • O método JFNK demonstrou robustez superior em casos hiperelásticos (grandes deformações/rotações) onde o método segregado falhou em convergir.
  • No entanto, o método JFNK divergiu em casos elastoplásticos (membrana de Cook com plasticidade). Os autores atribuem isso ao pré-condicionador de estêncil compacto, que é baseado em elasticidade de pequenas deformações, sendo uma aproximação pobre do Jacobiano para grandes deformações plásticas.
• Escolha do Pré-condicionador:
  • Solvers Diretos LU foram mais rápidos para problemas de tamanho pequeno a moderado e casos 2-D.
  • Multigrid (AMG) foi mais eficaz para problemas 3-D maiores, oferecendo melhor escalabilidade e menor crescimento de memória do que o LU.
  • ILU(k) foi geralmente o mais lento e menos robusto.
• Estabilização: O aumento do fator de estabilização de Rhie-Chow ($\alpha$) geralmente melhorou a velocidade e a robustez da convergência, embora ao custo da precisão pré-assintótica.

Significância e Alegações
O artigo afirma que os métodos Newton-Krylov Livres de Jacobiano são um avanço promissor para a mecânica dos sólidos de volumes finitos, particularmente para estruturas segregadas existentes onde modificações mínimas permitem sua adoção. O estudo demonstra que o JFNK pode alcançar acelerações significativas e melhor robustez para problemas elásticos sem o custo proibitivo de montagem do Jacobiano completo.

Os autores mantêm-se modestos quanto às limitações atuais do método, notando explicitamente que a abordagem ainda não é robusta para elastoplasticidade com a atual estratégia de pré-condicionamento. Eles identificam o desenvolvimento de pré-condicionadores aprimorados para plasticidade e a inclusão de correções não ortogonais como trabalhos futuros necessários. A implementação é disponibilizada publicamente no toolbox solids4foam para incentivar a exploração e extensão pela comunidade.