Programs in Physics & Physical Chemistry
|[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|
|Keywords: Plasma physics, Magnetic Confinement, Stability, Operator codes, Preprocessing, Linear.|
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.
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.
Heavy use is made of the IBM PL1 preprocessor.
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.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|