Adaptive spectral solution method for the Landau and Lenard-Balescu equations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:402

Article

Scullard, Christian R.¹  Hickok, Abigail^2,3  Sotiris, Justyna O.^2,4  Tzolova, Bilyana M.^2,5  Van Heyningen, R. Loek^2,6  Graziani, Frank R.¹ 

[1] Lawrence Livermore Natl Lab, Livermore, CA 94550 USA

[2] Univ Calif Los Angeles, Inst Pure & Appl Math, Los Angeles, CA 90095 USA

[3] Univ Calif Los Angeles, Dept Math, Los Angeles, CA 90095 USA

[4] Univ New Mexico, Dept Math & Stat, MSC01 1115, Albuquerque, NM 87131 USA

[5] Johns Hopkins Univ, Dept Math, Baltimore, MD 21218 USA

[6] Univ Calif Berkeley, Dept Math, Berkeley, CA 94720 USA

关键词: Kinetic theory;    Landau equation;    Lenard-Balescu equation;   

DOI  :  10.1016/j.jcp.2019.109110

来源: Elsevier

【 摘 要 】

We present an adaptive spectral method for solving the Landau/Fokker-Planck equation for electron-ion systems. The heart of the algorithm is an expansion in Laguerre polynomials, which has several advantages, including automatic conservation of both energy and particles without the need for any special discretization or time-stepping schemes. One drawback of such an expansion is the O (N-3) memory requirement, where N is the number of polynomials used. This can impose an inconvenient limit in cases of practical interest, such as when two particle species have widely separated temperatures. The algorithm we describe here addresses this problem by periodically re-projecting the solution onto a judicious choice of new basis functions that are still Laguerre polynomials but have arguments adapted to the current physical conditions. This results in a reduction in the number of polynomials needed, at the expense of increased solution time. Because the equations are solved with little difficulty, this added time is not of much concern compared to the savings in memory. To demonstrate the algorithm, we solve several relaxation problems that could not be computed with the spectral method without re-projection. Another major advantage of this method is that it can be used for collision operators more complicated than that of the Landau equation, and we demonstrate this here by using it to solve the non-degenerate quantum Lenard-Balescu (QLB) equation for a hydrogen plasma. We conclude with some comparisons of temperature relaxation problems solved with the latter equation and the Landau equation with a Coulomb logarithm inspired by the properties of the QLB operator. We find that with this choice of Coulomb logarithm, there is little difference between using the two equations for these particular systems. (C) 2019 Elsevier Inc. All rights reserved.

【 授权许可】

Free