Entity -- Hardware-agnostic Particle-in-Cell Code for Plasma Astrophysics. III: Higher-order shape functions & generalized field stencils

Resumo Técnico: Entidade – Código de Partículas em Células Agnóstico a Hardware para Astrofísica de Plasma. III: Funções de Forma de Alta Ordem e Estênceis de Campo Generalizados

Declaração do Problema
Simulações de Partículas em Células (PIC) são essenciais para modelar processos de plasma astrofísico, variando de turbulência e reconexão magnética à aceleração de raios cósmicos. No entanto, essas simulações enfrentam desafios significativos relacionados a ruído numérico, aquecimento artificial e instabilidade de Cherenkov numérica. Esses artefatos frequentemente surgem de esquemas de interpolação de baixa ordem e estênceis de diferenças finitas padrão, que podem obscurecer efeitos físicos ou levar a acelerações espúrias de partículas. Embora o aumento da resolução numérica (células da grade ou partículas por célula) seja uma estratégia comum de mitigação, ela é limitada por recursos computacionais. Os autores identificam a necessidade de métodos de alta ordem para suprimir ruído e melhorar a estabilidade sem incorrer necessariamente em custos computacionais proibitivos, particularmente no contexto do código PIC Entity, projetado para clusters de GPU na escala exa.

Metodologia
Este trabalho (Artigo III da série Entity) implementa e testa dois avanços primários dentro do framework Entity:

1. Deposição de Corrente de Alta Ordem e Interpolação de Campo:
   Os autores implementam o esquema de deposição de corrente conservador de carga proposto por Esirkepov (2001), estendendo-o para suportar funções de forma de ordem arbitrária até a 11ª ordem ($O=11$).

   • Funções de Forma: O código utiliza funções de forma derivadas de B-splines construídas via convolução sucessiva. Essas funções satisfazem condições de simetria, normalização, positividade, extensão finita e suavidade.
   • Consistência: Para garantir consistência com a deposição de corrente, a interpolação de campo da grade para as posições das partículas também é atualizada para usar as mesmas funções de forma de alta ordem, substituindo esquemas de interpolação indireta anteriores.
   • Nota de Implementação: Ordens $O=10$ e $O=11$ requerem aritmética de precisão dupla devido a valores numéricos pequenos nas fronteiras do estêncil, limitando seu uso em certas arquiteturas de GPU.

2. Estênceis de Campo Generalizados para o Solucionador de Maxwell:
   O solucionador de Faraday (evoluindo o campo magnético via $\partial_t B = -c\nabla \times E$) é atualizado para usar estênceis de diferenças finitas generalizados.

   • Otimização: Em vez do estêncil Yee padrão de segunda ordem, os autores implementam estênceis que utilizam pontos de grade adjacentes adicionais. Os coeficientes desses estênceis são otimizados para minimizar o desvio entre a velocidade de fase numérica e a velocidade física da luz no espaço de Fourier.
   • Controle de Dispersão: Essa abordagem visa suprimir a instabilidade de Cherenkov numérica, que ocorre quando partículas relativísticas viajam mais rápido que a velocidade numérica da luz na grade.
   • Preservação de Divergência: Os autores garantem que configurações específicas de estêncil satisfaçam a condição $\beta_{ij} = \delta_j$ para impor automaticamente $\nabla \cdot B = 0$ ao longo de toda a simulação.

Principais Contribuições e Resultados
Os autores realizam extensos testes de convergência para validar os novos esquemas contra implementações padrão (deposição Zigzag e estêncil Yee) e várias configurações de filtro.

• Conservação de Carga: O esquema Esirkepov atualizado demonstra excelente conservação de carga em todas as ordens testadas ($O=1$ a $O=11$), com precisão limitada principalmente pelo erro de truncamento e não pela ordem da função de forma.
• Poder de Parada: Em testes medindo a força de parada eletrostática em partículas relativísticas, funções de forma de alta ordem reduzem significativamente a perda de energia numérica. A transição da 1ª para a 9ª ordem reduz o poder de parada efetivo, equivalente a aumentar a resolução numérica por um fator de 3–4.
• Aquecimento Numérico e Estabilidade:
  • Plasma Frio: Para plasmas frios com comprimentos de Debye sub-resolvidos, esquemas de alta ordem melhoram drasticamente a estabilidade. A transição da 1ª para a 2ª ordem reduz a resolução de grade necessária por um fator de $\sim4$; a transição para a 3ª ordem reduz por $\sim8$.
  • Modos Oblíquos: Em testes envolvendo feixes de raios cósmicos, formas de alta ordem suprimem o crescimento de ruído de aliasing numérico. Embora a 6ª ordem pareça suficiente para convergência nesses testes específicos, ordens mais altas atrasam ainda mais o início da instabilidade.
  • Turbulência Decrescente: Em turbulência decrescente 2D, funções de forma de alta ordem melhoram a conservação de energia, reduzindo o erro de energia total em até duas ordens de grandeza em comparação com a deposição Zigzag padrão.
• Desempenho do Estêncil de Campo:
  • Partícula Única: Testes com um único elétron relativístico mostram que estênceis otimizados (por exemplo, Cowan, min3, min6) reduzem significativamente artefatos de Cherenkov numérico em comparação com o estêncil Yee padrão.
  • Plasma Derivante: Em simulações de plasma quente e derivante, os estênceis min3 e min6 suprimem efetivamente o crescimento da instabilidade de Cherenkov e o aquecimento numérico associado.
  • Choques Relativísticos: Em uma simulação de choque de pares relativísticos, a combinação do estêncil min3 e funções de forma de 3ª ordem estabiliza o plasma a montante, prevenindo o aquecimento artificial e a filamentação de corrente observados com o estêncil Yee. Isso leva a estruturas de rampa de choque mais precisas e espectros de aceleração de partículas.
• Compensações de Desempenho:
  • 1D: O impacto no desempenho de esquemas de alta ordem é negligenciável.
  • 2D: O custo aumenta moderadamente (por exemplo, $\sim2,5\times$ para 3ª ordem, $\sim24\times$ para 11ª ordem).
  • 3D: O custo torna-se proibitivo para ordens $O > 3$ (por exemplo, $\sim5\times$ para 3ª ordem).
  • Mitigação: Os autores demonstram que o custo computacional aumentado da interpolação de alta ordem pode ser compensado reduzindo o número de passadas de filtro de corrente ou o número de partículas por célula, uma vez que esquemas de alta ordem fornecem inerentemente maior estabilidade.

Significado e Alegações
O artigo alega que, embora nenhum método único resolva todos os desafios numéricos em simulações PIC, a integração de funções de forma de alta ordem e estênceis de campo generalizados no Entity oferece um conjunto robusto de ferramentas para melhorar a fidelidade da simulação. Os autores enfatizam que:

1. Estabilidade vs. Resolução: Métodos de alta ordem permitem simulações estáveis de plasmas frios e choques relativísticos em resoluções mais baixas do que anteriormente possível, potencialmente reduzindo custos computacionais totais apesar da sobrecarga por passo.
2. Supressão de Ruído: Esses métodos suprimem efetivamente efeitos de Cherenkov numérico e ruído de aliasing, o que é crítico para estudar instabilidades físicas como a instabilidade de Bell ou o fluxo de raios cósmicos.
3. Flexibilidade: Os estênceis generalizados permitem que os usuários ajustem o solucionador para relações de dispersão específicas ou simetrias direcionais (por exemplo, usando o estêncil min3 para estabilidade geral ou Cowan/min6 para grandes números CFL).
4. Recomendações Práticas: Os autores recomendam o uso de funções de forma de 2ª ou 3ª ordem como um equilíbrio prático entre precisão e desempenho, e a utilização dos estênceis min3 (2D) ou min3D2 (3D) para simulações de propósito geral.

O trabalho conclui que essas melhorias estão agora disponíveis no Entity v1.3.0, permitindo modelagem mais precisa e estável de fenômenos complexos de plasma astrofísico em hardware na escala exa.