
-
Previous Article
Linear degree growth in lattice equations
- JCD Home
- This Issue
-
Next Article
Chains of rigid bodies and their numerical simulation by local frame methods
Study of adaptive symplectic methods for simulating charged particle dynamics
1. | LSEC, ICMSEC, Academy of Mathematics and Systems Science, CAS, Beijing 100190, China |
2. | School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China |
3. | Department of Engineering and Applied Physics, USTC, Hefei, Anhui 230026, China |
4. | Key Laboratory of Geospace Environment, CAS, Hefei, Anhui 230026, China |
In plasma simulations, numerical methods with high computational efficiency and long-term stability are needed. In this paper, symplectic methods with adaptive time steps are constructed for simulating the dynamics of charged particles under the electromagnetic field. With specifically designed step size functions, the motion of charged particles confined in a Penning trap under three different magnetic fields is studied, and also the dynamics of runaway electrons in tokamaks is investigated. The numerical experiments are performed to show the efficiency of the new derived adaptive symplectic methods.
References:
[1] |
G. Benettin and P. Sempio,
Adiabatic invariants and trapping of a point charge in a strong nonuniform magnetic field, Nonlinearity, 7 (1994), 281-303.
doi: 10.1088/0951-7715/7/1/014. |
[2] |
M. P. Calvo and J. M. Sanz-Serna,
The Development of variable-step symplectic integrators, with application to the two-body problem, SIAM Journal on Scientific Computing, 14 (1993), 936-952.
doi: 10.1137/0914057. |
[3] |
H. Dreicer,
Electron and ion runaway in a fully ionized gas, Physical Review, 115 (1959), 238-249.
doi: 10.1103/PhysRev.115.238. |
[4] |
B. Gladman, M. Duncan and J. Candy,
Symplectic integrators for long-term integrations in celestial mechanics, Celestial Mechanics & Dynamical Astronomy, 52 (1991), 221-240.
doi: 10.1007/BF00048485. |
[5] |
E. Hairer,
Variable time step integration with symplectic methods, Applied Numerical Mathematics, 25 (1997), 219-227.
doi: 10.1016/S0168-9274(97)00061-5. |
[6] |
E. Hairer, C. Lubich and G. Wanner, Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, Second edition, Springer Series in Computational Mathematics, 31. Springer-Verlag, Berlin, 2006. |
[7] |
E. Hairer and C. Lubich, Long-term analysis of a variational integrator for charged-particle dynamics in a strong magnetic field, preprint, Available from: https://na.uni-tuebingen.de/pub/lubich/papers/adiabatic.pdf. |
[8] |
Y. He, Y. J. Sun, J. Liu and H. Qin,
Volume-preserving algorithms for charged particle dynamics, Journal of Computational Physics, 281 (2015), 135-147.
doi: 10.1016/j.jcp.2014.10.032. |
[9] |
Y. He, Y. J. Sun, R. L. Zhang, Y. L. Wang, J. Liu and H. Qin, High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields, Physics of Plasmas, 23 (2016), 092109.
doi: 10.1063/1.4962677. |
[10] |
Y. He, Z. Q. Zhou, Y. J. Sun, J. Liu and H. Qin,
Explicit $K$-symplectic algorithms for charged particle dynamics, Physics Letters A, 381 (2017), 568-573.
doi: 10.1016/j.physleta.2016.12.031. |
[11] |
W. Z. Huang and B. Leimkuhler,
The adaptive verlet method, SIAM Journal on Scientific Computing, 18 (1997), 239-256.
doi: 10.1137/S1064827595284658. |
[12] |
J. D. Jackson,
From Lorenz to Coulomb and other explicit gauge transformations, American Journal of Physics, 70 (2002), 917-928.
doi: 10.1119/1.1491265. |
[13] |
C. Knapp, A. Kendl, A. Koskela and A. Ostermann, Splitting methods for time integration of trajectories in combined electric and magnetic fields, Physical Review E, 92 (2015), 063310, 13 pp. |
[14] |
M. Kretzschmar,
Single particle motion in a Penning trap: Description in the classical canonical formalism, Physica Scripta, 46 (1992), 544-554.
doi: 10.1088/0031-8949/46/6/011. |
[15] |
P. Langevin,
Sur la théorie du mouvement brownien, C. R. Acad. Sci. (Paris), 146 (1908), 530-533.
|
[16] |
J. Liu, H. Qin, Y. Wang and et al., Largest particle simulations downgrade the runaway electron risk for Iter, arXiv: 1611.02362. |
[17] |
J. E. Marsden and T. S. Ratiu, Introduction to Mechanics and Symmetry. A Basic Exposition of Classical Mechanical Systems, Texts in Applied Mathematics, 17. Springer-Verlag, New York, 1994.
doi: 10.1007/978-1-4612-2682-6. |
[18] |
H. Qin and R. C. Davidson, An exact magnetic-moment invariant of charged-particle gyromotion, Physical Review Letters, 96 (2006), 085003.
doi: 10.1103/PhysRevLett.96.085003. |
[19] |
H. Qin and X. Y. Guan, Variational symplectic integrator for long-time simulations of the guiding-center motion of charged particles in general magnetic fields, Physical Review Letters, 100 (2008), 035006.
doi: 10.1103/PhysRevLett.100.035006. |
[20] |
H. Qin, X. Y. Guan and W. M. Tang, Variational symplectic algorithm for guiding center dynamics and its application in tokamak geometry, Physics of Plasmas, 16 (2009), 042510.
doi: 10.1063/1.3099055. |
[21] |
S. Reich,
Backward error analysis for numerical integrators, SIAM Journal on Numerical Analysis, 36 (1999), 1549-1570.
doi: 10.1137/S0036142997329797. |
[22] |
A. S. Richardson and J. M. Finn,
Symplectic integrators with adaptive time steps, Plasma Physics and Controlled Fusion, 54 (2012), 96-100.
doi: 10.1088/0741-3335/54/1/014004. |
[23] |
C. C. Rodegheri, K. Blaum, H. Kracke, S. Kreim, A. Mooser, W. Quint, S. Ulmer and J. Walz, An experiment for the direct determination of the g-factor of a single proton in a Penning trap, New Journal of Physics, 14 (2012), 063011.
doi: 10.1088/1367-2630/14/6/063011. |
[24] |
J. Schmitt and M. Leok, Adaptive variational integrators, arXiv: 1709.01975. |
[25] |
Y. Y. Shi, Y. J. Sun, Y. He, H. Qin and J. Liu,
Symplectic integrators with adaptive time step applied to runaway electron dynamics, Numerical Algorithms, 81 (2019), 1295-1309.
doi: 10.1007/s11075-018-0636-6. |
[26] |
M. Toggweiler, A. Adelmann, P. Arbenz and J. J. Yang,
A novel adaptive time stepping variant of the Boris–Buneman integrator for the simulation of particle accelerators with space charge, Journal of Computational Physics, 273 (2014), 255-267.
doi: 10.1016/j.jcp.2014.05.008. |
[27] |
Y. L. Wang, J. Liu and H. Qin, Lorentz covariant canonical symplectic algorithms for dynamics of charged particles, Physics of Plasmas, 23 (2016), 122513.
doi: 10.1063/1.4972824. |
[28] |
R. L. Zhang, Y. L. Wang, Y. He, J. Y. Xiao, J. Liu, H. Qin and Y. F. Tang, Explicit symplectic algorithms based on generating functions for relativistic charged particle dynamics in time-dependent electromagnetic field, Physics of Plasmas, 25 (2018), 022117.
doi: 10.1063/1.5012767. |
[29] |
Z. Q. Zhou, Y. He, Y. J. Sun, J. Liu and H. Qin, Explicit symplectic methods for solving charged particle trajectories, Physics of Plasmas, 24 (2017), 052507.
doi: 10.1063/1.4982743. |
show all references
References:
[1] |
G. Benettin and P. Sempio,
Adiabatic invariants and trapping of a point charge in a strong nonuniform magnetic field, Nonlinearity, 7 (1994), 281-303.
doi: 10.1088/0951-7715/7/1/014. |
[2] |
M. P. Calvo and J. M. Sanz-Serna,
The Development of variable-step symplectic integrators, with application to the two-body problem, SIAM Journal on Scientific Computing, 14 (1993), 936-952.
doi: 10.1137/0914057. |
[3] |
H. Dreicer,
Electron and ion runaway in a fully ionized gas, Physical Review, 115 (1959), 238-249.
doi: 10.1103/PhysRev.115.238. |
[4] |
B. Gladman, M. Duncan and J. Candy,
Symplectic integrators for long-term integrations in celestial mechanics, Celestial Mechanics & Dynamical Astronomy, 52 (1991), 221-240.
doi: 10.1007/BF00048485. |
[5] |
E. Hairer,
Variable time step integration with symplectic methods, Applied Numerical Mathematics, 25 (1997), 219-227.
doi: 10.1016/S0168-9274(97)00061-5. |
[6] |
E. Hairer, C. Lubich and G. Wanner, Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, Second edition, Springer Series in Computational Mathematics, 31. Springer-Verlag, Berlin, 2006. |
[7] |
E. Hairer and C. Lubich, Long-term analysis of a variational integrator for charged-particle dynamics in a strong magnetic field, preprint, Available from: https://na.uni-tuebingen.de/pub/lubich/papers/adiabatic.pdf. |
[8] |
Y. He, Y. J. Sun, J. Liu and H. Qin,
Volume-preserving algorithms for charged particle dynamics, Journal of Computational Physics, 281 (2015), 135-147.
doi: 10.1016/j.jcp.2014.10.032. |
[9] |
Y. He, Y. J. Sun, R. L. Zhang, Y. L. Wang, J. Liu and H. Qin, High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields, Physics of Plasmas, 23 (2016), 092109.
doi: 10.1063/1.4962677. |
[10] |
Y. He, Z. Q. Zhou, Y. J. Sun, J. Liu and H. Qin,
Explicit $K$-symplectic algorithms for charged particle dynamics, Physics Letters A, 381 (2017), 568-573.
doi: 10.1016/j.physleta.2016.12.031. |
[11] |
W. Z. Huang and B. Leimkuhler,
The adaptive verlet method, SIAM Journal on Scientific Computing, 18 (1997), 239-256.
doi: 10.1137/S1064827595284658. |
[12] |
J. D. Jackson,
From Lorenz to Coulomb and other explicit gauge transformations, American Journal of Physics, 70 (2002), 917-928.
doi: 10.1119/1.1491265. |
[13] |
C. Knapp, A. Kendl, A. Koskela and A. Ostermann, Splitting methods for time integration of trajectories in combined electric and magnetic fields, Physical Review E, 92 (2015), 063310, 13 pp. |
[14] |
M. Kretzschmar,
Single particle motion in a Penning trap: Description in the classical canonical formalism, Physica Scripta, 46 (1992), 544-554.
doi: 10.1088/0031-8949/46/6/011. |
[15] |
P. Langevin,
Sur la théorie du mouvement brownien, C. R. Acad. Sci. (Paris), 146 (1908), 530-533.
|
[16] |
J. Liu, H. Qin, Y. Wang and et al., Largest particle simulations downgrade the runaway electron risk for Iter, arXiv: 1611.02362. |
[17] |
J. E. Marsden and T. S. Ratiu, Introduction to Mechanics and Symmetry. A Basic Exposition of Classical Mechanical Systems, Texts in Applied Mathematics, 17. Springer-Verlag, New York, 1994.
doi: 10.1007/978-1-4612-2682-6. |
[18] |
H. Qin and R. C. Davidson, An exact magnetic-moment invariant of charged-particle gyromotion, Physical Review Letters, 96 (2006), 085003.
doi: 10.1103/PhysRevLett.96.085003. |
[19] |
H. Qin and X. Y. Guan, Variational symplectic integrator for long-time simulations of the guiding-center motion of charged particles in general magnetic fields, Physical Review Letters, 100 (2008), 035006.
doi: 10.1103/PhysRevLett.100.035006. |
[20] |
H. Qin, X. Y. Guan and W. M. Tang, Variational symplectic algorithm for guiding center dynamics and its application in tokamak geometry, Physics of Plasmas, 16 (2009), 042510.
doi: 10.1063/1.3099055. |
[21] |
S. Reich,
Backward error analysis for numerical integrators, SIAM Journal on Numerical Analysis, 36 (1999), 1549-1570.
doi: 10.1137/S0036142997329797. |
[22] |
A. S. Richardson and J. M. Finn,
Symplectic integrators with adaptive time steps, Plasma Physics and Controlled Fusion, 54 (2012), 96-100.
doi: 10.1088/0741-3335/54/1/014004. |
[23] |
C. C. Rodegheri, K. Blaum, H. Kracke, S. Kreim, A. Mooser, W. Quint, S. Ulmer and J. Walz, An experiment for the direct determination of the g-factor of a single proton in a Penning trap, New Journal of Physics, 14 (2012), 063011.
doi: 10.1088/1367-2630/14/6/063011. |
[24] |
J. Schmitt and M. Leok, Adaptive variational integrators, arXiv: 1709.01975. |
[25] |
Y. Y. Shi, Y. J. Sun, Y. He, H. Qin and J. Liu,
Symplectic integrators with adaptive time step applied to runaway electron dynamics, Numerical Algorithms, 81 (2019), 1295-1309.
doi: 10.1007/s11075-018-0636-6. |
[26] |
M. Toggweiler, A. Adelmann, P. Arbenz and J. J. Yang,
A novel adaptive time stepping variant of the Boris–Buneman integrator for the simulation of particle accelerators with space charge, Journal of Computational Physics, 273 (2014), 255-267.
doi: 10.1016/j.jcp.2014.05.008. |
[27] |
Y. L. Wang, J. Liu and H. Qin, Lorentz covariant canonical symplectic algorithms for dynamics of charged particles, Physics of Plasmas, 23 (2016), 122513.
doi: 10.1063/1.4972824. |
[28] |
R. L. Zhang, Y. L. Wang, Y. He, J. Y. Xiao, J. Liu, H. Qin and Y. F. Tang, Explicit symplectic algorithms based on generating functions for relativistic charged particle dynamics in time-dependent electromagnetic field, Physics of Plasmas, 25 (2018), 022117.
doi: 10.1063/1.5012767. |
[29] |
Z. Q. Zhou, Y. He, Y. J. Sun, J. Liu and H. Qin, Explicit symplectic methods for solving charged particle trajectories, Physics of Plasmas, 24 (2017), 052507.
doi: 10.1063/1.4982743. |



















Quantities | Symbols | Non-relativistic | Relativistic | |
Units | Units | |||
Time | ||||
Position | ||||
Velocity | ||||
Momentum | ||||
Canonical Momentum | ||||
Electric field | ||||
Magnetic field | ||||
Vector field | ||||
Scalar field | ||||
Hamiltonian | ||||
Quantities | Symbols | Non-relativistic | Relativistic | |
Units | Units | |||
Time | ||||
Position | ||||
Velocity | ||||
Momentum | ||||
Canonical Momentum | ||||
Electric field | ||||
Magnetic field | ||||
Vector field | ||||
Scalar field | ||||
Hamiltonian | ||||
period | 1st | 70th | 120th | 150th |
fixed | 100 | 451 | 2134 | 5470 |
adaptive | 100 | 100 | 100 | 100 |
period | 1st | 70th | 120th | 150th |
fixed | 100 | 451 | 2134 | 5470 |
adaptive | 100 | 100 | 100 | 100 |
[1] |
Yuezheng Gong, Jiaquan Gao, Yushun Wang. High order Gauss-Seidel schemes for charged particle dynamics. Discrete and Continuous Dynamical Systems - B, 2018, 23 (2) : 573-585. doi: 10.3934/dcdsb.2018034 |
[2] |
Mihai Bostan. On the Boltzmann equation for charged particle beams under the effect of strong magnetic fields. Discrete and Continuous Dynamical Systems - B, 2015, 20 (2) : 339-371. doi: 10.3934/dcdsb.2015.20.339 |
[3] |
Mathieu Lutz. Application of Lie transform techniques for simulation of a charged particle beam. Discrete and Continuous Dynamical Systems - S, 2015, 8 (1) : 185-221. doi: 10.3934/dcdss.2015.8.185 |
[4] |
Raffaele D'Ambrosio, Martina Moccaldi, Beatrice Paternoster. Numerical preservation of long-term dynamics by stochastic two-step methods. Discrete and Continuous Dynamical Systems - B, 2018, 23 (7) : 2763-2773. doi: 10.3934/dcdsb.2018105 |
[5] |
Lijin Wang, Jialin Hong. Generating functions for stochastic symplectic methods. Discrete and Continuous Dynamical Systems, 2014, 34 (3) : 1211-1228. doi: 10.3934/dcds.2014.34.1211 |
[6] |
Per Christian Moan, Jitse Niesen. On an asymptotic method for computing the modified energy for symplectic methods. Discrete and Continuous Dynamical Systems, 2014, 34 (3) : 1105-1120. doi: 10.3934/dcds.2014.34.1105 |
[7] |
Alexandre Mouton. Two-scale semi-Lagrangian simulation of a charged particle beam in a periodic focusing channel. Kinetic and Related Models, 2009, 2 (2) : 251-274. doi: 10.3934/krm.2009.2.251 |
[8] |
Tomáš Roubíček, Giuseppe Tomassetti. Dynamics of charged elastic bodies under diffusion at large strains. Discrete and Continuous Dynamical Systems - B, 2020, 25 (4) : 1415-1437. doi: 10.3934/dcdsb.2019234 |
[9] |
Andrew J. Steyer, Erik S. Van Vleck. Underlying one-step methods and nonautonomous stability of general linear methods. Discrete and Continuous Dynamical Systems - B, 2018, 23 (7) : 2859-2877. doi: 10.3934/dcdsb.2018108 |
[10] |
Sören Bartels, Jakob Keck. Adaptive time stepping in elastoplasticity. Discrete and Continuous Dynamical Systems - S, 2021, 14 (1) : 71-88. doi: 10.3934/dcdss.2020323 |
[11] |
Qifeng Cheng, Xue Han, Tingting Zhao, V S Sarma Yadavalli. Improved particle swarm optimization and neighborhood field optimization by introducing the re-sampling step of particle filter. Journal of Industrial and Management Optimization, 2019, 15 (1) : 177-198. doi: 10.3934/jimo.2018038 |
[12] |
Cesare Bracco, Annalisa Buffa, Carlotta Giannelli, Rafael Vázquez. Adaptive isogeometric methods with hierarchical splines: An overview. Discrete and Continuous Dynamical Systems, 2019, 39 (1) : 241-261. doi: 10.3934/dcds.2019010 |
[13] |
Giacomo Frassoldati, Luca Zanni, Gaetano Zanghirati. New adaptive stepsize selections in gradient methods. Journal of Industrial and Management Optimization, 2008, 4 (2) : 299-312. doi: 10.3934/jimo.2008.4.299 |
[14] |
Marie-Claude Arnaud. When are the invariant submanifolds of symplectic dynamics Lagrangian?. Discrete and Continuous Dynamical Systems, 2014, 34 (5) : 1811-1827. doi: 10.3934/dcds.2014.34.1811 |
[15] |
Guillermo Dávila-Rascón, Yuri Vorobiev. Hamiltonian structures for projectable dynamics on symplectic fiber bundles. Discrete and Continuous Dynamical Systems, 2013, 33 (3) : 1077-1088. doi: 10.3934/dcds.2013.33.1077 |
[16] |
Pablo G. Barrientos, Artem Raibekas. Robustly non-hyperbolic transitive symplectic dynamics. Discrete and Continuous Dynamical Systems, 2018, 38 (12) : 5993-6013. doi: 10.3934/dcds.2018259 |
[17] |
Fasma Diele, Carmela Marangi. Positive symplectic integrators for predator-prey dynamics. Discrete and Continuous Dynamical Systems - B, 2018, 23 (7) : 2661-2678. doi: 10.3934/dcdsb.2017185 |
[18] |
Yoonsang Lee, Bjorn Engquist. Variable step size multiscale methods for stiff and highly oscillatory dynamical systems. Discrete and Continuous Dynamical Systems, 2014, 34 (3) : 1079-1097. doi: 10.3934/dcds.2014.34.1079 |
[19] |
Tingting Wu, Yufei Yang, Huichao Jing. Two-step methods for image zooming using duality strategies. Numerical Algebra, Control and Optimization, 2014, 4 (3) : 209-225. doi: 10.3934/naco.2014.4.209 |
[20] |
Angelamaria Cardone, Dajana Conte, Beatrice Paternoster. Two-step collocation methods for fractional differential equations. Discrete and Continuous Dynamical Systems - B, 2018, 23 (7) : 2709-2725. doi: 10.3934/dcdsb.2018088 |
Impact Factor:
Tools
Metrics
Other articles
by authors
[Back to Top]