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] acdy_v1_0.gz(14 Kbytes)
Manuscript Title: TOPIC: a tokamak plasma impurities code.
Authors: T.A. Beu, F. Spineanu, M. Vlad, R.I. Campeanu, I.I. Popescu
Program title: TOPIC
Catalogue identifier: ACDY_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 36(1985)161
Programming language: Fortran.
Computer: CDC CYBER 170-720.
Operating system: NOS1P4 552/552.
RAM: 1K words
Word size: 64
Peripherals: disc.
Keywords: Tokamak, Plasma physics, Impurities, Diffusion equation, Coupled differential Equations, Three-diagonal block Matrix inversion, Magnetic confinement.
Classification: 19.9.

Nature of problem:
In high temperature low density tokamak plasmas, radiation cooling by impurity atoms can be an important energy loss mechanism, since radiation is not reabsorbed. The coupled set of time-dependent diffusion equations for ionized tokamak plasma impurities is solved in conjunction with a simple model for neutral impurities, on the assumption of a cylindrical symmetry. The resultant density distributions are used in the subsequent computation of the related energy losses.

Solution method:
The set of coupled second order differential equations, describing the diffusion of impurities in the MHD approach, is discretized using a finite difference scheme. The resultant linear system, with a three- diagonal block matrix structure, is then solved by a special forward elimination - backward substitution technique.

Running time:
About 6 s for the test run ouput.