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] aept_v1_0.tar.gz(1074 Kbytes)
Manuscript Title: GTROTA: A code for the solution of the coupled nonlinear extended neoclassical rotation equations in tokamak plasmas using successive over-relaxation and simulated annealing
Authors: C. Bae, W.M. Stacey, T.D. Morley
Program title: GTROTAv1
Catalogue identifier: AEPT_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 184(2013)2571
Programming language: Matlab.
Computer: Any workstation or PC where Matlab can be run.
Operating system: Any with Matlab available. Tested on Windows 7.
RAM: 1024 MB to run Matlab
Keywords: Nonlinear, Tokamak, Plasma, GTROTA.
Classification: 19.11.

Nature of problem:
Tokamak plasma rotation velocities and the poloidal in-out and up-down asymmetries of plasma densities are calculated from the coupled set of nonlinear equations, which shows a very instable iterative dynamics. To solve for the true solution for this nonlinearly coupled system that does not converge to a single solution, physics-based determination of the true solution becomes necessary and GTROTA provides the algorithm for users to find the true solution from the nonlinear topological maps.

Solution method:
The code decomposes the given system into three subsystems to stabilize the iterative dynamics and uses Successive Over-Relaxation and the concept of Simulated Annealing to determine the true solution.

Restrictions:
The code is designed for strong rotation analysis. An updated code for slow rotation and plasma edge analysis is to be developed in the future.

Unusual features:
For a new rotation calculation, the code requires several test runs to determine the iterations step for the true solution.

Additional comments:
Due to the extreme nonlinearity of the extended rotation calculation model, the signs of the inputs must be adjusted correctly to yield correct solutions. Refer to the User's Manual for this discussion.

Running time:
8 seconds (except test runs that depend on the complexity of the nonlinear dynamics)