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] acgt_v1_0.gz(41 Kbytes)
Manuscript Title: ECRCYL: a code for electron cyclotron radiation transport.
Authors: S. Chaturvedi, R.G. Mills
Program title: ECRCYL
Catalogue identifier: ACGT_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 70(1992)183
Programming language: Fortran.
Computer: VAX-8810.
Operating system: VMS v5.1.
RAM: 2578K words
Word size: 32
Keywords: Plasma physics, Cyclotron radiation, Transport radiation, Ray tracing, Synchrotron radiation, Fusion, Confinement magnetic.
Classification: 19.11.

External routines: (1) Routine ODE.
(2) VAX VMS run-time routine LIB$DATE_TIME.
(3) Function "Z" - Plasma dispersion function.
(4) IMSL routine MMDEI (Bessel function calculation).

Nature of problem:
The code performs electron cyclotron radiation transport calculations in a magnetized cylindrical plasma, taking account of plasma dispersion and modification of wave absorption due to finite plasma density. High accuracy calculations can be performed. With minor modifications, the code can handle the effects of spatial variations of wall reflectivity, and 'polarization scrambling' due to reflection at the wall.

Solution method:
Ray tracing is performed, using the cold plasma dispersion relation to determine ray trajectories. The expressions for absorption coefficients given by Bornatici et al. are used to calculate power deposition for different frequencies, angles of propagation and modes.

Restrictions:
(1) The plasma lies in the annular space between two concentric cylinders, and is bounded at both ends. The magnetic field is parallel to the circumference of the cylinders.
(2) The line-width of radiation about a harmonic must be significantly smaller than the electron cyclotron frequency omega ce.
(3) Reflection from the walls is ignored in the present form of the code.
(4) Electrons in the plasma are assumed to follow a Maxwellian distribution of velocities, so that the emission and absorption coefficients are related by Kirchhoff's law.

Running time:
Approximately 45 minutes on a VAX-8810 for the test job. The running time depends strongly on the mesh size used, and on the absorption mean free path for photons.