A Direct and Accurate Adaptive Semi-Lagrangian Scheme for the Vlasov-Poisson Equation

A Direct and Accurate Adaptive Semi-Lagrangian Scheme for the Vlasov-Poisson Equation This article aims at giving a simplified presentation of a new adaptive semi-Lagrangian scheme for solving the (1 + 1)-dimensional Vlasov-Poisson system, which was developed in 2005 with Michel Mehrenberger and first described in (Campos Pinto and Mehrenberger, 2007). The main steps of the analysis are also given, which yield the first error estimate for an adaptive scheme in the context of the Vlasov equation. This article focuses on a key feature of our method, which is a new algorithm to transport multiscale meshes along a smooth flow, in a way that can be said optimal in the sense that it satisfies both accuracy and complexity estimates which are likely to lead to optimal convergence rates for the whole numerical scheme. From the regularity analysis of the numerical solution and how it gets transported by the numerical flow, it is shown that the accuracy of our scheme is monitored by a prescribed tolerance parameter ε which represents the local interpolation error at each time step. As a consequence, the numerical solutions are proved to converge in L∞ towards the exact ones as ε and δt tend to zero, and in addition to the numerical tests presented in (Campos Pinto and Mehrenberger, 2007), some complexity bounds are established which are likely to prove the optimality of the meshes.

[1]  Wolfgang Dahmen Adaptive approximation by multivariate smooth splines , 1982 .

[2]  R. DeVore,et al.  Nonlinear approximation , 1998, Acta Numerica.

[3]  Jeffery Cooper,et al.  Boundary value problems for the Vlasov-Maxwell equation in one dimension , 1980 .

[4]  Jean,et al.  Développement et analyse de méthodes adaptatives pour les équations de transport , 2005 .

[5]  Francis Filbet,et al.  Vlasov simulations of beams with a moving grid , 2004, Comput. Phys. Commun..

[6]  Eric Sonnendrücker,et al.  Vlasov simulations on an adaptive phase-space grid , 2004, Comput. Phys. Commun..

[7]  Martin Campos Pinto,et al.  Convergence of an adaptive semi-Lagrangian scheme for the Vlasov-Poisson system , 2007, Numerische Mathematik.

[8]  P. Raviart An analysis of particle methods , 1985 .

[9]  Albert Cohen,et al.  Fully adaptive multiresolution finite volume schemes for conservation laws , 2003, Math. Comput..

[10]  E. Sonnendrücker,et al.  The Semi-Lagrangian Method for the Numerical Resolution of the Vlasov Equation , 1999 .

[11]  Francis Filbet,et al.  An adaptive numerical method for the Vlasov equation based on a multiresolution analysis , 2007, 0704.1595.

[12]  Harry Yserentant,et al.  Hierarchical bases , 1992 .

[13]  R. Glassey,et al.  The Cauchy Problem in Kinetic Theory , 1987 .

[14]  Wolfgang Dahmen,et al.  Multiresolution schemes for conservation laws , 2001 .

[15]  Olivier Roussel,et al.  A conservative fully adaptive multiresolution algorithm for parabolic PDEs , 2003 .

[16]  Martin Campos Pinto,et al.  Convergence of an Adaptive Scheme for the one dimensional Vlasov-Poisson system , 2005 .

[17]  G. Knorr,et al.  The integration of the vlasov equation in configuration space , 1976 .

[18]  E. Sonnendrücker,et al.  Numerical Methods for Hyperbolic and Kinetic Problems , 2005 .

[19]  Nicolas Besse Convergence of a Semi-Lagrangian Scheme for the One-Dimensional Vlasov-Poisson System , 2004, SIAM J. Numer. Anal..