No Arabic abstract
The neoclassical prescription to use an equation of motion to determine the electrostatic field within a tokamak plasma is fraught with difficulties. Herein we examine two popular expressions for the equilibrium electrostatic field so determined and show that one fails to withstand a formal scrutiny thereof while the other fails to respect the vector nature of the diamagnetic current. Reconsideration of the justification for the presence of the equilibrium electrostatic field indicates that no field is needed for a neutral plasma when considering the net bound current defined as the curl of the magnetization. With any shift in the toroidal magnetic flux distribution, a dynamic electric field is generated with both radial and poloidal components, providing an alternate explanation for any measurements thereof.
From a common expression for the poloidal electrostatic field of a tokamak, in the limit of large aspect ratio and concentric circular flux surfaces, one may determine the associated potential. This potential satisfies Poissons equation, which reduces to Laplaces equation when the medium has vanishing charge density, in axial geometry but not toroidal geometry. A simple transformation takes the potential over to the correct harmonic form for tokamak coordinates, and the resulting electrostatic field is calculated. From the radial field one may estimate the supporting charge density on the boundary, and from the poloidal field one may determine a prediction for the radial dependence of the electron temperature, which does not compare well with a rough estimate of the profile often seen in a tokamak.
Various MHD (magnetohydrodynamic) equilibrium tools, some of which being recently developed or considerably updated, are used on the COMPASS tokamak at IPP Prague. MHD equilibrium is a fundamental property of the tokamak plasma, whose knowledge is required for many diagnostics and modelling tools. Proper benchmarking and validation of equilibrium tools is thus key for interpreting and planning tokamak experiments. We present here benchmarks and comparisons to experimental data of the EFIT++ reconstruction code [L.C. Appel et al., EPS 2006, P2.184], the free-boundary equilibrium code FREEBIE [J.-F. Artaud, S.H. Kim, EPS 2012, P4.023], and a rapid plasma boundary reconstruction code VacTH [B. Faugeras et al., PPCF 56, 114010 (2014)]. We demonstrate that FREEBIE can calculate the equilibrium and corresponding poloidal field (PF) coils currents consistently with EFIT++ reconstructions from experimental data. Both EFIT++ and VacTH can reconstruct equilibria generated by FREEBIE from synthetic, optionally noisy diagnostic data. Hence, VacTH is suitable for real-time control. Optimum reconstruction parameters are estimated.
Uncertainties and errors in magnetic equilibrium reconstructions are a wide-spread problem in interpreting experimental data measured in the tokamak edge. This study demonstrates errors in EFIT++ reconstructions performed on the COMPASS tokamak by comparing the outer midplane separatrix position to the Velocity Shear Layer (VSL) position. The VSL is detected as the plasma potential peak measured by a reciprocating ball-pen probe. A subsequent statistical analysis of nearly 400 discharges shows a strong systematic trend in the reconstructed separatrix position relative to the VSL, where the primary factors are plasma triangularity and the magnetic axis radial position. This dependency is significantly reduced after the measuring coils positions as recorded in EFIT input are optimised to provide a closer match between the synthetic coil signal calculated by the Biot-Savart law in a vacuum discharge and the actual coil signal. In conclusion, we suggest that applying this optimisation may lead to more accurate and reliable reconstructions of the COMPASS equilibrium, which would have a positive impact on the accuracy of measurement analysis performed in the edge plasma.
A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and large deviations from thermal equilibrium. The formulation of the gyrokinetic model is based on a second order description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts, assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behavior of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function and includes a non-linear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a velocity-space Hermite-Laguerre polynomial basis and the gyrokinetic equation is projected onto the same basis, obtaining the spatial and temporal evolution of the Hermite-Laguerre expansion coefficients. The Hermite-Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwells equations derived from a variational principle, with the velocity integrals of the gyroaveraged distribution function explicitly evaluated.
Axisymmetric free-boundary equilibrium reconstructions of tokamak plasmas in the Lithium Tokamak eXperiment (LTX) are performed using the PSI-Tri equilibrium code. Reconstructions in LTX are complicated by the presence of long-lived non-axisymmetric eddy currents generated by vacuum vessel and first wall structures. To account for this effect, reconstructions are performed with additional toroidal current sources in these conducting regions. The source distributions are fixed poloidally, but their scale is adjusted as part of the full reconstruction. Eddy distributions are computed by toroidally averaging currents, generated by coupling to vacuum field coils, from a simplified 3D filament model of important conducting structures. The full 3D eddy current fields are also used to enable the inclusion of local magnetic field measurements, which have strong 3D eddy current pick-up, as reconstruction constraints. Using this method, equilibrium reconstruction yields good agreement with all available diagnostic signals. An accompanying field perturbation produced by 3D eddy currents on the plasma surface with primarily n=2, m=1 character is also predicted for these equilibria.