Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
10-15 May 2021
Nice, France
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Propagation of radio frequency waves through turbulent plasmas

11 May 2021, 14:00
4h 45m
Nice, France

Nice, France

Regular Poster Magnetic Fusion Theory and Modelling P2 Posters 2


Dr Abhay Ram (Massachusetts Institute of Technology)


The practical and economic viability of tokamak fusion reactors depends, in a significant way, on the efficiency of radio frequency (RF) waves to deliver energy and momentum to the plasma in the core of the reactor. The RF electromagnetic waves, excited by antenna structures placed near the wall of a tokamak, have to propagate through the turbulent edge plasma along their path to the core of the fusion device. In present day experiments, the radial width of the edge region and scrape-off layer is of the order of a few centimeters. In ITER, and in future fusion reactors, this width will be of the order of tens of centimeters. Any modifications to the properties of RF waves in the edge region can have deleterious effects on the efficiency and spatial distribution of power coupled to the plasma. This could affect, for example, the stabilization of the neoclassical tearing mode in ITER by electron cyclotron RF waves. Thus, it is imperative to properly understand the propagation of RF waves through plasma turbulence. This paper is on a multi-pronged, theoretical and computational, approach that is being pursued to quantify the effect of edge plasma turbulence on the propagation of RF waves. The theoretical and analytical models are based on solutions of the Faraday-Ampere equation in a magnetized plasma (1-3) and on the Kirchhoff tangent plane approximation (4). An effective medium approach has been developed so as to approximate the permittivity of a turbulent plasma by analytical expressions (5). The computations are being carried out with a newly developed code ScaRF that solves Maxwell's equations in three-dimensions for an arbitrarily assigned plasma permittivity (6).


A full-wave theoretical model, using Maxwell's equations, has been developed for scattering of RF plane waves by blobs and filaments embedded in a magnetized plasma (1-3). The spatial structure of the electric and magnetic fields is given by the Faraday-Ampere equation, and mathematical techniques developed for Mie scattering are used to solve the wave equations for a cold plasma permittivity. We have formulated a complimentary analytical model based on the Kirchhoff's tangent plane approximation [4]. The Kirchhoff method has commonly been used for studying scattering of electromagnetic waves from rough surfaces. The basis of the Kirchhoff approximation is built upon the theory of physical optics. In this approximation, the scattered field due to planar turbulence is determined by the wave fields on the surface of the turbulence separating two different plasma densities. Each point on this surface is assumed to be part of an infinitely extended plane that is along the local tangent at that point. Consequently, the Kirchhoff approximation, along with physical optics, leads to the study of reflection, refraction, and side-scattering of RF waves.This not only simplifies the modeling of RF scattering but, more importantly, gives physical insight into several, experimentally significant, aspects of scattering. We find that the spectrum of the RF fields is modified along, and across, the direction of the magnetic field. We have derived analytical expressions for the changes in the wave spectrum based on the spatial variation of the density fluctuations. The theory also predicts linear coupling of RF waves due to turbulence. For example, in the electron cyclotron range of frequencies, an incident ordinary wave can couple some of its power to the extraordinary wave in the presence of fluctuations. The theoretical results are being validated by full-wave simulations of RF scattering by different forms of plasma turbulence.


A full-wave electromagnetic computational code ScaRF, based on the finite difference frequency domain (FDFD) method, has been developed to study the effects of density turbulence on the propagation of RF waves (6). The anisotropic plasma permittivity used in the code is that for a magnetized, cold plasma. The code was initially used to study the propagation of an RF plane wave through a spatially modulated, periodic, density interface. Such an interface arises in the edge region due to magnetohydrodynamic instability and/or drift waves. Figure 1 shows the effect of a sinusoidal density perturbation on the propagation of a plane wave (incident from the top of the figure). Besides the transmitted wave through the density perturbation, there is a reflected wave as well. The Fourier spectrum of the transmitted RF electric field, Fig. 2, shows a broadening of the spectrum as the wave propagates away from the density fluctuation. The results from ScaRF are being used to validate the analytical results obtained from our theoretical models and for determining the limitations of the Kirchhoff approximation. While ScaRF has been used to study a periodic density fluctuation, the code is general enough to include different varieties of density fluctuations in the edge region - such as blobs and filaments, and spatially random fluctuations.
ScaRF simulation of scattering from a periodic fluctuation (black line). This shows the magnitude of the Poynting flux for an incident ordinary plane wave in the electron cyclotron range of frequencies. The magnetic field is in the z-direction and the wave is incident from the negative x-direction. The red lines and the arrows indicate the flow of RF power.

The Fourier spectrum of the x-component of the wave fields transmitted through the density fluctuation. This shows that the density perturbation can spread the RF power over a broader range of the spectrum along the magnetic field. The abscissa is the component of the wave vector normal to the magnetic field and pointing towards the core of the plasma. The incident wave has just a single value of the component defining its spectrum. Part of the broadening observed in the figure is due to the coupling of the incident ordinary wave with the transmitted extraordinary wave induced by the fluctuation.


The theoretical and computational studies on scattering, using either the full-wave approach or the physical optics model, require a proper description of the plasma permittivity in the turbulent edge region. The tenuous plasma in this region is a mix of fluctuations with a broad range of spatial correlations. Electromagnetic homogenization, or ``effective medium approximation'', is the process of estimating the effective electromagnetic properties of composite materials. We have generalized this approximation, following the Maxwell-Garnett approach to describe the plasma permittivity in the scrape-off layer (5). Our formulation is suitably adapted for magnetized plasmas and for RF waves that are commonly used for heating and current drive.The advantage of this approach is that we can describe the turbulent plasma as an effective homogeneous medium for which the permittivity is suitably approximated. This eliminates the need to have a detailed description of the edge plasma - an experimentally formidable task. The effective permittivity is being implemented in a computational code and will be validated against specific models of the plasma permittivity.


AKR is supported by the US Department of Energy Grant numbers DE-FG02-91ER-54109 and DE-FC02-01ER54648. All other authors are supported by EUROfusion and the National Program for Thermonuclear Fusion.


(1) A. K. Ram and K. Hizanidis, Phys. Plasmas 23, 022504 (2016).

(2) Z. Ioannidis, A. K. Ram, K. Hizanidis, and I. G. Tigelis, Phys. Plasmas 24, 102115 (2017).

(3) S.-I. Valvis et al., J. Plasma Phys. 84, 745840604 (2018).

(4) P. Beckmann and A. Spizzichino, ``The scattering of electromagnetic waves from rough surfaces,'' Artech House, Massachusetts (1987).

(5) F. Bairaktaris et al., J. Plasma Phys. 84, 45840501 (2018).

(6) A. D. Papadopoulos et al., J. Plasma Phys. 85, 905850309 (2019).

Country or International Organization United States
Affiliation Massachusetts Institute of Technology

Primary author

Dr Abhay Ram (Massachusetts Institute of Technology)


Prof. Kyriakos Hizanidis (National Technical University of Athens) Dr Fotis Bairaktaris (National Technical University of Athens) Dr Aristeides Papadopoulos (National Technical University of Athens) Dr S.-I. Valvis (National Technical University of Athens)

Presentation Materials

There are no materials yet.