Convective radial energy flux due to resonant magnetic perturbations and magnetic curvature at the tokamak plasma edge
With the resonant magnetic perturbations (RMPs) consolidating as an important tool to control the transport barrier relaxation, the mechanism on how they work is still a subject to be clearly understood. In this work we investigate the equilibrium states in the presence of RMPs for a reduced MHD model using 3D electromagnetic fluid numerical code (EMEDGE3D) with a single harmonic RMP (single magnetic island chain) and multiple harmonics RMPs in cylindrical and toroidal geometry. Two different equilibrium states were found in the presence of the RMPs with different characteristics for each of the geometries used. For the cylindrical geometry in the presence of a single RMP, the equilibrium state is characterized by a strong convective radial thermal flux and the generation of a mean poloidal velocity shear. In contrast, for toroidal geometry the thermal flux is dominated by the magnetic flutter. For multiple RMPs, the high amplitude of the convective flux and poloidal rotation are basically the same in cylindrical geometry, but in toroidal geometry the convective thermal flux and the poloidal rotation appear only with the islands overlapping of the linear coupling between neighbouring poloidal wavenumbers , , .
The key to improve plasma confinement in tokamaks at high temperature and density lies on the transport barriers achieved in so called H-mode regimes, first discovered on ASDEX Wagner_PRL82 () in the early 1980’s. This high energy confinement mode (H-mode) with a steep pressure gradient at the edge represents out of few percent of the plasma radius in toroidal magnetic systems. In this regime the plasma pressure drops sharply over a narrow layer in the middle of the plasma edge. Once the H-mode is settled, this pressure gradient tends to increase with time until the peeling-ballooning instability leads to a turbulent state triggering the Edge Localized Modes (ELMs), which release the transport barriers in a time shot compared to the the energy confinement time . Connor_PoP98 (); Groeb_NF09 (). ELMs are periodic fast bursts of hot dense plasma on a fast time scale (s) and low amplitude, followed by expulsion of edge plasma and MHD activity, leading to confinement degradation Connor_PPCF98 (). Along with the fast barrier relaxation event, a turbulent transport through the barrier strongly increases and the pressure gradient drops down reaching a new configuration and then the barrier builds up again on a slow collisional time scale.
Turbulence simulations of transport barrier relaxations at the tokamak plasma edge have revealed that the control of such relaxations by RMPs is attributed to a local erosion of the barrier LecB_PRL09 (); LecB_NF10 (); BeyS_PPCF11 (). This erosion at the resonance position is known to be caused - at least partly - by the enhancement of the radial heat flux in presence of the RMP due to the strong parallel heat flow along perturbed field lines as is the so called magnetic flutter flux Fitz_PoP95 (). However, in the presence of magnetic curvature, an additional transport mechanism exists and is linked to stationary convection cells associated with the magnetic islands induced by RMPsLecB_PRL09 (); LecB_NF10 (). In presence of a mean poloidal velocity shear, this additional convective transport can be considerably high and even larger than the thermal flux from the magnetic flutter fluxBeyS_PPCF11 (). This previous result has been obtained in the frame of an electrostatic model.
We show here by numerical simulations in the frame of a 3D electromagnetic fluid turbulence model (EMEDGE3D)FuhB_PRL08 () and in the basic situation without turbulent fluctuations and without imposed mean velocity shear that two different equilibrium plasma states exists in presence of a RMP. The first regime is characterized by the absence of mean poloidal flow and a low level of convective transport. The second regime shows mean poloidal rotation and large convective transport. Here we show how these two equilibria depend on magnetic curvature and RMP amplitude. Interestingly, in the case of cylindrical curvature, the simple equilibrium without mean poloidal rotation and without considerable convective transport is found to be unstable such that the plasma evolves self-consistently to the poloidally rotating state where large convective transport is present. In the case of toroidal curvature, the simple equilibrium is found to be stable. A detailed stability analysis is presented.
In section II, the dimensionless 3D model of partial differential equations is presented with the profiles used for the code EMEDGE3D. In section III, the helical equilibrium states of the plasma driven by the RMP coils is derived for the cylindrical case and compared with simulations, for then analyse the results obtained in terms of the islands width in section IV. We show the different equilibrium states obtained with the toroidal geometry with multiple harmonics RMPs in section V, then we detach the main results of the work in section VI as a conclusion.
Ii 3D Plasma-RMP Model
Following BeyerBeyB_PPCF07 () and FuhrFuhB_PRL08 (), we introduce three reduced normalized magneto-hydrodynamical PDEs with source for the strength of the RMP currents. The equations are used to study the spatio-temporal behaviour of the three-dimensional fields of plasma pressure , electrostatic potential and magnetic flux .
In toroidal coordinates and in a slab geometry in the vicinity of a reference surface at the plasma edge, i.e. , , , the normalized operators are
Here, is the safety factor at the reference surface, is the major radius of the magnetic axis and is the magnetic shear length. Lengths parallel () and perpendicular () to the unperturbed magnetic field are normalized by and , respectively, and time is normalized by , where the resistive ballooning radial correlation length and the interchange time are given by
where is the ratio of the electron to the ion mass, and , , , are reference values of the electron collision time, the sound speed, the pressure gradient length, and the ion Larmor radius at electron temperature, respectively. For a collisional tokamak plasma edge, one typically finds and .
The dimensionless ballooning pressure gradient coefficient is similar to the normalized pressure gradient widely used in MHD tokamak stability theory , where is the ratio of plasma pressure to magnetic pressure. The parallel and perpendicular heat conductivity coefficients and are normalized by and , respectively. Therefore, a ratio of the normalized coefficients of , corresponds to a ratio of the dimensional coefficients of . In the present simulations, we use , , and the curvature parameter is set to .
The main computational domain corresponds to the volume delimited by the toroidal surfaces characterized by and , respectively, and including the reference surface . Here, a linear profile is assumed, and , for the reference parameters. The complete computational domain is slightly larger and delimited by and . The source is located in the inner buffer zone and gives rise to a constant incoming energy flux, from the plasma center into the main computational domain. The helical current
is located in the outer buffer zone (Fig. 2a). Here, the RMP coil structure is a function of the helical angle
Here we choose the regime RMP coil with , which the external current induces a magnetic perturbation resonant at as shown in Fig. 2b.
Iii Equilibrium states in presence of RMP
In the following, for a given RMP coil, two different helical equilibrium states of the plasma are obtained by the numerical calculations we first deduce analytically their main properties in the case of cylindrical curvature. In this case, there is no linear coupling between neighbouring poloidal wavenumbers , , , etc., and the equilibrium can be assumed to be of the form
where the bar and the index ”1” designate respectively the axisymmetric and the components, in which
Here, and . Eqs. (8)–(12) admit a symmetric solution corresponding to a helical equilibrium where the induced magnetic perturbation and the associated pressure perturbation are in phase with the external current (4), i.e. , , and the potential variation is in phase quadrature, . In this equilibrium, the poloidal rotation is zero, . In summary, the solutions of Eq. (8)–(12) reduces to
where , and are real fields corresponding to the imaginary (I) and real (R) parts of , and , respectively. Typical radial profiles of , , and are shown in Fig. 3. Note that there is a non-vanishing convective flux associated with the equilibrium (13), as the pressure and potential perturbations are phase shifted by . However, in typical situations as the one illustrated in Fig. 3, this convective flux is found to be much smaller than the magnetic flutter flux .
Stable helical equilibrium states of the plasma in presence of the RMP can be calculated numerically by integrating Eq. (1)–(3) in time starting with low level noise for all fields, providing that the pressure gradient stays below the resistive ballooning instability limit. This is guaranteed here by choosing a sufficiently low value of the total energy flux . However, with this convergence method, the symmetric equilibrium state (13) can only be obtained when explicitly forcing no rotation. This means that the symmetric equilibrium is unstable. This property will be discussed in the next section.
When the temporal evolution is calculated self-consistently including poloidal rotation, the system evolves to a new helical stationary state with non vanishing sheared plasma rotation. We call this new equilibrium here rotating state and it is characterized by: A. a phase shift between the external current and the magnetic perturbation induced in the plasma and B. a large convective flux . Typical radial profiles of the axisymmetric components as well as the helical amplitudes and the fluxes are shown in Fig. 4. Note that the rotation velocity vanishes close to the resonant surface which is consistent with a nearly complete penetration of the magnetic perturbation MonF_NF14 (), i.e. the amplitudes of the at the resonance surface are similar in both equilibria and close to the maximum value corresponding to the vacuum case. We may summarise the fields in the rotating steady state as
Note also that due to the strong convective flux, the flattening of the pressure profile on the resonant surface is significantly more pronounced in the rotating state compared with the symmetric state (cf. Figs. 3a and 4a). Two dimensional maps of the pressure and the non-axisymmetric part of potential corresponding to the two different equilibrium states are shown in Fig. 5. The periodic structure that corresponds to the flattening of the pressure on the magnetic islands is symmetric in the symmetric state (Fig. 5a) and distorted in the rotation state (Fig. 5b) from the shear flow. The potential structure together with the pressure structure, as illustrated in Fig. 5d, is responsible for the strong convective flux in the rotation state. Note that with increasing amplitude of the external RMP current, both the magnetic flutter flux and the convective flux increase in the stable rotating equilibrium, but the convective flux is always larger than the thermal flux produced by magnetic flutters as shown in Fig. 6.
Iv Transition to the rotating state with strong convective flux
The instability of the symmetric equilibrium and the transition to the second helical state described above can be illustrated by performing a time integration in two successive phases and following the evolution of the convective and magnetic flutter fluxes at the resonant surface (Fig. 7a). In an early phase of the integration (from on), the rotation is forced to zero and the system is rapidly evolving to the stationary state obtained in Eq. (13), where the convective flux is much smaller than the flutter flux . Then, from on, we release the constraint on the poloidal rotation and the system evolves self-consistently to the rotating state characterized by .
The instability of the symmetric equilibrium can be characterized by the growth rate of a small perturbations of this equilibrium. The growth rate can be determined in the numerical experiment described above. A series of such numerical simulations for different amplitudes of the external magnetic perturbation reveals that the growth rate is nearly constant up to a critical value of the RMP perturbation level. Above this level, the growth rate is strongly increasing with the external perturbation amplitude. This is illustrated in Fig. 7b where the growth rate is plotted against the magnetic island width linked to the perturbation amplitude via
The time growth rate of strongly increases for island widths . In that case, the half-width of the island approaches the distance between the resonant surface and the external boundary of the main computational domain and the island likely is influenced by the boundary. Also, higher order harmonics become significant for , where the critical island with is given by Fitz_PoP95 () . As illustrated in Fig. 8, for the amplitude of the second order mode is one order of magnitude lower than the amplitude of the main harmonics but for , the second harmonic is only lower by a factor of .
V Stable rotating state in toroidal geometry with multiple RMPs
Toroidal curvature gives rise to linear coupling between , , and modes. In particular, in the vorticity equation (1) a harmonic of the pressure couples to and harmonics of the electrostatic potential . If a self-consistent state with rotation and convective transport similar to the one shown in the Fig. 4 exists in toroidal geometry, then the state involve multiple harmonics. We therefore induce a multiple harmonics RMP with poloidal wavenumbers and toroidal wavenumber . In order to compare with the results shown above, we first apply the multiple harmonics RMP in the cylindrical curvature case. As illustrated in Fig. 9, a stable equilibrium with rotation and important convective flux is recovered in this case. The plasma self-organizes such that the rotation velocity vanishes close to the resonant surfaces . The convective and magnetic flutter fluxes show local maxima close to these resonant surfaces (Fig. 9b). The width of the island is and the width of the closest neighbouring island is , so the corresponding Chirikov overlapping parameter Chirikov () is
The islands therefore are overlapping, but this overlapping is weak enough such that distinct local maxima are observed for the and fluxes, as shown in Fig. 9b.
When applying a multiple harmonics RMPs in the toroidal curvature case, the plasma also evolves to a stable equilibrium with rotation and strong convective flux. This is illustrated in Fig. 11, which may be compared with Fig. 9. Differently from the cylindrical case in which a poloidal rotation is induced for all values of RMPs amplitudes, the poloidal rotation in toroidal geometry is triggered only for multiple harmonics RMPs and when Chirikov overlapping parameter , which correspond to values of external resonant magnetic perturbation with island width . In Fig. 10 we see for the three first values of , stays practically constant as the value of increases. This corresponds to the situation where the Chirikov overlapping parameter is smaller than one.
Comparing figure 4, 9 and 11 we see that even for different shapes and intensities of the equilibrium electric potential and convective fluxes, the poloidal rotation in all the three figures are at the same order of magnitude and have the maximum shear at the resonant surface. This indicates that an important induced poloidal rotation can be controlled by the RMP currents.
In this work we investigate the dependence of the radial transport of the thermal energy from convection and magnetic flutter from the plasma response to a reference model of the resonant magnetic perturbations RMPs. The stationary states are studied running the 3D plasma edge turbulence code EMEDGE3D below to primary ballooning instability threshold. The simple static equilibrium in which the magnetic perturbation inside the plasma is in phase with the external perturbation and the plasma is not rotating is found to be unstable. The plasma is self-organized into a more complex state where the perturbation becomes phase shifted and the plasma rotates. This is due to the coupling between pressure and electrostatic potential perturbations induced by the magnetic curvature. In the stable equilibrium state, the phase between pressure and electrostatic potential is such that produces a rotating plasma with a significant thermal convective that exceeds the thermal flux from the magnetic flutter by a factor of 3–10 times. The instability of the simple equilibrium and the subsequent evolution to a new stable complex equilibrium has been first investigated in cylindrical geometry for a single harmonic RMP perturbation and latter for multiple RMP modes. We showed that for both situations the system reaches the same final equilibrium state with high convective flux and induced poloidal rotation. Then, we showed that the corresponding coupling mechanism between pressure and the electric potential also produces plasma rotation and a strong convective flux in the toroidal geometry with multiple RMPs. The convective thermal flux and the magnetic flutter flux are about the same order of magnitude, but the presence of the thermal convective flux condition is only achieved when Chirikov overlapping parameter is greater than one.
Acknowledgements.The authors gratefully acknowledge Wendell Horton and Iberê Caldas for useful discussions. This research was supported by the French National Research Agency, project ANR- 2010-BLAN-940-01, and Brazilian national research agencies CNPq/CAPES under project No. 163395/2013-6 and FAPESP No. 201119296-1.
- (1) F. Wagner, G. Becker, K. Behringer, D. Campbell, A. Eberhagen, W. Engelhardt, G. Fussmann, O. Gehre, J. Gernhardt, G. v. Gierke, G. Haas, M. Huang, F. Karger, M. Keilhacker, O. Klüber, M. Kornherr, K. Lackner, G. Lisitano, G. G. Lister, H. M. Mayer, D. Meisel, E. R. Müller, H. Murmann, H. Niedermeyer, W. Poschenrieder, H. Rapp, H. Röhr, F. Schneider, G. Siller, E. Speth, A. Stäbler, K. H. Steuer, G. Venus, O. Vollmer, and Z. Yü, Phys. Rev. Lett. 49, 1408 (Nov 1982), http://link.aps.org/doi/10.1103/PhysRevLett.49.1408.
- (2) J. W. Connor, R. J. Hastie, H. R. Wilson, and R. L. Miller, Physics of Plasmas 5, 2687 (1998), http://scitation.aip.org/content/aip/journal/pop/5/7/10.1063/1.872956% BibitemShutNoStop
- (3) R. J. Groebner, T. Osborne, A. Leonard, and M. Fenstermacher, Nuclear Fusion 49, 045013 (2009).
- (4) J. W. Connor, Plasma Physics and Controlled Fusion 40, 531 (1998), http://stacks.iop.org/0741-3335/40/i=5/a=002.
- (5) M. Leconte, P. Beyer, X. Garbet, and S. Benkadda, Physical Review Letters 102, 045006 (2009), http://link.aps.org/doi/10.1103/PhysRevLett.102.045006.
- (6) M. Leconte, P. Beyer, X. Garbet, and S. Benkadda, Nuclear Fusion 50, 054008 (2010).
- (7) P. Beyer, F. de Solminihac, M. Leconte, X. Garbet, F. L. Waelbroeck, A. I. Smolyakov, and S. Benkadda, Plasma Physics and Controlled Fusion 53, 054003 (2011).
- (8) R. Fitzpatrick, Physics of Plasmas 2, 825 (1995).
- (9) G. Fuhr, P. Beyer, and S. Benkadda, Physical Review Letters 101, 195001 (2008).
- (10) P. Beyer, S. Benkadda, G. Fuhr-Chaudier, X. Garbet, P. Ghendrih, and Y. Sarazin, Plasma Physics and Controlled Fusion 49, 507 (2007).
- (11) P. B. A. Monnier, G. Fuhr, Nuclear Fusion Accepted for publication, 000 (2014).
- (12) B. Chirikov, Physics Reports 52, 263 (1979).