Improving divergence cleaning in cosmological SPMHD simulations

Resumo Técnico: Melhorando a Limpeza de Divergência em Simulações de SPMHD Cosmológicas

Definição do Problema
Em simulações de magnetohidrodinâmica por partículas suavizadas (SPMHD) cosmológicas, manter a condição solenoidal ($\nabla \cdot \mathbf{B} = 0$) é um desafio numérico persistente. Os esquemas padrão em códigos como o Gadget dependem tipicamente da abordagem de "8 ondas" de Powell et al. (1999). Embora este método advecte os erros de divergência com o fluxo do fluido, ele não os remove ativamente. Em códigos baseados em SPH, os erros residuais de $\nabla \cdot \mathbf{B}$ podem atingir ordem unitária, levando à amplificação não física do campo magnético e à "instabilidade de tensão", onde termos de fonte espúrios na equação do momento dominam quando a pressão magnética excede a pressão do gás. Embora o esquema de limpeza hiperbólico/parabólico proposto por Dedner et al. (2002) ofereça uma solução, sua formulação original provou ser instável em SPMHD, particularmente em ambientes dinâmicos com passos de tempo individuais e saltos de densidade fortes, agindo frequentemente como uma fonte não física de energia magnética.

Metodologia
Os autores implementam uma versão restrita e conservadora de energia do algoritmo de limpeza do tipo Dedner (Tricco & Price 2012; Tricco et al. 2016a) no código de SPMHD cosmológica OpenGadget3. As principais adaptações metodológicas incluem:

• Conservação de Energia: As equações equacionadas são reformuladas para garantir a conservação da energia total. Isso envolve evoluir a variável de campo de limpeza como $\psi/c_h$ (onde $\psi$ é o escalar auxiliar e $c_h$ é a velocidade de limpeza) em vez de apenas $\psi$, e incluir um termo de fonte específico $-\frac{1}{2}\psi(\nabla \cdot \mathbf{v})$ na equação de evolução. Isso garante que o amortecimento parabólico do campo de limpeza remova a energia magnética sem introduzir crescimento ilimitado.
• Adaptação Cosmológica: As equações são adaptadas para coordenadas comoventes e passos de tempo individuais de partículas. A velocidade de limpeza $c_h$ é definida como a velocidade magnetossônica rápida local ($\sqrt{v_A^2 + c_s^2}$), garantindo que a limpeza seja intrinsecamente estável em termos de Courant dentro das restrições de passo de tempo padrão do SPH.
• Configuração da Simulação: Os autores realizam uma série de simulações em um aglomerado de galáxias massivo ($M_{200c} \approx 10^{15} M_\odot$) semeado com um campo magnético uniforme ($10^{-14}$ G comóvel). Cinco execuções distintas são comparadas:
  1. NIFTY-noclean: Execução de referência usando apenas o esquema de advecção de Powell.
  2. NIFTY-clean: Execução fiduciária com o esquema de limpeza restrito completo (incluindo retroação de entropia).
  3. NIFTY-clean2: O mesmo que o anterior, mas com um fator de sobrelimpeza $\zeta=2$ (limpeza a duas vezes a velocidade de onda local).
  4. NIFTY-advect: Esquema de limpeza com o termo de fonte de divergência de velocidade desligado.
  5. NIFTY-heat: Esquema de limpeza onde a energia de limpeza dissipada é descartada em vez de ser alimentada de volta na equação de entropia.

Resultos Principais

• Redução do Erro de Divergência: O resultado primário é uma redução nos erros de divergência relativa ($|\nabla \cdot \mathbf{B}|/B$) em 2 a 3 ordens de magnitude em todo o volume do aglomerado em comparação com o esquema apenas de Powell. Esta redução é uniforme tanto no núcleo quanto nas periferias do aglomerado.
• Amplificação do Campo Magnético: Nas periferias do aglomerado ($r \approx 1\text{--}3 \, h^{-1}$ Mpc), onde o campo magnético é mais fraco e a resolução é menor, o esquema de limpeza ativo permite uma amplificação do campo magnético de um fator de 5 a 10 em comparação com a execução apenas de Powell. Os autores atribuem isso à remoção de forças de Lorentz não físicas causadas por erros de divergência, que anteriormente suprimiam o dínamo de pequena escala nessas regiões.
• Estabilidade do Núcleo: A estrutura de densidade global e a geometria do campo magnético no núcleo do aglomerado permanecem amplamente inalteradas entre as execuções com e sem limpeza, indicando que a limpeza não altera artificialmente o estado saturado do dínamo do núcleo.
• Consistência Termodinâmica: A inclusão da retroação de entropia (alimentar a energia de limpeza dissipada de volta como calor) mostra-se crítica para a consistência termodinâmica. Execuções sem esta retroação (NIFTY-heat) resultam em menor entropia central e campos magnéticos centrais saturados mais altos devido à compressão adiabática de um núcleo mais frio.
• Custo Computacional: A implementação não introduz nenhum custo de tempo de execução mensurável; a limpeza melhorada é alcançada essencialmente "de graça" dentro das restrições de passo de tempo de Courant existentes.

Significância e Alegações
O artigo afirma que o controle preciso dos erros de divergência é essencial para capturar a amplificação do campo magnético nas regiões de baixa densidade e mal resolvidas das periferias de aglomerados de galáxias. Os autores argumentam que os erros de divergência não controlados nos esquemas padrão de Powell atuam como ruído numérico que suprime o dínamo em regiões sub-saturadas. Ao implementar o esquema de limpeza restrito e conservador de energia, a simulação alcança uma realização mais fiel da física de MHD sem penalidade computacional.

Os autores concluem que o esquema de limpeza restrito:

1. Reduz dramaticamente os erros de divergência numérica (2 a 3 ordens de magnitude).
2. Permite uma representação mais precisa da amplificação do campo em regiões sub-resolvidas.
3. Produz estruturas de campo magnético mais suaves e coerentes.
4. Estabiliza a evolução do campo magnético através de diferentes resoluções, melhorando a fidelidade física dos resultados mesmo em resoluções moderadas.

O trabalho recomenda que aplicações de SPMHD cosmológica adotem este esquema de limpeza por padrão, utilizando a velocidade de limpeza natural ($\zeta=1$) e mantendo a retroação de energia parabólica ativa para fechar o balanço de energia. A implementação está agora disponível publicamente no OpenGadget3.