Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

Resumo Técnico: Análise de Erro de um Esquema de Elementos Finitos Misto Preservador de Divergência para as Equações Hall–MHD Incompressíveis

Enunciado do Problema
O sistema de magnetohidrodinâmica de Hall (Hall–MHD) incompressível modela a interação entre o movimento de fluidos e campos magnéticos em fluidos condutores eletricamente, incorporando o efeito Hall, que se torna significativo em escalas comparáveis à profundidade de pele iônica. Este efeito introduz um termo rígido e altamente não linear ($\eta \text{curl}(\text{curl} B \times B)$) na equação de indução, levando a comportamento dispersivo, gradientes íngremes e oscilações de alta frequência (ondas de assobio). Essas características apresentam desafios analíticos e computacionais substanciais. Especificamente, a existência global de soluções fortes em três dimensões permanece um problema em aberto, e os métodos numéricos padrão frequentemente lutam com a estabilidade e a exigência estrita de que o campo magnético permaneça solenoidal ($\text{div } B = 0$). Violações desta condição livre de divergência podem induzir monopolos magnéticos não físicos e desencadear instabilidades numéricas severas. Além disso, abordagens numéricas existentes, como aquelas baseadas em métodos de Variável Auxiliar Escalar (SAV) ou técnicas espectrais, frequentemente dissipam funcionais de energia modificados em vez da energia física real, falham em impor a condição livre de divergência exatamente no nível discreto, ou carecem de análise de erro rigorosa e flexibilidade geométrica.

Metodologia
Para abordar esses desafios, os autores propõem um método de elementos finitos misto (FEM) totalmente discreto, linear e que preserva estrutura, para um sistema Hall–MHD regularizado por Voigt. A regularização de Voigt substitui a derivada temporal na equação do momento por uma versão filtrada espacialmente $(I - \alpha_1 \Delta)\partial_t u$ e retém um termo de inércia eletrônica $\alpha_2 \partial_t J$ na lei de Ohm generalizada. Esta modificação melhora a bem-postura e a estabilidade numérica, preservando ao mesmo tempo a estrutura física essencial e as leis de energia do sistema original.

O esquema numérico é construído utilizando os seguintes componentes:

1. Formulação Mista: O problema é reformulado utilizando o campo elétrico $E$ e a densidade de corrente $J$ como variáveis primárias, levando a um sistema acoplado de equações de momento, indução, lei de Ohm e lei de Ampère.
2. Discretização Espacial: O esquema emprega o Cálculo de Elementos Finitos Exterior (FEEC) para garantir que a condição livre de divergência seja satisfeita exatamente.
   • A velocidade $u$ e a pressão $p$ são aproximadas usando o par de elementos finitos MINI (inf-sup estável de ordem mais baixa).
   • O campo magnético $B$ e a densidade de corrente $J$ são aproximados usando os espaços de Raviart–Thomas ($RT_0$) e Nédélec ($X_0$) de ordem mais baixa, respectivamente.
   • Esses espaços formam um complexo de de Rham discreto, garantindo que $\text{div } B_h = 0$ seja preservado exatamente no nível discreto.
3. Discretização Temporal: Um esquema de passo de tempo linearmente implícito e antissimétrico é utilizado. Os termos convectivos e de Lorentz são tratados semi-implicitamente usando valores do passo de tempo anterior para manter a linearidade e a estabilidade energética incondicional.
4. Análise de Erro: Os autores estabelecem taxas de convergência ótimas decompondo o erro em erros de aproximação (usando projeções de Stokes e projeções ortogonais) e erros discretos. Eles utilizam desigualdades de Gaffney discretas e estimativas de energia para lidar com os termos de acoplamento não linear.

Principais Contribuições

• Esquema Preservador de Estrutura: O artigo apresenta um FEM totalmente discreto e linear que impõe a restrição solenoidal no campo magnético exatamente, evitando a necessidade de etapas de projeção ou métodos de penalidade que podem degradar a precisão ou a estabilidade.
• Estabilidade Energética Incondicional: O esquema proposto satisfaz uma lei de dissipação de energia discreta incondicionalmente, espelhando a lei de energia física contínua, mesmo na presença do termo Hall rígido.
• Estimativas de Erro Rigorosas: Os autores fornecem uma análise de erro completa para o problema regularizado por Voigt, estabelecendo taxas de convergência ótimas da ordem $O(h^s + \tau)$ sob suposições de regularidade adequadas.
• Extensão ao Sistema Não Regularizado: Crucialmente, a análise estende-se ao sistema Hall–MHD não regularizado (onde $\alpha_1 = \alpha_2 = 0$). Os autores demonstram que a convergência ainda pode ser provada para o sistema não regularizado, desde que a resistividade $\sigma > 0$. No regime não resistivo ($\sigma = 0$), a regularização de Voigt mostra-se crucial nas estimativas de erro.
• Implementação 2.5D e 3D: A estrutura é validada tanto em configurações 2.5D (invariância translacional) quanto totalmente 3D, demonstrando a aplicabilidade do método a geometrias complexas e fenômenos físicos.

Resultados
Simulações numéricas em 2.5D e 3D confirmam as descobertas teóricas:

• Convergência: O esquema atinge as taxas de convergência espacial e temporal previstas para a velocidade, campo magnético e densidade de corrente. Embora as taxas de convergência para velocidade e densidade de corrente sejam ligeiramente reduzidas em alguns casos (atribuídas à regularidade limitada da solução devido ao termo Hall), o campo magnético mantém convergência ótima.
• Preservação de Divergência: O erro máximo de divergência ($\max | \text{div } B_h |$) permanece no nível da tolerância do solucionador durante toda as simulações, confirmando a preservação exata da restrição solenoidal.
• Estabilidade Energética: A energia total do sistema decai monotonicamente ao longo do tempo, consistente com a lei de dissipação de energia teórica.
• Fidelidade Física: Simulações do fluxo de Arnold–Beltrami–Childress (ABC), do vórtice de Orszag–Tang e da reconexão impulsionada por Hall em uma folha de Harris capturam com sucesso fenômenos físicos complexos, incluindo reconexão magnética, formação de ilhas magnéticas, cordas de fluxo e os padrões de corrente quadrupolares característicos associados à dinâmica Hall–MHD. O esquema resolve essas características com precisão mesmo em malhas relativamente grosseiras.

Significado
O artigo afirma que este trabalho fornece uma estrutura matematicamente rigorosa e computacionalmente robusta para a simulação de sistemas Hall–MHD. Ao combinar a regularização de Voigt com uma formulação de elementos finitos mistos preservadora de divergência, os autores superam os desafios duplos da não linearidade do termo Hall e da exigência estrita de divergência nula. O método oferece uma aproximação fisicamente consistente que é estável, convergente e capaz de lidar com geometrias complexas, preenchendo uma lacuna deixada por métodos anteriores que careciam de análise de erro rigorosa ou preservação exata de divergência. Os resultados sugerem que esta abordagem é uma ferramenta viável para o estudo de fenômenos de plasma astrofísicos e laboratoriais onde os efeitos Hall são dominantes.