2 Properties of the anisotropic diffusion equation

The ADE can be expressed as

(2.1)\begin{equation} \boldsymbol{\nabla} \boldsymbol{\cdot} \left(\boldsymbol{\kappa} \boldsymbol{\cdot} \boldsymbol{\nabla} T \right) = 0, \end{equation}

with the diffusion tensor

(2.2)\begin{equation} \boldsymbol{\kappa} = \kappa_{{\perp}}\boldsymbol{I} + (\kappa_{\|} - \kappa_{{\perp}}) \hat{\boldsymbol{b}} \hat{\boldsymbol{b}}, \end{equation}

where the perpendicular, $\kappa _{\perp }$, and parallel, $\kappa _{\|}$, diffusion coefficients are taken to be constants. We compute solutions to (2.1) in an annular volume $\varOmega$ bounded by two toroidal surfaces, $S_{-}$ and $S_{+}$ with constant Dirichlet boundary conditions $T_{-}$ and $T_{+}$.

2.1 Ellipticity

We note that $\boldsymbol {\kappa }$ is positive–definite for $\kappa _{\perp } > 0$ and $\kappa _{\|} \ge \kappa _{\perp }$ since

(2.3)\begin{equation} \boldsymbol{\nabla} T \boldsymbol{\cdot} \boldsymbol{\kappa} \boldsymbol{\cdot} \boldsymbol{\nabla} T = \kappa_{{\perp}} |\boldsymbol{\nabla} T|^2 + (\kappa_{\|}- \kappa_{{\perp}}) (\hat{\boldsymbol{b}} \boldsymbol{\cdot} \boldsymbol{\nabla} T )^2 > 0, \end{equation}

for all $\boldsymbol {\nabla } T \ne 0$. (While the assumption $\kappa _{\|} \ge \kappa _{\perp }$ is not necessary for ellipticity, it is sufficient according to (2.3) and not too restrictive since we are interested in strong anisotropy, $\kappa _{\|}\gg \kappa _{\perp }$.) We can furthermore see that the anisotropic diffusion operator is elliptic, as it can be expressed as

(2.4)\begin{equation} \sum_{i,j} \partial_i \left(\kappa_{ij} \partial_j T\right) = \sum_{i,j} \kappa_{ij} \partial_i \partial_j T + \left(\partial_i \kappa_{ij}\right) \partial_j T. \end{equation}

Ellipticity of the anisotropic diffusion operator (2.1) follows from the symmetric positive– definiteness of the tensor which multiplies the second derivatives, $\kappa _{ij}$. Ellipticity implies the existence of a maximum principle (Appendix A),

(2.5)\begin{equation} \max_{\bar{\varOmega}} T = \max_{\partial \varOmega} T. \end{equation}

The corresponding statement also holds for the minimum. In words, the maximum (minimum) temperature must be obtained on the boundary. Furthermore, no local minima or maxima can occur within $\varOmega$. As a further consequence of ellipticity, the topology of the isotherms is constrained (Appendix B). Specifically, an isotherm cannot enclose a net volume, but instead must link the domain toroidal and poloidally as the boundaries (with the possible addition of a handle).

2.2 Variational principle

The ADE can be obtained (Helander, Hudson & Paul Reference Helander, Hudson and Paul2021) from a variational principle involving the functional

(2.6)\begin{equation} \mathcal{W}[T] = \frac{1}{2}\int_{\varOmega}\,{\rm d}^3 x \, \boldsymbol{\nabla} T \boldsymbol{\cdot} \boldsymbol{\kappa} \boldsymbol{\cdot} \boldsymbol{\nabla} T. \end{equation}

The first variation of $\mathcal {W}$ with respect to $T$ is computed to be

(2.7)\begin{equation} \delta \mathcal{W}[\delta T] = \int_{\varOmega} \,{\rm d}^3 x \, \boldsymbol{\nabla} \delta T \boldsymbol{\cdot} \boldsymbol{\kappa} \boldsymbol{\cdot} \boldsymbol{\nabla} T={-} \int_{\varOmega} \,{\rm d}^3 x \, \delta T \boldsymbol{\nabla} \boldsymbol{\cdot} \left( \boldsymbol{\kappa} \boldsymbol{\cdot} \boldsymbol{\nabla} T \right), \end{equation}

upon application of the boundary condition $\delta T \rvert _{\partial \varOmega } = 0$. Thus stationary points of $\mathcal {W}$ with respect to $T$ correspond to solutions of (2.1). We can now compute the second variations as

(2.8)\begin{equation} \delta^2\mathcal{W}[\delta T,\delta T'] = \int_{\varOmega} \,{\rm d}^3 x \, \boldsymbol{\nabla} \delta T \boldsymbol{\cdot} \boldsymbol{\kappa} \boldsymbol{\cdot} \boldsymbol{\nabla} \delta T'. \end{equation}

We note that $\delta ^2 \mathcal {W}[\delta T,\delta T] \ge 0$ under the assumption that $\kappa _{\perp }<\kappa _{\|}$. Therefore, $\mathcal {W}$ is a convex function, implying that any stationary point of $\mathcal {W}$ is a global minimum. We will refer to the functional $\mathcal {W}$ that appears in this variational principle as the diffusion integral

(2.9)\begin{equation} D = \frac{1}{2} \int_{\varOmega}\, {\rm d}^3 x \, \left( \kappa_{\|}| \boldsymbol{\nabla}_\| T |^2 + \kappa_\perp |\boldsymbol{\nabla}_\perp T |^2 \right), \end{equation}

as it is related to the entropy production functional (Hameiri & Bhattacharjee Reference Hameiri and Bhattacharjee1987). Thus the temperature profile minimizes the entropy production subject to the boundary conditions.

We note that the total diffusion is related to the heat flux through the boundaries

(2.10)\begin{equation} Q ={-}\int_{S_+} \,{\rm d}^2 x \, \hat{\boldsymbol{n}} \boldsymbol{\cdot} \boldsymbol{q} ={-}\int_{S_-} \,{\rm d}^2 x \, \hat{\boldsymbol{n}} \boldsymbol{\cdot}\boldsymbol{q}, \end{equation}

through

(2.11)\begin{equation} Q = \frac{2 D[T]}{T_{+} - T_{-}}, \end{equation}

where $\boldsymbol {q} = -\kappa _{\|} \boldsymbol {\nabla }_{\|}T - \kappa _{\perp } \boldsymbol {\nabla }_{\perp } T$ and $\hat {\boldsymbol {n}}$ is the unit normal oriented in the direction from $S_-$ to $S_+$. This relation follows from

(2.12)\begin{align} 0 & = \int_\varOmega \,{\rm d}^3 x \, T \boldsymbol{\nabla} \boldsymbol{\cdot} {\boldsymbol q} ={-} \int_\varOmega \,{\rm d}^3 x \, {\boldsymbol q} \boldsymbol{\cdot} \boldsymbol{\nabla} T + \int_{S_+} \,{\rm d}^2 x \, T \hat{\boldsymbol{n}} \boldsymbol{\cdot} {\boldsymbol q} - \int_{S_-} \,{\rm d}^2 x \, T \hat{\boldsymbol{n}} \boldsymbol{\cdot} {\boldsymbol q}\nonumber\\ & = 2 D[T] - Q (T_+{-} T_-). \end{align}

Thus the total diffusion quantifies the heat flux required to support a given temperature gradient across a volume. In this work we will consider $T_+$ and $T_-$ to be fixed such that $D$ measures the required heat flux to support the prescribed temperature gradient. Since the perpendicular transport is relatively insensitive to the structure of the field, we can assess the impact of non-integrability on the transport through the local parallel diffusion

(2.13)\begin{equation} d_{\|} = \tfrac{1}{2}\kappa_{\|} |\boldsymbol{\nabla}_{\|}T|^2, \end{equation}

in comparison with the local perpendicular diffusion,

(2.14)\begin{equation} d_{{\perp}} = \tfrac{1}{2}\kappa_{{\perp}} |\boldsymbol{\nabla}_{{\perp}}T|^2. \end{equation}

3 The effective volume of parallel diffusion

Motivated by comparing $d_{\|}$ to $d_{\perp }$ in order to determine the impact of non-integrability on transport, we propose the following metric:

(3.1)\begin{equation} \mathcal{V}_{\text{PD}} = \frac{\int_{\varOmega} \,{\rm d}^3 x \, \varTheta \left(\kappa_{\|}|\boldsymbol{\nabla}_{\|}T|^2-\kappa_{{\perp}} |\boldsymbol{\nabla}_{{\perp}}T|^2\right)}{\int_{\varOmega} \,{\rm d}^3 x}, \end{equation}

where $\varTheta$ is the Heaviside step function. In words, $\mathcal {V}_{\text {PD}}$ measures the fraction of volume over which the local parallel transport is larger than the perpendicular transport. While one could imagine many ways in which to compare the parallel and perpendicular diffusion, the above definition is easy to interpret as it is a positive scalar in $[0,1]$.

While the definition of $\mathcal {V}_{\text {PD}}$ is not intrinsic to the field structure as it depends on $\kappa _{\perp }/\kappa _{\|}$, we can infer general trends in the limit of small $\kappa _{\perp }/\kappa _{\|}$. In the limit of an integrable magnetic field with sufficiently small $\kappa _{\perp }$, the isotherms will largely coincide with magnetic surfaces, $|\hat {\boldsymbol {b}} \boldsymbol {\cdot } \boldsymbol {\nabla } T|$ will be relatively small, and thus the local parallel transport will be negligible. Thus $\mathcal {V}_{\text {PD}}$ will be small in the limit of integrability. Since the topology of the isotherms is constrained, isotherms will not be able to completely align to the structure of the field in regions of non-integrability, and $|\hat {\boldsymbol {b}} \boldsymbol {\cdot } \boldsymbol {\nabla } T|$ will be consequently larger, as measured by $\mathcal {V}_{\text {PD}}$.

We can consider the metric $\mathcal {V}_{\text {PD}}$ to be a ‘diffusive’ analogue of a metric based on converse KAM theory. With a converse KAM approach, the non-existence of an invariant torus at a given point in phase space is determined from the rotation of nearby trajectories with respect to a chosen foliation. This allows one to determine regions of phase space over which invariant surfaces do not exist with the topology of the chosen foliation, i.e. a given class of tori (MacKay & Percival Reference MacKay and Percival1985; MacKay Reference MacKay2018). Similarly, by comparing the local parallel and perpendicular diffusion, we can assess whether a given region of physical space is ‘effectively non-integrable’ (indicated by enhanced parallel diffusion) with respect to its impact on the transport. This assessment similarly depends on the topology of the isotherms, which foliate the volume. Since the isotherms must link the domain as the boundaries do, the ADE similarly provides an assessment of the non-integrability with respect to a limited class of tori. For instance, an island chain will impact the parallel transport if the boundary surfaces are chosen to have the same topology as the unperturbed Hamiltonian, given that isotherms cannot have the same linking as the island. If instead the boundary surfaces have the same linking as the surfaces within the island chain, the parallel transport will not be impacted.

To assess this metric, we consider two cases for which analytic predictions for the parallel transport can be employed: diffusion across a single island chain (§ 5) and across a strongly chaotic layer (§ 6). In § 7 we then consider a model field for which such prediction is more challenging, namely a critically chaotic field with remnants of island chains and partial barriers.

4 Numerical methods

4.2 Model magnetic field

We consider a model Hamiltonian,

(4.1)\begin{equation} \chi(\psi,\theta,\zeta) = \frac{\iota' \psi^2}{2} + \sum_{m,n}\epsilon_{m,n} \psi (\psi - \bar{\psi}) \cos(m \theta - n \zeta), \end{equation}

which yields the magnetic field,

(4.2)\begin{equation} \boldsymbol{B} = \boldsymbol{\nabla} \psi \times \boldsymbol{\nabla} \theta - \boldsymbol{\nabla} \chi \times \boldsymbol{\nabla} \zeta, \end{equation}

where $\psi = \bar {\psi } \rho /\bar {\rho }$ and $(\rho,\theta,\zeta )$ forms an orthogonal coordinate system ($\hat {\boldsymbol {\rho }} \times \hat {\boldsymbol {\theta }} \boldsymbol {\cdot } \hat {\boldsymbol {\zeta }} = 1$) with $\rho \in [0,\bar {\rho }]$, $\theta \in [0,2{\rm \pi} )$, $\zeta \in [0,2{\rm \pi} )$, $|\boldsymbol {\nabla } \rho | = 1$, $|\boldsymbol {\nabla } \theta | = 1/\bar {\rho }$ and $|\boldsymbol {\nabla } \zeta | = 1/L_{\zeta }$. We will use the notation $\tilde {\epsilon }_{m,n}(\psi ) = \epsilon _{m,n}\psi (\psi -\bar {\psi })$ to denote the strength of the perturbation modes. In the limit $\epsilon _{m,n} \rightarrow 0$ we recover an integrable magnetic field lying on surfaces of constant $\rho$ with rotational transform $\iota = \iota ' \psi$ and constant shear $\iota '(\psi ) = \iota '$. The perturbation is chosen to vanish at the upper ($\rho = \bar {\rho }$) and lower ($\rho = 0$) boundaries such that there is no parallel flux into the volume. For all of the following numerical calculations, we will take $\bar {\psi } = \bar {\rho } = L_{\zeta } = \iota ' = 1$.

5 Transport across an island chain

In this section, we consider the relative perpendicular and parallel transport in the presence of a single island chain generated by a resonance $\iota = n/m$. We consider the case in which $\kappa _{\perp }$ remains small but non-zero such that perpendicular diffusion competes with parallel diffusion.

To analyse the behaviour of the transport in such a field, in Appendix C we perform a perturbation analysis in the smallness of the amplitude $\epsilon _{m,n}$. Due to the smallness of $\kappa _{\perp }$, outside of a small boundary layer of width

(5.1)\begin{equation} \mathcal{W}_{{c}} = \kappa_{{\perp}}^{1/4}\left(1 + \frac{n^2\bar{\rho}^2}{m^2 L_{\zeta}^2}\right)^{1/4}\left(\frac{L_{\zeta} \bar{\rho}}{m \iota' \bar{\psi}}\right)^{1/2}, \end{equation}

the perpendicular diffusion becomes unimportant. The solution is evaluated both inside and outside the boundary layer. We conclude that for both the inner and outer solution, the parallel and perpendicular diffusion can only be comparable within the island separatrices when $\mathcal {W}_{m,n}/\mathcal {W}_{{c}} \gtrsim 1$, where the island half-width is

(5.2)\begin{equation} \mathcal{W}_{m,n} = \frac{2\bar{\rho}}{\bar{\psi}}\sqrt{\left.\frac{\tilde{\epsilon}_{m,n}(\psi)}{\iota'(\psi)}\right\rvert} _{\psi = n/(m\iota')} . \end{equation}

In the limit that $\kappa _{\perp } \rightarrow 0$, $\mathcal {W}_{m,n}/\mathcal {W}_{{c}} \gg 1$, and the parallel transport can compete with perpendicular transport throughout the volume within the separatrices. In this sense, an island can be considered unimportant to the transport if

(5.3)\begin{equation} \mathcal{W}_{m,n} \ll \mathcal{W}_{{c}}. \end{equation}

In figure 1 we present the numerically computed values of the effective volume (3.1) for a single island chain created with resonance $\epsilon _{2,1}$ with varying amplitudes and values of $\kappa _{\perp }$. The critical value, $\epsilon _{2,1}^{\text {crit}} = \sqrt {\kappa _{\perp }}/2$, predicted for the parallel transport to compete with the perpendicular transport, is indicated by the vertical dashed line. We see that, as expected, when the island is smaller than the critical width $\mathcal {W}_{{c}}$, the computed effective volume vanishes. The point at which $\mathcal {V}_{\text {PD}}$ becomes positive is within an order of magnitude of $\epsilon _{2,1}^{\text {crit}}$. For large values of $\epsilon _{2,1}$, the effective volume approaches the expected scaling of the total volume enclosed by the separatrices, $\mathcal {W}_{2,1} \sim \sqrt {\epsilon _{2,1}}$, shown in orange. At non-zero values of the perpendicular diffusion, the scaling of $\mathcal {V}_{\text {PD}}$ is shallower than the scaling of the total island width since the entire island chain is not dominated by parallel diffusion.

Figure 1. The numerically computed value of the effective volume (3.1) is shown (blue) for the model field (4.2) with the displayed values of $\epsilon _{2,1}$ and $\kappa _{\perp }$. Also displayed is the scaling of the island width with $\epsilon _{2,1}$ (orange), and the critical value of $\epsilon _{2,1}$ for which the parallel transport may compete with perpendicular transport (black). The numerically computed value is predicted to approach the $\epsilon _{2,1}^{1/2}$ scaling in the limit $\kappa _{\perp } \rightarrow 0$; (a) $\kappa _{\perp } = 10^{-6}$,(b) $\kappa _{\perp } = 10^{-7}$, (c) $\kappa _{\perp } = 10^{-8}$.

In figure 2 we present the isotherms (red) computed with the displayed values of $\kappa _{\perp }$ on the Poincaré section of the model magnetic field (4.2) with $\epsilon _{2,1} = 0.008$. As $\kappa _{\perp }$ decreases, the isotherms conform more strongly to the structure of the magnetic field. We note the presence of isotherms which cut across the island chain at smaller values of $\kappa _{\perp }$ due to the topological restrictions. On the right we display the Poincaré section with colour highlighting the region over which the local parallel diffusion is dominating. The colour scale indicates the ratio of the local parallel to perpendicular transport. The parallel transport is strongest near the X-points, where the misalignment of the isotherms with the local field is the strongest. The parallel transport is relatively weak at the top and bottom of the island chain, where the isotherms are more strongly aligned with the magnetic field and the resonant field is weaker. As $\kappa _{\perp }$ decreases, we see the area over which the parallel transport dominates increase as well as its magnitude in comparison with the perpendicular transport.

Figure 2. The Poincaré section of the model field (4.2) with $\epsilon _{2,1} = 0.008$ is displayed with isotherms (red) on the left computed with the three values of $\kappa _{\perp }$ displayed. On the right, the Poincaré section is displayed with the ratio of the parallel (2.13) to perpendicular (2.14) diffusion. The colour scale is set to white where the parallel diffusion is smaller than the perpendicular diffusion; (a,b) $\kappa _{\perp } = 10^{-6}$, (c,d) $\kappa _{\perp } = 10^{-7}$, (e, f) $\kappa _{\perp } = 10^{-8}$.

6 Transport across a strongly chaotic layer

In this section we consider the relative impact of perpendicular and parallel transport in the presence of strong chaos. We will assume a quasilinear model (Chirikov Reference Chirikov1979; Hazeltine & Meiss Reference Hazeltine and Meiss2003; Lichtenberg & Lieberman Reference Lichtenberg and Lieberman2013) to obtain a parallel diffusion coefficient. We will assume that the parallel transport at a given location will be dominated by nearby resonances within an island width and that the phase of the resonances is decorrelated after one toroidal transit. We assume strong chaos such that the Chirikov island overlap parameter is large,

(6.1)\begin{equation} \frac{\mathcal{W}_{m,n} + \mathcal{W}_{m',n'}}{(\Delta \rho)_{m,n;m',n'}} \gg 1, \end{equation}

where $\mathcal {W}_{m,n}$ is the island half-width associated with resonances (5.2) and $(\Delta \rho )_{m,n;m',n'}$ is the distance between the resonances.

Under the assumption that the resonances are strongly overlapping but $\epsilon _{m,n}$ is sufficiently small that a perturbation analysis is appropriate, the motion in the field can be treated as a random walk

(6.2)\begin{equation} \frac{\partial T}{\partial \zeta} \sim D_{\|} \frac{\partial ^2T}{\partial \psi^2}, \end{equation}

with a quasi-linear diffusion coefficient given by

(6.3)\begin{equation} D_{\|} \approx \frac{\rm \pi}{2} \sum_{m} m^2 \left(\sum_n \varTheta(\mathcal{W}_{m,n}/2 -|\rho-\rho_{{r}}|) \tilde{\epsilon}_{m,n}(\psi_0) \right)^2. \end{equation}

See Appendix D for details. The local parallel transport then scales as

(6.4)\begin{equation} |\boldsymbol{\nabla}_{\|} T|^2 \sim D_{\|}^2 |\boldsymbol{\nabla}_{{\perp}} ^2T|^2 \frac{\bar{\rho}^4}{\bar{\psi}^4 L_{\zeta}^2}. \end{equation}

We can estimate $|\boldsymbol {\nabla }_{\perp }^2 T| \sim |\boldsymbol {\nabla }_{\perp }T|/L_{\perp }$ for some perpendicular length scale $L_{\perp }$. Given the analysis in § 5 and Appendix C, we can note several relevant length scales. Within the boundary layer, $L_{\perp } \sim \mathcal {W}_{{c}}$. Outside a width $\mathcal {W}_{m,n}$ of a given resonance and outside the boundary layer, $L_{\perp } \sim \bar {\rho }$. Within a width $\mathcal {W}_{m,n}$ and outside the boundary layer, $L_{\perp } \sim |\rho -\rho _{{r}}|$. Given that the resonances are strongly overlapping, the perpendicular temperature gradient will be dominated by nearby resonances. Thus we can estimate $L_{\perp } \lesssim \mathcal {W}_{m,n}$ and,

(6.5)\begin{equation} |\boldsymbol{\nabla}_{\|} T|^2\gtrsim D_{\|}^2 \frac{|\boldsymbol{\nabla}_{{\perp}} T|^2}{\mathcal{W}_{m,n}^2} \frac{\bar{\rho}^4}{\bar{\psi}^4L_{\zeta}^2}. \end{equation}

We conclude that the parallel transport may compete with the perpendicular transport when,

(6.6)\begin{equation} \frac{D_{\|}^2 \bar{\rho}^4}{\mathcal{W}_{m,n}^2 \bar{\psi}^4L_{\zeta}^2} \sim \kappa_{{\perp}}. \end{equation}

To test the expected scaling, we consider the model field (4.2) with $m = 12$ and $n = [2,\ldots,10]$. We evaluate three strongly chaotic fields with amplitudes $\epsilon _{m,n}$ chosen such that the Chirikov parameter (6.1) is 2, $2 \sqrt {2}$ and 4. We compute solutions to the ADE for a range of $\kappa _{\perp }$. The effective volume (3.1) is shown in figure 3 along with the critical value of $\kappa _{\perp }$ (vertical dashed) at which the parallel diffusion may compete (6.6). We note that the point at which the effective volume becomes non-zero is within an order of magnitude of that expected from quasilinear theory for each value of the overlap parameter. As $\kappa _{\perp }$ decreases, the effective volume computed from the three model fields approach each other.

Figure 3. The numerically computed value of the effective volume (3.1) is shown for the model field (4.2) with $m=12$ resonances with amplitudes chosen to provide strong island overlap (values of the overlap parameter (6.1) are indicated). Also displayed is the critical value of $\kappa _{\perp }$ for which the parallel transport is predicted to compete with perpendicular transport (vertical dashed).

In figure 4 we display the Poincaré section with isotherms (left) and with the parallel-transport-dominated regions highlighted in colour (right). Although the field is strongly chaotic, we note several isotherms near the boundaries which adapt to the $m = 12$ structure of the residual island chains. For $\kappa _{\perp } = 10^{-4}$, we note only a very small volume over which the parallel transport dominates. As we decrease the value of $\kappa _{\perp }$, we see that the parallel transport dominates over much of the volume contained within the outermost island separatrices and the temperature gradient across the chaotic layer is significantly reduced. We note an $n = 6$ variation in the structure of the parallel diffusion due to the poloidal dependence of the radial field (4.2). Near the edges of the chaotic layer, the isotherms become strongly shaped in comparison with the relatively flat isotherms in the chaotic layer. Given that the magnitude of the radial field is smaller than the toroidal field by the amplitude of the perturbation, this feature enhances the parallel diffusion near the boundaries of the chaotic layer.

Figure 4. The Poincaré section of the model field (4.2) with $m = 12$ resonances chosen to provide an overlap parameter (6.1) of 4 is displayed with isotherms (red) on the left for the three values of $\kappa _{\perp }$ displayed. On the right, the Poincaré section is displayed with the ratio of the parallel (2.13) to perpendicular (2.14) diffusion. The colour scale is set to white where the parallel diffusion is smaller than the perpendicular diffusion; (a,b) $\kappa _{\perp } = 10^{-4}$, (c,d) $\kappa _{\perp } = 10^{-5}$, (e, f) $\kappa _{\perp } = 10^{-6}$.

7 Transport across a critically chaotic layer

Although simple analytic scalings can be derived in the case of a single island chain or strong chaos, in reality a magnetic field may possess a complicated mixture of KAM surfaces, partial barriers, chaos and island chains. In this case, we can use the metric defined in § 3 to assess the relative impact of non-integrability on the transport through the volume of parallel diffusion.

As a point of comparison, we will consider four marginally chaotic magnetic fields generated by resonances with different poloidal mode numbers. The amplitude of the perturbation will be chosen to achieve critical island overlap with a comparable volume enclosed by the outermost separatrices. For the first field, we choose $n/m = [1/4, 2/4, 3/4]$. Since the resonances are equally spaced with $\Delta \psi = 1/4$, we choose the perturbation amplitudes $\epsilon _{m,n}$ such that the associated half-island widths (5.2) are $1/8$. We similarly consider fields with $n/m = [2/12,\ldots,10/12]$, $n/m =[3/20,\ldots,17/20]$, and $n/m = [5/36,\ldots,31/36]$ with critical island overlap. The separatrices associated with the outermost island chains for each of these model fields are located at $\rho =1/8$ and $\rho = 7/8$ such that they have a comparable ‘volume’ of chaos.

If the remnants of the island chains are dominating the transport, we anticipate the parallel transport to dominate when

(7.1)\begin{equation} \epsilon_{m,n}^2 m^2 \sim \kappa_{{\perp}}, \end{equation}

from (5.3), while if strong chaos is dominating the transport

(7.2)\begin{equation} \epsilon_{m,n}^3 m^4 \sim \kappa_{{\perp}}, \end{equation}

from (6.6). Since we have chosen $\epsilon _{m,n} \sim 1/m^2$ to achieve critical overlap, we expect the parallel transport to be relatively more impactful when $m$ is small. In figure 5 we compare the effective volume between the model fields. Since the volume contained within the outermost separatrices differs slightly between the two model fields, we compute the effective volume only for $\rho \in [0.25, 0.75]$. As anticipated, the $m = 4$ resonances produce a larger effective volume of parallel diffusion, and the relative difference between the model fields increases with increasing $\kappa _{\perp }$. To quantify the effect of the field structure on temperature flattening, we also compare the total temperature differential between $\rho = 0.75$ and $\rho = 0.25$. At large values of $\kappa _{\perp }$, each of the model fields can support a significant temperature gradient, although it is reduced for the $m =4$ field.

Figure 5. (a) The effective volume (3.1) is computed for the model field (4.2) with the displayed mode number perturbations with amplitudes $\epsilon _{m,n}$ chosen for critical island overlap. The effective volume is computed over a subset of the entire volume, $\rho \in [0.25, 0.75]$. (b) The total temperature difference between $\rho = 0.75$ and $\rho = 0.25$ upon averaging over the angles.

As can be seen from figure 6, each of the fields are visually chaotic with secondary island chains as well as stochastic regions. We furthermore note that the magnitude of the parallel transport, in addition to the volume over which it dominates, is considerably larger in the $m = 4$ field. For the $m = 36$ perturbation, we note that the temperature gradient is not markedly reduced in the chaotic layer when $\kappa _{\perp } = 10^{-6}$, although all of the flux surfaces are destroyed. Thus while all of the fields would be classified as very far from integrability by traditional measures such as the Chirikov overlap criterion or converse KAM theory, the impact of the magnetic field structure on the overall transport at non-zero $\kappa _{\perp }$ is remarkably distinct.

Figure 6. The ADE is solved for the model field (4.2) with critically overlapping resonances of $m = 4$, $m = 12$ and $m=36$ with $\kappa _{\perp } = 10^{-6}$. On the left are the isotherms (red), and on the right is the ratio of the parallel (2.13) to perpendicular (2.14) diffusion. The colour scale is set to white where the parallel diffusion is smaller than the perpendicular diffusion; (a,b) $m=4$, (c,d) $m = 12$, (e, f) $m = 36$.