AP-Cloud: Adaptive Particle-in-Cloud method for optimal solutions to Vlasov-Poisson equation

We propose a new adaptive Particle-in-Cloud (AP-Cloud) method for obtaining optimal numerical solutions to the Vlasov-Poisson equation. Unlike the traditional particle-in-cell (PIC) method, which is commonly used for solving this problem, the AP-Cloud adaptively selects computational nodes or particles to deliver higher accuracy and efficiency when the particle distribution is highly non-uniform. Unlike other adaptive techniques for PIC, our method balances the errors in PDE discretization and Monte Carlo integration, and discretizes the differential operators using a generalized finite difference (GFD) method based on a weighted least square formulation. As a result, AP-Cloud is independent of the geometric shapes of computational domains and is free of artificial parameters. Efficient and robust implementation is achieved through an octree data structure with 2:1 balance. We analyze the accuracy and convergence order of AP-Cloud theoretically, and verify the method using an electrostatic problem of a particle beam with halo. Simulation results show that the AP-Cloud method is substantially more accurate and faster than the traditional PIC, and it is free of artificial forces that are typical for some adaptive PIC techniques.

[1]  E. Oñate,et al.  A FINITE POINT METHOD IN COMPUTATIONAL MECHANICS. APPLICATIONS TO CONVECTIVE TRANSPORT AND FLUID FLOW , 1996 .

[2]  Georges-Henri Cottet,et al.  Particle Methods for the One-Dimensional Vlasov–Poisson Equations , 1984 .

[3]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[4]  Phillip Colella,et al.  A Particle-in-cell Method with Adaptive Phase-space Remapping for Kinetic Plasmas , 2011, SIAM J. Sci. Comput..

[5]  T. Liszka,et al.  The finite difference method at arbitrary irregular grids and its application in applied mechanics , 1980 .

[6]  Luis Gavete,et al.  An h-adaptive method in the generalized finite differences , 2003 .

[7]  Jean-Luc Vay,et al.  Mesh Refinement for Particle-In-Cell Plasma Simulations: Applications to - and benefits for - Heavy-Ion-Fusion , 2002 .

[8]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[9]  Rolf Dieter Grigorieff,et al.  Supraconvergence of a finite difference scheme for solutions in Hs(0, L) , 2005 .

[10]  Kun Zhou,et al.  Data-Parallel Octrees for Surface Reconstruction. , 2011, IEEE transactions on visualization and computer graphics.

[11]  Phillip Colella,et al.  Application of adaptive mesh refinement to particle-in-cell simulations of plasmas and beams , 2003 .

[12]  Roman Samulyak,et al.  An embedded boundary method for elliptic and parabolic problems with interfaces and application to multi-material systems with phase transitions , 2010 .

[13]  Luis Gavete,et al.  Influence of several factors in the generalized finite difference method , 2001 .

[14]  Phillip Colella,et al.  Controlling self-force errors at refinement boundaries for AMR-PIC , 2010, J. Comput. Phys..