TY - JOUR

T1 - Accurate spectral numerical schemes for kinetic equations with energy diffusion

AU - Wilkening, Jon

AU - Cerfon, Antoine J.

AU - Landreman, Matt

N1 - Funding Information:
J.W. was supported in part by the U.S. Department of Energy , Office of Science, Office of Advanced Scientific Computing Research, Applied Mathematics program under contract number DE-AC02-05CH11231 , and by the National Science Foundation under Grant No. DMS-0955078 . A.J.C. was supported by the U.S. Department of Energy , Office of Science, Fusion Energy Sciences under Award Nos. DE-FG02-86ER53223 and DE-SC0012398 . M.L. was supported in part by the U.S. Department of Energy , Office of Science, Office of Fusion Energy Science, under award numbers DE-FG02-93ER54197 and DE-FC02-08ER54964 .
Publisher Copyright:
© 2015.

PY - 2015/8/1

Y1 - 2015/8/1

N2 - We examine the merits of using a family of polynomials that are orthogonal with respect to a non-classical weight function to discretize the speed variable in continuum kinetic calculations. We consider a model one-dimensional partial differential equation describing energy diffusion in velocity space due to Fokker-Planck collisions. This relatively simple case allows us to compare the results of the projected dynamics with an expensive but highly accurate spectral transform approach. It also allows us to integrate in time exactly, and to focus entirely on the effectiveness of the discretization of the speed variable. We show that for a fixed number of modes or grid points, the non-classical polynomials can be many orders of magnitude more accurate than classical Hermite polynomials or finite-difference solvers for kinetic equations in plasma physics. We provide a detailed analysis of the difference in behavior and accuracy of the two families of polynomials. For the non-classical polynomials, if the initial condition is not smooth at the origin when interpreted as a three-dimensional radial function, the exact solution leaves the polynomial subspace for a time, but returns (up to roundoff accuracy) to the same point evolved to by the projected dynamics in that time. By contrast, using classical polynomials, the exact solution differs significantly from the projected dynamics solution when it returns to the subspace. We also explore the connection between eigenfunctions of the projected evolution operator and (non-normalizable) eigenfunctions of the full evolution operator, as well as the effect of truncating the computational domain.

AB - We examine the merits of using a family of polynomials that are orthogonal with respect to a non-classical weight function to discretize the speed variable in continuum kinetic calculations. We consider a model one-dimensional partial differential equation describing energy diffusion in velocity space due to Fokker-Planck collisions. This relatively simple case allows us to compare the results of the projected dynamics with an expensive but highly accurate spectral transform approach. It also allows us to integrate in time exactly, and to focus entirely on the effectiveness of the discretization of the speed variable. We show that for a fixed number of modes or grid points, the non-classical polynomials can be many orders of magnitude more accurate than classical Hermite polynomials or finite-difference solvers for kinetic equations in plasma physics. We provide a detailed analysis of the difference in behavior and accuracy of the two families of polynomials. For the non-classical polynomials, if the initial condition is not smooth at the origin when interpreted as a three-dimensional radial function, the exact solution leaves the polynomial subspace for a time, but returns (up to roundoff accuracy) to the same point evolved to by the projected dynamics in that time. By contrast, using classical polynomials, the exact solution differs significantly from the projected dynamics solution when it returns to the subspace. We also explore the connection between eigenfunctions of the projected evolution operator and (non-normalizable) eigenfunctions of the full evolution operator, as well as the effect of truncating the computational domain.

KW - Continuous spectrum

KW - Continuum kinetic calculations

KW - Fokker-Plank collisions

KW - Orthogonal polynomials

KW - Sturm-Liouville theory

UR - http://www.scopus.com/inward/record.url?scp=84926454980&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84926454980&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2015.03.039

DO - 10.1016/j.jcp.2015.03.039

M3 - Article

AN - SCOPUS:84926454980

VL - 294

SP - 58

EP - 77

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -