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. Following N. Crouseilles et al. [A parallel Vlasov solver based on local cubic spline interpolation on patches. Journal of Computational Physics (2009)], 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.

References in corpus (6)

• Stable interpolation with isotropic and anisotropic Gaussians using Hermite generating function
• A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation
• Jet schemes for advection problems
• Non-linear evolution of the diocotron instability in a pulsar electrosphere: 2D PIC simulations
• The diocotron instability in a pulsar cylindrical electrosphere
• Solving hyperbolic-elliptic problems on singular mapped disk-like domains with the method of characteristics and spline finite elements