Mesh-free Hamiltonian implementation of two dimensional Darwin model

Siddi, Lorenzo; Lapenta, Giovanni; Gibbon, Paul

doi:10.1063/1.4994705

Journal Article

FZJ-2017-05054

Mesh-free Hamiltonian implementation of two dimensional Darwin model

Siddi, L.Extern* ; Lapenta, G. ; Gibbon, P. (Corresponding author)FZJ*

2017
American Institute of Physics [S.l.]

Physics of plasmas 24(8), 082103 (2017) [10.1063/1.4994705]

This record in other databases:

Please use a persistent id in citations: http://hdl.handle.net/2128/18130 doi:10.1063/1.4994705

Abstract: 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.

Classification: