%0 Journal Article
%A Siddi, Lorenzo
%A Lapenta, Giovanni
%A Gibbon, Paul
%T Mesh-free Hamiltonian implementation of two dimensional Darwin model
%J Physics of plasmas
%V 24
%N 8
%@ 1089-7674
%C [S.l.]
%I American Institute of Physics
%M FZJ-2017-05054
%P 082103
%D 2017
%X A new approach to Darwin or magnetoinductive plasma simulation is presented, which combines amesh-free field solver with a robust time-integration scheme avoiding numerical divergence errorsin the solenoidal field components. The mesh-free formulation employs an efficient parallelBarnes-Hut tree algorithm to speed up the computation of fields summed directly from the particles,avoiding the necessity of divergence cleaning procedures typically required by particle-incellmethods. The time-integration scheme employs a Hamiltonian formulation of the Lorentzforce, circumventing the development of violent numerical instabilities associated with time differentiationof the vector potential. It is shown that a semi-implicit scheme converges rapidly and isrobust to further numerical instabilities which can develop from a dominant contribution of the vectorpotential to the canonical momenta. The model is validated by various static and dynamicbenchmark tests, including a simulation of the Weibel-like filamentation instability in beam-plasmainteractions.
%F PUB:(DE-HGF)16
%9 Journal Article
%U <Go to ISI:>//WOS:000409232200015
%R 10.1063/1.4994705
%U https://juser.fz-juelich.de/record/835922