Numerical integration of particle orbits in discontinuous fields using VENUS-LEVIS and SPEC

The orbit code VENUS-LEVIS is adapted to follow particles in plasma equilibria with discontinuous fields generated by the Stepped Pressure Equilibrium Code (SPEC). The latter is an implementation of the Multi-Region relaxed MHD model, which efficiently computes Taylor states in a series of nested toroidal volumes and supports the formation of magnetic islands and chaotic field regions. To adapt VENUS-LEVIS, an event location procedure is implemented in the existing numerical integrator, which ensures the particle sees the correct field along its trajectory, regardless of the discontinuities present in the Stepped Pressure Equilibrium model. The algorithm is tested in the case where the magnetic field is uniform in the upper and lower half-spaces but has a discontinuity in its direction (shear) on the plane z=0. Particle drifts due to the discontinuity are studied. The convergence properties of the numerical scheme are highlighted by the numerical accuracy, and conservation of the system's invariants, such as energy and momentum. Simulations and convergence studies using the SPEC-LEVIS interface in axisymmetric geometry are then presented. Finally, illustrative particle drifts due to the discontinuity are studied and explained: we examine drifts associated with the change in Larmor radius of passing particles with small excursion from flux surfaces.