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] aafz_v1_0.gz(104 Kbytes)
Manuscript Title: ODRIC: a one-dimensional linear resistive MHD code in cylindrical geometry.
Authors: A.A. Mirin, R.J. Bonugli, N.J. O'Neill, J. Killeen
Program title: ODRIC
Catalogue identifier: AAFZ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 41(1986)85
Programming language: Fortran.
Computer: CRAY-1, X-MP.
Operating system: CTSS.
Word size: 64
Keywords: Plasma physics, One-dimensional, Resistive, Implicit, Linear, Stability.
Classification: 19.6.

Nature of problem:
The primitive linearized resistive MHD equations in cylindrical geometry are advanced in time. Separate equations for the electron and ion temperature perturbations are solved. Hall terms and the thermal force vector are included in Ohm's law. Anisotropic thermal conductivity and viscosity are included in the model. The plasma is assumed to be quasineutral. Fourier analysis in the poloidal and toroidal directions is performed, resulting in a 1-D (radial) system for the perturbed quantities.

Solution method:
A finite difference method is used. Spatial derivatives are approximated by central differences. The mesh need not be uniform. Temporal differences are fully implicit.

The user must adjust the size and type of radial mesh, in particular in the singular layer, in order to model the physics accurately. For nonaxisymmetric perturbations (m ^= 0; see Eq. (21)) a small amount of viscosity is necessary to prevent numerical instability. This is discussed further.

Running time:
Two versions of the code are provided (see Sec.1 of long write-up). The real version takes 40 ms per meshpoint per timestep on the Cray-1. The complex version takes 113 ms per meshpoint per timestep.