The current response of a hot magnetized plasma to a radio-frequency wave is non-local, turning the electromagnetic wave equation into an integro-differential equation. Non-local physics gives rise to wave physics and absorption processes not observed in local media. Furthermore, non-local physics alters wave propagation and absorption properties of the plasma. In this work, an iterative method that accounts for parallel non-local effects in 2D axisymmetric tokamak plasmas is developed, implemented, and verified. The iterative method is based on the finite element method and Fourier decomposition, with the advantage that this numerical scheme can describe non-local effects while using a high-fidelity antenna and wall representation, as well as limiting memory usage. The proposed method is implemented in the existing full wave solver FEMIC and applied to a minority heating scenario in ITER to quantify how parallel non-local physics affect wave propagation and dissipation in the ion cyclotron range of frequencies (ICRF). The effects are then compared to a reduced local plane wave model, both verifying the physics implemented in the model, as well as estimating how well a local plane wave approximation performs in scenarios with high single pass damping. Finally, the new version of FEMIC is benchmarked against the ICRF code TORIC.
This paper presents an iterative method to account for parallel non-local effects in 2D axisymmetric tokamak plasmas during radio-frequency wave propagation and absorption. The method is implemented in the FEMIC code and applied to study minority heating in ITER, comparing the results to a local plane wave model. The new approach can describe non-local physics while using a high-fidelity representation of the antenna and wall, and requires less memory.