Magnetohydrodynamic stability of plasmas with ideal and relaxed regions

Abstract

AbstractA unified energy principle approach is presented for analysing the magnetohydrodynamic (MHD) stability of plasmas consisting of multiple ideal and relaxed regions. The gauge a = ξ × B for the vector potential, a, of linearized perturbations is used, with the equilibrium magnetic field B obeying a Beltrami equation, ∇ × B = αB, in relaxed regions. In a region with such a force-free equilibrium Beltrami field we show that ξ obeys the same Euler–Lagrange equation whether ideal or relaxed MHD is used for perturbations, except in the neighbourhood of the magnetic surfaces where B · ∇ is singular. The difference at singular surfaces is analysed in cylindrical geometry: in ideal MHD only Newcomb's small solutions are allowed, whereas in relaxed MHD only the odd-parity large solution and even-parity small solution are allowed. A procedure for constructing global multi-region solutions in cylindrical geometry is presented. Focusing on the limit where the two interfaces approach each other arbitrarily closely, it is shown that the singular-limit problem encountered previously by Hole et al. in multi-region relaxed MHD is stabilized if the relaxed-MHD region between the coalescing interfaces is replaced by an ideal-MHD region. We then present a stable (k, pressure) phase-space plot, which allows us to determine the form a stable pressure and field profile must take in the region between the interfaces. From this knowledge, we conclude that there exists a class of single-interface plasmas that were found to be stable by Kaiser and Uecker, but are shown to be unstable when the interface is resolved.