Numerical thermalization in $n$-D particle-in-cell simulations

Resumen Técnico: Termalización Numérica en Simulaciones de Partícula en Celda de n-D

Planteamiento del Problema
Las simulaciones de partícula en celda (PIC) se utilizan ampliamente para modelar el comportamiento cinético de plasmas mediante la representación de partículas físicas como macropartículas con funciones de forma finitas. Aunque a menudo se asume que resuelven la ecuación de Vlasov sin colisiones debido a los efectos de suavizado de la interpolación de la malla, las simulaciones PIC sufren frecuentemente de una colisionalidad artificialmente alta. Esta "termalización numérica" surge cuando el peso de la macropartícula (la relación entre partículas físicas y macropartículas) es grande, amplificando las fuerzas interparticulares a pesar del suavizado de la malla. En dimensiones superiores (2D y 3D), lograr que un tiempo de termalización numérica ($\tau_{num}$) sea más largo que el tiempo de termalización física ($\tau_{phys}$) es a menudo intratable, incluso cuando la malla resuelve la longitud de Debye. Esta colisionalidad artificial distorsiona las distribuciones de velocidad y afecta magnitudes clave de interés (por ejemplo, transporte, emisión, ruptura) durante escalas de tiempo de simulación $t \gtrsim \tau_{num}$, lo que puede conducir potencialmente a una falsa convergencia o a predicciones físicas incorrectas.

Metodología
Los autores validan una teoría cinética para tiempos de colisión numéricos frente a simulaciones PIC electrostáticas en dimensionalidades 1D, 2D y 3D.

1. Simulaciones: Se realizaron simulaciones PIC electrostáticas de una sola especie con condiciones iniciales uniformes. Se calculó la función de autocorrelación de velocidad (VACF), $Z(t)$, a partir de los datos de trayectoria de las partículas. El tiempo de colisión numérica se definió como el tiempo de decaimiento exponencial de $Z(t)$.
2. Parámetros: Las simulaciones variaron el peso de la macropartícula ($W$), el espaciamiento de la malla ($\Delta$) y el radio de la forma de la macropartícula ($R$). Tanto se realizaron simulaciones sin filtrar (donde $R = \Delta/2$) como simulaciones filtradas (donde $R > \Delta/2$ mediante filtrado en el espacio $k$) para aislar los efectos de la forma de los efectos de calentamiento de la malla.
3. Teoría: El estudio emplea la teoría cinética desarrollada por Okuda, Birdsall y Langdon. Este marco incorpora la función de forma y el peso de la macropartícula en una ecuación cinética de tipo Lenard-Balescu. La teoría deriva operadores de colisión para macropartículas con forma, contabilizando la modificación del núcleo de fuerza interparticular por la función de forma $S(k)$ y el peso $W$.
4. Comparación: Las predicciones teóricas para el tiempo de colisión se derivaron utilizando aproximaciones dinámica, estática y de baja velocidad del integral de arrastre. Estas se compararon directamente contra las tasas de decaimiento observadas en las simulaciones PIC.

Contribuciones Clave y Resultados

• Validación de la Teoría Cinética: El artículo demuestra que la teoría cinética de Okuda, Birdsall y Langdon predice con precisión el tiempo de colisión numérica en simulaciones PIC 1D, 2D y 3D. La teoría captura correctamente la dependencia de $\tau_{num}$ respecto al peso de la macropartícula, el tamaño de la forma y la dimensionalidad.
• Interpretación del Núcleo de Fuerza: El estudio aclara que las colisiones numéricas son impulsadas por una fuerza interparticular efectiva donde la función de forma finita debilita y suaviza la fuerza, mientras que el peso de la macropartícula la amplifica. La teoría muestra que la relación entre los tiempos de colisión numéricos y físicos es independiente de las razones de masa y carga de las especies, dependiendo principalmente de $W$ y $R$.
• Efectos de Dimensionalidad: Los resultados resaltan que alcanzar el tiempo de termalización física es significativamente más difícil en 3D que en 1D o 2D para tamaños de macropartícula que resuelven la longitud de Debye. En 3D, el tiempo de colisión escala logarítmicamente con el parámetro de plasma $\Lambda$ en el límite de un radio de partícula pequeño ($R \ll \lambda_D$).
• Análisis de Compromiso (Trade-off): El artículo cuantifica el compromiso entre reducir la termalización numérica y mantener la exactitud física. Aumentar el tamaño de la macropartícula ($R$) para reducir las colisiones (mediante el suavizado de fuerzas) debilita simultáneamente el campo medio de Vlasov y altera los procesos sin colisiones (por ejemplo, la dispersión y amortiguación de ondas). Existe un límite de "campo distorsionado" donde aumentar $R$ viola las restricciones físicas sobre la resolución espacial.
• Formas de Orden Superior: El análisis se extiende a funciones de forma de orden superior (por ejemplo, B-splines). Se muestra que aumentar el orden de la spline es matemáticamente equivalente a aumentar el ancho efectivo de la partícula. Por lo tanto, el uso de formas de orden superior no altera fundamentalmente la escala de la termalización numérica, sino que aumenta efectivamente el tamaño de la partícula.
• Aplicación a Sistemas Reales: Los autores aplican las restricciones derivadas a varios estudios PIC existentes (por ejemplo, puntos calientes de ICF, reconexión magnética, descargas de Penning, propulsores Hall, descargas de streamers e inyectores de haces neutros de ITER). En muchos casos 3D, se encuentra que el tiempo de colisión numérica es más corto que el tiempo de colisión física, lo que sugiere que la termalización numérica puede dominar o influir significativamente en los resultados, a menos que se empleen estrategias de mitigación específicas (como colisiones de Monte Carlo o un refinamiento extremo de la malla).

Significancia
La principal significancia de este trabajo es la provisión de un marco predictivo y validado para estimar la escala de tiempo de las interacciones de Coulomb autoconsistentes en simulaciones PIC. Al establecer que la teoría cinética predice con precisión $\tau_{num}$ como una función del espaciamiento de la malla y el peso de la macropartícula, el artículo proporciona a los usuarios de PIC los medios para:

1. Caracterizar la Relevancia: Determinar a priori si la termalización numérica impactará la física de interés específica (por ejemplo, si $\tau_{num}$ es más corto que el tiempo de simulación o el tiempo de colisión física).
2. Establecer Restricciones: Definir regímenes operativos permisibles para el peso de la macropartícula y el tamaño de la forma para asegurar que los artefactos numéricos no oscurezcan los fenómenos físicos.
3. Evitar la Falsa Convergencia: Reconocer que la convergencia con respecto a partículas por celda no garantiza la eliminación de la termalización numérica; se requieren conteos de macropartículas suficientes para resolver la tasa de colisión física, no solo las características de ruido.

El artículo concluye que, si bien los métodos PIC son poderosos, la suposición "sin colisiones" es a menudo inválida en simulaciones 3D con parámetros estándar. Los usuarios deben contabilizar explícitamente la termalización numérica para asegurar la fidelidad de los procesos de transporte y colisionales de larga escala temporal.