The anatomy of Boris type solvers and the Lie operator formalism for deriving large time-step magnetic field integrators
| Autor: | Chin, Siu A., Cator, Durward |
|---|---|
| Rok vydání: | 2021 |
| Předmět: | |
| Druh dokumentu: | Working Paper |
| DOI: | 10.1016/j.jcp.2022.111422 |
| Popis: | This work gives a Lie operator derivation of various Boris solvers via a detailed study of trajectory errors in a constant magnetic field. These errors in the gyrocenter location and the gyroradius are the foundational basis for why Boris solvers existed, independent of any finite-difference schemes. This work shows that there are two distinct ways of eliminating these errors so that the trajectory of a charged particle in a constant magnetic field is exactly on the cyclotron orbit. One way reproduces the known second-order symmetric Boris solver. The other yields a previously unknown, but also on-orbit solver, not derivable from finite-difference schemes. By revisiting some historical calculations, it is found that many publications do not distinguish the poorly behaved leap-frog Boris solver from the symmetric second-order Boris algorithm. This symmetric second-order Boris solver's trajectory is much more accurate and remains close to the exact orbit in a combined $nonuniform$ electric and magnetic field at time-steps greater than the cyclotron period. Finally, this operator formalism showed that Buneman's cycloid fitting scheme is mathematically identical to Boris' on-orbit solver and that Boris' E-B splitting is unnecessary. Comment: Completely revised with 32 pages and 18 figures |
| Databáze: | arXiv |
| Externí odkaz: |
|