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

Resumen Técnico: Entidad – Código de Partículas en Celda agnóstico al hardware para Astrofísica de Plasmas. III: Funciones de forma de orden superior y plantillas de campo generalizadas

Enunciado del Problema
Las simulaciones de Partículas en Celda (PIC) son esenciales para modelar procesos de plasma astrofísico, que van desde la turbulencia y la reconexión magnética hasta la aceleración de rayos cósmicos. Sin embargo, estas simulaciones enfrentan desafíos significativos relacionados con el ruido numérico, el calentamiento artificial y la inestabilidad de Cherenkov numérica. Estos artefactos a menudo surgen de esquemas de interpolación de bajo orden y plantillas de diferencias finitas estándar, que pueden oscurecer efectos físicos o conducir a una aceleración de partículas espuria. Aunque aumentar la resolución numérica (celdas de la cuadrícula o partículas por celda) es una estrategia común de mitigación, está limitada por los recursos computacionales. Los autores identifican la necesidad de métodos de orden superior para suprimir el ruido y mejorar la estabilidad sin incurrir necesariamente en costos computacionales prohibitivos, particularmente en el contexto del código PIC Entity, diseñado para clústeres de GPU de exaescala.

Metodología
Este trabajo (Papel III de la serie Entity) implementa y prueba dos avances principales dentro del marco de Entity:

1. Deposición de Corriente de Alto Orden e Interpolación de Campo:
   Los autores implementan el esquema de deposición de corriente que conserva la carga, propuesto por Esirkepov (2001), extendiéndolo para soportar funciones de forma de orden arbitrario hasta el undécimo orden ($O=11$).

   • Funciones de Forma: El código utiliza funciones de forma derivadas de B-splines, construidas mediante convolución sucesiva. Estas funciones satisfacen las condiciones de simetría, normalización, positividad, extensión finita y suavidad.
   • Consistencia: Para garantizar la consistencia con la deposición de corriente, la interpolación de campo desde la cuadrícula hacia las posiciones de las partículas también se actualiza para utilizar las mismas funciones de forma de alto orden, reemplazando los esquemas de interpolación indirecta anteriores.
   • Nota de Implementación: Los órdenes $O=10$ y $O=11$ requieren aritmética de doble precisión debido a los pequeños valores numéricos en los límites de la plantilla, lo que limita su uso en ciertas arquitecturas de GPU.

2. Plantillas de Campo Generalizadas para el Solucionador de Maxwell:
   El solucionador de Faraday (que evoluciona el campo magnético mediante $\partial_t B = -c\nabla \times E$) se actualiza para utilizar plantillas de diferencias finitas generalizadas.

   • Optimización: En lugar de la plantilla Yee estándar de segundo orden, los autores implementan plantillas que utilizan puntos de cuadrícula adyacentes adicionales. Los coeficientes de estas plantillas se optimizan para minimizar la desviación entre la velocidad de fase numérica y la velocidad física de la luz en el espacio de Fourier.
   • Control de Dispersión: Este enfoque busca suprimir la inestabilidad de Cherenkov numérica, que ocurre cuando partículas relativistas viajan más rápido que la velocidad numérica de la luz en la cuadrícula.
   • Preservación de la Divergencia: Los autores aseguran que configuraciones específicas de la plantilla satisfagan la condición $\beta_{ij} = \delta_j$ para imponer automáticamente $\nabla \cdot B = 0$ durante toda la simulación.

Contribuciones y Resultados Clave
Los autores realizan pruebas exhaustivas de convergencia para validar los nuevos esquemas frente a implementaciones estándar (deposición Zigzag y plantilla Yee) y varias configuraciones de filtros.

• Conservación de Carga: El esquema de Esirkepov actualizado demuestra una excelente conservación de la carga en todos los órdenes probados ($O=1$ a $O=11$), con una precisión limitada principalmente por el error de truncamiento más que por el orden de la función de forma.
• Poder de Frenado: En pruebas que miden la fuerza de frenado electrostática sobre partículas relativistas, las funciones de forma de orden superior reducen significativamente la pérdida de energía numérica. Pasar del primer al noveno orden reduce el poder de frenado efectivo, equivalente a aumentar la resolución numérica por un factor de 3–4.
• Calentamiento Numérico y Estabilidad:
  • Plasma Frío: Para plasmas fríos con longitudes de Debye subresueltas, los esquemas de orden superior mejoran drásticamente la estabilidad. Pasar del primer al segundo orden reduce la resolución de cuadrícula requerida por un factor de $\sim4$; pasar al tercer orden la reduce en $\sim8$.
  • Modos Oblicuos: En pruebas que involucran haces de rayos cósmicos, las formas de orden superior suprimen el crecimiento del ruido de aliasing numérico. Aunque el sexto orden parece suficiente para la convergencia en estas pruebas específicas, órdenes superiores retrasan aún más el inicio de la inestabilidad.
  • Turbulencia Decreciente: En turbulencia 2D decreciente, las funciones de forma de orden superior mejoran la conservación de la energía, reduciendo el error total de energía en hasta dos órdenes de magnitud en comparación con la deposición Zigzag estándar.
• Rendimiento de la Plantilla de Campo:
  • Partícula Única: Las pruebas con un solo electrón relativista muestran que las plantillas optimizadas (por ejemplo, Cowan, min3, min6) reducen significativamente los artefactos de Cherenkov numérico en comparación con la plantilla Yee estándar.
  • Plasma Derivante: En simulaciones de plasma cálido y derivante, las plantillas min3 y min6 suprimen eficazmente el crecimiento de la inestabilidad de Cherenkov y el calentamiento numérico asociado.
  • Choques Relativistas: En una simulación de choque de pares relativista, la combinación de la plantilla min3 y las funciones de forma de tercer orden estabiliza el plasma aguas arriba, previniendo el calentamiento artificial y la filamentación de corriente observados con la plantilla Yee. Esto conduce a estructuras de rampa de choque más precisas y espectros de aceleración de partículas.
• Compensaciones de Rendimiento:
  • 1D: El impacto en el rendimiento de los esquemas de orden superior es insignificante.
  • 2D: El costo aumenta moderadamente (por ejemplo, $\sim2.5\times$ para tercer orden, $\sim24\times$ para undécimo orden).
  • 3D: El costo se vuelve prohibitivo para órdenes $O > 3$ (por ejemplo, $\sim5\times$ para tercer orden).
  • Mitigación: Los autores demuestran que el aumento del costo computacional de la interpolación de orden superior puede compensarse reduciendo el número de pasos de filtro de corriente o el número de partículas por celda, ya que los esquemas de orden superior proporcionan inherentemente mayor estabilidad.

Significado y Afirmaciones
El artículo afirma que, aunque ningún método único resuelve todos los desafíos numéricos en las simulaciones PIC, la integración de funciones de forma de orden superior y plantillas de campo generalizadas en Entity ofrece un conjunto de herramientas robusto para mejorar la fidelidad de la simulación. Los autores enfatizan que:

1. Estabilidad vs. Resolución: Los métodos de orden superior permiten simulaciones estables de plasmas fríos y choques relativistas a resoluciones más bajas que las posibles anteriormente, reduciendo potencialmente los costos computacionales totales a pesar de la sobrecarga por paso.
2. Supresión de Ruido: Estos métodos suprimen eficazmente los efectos de Cherenkov numérico y el ruido de aliasing, lo cual es crítico para estudiar inestabilidades físicas como la inestabilidad de Bell o el flujo de rayos cósmicos.
3. Flexibilidad: Las plantillas generalizadas permiten a los usuarios ajustar el solucionador para relaciones de dispersión específicas o simetrías direccionales (por ejemplo, usando la plantilla min3 para estabilidad general o Cowan/min6 para números CFL grandes).
4. Recomendaciones Prácticas: Los autores recomiendan usar funciones de forma de segundo o tercer orden como un equilibrio práctico entre precisión y rendimiento, y utilizar las plantillas min3 (2D) o min3D2 (3D) para simulaciones de propósito general.

El trabajo concluye que estas mejoras están ahora disponibles en Entity v1.3.0, permitiendo un modelado más preciso y estable de fenómenos complejos de plasma astrofísico en hardware de exaescala.