BORAY: An Axisymmetric Ray Tracing Code Supports Both Closed and Open Field Lines Plasmas

Hua-sheng Xie,Banerjee Debabrata,Yu-kun Bai,Han-yue Zhao,Jing-chun Li

DOI: https://doi.org/10.1016/j.cpc.2022.108363

2021-05-25

Abstract:Ray tracing codes are useful to study the electromagnetic wave propagation and absorption in the geometrical optics approximation. In magnetized fusion plasma community, most ray tracing codes assume the plasma density and temperature be functions of the magnetic flux and study waves only inside the last closed flux surface, which are sufficient for the present day tokamak. However, they are difficult to be used for configurations with open magnetic field line plasmas, such as mirror machine and field-reversed-configuration (FRC). We develop a ray tracing code in cylindrical coordinates $(r,\phi,z)$ to support arbitrary axisymmetric configurations with both closed and open field lines plasmas. For wave propagation, the cold plasma dispersion relation is usually sufficient, and we require the magnetic field ${\bf B}(r,z)$ and species densities $n_{s0}(r,z)$ profiles as input. For wave absorption, we require a further temperature $T_{s0}(r,z)$ profile to solve a hot kinetic plasma dispersion relation. In difference to other ray tracing codes which calculate the imaginary part of wave vector ${\bf k}_{\perp,i}$ for wave absorption, we calculate the imaginary part of wave frequency $\omega_i$, which is shown to be equivalent with the former technique under weak damping approximation. The code can use either numerical or analytical equilibrium. Examples and benchmarks with electron cyclotron wave, lower hybrid wave and ion cyclotron wave for tokamak, spherical tokamak (ST), FRC and mirror machine are shown.

Plasma Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is that in magnetically confined plasmas, most of the existing ray - tracing codes assume that the plasma density and temperature are functions of magnetic flux, and mainly study wave propagation within closed magnetic field lines. This is sufficient for use in current tokamak devices. However, these assumptions are difficult to apply to plasma configurations with open magnetic field lines (such as mirror machines and field - reversed configurations FRC). Therefore, the author has developed a new ray - tracing code, BORAY, which can support closed and open magnetic field line plasmas in any axisymmetric configuration in cylindrical coordinate systems, thus overcoming the limitations of existing codes for such configurations.

Specifically, the BORAY code aims to:

1. **Support open magnetic field line plasmas**: Unlike most existing codes, BORAY can handle plasma configurations with open magnetic field lines, such as mirror machines and field - reversed configurations (FRC).
2. **Provide broader applicability**: By not relying on the assumption that the plasma density and temperature must be functions of magnetic flux, BORAY can be applied to more types of plasma configurations, including tokamaks, spherical tokamaks, FRCs, and mirror machines.
3. **Improve the calculation method of wave absorption**: BORAY determines wave absorption by calculating the imaginary part of the wave frequency. This method is equivalent to calculating the imaginary part of the wave vector under the weak - damping approximation, but provides more accurate results.

In summary, the main objective of this paper is to develop a more general and accurate ray - tracing code to support the study of electromagnetic wave propagation and absorption in different plasma configurations, especially those with open magnetic field lines. This will be helpful for designing and optimizing plasma heating systems, selecting appropriate wave parameters, and improving plasma heating efficiency.