Efficient lower hybrid current drive (LHCD) is demonstrated at densities up to nover ¯_{e}≈1.5×10^{20} m^{-3} in diverted plasmas on the Alcator C-Mod tokamak by operating at increased plasma ...current-and therefore reduced Greenwald density fraction. This density exceeds the nominal "LH density limit" at nover ¯_{e}≈1.0×10^{20} m^{-3} reported previously, above which an anomalous loss of current drive efficiency was observed. The recovery of current drive efficiency to a level consistent with engineering scalings is correlated with a reduction in density shoulders and turbulence levels in the far scrape-off layer. Concurrently, rf wave interaction with the edge and/or scrape-off-layer plasma is reduced, as indicated by a minimal broadening of the wave frequency spectrum measured at the plasma edge. These results have important implications for sustaining steady-state tokamak operation and indicate a pathway forward for implementing efficient LHCD in a reactor.
Abstract
This paper reports a new numerical scheme to simulate the radio-frequency (RF) induced RF sheath, which is suitable for a large 3D simulation. In the RF sheath boundary model, the tangential ...component of the electric field (
E
t
) is given by the gradient of a scalar electric field potential. We introduce two additional scalar potentials for the tangential components of the magnetic field, which effectively impose the normal electric displacement (
D
n
) on the plasma sheath boundary condition via in-homogeneous Neumann boundary condition and constrain the tangential electric field on the surface as curl-free (
∇
×
E
t
=
0
). In our approach, the non-linear sheath impedance is formulated as a natural extension of the large thickness (or asymptotic) sheath limit (
D
n
=
0
), allowing for handling both asymptotic and non-linear regimes seamlessly. The new scheme is implemented using the Petra-M finite element method analysis framework and is verified with simulations in the literature. The significance of non-linearity is discussed in various plasma conditions. An application of this scheme to asymptotic RF sheath simulation on the WEST ICRF antenna side limiters is also discussed.
A newly developed generic electro-magnetic (EM) simulation tool for modeling RF wave propagation in SOL plasmas is presented. The primary motivation of this development is to extend the domain ...partitioning approach for incorporating arbitrarily shaped SOL plasmas and antenna to the TORIC core ICRF solver, which was previously demonstrated in the 2D geometry S. Shiraiwa, et. al., “HISTORIC: extending core ICRF wave simulation to include realistic SOL plasmas”, Nucl. Fusion in press, to larger and more complicated simulations by including a 3D realistic antenna and integrating RF rectified sheath potential model. Such an extension requires a scalable high fidelity 3D edge plasma wave simulation. We used the MFEM http://mfem.org, open source scalable C++ finite element method library, and developed a Python wrapper for MFEM (PyMFEM), and then a radio frequency (RF) wave physics module in Python. This approach allows for building a physics layer rapidly, while separating the physics implementation being apart from the numerical FEM implementation. An interactive modeling interface was built on pScope S Shiraiwa, et. al. Fusion Eng. Des. 112, 835 to work with an RF simulation model in a complicated geometry.
A new technique for the calculation of RF waves in toroidal geometry enables the simultaneous incorporation of antenna geometry, plasma facing components (PFCs), the scrape off-layer (SOL) and core ...propagation Shiraiwa et al, “HIS-TORIC: Extending core ICRF wave simulation to include realistic SOL plasmas”, Nucl. Fusion in press 10 (2017). The calculation of antenna coupling has been traditionally performed separately from those of core propagation whereas this method allows them to be combined permitting the modelling of such phenomena as waves reflected from plasma-vacuum boundary or cyclotron resonances back to the antenna. Calculations with this technique also naturally capture wave propagation in the SOL and its interactions with non-conforming PFCs permitting self-consistent calculation of core absorption and edge power loss.
Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D ...simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ($n_{e0}$, $T_{e0}$, $I_p$, $B_t$, $R_0$, $n_{\|}$, $Z_{{\rm eff}}$, $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ms with high accuracy across the input parameter space.