Optimal error estimates for a fully discrete, highly efficient decoupled scheme for the 2D/3D diffuse interface two-phase MHD flows

Resumo Técnico: Estimativas de Erro Ótimas para um Esquema Decoplado de Alta Eficiência para Fluxos MHD de Duas Fases com Interface Difusa

1. O Problema
O artigo aborda a modelagem e simulação numérica de sistemas de Magnetohidrodinâmica (MHD) de duas fases, que descrevem o comportamento dinâmico de dois fluidos condutores incompressíveis e imiscíveis sob a influência de um campo eletromagnético externo. O modelo governante acopla três conjuntos de equações:

• Equações de Cahn-Hilliard: Descrevem a evolução da interface difusa entre os fluidos.
• Equações de Navier-Stokes: Descrevem a hidrodinâmica e o transporte de momento.
• Equações de Maxwell: Descrevem a evolução do campo magnético.

O desafio principal reside na forte não-linearidade, no acoplamento complexo entre os campos (forças de Lorentz, tensões de interface) e na necessidade de desenvolver esquemas numéricos que sejam simultaneamente eficientes, estáveis energeticamente e que ofereçam estimativas de erro precisas (ótimas) para todas as variáveis, especialmente em esquemas decoplados.

2. Metodologia
Os autores propõem um esquema de elementos finitos totalmente discreto e decoplado baseado na técnica de divisão convexa (convex-splitting). As principais características metodológicas incluem:

• Discretização Temporal: Utilização do esquema semi-implícito de Euler para garantir estabilidade incondicional.
• Discretização Espacial: Emprego de elementos finitos de Taylor-Hood (para velocidade e pressão) e elementos Mini, garantindo a estabilidade inf-sup.
• Estratégia de Decoplamento: Adoção de um esquema de correção de pressão (pressure-correction scheme). Isso permite resolver as equações de velocidade, pressão e campo magnético de forma desacoplada em cada passo de tempo, reduzindo significativamente o custo computacional comparado a esquemas totalmente acoplados.
• Técnicas de Análise de Erro Inovadoras: Para superar as limitações de precisão causadas pela "poluição artificial" da aproximação do campo de fase ($\phi$) em estimativas de erro tradicionais, os autores utilizam projeções elípticas de Ritz e de Stokes quasi-projeções (quasi-projections), além de uma projeção de Maxwell. Essas ferramentas permitem isolar os erros de aproximação e obter taxas de convergência ótimas.

3. Contribuições Principais

• Estimativas de Erro Ótimas: O trabalho deriva pela primeira vez, para este tipo de esquema decoplado, estimativas de erro ótimas nas normas $L^2$ e $H^1$ para o campo de fase, velocidade, pressão e campo magnético.
  • Para elementos de Taylor-Hood de ordem $r \ge 2$, obtém-se convergência ótima de ordem $O(\Delta t + h^{r+1})$ para o campo de fase e $O(\Delta t + h^{r+2})$ para velocidade e campo magnético (em norma $L^2$).
  • Para o caso $r=1$ (elementos lineares), os resultados numéricos confirmam a optimalidade, embora a análise teórica para velocidade e campo magnético mostre uma ordem ligeiramente inferior em certos casos, o que é consistente com a literatura, mas aqui refinado.
• Estabilidade Energética Incondicional: O esquema proposto é provado ser incondicionalmente estável energeticamente, preservando a lei de dissipação de energia do sistema contínuo no nível discreto, independentemente do tamanho do passo de tempo ($\Delta t$).
• Tratamento da Poluição de Aproximação: A introdução e aplicação rigorosa das quasi-projeções de Ritz e Stokes para lidar com o termo de acoplamento entre o gradiente do campo de fase e a velocidade, eliminando a perda de ordem de convergência observada em análises anteriores.

4. Resultados

• Análise Teórica: O Teorema 2.1 estabelece as estimativas de erro rigorosas. Por exemplo, para elementos de Taylor-Hood com $r \ge 2$, o erro na norma $L^2$ para velocidade e campo magnético é da ordem de $O(\Delta t + h^{r+2})$, e para o campo de fase é $O(\Delta t + h^{r+1})$.
• Validação Numérica:
  • Convergência: Simulações em 2D e 3D utilizando o software FreeFem confirmam as taxas de convergência teóricas para diferentes pares de elementos finitos (Casos I e II).
  • Decomposição Espinodal: Simulações mostram a evolução da separação de fases, demonstrando a conservação de massa e a estabilidade da energia do sistema ao longo do tempo.
  • Instabilidade de Kelvin-Helmholtz: O método foi aplicado para simular a instabilidade de K-H em 2D e 3D, tanto com perturbações de modo único quanto duplo. Os resultados capturam corretamente a formação de vórtices em espiral ("olho de gato") e a interação com o campo magnético, validando a capacidade do esquema de lidar com fenômenos físicos complexos e acoplados.

5. Significância
Este trabalho é significativo porque preenche uma lacuna importante na literatura numérica de MHD de duas fases. Até então, a maioria das análises de erro para esquemas decoplados não alcançava a optimalidade nas normas $L^2$ para todas as variáveis devido aos efeitos de acoplamento não lineares. Ao fornecer estimativas de erro ótimas e provar a estabilidade incondicional, o artigo oferece uma base teórica sólida para o uso de esquemas computacionalmente eficientes (decoplados) em aplicações de alta fidelidade, como fusão nuclear, metalurgia e bombas de metal líquido, onde a precisão e a estabilidade são críticas. Além disso, a metodologia de quasi-projeções pode ser adaptada para outros problemas de acoplamento de múltiplos campos físicos.