A numerical instability in an ADI algorithm for gyrokinetics

We explore the implementation of an Alternating Direction Implicit (ADI) algorithm for a gyrokinetic plasma problem and its resulting numerical stability properties. This algorithm, which uses a standard ADI scheme to divide the field solve from the particle distribution function advance, has previously been found to work well for certain plasma kinetic problems involving one spatial and two velocity dimensions, including collisions and an electric field. However, for the gyrokinetic problem we find a severe stability restriction on the time step. Furthermore, we find that this numerical instability limitation also affects some other algorithms, such as a partially implicit Adams-Bashforth algorithm, where the parallel motion operator v{sub {parallel}} {partial_derivative}/{partial_derivative}z is treated implicitly and the field terms are treated with an Adams-Bashforth explicit scheme. Fully explicit algorithms applied to all terms can be better at long wavelengths than these ADI or partially implicit algorithms.

[1]  Frank Jenko,et al.  Electron temperature gradient driven turbulence , 1999 .

[2]  Frank Jenko,et al.  Massively parallel Vlasov simulation of electromagnetic drift-wave turbulence , 2000 .

[3]  B. Cohen,et al.  Kinetic electron closures for electromagnetic simulation of drift and shear-Alfvén waves. II , 2001 .

[4]  O. Sauter,et al.  Kinetic modeling of scrape‐off layer plasmas , 1996 .

[5]  E. Frieman,et al.  Nonlinear gyrokinetic equations for low-frequency electromagnetic waves in general plasma equilibria , 1981 .

[6]  Scott E. Parker,et al.  A δ f particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations , 2003 .

[7]  H. H. Rachford,et al.  On the numerical solution of heat conduction problems in two and three space variables , 1956 .

[8]  T. Antonsen,et al.  Kinetic equations for low frequency instabilities in inhomogeneous plasmas , 1980 .

[9]  Frank Jenko,et al.  Vlasov simulation of kinetic shear Alfvén waves , 2004, Comput. Phys. Commun..

[10]  D. Durran.  Numerical methods for wave equations in geophysical fluid dynamics , 1999 .

[11]  F. Jenko,et al.  Electron temperature gradient turbulence. , 2000, Physical review letters.

[12]  R. E. Waltz,et al.  An Eulerian gyrokinetic-Maxwell solver , 2003 .

[13]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[14]  R. Waltz,et al.  Anomalous transport scaling in the DIII-D tokamak matched by supercomputer simulation. , 2003, Physical review letters.

[15]  B. M. Fulk.  MATH , 1992 .

[16]  W. Dorland,et al.  Fluid models of phase mixing, Landau damping, and nonlinear gyrokinetic dynamics , 1992 .

[17]  W. Lee,et al.  Gyrokinetic Particle Simulation Model , 1987 .

[18]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[19]  Willem Hundsdorfer,et al.  A note on stability of the Douglas splitting method , 1998, Math. Comput..

[20]  Yang Chen,et al.  A gyrokinetic ion zero electron inertia fluid electron model for turbulence simulations , 2001 .