Efficient dual space source interpolation method for the numerical solution of self-consistent static beam-wave interactions

Many numerical solvers model the general nonlinear problem of particle motion in self-generated fields through the computation of particle trajectories and fields in a Lagrangian and Eulerian framework, respectively. The particle-weighting step mapping the source terms between these frameworks is a primary source of noise and computation time. Herein we eliminate this costly step by treating the trajectories as bounds on the particle flow, applying charge conservation to interpolate the source distribution. The result is an order of magnitude improvement in accuracy and an 80 to 120 fold speed-up.