DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

Abstract

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

Authors:
; ORCiD logo
Publication Date:
Research Org.:
Univ. of Nebraska, Lincoln, NE (United States)
Sponsoring Org.:
USDOE; National Science Foundation (NSF)
OSTI Identifier:
1893545
Alternate Identifier(s):
OSTI ID: 1894688
Grant/Contract Number:  
FG02-08ER55000; SC0008382
Resource Type:
Published Article
Journal Name:
Journal of Plasma Physics
Additional Journal Information:
Journal Name: Journal of Plasma Physics Journal Volume: 88 Journal Issue: 2; Journal ID: ISSN 0022-3778
Publisher:
Cambridge University Press (CUP)
Country of Publication:
United Kingdom
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; plasma simulation; plasma instabilities; plasma dynamics

Citation Formats

Carrié, M., and Shadwick, B. A. An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system. United Kingdom: N. p., 2022. Web. doi:10.1017/S0022377821001124.
Carrié, M., & Shadwick, B. A. An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system. United Kingdom. https://doi.org/10.1017/S0022377821001124
Carrié, M., and Shadwick, B. A. Tue . "An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system". United Kingdom. https://doi.org/10.1017/S0022377821001124.
@article{osti_1893545,
title = {An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system},
author = {Carrié, M. and Shadwick, B. A.},
abstractNote = {The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.},
doi = {10.1017/S0022377821001124},
journal = {Journal of Plasma Physics},
number = 2,
volume = 88,
place = {United Kingdom},
year = {Tue Mar 08 00:00:00 EST 2022},
month = {Tue Mar 08 00:00:00 EST 2022}
}

Works referenced in this record:

Principles of Plasma Physics
journal, December 1973

  • Krall, Nicholas A.; Trivelpiece, Alvin W.; Gross, Robert A.
  • American Journal of Physics, Vol. 41, Issue 12
  • DOI: 10.1119/1.1987587

Conservative Numerical Schemes for the Vlasov Equation
journal, September 2001

  • Filbet, Francis; Sonnendrücker, Eric; Bertrand, Pierre
  • Journal of Computational Physics, Vol. 172, Issue 1
  • DOI: 10.1006/jcph.2001.6818

On the theory of stationary waves in plasmas
journal, January 1955


On the velocity space discretization for the Vlasov–Poisson system: Comparison between implicit Hermite spectral and Particle-in-Cell methods
journal, January 2016


Vlasov Simulations Using Velocity-Scaled Hermite Representations
journal, August 1998

  • Schumer, Joseph W.; Holloway, James Paul
  • Journal of Computational Physics, Vol. 144, Issue 2
  • DOI: 10.1006/jcph.1998.5925

Iterative Refinement in Floating Point
journal, April 1967


Thermal effects in plasma-based accelerators
journal, May 2007

  • Esarey, E.; Schroeder, C. B.; Cormier-Michel, E.
  • Physics of Plasmas, Vol. 14, Issue 5
  • DOI: 10.1063/1.2714022

A time-implicit numerical method and benchmarks for the relativistic Vlasov–Ampere equations
journal, January 2016

  • Carrié, Michael; Shadwick, B. A.
  • Physics of Plasmas, Vol. 23, Issue 1
  • DOI: 10.1063/1.4938035

Comparison of two Eulerian solvers for the four-dimensional Vlasov equation: Part I
journal, February 2008

  • Crouseilles, N.; Gutnic, M.; Latu, G.
  • Communications in Nonlinear Science and Numerical Simulation, Vol. 13, Issue 1
  • DOI: 10.1016/j.cnsns.2007.03.010

A finite element code for the simulation of one-dimensional vlasov plasmas. II. Applications
journal, November 1988


Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004


On the Stability of Plasma in Static Equilibrium
journal, January 1958

  • Kruskal, M. D.; Oberman, C. R.
  • Physics of Fluids, Vol. 1, Issue 4
  • DOI: 10.1063/1.1705885

Quadratic conservative scheme for relativistic Vlasov–Maxwell system
journal, February 2019

  • Shiroto, Takashi; Ohnishi, Naofumi; Sentoku, Yasuhiko
  • Journal of Computational Physics, Vol. 379
  • DOI: 10.1016/j.jcp.2018.10.041

Composition constants for raising the orders of unconventional schemes for ordinary differential equations
journal, July 1997


Variational formulation of macro-particle plasma simulation algorithms
journal, May 2014

  • Shadwick, B. A.; Stamm, A. B.; Evstatiev, E. G.
  • Physics of Plasmas, Vol. 21, Issue 5
  • DOI: 10.1063/1.4874338

Geometric Numerical Integration
book, January 2002

  • Hairer, Ernst; Wanner, Gerhard; Lubich, Christian
  • Springer Series in Computational Mathematics
  • DOI: 10.1007/978-3-662-05018-7

Radiotherapy using a laser proton accelerator
conference, January 2008

  • Murakami, Masao; Hishikawa, Yoshio; Miyajima, Satoshi
  • LASER-DRIVEN RELATIVISTIC PLASMAS APPLIED FOR SCIENCE, INDUSTRY, AND MEDICINE: The 1st International Symposium, AIP Conference Proceedings
  • DOI: 10.1063/1.2958203

Updating the Inverse of a Matrix
journal, June 1989

  • Hager, William W.
  • SIAM Review, Vol. 31, Issue 2, p. 221-239
  • DOI: 10.1137/1031049

Gamma-rays from harmonically resonant betatron oscillations in a plasma wake
journal, September 2011

  • Cipiccia, Silvia; Islam, Mohammad R.; Ersfeld, Bernhard
  • Nature Physics, Vol. 7, Issue 11
  • DOI: 10.1038/nphys2090

A discontinuous Galerkin method for the Vlasov–Poisson system
journal, February 2012

  • Heath, R. E.; Gamba, I. M.; Morrison, P. J.
  • Journal of Computational Physics, Vol. 231, Issue 4
  • DOI: 10.1016/j.jcp.2011.09.020

Eulerian codes for the numerical solution of the Vlasov equation
journal, February 2008


Exactly Conservative Integrators
journal, January 1998

  • Shadwick, B. A.; Bowman, John C.; Morrison, P. J.
  • SIAM Journal on Applied Mathematics, Vol. 59, Issue 3
  • DOI: 10.1137/S0036139995289313

VALIS: A split-conservative scheme for the relativistic 2D Vlasov–Maxwell system
journal, July 2009


A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type
journal, January 1947

  • Crank, J.; Nicolson, P.
  • Mathematical Proceedings of the Cambridge Philosophical Society, Vol. 43, Issue 1
  • DOI: 10.1017/S0305004100023197

Numerical solution of the Vlasov–Poisson system using generalized Hermite functions
journal, October 2006

  • Le Bourdiec, S.; de Vuyst, F.; Jacquet, L.
  • Computer Physics Communications, Vol. 175, Issue 8
  • DOI: 10.1016/j.cpc.2006.07.004

The free energy of Maxwell–Vlasov equilibria
journal, June 1990

  • Morrison, P. J.; Pfirsch, D.
  • Physics of Fluids B: Plasma Physics, Vol. 2, Issue 6
  • DOI: 10.1063/1.859246

Physics of laser-driven plasma-based electron accelerators
journal, August 2009


On the Construction and Comparison of Difference Schemes
journal, September 1968

  • Strang, Gilbert
  • SIAM Journal on Numerical Analysis, Vol. 5, Issue 3
  • DOI: 10.1137/0705041

The Semi-Lagrangian Method for the Numerical Resolution of the Vlasov Equation
journal, March 1999

  • Sonnendrücker, Eric; Roche, Jean; Bertrand, Pierre
  • Journal of Computational Physics, Vol. 149, Issue 2
  • DOI: 10.1006/jcph.1998.6148

Construction of higher order symplectic integrators
journal, November 1990


Comparison of Eulerian Vlasov solvers
journal, February 2003


The integration of the vlasov equation in configuration space
journal, November 1976


Electrostatic Instabilities of a Uniform Non-Maxwellian Plasma
journal, January 1960


A Numerical Scheme for the Integration of the Vlasov–Maxwell System of Equations
journal, July 2002

  • Mangeney, A.; Califano, F.; Cavazzoni, C.
  • Journal of Computational Physics, Vol. 179, Issue 2
  • DOI: 10.1006/jcph.2002.7071

Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime
journal, June 2012


Plasma Physics via Computer Simulation
book, January 1991


Comparison of two Eulerian solvers for the four-dimensional Vlasov equation: Part II
journal, February 2008

  • Crouseilles, N.; Gutnic, M.; Latu, G.
  • Communications in Nonlinear Science and Numerical Simulation, Vol. 13, Issue 1
  • DOI: 10.1016/j.cnsns.2007.03.017

On the Numerical Integration of Ordinary Differential Equations by Symmetric Composition Methods
journal, January 1995

  • McLachlan, Robert I.
  • SIAM Journal on Scientific Computing, Vol. 16, Issue 1
  • DOI: 10.1137/0916010

Bound on the Energy Available from a Plasma
journal, January 1963


Unphysical kinetic effects in particle-in-cell modeling of laser wakefield accelerators
journal, July 2008


A Critical Comparison of Eulerian-Grid-Based Vlasov Solvers
journal, July 2002

  • Arber, T. D.; Vann, R. G. L.
  • Journal of Computational Physics, Vol. 180, Issue 1
  • DOI: 10.1006/jcph.2002.7098

A numerical solution to the Vlasov equation
journal, February 1999