Reducing Noise for PIC Simulations Using Kernel Density Estimation Algorithm

TL;DR

This paper introduces a Kernel Density Estimation-based framework to significantly reduce noise in Particle-In-Cell simulations, improving accuracy and efficiency by optimizing shape functions and particle widths.

Contribution

It develops a novel KDE-based approach with optimal shape functions and width adjustment algorithms to reduce noise and improve PIC simulation accuracy.

Findings

Noise reduced by 98% in density estimation.

Simulation accuracy improved for linear damping rate.

Computational efficiency increased by 40%.

Abstract

Noise is a major concern for Particle-In-Cell (PIC) simulations. We propose a new theoretical and algorithmic framework to evaluate and reduce the noise level for PIC simulations based on the Kernel Density Estimation (KDE) theory, which has been widely adopted in machine learning and big data science. According to this framework, the error on particle density estimation for PIC simulations can be characterized by the Mean Integrated Square Error (MISE), which consists of two parts, systematic error and noise. A careful analysis shows that in the standard PIC methods noise is the dominate error, and the noise level can be reduced if we select different shape functions that are capable of balancing the systematic error and the noise. To improve performance, we use the von Mises distribution as the shape function and seek an optimal particle width that minimizes the MISE, represented by a Cross-Validation (CV) function. This procedure significantly reduces both the noise and the MISE for PIC simulations. A particle-wise width adjustment algorithm and a width update algorithm are further developed to reduce the MISE. Simulations using the examples of Langmuir wave and Landau Damping demonstrate that the KDE algorithm developed in the present study reduces the noise level on density estimation by 98%, and gives a much more accurate result on the linear damping rate compared to the standard PIC methods. Meanwhile, it is computational efficient that can save 40% time to achieve the same accuracy.