Multi-dimensional, fully-implicit, spectral method for the Vlasov-Maxwell equations with exact conservation laws in discrete form

JOURNAL OF COMPUTATIONAL PHYSICS

卷:301

Article

Delzanno, G. L.

关键词: Vlasov-Maxwell equations;    Kinetic plasma theory;    Spectral method;    Fourier-Hermite basis functions;    Fully implicit Crank-Nicolson time discretization;    Jacobian-Free Newton-Krylov solver;   

DOI  :  10.1016/j.jcp.2015.07.028

来源: Elsevier

【 摘 要 】

A spectral method for the numerical solution of the multi-dimensional Vlasov-Maxwell equations is presented. The plasma distribution function is expanded in Fourier (for the spatial part) and Hermite (for the velocity part) basis functions, leading to a truncated system of ordinary differential equations for the expansion coefficients (moments) that is discretized with an implicit, second order accurate Crank-Nicolson time discretization. The discrete non-linear system is solved with a preconditioned Jacobian-Free Newton-Krylov method. It is shown analytically that the Fourier-Hermite method features exact conservation laws for total mass, momentum and energy in discrete form. Standard tests involving plasma waves and the whistler instability confirm the validity of the conservation laws numerically. The whistler instability test also shows that we can step over the fastest time scale in the system without incurring in numerical instabilities. Some preconditioning strategies are presented, showing that the number of linear iterations of the Krylov solver can be drastically reduced and a significant gain in performance can be obtained. (C) 2015 Elsevier Inc. All rights reserved.

【 授权许可】

Free