This title appears in the Scientific Report :
2017
Please use the identifier:
http://hdl.handle.net/2128/18130 in citations.
Please use the identifier: http://dx.doi.org/10.1063/1.4994705 in citations.
Mesh-free Hamiltonian implementation of two dimensional Darwin model
Mesh-free Hamiltonian implementation of two dimensional Darwin model
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 alg...
Saved in:
Personal Name(s): | Siddi, Lorenzo |
---|---|
Lapenta, Giovanni / Gibbon, Paul (Corresponding author) | |
Contributing Institute: |
JARA - HPC; JARA-HPC Jülich Supercomputing Center; JSC |
Published in: | Physics of plasmas, 24 (2017) 8, S. 082103 |
Imprint: |
[S.l.]
American Institute of Physics
2017
|
DOI: |
10.1063/1.4994705 |
Document Type: |
Journal Article |
Research Program: |
Kinetic Plasma Simulation with Highly Scalable Particle Codes Computational Science and Mathematical Methods |
Link: |
Published on 2017-07-20. Available in OpenAccess from 2018-07-20. Published on 2017-07-20. Available in OpenAccess from 2018-07-20. Published on 2017-07-20. Available in OpenAccess from 2018-07-20. |
Publikationsportal JuSER |
Please use the identifier: http://dx.doi.org/10.1063/1.4994705 in citations.
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. |