Hostname: page-component-586b7cd67f-l7hp2 Total loading time: 0 Render date: 2024-11-23T10:52:15.793Z Has data issue: false hasContentIssue false

A one-dimensional tearing mode equation for pedestal stability studies in tokamaks

Published online by Cambridge University Press:  14 June 2018

J. W. Connor
Affiliation:
UKAEA-CCFE, Culham Science Centre, Abingdon, Oxon OX14 3DB, UK Imperial College of Science, Technology and Medicine, London SW7 2BZ, UK
R. J. Hastie*
Affiliation:
UKAEA-CCFE, Culham Science Centre, Abingdon, Oxon OX14 3DB, UK
C. Marchetto
Affiliation:
ISTITUTO di FISICA del Plasma, CNR, Via Roberto Cozzi, 53, 20125 Milano, Italy
C. M. Roach
Affiliation:
UKAEA-CCFE, Culham Science Centre, Abingdon, Oxon OX14 3DB, UK
*
Email address for correspondence: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Starting from expressions in Connor et al. (Phys. Fluids, vol. 31, 1988, p. 577), we derive a one-dimensional tearing equation similar to the approximate equation obtained by Hegna & Callen (Phys. Plasmas, vol. 1, 1994, p. 2308) and Nishimura et al. (Phys. Plasmas, vol. 5, 1998, p. 4292), but for more realistic toroidal equilibria. The intention is to use this approximation to explore the role of steep profiles, bootstrap currents and strong shaping in the vicinity of a separatrix, on the stability of tearing modes which are resonant in the H-mode pedestal region of finite aspect ratio, shaped cross-section tokamaks, e.g. the Joint European Torus (JET). We discuss how this one-dimensional model for tearing modes, which assumes a single poloidal harmonic for the perturbed poloidal flux, compares with a model that includes poloidal coupling Fitzpatrick et al. (Nucl. Fusion, vol. 33, 1993, p. 1533).

Type
Research Article
Copyright
© Cambridge University Press 2018 

1 Introduction

Edge localised modes (ELMs) are a ubiquitous feature of H-mode tokamak plasmas with important consequences for confinement and for transient heat loads on divertor target plates. Most theoretical models appeal to ideal magnetohydrodynamic (MHD) ballooning and peeling modes (Hegna et al. Reference Hegna, Connor, Hastie and Wilson1996; Connor et al. Reference Connor, Hastie, Wilson and Miller1998; Wilson et al. Reference Wilson, Connor, Field, Fielding, Miller, Lao, Ferron and Turnbull1999; Snyder et al. Reference Snyder, Wilson, Ferron, Lao, Leonard, Mossessian, Murakami, Osborne, Turnbull and Xu2004) as the trigger for ELMs. While this may well be the case for larger Type I ELMs, the smaller Type III may involve resistive ballooning modes (Connor Reference Connor1998). Furthermore it is unclear whether ideal peeling modes are ever unstable due to the presence of a separatrix in divertor tokamaks (Huysmans Reference Huysmans2005; Webster & Gimblett Reference Webster and Gimblett2009) or can lead to the required destruction of magnetic surfaces seen in resistive MHD simulations, e.g. ASDEX Team (Reference ASDEX Team1989). However, concerning the need for resistivity, one should mention a model for ELMs in which an unstable ideal peeling mode does play a part, triggering a Taylor relaxation in the edge plasma, thus involving reconnection. The relaxation region grows in size until the ideal mode becomes stable (Gimblett, Hastie & Helander Reference Gimblett, Hastie and Helander2006). An alternative possible explanation is that ELMs might be triggered by tearing modes being driven unstable by the large bootstrap current density that results from the pressure gradients in the H-mode pedestal.

The theory of tearing modes utilises asymptotic matching techniques (Furth, Killeen & Rosenbluth Reference Furth, Killeen and Rosenbluth1963). Thus solutions of the resistive equations (or those corresponding to more complex plasma models e.g. Antonsen & Coppi (Reference Antonsen and Coppi1981), Drake et al. (Reference Drake, Antonsen, Hassam and Gladd1983), Cowley, Kulsrud & Hahm (Reference Cowley, Kulsrud and Hahm1986), Pegoraro & Schep (Reference Pegoraro and Schep1986), Porcelli (Reference Porcelli1987), Fitzpatrick (Reference Fitzpatrick1989), Connor, Hastie & Zocco (Reference Connor, Hastie and Zocco2012)) that pertain near resonant surfaces, $m=nq(\unicode[STIX]{x1D70C}_{s})$ , are matched to solutions of ideal MHD equations that describe the regions away from the resonance to obtain a dispersion relation determining their stability. Here $m$ and $n$ are poloidal and toroidal mode numbers of the perturbation, $q(\unicode[STIX]{x1D70C})$ is the safety factor, $\unicode[STIX]{x1D70C}$ is a flux surface label with dimensions of length and $\unicode[STIX]{x1D70C}_{s}$ is the resonance position. This matching procedure involves obtaining the asymptotic forms of the ideal MHD solutions as $\unicode[STIX]{x1D70C}\rightarrow \unicode[STIX]{x1D70C}_{s}$ from both left and right, and the matching is characterised by a quantity $\unicode[STIX]{x1D6E5}^{\prime }$ . Stability of a mode is determined by comparing $\unicode[STIX]{x1D6E5}^{\prime }$ with $\unicode[STIX]{x1D6E5}_{\text{crit}}^{\prime }$ , a parameter that is determined from the solution of the equation describing the narrow layer around the resonance. The quantity $\unicode[STIX]{x1D6E5}_{\text{crit}}^{\prime }$ is usually a large positive number (Glasser, Greene & Johnson Reference Glasser, Greene and Johnson1975; Drake et al. Reference Drake, Antonsen, Hassam and Gladd1983; Cowley et al. Reference Cowley, Kulsrud and Hahm1986), but physics close to the resonance can make $\unicode[STIX]{x1D6E5}_{\text{crit}}^{\prime }$ negative: e.g. when microtearing modes are unstable, as has been reported for the plasma region around the H-mode pedestal in the Mega Ampère Spherical Tokamak (MAST) (Dickinson et al. Reference Dickinson, Roach, Saarelma, Scannell, Kirk and Wilson2012) and in the Joint European Torus (JET) (Hatch et al. Reference Hatch, Kotschenreuther, Mahajan, Valanju, Jenko, Told, Görler and Saarelma2016).

The linear theory of tearing instability in toroidal geometry (Connor et al. Reference Connor, Cowley, Hastie, Hender, Hood and Martin1988) is a complex problem, raising issues associated with the coupling of different poloidal harmonics and with the decoupling of resonances at different rational surfaces due to differing diamagnetic frequencies at such surfaces. Hegna & Callen (Reference Hegna and Callen1994) proposed a simple approximation that the perturbed poloidal flux has a single poloidal harmonic, of admittedly uncertain accuracy, to obtain a master equation for tearing instability, with similar one-dimensional (1-D) character to that holding in a straight cylinder. This equation was derived for equilibria with weakly shaped poloidal cross-section, and under the additional assumptions of large aspect ratio, low $\unicode[STIX]{x1D6FD}$ (where $\unicode[STIX]{x1D6FD}$ is the ratio of plasma pressure, $p$ , to the magnetic field energy density, $\unicode[STIX]{x1D6FD}=2\unicode[STIX]{x1D707}_{0}p/B^{2}$ ) and with the toroidal magnetic field greatly exceeding the poloidal field:

(1.1) $$\begin{eqnarray}\frac{I}{q}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\frac{q}{I}\langle g^{\unicode[STIX]{x1D713}\unicode[STIX]{x1D713}}\rangle \frac{\text{d}\tilde{A}}{\text{d}\unicode[STIX]{x1D713}}-\left[m^{2}\langle g^{\unicode[STIX]{x1D703}\unicode[STIX]{x1D703}}\rangle +\frac{m}{m-nq}I\langle \unicode[STIX]{x1D70E}\rangle ^{\prime }+\frac{m^{2}}{(m-nq)^{2}}\unicode[STIX]{x1D707}_{0}Ip^{\prime }\langle J\rangle ^{\prime }\right]\tilde{A}=0,\end{eqnarray}$$

where $\tilde{A}$ and $\unicode[STIX]{x1D713}$ are respectively the perturbed and equilibrium poloidal flux, the magnetic field is $\boldsymbol{B}=I\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}+\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}\times \unicode[STIX]{x1D735}\unicode[STIX]{x1D713}$ , $I$ is the toroidal field function $I=RB_{\unicode[STIX]{x1D719}}$ , $\unicode[STIX]{x1D719}$ is the toroidal coordinate, $\unicode[STIX]{x1D70E}=\unicode[STIX]{x1D707}_{0}j_{\Vert }/B$ , with the parallel current density $j_{\Vert }=\boldsymbol{j}\boldsymbol{\cdot }\boldsymbol{B}/B$ , $q$ is the safety factor, $^{\prime }$ denotes the radial gradient with respect to $\unicode[STIX]{x1D713}$ , $\langle Y\rangle$ is the flux surface average of $Y$ for any quantity $Y(\unicode[STIX]{x1D713},\unicode[STIX]{x1D703})$ ,

(1.2) $$\begin{eqnarray}\langle Y\rangle =\frac{1}{2\unicode[STIX]{x03C0}}\oint Y\,\text{d}\unicode[STIX]{x1D703}\end{eqnarray}$$

and the metric elements are $g^{\unicode[STIX]{x1D713}\unicode[STIX]{x1D713}}=|\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}|^{2}$ and $g^{\unicode[STIX]{x1D703}\unicode[STIX]{x1D703}}=|\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}|^{2}$ . The $\unicode[STIX]{x1D703}$ coordinate is a straight field line poloidal angle and $J=(\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}\times \unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D719})^{-1}=R^{2}q/I$ is the Jacobian. Nishimura, Callen & Hegna (Reference Nishimura, Callen and Hegna1998) presented numerical solutions of a similar equation, for a family of equilibrium profiles resembling those studied previously by Furth, Rutherford & Selberg (Reference Furth, Rutherford and Selberg1973) in a cylindrical geometry.

To assist the tearing mode stability analysis of the H-mode pedestal, in this paper we develop a 1-D ideal MHD equation for application to realistic, fully toroidal tokamak equilibria at high $\unicode[STIX]{x1D6FD}$ , thus generalising the earlier seminal works by Hegna & Callen (Reference Hegna and Callen1994) and Nishimura et al. (Reference Nishimura, Callen and Hegna1998). This contrasts with alternative approximate treatments described in Fitzpatrick et al. (Reference Fitzpatrick, Hastie, Martin and Roach1993), where the effect of poloidal mode coupling was calculated for toroidal equilibria of large aspect ratio, low $\unicode[STIX]{x1D6FD}$ and weak shaping, and approximate solutions with seven poloidal harmonics were used to obtain  $\unicode[STIX]{x1D6E5}^{\prime }$ .

2 A 1-D tearing mode equation

We start from equations (A5) and (A6) of Connor et al. (Reference Connor, Cowley, Hastie, Hender, Hood and Martin1988), which respectively govern the radial component of the displacement ( $\unicode[STIX]{x1D743}$ ), and the perturbed toroidal magnetic field: these quantities manifest themselves in Connor et al.’s variables $y=R_{0}f\unicode[STIX]{x1D743}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}$ and $z=R^{2}\unicode[STIX]{x1D6FF}\boldsymbol{B}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}/B_{0}$ . Here the equilibrium magnetic field is written as $\boldsymbol{B}=R_{0}B_{0}[g\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}+f\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}\times \unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}]$ , where $\unicode[STIX]{x1D70C}$ is a flux surface label with dimension of length, $B_{0}g(\unicode[STIX]{x1D70C})R_{0}/R$ is the full toroidal magnetic field, $B_{0}$ is the vacuum toroidal field at the major radius, $R_{0}$ , of the magnetic axis and $q=(\unicode[STIX]{x1D70C}/R_{0})(g/f)$ . The variable $y$ , which is related to the perturbed poloidal flux, now denoted $\unicode[STIX]{x1D713}$ , by $y=\unicode[STIX]{x1D713}/(m-nq)$ ,Footnote 1 is assumed to contain only a single poloidal harmonic, $\text{e}^{\text{i}m\unicode[STIX]{x1D703}}$ , where $\unicode[STIX]{x1D703}$ is the poloidal angle in straight field line coordinates. These equations can be used to generate the 1-D ideal MHD equation for $\unicode[STIX]{x1D713}$ .

Then equations (A5) and (A6) of Connor et al. (Reference Connor, Cowley, Hastie, Hender, Hood and Martin1988) take the form

(2.1) $$\begin{eqnarray}\text{i}\frac{\text{d}\unicode[STIX]{x1D713}}{\text{d}\unicode[STIX]{x1D70C}}\text{e}^{\text{i}m\unicode[STIX]{x1D703}}=-\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\left[\unicode[STIX]{x1D713}\text{e}^{\text{i}m\unicode[STIX]{x1D703}}\left(\text{i}T+\frac{U}{m-nq}\right)\right]+\left(Sz+\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}Q\frac{\unicode[STIX]{x2202}z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\right)\end{eqnarray}$$
(2.2) $$\begin{eqnarray}\displaystyle \left(\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}-\text{i}nq\right)\frac{\unicode[STIX]{x2202}z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}} & = & \displaystyle \unicode[STIX]{x1D713}\text{e}^{\text{i}m\unicode[STIX]{x1D703}}\left[\text{i}W+\frac{X}{m-nq}+(m-nq)V\right]-\text{i}\unicode[STIX]{x1D713}\text{e}^{\text{i}m\unicode[STIX]{x1D703}}\frac{\unicode[STIX]{x2202}V}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\nonumber\\ \displaystyle & & \displaystyle +\,U\frac{\unicode[STIX]{x2202}z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}-\left(\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}-\text{i}nq\right)\left[T^{\ast }\frac{\unicode[STIX]{x2202}z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\right],\end{eqnarray}$$

where the equilibrium quantities, $Q,S,T,U,V,W,X$ are defined in equation (A7) of Connor et al. (Reference Connor, Cowley, Hastie, Hender, Hood and Martin1988),

(2.3) $$\begin{eqnarray}\displaystyle & \displaystyle S=\text{i}n\unicode[STIX]{x1D70C}/R_{0}, & \displaystyle\end{eqnarray}$$
(2.4) $$\begin{eqnarray}\displaystyle & \displaystyle Q=-\text{i}\frac{R_{0}}{n\unicode[STIX]{x1D70C}}\frac{1}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}, & \displaystyle\end{eqnarray}$$
(2.5) $$\begin{eqnarray}\displaystyle & \displaystyle T=\frac{\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}+\text{i}\frac{R_{0}g^{\prime }}{n\unicode[STIX]{x1D70C}f}\frac{1}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}, & \displaystyle\end{eqnarray}$$
(2.6) $$\begin{eqnarray}\displaystyle & \displaystyle U=\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}\frac{R^{2}}{R_{0}^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}, & \displaystyle\end{eqnarray}$$
(2.7) $$\begin{eqnarray}\displaystyle & \displaystyle V=\text{i}\frac{n}{\unicode[STIX]{x1D70C}R_{0}}\frac{R_{0}^{2}}{R^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}-\text{i}\frac{R_{0}}{\unicode[STIX]{x1D70C}n}\left(\frac{g^{\prime }}{f}\right)^{2}\frac{1}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}} & \displaystyle\end{eqnarray}$$
(2.8) $$\begin{eqnarray}\displaystyle & \displaystyle W=\frac{2\unicode[STIX]{x1D707}_{0}p^{\prime }g^{\prime }}{B_{0}^{2}f^{3}}\frac{R^{2}}{R_{0}^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}-\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\frac{g^{\prime }}{f} & \displaystyle\end{eqnarray}$$
(2.9) $$\begin{eqnarray}\displaystyle & \displaystyle X=\text{i}\frac{n\unicode[STIX]{x1D707}_{0}p^{\prime }\unicode[STIX]{x1D70C}}{B_{0}^{2}f^{2}R_{0}}\left[\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\left(T^{\ast }\frac{R^{2}}{R_{0}^{2}}\right)+\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\frac{R^{2}}{R_{0}^{2}}-\frac{R^{2}}{R_{0}^{2}}\left(\frac{f^{\prime }}{f}-\frac{1}{\unicode[STIX]{x1D70C}}\right)-\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}\frac{R^{4}}{R_{0}^{4}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right], & \displaystyle\end{eqnarray}$$

with $^{\prime }$ now representing the radial derivative with respect to $\unicode[STIX]{x1D70C}$ . We note here that the above expressions were derived for equilibria of arbitrary aspect ratio, cross-sectional shape and $\unicode[STIX]{x1D6FD}$ . The method employed in the following analysis is rather general and does not assume that the second dependent variable, $z$ , is also of single poloidal harmonic structure.

To simplify the analysis we neglect the term involving $S$ , relative to $m^{2}Q$ in (2.1). This is equivalent to reducing the field line bending energy in a circular cylinder from the $(m^{2}+k_{z}^{2}r^{2})$ of the Newcomb (Reference Newcomb1960) analysis of stability in a linear pinch, to $m^{2}$ . In a torus this is equivalent to an assumption that $(\unicode[STIX]{x1D716}/q_{s})^{2}\ll 1$ , where $\unicode[STIX]{x1D716}$ is the local aspect ratio and $q_{s}=m/n$ is the value of the safety factor at the resonance. Since our focus will be on tearing modes which are resonant in the pedestal region of a tokamak of aspect ratio around $1/3$ , $q_{s}$ may be of order $4$ or greater, so this approximation would appear to introduce errors of only approximately $1\,\%$ .

The required 1-D tearing equation is now obtained by solving (2.1) for $\unicode[STIX]{x2202}z/\unicode[STIX]{x2202}\unicode[STIX]{x1D703}$ , inserting the result in (2.2) and taking the flux surface average. Thus:

(2.10) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}=\frac{\text{e}^{\text{i}m\unicode[STIX]{x1D703}}}{Q}\left[\frac{1}{m}\frac{\text{d}\unicode[STIX]{x1D713}}{\text{d}\unicode[STIX]{x1D70C}}+\unicode[STIX]{x1D713}\left(\text{i}T+\frac{U}{m-nq}\right)\right]-\frac{K(\unicode[STIX]{x1D70C})}{Q},\end{eqnarray}$$

where $K(\unicode[STIX]{x1D70C})$ is a flux surface dependent constant of integration to be determined by a periodicity constraint on $z(\unicode[STIX]{x1D70C},\unicode[STIX]{x1D703})$ . Thus

(2.11) $$\begin{eqnarray}K(\unicode[STIX]{x1D70C})=\frac{1}{m}\frac{\text{d}\unicode[STIX]{x1D713}}{\text{d}\unicode[STIX]{x1D70C}}\unicode[STIX]{x1D6FC}_{m}+\unicode[STIX]{x1D713}\left[\unicode[STIX]{x1D6FE}_{m}+\frac{\unicode[STIX]{x1D6FF}_{m}}{m-nq}\right],\end{eqnarray}$$

with

(2.12) $$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D6FC}_{m}=\frac{\langle \text{e}^{\text{i}m\unicode[STIX]{x1D703}}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }{\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }=\frac{\langle \cos (m\unicode[STIX]{x1D703})|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }{\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle } & \displaystyle\end{eqnarray}$$
(2.13) $$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D6FE}_{m}=\text{i}\frac{\langle \text{e}^{\text{i}m\unicode[STIX]{x1D703}}T|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }{\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }=-\frac{\langle \sin (m\unicode[STIX]{x1D703})\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}\rangle }{\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle } & \displaystyle\end{eqnarray}$$
(2.14) $$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D6FF}_{m}=\frac{\langle \text{e}^{\text{i}m\unicode[STIX]{x1D703}}U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }{\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }=\frac{\langle \cos (m\unicode[STIX]{x1D703})R^{2}\rangle }{R_{0}^{2}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}, & \displaystyle\end{eqnarray}$$

where the second form in (2.12)–(2.14) applies for equilibria which are symmetric above and below the median plane. Now, since the $m$ number for tearing modes which are resonant in the pedestal region of a tokamak is likely to be moderately large, the coefficients defined by $\unicode[STIX]{x1D6FC}_{m},\unicode[STIX]{x1D6FE}_{m}$ and $\unicode[STIX]{x1D6FF}_{m}$ in (2.12)–(2.14) may be very small unless there is strong shaping. Consequently, we can normally neglect the integration constant, $K(\unicode[STIX]{x1D70C})$ defined in (2.11). In § 2.1 we will investigate the consequences of retaining finite $K(\unicode[STIX]{x1D70C})$ .

Inserting the expression for $(\unicode[STIX]{x2202}z/\unicode[STIX]{x2202}\unicode[STIX]{x1D703})$ (in the $K(\unicode[STIX]{x1D70C})=0$ limit) into (2.2) and multiplying by the factor $\text{e}^{-\text{i}m\unicode[STIX]{x1D703}}$ , we take the flux surface average to obtain a 1-D tearing mode equation. Expressed in terms of the equilibrium quantities,  $Q,T,U,V,W$ and $X$ , this takes the form:

(2.15) $$\begin{eqnarray}\displaystyle & & \displaystyle \frac{(m-nq)}{m^{2}}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\left\langle \frac{1}{Q}\right\rangle \frac{\text{d}\unicode[STIX]{x1D713}}{\text{d}\unicode[STIX]{x1D70C}}\right]+\unicode[STIX]{x1D713}\frac{(m-nq)}{m}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\text{i}\left\langle \frac{T}{Q}\right\rangle +\frac{1}{(m-nq)}\left\langle \frac{U}{Q}\right\rangle \right]\nonumber\\ \displaystyle & & \displaystyle \qquad =\unicode[STIX]{x1D713}\left[(m-nq)\langle V\rangle +\text{i}\langle W\rangle +\frac{\langle X\rangle }{(m-nq)}+(m-nq)\left\langle \frac{T~T^{\ast }}{Q}\right\rangle +\frac{1}{(m-nq)}\left\langle \frac{U^{2}}{Q}\right\rangle \right.\nonumber\\ \displaystyle & & \displaystyle \qquad \quad \left.+\,\text{i}\left\langle \frac{U(T-T^{\ast })}{Q}\right\rangle \right].\end{eqnarray}$$

Now, writing $1/Q=\unicode[STIX]{x1D706}\unicode[STIX]{x1D70C}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}$ , where $\unicode[STIX]{x1D706}=\text{i}n/R_{0}$ and dividing through by $\unicode[STIX]{x1D706}(m-nq)/m^{2}$ , equation (2.15) takes the form of the second order differential equation:

(2.16) $$\begin{eqnarray}\displaystyle & & \displaystyle \frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\unicode[STIX]{x1D70C}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle \frac{\text{d}\unicode[STIX]{x1D713}}{\text{d}\unicode[STIX]{x1D70C}}\right]+m\unicode[STIX]{x1D713}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\text{i}\unicode[STIX]{x1D70C}\langle T|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle +\frac{\unicode[STIX]{x1D70C}}{(m-nq)}\langle U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle \right]\nonumber\\ \displaystyle & & \displaystyle \qquad =\unicode[STIX]{x1D713}\frac{m^{2}}{\unicode[STIX]{x1D706}}\left[\langle V\rangle +\text{i}\frac{\langle W\rangle }{(m-nq)}+\frac{\langle X\rangle }{(m-nq)^{2}}\right]\nonumber\\ \displaystyle & & \displaystyle \qquad \quad +\,\unicode[STIX]{x1D713}m^{2}\unicode[STIX]{x1D70C}\left[\langle TT^{\ast }|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle +\frac{\text{i}\langle U(T-T^{\ast })|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }{(m-nq)}+\frac{\langle U^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle }{(m-nq)^{2}}\right],\end{eqnarray}$$

which is of the same structure as the equation derived by Hegna & Callen (Reference Hegna and Callen1994), namely

(2.17) $$\begin{eqnarray}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[A(\unicode[STIX]{x1D70C})\frac{\text{d}\unicode[STIX]{x1D713}}{\text{d}\unicode[STIX]{x1D70C}}\right]-\left[B(\unicode[STIX]{x1D70C})+\frac{mC(\unicode[STIX]{x1D70C})}{(m-nq)}+\frac{m^{2}D(\unicode[STIX]{x1D70C})}{(m-nq)^{2}}\right]\unicode[STIX]{x1D713}=0,\end{eqnarray}$$

where, on inserting the definitions (2.4)–(2.9),

(2.18) $$\begin{eqnarray}\displaystyle & \displaystyle A=\unicode[STIX]{x1D70C}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle & \displaystyle\end{eqnarray}$$
(2.19) $$\begin{eqnarray}\displaystyle & \displaystyle B=\frac{m^{2}}{\unicode[STIX]{x1D70C}}\left[\left\langle \frac{R_{0}^{2}}{R^{2}}\frac{1}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle +\unicode[STIX]{x1D70C}^{2}\left\langle \frac{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle \right]=m^{2}\unicode[STIX]{x1D70C}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}|^{2}\rangle & \displaystyle\end{eqnarray}$$
(2.20) $$\begin{eqnarray}\displaystyle & \displaystyle C=-q\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\frac{R_{0}g^{\prime }}{f}+\frac{R_{0}}{fg}\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}}\left\langle \frac{R^{2}}{R_{0}^{2}}\right\rangle \right] & \displaystyle\end{eqnarray}$$
(2.21) $$\begin{eqnarray}\displaystyle & \displaystyle D=\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}\left[\unicode[STIX]{x1D70C}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left\langle \frac{R^{2}}{R_{0}^{2}}\right\rangle -\left\langle \frac{R^{2}}{R_{0}^{2}}\right\rangle \left(\frac{\unicode[STIX]{x1D70C}g^{\prime }}{g}\right)\right]. & \displaystyle\end{eqnarray}$$

Some details of the derivation of (2.18)–(2.21) are given in appendix A, and in appendix B we express equation (2.17) in terms of the variables used in equation (26) of Hegna & Callen (Reference Hegna and Callen1994).

2.1 Consequences of finite K

We now return to (2.10) and (2.11) and construct the additional terms that will appear in the tearing equation when we retain terms with finite $K(\unicode[STIX]{x1D70C})$ . After lengthy, but straightforward, further analysis, we find that each of the coefficients $A(\unicode[STIX]{x1D70C}),B(\unicode[STIX]{x1D70C}),C(\unicode[STIX]{x1D70C})$ and $D(\unicode[STIX]{x1D70C})$ is modified by an additional contribution, which we shall denote by a circumflex. Thus

(2.22) $$\begin{eqnarray}\displaystyle \left.\begin{array}{@{}c@{}}A\rightarrow A(\unicode[STIX]{x1D70C})-\hat{A}(\unicode[STIX]{x1D70C}),\\ B\rightarrow B(\unicode[STIX]{x1D70C})-\hat{B}(\unicode[STIX]{x1D70C}),\\ C\rightarrow C(\unicode[STIX]{x1D70C})-{\hat{C}}(\unicode[STIX]{x1D70C}),\\ D\rightarrow D(\unicode[STIX]{x1D70C})-\hat{D}(\unicode[STIX]{x1D70C}),\end{array}\right\} & & \displaystyle\end{eqnarray}$$

with

(2.23) $$\begin{eqnarray}\displaystyle & \displaystyle \hat{A}(\unicode[STIX]{x1D70C})=A(\unicode[STIX]{x1D70C})|\unicode[STIX]{x1D6FC}_{m}|^{2}, & \displaystyle\end{eqnarray}$$
(2.24) $$\begin{eqnarray}\displaystyle & \displaystyle \hat{B}(\unicode[STIX]{x1D70C})=m^{2}A(\unicode[STIX]{x1D70C})\unicode[STIX]{x1D6FE}_{m}^{2}-m\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}[A(\unicode[STIX]{x1D70C})\unicode[STIX]{x1D6FE}_{m}\unicode[STIX]{x1D6FC}_{m}], & \displaystyle\end{eqnarray}$$
(2.25) $$\begin{eqnarray}\displaystyle & \displaystyle {\hat{C}}(\unicode[STIX]{x1D70C})=-q\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\frac{\unicode[STIX]{x1D6FC}_{m}\unicode[STIX]{x1D6FF}_{m}A(\unicode[STIX]{x1D70C})}{q}\right], & \displaystyle\end{eqnarray}$$
(2.26) $$\begin{eqnarray}\displaystyle & \displaystyle \hat{D}(\unicode[STIX]{x1D70C})=A(\unicode[STIX]{x1D70C})\unicode[STIX]{x1D6FF}_{m}\left(\unicode[STIX]{x1D6FF}_{m}-\frac{s\unicode[STIX]{x1D6FC}_{m}}{\unicode[STIX]{x1D70C}}\right) & \displaystyle\end{eqnarray}$$

where $s=\unicode[STIX]{x1D70C}q^{\prime }/qs=\unicode[STIX]{x1D70C}q^{\prime }/q$ is the magnetic shear. $\hat{A}$ , $\hat{B}$ , ${\hat{C}}$ and $\hat{D}$ are small in the large $m$ limit because the numerators in the definitions of $\unicode[STIX]{x1D6FC}_{m}$ , $\unicode[STIX]{x1D6FF}_{m}$ , $\unicode[STIX]{x1D6FE}_{m}$ (see (2.12)–(2.14)) must vanish both at high $m$ , or with weak shaping. At a fixed finite $m$ these terms can, however, become more important with stronger shaping (e.g. as one approaches the separatrix).

2.2 Comparison with earlier results

The Hegna–Callen equation represented a significant advance on earlier work by making possible a simple 1-D tearing analysis of large aspect ratio toroidal equilibria with weakly shaped poloidal cross-sections. Our derivation has not only extended the validity of the 1-D equation to finite aspect ratio equilibria, subject to $(\unicode[STIX]{x1D716}/q_{s})^{2}\ll 1$ , with arbitrary poloidal shaping, but it has also revealed the presence of new terms arising from finite values of the integration constant $K(\unicode[STIX]{x1D70C})$ . These additional terms of (2.23)–(2.26) have no counterpart in Hegna & Callen (Reference Hegna and Callen1994) or Nishimura et al. (Reference Nishimura, Callen and Hegna1998), but they are small unless there is strong shaping containing poloidal harmonics that couple to the mode number, $m$ .

We now compare our tearing equation (2.17) with Hegna & Callen (Reference Hegna and Callen1994) and Nishimura et al. (Reference Nishimura, Callen and Hegna1998). We begin by transforming from the Hegna–Callen equilibrium variables, $I$ and $\unicode[STIX]{x1D713}$ , to the $g,f,\unicode[STIX]{x1D70C}$ variables of the present work. Thus:

(2.27) $$\begin{eqnarray}\displaystyle & \displaystyle I\rightarrow R_{0}B_{0}g(\unicode[STIX]{x1D70C}), & \displaystyle\end{eqnarray}$$
(2.28) $$\begin{eqnarray}\displaystyle & \displaystyle \frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\rightarrow \frac{1}{\unicode[STIX]{x1D713}^{\prime }}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}, & \displaystyle\end{eqnarray}$$
(2.29) $$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D713}^{\prime }(\unicode[STIX]{x1D70C})\rightarrow R_{0}B_{0}f(\unicode[STIX]{x1D70C}). & \displaystyle\end{eqnarray}$$

The coefficients $A,B,C$ and $D$ can then be identified in equation (26) of Hegna & Callen (Reference Hegna and Callen1994) and compared to (2.18)–(2.21). This shows agreement in the expressions for $A$ and $B$ , close agreement on $C$ , but not for $D$ . Since

(2.30) $$\begin{eqnarray}\unicode[STIX]{x1D70E}=\frac{1}{f}\left(g^{\prime }+\frac{g\unicode[STIX]{x1D707}_{0}p^{\prime }}{B^{2}}\right),\end{eqnarray}$$

one can indeed write $C\propto (\unicode[STIX]{x2202}\langle \unicode[STIX]{x1D70E}\rangle /\unicode[STIX]{x2202}\unicode[STIX]{x1D70C})$ if $B\simeq B_{\unicode[STIX]{x1D719}}$ , as in Hegna & Callen (Reference Hegna and Callen1994). There is some similarity with the expression for $D$ that appears in equation (19) of Nishimura et al. (Reference Nishimura, Callen and Hegna1998), where special equilibria with $g=$ constant were studied so that the last term in (2.21) is absent, but nevertheless their $D\propto n^{2}q^{2}$ rather than $m^{2}$ , and so it differs away from the resonance.

As noted by Hegna, Callen and Nishimura, there is an important comparison for the expression given in (2.21) for $D(\unicode[STIX]{x1D70C})$ . This is associated with the Mercier stability criterion, $D_{M}<0$ , for the ideal MHD stability of a mode localised around a rational surface (Mercier Reference Mercier1960). Glasser et al. (Reference Glasser, Greene and Johnson1975) showed $D_{M}$ Footnote 2 plays an important role in the theory of tearing mode stability in a torus. They found the asymptotic form of the ideal MHD solutions as $\unicode[STIX]{x1D70C}\rightarrow \unicode[STIX]{x1D70C}_{s}$ is

(2.31) $$\begin{eqnarray}\unicode[STIX]{x1D713}\sim c_{0}\left|x-1\right|^{\unicode[STIX]{x1D708}_{-}}+c_{1}\left|x-1\right|^{\unicode[STIX]{x1D708}_{+}},\end{eqnarray}$$

where $x=\unicode[STIX]{x1D70C}/\unicode[STIX]{x1D70C}_{s}$ , constants $c_{0}$ and $c_{1}$ have different values to the left and right of the resonance, and the Mercier indices $\unicode[STIX]{x1D708}_{\pm }$ have values:

(2.32) $$\begin{eqnarray}\unicode[STIX]{x1D708}_{\pm }={\textstyle \frac{1}{2}}\pm \sqrt{-D_{M}}.\end{eqnarray}$$

This serves to define a generalised $\unicode[STIX]{x1D6E5}^{\prime }$

(2.33) $$\begin{eqnarray}\unicode[STIX]{x1D6E5}^{\prime }=\left.\frac{c_{1}}{c_{0}}\right|_{R}+\left.\frac{c_{1}}{c_{0}}\right|_{L},\end{eqnarray}$$

where $R$ and $L$ denote locations immediately to the right and left of the resonance, respectively. This expression, obtained from the ideal MHD solution, must be matched to the analogous quantity arising from the inner resonant layer solution, to obtain the tearing mode dispersion relation.

Using the results in Glasser et al. (Reference Glasser, Greene and Johnson1975) and Connor et al. (Reference Connor, Cowley, Hastie, Hender, Hood and Martin1988)Footnote 3 we find that, at the tearing mode resonance, $D$ of (2.21) should be compared to $-(As^{2}/\unicode[STIX]{x1D70C}^{2})((1/4)+D_{M})$ , where

(2.34) $$\begin{eqnarray}\displaystyle D_{M} & = & \displaystyle -\frac{1}{4}+E+F+H\nonumber\\ \displaystyle & = & \displaystyle -\frac{1}{4}+\frac{q}{q^{\prime }}\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}\left\langle \frac{R^{2}}{R_{0}^{2}}\frac{1}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle -\frac{q^{2}}{q^{\prime 2}}\left(\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}\right)^{2}\left\langle \frac{R^{2}}{R_{0}^{2}}\frac{1}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle ^{2}\nonumber\\ \displaystyle & & \displaystyle -\,\left(\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}}\right)\frac{1}{q^{\prime 2}}\left(\frac{\unicode[STIX]{x1D70C}}{R_{0}f}\right)^{2}\times \left\langle \frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\left(\frac{R^{2}}{R_{0}^{2}}\right)-\frac{R^{2}}{R_{0}^{2}}\frac{\unicode[STIX]{x1D70C}}{f}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left(\frac{f}{\unicode[STIX]{x1D70C}}\right)\right\rangle \left\langle \frac{B^{2}R^{2}}{B_{0}^{2}R_{0}^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle \nonumber\\ \displaystyle & & \displaystyle +\,\left(\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}f^{2}q^{\prime }}\right)^{2}\left(\frac{\unicode[STIX]{x1D70C}}{R_{0}f}\right)^{2}\left\langle \frac{R^{4}}{R_{0}^{4}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle \left\langle \frac{B^{2}R^{2}}{B_{0}^{2}R_{0}^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}}\right\rangle ,\end{eqnarray}$$

and the quantities $E,~F$ and $H$ are defined in Glasser et al. (Reference Glasser, Greene and Johnson1975). (In a later paper, Glasser, Greene & Johnson (Reference Glasser, Greene and Johnson1976) showed that for a large aspect ratio circular cross-section plasma:

(2.35) $$\begin{eqnarray}E+F+H=\frac{2\unicode[STIX]{x1D70C}\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}}\frac{q^{2}-1}{s^{2}},\end{eqnarray}$$

where the important factor $q^{2}-1$ removes, for $q>1$ , the possibility of the instability predicted by Suydam (Reference Suydam1958) in a straight cylinder.) Thus we can write:

(2.36) $$\begin{eqnarray}\frac{1}{4}+D_{M}\propto \frac{\unicode[STIX]{x1D70C}\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}}\frac{\unicode[STIX]{x1D705}_{\text{eff}}}{s^{2}},\end{eqnarray}$$

with the ‘effective’ curvature, $\unicode[STIX]{x1D705}_{\text{eff}}$ , deduced from (2.34). However, Hegna & Callen (Reference Hegna and Callen1994), perhaps seeking a $D$ consistent with this argument, assumed $\unicode[STIX]{x1D705}_{\text{eff}}$ was the surface-averaged normal curvature, $\unicode[STIX]{x1D705}_{n}$ , and, furthermore, that $\unicode[STIX]{x1D705}_{n}\propto V^{\prime \prime }=(\text{d}\langle J\rangle /\text{d}\unicode[STIX]{x1D70C})$ , where $\langle J\rangle =(\langle R^{2}\rangle q/R_{0}B_{0}g)$ , to obtain the following result for $D$ :

(2.37) $$\begin{eqnarray}D_{HC}\propto \frac{\unicode[STIX]{x1D70C}\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}s^{2}}\frac{\text{d}\langle J\rangle }{\text{d}\unicode[STIX]{x1D70C}}\propto \frac{\unicode[STIX]{x1D70C}\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}s^{2}}\frac{1}{R_{0}}\frac{q}{g}\left(\frac{\text{d}\langle R^{2}\rangle }{\text{d}\unicode[STIX]{x1D70C}}+\langle R^{2}\rangle \left(\frac{q^{\prime }}{q}-\frac{g^{\prime }}{g}\right)\right).\end{eqnarray}$$

However, at low $\unicode[STIX]{x1D6FD}$ and with $B_{\unicode[STIX]{x1D719}}\simeq B$ (e.g. at large aspect ratio),

(2.38) $$\begin{eqnarray}\unicode[STIX]{x1D705}_{n}\propto V^{\prime \prime }-\frac{\langle R^{2}\rangle }{R_{0}B_{0}}\frac{q^{\prime }}{g}\end{eqnarray}$$

(Connor, Hastie & Helander Reference Connor, Hastie and Helander2009), so that their argument should have implied

(2.39) $$\begin{eqnarray}D_{HC}\rightarrow D\propto \frac{\unicode[STIX]{x1D70C}\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}s^{2}}\frac{1}{R_{0}}\frac{q}{g}\left(\frac{\text{d}\langle R^{2}\rangle }{\text{d}\unicode[STIX]{x1D70C}}-\langle R^{2}\rangle \frac{g^{\prime }}{g}\right).\end{eqnarray}$$

Equation (2.39) is indeed consistent with our expression for $D$ in (2.21), and also with the work of Nishimura et al. (Reference Nishimura, Callen and Hegna1998) in the special case $g^{\prime }=0$ that they considered. Equations (2.21), (2.39) are not, however, consistent with $D=-(As^{2}/\unicode[STIX]{x1D70C}^{2})((1/4)+D_{M})$ , since $\unicode[STIX]{x1D705}_{\text{eff}}\neq \unicode[STIX]{x1D705}_{n}$ .Footnote 4 We should not expect $D$ to be exactly equal to $-(As^{2}/\unicode[STIX]{x1D70C}^{2})((1/4)+D_{M})$ , because the ideal instability investigated by Mercier, and later by Greene & Johnson (Reference Greene and Johnson1962) using Hamada coordinates, is a mode with a range of coupled poloidal harmonics, whereas $\unicode[STIX]{x1D713}$ of the envisaged tearing mode, has an isolated single poloidal harmonic.

It would be inconsistent with the ‘single poloidal harmonic’ assumption to simply replace $D(\unicode[STIX]{x1D70C})$ by the value corresponding to $D_{M}$ in (2.17); although the use of $D_{M}$ would capture the poloidal mode coupling effects close to the singular surface that can have a profound effect on the Mercier indices, which in turn influence the value of the generalised $\unicode[STIX]{x1D6E5}^{\prime }$ stability parameter (Glasser et al. Reference Glasser, Greene and Johnson1975).

3 Conclusions

Within the foregoing sections we have assumed that the perturbed poloidal flux function, $\unicode[STIX]{x1D713}(\unicode[STIX]{x1D70C},\unicode[STIX]{x1D703})$ , contained only one poloidal harmonic, $\text{e}^{\text{i}m\unicode[STIX]{x1D703}}$ . However our solution for the variable $z$ , equation (2.10), contains a full spectrum of poloidal harmonics. Under these assumptions we have extended the validity of the tearing equation proposed by Hegna & Callen (Reference Hegna and Callen1994) to axisymmetric equilibria of arbitrary aspect ratio and arbitrary $\unicode[STIX]{x1D6FD}$ . In doing so we have only made use of the approximation, $\unicode[STIX]{x1D716}^{2}/q_{s}^{2}\ll 1$ . This would certainly rule out the use of the resulting 1-D equation for studying internal kink type disruptions in tokamaks (where the $m=n=1$ harmonic plays a crucial role), but should prove to be an accurate approximation for modes which are resonant in the pedestal region of a tokamak in H-mode. An unexpected result of this calculation has been the appearance of a new set of terms arising from the effect of the integration constant $K$ (denoted by $\hat{A},~\hat{B},~{\hat{C}}$ and $\hat{D}$ ). However, it appears unlikely that such terms will play a significant role in determining tearing stability since they are normally negligibly small, except perhaps in very strongly shaped cross-sections or, e.g. in the vicinity of a separatrix boundary. For simplicity we ignored these extra terms in (2.17).

It is also clear from the foregoing derivation of a 1-D equation that the pressure gradient term, $D(\unicode[STIX]{x1D70C})$ of (2.17), differs from the quantity $-(As^{2}/\unicode[STIX]{x1D70C}^{2})((1/4)+D_{M})$ that would be expected in general tearing mode theory, as the singular surface is approached. The difference arises because the derivation of (2.17) is based on a single poloidal harmonic assumption, whereas retention of the coupled poloidal harmonics is required to capture the true value in the limit as $\unicode[STIX]{x1D70C}\rightarrow \unicode[STIX]{x1D70C}_{s}$ . The approach outlined in Fitzpatrick et al. (Reference Fitzpatrick, Hastie, Martin and Roach1993) retains seven coupled poloidal harmonics, but its restrictions to weak shaping and low $\unicode[STIX]{x1D6FD}$ severely impede its application to the pedestal. The single poloidal harmonic approach outlined in this paper accommodates strong shaping and $\unicode[STIX]{x1D6FD}$ effects, but neglects poloidal mode coupling that is needed to describe $D_{M}$ at the resonance and that may be important more globally. Nevertheless, for $\unicode[STIX]{x1D6E5}^{\prime }$ calculations at the foot of the pedestal where $s^{2}$ becomes large near a separatrix boundary, both the exact Mercier indices and the approximate (1-D) ones return to similar, low $\unicode[STIX]{x1D6FD}$ , values (of 0 and 1), and the 1-D approximation may give a good indication of tearing instability in a rather simple manner.

Numerical investigations of H-mode equilibria are presently underway.

Acknowledgements

This work has been carried out under the RCUK Energy Programme (grant number EP/P012450/1) with partial funding through the framework of the EUROfusion Consortium from the Euratom research and training programme 2014-2018 under grant agreement no. 633053. To obtain further information on the data and models underlying this paper please contact Email address for correspondence: . The views and opinions expressed herein do not necessarily reflect those of the European Commission. The authors are grateful to C. Ham for helpful discussions.

Appendix A

We can generate unique expressions for the coefficients $B,~C$ and $D$ , by exploiting the fact that all toroidal mode number dependencies in the 1-D tearing (2.17) can be expressed as powers, up to quadratic, of $m/(m-nq)$ .

First, we collect all three terms in (2.16), that include parts proportional to $(m/(m-nq))^{2}$ and contribute to the coefficient $D(\unicode[STIX]{x1D70C})$ in (2.17), namely;

(A 1) $$\begin{eqnarray}\frac{1}{\unicode[STIX]{x1D706}}\langle X\rangle +\unicode[STIX]{x1D70C}\langle U^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle -\unicode[STIX]{x1D70C}\frac{n}{m}q^{\prime }\langle U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle .\end{eqnarray}$$

Now replacing $n$ by the identity $(nq-m+m)/q$ and using $s=\unicode[STIX]{x1D70C}q^{\prime }/q$ , this expression becomes

(A 2) $$\begin{eqnarray}\frac{1}{\unicode[STIX]{x1D706}}\langle X\rangle +\unicode[STIX]{x1D70C}\langle U^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle -s\langle U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle +\frac{(m-nq)}{m}s\langle U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle ,\end{eqnarray}$$

where the first three terms yield (2.21) for $D$ and the last term now contributes to the expression for the coefficient $C$ , rather than $D$ . Three different terms from (2.16) and the final term of equation (A 2), contribute the term in (2.17) that is proportional to $m/(m-nq)$ , with the following factor in the coefficient:

(A 3) $$\begin{eqnarray}\frac{\text{i}m}{\unicode[STIX]{x1D706}}\langle W\rangle +\text{i}m\unicode[STIX]{x1D70C}\langle U(T-T^{\ast })|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle -\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left(\unicode[STIX]{x1D70C}\langle U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle \right)+s\langle U|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}\rangle ,\end{eqnarray}$$

where the last term is the contribution from equation (A 2) above. Using (2.5), (2.6) and (2.8) for $T$ , $U$ and $W$ , the expression in (A 3) becomes:

(A 4) $$\begin{eqnarray}-\frac{m}{n}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\frac{R_{0}g^{\prime }}{f}\right]-q\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\frac{R_{0}}{fg}\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}}\left\langle \frac{R^{2}}{R_{0}^{2}}\right\rangle \right].\end{eqnarray}$$

Now, on replacing $m$ by the identity $m-nq+nq$ , we obtain the following expression:

(A 5) $$\begin{eqnarray}C=-q\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left[\frac{R_{0}g^{\prime }}{f}+\frac{R_{0}}{fg}\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{B_{0}^{2}}\left\langle \frac{R^{2}}{R_{0}^{2}}\right\rangle \right]-\frac{(m-nq)}{n}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D70C}}\left(\frac{R_{0}g^{\prime }}{f}\right),\end{eqnarray}$$

where the first two terms coincide with (2.20) for $C(\unicode[STIX]{x1D70C})$ , and the third term contributes to the coefficient $B(\unicode[STIX]{x1D70C})$ and exactly cancels the remaining $n$ dependence in $B$ , leading to (2.19) for $B(\unicode[STIX]{x1D70C})$ .

To demonstrate the second equality in (2.19) we consider cylindrical toroidal coordinates $R,Z,\unicode[STIX]{x1D719}$ . The Jacobian for the transformation $(R,Z)\rightarrow (\unicode[STIX]{x1D70C},\unicode[STIX]{x1D703})$ is:

(A 6) $$\begin{eqnarray}J=\frac{\unicode[STIX]{x1D70C}R}{R_{0}}=\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}-\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}.\end{eqnarray}$$

We can obtain $\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}$ and $\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}$ , using

(A 7) $$\begin{eqnarray}\displaystyle \left.\begin{array}{@{}c@{}}\displaystyle \unicode[STIX]{x1D735}R=\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}+\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\\ \displaystyle \unicode[STIX]{x1D735}Z=\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}+\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\end{array}\right\} & & \displaystyle\end{eqnarray}$$

and deduce:

(A 8) $$\begin{eqnarray}\displaystyle \left.\begin{array}{@{}c@{}}\displaystyle J^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}|^{2}=\left(\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\right)^{2}+\left(\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\right)^{2}\\ \displaystyle J^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}|^{2}=\left(\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\right)^{2}+\left(\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\right)^{2}\\ \displaystyle J^{2}\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}=-\left[\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\frac{\unicode[STIX]{x2202}R}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}+\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D703}}\frac{\unicode[STIX]{x2202}Z}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}\right].\end{array}\right\} & & \displaystyle\end{eqnarray}$$

Squaring and adding (A 6) and (A 7) using (A 8) one finds:

(A 9) $$\begin{eqnarray}\frac{R^{2}}{R_{0}^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}^{2}|}+\unicode[STIX]{x1D70C}^{2}\frac{\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}}{|\unicode[STIX]{x1D735}\unicode[STIX]{x1D70C}^{2}|}=\unicode[STIX]{x1D70C}^{2}|\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}|^{2}.\end{eqnarray}$$

Appendix B

The 1-D tearing (2.17) is expressed in terms of equilibrium variables $\unicode[STIX]{x1D70C},g$ and $f$ . More familiar variables are the equilibrium poloidal flux $\unicode[STIX]{x1D713}$ and $I(\unicode[STIX]{x1D713})$ as used by Hegna & Callen (Reference Hegna and Callen1994). These are related by (2.27)–(2.29). In this appendix we give the form that (2.17) takes when expressed in these Hegna–Callen variables. Of the four terms in (2.17) we find:

(B 1) $$\begin{eqnarray}\displaystyle & \displaystyle Term~1\rightarrow \frac{I\unicode[STIX]{x1D70C}}{q}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\left[\frac{q}{I}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}|^{2}\rangle \frac{\text{d}\tilde{A}}{\text{d}\unicode[STIX]{x1D713}}\right], & \displaystyle\end{eqnarray}$$
(B 2) $$\begin{eqnarray}\displaystyle & \displaystyle Term~2\rightarrow -m^{2}\unicode[STIX]{x1D70C}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}|^{2}\rangle \tilde{A}, & \displaystyle\end{eqnarray}$$
(B 3) $$\begin{eqnarray}\displaystyle & \displaystyle Term~3\rightarrow \frac{m\unicode[STIX]{x1D70C}I}{(m-nq)}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\left[I^{\prime }(\unicode[STIX]{x1D713})+\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }(\unicode[STIX]{x1D713})}{I(\unicode[STIX]{x1D713})}\langle R^{2}\rangle \right]\tilde{A}, & \displaystyle\end{eqnarray}$$
(B 4) $$\begin{eqnarray}\displaystyle & \displaystyle Term~4\rightarrow -\frac{\unicode[STIX]{x1D70C}Im^{2}\unicode[STIX]{x1D707}_{0}p^{\prime }}{(m-nq)^{2}}\left[\frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\left\langle \frac{R^{2}}{I}\right\rangle \right]\tilde{A}, & \displaystyle\end{eqnarray}$$

where, as in the work of Hegna and Callen, the dependent variable $\tilde{A}$ is the, single poloidal harmonic, tearing mode eigenfunction and $^{\prime }$ denotes the radial derivative with respect to $\unicode[STIX]{x1D713}$ . Finally, on multiplying through by the factor $q/\unicode[STIX]{x1D70C}I$ we obtain the 1-D tearing equation in a rather simple form:

(B 5) $$\begin{eqnarray}\displaystyle & & \displaystyle \frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\left[\frac{q}{I}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}|^{2}\rangle \frac{\text{d}\tilde{A}}{\text{d}\unicode[STIX]{x1D713}}\right]-\left\{\frac{m^{2}q}{I}\langle |\unicode[STIX]{x1D735}\unicode[STIX]{x1D703}|^{2}\rangle -\frac{mq}{(m-nq)}\frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\left[I^{\prime }+\frac{\unicode[STIX]{x1D707}_{0}p^{\prime }}{I}\langle R^{2}\rangle \right]\right.\nonumber\\ \displaystyle & & \displaystyle \qquad \quad \left.+\,\frac{m^{2}q\unicode[STIX]{x1D707}_{0}p^{\prime }}{(m-nq)^{2}}\left[\frac{\text{d}}{\text{d}\unicode[STIX]{x1D713}}\left\langle \frac{R^{2}}{I}\right\rangle \right]\right\}\tilde{A}=0.\end{eqnarray}$$

Footnotes

1 The behaviour of $\unicode[STIX]{x1D713}$ near the resonant surface at $m=nq(\unicode[STIX]{x1D70C})$ , is a combination of large and small solutions (Glasser et al. Reference Glasser, Greene and Johnson1975), and it is this combination that must be matched to the solution in the inner resonant layer.

2 The quantity labelled $D_{M}$ here is precisely the object denoted by $D_{I}$ in Glasser et al. (Reference Glasser, Greene and Johnson1975).

3 A factor $1/f^{2}$ was missed from the final term of equation (B3) of Connor et al. (Reference Connor, Cowley, Hastie, Hender, Hood and Martin1988) due to a typographical error, and this is correctly included here.

4 A more detailed discussion of the relation of $D_{M}$ to $\unicode[STIX]{x1D705}_{n}$ is given by Johnson & Greene (Reference Johnson and Greene1967).

References

Antonsen, T. M. & Coppi, B. 1981 Non-asymptotic theory of collisionless reconnecting modes. Phys. Lett. A 81, 335.Google Scholar
ASDEX Team, 1989 The H-Mode of ASDEX. Nucl. Fusion 29, 1959.Google Scholar
Connor, J. W. 1998 Edge-localized modes – physics and theory. Plasma Phys. Control. Fusion 40, 531.Google Scholar
Connor, J. W., Cowley, S. C., Hastie, R. J., Hender, T. C., Hood, A. & Martin, T. J. 1988 Tearing modes in toroidal geometry. Phys. Fluids 31, 577.Google Scholar
Connor, J. W., Hastie, R. J. & Helander, P. 2009 Linear tearing mode stability equations for a low collisionality toroidal plasma. Plasma Phys. Control. Fusion 51, 015009.Google Scholar
Connor, J. W., Hastie, R. J., Wilson, H. R. & Miller, R. L. 1998 Magnetohydrodynamic stability of tokamak edge plasmas. Phys. Plasmas 5, 2687.Google Scholar
Connor, J. W., Hastie, R. J. & Zocco, A. 2012 Unified theory of the semi-collisional tearing mode and internal kink mode in a hot tokamak: implications for sawtooth modelling. Plasma Phys. Control. Fusion 54, 035003.Google Scholar
Cowley, S. C., Kulsrud, R. M. & Hahm, T. S. 1986 Linear stability of tearing modes. Phys. Fluids 29, 3230.Google Scholar
Dickinson, D., Roach, C. M., Saarelma, S., Scannell, R., Kirk, A. & Wilson, H. R. 2012 Kinetic instabilities that limit $\unicode[STIX]{x1D6FD}$ in the edge of a tokamak plasma: a picture of an H-mode pedestal. Phys. Rev. Lett. 108, 135002.Google Scholar
Drake, J. F., Antonsen, T. M. Jr., Hassam, A. B. & Gladd, N. T. 1983 Stabilization of the tearing mode in high-temperature plasma. Phys. Fluids 26, 2509.Google Scholar
Fitzpatrick, R. 1989 Linear stability of low mode number tearing modes in the banana collisionality regime. Phys. Plasmas 1, 2381.Google Scholar
Fitzpatrick, R., Hastie, R. J., Martin, T. J. & Roach, C. M. 1993 Stability of coupled tearing modes in tokamaks. Nucl. Fusion 33, 1533.Google Scholar
Furth, H. P., Killeen, J. & Rosenbluth, M. N. 1963 Finite-resistivity instabilities of a sheet pinch. Phys. Fluids 6, 459.CrossRefGoogle Scholar
Furth, H. P., Rutherford, P. H. & Selberg, H. 1973 Tearing mode in cylindrical plasmas. Phys. Fluids 16, 1054.Google Scholar
Gimblett, C. G., Hastie, R. J. & Helander, P. 2006 Model for current-driven edge-localized modes. Phys. Rev. Lett. 96, 035006.Google Scholar
Glasser, A. H., Greene, J. M. & Johnson, J. L. 1975 Resistive instabilities in general toroidal plasma configurations. Phys. Fluids 18, 875.CrossRefGoogle Scholar
Glasser, A. H., Greene, J. M. & Johnson, J. L. 1976 Resistive instabilities in a tokamak. Phys. Fluids 19, 567.Google Scholar
Greene, J. M. & Johnson, J. L. 1962 Stability criterion for arbitrary hydromagnetc equilibria. Phys. Fluids 5, 510.Google Scholar
Hatch, D. R., Kotschenreuther, M., Mahajan, S., Valanju, P., Jenko, F., Told, D., Görler, T. & Saarelma, S. 2016 Microtearing turbulence limiting the JET-ILW pedestal. Nucl. Fusion 56, 104003.Google Scholar
Hegna, C. C. & Callen, J. D. 1994 Stability of tearing modes in tokamak plasmas. Phys. Plasmas 1, 2308.Google Scholar
Hegna, C. C., Connor, J. W., Hastie, R. J. & Wilson, H. R. 1996 Toroidal coupling of ideal magnetohydrodynamic instabilities in tokamak plasmas. Phys. Plasmas 3, 584.Google Scholar
Huysmans, G. T. A. 2005 External kink (peeling) modes in x-point geometry. Plasma Phys. Control. Fusion 47, 2107.Google Scholar
Johnson, J. L. & Greene, J. M. 1967 Resistive interchanges and the negative v” criterion. Plasma Phys. 9, 611.Google Scholar
Mercier, C. 1960 A necessary condition for hydromagnetic stability of plasma with axial symmetry. Nucl. Fusion 1, 47.Google Scholar
Newcomb, W. A. 1960 Hydromagnetic stability of a diffuse linear pinch. Ann. Phys. NY 10, 232.Google Scholar
Nishimura, Y., Callen, J. D. & Hegna, C. C. 1998 Tearing mode analysis in tokamaks, revisited. Phys. Plasmas 5, 4292.Google Scholar
Pegoraro, F. & Schep, T. J. 1986 Theory of resistive modes in the ballooning representation. Plasma Phys. Control. Fusion 28, 647.Google Scholar
Porcelli, F. 1987 Viscous resistive magnetic reconnection. Phys. Fluids 30, 1734.Google Scholar
Snyder, P. B., Wilson, H. R., Ferron, J. R., Lao, L. L., Leonard, A. W., Mossessian, D., Murakami, M., Osborne, T. H., Turnbull, A. D. & Xu, X. Q. 2004 ELMs and constraints on the H-mode pedestal: peeling-ballooning stability calculation and comparison with experiment. Nucl. Fusion 44, 320.Google Scholar
Suydam, B. R. 1958 Stability of a linear pinch. In Proc. 2nd United Nations Conference PAUE, United Nations, Geneva, vol. 31, p. 354.Google Scholar
Webster, A. J. & Gimblett, C. G. 2009 Magnetohydrodynamic stability of a toroidal plasmas separatrix. Phys. Rev. Lett. 102, 035003.Google Scholar
Wilson, H. R., Connor, J. W., Field, A. R., Fielding, S. J., Miller, R. L., Lao, L. L., Ferron, J. R. & Turnbull, A. D. 1999 Ideal magnetohydrodynamic stability of the tokamak high-confinement-mode edge region. Phys. Plasmas 6, 1925.Google Scholar