Spectral broadening from turbulence in multiscale lower hybrid current drive simulations

The scattering of lower hybrid (LH) waves due to scrape-off layer (SOL) filaments is investigated. It is revealed that scattering can account for the LH spectral gap without any ad hoc modification to the wave-spectrum. This is shown using a multiscale simulation approach which allows, for the first time, the inclusion of full-wave scattering physics in ray-tracing/Fokker-Planck calculations. In this approach, full-wave scattering probabilities are calculated for a wave interacting with a statistical ensemble of filaments. These probabilities are coupled to ray-tracing equations using radiative transfer (RT) theory. This allows the modeling of scattering along the entire ray-trajectory, which can be important in the multi-pass regime. Simulations are conducted for lower hybrid current drive (LHCD) in Alcator C-Mod, resulting in excellent agreement with experimental current and hard X-ray (HXR) profiles. A region in filament parameter space is identified in which the impact of scattering on LHCD is saturated. Such a state coincides with experimental LHCD measurements, suggesting saturation indeed occurs in C-Mod, and therefore the exact statistical properties of the filaments are not important.Comment: 23 pages, 13 figures. Revised to fix reference formattin

Application of linear electron Bernstein current drive models in reactor-relevant spherical tokamaks

Electron Bernstein current drive (EBCD) systems in spherical tokamaks are sensitive to plasma and launch conditions, and therefore require large parametric scans to optimise their design. One particular bottleneck in the simulation workflow is quasilinear modelling of current drive efficiency. Linear adjoint models are an attractive alternative, offering a ∼103 × speed-up compared to quasilinear codes. While linear models are well-tested and commonly used for electron cyclotron current drive (ECCD), they have seen little use in EBCD modelling. In this work, variants of the linear model are applied to EBCD and compared to quasilinear results in a reactor-relevant plasma, i.e. Spherical Tokamak for Energy Production (STEP). This comparison reveals it is important to accurately model the collision operator and finite Larmor radius effects in the linear model. When done properly, good agreement is found with quasilinear calculations, at least for normalised minor radii ρ < 0.7 and at low power densities. The power density threshold for quasilinear effects during EBCD is found to be significantly lower than that of ECCD. This is attributed to the much lower group velocity of the electron Bernstein wave (EBW). Thus, the linear model is only valid for EBCD modelling at low power densities (e.g. ≲ 1 MW launched EBW power in STEP). This may be satisfied in present-day experimental devices, but certainly not in reactors targeting non-inductive operation

A hybrid full-wave Markov chain approach to calculating radio-frequency wave scattering from scrape-off layer filaments

The interaction of radio-frequency (RF) waves with edge turbulence modifies the incident wave spectrum, and can significantly affect RF heating and current drive in tokamaks. Previous lower hybrid (LH) scattering models have either used the weak-turbulence approximation, or treated more realistic, filamentary turbulence in the ray tracing limit. In this work, a new model is introduced which retains full-wave effects of RF scattering in filamentary turbulence. First, a Mie-scattering technique models the interaction of an incident wave with a single Gaussian filament. Next, an effective differential scattering width is derived for a statistical ensemble of filaments. Lastly, a Markov chain solves for the transmitted wave spectrum in slab geometry. This model is applied to LH launching for current drive. The resulting wave spectrum is asymmetrically broadened in angular wavenumber space. This asymmetry is not accounted for in previous LH scattering models. The modified wave spectrum is coupled to a ray tracing/Fokker–Planck solver (GENRAY/CQL3D) to study its impact on current drive. The resulting current profile is greatly altered, and there is significant increase in the on-axis current and decrease in the off-axis peaks. This is attributed to a portion of the modified wave spectrum that is strongly dampened on-axis during the first pass.</jats:p