EDP Sciences, ESAIM: Proceedings, (16), p. 195-210
DOI: 10.1051/proc:2007013
Full text: Download
We aim at improving the accuracy of a Vlasov-Poisson solver using a 2D Particle-In-Cell (PIC) scheme, by denoising its charge density field. To this end, we have used an improvement of Donoho and Johnstone’s wavelet denoising technique. To some extent, our work is a continuation of that performed by J.-P. Chehab et al. [CEMRACS 2003, IRMA Lectures in Mathematics and Theoretical Physics 7, 29–42 (2005), see Zbl 1062.76002 for the entire collection]. Indeed, they made such a study in the one dimensional case and validated their analysis by considering the simulation of the Landau damping phenomenon. They concluded on the efficiency of the method in reducing the number of particles. However, our approach is quite different, since we do not use wavelets to directly interpolate the charge density, but we smooth the density field calculated by the PIC code. This is carried out via an iterative wavelet denoising technique introduced by A. Azzalini, M. Farge and K. Schneider [Appl. Comput. Harmon. Anal. 18, No. 2, 177–185 (2005; Zbl 1061.42022)]. Our work consists in studying the application of the method as a post-processing tool, in view of a future embedding into the PIC code. The results are the following: first, we showed that the hypotheses underlying the application of this method are valid. Secondly, we can infer from this study that it is possible to significantly reduce the amount of data needed for a simulation.