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] aafs_v1_0.gz(14 Kbytes)
Manuscript Title: STRIMP: program for studying the impurity evolution in tokamak plasma.
Authors: F. Spineanu, M. Vlad, I.I. Popescu
Program title: STRIMP
Catalogue identifier: AAFS_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 41(1986)155
Programming language: Fortran.
Computer: CDC CYBER 170-720.
Operating system: NOS1P4 552/552.
RAM: 25K words
Word size: 60
Keywords: Plasma physics, Tokamak plasma, Transport equations, Impurity, Cubic b-spline, Representation.
Classification: 19.11.

Nature of problem:
The code computes the stationary radial profiles for each impurity stage of ionization and/or their time evolution, taking into account diffusion, ionization and recombination processes occuring in each point of the plasma of a tokamak discharge. This code can be connected with a transport code in order to study the influence of the impurities (e.g. C,O,Fe,...) on the plasma parameters or it can be used to interpret the results of some diagnostic experiments.

Solution method:
The system of partial differential equations representing the transport equations for the density of the impurity ions in each ionization stage is solved using the cubic B-spline approximation of the unknown functions. Restrictions: The tomkamak plasma is assumed to be axisymmetric and to have a circular meridional cross section.

Unusual features:
Due to the uncertainities in the physical theory of impurity diffusion, the code is specially structured to allow an easy implementation of any new transport model. The formulas of the rates of the atomic processes can be modified as well. The code provides the adaptation of the necessary input information and can use different types of boundary conditions.

Running time:
The main computational sequence of the program consisting in loading the matrix of the system and finding the solution needs about 1 second.