First principles model of internal transport barriers in negative central shear discharges

At and around minima or maxima of thc safety factor q, overlapping, and therefore toroidicity induced coupling of drift eigenmodes centred on neighbouring rational surfaces r(l,m) defined by q(r(l,m))= m/l and r(l,m+/-1) defined by q(r(l,m+/-1)) = (m+/-1)/l, is generally negligible; l and m(m+/-1) are the toroidal and (main) poloidal mode numbers, respectively. Under these conditions, a careful analysis shows that the growth rate of the ion temperature gradient (ITC) mode is proportional to thc absolute value of the magnetic shear parameter; its radial width and Frequency (in the E x B rotating frame) arc, larger in the axisymmetric torus than in the plasma slab or cylinder, especially at small values of k(theta)a, (k(theta) is the poloidal mode number and a(i) the characteristic ion Larmor radius). These results provide a straightforward theoretical interpretation of the origin of internal ion transport barriers (ITBs) in discharges with negative central shear if the ITC instability is indeed the main channel for ion energy transport. In agreement with recent experiments, the mechanism does not degrade for T-e similar to T-i, an important result for extrapolation to reactors. Once initiated, an ITB will locally increase the curvature of the ion temperature profile and, consequently, the gradients of the radial electric field and E x B rotation profiles (until neoclassical (or sub-neoclassical) transport conics into play). The damping associated with velocity shear: yields an upper bound on the unstable mode number range, whereas Landau damping is expected to introduce a lower bound. In cases where toroidal effects arc subdominant, the ratio of the damping rate associated with velocity shear to the growth rate associated with magnetic shear leads naturally to the Hamaguchi-Horton parameter (divided by the factor eta (i) - 2/3 which characterizes the departure from the instability threshold (eta (i) is the ratio of the density and temperature length scales)) when introducing the most unfavourable mode number. Toroidal effects are always dominant for sufficiently weak magnetic shear or sufficiently large values of \L-N\/R (where L-N is the density gradient length), and therefore before steepening of the profiles has occurred near the zero magnetic shear layer; the counterpart of the Hamaguchi-Horton characteristic parameter appropriate to that case is also obtained, considering again the most unfavourable mode number.