An Asymptotically Stable Semi-Lagrangian scheme in the Quasi-neutral Limit

This paper deals with the numerical simulations of the Vlasov-Poisson equation using a phase space grid in the quasi-neutral regime. In this limit, explicit numerical schemes suffer from numerical constraints related to the small Debye length and large plasma frequency. Here, we propose a semi-Lagrangian scheme for the Vlasov-Poisson model in the quasi-neutral limit. The main ingredient relies on a reformulation of the Poisson equation derived in (Crispel et al. in C. R. Acad. Sci. Paris, Ser. I 341:341–346, 2005) which enables asymptotically stable simulations. This scheme has a comparable numerical cost per time step to that of an explicit scheme. Moreover, it is not constrained by a restriction on the size of the time and length step when the Debye length and plasma period go to zero. A stability analysis and numerical simulations confirm this statement.

[1]  R. J. Mason,et al.  Implicit moment PIC-hybrid simulation of collisional plasmas , 1983 .

[2]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .

[3]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[4]  J. Delcroix.  Physique des plasmas , 1963 .

[5]  Bruce I. Cohen,et al.  Implicit time integration for plasma simulation , 1982 .

[6]  P. Rambo.  Finite-grid instability in quasineutral hybrid simulations , 1995 .

[7]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[8]  Jian-Guo Liu,et al.  Analysis of an Asymptotic Preserving Scheme for the Euler-Poisson System in the Quasineutral Limit , 2008, SIAM J. Numer. Anal..

[9]  T. Yabe,et al.  Cubic interpolated propagation scheme for solving the hyper-dimensional Vlasov-Poisson equation in phase space , 1999 .

[10]  Eric Sonnendrücker,et al.  The Semi-Lagrangian Method for the Numerical Resolution of Vlasov Equations , 1998 .

[11]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[12]  D. Hewett,et al.  Electromagnetic direct implicit plasma simulation , 1987 .

[13]  Francis Filbet,et al.  Numerical approximation of collisional plasmas by high order methods , 2004 .

[14]  C. W. Nielson,et al.  A multidimensional quasineutral plasma simulation model , 1978 .

[15]  P. Bertrand,et al.  Conservative numerical schemes for the Vlasov equation , 2001 .

[16]  Céline Parzani,et al.  Plasma Expansion in Vacuum: Modeling the Breakdown of Quasi Neutrality , 2003, Multiscale Model. Simul..

[17]  M. Shoucri.  Nonlinear evolution of the bump‐on‐tail instability , 1979 .

[18]  Guillaume Latu,et al.  Hermite Spline Interpolation on Patches for Parallelly Solving the Vlasov-Poisson Equation , 2007, Int. J. Appl. Math. Comput. Sci..

[19]  Y. Brenier,et al.  convergence of the vlasov-poisson system to the incompressible euler equations , 2000 .

[20]  Pierre Degond,et al.  An asymptotically stable discretization for the Euler–Poisson system in the quasi-neutral limit , 2005 .

[21]  R. J. Mason,et al.  Implicit moment particle simulation of plasmas , 1981 .

[22]  Fabrice Deluzet,et al.  An asymptotically stable Particle-in-Cell (PIC) scheme for collisionless plasma simulations near quasineutrality , 2006 .

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

[24]  E. Sonnendrücker,et al.  Comparison of Eulerian Vlasov solvers , 2003 .

[25]  Burton D. Fried,et al.  The Plasma Dispersion Function , 1961 .

[26]  Pierre Degond,et al.  An asymptotic preserving scheme for the two-fluid Euler-Poisson model in the quasineutral limit , 2007, J. Comput. Phys..

[27]  J. Denavit,et al.  Hybrid simulation of ion beams in background plasma , 1987 .

[28]  Bruce I. Cohen,et al.  Direct implicit large time-step particle simulation of plasmas , 1983 .

[29]  Seung-Yeal Ha,et al.  Global Existence of Plasma Ion-Sheaths and Their Dynamics , 2003 .

[30]  Sylvie Fabre.  Stability analysis of the Euler-poisson equations , 1992 .