2D-2P multiscale semi-Lagrangian algorithm for fast electron transport in the relativistic Vlasov-Fokker-Planck equation

For sufficiently strong electric fields, electrons may break away from thermal equilibrium and approach relativistic speeds. Such “runaway” electrons, common in tokamak disruptions, traverse orbits at much faster time scales than collisional ones, while dynamics of interest saturate on time scales much longer than these. In this study, we propose a 2D-2P semi-Lagrangian scheme to efficiently bridge these temporal scales. The approach reformulates the Vlasov equation as an integro-differential operator using Green’s functions along electron orbits, and employs operator splitting to decouple integrals over the relativistic collisional source [1,2]. We consider 2D magnetic fields, but the formulation generalizes to arbitrary 3D magnetic fields. The proposed 2D-2P treatment is formally first-order accurate in time, but (i) preserves asymptotic properties associated with stiff Vlasov term, (ii) is uniformly accurate in Δt/ε, where Δt is the timestep and ε is the ratio of advection to collisional time scales, and (iii) is optimal (i.e., scalable with the total number of mesh points in the domain). We will demonstrate the algorithm in circular tokamak geometries [3].

[1] Chacón et al., arXiv:2105.01623

[2] Daniel et al., CPC (2020) 

[3] McDevitt et al., EPL (2019)