Multigrid Preconditioning for FEEC using Mass-Lumping and Transforming Smoothers

Resumo Técnico: Précondicionamento de Multigrid para FEEC usando Mass-Lumping e Smoothers de Transformação

Definição do Problema
O artigo aborda o desafio de resolver sistemas lineares decorrentes de discretizações de Cálculo Exterior de Elementos Finitos (FEEC) de equações diferenciais parciais (EDPs) formuladas no complexo de de Rham. Especificamente, o trabalho visa formulações mistas e de ponto de sela para o operador de Hodge-Dirac, o Laplaciano de Hodge e sistemas de magnetostática (equações de Maxwell no gauge de Coulomb). Essas discretizações tipicamente resultam em sistemas lineares indefinidos com estruturas de ponto de sela. Smoothers iterativos padrão, como Gauss-Seidel ou Jacobi, são ineficazes para esses sistemas indefinidos porque os operadores carecem da definição positiva necessária para esses esquemas de relaxação. Embora trabalhos anteriores tenham proposto abordagens de multigrid baseadas em complementos de Schur, este artigo busca uma estratégia alternativa que evite a complexidade da inversão do complemento de Schur.

Metodologia
O framework proposto combina mass-lumping com smoothers de transformação para construir um précondicionador eficaz para os sistemas consistentes de FEEC.

1. Mass-Lumping e Equivalência Espectral:
   Os autores introduzem uma versão de mass-lumped dos operadores de FEEC onde as matrizes de massa consistentes são substituídas por matrizes diagonais (ou explicitamente invertíveis). Isso permite a inversão explícita do isomorfismo de Riesz, um passo que, de outra forma, seria proibitivamente caro com matrizes de massa consistentes densas.

   • Sob suposições suaves de equivalência de norma $h$-uniforme entre os produtos internos consistentes e de mass-lumped, o artigo prova uma desigualdade de Poincaré discreta para as normas de mass-lumped.
   • Este resultado de estabilidade leva a uma prova de equivalência espectral entre os operadores FEEC consistentes originais e os operadores de mass-lumped. Essa equivalência justifica o uso do sistema de mass-lumped como um proxy para o précondicionamento do sistema original.

2. Smoothers de Transformação:
   Para lidar com a indefinição dos operadores de mass-lumped, o método emprega smoothers de transformação. A ideia central é transformar algebricamente o operador em uma forma de bloco onde os blocos diagonais são definidos positivos, permitindo o uso de smoothers padrão como o Gauss-Seidel.

   • Para os problemas de Hodge-Dirac, Laplaciano de Hodge e magnetostática, matrizes de transformação à esquerda e à direita específicas ($T_L, T_R$) são construídas.
   • Essas transformações mapeiam o sistema indefinido original para um sistema triangular inferior de blocos (ou um sistema com blocos diagonais definidos positivos).
   • Uma relaxação de Gauss-Seidel é aplicada a esse sistema transformado, e o resultado é mapeado de volta para o espaço original.

3. Précondicionamento de Multigrid:
   O algoritmo de multigrid é projetado especificamente para o sistema de mass-lumped transformado. Este ciclo de multigrid ($M_{MG}$) serve como uma aproximação inversa para o operador de mass-lumped.

   • Para resolver o sistema FEEC original $G_h u = g$, o précondicionador $B$ é construído como $B \approx J_h^{-1} M_{MG} \tilde{H}_h (J_h^*)^{-1}$, onde $J_h$ é o isomorfismo entre os espaços padrão e de mass-lumped (frequentemente a identidade na prática) e $\tilde{H}_h$ é o mapa de Riesz de mass-lumped.
   • Como o sistema précondicionado resultante permanece indefinido, os autores utilizam o GMRES (método de Resíduos Mínimos Generalizados) como o solver externo, já que métodos padrão como PCG ou MINRES requerem précondicionadores definidos positivos.

Principais Contribuições

• Design Algorítmico: O artigo propõe um framework unificado de précondicionamento de multigrid que integra mass-lumping com smoothers de transformação para problemas de FEEC.
• Estabilidade Teórica: Estabelece a estabilidade dos sistemas de mass-lumped e prova a equivalência espectral com os operadores FEEC originais sob suposições suaves, fornecendo uma base teórica para a estratégia de précondicionamento.
• Tratamento de Indefinição: O trabalho demonstra como os smoothers de transformação podem converter operadores de bloco indefinidos em formas passíveis de técnicas de relaxação clássicas sem a necessidade de aproximações de complemento de Schur.
• Robustez em Topologias Complexas: O método mostra-se capaz de lidar efetivamente com domínios de topologias não triviais (formas harmônicas), onde a deficiência de posto do sistema é determinada pelos números de Betti.

Resultados Numéricos
Experimentos numéricos extensos foram conduzidos em 2D e 3D em várias malhas, incluindo aquelas com cantos reentrantes e buracos (topologia não trivial).

• Convergência: Os ciclos de multigrid (V-cycle e W-cycle) combinados com os smoothers de transformação demonstraram convergência robusta. Os raios espectrais dos operadores de propagação de erro permaneceram limitados longe da unidade, independentemente do tamanho da malha $h$.
• Estratégias de Mass-Lumping: Três estratégias de mass-lumping foram testadas: identidade escalonada, somas de linhas e duais baricêntricos (inspirados em Cálculo Exterior Discreto). Os resultados indicaram que essas diferentes estratégias de mass-lumping não levaram a diferenças qualitativas pronunciadas no comportamento de convergência.
• Desempenho do Solver: Quando usado como précondicionador para GMRES, o método resultou em contagens de iteração que foram amplamente independentes do nível de refinamento da malha. O W-cycle geralmente proporcionou uma convergência mais robusta do que o V-cycle.
• Formas Harmônicas: Os experimentos confirmaram que o método lida corretamente com a deficiência de posto associada às formas harmônicas, onde o número de modos autovalores não convergentes corresponde à dimensão do espaço harmônico.

Significância e Alegações
Os autores posicionam este trabalho primariamente como uma contribuição algorítmica, em vez de uma prova formal de teoria de convergência. Eles declaram explicitamente que, embora provem a equivalência espectral e a estabilidade, uma teoria de convergência rigorosa que estabeleça limites uniformes nas contagens de iteração do GMRES permanece como um problema aberto para investigação futura.

A significância da abordagem reside em sua capacidade de fornecer um précondicionador prático e robusto para sistemas FEEC indefinidos usando operações explícitas e baratas (mass-lumping) e transformações algébricas, evitando as complexidades dos métodos de complemento de Schur. A evidência numérica sugere que a abordagem é uma direção válida e promissora para resolver problemas de elementos finitos mistos em eletromagnetismo e áreas relacionadas, justificando novos estudos teóricos e extensões para outros problemas, como fluxo de Stokes ou sistemas com saltos de coeficientes.