Development of momentum and energy conserving PIC code with nonuniform 2D and 3D Cartesian grids

  • Vedoucí práce / Supervisor: Dr. David Tskhakaya
  • Pracoviště / Workplace: ÚFP AV ČR
  • Kontakt / Contact:


The PIC (Particle-in Cell) codes represent powerful tools for plasma transport study [1,2]. They are frequently used in fusion plasma applications and contribute to design of next generation fusion devices. As it was demonstrated in [3], PIC simulations with Cartesian grid do conserve momentum and energy simultaneously only for uniform grids. This requirement represents very heavy condition for plasma edge modelling, when nonuniform space resolution is favorable. Some time ago a solution for 1D codes has been found [4], but so far up to now no multidimensional generalization of this approach has been done.

The aim of this PhD work is to develop a multidimensional energy and momentum conserving PIC code with nonuniform Cartesian grid. We expected significant speedup, by order of magnitude, of such code against the conventional codes, which will save tens of millions of CPU hours per year for average PIC user on high performance computer facilities.



[1] D. Tskhakaya, Plasma Phys. Cont. Fus., 59 (2017) 114001

[2] T. Takizuka, et al., Contrib. Plasma Phys., 59 (2017) 034008

[3] D. Tskhakaya, et al., Contrib. Plasma Phys., 47 (2007) 563

[4] J. Duras, et al., Contrib. Plasma Phys., 54 (2014) 697