Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators
Abstract
Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy tomore »
- Authors:
-
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Plasma Science and Fusion Center; Dartmouth College, Hanover, NH (United States)
- Univ. of Maryland, College Park, MD (United States)
- Publication Date:
- Research Org.:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC); US Air Force Office of Scientific Research (AFOSR); National Aeronautics and Space Administration (NASA); National Science Foundation (NSF)
- OSTI Identifier:
- 1668773
- Grant/Contract Number:
- AC02-06CH11357; FG02-91-ER54109; AC02-09CH11466; SC-0010508; FC02-08ER54966; 80NSSC17K0428; ACI-1548562; FA9550-15-1-0193
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Plasma Physics
- Additional Journal Information:
- Journal Volume: 86; Journal Issue: 4; Journal ID: ISSN 0022-3778
- Publisher:
- Cambridge University Press
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; plasma dynamics; plasma nonlinear phenomena; plasma simulation
Citation Formats
Hakim, Ammar, Francisquez, Manaure, Juno, James, and Hammett, Gregory W. Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators. United States: N. p., 2020.
Web. doi:10.1017/s0022377820000586.
Hakim, Ammar, Francisquez, Manaure, Juno, James, & Hammett, Gregory W. Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators. United States. https://doi.org/10.1017/s0022377820000586
Hakim, Ammar, Francisquez, Manaure, Juno, James, and Hammett, Gregory W. Fri .
"Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators". United States. https://doi.org/10.1017/s0022377820000586. https://www.osti.gov/servlets/purl/1668773.
@article{osti_1668773,
title = {Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators},
author = {Hakim, Ammar and Francisquez, Manaure and Juno, James and Hammett, Gregory W.},
abstractNote = {Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.},
doi = {10.1017/s0022377820000586},
journal = {Journal of Plasma Physics},
number = 4,
volume = 86,
place = {United States},
year = {Fri Jul 17 00:00:00 EDT 2020},
month = {Fri Jul 17 00:00:00 EDT 2020}
}
Works referenced in this record:
The Mathematical Theory of Finite Element Methods
book, January 2008
- Brenner, Susanne C.; Scott, L. Ridgway
- Texts in Applied Mathematics
Discontinuous Galerkin algorithms for fully kinetic plasmas
journal, January 2018
- Juno, J.; Hakim, A.; TenBarge, J.
- Journal of Computational Physics, Vol. 353
Implications of advanced collision operators for gyrokinetic simulation
journal, February 2017
- Belli, E. A.; Candy, J.
- Plasma Physics and Controlled Fusion, Vol. 59, Issue 4
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
Model Fokker-Planck Equation for a Plasma and Its Solution
journal, January 1964
- Dougherty, J. P.
- Physics of Fluids, Vol. 7, Issue 11
Effect of a weak ion collisionality on the dynamics of kinetic electrostatic shocks
journal, February 2019
- Sundström, Andréas; Juno, James; TenBarge, Jason M.
- Journal of Plasma Physics, Vol. 85, Issue 1
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
Improved Bhatnagar-Gross-Krook model of electron-ion collisions
journal, January 1973
- Greene, John M.
- Physics of Fluids, Vol. 16, Issue 11
Gyrokinetic continuum simulations of plasma turbulence in the Texas Helimak
journal, April 2019
- Bernard, T. N.; Shi, E. L.; Gentle, K. W.
- Physics of Plasmas, Vol. 26, Issue 4
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
Improved collision operator for plasma kinetic simulations with multi-species ions and electrons
journal, December 2015
- Nakata, Motoki; Nunami, Masanori; Watanabe, Tomo-Hiko
- Computer Physics Communications, Vol. 197
Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma
journal, May 2017
- Shi, E. L.; Hammett, G. W.; Stoltzfus-Dueck, T.
- Journal of Plasma Physics, Vol. 83, Issue 3
Eigenfunctions and eigenvalues of the Dougherty collision operator
journal, May 2007
- Anderson, M. W.; O’Neil, T. M.
- Physics of Plasmas, Vol. 14, Issue 5
Heating of a Confined Plasma by Oscillating Electromagnetic Fields
journal, January 1958
- Berger, J. M.; Newcomb, W. A.; Dawson, J. M.
- Physics of Fluids, Vol. 1, Issue 4
Conservative discretization of the Landau collision integral
journal, March 2017
- Hirvijoki, E.; Adams, M. F.
- Physics of Plasmas, Vol. 24, Issue 3
Magnetic Pumping as a Source of Particle Heating and Power-law Distributions in the Solar Wind
journal, November 2017
- Lichko, E.; Egedal, J.; Daughton, W.
- The Astrophysical Journal, Vol. 850, Issue 2
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
Transition between Landau and Van Kampen Treatments of the Vlasov Equation
journal, January 1967
- Grant, Frederick C.
- Physics of Fluids, Vol. 10, Issue 6
The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems
journal, December 1998
- Cockburn, Bernardo; Shu, Chi-Wang
- SIAM Journal on Numerical Analysis, Vol. 35, Issue 6
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
Full- f gyrokinetic simulation of turbulence in a helical open-field-line plasma
journal, January 2019
- Shi, E. L.; Hammett, G. W.; Stoltzfus-Dueck, T.
- Physics of Plasmas, Vol. 26, Issue 1
The Serendipity Family of Finite Elements
journal, March 2011
- Arnold, Douglas N.; Awanou, Gerard
- Foundations of Computational Mathematics, Vol. 11, Issue 3
A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma
journal, June 2016
- Hager, Robert; Yoon, E. S.; Ku, S.
- Journal of Computational Physics, Vol. 315
A drift-kinetic analytical model for scrape-off layer plasma dynamics at arbitrary collisionality
journal, November 2017
- Jorge, R.; Ricci, P.; Loureiro, N. F.
- Journal of Plasma Physics, Vol. 83, Issue 6
Theory of first‐order plasma heating by collisional magnetic pumping
journal, May 1989
- Laroussi, M.; Roth, J. Reece
- Physics of Fluids B: Plasma Physics, Vol. 1, Issue 5
Collisional relaxation: Landau versus Dougherty operator
journal, October 2014
- Pezzi, Oreste; Valentini, F.; Veltri, P.
- Journal of Plasma Physics, Vol. 81, Issue 1
Plasma Instabilities and Magnetic Field Growth in Clusters of Galaxies
journal, August 2005
- Schekochihin, A. A.; Cowley, S. C.; Kulsrud, R. M.
- The Astrophysical Journal, Vol. 629, Issue 1
Plasma Oscillations with Diffusion in Velocity Space
journal, December 1958
- Lenard, A.; Bernstein, Ira B.
- Physical Review, Vol. 112, Issue 5
Developments in the gyrofluid approach to Tokamak turbulence simulations
journal, August 1993
- Hammett, G. W.; Beer, M. A.; Dorland, W.
- Plasma Physics and Controlled Fusion, Vol. 35, Issue 8
Linearized model Fokker–Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests
journal, July 2009
- Barnes, M.; Abel, I. G.; Dorland, W.
- Physics of Plasmas, Vol. 16, Issue 7