14.5 Safety factor profile

The magnetic shear for a concentric-circular configuration is defined by

ˆs = dq r,
   dr q
(385)

where r is the minor radius of a magnetic surface. The above expression can be re-arranged as

 dr   dq
ˆsr-=  q-
(386)

Integrating the above equation over r and assuming ŝ is a constant, we obtain

 ∫       ∫
   r dr    r dq
ˆs r0 r =  r0 q ,
(387)

Performing the integration, the above equation is written as

ˆs(ln r− ln r0) = ln q− ln q0,
(388)

where q0 = q(r0). Equation (388) can be finally written as

     ( r )ˆs
q = q0 r-  .
        0
(389)

This is a profile with a constant magnetic shear s. In Ben’s toroidal ITG simulation, the following q profile is used:

               ′
q = q0 + (r − r0)q(r0),
(390)

with q(r0) = ŝq0∕r0. This is a linear profile over r, with the values of q and the shear at r = r0 being q0 and ŝ, respectively.