You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Plasma Phys. Control. Fusion 53 (2011) 093001<br />
Topical Review<br />
Figure 21. Eigenfunctions for <strong>the</strong> m = 1 mode for a Kadomtsev pr<strong>of</strong>ile for ka = 4[149, figure 8].<br />
3.3. Resistive MHD stability<br />
Because <strong>of</strong> <strong>the</strong> absence <strong>of</strong> singular surfaces where k · B o is zero <strong>the</strong> introduction <strong>of</strong> resistivity<br />
does not trigger tearing modes (with subsequent magnetic islands and reconnection), which<br />
would occur even at high values <strong>of</strong> Lundqvist number S (equation (3.4)). Instead both resistive<br />
diffusion and Joule heating are <strong>the</strong> effects that would modify ideal MHD, and values <strong>of</strong> S <strong>of</strong><br />
order 1 to 10 2 will be considered. This does not imply very cold plasmas, because <strong>the</strong> <strong>pinch</strong><br />
radius for a Z-<strong>pinch</strong> can be small ∼1 mm. Ra<strong>the</strong>r <strong>the</strong> motivation was to explain <strong>the</strong> apparently<br />
anomalous stability <strong>of</strong> initially cryogenic deuterium fibre experiments [151, 152] where <strong>the</strong><br />
small plasma radius leads to small S.<br />
Because <strong>of</strong> <strong>the</strong> perceived importance <strong>of</strong> Joule heating both in <strong>the</strong> equilibrium and in <strong>the</strong><br />
perturbed state we have to consider <strong>the</strong> stability <strong>of</strong> a time-evolving Joule heated Z-<strong>pinch</strong>.<br />
Culverwell and Coppins [133, 153] took an equilibrium constant radius <strong>pinch</strong> in which <strong>the</strong><br />
current varied as t 1/3 , satisfying <strong>the</strong> Haines–Hammel (see section 2.5) curve. S commenced at<br />
a low value <strong>of</strong> ∼1 and rose in time to 100. An initial perturbation is applied and its growth is<br />
followed by <strong>the</strong> linearized equations. At high S it becomes an ideal MHD mode; exponential<br />
growth <strong>of</strong> this mode implies that <strong>the</strong> logarithm <strong>of</strong> <strong>the</strong> kinetic energy <strong>of</strong> <strong>the</strong> perturbation will<br />
grow as v A t/a; since v a ∝ B ∝ t 1/3 it follows that at large S linear growth will be proportional<br />
to t 4/3 or S. Figure 22 shows how <strong>the</strong> growth <strong>of</strong> <strong>the</strong> MHD mode is delayed; specifically until<br />
a critical S is reached <strong>of</strong> about 50(ka) −0.86 . Concurrent work by Cochran and Robson [154]<br />
employed a nonlinear 2D time-dependent resistive MHD code, and similar conclusions were<br />
reached.<br />
However, a simplified analytic model by Lampe [155] claimed that growth rates were only<br />
reduced by a factor <strong>of</strong> 2 from ideal MHD and that an error occurred in Culverwell and Coppins<br />
arising from <strong>the</strong> neglect <strong>of</strong> <strong>pinch</strong> radius change in <strong>the</strong> perturbed Joule heating.<br />
38