TY - JOUR
AU - Siddi, Lorenzo
AU - Lapenta, Giovanni
AU - Gibbon, Paul
TI - Mesh-free Hamiltonian implementation of two dimensional Darwin model
JO - Physics of plasmas
VL - 24
IS - 8
SN - 1089-7674
CY - [S.l.]
PB - American Institute of Physics
M1 - FZJ-2017-05054
SP - 082103
PY - 2017
AB - 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.
LB - PUB:(DE-HGF)16
UR - <Go to ISI:>//WOS:000409232200015
DO - DOI:10.1063/1.4994705
UR - https://juser.fz-juelich.de/record/835922
ER -