Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration
Abstract
Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this study, we propose an implicit–explicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithm than those using explicit time integrators, where the time step sizes are constrained by the stiff collisional time scales. Finally, we implement semi-implicit additive Runge–Kutta methods in COGENT, a high-order finite-volume gyrokinetic code and test the accuracy, convergence, and computational cost of these semi-implicit methods for test cases with highly-collisional plasmas.
- Authors:
-
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientific Computing
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Physics Division
- Publication Date:
- Research Org.:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
- OSTI Identifier:
- 1458696
- Report Number(s):
- LLNL-JRNL-735522
Journal ID: ISSN 0885-7474; 887762
- Grant/Contract Number:
- AC52-07NA27344
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Scientific Computing
- Additional Journal Information:
- Journal Volume: 77; Journal ID: ISSN 0885-7474
- Publisher:
- Springer
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; IMEX time integration; plasma physics; gyrokinetic simulations; Vlasov–Fokker–Planck equations
Citation Formats
Ghosh, Debojyoti, Dorf, Mikhail A., Dorr, Milo R., and Hittinger, Jeffrey A. F. Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration. United States: N. p., 2018.
Web. doi:10.1007/s10915-018-0726-6.
Ghosh, Debojyoti, Dorf, Mikhail A., Dorr, Milo R., & Hittinger, Jeffrey A. F. Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration. United States. https://doi.org/10.1007/s10915-018-0726-6
Ghosh, Debojyoti, Dorf, Mikhail A., Dorr, Milo R., and Hittinger, Jeffrey A. F. Sat .
"Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration". United States. https://doi.org/10.1007/s10915-018-0726-6. https://www.osti.gov/servlets/purl/1458696.
@article{osti_1458696,
title = {Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration},
author = {Ghosh, Debojyoti and Dorf, Mikhail A. and Dorr, Milo R. and Hittinger, Jeffrey A. F.},
abstractNote = {Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this study, we propose an implicit–explicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithm than those using explicit time integrators, where the time step sizes are constrained by the stiff collisional time scales. Finally, we implement semi-implicit additive Runge–Kutta methods in COGENT, a high-order finite-volume gyrokinetic code and test the accuracy, convergence, and computational cost of these semi-implicit methods for test cases with highly-collisional plasmas.},
doi = {10.1007/s10915-018-0726-6},
journal = {Journal of Scientific Computing},
number = ,
volume = 77,
place = {United States},
year = {2018},
month = {5}
}
Web of Science
Works referenced in this record:
Progress in Kinetic Simulation of Edge Plasmas
journal, March 2008
- Cohen, R. H.; Xu, X. Q.
- Contributions to Plasma Physics, Vol. 48, Issue 1-3
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
Tempest Simulations of Collisionless Damping of the Geodesic-Acoustic Mode in Edge-Plasma Pedestals
journal, May 2008
- Xu, X. Q.; Xiong, Z.; Gao, Z.
- Physical Review Letters, Vol. 100, Issue 21
Detailed comparison of simulated and measured plasma profiles in the scrape-off layer and edge plasma of DIII-D
journal, September 2000
- Porter, G. D.; Isler, R.; Boedo, J.
- Physics of Plasmas, Vol. 7, Issue 9
Discretization methods for one-dimensional Fokker-Planck operators
journal, December 1985
- Larsen, E. W.; Levermore, C. D.; Pomraning, G. C.
- Journal of Computational Physics, Vol. 61, Issue 3
Entropy Stable Numerical Schemes for Two-Fluid Plasma Equations
journal, November 2011
- Kumar, Harish; Mishra, Siddhartha
- Journal of Scientific Computing, Vol. 52, Issue 2
An IMEX-Scheme for Pricing Options under Stochastic Volatility Models with Jumps
journal, January 2014
- Salmi, Santtu; Toivanen, Jari; von Sydow, Lina
- SIAM Journal on Scientific Computing, Vol. 36, Issue 5
Gyrokinetic study of the edge shear layer
journal, April 2006
- Scott, B.
- Plasma Physics and Controlled Fusion, Vol. 48, Issue 5A
Conservative gyrokinetic Vlasov simulation
journal, February 2008
- Idomura, Yasuhiro; Ida, Masato; Tokuda, Shinji
- Communications in Nonlinear Science and Numerical Simulation, Vol. 13, Issue 1
A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov–Poisson equations
journal, July 2011
- Rossmanith, James A.; Seal, David C.
- Journal of Computational Physics, Vol. 230, Issue 16
Progress with the COGENT Edge Kinetic Code: Collision Operator Options
journal, June 2012
- Dorf, M. A.; Cohen, R. H.; Compton, J. C.
- Contributions to Plasma Physics, Vol. 52, Issue 5-6
Relaxation of a system of charged particles
journal, July 1985
- Kho, T. H.
- Physical Review A, Vol. 32, Issue 1
TEMPEST simulations of the plasma transport in a single-null tokamak geometry
journal, May 2010
- Xu, X. Q.; Bodi, K.; Cohen, R. H.
- Nuclear Fusion, Vol. 50, Issue 6
A drift-kinetic Semi-Lagrangian 4D code for ion turbulence simulation
journal, September 2006
- Grandgirard, V.; Brunetti, M.; Bertrand, P.
- Journal of Computational Physics, Vol. 217, Issue 2
Fast implicit schemes for the Fokker–Planck–Landau equation
journal, May 2004
- Lemou, Mohammed; Mieussens, Luc
- Comptes Rendus Mathematique, Vol. 338, Issue 10
A mass, momentum, and energy conserving, fully implicit, scalable algorithm for the multi-dimensional, multi-species Rosenbluth–Fokker–Planck equation
journal, September 2015
- Taitano, W. T.; Chacón, L.; Simakov, A. N.
- Journal of Computational Physics, Vol. 297
Nonlinear gyrokinetic equations for turbulence in core transport barriers
journal, December 1996
- Hahm, T. S.
- Physics of Plasmas, Vol. 3, Issue 12
Additive Runge–Kutta schemes for convection–diffusion–reaction equations
journal, January 2003
- Kennedy, Christopher A.; Carpenter, Mark H.
- Applied Numerical Mathematics, Vol. 44, Issue 1-2
Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004
- Knoll, D. A.; Keyes, D. E.
- Journal of Computational Physics, Vol. 193, Issue 2
Edge gyrokinetic theory and continuum simulations
journal, July 2007
- Xu, X. Q.; Xiong, Z.; Dorr, M. R.
- Nuclear Fusion, Vol. 47, Issue 8
NITSOL: A Newton Iterative Solver for Nonlinear Systems
journal, January 1998
- Pernice, Michael; Walker, Homer F.
- SIAM Journal on Scientific Computing, Vol. 19, Issue 1
Implicit–explicit time integration of a high-order particle-in-cell method with hyperbolic divergence cleaning
journal, October 2009
- Jacobs, G. B.; Hesthaven, J. S.
- Computer Physics Communications, Vol. 180, Issue 10
High-order finite-volume methods for hyperbolic conservation laws on mapped multiblock grids
journal, May 2015
- McCorquodale, P.; Dorr, M. R.; Hittinger, J. A. F.
- Journal of Computational Physics, Vol. 288
Fast Algorithms for Numerical, Conservative, and Entropy Approximations of the Fokker–Planck–Landau Equation
journal, May 1997
- Buet, C.; Cordier, S.; Degond, P.
- Journal of Computational Physics, Vol. 133, Issue 2
A forward semi-Lagrangian method for the numerical solution of the Vlasov equation
journal, October 2009
- Crouseilles, Nicolas; Respaud, Thomas; Sonnendrücker, Eric
- Computer Physics Communications, Vol. 180, Issue 10
Fokker-Planck Equation for an Inverse-Square Force
journal, July 1957
- Rosenbluth, Marshall N.; MacDonald, William M.; Judd, David L.
- Physical Review, Vol. 107, Issue 1
New velocity-space discretization for continuum kinetic calculations and Fokker–Planck collisions
journal, June 2013
- Landreman, Matt; Ernst, Darin R.
- Journal of Computational Physics, Vol. 243
An implicit Vlasov–Fokker–Planck code to model non-local electron transport in 2-D with magnetic fields
journal, February 2004
- Kingham, R. J.; Bell, A. R.
- Journal of Computational Physics, Vol. 194, Issue 1
Implicit-explicit runge-kutta schemes and applications to hyperbolic systems with relaxation
journal, November 2005
- Pareschi, Lorenzo; Russo, Giovanni
- Journal of Scientific Computing, Vol. 25, Issue 1-2
Kinetic simulation of a plasma collision experiment
journal, August 1993
- Larroche, Olivier
- Physics of Fluids B: Plasma Physics, Vol. 5, Issue 8
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986
- Saad, Youcef; Schultz, Martin H.
- SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
Positivity preserving semi-Lagrangian discontinuous Galerkin formulation: Theoretical analysis and application to the Vlasov–Poisson system
journal, September 2011
- Qiu, Jing-Mei; Shu, Chi-Wang
- Journal of Computational Physics, Vol. 230, Issue 23
Implicit–explicit Runge–Kutta methods applied to atmospheric chemistry-transport modelling
journal, September 2000
- Wolke, Ralf; Knoth, Oswald
- Environmental Modelling & Software, Vol. 15, Issue 6-7
An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov–Darwin particle-in-cell algorithm
journal, October 2014
- Chen, G.; Chacón, L.
- Computer Physics Communications, Vol. 185, Issue 10
Particle Simulation of the Neoclassical Plasmas
journal, November 2001
- Heikkinen, J. A.; Kiviniemi, T. P.; Kurki-Suonio, T.
- Journal of Computational Physics, Vol. 173, Issue 2
Conservative finite-difference schemes for the Fokker-Planck equation not violating the law of an increasing entropy
journal, March 1987
- Berezin, Yu. A.; Khudick, V. N.; Pekker, M. S.
- Journal of Computational Physics, Vol. 69, Issue 1
Simulation of neoclassical transport with the continuum gyrokinetic code COGENT
journal, January 2013
- Dorf, M. A.; Cohen, R. H.; Dorr, M.
- Physics of Plasmas, Vol. 20, Issue 1
Fully Implicit Kinetic Solution of Collisional Plasmas
journal, September 1997
- Mousseau, V. A.; Knoll, D. A.
- Journal of Computational Physics, Vol. 136, Issue 2
Spontaneous rotation sources in a quiescent tokamak edge plasma
journal, June 2008
- Chang, C. S.; Ku, S.
- Physics of Plasmas, Vol. 15, Issue 6
Full linearized Fokker–Planck collisions in neoclassical transport simulations
journal, December 2011
- Belli, E. A.; Candy, J.
- Plasma Physics and Controlled Fusion, Vol. 54, Issue 1
Progress with the COGENT Edge Kinetic Code: Implementing the Fokker-Planck Collision Operator: Progress with the COGENT Edge Kinetic Code: Implementing the Fokker-Planck Collision Operator
journal, June 2014
- Dorf, M. A.; Cohen, R. H.; Dorr, M.
- Contributions to Plasma Physics, Vol. 54, Issue 4-6
Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator
journal, September 2011
- Pataki, Andras; Greengard, Leslie
- Journal of Computational Physics, Vol. 230, Issue 21
Kinetic simulation of a collisional shock wave in a plasma
journal, October 1991
- Casanova, Michel; Larroche, Olivier; Matte, Jean-Pierre
- Physical Review Letters, Vol. 67, Issue 16
Conservative and Entropy Decaying Numerical Scheme for the Isotropic Fokker–Planck–Landau Equation
journal, September 1998
- Buet, C.; Cordier, S.
- Journal of Computational Physics, Vol. 145, Issue 1
Semi-Implicit Formulations of the Navier–Stokes Equations: Application to Nonhydrostatic Atmospheric Modeling
journal, January 2010
- Giraldo, F. X.; Restelli, M.; Läuter, M.
- SIAM Journal on Scientific Computing, Vol. 32, Issue 6
An adaptive, conservative 0D-2V multispecies Rosenbluth–Fokker–Planck solver for arbitrarily disparate mass and temperature regimes
journal, August 2016
- Taitano, W. T.; Chacón, L.; Simakov, A. N.
- Journal of Computational Physics, Vol. 318
An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm
journal, August 2011
- Chen, G.; Chacón, L.; Barnes, D. C.
- Journal of Computational Physics, Vol. 230, Issue 18
Kinetic simulations of fuel ion transport in ICF target implosions
journal, November 2003
- Larroche, O.
- The European Physical Journal D - Atomic, Molecular and Optical Physics, Vol. 27, Issue 2
FPPAC: A two-dimensional multispecies nonlinear Fokker-Planck package
journal, September 1981
- McCoy, M. G.; Mirin, A. A.; Killeen, J.
- Computer Physics Communications, Vol. 24, Issue 1
A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm
journal, December 2015
- Chen, G.; Chacón, L.
- Computer Physics Communications, Vol. 197
The solution of poisson's equation for isolated source distributions
journal, October 1977
- James, R. A.
- Journal of Computational Physics, Vol. 25, Issue 2
hypre: A Library of High Performance Preconditioners
book, April 2002
- Falgout, Robert D.; Yang, Ulrike Meier; Goos, Gerhard
- Computational Science — ICCS 2002: International Conference Amsterdam, The Netherlands, April 21–24, 2002 Proceedings, Part III
Implicit and Conservative Difference Scheme for the Fokker-Planck Equation
journal, June 1994
- Epperlein, E. M.
- Journal of Computational Physics, Vol. 112, Issue 2
Gyrokinetic Simulation of Particle and Heat Transport in the Presence of Wide Orbits and Strong Profile Variations in the Edge Plasma
journal, September 2006
- Heikkinen, J. A.; Henriksson, S.; Janhunen, S.
- Contributions to Plasma Physics, Vol. 46, Issue 7-9
Gyrokinetic approach in particle simulation
journal, January 1983
- Lee, W. W.
- Physics of Fluids, Vol. 26, Issue 2
A hybrid method of semi-Lagrangian and additive semi-implicit Runge–Kutta schemes for gyrokinetic Vlasov simulations
journal, September 2012
- Maeyama, S.; Ishizawa, A.; Watanabe, T. -H.
- Computer Physics Communications, Vol. 183, Issue 9
Implicit–Explicit Multistep Methods for Fast-Wave–Slow-Wave Problems
journal, April 2012
- Durran, Dale R.; Blossey, Peter N.
- Monthly Weather Review, Vol. 140, Issue 4
Conservative global gyrokinetic toroidal full-f five-dimensional Vlasov simulation
journal, September 2008
- Idomura, Yasuhiro; Ida, Masato; Kano, Takuma
- Computer Physics Communications, Vol. 179, Issue 6
An Implicit Energy-Conservative 2D Fokker–Planck Algorithm
journal, January 2000
- Chacón, L.; Barnes, D. C.; Knoll, D. A.
- Journal of Computational Physics, Vol. 157, Issue 2
Implicit-Explicit Formulations of a Three-Dimensional Nonhydrostatic Unified Model of the Atmosphere (NUMA)
journal, January 2013
- Giraldo, F. X.; Kelly, J. F.; Constantinescu, E. M.
- SIAM Journal on Scientific Computing, Vol. 35, Issue 5
Two-Dimensional Nonlocal Electron Transport in Laser-Produced Plasmas
journal, November 1988
- Epperlein, E. M.; Rickard, G. J.; Bell, A. R.
- Physical Review Letters, Vol. 61, Issue 21
Vlasov simulations of electron-ion collision effects on damping of electron plasma waves
journal, March 2016
- Banks, J. W.; Brunner, S.; Berger, R. L.
- Physics of Plasmas, Vol. 23, Issue 3
A second order self-consistent IMEX method for radiation hydrodynamics
journal, November 2010
- Kadioglu, Samet Y.; Knoll, Dana A.; Lowrie, Robert B.
- Journal of Computational Physics, Vol. 229, Issue 22
Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations
journal, November 1997
- Ascher, Uri M.; Ruuth, Steven J.; Spiteri, Raymond J.
- Applied Numerical Mathematics, Vol. 25, Issue 2-3
A review of Vlasov–Fokker–Planck numerical modeling of inertial confinement fusion plasma
journal, February 2012
- Thomas, A. G. R.; Tzoufras, M.; Robinson, A. P. L.
- Journal of Computational Physics, Vol. 231, Issue 3
Implicit Schemes for the Fokker--Planck--Landau Equation
journal, January 2005
- Lemou, Mohammed; Mieussens, Luc
- SIAM Journal on Scientific Computing, Vol. 27, Issue 3
Rapid self-magnetization of laser speckles in plasmas by nonlinear anisotropic instability
journal, March 2009
- Thomas, A. G. R.; Kingham, R. J.; Ridgers, C. P.
- New Journal of Physics, Vol. 11, Issue 3
Semi-implicit Time Integration of Atmospheric Flows with Characteristic-Based Flux Partitioning
journal, January 2016
- Ghosh, Debojyoti; Constantinescu, Emil M.
- SIAM Journal on Scientific Computing, Vol. 38, Issue 3
The integration of the vlasov equation in configuration space
journal, November 1976
- Cheng, C. Z.; Knorr, Georg
- Journal of Computational Physics, Vol. 22, Issue 3
Multipole expansions for the Fokker-Planck-Landau operator
journal, February 1998
- Lemou, M.
- Numerische Mathematik, Vol. 78, Issue 4
A Numerical Method for the Accurate Solution of the Fokker–Planck–Landau Equation in the Nonhomogeneous Case
journal, June 2002
- Filbet, Francis; Pareschi, Lorenzo
- Journal of Computational Physics, Vol. 179, Issue 1
High-order, finite-volume methods in mapped coordinates
journal, April 2011
- Colella, P.; Dorr, M. R.; Hittinger, J. A. F.
- Journal of Computational Physics, Vol. 230, Issue 8
A conservative high order semi-Lagrangian WENO method for the Vlasov equation
journal, February 2010
- Qiu, Jing-Mei; Christlieb, Andrew
- Journal of Computational Physics, Vol. 229, Issue 4
A practical difference scheme for Fokker-Planck equations
journal, August 1970
- Chang, J. S.; Cooper, G.
- Journal of Computational Physics, Vol. 6, Issue 1
An Implicit Energy-Conservative 2D Fokker–Planck Algorithm
journal, January 2000
- Chacón, L.; Barnes, D. C.; Knoll, D. A.
- Journal of Computational Physics, Vol. 157, Issue 2
FPPAC: A two-dimensional multispecies nonlinear fokker-planck package
journal, January 1984
- McCoy, M. G.; Mirin, A. A.; Killeen, J.
- Computer Physics Communications, Vol. 35
A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
text, January 2010
- Rossmanith, J. A.; Seal, D. C.
- arXiv
An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
text, January 2011
- Chen, Guangye; Chacón, Luis; Barnes, Daniel C.
- arXiv
New velocity-space discretization for continuum kinetic calculations and Fokker-Planck collisions
text, January 2012
- Landreman, Matt; Ernst, Darin R.
- arXiv
An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm
text, January 2013
- Chen, Guangye; Chacon, Luis
- arXiv
A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
text, January 2015
- Chen, Guangye; Chacon, Luis
- arXiv