Steady Incremental Viscosity Splitting Method for solving the stationary Navier-Stokes equation

Resumo Técnico: Método de Divisão Incremental de Viscosidade Estacionária para Equações de Navier-Stokes Estacionárias

Formulação do Problema
O artigo aborda a solução numérica das equações de Navier-Stokes estacionárias e incompressíveis com estabilização grad-div. As equações governantes são:
$$-\nu\Delta \mathbf{u} + \mathbf{u} \cdot \nabla \mathbf{u} + \nabla p - \gamma\nabla(\nabla \cdot \mathbf{u}) = \mathbf{f}, \quad \nabla \cdot \mathbf{u} = 0$$
em um domínio $\Omega$, com condições de contorno de Dirichlet homogêneas. Os autores observam que a resolução dos sistemas lineares de ponto de sela não simétricos resultantes (decorrentes da linearização padrão de Picard) é computacionalmente dispendiosa, particularmente para problemas 3D, onde solutores diretos falham devido a restrições de memória e solutores iterativos lutam com a estrutura de ponto de sela. Embora métodos de divisão estejam bem estabelecidos para fluxos não estacionários, sua aplicação a problemas genuinamente estacionários é menos desenvolvida.

Metodologia
Os autores propõem um Método de Divisão Incremental de Viscosidade Estacionária (SIVS). Esta abordagem adapta o esquema de Divisão Incremental de Viscosidade (IVS), anteriormente utilizado para fluxos não estacionários, ao contexto de estado estacionário. O método opera como um esquema iterativo de ponto fixo:

1. Etapa 1 (Previsão de Velocidade): Dada a iteração anterior $(\mathbf{u}^{k-1}, p^{k-1})$, resolve-se uma EDP elíptica para uma velocidade intermediária $\tilde{\mathbf{u}}^k$:
   $$a(\tilde{\mathbf{u}}^k, \mathbf{v}) + c^*(\mathbf{u}^{k-1}, \tilde{\mathbf{u}}^k, \mathbf{v}) + b(p^{k-1}, \mathbf{v}) = \langle \mathbf{f}, \mathbf{v} \rangle$$
   Aqui, a forma bilinear $a(\cdot, \cdot)$ inclui o termo viscoso e o termo de estabilização grad-div ($\gamma$).

2. Etapa 2 (Correção de Pressão e Projeção): Resolve-se para a velocidade atualizada $\mathbf{u}^k$ e a pressão $p^k$, corrigindo a solução intermediária:
   $$a(\mathbf{u}^k - \tilde{\mathbf{u}}^k, \mathbf{v}) + b(p^k - p^{k-1}, \mathbf{v}) = 0$$
   $$b(q, \mathbf{u}^k) = 0$$
   Esta etapa garante que a velocidade final seja livre de divergência e satisfaça as condições de contorno.

Interpretação Algébrica
O artigo interpreta o SIVS como um método de divisão algébrica. Ao realizar uma decomposição LU incompleta (ILU) na matriz do sistema, os autores mostram que o SIVS corresponde a uma fatoração de blocos específica onde a matriz do complemento de Schur $S = -B\tilde{A}^{-1}B^T$ (com $\tilde{A}$ representando o operador viscoso linear) é espectralmente equivalente à matriz de massa da pressão. Crucialmente, ao contrário das iterações de Picard padrão que exigem a inversão de um complemento de Schur diferente a cada passo, o esquema SIVS requer a resolução de um sistema com uma matriz fixa SPD para a variável de pressão em todas as iterações não lineares. Isso permite o pré-cálculo e a reutilização de pré-condicionadores, reduzindo significativamente o custo computacional.

Principais Contribuições e Resultados Teóricos

• Esquema Novel: Introdução do método SIVS para equações de Navier-Stokes estacionárias com estabilização grad-div, garantindo que o campo de velocidade discreto seja livre de divergência e satisfaça condições de contorno de Dirichlet completas.
• Análise de Convergência: Os autores provam o limitamento uniforme das iterações e a convergência geométrica do esquema sob uma condição de pequenez nos dados ($\Lambda < 1/\sqrt{2}$).
• Eficiência Algébrica: O método é mostrado como uma abordagem de divisão algébrica que desacopla as resoluções de velocidade e pressão em sistemas lineares mais simples, comparados ao sistema de ponto de sela acoplado padrão.

Resultados Numéricos
Os autores validam o método usando FreeFem++ com elementos Taylor-Hood ($P_2, P_1$) em três casos de teste:

1. Teste de Convergência: Usando uma solução fabricada, o método alcançou as taxas de convergência ótimas esperadas para erros de velocidade, pressão e divergência. O número de iterações não lineares permaneceu uniforme através de refinamentos de malha.
2. Cavidade Acionada por Tampo 2D: O esquema foi testado em números de Reynolds até $Re=5000$. Os resultados para perfis de velocidade e linhas de corrente corresponderam com precisão aos dados de referência (Ghia et al.). As contagens de iteração foram comparáveis ao método Incremental Picard-Yosida (IPY), mas os autores antecipam menor tempo de CPU devido aos sistemas lineares mais simples.
3. Fluxo 3D ao Redor de um Cilindro: Um problema de referência foi resolvido com $Re \approx 450$ (baseado em $\nu=10^{-3}$). O esquema convergiu em 11 iterações não lineares, produzindo coeficientes de arrasto, sustentação e queda de pressão dentro dos intervalos de referência fornecidos na literatura.
4. Influência de Parâmetros: Testes sobre o parâmetro grad-div $\gamma$ indicaram que valores maiores de $\gamma$ aceleram a convergência para o estado estacionário, embora possam deteriorar o condicionamento do sistema linear na Etapa 1.

Significado
O artigo afirma que o método SIVS oferece uma alternativa computacionalmente eficiente às linearizações de Picard padrão para fluxos incompressíveis estacionários. Ao aproveitar a estabilização grad-div e a divisão algébrica, o método produz sistemas lineares mais simples que são mais fáceis de pré-condicionar e resolver iterativamente. Isso torna a abordagem particularmente adequada para problemas 3D de grande escala onde solutores diretos são inviáveis. Os autores concluem que o método mantém baixas contagens de iteração enquanto reproduz com precisão benchmarks padrão, sugerindo sua viabilidade para aplicações complexas de multipísica.