Issues in measure-preserving three dimensional flow integrators: Self-adjointness, reversibility, and non-uniform time stepping

Abstract

Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a "special divergence-free" (SDF) property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volumepreserving flows. We investigate two potential schemes, both based on the general method of Feng and Shang [Numer. Math. 71, 451 (1995)], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Richardson and Finn [Plasma Phys. Controlled Fusion 54, 014004 (2012)], appears to work very well. V C 2015 AIP Publishing LLC. [http://dx.

Improved accuracy in degenerate variational integrators for guiding center and magnetic field line flow

“…If a discrete scheme T h also satisfies this map reversibility, it should have the favorable properties due to reversibility 16 . In fact, neither T h nor T † h (where adjoint is defined as T † h = T −1 −h ; see Ref.…”

“…12) satisfy this map reversibility property. However, T h does satisfy the related property weak reversibility 16 ,…”

“…) is also weakly reversible, and because it is self-adjoint, it is also reversible 16 . This map reversibility appears to be responsible for the observed good long-time behavior.…”

Analysis of the Chaotic Behavior of the Lower Hybrid Wave Propagation in Magnetised Plasma by Hamiltonian Theory

“…In particular, Figure 1 of [34] is analogous to our Figure 10. More recently, several different algorithms have been developed for general applications [35,36], and in particular, they have been applied to plasma-physics-related problems (see [37]). For our purposes, we chose a second order implicit symplectic algorithm inspired by Channel and Scovel.…”

“…Numerical analysis by means of a symplectic algorithm [34], implemented in a ray-tracing code, support the analytical considerations. It is worth remarking here that accurate symplectic algorithms have been recently implemented (see, for example, [35,36]), and they have been applied to the study of charged particles' motion in a plasma (see [37]). …”