Cosmic Ray Magnetohydrodynamics: A New Two-Moment Framework with Numerical Implementation

Resumo Técnico: Magnetohidrodinâmica de Raios Cósmicos: Uma Nova Estrutura de Dois Momentos com Implementação Numérica

Problema
Os raios cósmicos (RC) são componentes críticos de sistemas astrofísicos, fornecendo feedback em uma ampla gama de escalas, desde a regulação da formação estelar até o acionamento de ventos galácticos. No entanto, modelar o transporte de RC permanece desafiador devido à natureza multiescalar do problema e às complexas microfísicas envolvidas. Embora abordagens recentes de hidrodinâmica de RC de dois momentos tenham melhorado os métodos anteriores de um momento (que sofrem com instabilidades numéricas em extremos de pressão), as estruturas existentes possuem limitações. Especificamente, formulações atuais frequentemente assumem que as ondas de Alfvén se propagam em apenas uma direção ao longo do campo magnético, falhando em capturar configurações de ondas mais gerais, como aquelas impulsionadas pela Instabilidade de Anisotropia de Pressão de Raios Cósmicos (CRPAI). Além disso, desafios numéricos persistem em relação a instabilidades decorrentes de termos de fonte de streaming e ambiguidades no design de solvers para sistemas acoplados de CR-MHD.

Metodologia
Os autores desenvolvem uma nova descrição de dois momentos para a hidrodinâmica de RC derivada de princípios fundamentais, partindo da equação de Boltzmann colisional para partículas relativísticas. Etapas metodológicas principais incluem:

1. Derivação da Teoria Cinética: A estrutura começa com a função de distribuição de RC em um referencial com velocidade de gás perpendicular ($\vec{u}_\perp$). Ao tirar os momentos da equação cinética, os autores derivam equações de evolução para a densidade de energia de RC ($E_{cr}$) e o fluxo de energia ($\vec{F}_{cr}$).
2. Incorporação de Anisotropia e Ondas Bidirecionais: Diferente de modelos anteriores que assumem implicitamente a propagação unidirecional de ondas, esta formulação contabiliza explicitamente a anisotropia de pressão de RC ($\Delta P_{cr}$) e permite ondas de Alfvén propagando-se tanto para frente quanto para trás ao longo do campo magnético. Isso é alcançado expandindo a função de distribuição usando polinômios de Legendre até a segunda ordem.
3. Prescrições de Física de Subgrade: A estrutura incorpora coeficientes de espalhamento ($\sigma$) derivados de processos microfísicos. Os autores fornecem prescrições para:
   • Instabilidade de Streaming de RC (CRSI): Onde gradientes de pressão impulsionam ondas em uma direção específica.
   • Instabilidade de Anisotropia de Pressão de RC (CRPAI): Onde a invariância adiabática em campos magnéticos variáveis impulsiona ondas em ambas as direções.
   • Turbulência Extrínseca: Onde as ondas não possuem uma direção preferencial.
4. Implementação Numérica no Athena++: As equações são implementadas como uma espécie de fluido acoplada ao gás MHD no código Athena++.
   • Solver de Riemann: Um novo solver de Riemann HLLE é empregado. Para mitigar a difusão numérica excessiva, os autores derivam a relação de dispersão completa do sistema CR-MHD para identificar velocidades de sinal realistas (incluindo ondas acústicas de CR) em vez de depender apenas do limite reduzido da velocidade da luz.
   • Integração Temporal: Integradores implícitos (Backward Euler e Implicit van Leer 2) são usados para lidar com termos de fonte rígidos (stiff) provenientes do espalhamento.
   • Critério de Estabilidade: Uma análise de estabilidade de von Neumann é realizada para derivar uma diretriz de passo de tempo especificamente desenhada para suprimir oscilações de escala de grade causadas por termos de streaming com inversão de sinal.

Resultos Principais
Os autores validam o novo módulo através de uma série de testes de benchmark:

• Testes de Streaming 1D: O código reproduz com sucesso soluções analíticas para o streaming de RC sob regimes de amortecimento linear (íon-neutro) e não linear (Landau). Os testes confirmam que o critério de estabilidade derivado mitiga efetivamente instabilidades numéricas quando os coeficientes de espalhamento são grandes.
• Testes de Difusão 1D: O módulo modela com precisão o transporte difusivo impulsionado por turbulência extrínseca, correspondendo a soluções de função de Green para perfis Gaussianos em fundos de gás estáticos e em movimento.
• Transporte Anisotrópico: Em testes 2D com campos magnéticos circulares, o código lida corretamente com o alinhamento do fluxo de RC com as linhas de campo, demonstrando que o termo de força de Lorentz consegue contrapor a inércia perpendicular de RC.
• Coordenadas Curvilíneas: A estrutura funciona corretamente em coordenadas esféricas e cilíndricas, lidando com descontinuidades agudas e efeitos geométricos (ex: campos monopolares e anisotropia de pressão cilíndrica) consistentes com as expectativas teóricas.
• Relações de Dispersão: A derivação da relação de dispersão completa identifica a instabilidade acústica de CR e outras ramificações de ondas, fornecendo benchmarks rigorosos para o código e permitindo o uso de velocidades de sinal fisicamente realistas no solver de Riemann.

Significância e Alegações
O artigo afirma fornecer uma estrutura geral e extensível para a modelagem de feedback de RC que supera as limitações direcionais de esquemas anteriores de dois momentos. Ao acomodar configurações arbitrárias de ondas de Alfvén e incluir explicitamente a anisotropia de pressão de RC, a estrutura pode incorporar diversos mecanismos de transporte, incluindo CRSI, CRPAI e suas combinações. Os autores enfatizam que a inclusão de velocidades de sinal realistas no solver de Riemann e a derivação de um critério de passo de tempo específico aumentam significativamente a robustez numérica. Este trabalho visa facilitar estudos mais precisos de fenômenos impulsionados por RC, como ventos galácticos, interações de nuvens e a estabilidade térmica do meio circumgaláctico, ao fornecer uma ferramenta numericamente estável que faz a ponte entre processos de espalhamento microfísico e dinâmica de fluidos macroscópica.