TORBEAM 2.0, a paraxial beam tracing code for electron-cyclotron beams in fusion plasmas for extended physics applications

• Published in: Computer physics communications Vol. 225; no. C; pp. 36 - 46
• Main Authors: Poli, E., Bock, A., Lochbrunner, M., Maj, O., Reich, M., Snicker, A., Stegmeir, A., Volpe, F., Bertelli, N., Bilato, R., Conway, G.D., Farina, D., Felici, F., Figini, L., Fischer, R., Galperti, C., Happel, T., Lin-Liu, Y.R., Marushchenko, N.B., Mszanowski, U., Poli, F.M., Stober, J., Westerhof, E., Zille, R., Peeters, A.G., Pereverzev, G.V.
• Format: Journal Article
• Language: English
• Published: United States Elsevier B.V 01.04.2018; Elsevier
• Subjects: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Electron cyclotron waves; Magnetic confinement; MATHEMATICS AND COMPUTING; Paraxial beam tracing; Plasma physics; Wave–plasma interactions; Wave–plasma interactions; Plasma physics; Paraxial beam tracing; Magnetic confinement; Electron cyclotron waves
• ISSN: 0010-4655; 1879-2944
• DOI: 10.1016/j.cpc.2017.12.018

The paraxial WKB code TORBEAM (Poli, 2001) is widely used for the description of electron-cyclotron waves in fusion plasmas, retaining diffraction effects through the solution of a set of ordinary differential equations. With respect to its original form, the code has undergone significant transformations and extensions, in terms of both the physical model and the spectrum of applications. The code has been rewritten in Fortran 90 and transformed into a library, which can be called from within different (not necessarily Fortran-based) workflows. The models for both absorption and current drive have been extended, including e.g. fully-relativistic calculation of the absorption coefficient, momentum conservation in electron–electron collisions and the contribution of more than one harmonic to current drive. The code can be run also for reflectometry applications, with relativistic corrections for the electron mass. Formulas that provide the coupling between the reflected beam and the receiver have been developed. Accelerated versions of the code are available, with the reduced physics goal of inferring the location of maximum absorption (including or not the total driven current) for a given setting of the launcher mirrors. Optionally, plasma volumes within given flux surfaces and corresponding values of minimum and maximum magnetic field can be provided externally to speed up the calculation of full driven-current profiles. These can be employed in real-time control algorithms or for fast data analysis.