J. W. Burby

2papers

2 Papers

81.0NAMay 20
A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

Mandela B. Quashie, J. W. Burby, Andrew J. Christlieb et al.

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

COMP-PHSep 9, 2015
Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational

C. L. Ellison, J. W. Burby, H. Qin

The Boris algorithm for integrating charged particle trajectories in electric and magnetic fields is popular due to its simple implementation, rapid iteration, and observed long-term numerical fidelity. The underlying cause of this long-term fidelity has become a matter of controversy, with one article claiming the method to be symplectic [S. D. Webb, J. Comput. Phys. 270 (2014) 570], and others claiming the method to be volume preserving but not symplectic [e.g. H. Qin et al., Phys. Plasmas 20 (2013) 084503]. To resolve the discrepancy, this letter leverages a discrete Helmholtz condition to demonstrate that no variational formulation of the Boris algorithm exists, indicating that the long-term fidelity should be attributed to the volume-preserving properties of the algorithm.