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] abuq_v1_0.gz(18 Kbytes)
Manuscript Title: Linear and nonlinear ideal MHD codes - V103.
Authors: H.R. Hicks, J.W. Wooten
Program title: 2LDV103
Catalogue identifier: ABUQ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 13(1977)117
Programming language: PL/1.
Computer: IBM 360/91.
Operating system: OS/360 WITH HASP II.
RAM: 224K words
Word size: 8
Peripherals: disc.
Keywords: Plasma physics, Magnetic Confinement, Stability, Operator codes, Preprocessing, Linear.
Classification: 19.6.

Nature of problem:
The linearized ideal MHD equations are advanced in time. Given an unstable MHD equilibrium, the result is the largest growth rate and the corresponding eigenmode. The geometry is a rectangular cylinder with conducting walls and periodic on the ends.

Solution method:
An explicit leap-frog time advancement scheme is used on a cartesian grid. The finite difference scheme is space and time centered. Eventually, the fastest growing eigenmode dominates.

Only large scale instabilities can be well represented on the finite difference grid. Lower growth rates may be missed due to faster growing numerical instalbilities. The user should realize that large areas of parameter space have not been tested. Although we have run many cases which we believe are valid, we make no representation that results are valid for other than the case reported here.

Unusual features:
Heavy use is made of the IBM PL1 preprocessor.

Running time:
The run time is strongly dependent on the data. The case given in the test deck takes 2.14 min to compile, link edit, and execute on a 360/91.