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] acpk_v1_0.gz(6 Kbytes)
Manuscript Title: A Fokker-Planck code for the numerical solution of the plasma heating and the current drive problems with synergetic effects (FW/EC - FW/LH - LH/EC).
Authors: M. Shoucri, I. Shkarofsky
Program title: SOLV2DFP
Catalogue identifier: ACPK_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 78(1993)199
Programming language: Fortran, LaTeX.
Computer: IBM 3092.
Operating system: VAX.
RAM: 9999K words
Keywords: Plasma physics, Kinetic model, Fokker-planck, Pde-protran, Pde2d, Current drive, Synergetic effects.
Classification: 19.8.

Nature of problem:
A code is presented that solves numerically the linearized relativistic 2-D Fokker-Planck equation for the problem of plasma heating or current drive using lower hybrid waves, fast waves, electron cyclotron waves, as well as the synergetic effects of these waves. The code uses the IMSL finite element library PDE-PROTRAN, or its new version PDE2D. Very good precision for the solution of the distribution function can be obtained for high relativistic velocities, which allows the hot population to be calculated accurately. In the case of the fast wave current drive problem, 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 code solves for the equilibrium solution. The code defining the input parameters and the 2-D differential equation for the library used requires about 170 FORTRAN lines. This application also provides a powerful example on how to use PDE-PROTRAN or PDE2D to solve a partial differential equation in a 2-D domain.

Running time:
An execution on an IBM 3092 using 20 000 quadratic elements (which gives very high precision to high relativistic velocities), would require an execution CPU time of about 30 minutes. A few minutes would be necessary on a CRAY machine.