A family of second order, linear, unconditionally stable methods for the Cahn-Hilliard-Navier-Stokes equations

Resumo Técnico: Uma Família de Métodos de Segunda Ordem, Lineares e Incondicionalmente Estáveis para as Equações de Cahn-Hilliard-Navier-Stokes

Definição do Problema
A simulação numérica precisa de dinâmicas de fluidos multifásicos complexos, especificamente aqueles modelados pelas equações de Cahn-Hilliard-Navier-Stokes (CHNS), apresenta desafios significativos devido à rigidez matemática do operador de quarta ordem de Cahn-Hilliard e sua acoplagem não linear com as equações de Navier-Stokes incompressíveis. Esquemas padrão de passo de tempo explícito enfrentam restrições severas de passo de tempo (por exemplo, $\Delta t \sim O(\Delta x^4)$) para manter a estabilidade. Embora métodos totalmente implícitos ou de divisão convexa (convex-splitting) ofereçam estabilidade de energia incondicional, eles frequentemente exigem a resolução de grandes sistemas não lineares acoplados em cada passo de tempo, aumentando o custo computacional. Há uma necessidade de algoritmos que combinem precisão temporal de segunda ordem, estabilidade incondicional e solvabilidade linear sem sacrificar as propriedades de conservação física do sistema.

Metodologia
Os autores propõem uma família de métodos implícito-explícitos (IMEX) de segunda ordem, lineares e incondicionalmente estáveis para fluxos de duas fases com densidade correspondente (matched-density). O esquema semidiscreto integra três componentes primários:

1. Extrapolação de Termos Não Lineares: Termos convectivos e de acoplamento não lineares são tratados explicitamente usando quantidades extrapoladas (denotadas por $H_{n+\theta}$), evitando a necessidade de iterações não lineares.
2. Formulação de Variável Auxiliar: O termo de energia livre não linear é tratado via uma variável auxiliar $q$ (derivada do potencial de poço duplo $F(\phi)$), reformulando o problema para manter a linearidade.
3. Regularização de Curvatura Temporal: Um mecanismo de estabilização controlado por um parâmetro $\epsilon \geq 0$ é introduzido. Isso envolve um operador de interpolação específico $J^\epsilon_{n+\theta}$ que incorpora a curvatura temporal discreta da solução, um conceito adaptado de trabalhos anteriores sobre regularização de Navier-Stokes.

O algoritmo resultante (Algoritmo 2.1) requer apenas a resolução de sistemas lineares com coeficientes constantes em cada passo de tempo. O esquema é parametrizado por $\theta \in (1/2, 1]$ e $\epsilon \geq 0$.

Principais Contribuições e Resultados Teóricos

• Estabilidade Incondicional: O artigo estabelece uma estimativa de energia discreta rigorosa provando que o esquema proposto é incondicionalmente estável a longo prazo para $\theta \in (1/2, 1]$ e $\epsilon \geq 0$. A prova utiliza matrizes simétricas definidas positivas para definir normas discretas e demonstra que o esquema satisfaz uma lei de dissipação de energia discreta sem restrições de passo de tempo.
• Solvabilidade Linear: Diferente de abordagens totalmente implícitas ou de divisão convexa que podem exigir solvers não lineares, este método produz sistemas lineares em cada passo de tempo, reduzindo significativamente a complexidade computacional.
• Precisão de Segunda Ordem: A análise teórica e os experimentos numéricos confirmam que o método mantém a precisão temporal de segunda ordem.

Resultados Numéricos
Os autores validam o método através de uma série de computações de referência:

• Análise de Convergência: Usando o Método de Soluções Fabricadas (MMS), o esquema demonstra convergência temporal aproximadamente de segunda ordem para velocidade, pressão e variáveis de campo de fase através de diferentes configurações de elementos finitos ($P2-P1-P2-P2$ e $P3-P2-P3-P3$).
• Robustez e Regularização: No problema de cavidade impulsionada por tampa (lid-driven cavity) de duas fases, o estudo examina a interação entre $\theta$ e $\epsilon$. Observa-se que para valores de $\theta$ próximos ao limite de Crank-Nicolson (por exemplo, $\theta=0.51$) com $\epsilon=0$, o esquema pode tornar-se instável devido a oscilações temporais fracamente amortecidas. A introdução de regularização de curvatura positiva ($\epsilon > 0$) amortece efetivamente essas oscilações, restaurando a limitação e a robustez mesmo para passos de tempo maiores.
• Benchmarks Físicos:
  • Decomposição Espinodal: O método captura com precisão a dissipação de energia e a conservação de massa em simulações de longo prazo ($T=100$), mostrando separação de fase e coalescência eficazes.
  • Relaxação de Gotícula: O algoritmo resolve com sucesso o movimento da interface impulsionado pela tensão superficial, transicionando uma gota quadrada para um equilíbrio circular enquanto preserva a massa.
  • Cavidade Impulsionada por Tampa e Instabilidade de Rayleigh-Taylor: O esquema lida com cisalhamento forte, deformação interfacial e instabilidades impulsionadas pela densidade (incluindo o enrolamento de Kelvin-Helmholtz em regimes de baixa viscosidade) sem oscilações espúrias visíveis.

Significância
O artigo afirma que a família de métodos proposta oferece uma alternativa prática e eficiente para sistemas CHNS. Ao desacoplar as não linearidades através de extrapolação e variáveis auxiliares, enquanto estabiliza a discretização temporal via regularização de curvatura, o método alcança um equilíbrio entre precisão de segunda ordem, estabilidade incondicional e custo computacional linear. Os autores enfatizam que o esquema preserva as principais estruturas físicas do modelo, incluindo a dissipação de energia e a conservação de massa, tornando-o adequado para simular fluxos interfaciais complexos em diversos regimes.