Full-wave electromagnetic field Simulations of lower hybrid waves in tokamaks

The most common method for treating wave propagation in tokamaks in the lower hybrid range of frequencies (LHRF) has been toroidal ray tracing, owing to the short wavelengths (relative to the system size) found in this regime. Although this technique provides an accurate description of 2D and 3D plasma inhomogeneity effects on wave propagation, the approach neglects important effects related to focusing, diffraction, and finite extent of the RF launcher. Also, the method breaks down at plasma cutoffs and caustics. Recent adaptation of full-wave electromagnetic field solvers to massively parallel computers [1] has made it possible to accurately resolve wave phenomena in the LHRF. One such solver, the TORIC code, has been modified to simulate LH waves by implementing boundary conditions appropriate for coupling the fast electromagnetic and the slow electrostatic waves in the LHRF. In this frequency regime the plasma conductivity operator can be formulated in the limits of unmagnetized ions and strongly magnetized electrons, resulting in a relatively simple and explicit form. Simulations have been done for parameters typical of the planned LHRF experiments on Alcator C-Mod, demonstrating fully resolved fast and slow LH wave fields using a Maxwellian non-relativistic plasma dielectric. Significant spectral broadening of the injected wave spectrum and focusing of the wave fields have been found, especially at caustic surfaces. Comparisons with toroidal ray tracing have also been done and differences between the approaches have been found, especially for cases where wave caustics form. The possible role of this diffraction-induced spectral broadening [2] in filling the spectral gap in LH heating and current drive will be discussed.