Conservative integrators for vortex blob methods on the plane

Conservative integrators have recently been derived for many-body systems using the Discrete Multiplier Method, specifically for the point vortex problems in [1]. In the present paper, we extend this work to derive conservative symmetric second-order one-step integrators for a family of vortex blob models. A rational function approximation was used to approximate the exponential integral which appears in the Hamiltonian. Conservative properties and second order convergence are proved. Numerical experiments are shown to verify the conservative property of these integrators, their second-order accuracy, and as well as the resulting spatial and temporal accuracy of the vortex blob method. Moreover, the derived implicit conservative integrators are shown to better preserve conserved quantities than standard higher-order explicit integrators on comparable computation times.