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] acjo_v1_0.gz(867 Kbytes)
Manuscript Title: A computer model for resistive MHD analysis.
Authors: M.H. Hughes, M.W. Phillips, R.G. Storer
Program title: ARES
Catalogue identifier: ACJO_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 72(1992)76
Programming language: Fortran.
Computer: IBM RS 6000, MODEL 550.
Operating system: UNIX.
RAM: 5200K words
Word size: 64
Peripherals: disc.
Keywords: Plasma physics, Tokamaks, Stability, Magnetohydrodynamics.
Classification: 19.6.

Subprograms used:
Cat Id Title Reference
ABUF_v2_0 OLYMPUS FOR IBM 370/165 CPC 9(1975)51

Nature of problem:
This paper describes the formulation of a linear MHD stability code in toroidal geometry. The only assumption involved is that the perturb- ations are incompressible. The emphasis is to construct a code capable of calculating the effects of finite plasma resistivity. Thus, a magnetic flux coordinate system is employed since resistivity is known to be important in narrow regions in the vicinity of resonant magnetic surfaces. The code is constructed as an Initial Value problem, using fully implicit temporal differencing to eliminate shortest timescales. This permits 'timesteps' comparable to the mode growth times of interest. Illustrative exapmles of both ideal and resitive instabilities are presented.

Running time:
Problem dependent; the test run required 59 seconds cpu time, including a 25 second overhead for setting up the equilibrium data.