An arbitrary order time-stepping algorithm for tracking particles in inhomogeneous magnetic fields
From MaRDI portal
Publication:6145384
Abstract: The Lorentz equations describe the motion of electrically charged particles in electric and magnetic fields and are used widely in plasma physics. The most popular numerical algorithm for solving them is the Boris method, a variant of the St"ormer-Verlet algorithm. Boris' method is phase space volume conserving and simulated particles typically remain near the correct trajectory. However, it is only second order accurate. Therefore, in scenarios where it is not enough to know that a particle stays on the right trajectory but one needs to know where on the trajectory the particle is at a given time, Boris method requires very small time steps to deliver accurate phase information, making it computationally expensive. We derive an improved version of the high-order Boris spectral deferred correction algorithm (Boris-SDC) by adopting a convergence acceleration strategy for second order problems based on the Generalised Minimum Residual (GMRES) method. Our new algorithm is easy to implement as it still relies on the standard Boris method. Like Boris-SDC it can deliver arbitrary order of accuracy through simple changes of runtime parameter but possesses better long-term energy stability. We demonstrate for two examples, a magnetic mirror trap and the Solev'ev equilibrium, that the new method can deliver better accuracy at lower computational cost compared to the standard Boris method. While our examples are motivated by tracking ions in the magnetic field of a nuclear fusion reactor, the introduced algorithm can potentially deliver similar improvements in efficiency for other applications.
Cites work
- A high-order boris integrator
- Accelerating the convergence of spectral deferred correction methods
- Achieving Brouwer's law with implicit Runge-Kutta methods
- Arbitrary order Krylov deferred correction methods for differential algebraic equations
- Common molecular dynamics algorithms revisited: Accuracy and optimal time steps of Störmer-Leapfrog integrators
- Energy behaviour of the Boris method for charged-particle dynamics
- Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
- Faster SDC convergence on non-equidistant grids by DIRK sweeps
- Higher order volume-preserving schemes for charged particle dynamics
- scientific article; zbMATH DE number 3143576 (Why is no real title available?)
- scientific article; zbMATH DE number 1745051 (Why is no real title available?)
- scientific article; zbMATH DE number 940566 (Why is no real title available?)
- scientific article; zbMATH DE number 781821 (Why is no real title available?)
- Parallel methods for ordinary differential equations
- Parallelizing spectral deferred corrections across the method
- Short pulse laser interactions with matter. An introduction.
- Solving Ordinary Differential Equations I
- Spectral deferred correction methods for ordinary differential equations
- Symmetric multistep methods for charged-particle dynamics
- Volume-preserving algorithms for charged particle dynamics
Cited in
(3)
This page was built for publication: An arbitrary order time-stepping algorithm for tracking particles in inhomogeneous magnetic fields
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6145384)
