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

Abstract 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.

[1]  Nikolai B. Marushchenko,et al.  Corrections to the paper ‘electron cyclotron current drive calculated for ITER conditions using different models’ (2008 Nucl. Fusion 48 054002) , 2009 .

[2]  M. Reich,et al.  Real-time beam tracing for control of the deposition location of electron cyclotron waves , 2015 .

[3]  Y. Kravtsov,et al.  Geometrical optics of inhomogeneous media , 2019, Geometrical Optics of Weakly Anisotropic Media.

[4]  E. Strumberger,et al.  Analysis of the ITER low field side reflectometer employing the Beam Tracing Method , 2011 .

[5]  Timothy Goodman,et al.  On recent results in the modelling of neoclassical-tearing-mode stabilization via electron cyclotron current drive and their impact on the design of the upper EC launcher for ITER , 2015 .

[6]  R. W. Harvey,et al.  Benchmarking of codes for electron cyclotron heating and electron cyclotron current drive under ITER conditions , 2008 .

[7]  T. Görler,et al.  Core turbulence behavior moving from ion-temperature-gradient regime towards trapped-electron-mode regime in the ASDEX Upgrade tokamak and comparison with gyrokinetic simulation , 2015 .

[8]  Kjell Rönnmark Kinetic theory of plasma waves , 1985 .

[9]  Jet Contributors,et al.  Stationary Zonal Flows during the Formation of the Edge Transport Barrier in the JET Tokamak. , 2016, Physical review letters.

[10]  Stefano Coda,et al.  Distributed digital real-time control system for the TCV tokamak and its applications , 2017 .

[11]  H. Kogelnik,et al.  Coupling and Conversion Coefficients for Optical Modes , 1964 .

[12]  M. Brambilla,et al.  Feasibility of electron Bernstein wave coupling via O-X-B mode conversion in the RFX-mod reversed field pinch device , 2009 .

[13]  Lorenzo Figini,et al.  Electron-cyclotron-current-drive efficiency in DEMO plasmas , 2012 .

[14]  O. Sauter,et al.  EU developments of the ITER ECRH system , 2007 .

[15]  J. Stober,et al.  Coupling of the Flux Diffusion Equation with the Equilibrium Reconstruction at ASDEX Upgrade , 2016 .

[16]  E. Westerhof,et al.  Relax a computer code for the study of collisional and wave driven relaxation of the electron distribution function in toroidal geometry , 1992 .

[17]  G. V. Pereverzev,et al.  TORBEAM, a beam tracing code for electron-cyclotron waves in tokamak plasmas , 2001 .

[18]  Jet Efda Contributors,et al.  The European Integrated Tokamak Modelling (ITM) effort: achievements and first physics results , 2014 .

[19]  V. Erckmann,et al.  Electron cyclotron resonance heating and current drive in the W7-X stellarator , 1998 .

[20]  Nikolai B. Marushchenko,et al.  Electron cyclotron current drive calculated for ITER conditions using different models , 2008 .

[21]  H. Zohm,et al.  Recent ASDEX Upgrade research in support of ITER and DEMO , 2015 .

[22]  Roddy G. L. Vann,et al.  2D Doppler backscattering using synthetic aperture microwave imaging of MAST edge plasmas , 2016 .

[23]  Francis F. Chen Plasma Physics and Controlled Fusion , 1984 .

[24]  Winfried Kernbichler,et al.  Electron cyclotron current drive in low collisionality limit: On parallel momentum conservation , 2011 .

[25]  Omar Maj,et al.  The wave energy flux of high frequency diffracting beams in complex geometrical optics , 2013 .

[26]  G. V. Pereverzev,et al.  Validation of the paraxial beam-tracing method in critical cases , 2009 .

[27]  S. Zweben Principles of Plasma Diagnostics , 1990 .

[28]  Omar Maj,et al.  Scattering of diffracting beams of electron cyclotron waves by random density fluctuations in inhomogeneous plasmas , 2015 .

[29]  D. G. Swanson,et al.  Chapter 4 – Kinetic Theory of Plasma Waves , 1989 .

[30]  Henrik Bindslev,et al.  Relativistic expressions for plasma cutoffs , 1993 .

[31]  Lorenzo Figini,et al.  Assessment of the ITER electron cyclotron upper launcher capabilities in view of an optimized design , 2015 .

[32]  F. Felici,et al.  Real-time physics-model-based simulation of the current density profile in tokamak plasmas , 2011 .

[33]  G. V. Pereverzev,et al.  Use of the multidimensional WKB method to describe propagation of lower hybrid waves in tokamak plasma , 1992 .

[34]  T. C. Luce,et al.  Electron cyclotron current drive efficiency in general tokamak geometry , 2003 .

[35]  P. T. Bonoli,et al.  Paraxial Wentzel–Kramers–Brillouin method applied to the lower hybrid wave propagation , 2012 .

[36]  C. Fuchs,et al.  Integrated Data Analysis of Profile Diagnostics at ASDEX Upgrade , 2010 .

[37]  Thomas M. Antonsen,et al.  Radio frequency current generation by waves in toroidal geometry , 1982 .

[38]  W. Treutterer,et al.  Performance, heating and current drive scenarios of ASDEX Upgrade advanced tokamak discharges , 2001 .

[39]  F. Felici,et al.  Non-linear model-based optimization of actuator trajectories for tokamak plasma profile control , 2012 .

[40]  Omar Maj,et al.  Effects of aberration on paraxial wave beams: beam tracing versus quasi-optical solutions , 2010 .

[41]  Daniela Farina,et al.  A Quasi-Optical Beam-Tracing Code for Electron Cyclotron Absorption and Current Drive: GRAY , 2007 .

[42]  E Egbert Westerhof,et al.  Wave propagation through an electron cyclotron resonance layer , 1997 .