Using linear multistep methods for the time stepping in Vlasov–Poisson simulations

2021 ◽  
Vol 40 (8) ◽  
Author(s):  
Denis Lorenzon ◽  
Sergio Elaskar
Author(s):  
Giacomo Albi ◽  
Lorenzo Pareschi

AbstractWe consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form, typically used in implicit-explicit (IMEX) methods, is not possible. As shown in Boscarino et al. (J. Sci. Comput. 68: 975–1001, 2016) for Runge-Kutta methods, these semi-implicit techniques give a great flexibility, and allow, in many cases, the construction of simple linearly implicit schemes with no need of iterative solvers. In this work, we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype linear advection-diffusion equation and in the setting of strong stability preserving (SSP) methods. Our findings are demonstrated on several examples, including nonlinear reaction-diffusion and convection-diffusion problems.


Mathematics ◽  
2019 ◽  
Vol 7 (12) ◽  
pp. 1197
Author(s):  
Changbum Chun ◽  
Beny Neta

Numerical methods for the solution of ordinary differential equations are based on polynomial interpolation. In 1952, Brock and Murray have suggested exponentials for the case that the solution is known to be of exponential type. In 1961, Gautschi came up with the idea of using information on the frequency of a solution to modify linear multistep methods by allowing the coefficients to depend on the frequency. Thus the methods integrate exactly appropriate trigonometric polynomials. This was done for both first order systems and second order initial value problems. Gautschi concluded that “the error reduction is not very substantial unless” the frequency estimate is close enough. As a result, no other work was done in this direction until 1984 when Neta and Ford showed that “Nyström’s and Milne-Simpson’s type methods for systems of first order initial value problems are not sensitive to changes in frequency”. This opened the flood gates and since then there have been many papers on the subject.


Sign in / Sign up

Export Citation Format

Share Document