arrow
Volume 24, Issue 3
Comparison of Preconditioning Strategies in Energy Conserving Implicit Particle in Cell Methods

Lorenzo Siddi, Emanuele Cazzola & Giovanni Lapenta

Commun. Comput. Phys., 24 (2018), pp. 672-694.

Published online: 2018-05

Export citation
  • Abstract

This work presents a set of preconditioning strategies able to significantly accelerate the performance of fully implicit energy-conserving Particle-in-Cell methods to a level that becomes competitive with semi-implicit methods. We compare three different preconditioners. We consider three methods and compare them with a straight unpreconditioned Jacobian Free Newton Krylov (JFNK) implementation. The first two focus, respectively, on improving the handling of particles (particle hiding) or fields (field hiding) within the JFNK iteration. The third uses the field hiding preconditioner within a direct Newton iteration where a Schwarz-decomposed Jacobian is computed analytically. Clearly, field hiding used with JFNK or with the direct Newton-Schwarz (DNS) method outperforms all method. We compare these implementations with a recent semi-implicit energy conserving scheme. Fully implicit methods are still lag behind in cost per cycle but not by a large margin when proper preconditioning is used. However, for exact energy conservation, preconditioned fully implicit methods are significantly easier to implement compared with semi-implicit methods and can be extended to fully relativistic physics.

  • AMS Subject Headings

65Z05, 65F08, 65M75

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-24-672, author = {Lorenzo Siddi, Emanuele Cazzola and Giovanni Lapenta}, title = {Comparison of Preconditioning Strategies in Energy Conserving Implicit Particle in Cell Methods}, journal = {Communications in Computational Physics}, year = {2018}, volume = {24}, number = {3}, pages = {672--694}, abstract = {

This work presents a set of preconditioning strategies able to significantly accelerate the performance of fully implicit energy-conserving Particle-in-Cell methods to a level that becomes competitive with semi-implicit methods. We compare three different preconditioners. We consider three methods and compare them with a straight unpreconditioned Jacobian Free Newton Krylov (JFNK) implementation. The first two focus, respectively, on improving the handling of particles (particle hiding) or fields (field hiding) within the JFNK iteration. The third uses the field hiding preconditioner within a direct Newton iteration where a Schwarz-decomposed Jacobian is computed analytically. Clearly, field hiding used with JFNK or with the direct Newton-Schwarz (DNS) method outperforms all method. We compare these implementations with a recent semi-implicit energy conserving scheme. Fully implicit methods are still lag behind in cost per cycle but not by a large margin when proper preconditioning is used. However, for exact energy conservation, preconditioned fully implicit methods are significantly easier to implement compared with semi-implicit methods and can be extended to fully relativistic physics.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2017-0171}, url = {http://global-sci.org/intro/article_detail/cicp/12276.html} }
TY - JOUR T1 - Comparison of Preconditioning Strategies in Energy Conserving Implicit Particle in Cell Methods AU - Lorenzo Siddi, Emanuele Cazzola & Giovanni Lapenta JO - Communications in Computational Physics VL - 3 SP - 672 EP - 694 PY - 2018 DA - 2018/05 SN - 24 DO - http://doi.org/10.4208/cicp.OA-2017-0171 UR - https://global-sci.org/intro/article_detail/cicp/12276.html KW - Preconditioners in Particle-in-Cell, energy conserving Particle-in-Cell, plasma numerics. AB -

This work presents a set of preconditioning strategies able to significantly accelerate the performance of fully implicit energy-conserving Particle-in-Cell methods to a level that becomes competitive with semi-implicit methods. We compare three different preconditioners. We consider three methods and compare them with a straight unpreconditioned Jacobian Free Newton Krylov (JFNK) implementation. The first two focus, respectively, on improving the handling of particles (particle hiding) or fields (field hiding) within the JFNK iteration. The third uses the field hiding preconditioner within a direct Newton iteration where a Schwarz-decomposed Jacobian is computed analytically. Clearly, field hiding used with JFNK or with the direct Newton-Schwarz (DNS) method outperforms all method. We compare these implementations with a recent semi-implicit energy conserving scheme. Fully implicit methods are still lag behind in cost per cycle but not by a large margin when proper preconditioning is used. However, for exact energy conservation, preconditioned fully implicit methods are significantly easier to implement compared with semi-implicit methods and can be extended to fully relativistic physics.

Lorenzo Siddi, Emanuele Cazzola and Giovanni Lapenta. (2018). Comparison of Preconditioning Strategies in Energy Conserving Implicit Particle in Cell Methods. Communications in Computational Physics. 24 (3). 672-694. doi:10.4208/cicp.OA-2017-0171
Copy to clipboard
The citation has been copied to your clipboard