Structure-preserving algorithms for guiding center dynamics based on the slow manifold of classical Pauli particle

作     者：张若涵 王正汹 肖建元 王丰 Ruohan ZHANG;Zhengxiong WANG;Jianyuan XIAO;Feng WANG

作者机构：Key Laboratory of Materials Modification by LaserIon and Electron BeamsSchool of PhysicsDalian University of TechnologyDalian 116024People’s Republic of China Department of Plasma Physics and Fusion EngineeringUniversity of Science and Technology of ChinaHefei 230026People’s Republic of China 

出 版 物：《Plasma Science and Technology》 (等离子体科学和技术（英文版）)

年 卷 期：2024年第26卷第6期

页      面：88-102页

核心收录：

学科分类：07[理学] 070204[理学-等离子体物理] 0702[理学-物理学] 

基　　金：supported by National Natural Science Foundation of China (Nos. 11975068 and 11925501) the National Key R&D Program of China (No. 2022YFE03090000) the Fundamental Research Funds for the Central Universities (No. DUT22ZD215) 

主　　题：structure-preserving algorithm averaged vector field classical Pauli particle guiding center dynamics 

摘      要：The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbits. Demonstrating significantly higher efficiency, this advanced method is capable of accomplishing the simulation in less than one-third of the time of directly computing the guiding center motion. In contrast to the CPP-based Boris algorithm, this approach inherits the advantages of the AVF method, yielding stable trajectories even achieved with a tenfold time step and reducing the energy error by two orders of magnitude. By comparing these two CPP algorithms with the traditional RK4 method, the numerical results indicate a remarkable performance in terms of both the computational efficiency and error elimination. Moreover, we verify the properties of slow manifold integrators and successfully observe the bounce on both sides of the limiting slow manifold with deliberately chosen perturbed initial conditions. To evaluate the practical value of the methods, we conduct simulations in non-axisymmetric perturbation magnetic fields as part of the experiments,demonstrating that our CPP-based AVF method can handle simulations under complex magnetic field configurations with high accuracy, which the CPP-based Boris algorithm lacks. Through numerical experiments, we demonstrate that the CPP can replace guiding center dynamics in using energy-preserving algorithms for computations, providing a new, efficient, as well as stable approach for applying structure-preserving algorithms in plasma simulations.