Island Width Evolution

Let us write the generalized Rutherford island width evolution equation, (13.7), in the form

$\displaystyle \tilde{\tau}_R\,\frac{d\hat{W}}{dt}= {\mit\Delta} - f_{sat}\,\hat...
...W}} -\frac{f_{curv}}{\hat{W}}-\frac{f_{polz}}{\hat{W}^3}
+F(\hat{W}_c,\hat{W}),$ (13.63)

where

$\displaystyle F(\hat{W}_c,\hat{W})=f_{polz1}\,\hat{W}_c^{\,2}\,\sin\varphi+ f_{...
...\,4}\,\hat{W}^3\,\sin^2\varphi + \frac{\hat{W}_c^{\,2}}{\hat{W}^2}\,\cos\varphi$ (13.64)

represents the three additional terms that appear on the right-hand side of the equation when an island chain interacts with a static resonant magnetic perturbation.

First, let us consider a rotating magnetic island chain. Given that the island width evolves on a much longer timescale than the island phase (see Table 13.2), it is a good approximation to replace $F(\hat{W}_c,\hat{W})$ by its period average in the generalized Rutherford equation [9,10]. [See Equation (13.49).] Now,

$\displaystyle \langle F(\hat{W}_c,\hat{W})\rangle = f_{polz1}\,\hat{W}_c^{\,2}\...
...2\varphi\rangle + \frac{\hat{W}_c^{\,2}}{\hat{W}^2}\,\langle\cos\varphi\rangle.$ (13.65)

However, to lowest order, $\langle \sin^2\varphi\rangle =1/2$, $\langle\sin\varphi\rangle = (1/2)\,\delta\,S_1$, and $\langle \cos\varphi\rangle = -(1/2)\,\delta\,C_1$. [See Equations (13.51) and (13.52).] Thus, we obtain

$\displaystyle \langle F(\hat{W}_c,\hat{W})\rangle = \frac{1}{2}\,f_{polz1}\,\ha...
...c^{\,4}\,\hat{W}^3-\frac{1}{2}\,\frac{\hat{W}_c^{\,2}}{\hat{W}^2}\,\delta\,C_1.$ (13.66)

As is clear from Table 13.1 and Equations (13.29)–(13.31), the first two terms on the right-hand side of the previous equation, which are ion polarization terms that arise from the modified ion flow in the immediate vicinity of the island chain when the chain is subject to an electromagnetic locking torque, are destabilizing (i.e., positive). On the other hand, the third term, which is due to the direct interaction between the island chain and the static resonant magnetic perturbation, is stabilizing (i.e., negative).

If we use the data given in Tables 13.1 and 13.2 to estimate the sizes of the three terms on the right-hand side of the previous equation in a tokamak fusion reactor then we discover that the first two terms are completely negligible compared with the third. In other words, the polarization terms associated with the modified ion flow in the vicinity of the island chain are much smaller than the term due to the direct interaction between the island chain and the static resonant magnetic perturbation. In this limit, we can use Equation (13.29) to write

$\displaystyle \langle F(\hat{W}_c,\hat{W})\rangle \simeq -\frac{1}{2^{11}}\left...
...\tau_\theta\,\tau_\varphi}}{\omega_0\,\tau_H^{\,2}}\right)C_1\,\hat{W}_c^{\,4}.$ (13.67)

Thus, we conclude that the interaction between a rotating island chain and a static magnetic perturbation gives rise to a stabilizing term in the chain's Rutherford equation that is independent of the island width, and scales as the square of the current flowing in the magnetic field-coil that generates the perturbation. [See Equation (13.3).] Making use of the data given in Table 13.2, we deduce that the stabilizing term takes the form

$\displaystyle \langle F(\hat{W}_c,\hat{W})\rangle \simeq -\left(\frac{\hat{W_c}}{0.109}\right)^4$ (13.68)

in a low-field tokamak fusion reactor (with $\omega_0=\omega_\ast$), and

$\displaystyle \langle F(\hat{W}_c,\hat{W})\rangle \simeq -\left(\frac{\hat{W_c}}{7.33\times 10^{-2}}\right)^4$ (13.69)

in a high-field tokamak fusion reactor. Of course, these estimates are only valid as long as the locking criteria (13.53) and (13.54) are not satisfied. Equations (13.53), (13.54), (13.68), and (13.69) suggest that the stabilizing effect of the interaction of a rotating magnetic island chain with a static resonant magnetic perturbation could completely stabilize a relatively narrow island chain, but would not be able to completely stabilize a wide island chain, because the chain would lock to the perturbation before complete stabilization was acheived. The complete stabilization of a narrow rotating magnetic island chain by a static resonant magnetic perturbation has been observed experimentally [16].

Let us now consider a locked magnetic island chain. Equations (13.59), (13.60), and (13.64) yield

$\displaystyle F(\hat{W}_c,\hat{W}) = f_{polz1}\left(\frac{\delta_{unlock}}{\del...
...delta_{unlock}}{\delta}\right)^2\right]^{1/2}\frac{\hat{W}_c^{\,2}}{\hat{W}^2}.$ (13.70)

Of course, the previous equation is only valid when $\delta>\delta_{unlock}$ (i.e., provided that the island chain remains locked to the static resonant magnetic perturbation). Assuming that $\delta\sim\delta_{unlock}$, $\hat{W}_c\lesssim 0.1$, and $\hat{W}\lesssim 0.1$, it is clear from the data given in Table 13.1 that the first two terms on the right-hand side of the previous equation are completely negligible compared to the third. In other words, just as is the case for a rotating island chain, the polarization terms in the Rutherford equation associated with the modified ion flow in the vicinity of the island chain are much smaller than the term due to the direct interaction between the island chain and the static resonant magnetic perturbation. In this case, we can write

$\displaystyle F(\hat{W}_c,\hat{W}) \simeq \left[1- \left(\frac{\delta_{unlock}}{\delta}\right)^2\right]^{1/2}\left(\frac{\hat{W}_c}{\hat{W}}\right)^2.$ (13.71)

Thus, it is clear that the destabilizing effect of a static resonant magnetic perturbation on locked island chain can be significant, especially in the limit in which $\hat{W}\ll \hat{W}_c$. This destabilizing effect has been observed experimentally [14].

Finally, experiments have been performed in which the rotation frequency of a magnetic island chain in a tokamak plasma is controlled by locking the chain to a rotating resonant magnetic perturbation [23]. More ambitious experiments have been performed in which the helical phase of a rotating resonant magnetic perturbation is actively controlled in such a manner that a magnetic island chain in a tokamak plasma is held in a permanent stabilizing phase relation with respect to the perturbation [19,21]. Unfortunately, neither of these control schemes is practical in a tokamak fusion reactor because the magnetic field-coils that generate the rotating magnetic perturbation would have to be placed inside the vacuum vessel (otherwise the oscillating magnetic field generated by the coils would be shielded from the plasma by eddy currents excited in the vessel), where they would be unacceptably vulnerable to neutrons produced by fusion reactions within the plasma, as well as large electromagnetic forces generated by plasma disruptions.