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.
- Research Organization:
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC) (US)
- DOE Contract Number:
- AC02-76CH03073
- OSTI ID:
- 836494
- Report Number(s):
- PPPL-4034; TRN: US0500615
- Resource Relation:
- Other Information: PBD: 17 Dec 2004
- Country of Publication:
- United States
- Language:
- English
Similar Records
Center for Gyrokinetic/MHD Hybrid Simulation of Energetic Particle Physics in Toroidal Plasmas (CSEPP). Final report
Numerical Stability and Accuracy of Temporally Coupled Multi-Physics Modules in Wind-Turbine CAE Tools