中国物理B ›› 2021, Vol. 30 ›› Issue (5): 55203-055203.doi: 10.1088/1674-1056/abd161
Energy behavior of Boris algorithm
Abdullah Zafar1,† and Majid Khan2
- 1 Department of Engineering and Applied Physics, University of Science and Technology of China, Hefei 230026, China;
2 Department of Physics, Quaid-i-Azam University, Islamabad 45320, Pakistan
Energy behavior of Boris algorithm
Abdullah Zafar1,? and Majid Khan2
- 1 Department of Engineering and Applied Physics, University of Science and Technology of China, Hefei 230026, China;
2 Department of Physics, Quaid-i-Azam University, Islamabad 45320, Pakistan
摘要： Boris numerical scheme due to its long-time stability, accuracy and conservative properties has been widely applied in many studies of magnetized plasmas. Such algorithms conserve the phase space volume and hence provide accurate charge particle orbits. However, this algorithm does not conserve the energy in some special electromagnetic configurations, particularly for long simulation times. Here, we empirically analyze the energy behavior of Boris algorithm by applying it to a 2D autonomous Hamiltonian. The energy behavior of the Boris method is found to be strongly related to the integrability of our Hamiltonian system. We find that if the invariant tori is preserved under Boris discretization, the energy error can be bounded for an exponentially long time, otherwise the said error will show a linear growth. On the contrary, for a non-integrable Hamiltonian system, a random walk pattern has been observed in the energy error.
中图分类号: (Magnetized plasmas)
- 52.25.Xz