IAEA0909 TOPIC-RUM.
TOPIC-RUM, Plasma Impurities in Tokamak Reactor by MHD Method
NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES
[ top ]
[ top ]
To submit a request, click below on the link of the version you wish to order.
Only liaison officers are authorised to submit online requests. Rules for requesters are
available here.
| Program name | Package id | Status | Status date |
|---|---|---|---|
| TOPIC-RUM | IAEA0909/01 | Tested | 09-MAY-1989 |
Machines used:
| Package ID | Orig. computer | Test computer |
|---|---|---|
| IAEA0909/01 | CDC CYBER 720 | CDC CYBER 830 |
[ top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION
In high temperature-low density tokamak plasmas, radiation cooling by impurity atoms can be an important energy loss mechanism, since radiation is not reabsorbed. The coupled set of time-dependent diffusion equations for ionized plasma impurities is solved in conjunction with a simple model for neutral impurities, on the assumption of a cylindrical symmetry. The resultant density distributions are used in the subsequent computation of the related energy losses.
In high temperature-low density tokamak plasmas, radiation cooling by impurity atoms can be an important energy loss mechanism, since radiation is not reabsorbed. The coupled set of time-dependent diffusion equations for ionized plasma impurities is solved in conjunction with a simple model for neutral impurities, on the assumption of a cylindrical symmetry. The resultant density distributions are used in the subsequent computation of the related energy losses.
[ top ]
4. METHOD OF SOLUTION
The set of coupled second order equations, describing the diffusion of impurities in the MHD approach, is discretized using a finite difference scheme. The resultant linear system, with a three-diagonal block matrix structure, is then solved by a special forward elimination backward substitution technique.
The set of coupled second order equations, describing the diffusion of impurities in the MHD approach, is discretized using a finite difference scheme. The resultant linear system, with a three-diagonal block matrix structure, is then solved by a special forward elimination backward substitution technique.
[ top ]
[ top ]
IAEA0909/01
NEA-DB ran the two test cases included in this package on a CDC CYBER 830 computer in less than 13 seconds of CPU time.[ top ]
[ top ]
[ top ]
IAEA0909/01, included references:
- T.A. Beu, F. Spineanu, M. Vlad, R.I. Campeanu and I.I. Popescu:TOPIC - A Tokamak Plasma Impurities Code.
Reprint from "Computer Physics Communication", Vol. 36(2),
pp. 161-176 (April 1985)
- User's Guide of the TOPIC Program (Detailed Input Specification).
[ top ]
11. MACHINE REQUIREMENTS
High speed storage required : 0.7 Kwords for the test run. Peripherals used: card reader, line printer, magnetic disk.
High speed storage required : 0.7 Kwords for the test run. Peripherals used: card reader, line printer, magnetic disk.
IAEA0909/01
55,000 (octal) words of CDC CYBER 830 main storage.[ top ]
IAEA0909/01
NOS2.5.1 (CDC CYBER 830).[ top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS
The number of words required by the arrays used in the program, NW, is used by the system subprograms SYMLIB, CNMALF and LOCF for dynamic memory allocation, relative to the starting adress of the array A dimensioned in the main program. On computers which do not have dynamic memory allocation facilities the sequence
CALL SYMLIB
CALL CMMALF (NW+1,0,0,IA)
LA=IA-LOCF (A(1))+1
should be replaced by
LA=1
while the array A should be dimensioned according to NW.
The number of words required by the arrays used in the program, NW, is used by the system subprograms SYMLIB, CNMALF and LOCF for dynamic memory allocation, relative to the starting adress of the array A dimensioned in the main program. On computers which do not have dynamic memory allocation facilities the sequence
CALL SYMLIB
CALL CMMALF (NW+1,0,0,IA)
LA=IA-LOCF (A(1))+1
should be replaced by
LA=1
while the array A should be dimensioned according to NW.
[ top ]
15. NAME AND ESTABLISHMENT OF AUTHORS
T.A. Beu, F. Spineanu, M. Vlad
Central Institute of Physics
Bucharest-Magurele, Romania
R.I. Campeanu
University of Cluj-Napoca
Department of physics
3400 Cluj-Napoca, Romania
I.I. Popescu
University of Bucharest
Department of Physics
Bucharest-Magurele, Romania
T.A. Beu, F. Spineanu, M. Vlad
Central Institute of Physics
Bucharest-Magurele, Romania
R.I. Campeanu
University of Cluj-Napoca
Department of physics
3400 Cluj-Napoca, Romania
I.I. Popescu
University of Bucharest
Department of Physics
Bucharest-Magurele, Romania
[ top ]
IAEA0909/01
| File name | File description | Records |
|---|---|---|
| IAEA0909_01.001 | Information file | 46 |
| IAEA0909_01.002 | JCL and control information | 11 |
| IAEA0909_01.003 | TOPIC-RUM source program | 1128 |
| IAEA0909_01.004 | Sample case input 1 | 35 |
| IAEA0909_01.005 | Sample case input 2 | 44 |
| IAEA0909_01.006 | Sample output | 443 |
Keywords: diffusion equations, magnetohydrodynamics, plasma, tokamak reactors.