Mitigation of ELMs by Electrostatic Field in Tokamaks

Mitigation of ELMs by electrostatic field is studied. The perpendicular heating in cyclotron waves tends to pile up the resonant particles toward the low magnetic field side in which a electrostatic field may result [J. Y. Hsu, V. S. Chan, R. W. Harvey, R. Prater, and S. K. Wong, Phys. Rev. Lett. 53, 564 (1984)]. The electrostatic field can make circulating particles trapped or make trapped particles circulating depending on the field direction. The trappedparticle population and bootstrap current change accordantly. Modulating bootstrap current, mitigation of type-1 ELM by the electrostatic field is possible. The electrostatic potential needed for the mitigation is quantitatively estimated. Experiments by either ECRH or biasing are being prepared to verify the theory.


Introduction
In present tokamaks operating in high-confinement regimes (H-modes), the steep pressure gradients at edge are often observed to relax through frequent intermittent discharges of energy, known as ELMs.The physics of ELMs is a key issue for ITER operation.The onset of ELMs constrains the pressure at top of edge transport barrier (pedestal height).The ELMs events transport sub-Journal of High Energy Physics, Gravitation and Cosmology stantial heat and particle loads to plasma-facing materials.A predictive understanding of the onset of type-I ELMs has been gained via the development of peeling-ballooning modes [1] in which EL Ms are triggered by instabilities driven by the large pressure gradient and bootstrap current in the edge.High pressure is important for fusion efficiency.The bootstrap current can be changed.
The perpendicular heating in cyclotron waves tends to pile up the resonant particles toward the low magnetic field side.An electrostatic field may result [2].
Variations of the electrostatic potential at plasma edge are observed in HL-2A [3].Full particle simulation is performed using the Boris algorithm [4].The electrostatic field can make circulating particles trapped or make trapped particles circulating depending on the field direction.With the assumption neoclassical transport the population of the trapped particles and bootstrap current change accordantly.Modulating bootstrao current by changing the electrostatic field, mitigation of type-1 ELM is possible.The electrostatic potential needed for the mitigation is quantitatively calculated.Experiments by either ECRH or biasing [5] are being prepared to verify the theory.

Full Particle Orbit Simulation in Tokamaks
In particle simulations of magnetized plasmas, the Boris algorithm [4] is the standard for advancing a charged particle in an electromagnetic field in accordance with the equation of motion associated with the Lorentz force, where the magnetic field and electric field are given respectively by where Ψ is the poloidal magnetic flux, where ( ) We use ITER's parameters: . So the tokamak magnetic field is well-determined.Full orbit simulations find electric trapping and de-trapping seen in Figure 1 and Figure 2 respectively.

Bootstrap Current
The gyro-averaged Hamiltonian has been given in Ref. [6], ( ) where the momenta p Rv e are conjugate to α , the gyrophase, φ, the toroidal angle, and x, expressed as where R and Z are the coordinates of the guiding center in a cylindrical system, ρ is the Larmor radius, Ω is the toroidal gyro-frequency.The particle mass is taken to be unity for simplicity.The electrostatic potential is assumed in a form, which is like the dipole potential produced by two close-point-charges, where ε is the inverse aspect ratio.
From Equation ( 6) we have where For the large aspect-ratio approximation we have, where The ions with are trapped, however, they are circulating without the electrostatic field.That is electrostatic trapping.For electrons the trapping condition is There is minimum of ( ) where p B is the poloidal magnetic field, n is the density, T is plasma temperature, n L is the density scale length.The gradients in the electron profiles con- tribute to typically 70% -90% of the total bootstrap current [9].

Peeling-Ballooning Modes
The criterion of peeling-ballooning modes can be expressed by the following formula [10], where m D is the Mercier coefficient, 1 4 m D < is the Mercier stability criterion, finite (positive) bootstrap current, j , is destabilizing and q′ is the de- rivative of the safety factor with respect to the poloidal magnetic flux.At pedestal the temperature is low, therefore, from Equation ( 16) bootstrap current is sensitive to the electrostatic potential.Now we use Equation ( 17) to calculate the criterion.For a large aspect ratio and low β ordering Equation ( 17) can be written [1]   || where  and e is the elongation [11].We neglect triangularity, δ , then Equation (17) becomes which can be produced in the practical experiments [5].
Electrostatic field, hopefully, can realize ELM-control like that in Ref. [12] and show synchronization of the ELM cycle with added electrostatic field.Electrostatic field, hopefully, can realize ELM-ree discharge which appears in I-mode of Alcator C-mod [13].

Summary
Full particle simulation is suitable to a multi-scale problem.The Boris algorithm makes long-time simulation accurate.The perpendicular heating in cyclotron waves tends to pile up the resonant particles toward the low magnetic field side in which electrostatic field may result [2].The electrostatic field can make circulating particles trapped or make trapped particles circulating depending on the field direction.The trapped-particle population and bootstrap current change accordingly in the process.Modulating bootstrap current, mitigation of type-1 ELM or ELM-free discharge is possible.Experiments by either ECRH or biasing [5] are being prepared to verify the theory in HL-2A Tokamak [14].

Figure 1 .
Figure 1.(a) The electrostatic field E is zero, particle is circulating; (b) The electrostatic field E is 35 kv/m, particle is trapped.Ion energy is 60 kev, pitch angle is 72˚.

Figure 2 .
Figure 2. (a) The electrostatic field E is zero, particle is trapped; (b) The electrostatic field E is −35 kv/m, particle is turned to be circulating which is called as de-trapping.Ion energy is 60 kev, pitch angle 151˚.
equilibrium distribution-function is Maxwellian it is easy to calculate trapped-electron population.The fraction of trapped electrons is