Local cubic spline interpolation for Vlasov-type equations on a multi-patch geometry

Abstract: We present a semi-Lagrangian method for the numerical resolution of Vlasov-type equations on multi-patch meshes. We employ a local cubic spline interpolation with Hermite boundary conditions between the patches. The derivative reconstruction is adapted to cope with non-uniform meshes as well as non-conforming situations. In the conforming case, there are no longer any constraints on the number of points for each patch; however, a small global system must now be solved. In that case, the local spline representations coincide with the corresponding global spline reconstruction. Alternatively, we can choose not to apply the global system and the derivatives can be approximated. The influence of the most distant points diminishes as the number of points per patch increases. For uniform per patch configurations, a study of the explicit and asymptotic behavior of this influence has been led. The method is validated using a two-dimensional guiding-center model with an O-point. All the numerical results are carried out in the Gyselalib++ library.