Monolithic Multi-level Overlapping Schwarz Solvers for Fluid Problems

Resumo Técnico: Solucionadores de Schwarz Sobrepostos Monolíticos Multinível para Problemas de Fluidos

Enunciado do Problema
O artigo aborda a simulação numérica da dinâmica de fluidos incompressíveis, especificamente a resolução das equações de Stokes discretizadas, que resultam em sistemas de ponto de sela em grande escala. À medida que a complexidade geométrica aumenta e as resoluções de malha tornam-se mais finas para resolver camadas limite, solucionadores diretos tornam-se inviáveis devido a custos computacionais proibitivos e consumo de memória. Os autores concentram-se nas limitações de escalabilidade dos métodos padrão de decomposição de domínio de dois níveis, especificamente na abordagem GDSW (Generalized Dryja–Smith–Widlund). Em um esquema de dois níveis, à medida que o número de subdomínios aumenta, o problema grosseiro ($A_0$) cresce em tamanho, tornando eventualmente sua fatoração direta computacionalmente intratável.

Metodologia
Os autores propõem e implementam um precondicionador de Schwarz sobreposto monolítico multinível para superar os limites de escalabilidade da abordagem GDSW de dois níveis.

• Formulação Matemática: O método estende o precondicionador GDSW de dois níveis aplicando-o recursivamente. Em vez de inverter diretamente o problema grosseiro $A_0$, os autores substituem $A_0^{-1}$ por uma aproximação derivada de um nível mais grosseiro (uma hierarquia de três níveis ou multinível). O precondicionador é definido como:
  $$B_{GDSW, 3-level} = \Phi^{(2)} \left( \Phi^{(3)} A_0^{(3)-1} \Phi^{(3)T} + \sum_{i=1}^{N^{(2)}} R_i^{(2)T} A_i^{(2)-1} R_i^{(2)} \right) \Phi^{(2)T} + \sum_{i=1}^{N^{(1)}} R_i^{(1)T} A_i^{(1)-1} R_i^{(1)}$$
  Esta estrutura permite que o problema grosseiro seja resolvido por meio de um nível adicional de decomposição de domínio, em vez de um solucionador direto.

• Implementação: O solucionador é implementado dentro da biblioteca FROSch (Fast and Robust Overlapping Schwarz), parte do pacote ShyLU no ecossistema Trilinos.

  • Acoplamento: Uma interface escalável conecta o FROSch à biblioteca de elementos finitos FEATFLOW. O FEATFLOW trata da montagem da matriz do sistema e dos vetores do lado direito, utilizando elementos $Q_2$ para a velocidade e elementos $P_1$ descontínuos para a pressão.
  • Natureza Algébrica: A implementação é puramente algébrica, exigindo apenas informações sobre a distribuição paralela dos graus de liberdade (dofs). Ela constrói automaticamente as matrizes locais de subdomínio, os operadores de restrição/interpolção e o problema grosseiro.
  • Construção do Espaço Grosseiro: Para a variável de pressão, que é descontínua através dos subdomínios, os autores utilizam funções constantes sobre todo o subdomínio como funções de base grosseiras, evitando a necessidade de uma interface de pressão $\Gamma_p$. Para a velocidade, são utilizadas extensões que minimizam a energia das restrições do espaço nulo.

Principais Contribuições

1. Extensão para Multinível: O trabalho estende com sucesso o precondicionador GDSW monolítico de dois níveis (anteriormente introduzido em [7]) para uma estrutura de três níveis (e potencialmente multinível) especificamente para problemas de fluidos de ponto de sela.
2. Acoplamento Escalável: O desenvolvimento de uma interface entre a biblioteca FEATFLOW e o framework Trilinos/FROSch, permitindo o acoplamento eficiente da montagem geométrica com solucionadores algébricos de decomposição de domínio.
3. Resultados Paralelos em Grande Escala: Os autores demonstram resultados de escalabilidade fraca em arquiteturas de computação de alto desempenho (HPC), alcançando escalabilidade até 32.768 ranks MPI para um escoamento de Poiseuille em um cubo unitário e até 4.096 ranks MPI para uma geometria complexa de matriz de extrusão.

Resultados Numéricos
Os autores avaliaram o solucionador usando dois casos de teste: um escoamento de Poiseuille em um cubo unitário e um escoamento de Poiseuille em uma geometria complexa de matriz de extrusão.

• Cubo Unitário (Decomposição Estruturada):

  • Solucionador de Dois Níveis: Em 32.768 subdomínios, o solucionador de dois níveis falhou porque o problema grosseiro (407.837 dofs) não pôde ser fatorado. O tempo de configuração para a fatoração do problema grosseiro tornou-se o gargalo.
  • Solucionador de Três Níveis: Resolveu com sucesso o problema em 32.768 subdomínios. Embora o número de iterações de Krylov tenha aumentado de 41 para 66 à medida que os subdomínios cresceram por um fator de 64, o solucionador permaneceu robusto. A abordagem de três níveis reduziu significativamente o tempo de configuração em comparação com a abordagem de dois níveis para grandes contagens de subdomínios, evitando a fatoração inviável do problema grosseiro massivo.

• Matriz de Extrusão (Decomposição Não Estruturada):

  • Solucionador de Dois Níveis: Em 4.096 subdomínios, o tempo de configuração foi dominado pela fatoração do problema grosseiro (84,3s).
  • Solucionador de Três Níveis: Demonstrou um ganho de desempenho. O tempo de configuração numérica diminuiu por um fator de aproximadamente 6 em comparação com o solucionador de dois níveis, e o tempo de solução diminuiu por quase um fator de 2. A abordagem de três níveis lidou com sucesso com a decomposição não estruturada onde a abordagem de dois níveis enfrentou limites de escalabilidade.

Significância e Afirmações
O artigo afirma que a principal significância deste trabalho é a demonstração bem-sucedida de que precondicionadores de Schwarz sobrepostos monolíticos multinível podem superar os gargalos de escalabilidade inerentes aos métodos GDSW de dois níveis para problemas de fluidos. Especificamente:

• A abordagem de três níveis permite a resolução de problemas com números significativamente maiores de subdomínios (até 32.768), onde o problema grosseiro de dois níveis seria, de outra forma, inviável de fatorar.
• O método fornece um ganho de desempenho tanto nos tempos de configuração quanto de solução para problemas em grande escala, particularmente quando o espaço grosseiro cresce.
• O trabalho representa um passo em direção a simulações eficientes de Dinâmica dos Fluidos Computacional (CFD) em futuras arquiteturas exascale com sistemas heterogêneos CPU-GPU, como parte do projeto alemão "StroemungsRaum".

Os autores permanecem modestos quanto ao aumento nas iterações de Krylov em contagens de subdomínios muito altas, observando que o precondicionador pode ainda não estar em sua faixa assintótica, mas enfatizam que o crescimento das iterações é aceitável dado o aumento massivo no paralelismo.