Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] acvl_v1_0.gz(24 Kbytes)
Manuscript Title: A fast 2-D Fokker-Planck solver with synergetic effects.
Authors: M. Shoucri, I. Shkarofsky
Program title: FASTFP
Catalogue identifier: ACVL_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 82(1994)287
Programming language: Fortran, Postscript.
Computer: SUN SPARC10.
Operating system: UNIX.
RAM: 21K words
Word size: 32
Keywords: Plasma physics, Kinetic model, Fokker-planck, D03uaf, Current drive, Synergetic effects.
Classification: 19.8.

Nature of problem:
We present in this work a fast solver for the linearized 2-D relativistic Fokker-Planck equation for current drive problems. The code makes use of the accurate relativistic collision operators presented in Ref. [14]. Convergence to the steady state can be reached in about 30 iterations. Very good precision for the solution of the distribution function can be obtained for high relativistic velocities. In the case of fast wave current drive problems, the damping of the wave is calculated self-consistently from the solution of the 2-D Fokker- Planck equation.

Solution method:
The linearized relativistic Fokker-Planck equation is solved on a 2-D domain. The 2-D equation is discretized using an implicit 5 point stencil, and solved using the NAG library subroutine D03UAF. This application also provides a powerful example on how to solve a partial differential equation on a 2-D domain with the subroutine D03UAF.

Running time:
An execution on a SPARC10 workstation requires about 70 microsec/grid point/time step. (The solution of a typical problem requires 200 x 100 grid points and about 30 iterations.)