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,ϕ,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 B(r,z) and species densities ns0(r,z) profiles as input. For wave absorption, we require a further temperature Ts0(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 k⊥,i for wave absorption, we calculate the imaginary part of wave frequency ω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.Program summaryProgram Title: BORAYCPC Library link to program files: https://doi.org/10.17632/tnkrjdbcz8.1Code Ocean capsule: https://codeocean.com/capsule/6205646Licensing provisions: BSD 3-clauseProgramming language: MatlabNature of problem: Solve the plasmas electromagnetic wave propagation and absorption in the geometrical optics approximation for magnetized plasmas based on ray tracing of plasma dispersion relation. In axisymmetric (r,z) coordinates, the code can be used for both closed and open field lines plasmas of various configurations such as tokamak, spherical tokamak, FRC and mirror machine.Solution method: Runge-Kutta time integral to solve ray tracing equations for wave propagation, and integral the imaginary part of the wave frequency in hot kinetic dispersion relation for wave absorption.Additional comments including restrictions and unusual features: Kinetic relativistic effects and collisional damping are not included in the present version yet. Only axisymmetric two-dimensional (2D) profiles are support in present version.

Quasi-Two-Dimensional Code for the Calculation of Antenna Impedance of the ICR Heating System

3D simulations of the RF wave propagation and power coupling in cold magnetized plasma near an ICRF antenna

Design of an optimized load-resilient Conjugate T for the ICRH system in the LHD using a novel hybrid circuit/3DLHDAP code and experimental results

Symplectic Simulations of Radial Diffusion of Fast Alpha Particles in the Presence of Low-Frequency Modes in Rippled Tokamaks

Modelling of Three-Ion ICRF Heating Scenario for Tokamak Globus-M2

Development of Impedance Sheath Boundary Condition in Stix Finite Element RF Code

Design and operations of a load-tolerant ICRH system in Experimental Advanced Superconducting Tokamak

Ion Cyclotron Resonance Heating (icrh) Systems for the Keda Mirror with Axisymmetry (Kmax)

Theoretical Study of the Impedance Matching of the Ferrite Tuner System for ICRF Heating

Preliminary Conceptual Design of the ICRH Antenna for DTT: A Systems Engineering Approach

A full wave solver[1] integrated with a Fokker–Planck code for optimizing ion heating with ICRF waves for the ITER deuterium-tritium plasma

ANSYS HFSS as a new numerical tool to study wave propagation inside anisotropic magnetized plasmas in the Ion Cylotron Range of Frequencies

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Integral electrical characteristics and local plasma parameters of a RF ion thruster

Ion Radial Transport Induced by ICRF Waves in Tokamaks

BORAY: A ray tracing code for various magnetized plasma configurations

LCODE: Quasistatic code for simulating long-term evolution of three-dimensional plasma wakefields

Impedance and voltage power spectra of a monopole antenna in a warm plasma -- derivation and application to CubeSats

Predicting high harmonic ion cyclotron heating efficiency in Tokamak plasmas

Reinterpreting radio frequency heating and current drive theory in a tokamak