Abstract
A numerical algorithm for solving full gyro orbit of relativistic charged particle motion in magnetized plasmas is presented. The algorithm developed here achieves the following features simultaneously. (1) The time advancement is explicit, and (2) the integration is performed with respect to the observation time in a laboratory frame. (3) It is suitable for accurate long time integration with its volume preserving properties in the phase space. (4) The algorithm can properly treat the E×B drift velocity in electromagnetic fields for large relativistic factors, and (5) can be extended to arbitrary high orders with the aid of symmetric composition methods. Because our algorithm is formulated in the Lorentz covariant form, explicit conservation of the Minkowski norm is not assumed. Nevertheless, no secular growth of the numerical errors in the norm does occur, and its influence can be minimized up to the levels of round-off errors when a high-order scheme is applied. Numerical results are compared with explicit and implicit Runge–Kutta methods. The numerical accuracy and the computational efficiency are discussed for long time integration in toroidal magnetic field configuration.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.