With cedram.org   version française
Search for an article
Table of contents for this volume | Previous article | Next article
Martin Campos Pinto; Frédérique Charles
From particle methods to forward-backward Lagrangian schemes
SMAI-Journal of computational mathematics, 4 (2018), p. 121-150, doi: 10.5802/smai-jcm.31
Article PDF


In this article we study a novel method for improving the accuracy of density reconstructions based on markers pushed forward by some available particle code. The method relies on the backward Lagrangian representation of the transported density, and it evaluates the backward flow using the current position of point particles seen as flow markers. Compared to existing smooth particle methods with either fixed or transformed shapes, the proposed reconstruction achieves higher locality and accuracy. This is confirmed by our error analysis which shows a theoretical gain of one convergence order compared to the LTP/QTP methods introduced in [8], and by numerical experiments that demonstrate significant CPU gains and an improved robustness relative to the remapping period.


[1] C. Alard & S. Colombi, “A cloudy Vlasov solution”, Monthly Notices of the Royal Astronomical Society 359 (2005) no. 1, p. 123-163
[2] W.B. Bateson & D.W. Hewett, “Grid and Particle Hydrodynamics”, Journal of Computational Physics 144 (1998), p. 358-378
[3] J.T. Beale & A. Majda, “Vortex methods. II. Higher order accuracy in two and three dimensions”, Mathematics of Computation 39 (1982) no. 159, p. 29-52
[4] M. Bergdorf, G.-H. Cottet & P. Koumoutsakos, “Multilevel Adaptive Particle Methods for Convection-Diffusion Equations”, Multiscale Modeling & Simulation 4 (2005) no. 1, p. 328-357
[5] M. Bergdorf & P. Koumoutsakos, “A Lagrangian particle-wavelet method”, Multiscale Modeling & Simulation 5 (2006) no. 3, p. 980-995
[6] A. Biancalani, A. Bottino, S. Briguglio, A. Koenies, Ph. Lauber, A. Mishchenko, E. Poli, B.D. Scott & F. Zonca, “Linear gyrokinetic particle-in-cell simulations of Alfven instabilities in tokamaks”, arXiv:1510.01945, 2015
[7] O. Bokanowski, J. Garcke, M. Griebel & I. Klompmaker, “An adaptive sparse grid semi-Lagrangian scheme for first order Hamilton-Jacobi Bellman equations”, Journal of Scientific Computing 55 (2013) no. 3, p. 575-605
[8] M. Campos Pinto, “Towards smooth particle methods without smoothing”, Journal of Scientific Computing (2014)
[9] M. Campos Pinto & F. Charles, “Uniform Convergence of a Linearly Transformed Particle Method for the Vlasov–Poisson System”, SIAM Journal on Numerical Analysis 54 (2016) no. 1, p. 137-160
[10] M. Campos Pinto, E. Sonnendrücker, A. Friedman, D.P. Grote & S.M. Lund, “Noiseless Vlasov–Poisson simulations with linearly transformed particles”, Journal of Computational Physics 275 (2014) no. C, p. 236-256
[11] Y. Chen, S.E. Parker, G. Rewoldt, S.-H. Ku, G.-Y. Park & C.-S. Chang, “Coarse-graining the electron distribution in turbulence simulations of tokamak plasmas”, Physics of Plasmas 15 (2008) no. 5
[12] A. Cohen & B. Perthame, “Optimal Approximations of Transport Equations by Particle and Pseudoparticle Methods”, SIAM Journal on Mathematical Analysis 32 (2000) no. 3, p. 616-636
[13] S. Colombi & C. Alard, “A “metric” semi-Lagrangian Vlasov-Poisson solver”, Submitted, 2016
[14] C.J. Cotter, J. Frank & S. Reich, “The remapped particle-mesh semi-Lagrangian advection scheme”, Quarterly Journal of the Royal Meteorological Society 133 (2007) no. 622, p. 251-260
[15] G.-H. Cottet, P. Koumoutsakos & M.L.O. Salihi, “Vortex Methods with Spatially Varying Cores”, Journal of Computational Physics 162 (2000) no. 1, p. 164-185
[16] G.H. Cottet & P. Koumoutsakos, Vortex Methods: Theory and Practice, Cambridge University Press, Cambridge, 2000
[17] N. Crouseilles, T. Respaud & E. Sonnendrücker, “A forward semi-Lagrangian method for the numerical solution of the Vlasov equation”, Computer Physics Communications 180 (2009) no. 10, p. 1730-1745
[18] J. Denavit, “Numerical Simulation of Plasmas with Periodic Smoothing in Phase Space”, Journal of Computational Physics 9 (1972), p. 75-98
[19] F. Filbet & E. Sonnendrücker, “Comparison of Eulerian Vlasov solvers”, Computer Physics Communications 150 (2003), p. 247-266
[20] Ch. Gauger, P. Leinen & H. Yserentant, “The finite mass method”, SIAM J. Numer. Anal. 37 (2000) no. 6, p. 1768-1799
[21] O.H. Hald, “Convergence of Vortex Methods for Euler’s Equations. II”, SIAM Journal on Numerical Analysis 16 (1979) no. 5, p. 726-755
[22] R.W. Hockney & J.W. Eastwood, Computer simulation using particles, Taylor & Francis, Inc, Bristol, PA, USA, 1988
[23] T.Y. Hou, “Convergence of a Variable Blob Vortex Method for the Euler and Navier-Stokes Equations”, SIAM Journal on Numerical Analysis 27 (1990) no. 6, p. 1387-1404
[24] P. Koumoutsakos, “Inviscid Axisymmetrization of an Elliptical Vortex”, Journal of Computational Physics 138 (1997), p. 821-857
[25] P. Koumoutsakos, “Inviscid axisymmetrization of an elliptical vortex”, Journal of Computational Physics 138 (1997) no. 2, p. 821-857
[26] R.J. LeVeque, “High-resolution conservative algorithms for advection in incompressible flow”, SIAM Journal on Numerical Analysis (1996), p. 627-665
[27] A. Magni & G.-H. Cottet, “Accurate, non-oscillatory, remeshing schemes for particle methods”, Journal of Computational Physics 231 (2012) no. 1, p. 152-172
[28] J.J. Monaghan, “Extrapolating B. Splines for Interpolation”, Journal of Computational Physics 60 (1985)
[29] R.D. Nair, J.S. Scroggs & F.H.M. Semazzi, “A forward-trajectory global semi-Lagrangian transport scheme”, Journal of Computational Physics 190 (2003) no. 1, p. 275-294
[30] W.M. Nevins, G.W. Hammett, A.M. Dimits, W. Dorland & D.E. Shumaker, “Discrete particle noise in particle-in-cell simulations of plasma microturbulence”, Physics of Plasmas 12 (2005)
[31] J.-M. Qiu & A. Christlieb, “A conservative high order semi-Lagrangian WENO method for the Vlasov equation”, Journal of Computational Physics 229 (2010), p. 1130-1149
[32] P.-A. Raviart, An analysis of particle methods, Numerical methods in fluid dynamics (Como, 1983), Lecture Notes in Mathematics, 1985, p. 243–324
[33] Selalib, “Semi-Lagrangian Library”, http://selalib.gforge.inria.fr/
[34] E. Sonnendrücker, J. Roche, P. Bertrand & A. Ghizzo, “The semi-Lagrangian method for the numerical resolution of the Vlasov equation”, Journal of Computational Physics 149 (1999) no. 2, p. 201-220
[35] A. Staniforth & J. Côté, “Semi-Lagrangian integration schemes for atmospheric models-a review”, Monthly weather review 119 (1991), p. 2206-2223
[36] M. Unser & I. Daubechies, “On the approximation power of convolution-based least squares versus interpolation”, Signal Processing, IEEE Transactions on 45 (1997) no. 7, p. 1697-1711
[37] B. Wang, G.H. Miller & P. Colella, “A Particle-In-Cell method with adaptive phase-space remapping for kinetic plasmas”, SIAM Journal on Scientific Computing 33 (2011), p. 3509-3537