Regardless of how small non-ideal effects may be, phenomena associated with changes in magnetic field line connections are frequently observed to occur on an Alfvénic time scale. Since it is mathematically impossible for magnetic field line connections to change when non-ideal effects are identically zero, an ideal evolution must naturally lead to states of unbound sensitivity to non-ideal effects. That such an evolution is natural is demonstrated by the use of Lagrangian coordinates based on the flow velocity of the magnetic field lines. The Lagrangian representation of an evolving magnetic field is highly constrained when neither the magnetic field strength nor the forces exerted by the magnetic field increase exponentially with time. The development of a state of fast reconnection consistent with these constraints (1) requires a three-dimensional evolution, (2) has an exponentially increasing sensitivity to non-ideal effects, and (3) has a parallel current density, which lies in exponentially thinning but exponentially widening ribbons, with a magnitude that is limited to a slow growth. The implication is that exponential growth in sensitivity is the cause of fast magnetic reconnection when non-ideal effects are sufficiently small. The growth of the non-ideal effect of the resistivity multiplied by the parallel current density is far too slow to be competitive.

## I. INTRODUCTION

Laboratory, space, solar, and astrophysical magnetic fields are commonly embedded in highly conducting, near-ideal, plasmas. In an ideal plasma, magnetic field lines do not change their connections or other topological properties, and the magnetic field evolves as

with $u\u2192$ being the velocity of the magnetic field lines,^{1} which need not be the velocity of the plasma $v\u2192$. As will be shown, Eq. (1) naturally leads to states of unbound sensitivity to non-ideal effects. This unbound sensitivity eventually produces a fast change in magnetic field line connections when non-ideal effects are sufficiently small but non-zero.

Equation (1) for an ideal evolution can be solved exactly for both the magnetic field, $B\u2192(x\u2192,t)$, and the current density, $j\u2192(x\u2192,t)=\u2207\u2192\xd7B\u2192/\mu 0$ in terms of the initial magnetic field $B\u21920(x\u21920)$ and the flow velocity of the magnetic field lines $u\u2192$. The solution is based on Lagrangian coordinates, $x\u21920$, which have a position vector $x\u2192(x\u21920,t)$ that is defined by

The position vector in ordinary Cartesian coordinates is $x\u2192=x(x0,y0,z0,t)x\u0302+y(x0,y0,z0,t)y\u0302+z(x0,y0,z0,t)z\u0302$.

How can the introduction of an unknown velocity of the magnetic field lines $u\u2192(x\u2192,t)$ provide important information on the evolution of magnetic fields? The answer is that the mathematical properties of Lagrangian coordinates allow the constraints of qualitative information to be included in the analysis.

The most important property of Lagrangian coordinates is exponentiation. What is meant by this? The Jacobian matrix of Lagrangian coordinates is defined as

The second matrix expression for $\u2202x\u2192/\u2202x\u21920$ is the Singular Value Decomposition (SVD) of the first. Any three-by-three matrix with real coefficients can be written in this form, where $U\u2194$ and $V\u2194$ are orthogonal matrices, $U\u2194\u2020\xb7U\u2194=1\u2194$. The coefficients Λ_{u} ≥ Λ_{m} ≥ Λ_{s} ≥ 0 of a singular value decomposition are called singular values and are positive real numbers. For all but exceptional flow velocities $u\u2192(x\u2192,t)$, the largest singular value Λ_{u} increases exponentially with time,

is called the Lyapunov exponent in the theory of dynamical systems. The singular value Λ_{s} generally decreases exponentially with time, $\Lambda s=e\u2212\sigma s$, while Λ_{m} has a much weaker evolution, sometimes a power law of time. The time dependence of the singular values is discussed in Sec. III.

It should be emphasized that the velocity $u\u2192(x\u2192,t)$ must have a very special form to avoid exponentiation, not the other way around.

The ideal approximation can become inadequate for two independent reasons: (1) The non-ideal part of the electric field, typically $\eta j||$, can become too large, where *η* is the plasma resistivity and $j||\u2261j\u2192\xb7B\u2192/B$. (2) The exponentiation of the separation between neighboring magnetic field lines with the distance along the lines amplifies an arbitrarily small non-ideal effect exponentially. The trajectory of a magnetic field line can be eventually affected on the scale of the system by an arbitrarily small non-ideal effect.

Appendix A shows the importance of exponentiation more directly. While non-ideal effects on the magnetic field remain small, the direction of the magnetic field lines is perturbed by the exponentially large factor Λ_{u} times the strength of the non-ideal part of the electric field.

Without constraints on $u\u2192$, the magnetic field strength would increase as Λ_{u}, which increases exponentially in time, and the force exerted by the field would increase as 1/Λ_{s}, which also increases exponentially in time. If and only if the evolution is in three-dimensional space, constraints can be imposed to eliminate these exponential increases. When these constraints are applied, the increase in the non-ideal effect $\eta j||/B$ is proportional to the increase in $\Lambda m2$, which is slow. But, the sensitivity to nonideal effects due to the separation of neighboring magnetic field lines increases as Λ_{u}.

The implication is that the effect on reconnection of the increase in the nonideal part of the electric field becomes subdominant to the effect of exponentiation as the nonideal effects become small. As will be discussed, singular current densities can arise in the infinite-time limit, but even when singular currents are predicted, as in the Parker conjecture,^{2,3} their effect on reconnection is subdominant to the effect of exponentiation when nonideal effects are sufficiently small.

A related argument on the subdominance of current singularities to exponentiation for magnetic reconnection was given in Ref. 4. There it was noted that the parallel current density $j||$ required to achieve a certain number of exponentiations in the separation between neighboring magnetic field lines scales linearly in the number of exponentiations.

Although the amplitude of $j||/B$ can increase only slowly, the gradient of $j||/B$ across the magnetic field lines increases exponentially in one direction and decreases exponentially in the other. That is, $j||/B$ lies in increasingly thin but wide ribbons. Although the thinning and widening of the current sheet and fast magnetic reconnection are correlated, both are due to large scale properties of the magnetic field evolution, not a local cause-and-effect relationship. The natural formation of current sheets in an evolving magnetic field should not be confused with the existence of Harris sheets,^{5} which are not consistent with three-dimensional plasma states.^{6} Nevertheless, a Harris sheet is a standard initial condition for two-dimensional studies of magnetic reconnection.

The requirement for a three-dimensional flow for developing states of fast magnetic reconnection may at first be surprising. Streamlines of a flow generally separate exponentially even when a time dependent velocity $u\u2192$ causes only a two-dimensional evolution. However, in a two-dimensional evolution, the only way to exponentially enhance reconnection is by an exponential increase in the magnetic field strength. When the evolution is two-dimensional, say in the *x* and *y* coordinates, the matrix

which has only the Λ_{u} and Λ_{s} singular values. In a three-dimensional evolution, the magnetic field can be proportional to the Λ_{m} singular value, but that singular value is missing for a two-dimensional evolution. The implication is that in two-dimensions, the only way to make a exponentially large change in the importance of nonideal terms is to exponentiate the magnetic field strength. Longcope and Strauss^{8} used the two-dimensional Jacobian matrix in their 1994 discussion of the formation of current layers.

Once magnetic field line connections are changed, large forces arise and must be relaxed. The reason for the large forces is each of the two parts of a newly connected field line will initially have a different parallel current. This implies a strong gradient in $j||/B$ with the distance along the line. When the Debye length is small, the current density must be divergence free,^{9} and $\u2207\u2192\xb7j\u2192=0$ can be written as

The Lorentz force, $f\u2192L$, is the force of the magnetic field on the plasma. The sudden change in the force on the plasma is balanced by the plasma viscosity, $f\u2192L=\u2207\u2192\xb7P\u2194$, or by the plasma inertia $f\u2192L=\rho (\u2202v\u2192/\u2202t+v\u2192\xb7\u2207\u2192v\u2192)$. The inertial term implies a relaxation by Alfvén waves. An Alfvénic limit on the speed of reconnection effects is not surprising—that is the fastest speed at which magnetic fields can transmit information either along or across magnetic field lines in the standard MHD approximation. Nevertheless, the time required for an Alfvénic relaxation is subtle because Alfvén waves propagating along exponentially separating field lines have enhanced damping.^{10,11}

Magnetic field lines are commonly observed to change their connections on a time scale consistent with the Alfvén speed *V _{A}*. A review of the observations and theory of reconnection proceeding at 0.1

*V*has been given by Cassak

_{A}*et al.*

^{12}Reconnection at a rate closely associated with Alfvénic rather than nonideal effects is called fast magnetic reconnection. Although Eq. (1) for an ideal evolution is inconsistent with magnetic field lines changing their topology, the equation predicts an exponentially increasing sensitivity to nonideal effects, Appendix A, which leads to fast magnetic reconnection.

In the limit as nonideal effects are very small, a fast reconnection event will occur after an adequate time, called a trigger time, for the ideal evolution to produce a sufficiently large Λ_{u}.

Assuming that the resistivity *η* is the most important nonideal term, the degree of exponentiation required for a reconnection trigger is $\Lambda u\u223c\u2111$. The dimensionless coefficient that measures the closeness to ideality, a Gothic “I,” is

is the time required for resistive diffusion across of the magnetic field lines on the scale *a* of the region that undergoes reconnection, and $\tau ev\u22611/|\u2207\u2192u\u2192|$ is the characteristic time scale for the magnetic evolution. When $\u2111\u2009\u2272\u20091$, resistive diffusion is so rapid that the difference between two and three-dimensional reconnection is of limited importance. As $\u2111$ becomes larger, the physics of reconnection becomes ever more sensitive to any breaking of the continuous symmetry that is assumed in two-dimensional models. The derivations of this paper are in the limit as $\u2111\u2192\u221e$, and important work remains to be done in how the validity of two dimensional models breaks down as $\u2111$ increases.

Remarkably, the importance of exponentiation has largely escaped notice in the literature on magnetic reconnection. For example, major recent reviews^{13,14} focused on two-dimensional plasmoid models of magnetic reconnection, which are descendants of the 1957 models of Sweet and Parker.^{15,16} An explanation for the speed of magnetic reconnection in a two dimensional plasmoid model based on a Harris current sheet warranted a 2017 Physical Review Letter.^{17}

When nonideal effects are small, the mathematical properties of an evolving magnetic field are fundamentally different between two and three dimensions. These differences bring the relevance of two-dimensional plasmoid models to three-dimensional problems into question.

Plasmoids are analogous to magnetic islands in topologically toroidal plasmas. Islands arise when perfect toroidal magnetic surfaces are perturbed. But islands are highly localized in toroidal plasmas; they only split magnetic surfaces that are rational surfaces, surfaces on which magnetic field lines close on themselves. Otherwise, Alfvén waves spread the effect of a perturbation over the volume of space covered by a single field line, and an island is not formed. How this spreading is consistent with plasmoid formation in three-dimensional space remains to be explained.

Delta function current densities are mathematically required when a rational magnetic surface in an ideal, steady-state, toroidal plasma is resonantly perturbed.^{18} Nevertheless, as shown by Hahm and Kulsrud,^{19} the current density in the vicinity of the rational surface increases only linearly in time after a resonant perturbation is applied. This arises from the time required for a shear-Alfvén wave propagating along the magnetic field to cover a near-rational surface and adequately sample its topological properties. The Hahm and Kulsrud time is proportional to $1/(N\u2212\iota M)$, where the magnetic field lines on the rational surface close on themselves after *M* toroidal and *N* poloidal transits. The rotational transform *ι* or twist of the magnetic field lines at the rational surface is *N*/*M*. The applicability of plasmoid models to reconnection in three-dimensional plasmas cannot be understood unless the analogue of the Hahm-Kulsrud time for the bounding surface of a plasmoid is obtained. Appendix B gives another example of a linear increase in the current density as a magnetic field is evolved ideally but slowly compared to the Alfvén speed.

What will be studied in this paper is the nature of magnetic reconnection in the limit as nonideal effects go to zero, the limit as $1/\u2111$ goes to zero. Turbulence produces three-dimensionality^{20} and can give the exponential sensitivity required for fast reconnection as nonideal effects become arbitrarily small, but turbulence is not required. Indeed, three-dimensionality on a large spatial scale causes the reconnection to occur on a large scale. Fast magnetic reconnection works much like stirring, Sec. VII. Large-scale stirring is a more effective way to mix a can of paint than small-scale stirring though both can cause mixing.

The sensitivity of the ideal constraint on magnetic evolution near places where current density becomes large has been discussed by a number of authors, including Low^{21} and Dewar *et al.*^{22} Here, it will be shown that in a highly ideal evolution the sensitivity is given by the magnitude of Λ_{u} and not directly by the current density.

A motivation for writing this paper was the development of an understanding of the rapid loss, $\u22721\u2009ms$, of magnetic surfaces that is commonly observed during the thermal quench phase of tokamak disruptions. For example, in JET^{23} the growth of resonant magnetic perturbations was observed to occur over 100's of millisecond before the current profile suddenly broadened reducing the internal inductance *ℓ _{i}* by approximately a factor of two. The relation between these observations and fast magnetic reconnection is discussed in Refs. 7, 24, and 25. In JET, the growth of islands over 100's of millisecond is consistent with the resistive opening of islands, but the sudden and large change in the current profile is not.

Section II derives the conditions for an ideal evolution. Section III derives the properties of Lagrangian coordinates required to determine the properties of ideally evolving magnetic fields. Section IV obtains the magnetic field evolution in Lagrangian coordinates. Section V gives expressions for the current density and the Lorentz, $j\u2192\xd7B\u2192$, force. Section VI obtains the separation of the magnetic field lines. Section VII discusses applications of Lagrangian coordinates in other areas of classical physics. Section VIII is a discussion, which focuses on the distinction between standard two-dimensional models of reconnection and the theory developed here. There are two appendixes. Appendix A derives the exponentially increasing departure of a magnetic field from its ideal form under the assumption that the departure is small. Appendix B derives the current density in the early stage of an ideal evolution when the system is driven by a slow flow in the perfectly conducting wall.

## II. CONDITIONS FOR AN IDEAL EVOLUTION

### A. Ohm's law

Equation (1) for the ideal evolution of a magnetic field holds in regions in which a magnetic field has nulls, $B\u2192=0$ points, when the plasma moving with a velocity $v\u2192$ obeys an Ohm's law $E\u2192+v\u2192\xd7B\u2192=0$. Then, the field line velocity is the plasma velocity $u\u2192\u2261v\u2192$.

The ideal evolution equation has a more general validity in regions of space in which there are no nulls of $B\u2192$. The magnetic field line velocity $u\u2192$ can be obtained from the generalized Ohm's law

which implies that the velocity of the fluid mass $v\u2192$ perpendicular to magnetic field lines that move with a velocity $u\u2192\u22a5$ satisfies

where Φ is an arbitrary single-valued function of position. When the inertial force $\rho dv\u2192/dt$ is negligible, the plasma velocity $v\u2192$ does not directly enter the evolution of the magnetic field. Nevertheless, the plasma flow can have a peculiar form relative to that of the field because of what are called Hall terms in the generalized Ohm's law.

The only physical effects that break the ideal evolution of the magnetic field are on the right hand side of the equation

The potential Φ can be chosen to minimize the field-line constant $Eni$ with

The integration limits on the two integrals are the same with $d\u2113\u2192$ the differential vector distance along a magnetic field line and *dℓ* the differential scalar distance. The integration limits can be plus and minus *L* → *∞* or from one perfectly conducting boundary to another. The ideality of the evolution can be measured by $\u2111=uB/Eni$.

### B. The Clebsch representation

A magnetic field can always be given in the Clebsch representation

A particular magnetic field line is specified by constant values of *α* and *β*. Stern^{26} has reviewed the history of the Clebsch representation.

The Clebsch potentials *α* and *β* can be used with the distance along the magnetic field lines *ℓ* to form a spatial coordinate system, which means that the three Cartesian coordinates (*x*, *y*, *z*) with $x\u2192=xx\u0302+yy\u0302+zz\u0302$ are given as functions of (*α*, *β*, *ℓ*). Clebsch coordinates are subtle when going from Cartesian space to Clebsch space, but not the other way around. When *α* and *β* are held fixed, $x\u2192(\alpha ,\beta ,\u2113)$ gives every point in Cartesian coordinates that a particular magnetic field line can reach as *ℓ* is varied. It is this property that makes the Clebsch variables (*α*, *β*, *ℓ*) useful for determining the deviation of a magnetic field from its ideal form due to a nonzero $Eni$, Appendix A.

Faraday's law, $\u2202B\u2192/\u2202t=\u2212\u2207\u2192\xd7E\u2192$, can be written using Eq. (14) for the electric field as $\u2202B\u2192/\u2202t=\u2207\u2192\xd7(u\u2192\xd7B\u2192\u2212EniB\u2192/B)$. Since $u\u2192\xd7B\u2192=(u\u2192\xb7\u2207\u2192\beta )\u2207\u2192\alpha \u2212(u\u2192\xb7\u2207\u2192\alpha )\u2207\u2192\beta $

Since $B\u2192\xb7\u2207\u2192\u2113=B$, the component of this equation along the magnetic field is

The evolution equations for *α* and *β* are

A coordinate transformation to Lagrangian coordinates gives the equations

The arbitrary function *g _{a}*(

*α*,

*β*,

*t*) just couples

*α*and

*β*but that does not change the identification of a magnetic field line with fixed

*α*and

*β*as the line is carried by the flow;

*g*can be chosen arbitrarily. The term $Eni\u2113$ in

_{a}*g*breaks the magnetic field lines from the flow and can change their topology. When $Eni=0$, the Clebsch coordinates are functions of the Lagrangian coordinates alone $\alpha (x\u21920)$ and $\beta (x\u21920)$.

A nontrivial example of *α* and *β* either satisfying the required equations for an ideal evolution, or not, is given by a toroidal plasma with magnetic surfaces. The toroidal magnetic flux enclosed by a surface, *ψ _{t}*, the poloidal angle,

*θ*, and the toroidal angle,

*φ*, can be chosen so

*α*=

*ψ*and $2\pi \beta =\theta \u2212\iota \phi $, where $\iota (\psi t,t)$ is called the rotational transform. When the toroidal field is strong

_{t}*ℓ*=

*R*

_{0}

*φ*and one can assume $(\u2202\u2113/\u2202t)x\u2192=0$. The nonideal part of $g=R0Eni(\psi t,t)\phi $. The evolution equations for

*α*and

*β*, then imply $(\u2202\psi t/\u2202t)x\u2192=\u2212u\u2192\xb7\u2207\u2192\psi t$ and $(\u2202\theta /\u2202t)x\u2192=\u2212u\u2192\xb7\u2207\u2192\theta $. When $(\u2202\iota /\u2202t)\psi t\u22600$, the evolution equation for

*β*implies $\phi (\u2202\iota /\u2202t)\psi t=2\pi \u2202g/\u2202\psi t$, or $(\u2202\iota /\u2202t)\psi t=2\pi R0(\u2202Eni/\u2202\psi t)$. The loop voltage is $V\u2113=2\pi R0Eni$, so $(\u2202\iota /\u2202t)\psi t=\u2202V\u2113/\u2202\psi t$, which is a well-known result in the physics of toroidal plasmas.

^{9}

## III. LAGRANGIAN COORDINATES

Lagrangian coordinates $x\u21920$, which were defined in Eq. (2), have a number of general properties that are required to understand their implications for magnetic evolution. These properties are derived in this section.

### A. The Jacobian matrix

Equation (4) for the SVD decomposition of the three-by-three matrix $\u2202x\u2192/\u2202x\u21920$ is equivalent to

The left eigenvectors which are defined by the orthogonal matrix $U\u2194$ obey

The right eigenvectors, which are defined by the orthogonal matrix $V\u2194$, obey analogous relations

The Jacobian of Lagrangian coordinates which is the determinant of the matrix $\u2202x\u2192/\u2202x\u21920$ obeys

The time derivative of the Jacobian can be found by considering the evolution of an arbitrary volume defined in fixed Lagrangian coordinates. Then

The component of the magnetic field line velocity $u\u2192$ in the direction along $B\u2192$ can be defined freely. Letting $b\u0302=B\u2192/B$, three attractive choices are: (1) $b\u0302\xb7u\u2192=0$, which is the simplest and the choice used in this paper, (2) $u\u2192\xb7\u2207\u2192\u2113=0$, which makes $u||=\u2212u\u2192\u22a5\xb7\u2207\u2192\u2113$, and (3) $b\u0302\xb7\u2207\u2192(b\u0302\xb7u\u2192)=u\u2192\xb7(b\u0302\xb7\u2207\u2192b\u0302)$, which makes *JB* = *B*_{0}, The requirement for the *JB* = *B*_{0} condition was incorrectly stated in Ref. 7 as $b\u0302\xb7u\u2192=0$, which is only true where the curvature of the magnetic field lines $b\u0302\xb7\u2207\u2192b\u0302$ vanishes.

### B. Evolution of the Jacobian matrix

To clarify the meaning of complicated dot products, Greek letters $\alpha ,\beta ,\gamma ,\u2026$ will be used to denote ordinary Cartesian coordinates $x\u2192$ and Latin letters $i,j,k,\u2026$ for Lagrangian coordinates $x\u21920$. Equation (2) for the definition of Lagrangian coordinates implies

In matrix notation

### C. The metric tensor

The exponential increase in Λ_{u} is made plausible by the expression for the evolution of the distance $\delta \u2192$ between neighboring (infinitesimally separated) streamlines. As $|\delta \u2192|\u21920$

Equation (37) is a linear equation for $\delta \u2192$, which generally has an exponentially increasing solution.

The distance between neighboring stream lines is

The symmetric tensor that appears in Eq. (39) is the metric tensor,

### D. Evolution of the metric tensor

The time derivative of the metric tensor can be written as

### E. Evolution of $\Lambda u,m,s2,\u2009u\u0302,\u2009m\u0302$, and $s\u0302$

The evolution of the basis vectors of Lagrangian coordinates was derived in 1987 by Goldhirsch *et al.*^{27} Contributions were also made by Tang and Boozer,^{28,29} by Thiffeault,^{31} and by Theffeault and Boozer.^{32} These techniques were applied to the evolution of the magnetic field in a dynamo in 2000 by Tang and Boozer^{30} and in 2003 by Thiffeault and Boozer.^{32}

To determine the time dependence of the $u\u0302$ basis vector holding Lagrangian coordinates constant, calculate

Using Eq. (41) for $g\u2194$ expanded in its basis vectors

Since $u\u0302\xb7u\u0302=1,\u2009u\u0302\xb7\u2202u\u0302/\u2202t=0$, the dot products of the unit vectors $u\u0302,\u2009m\u0302$, and $s\u0302$ with Eq. (45) imply

When the dot products of $G\u2194+G\u2194\u2020$ are not increasing as rapidly as Λ_{u} or 1/Λ_{s}, the unit vector $u\u0302$ has no further evolution once $\Lambda u\u226b\Lambda m\u226b\Lambda s$. An exponential increase in the $G\u2194+G\u2194\u2020$ dot products requires an exponential increase in $|\u2207\u2192u\u2192|2\u2261\u2211\alpha \beta (\u2202u\alpha /\u2202x\beta )2$.

The evolution of $m\u0302$ and $s\u0302$ can be found by appropriate changes in the letters that appear in the formula for the evolution of $u\u0302$. For example,

The orthogonality of the unit vectors $u\u0302\xb7m\u0302=0$ implies $u\u0302\xb7\u2202m\u0302/\u2202t=\u2212m\u0302\xb7\u2202u\u0302/\u2202t$

The evolution of the large singular value, Λ_{u}, is determined by Eq. (47). The properties of this evolution can be understood by noting that $G\u2194+G\u2194\u2020$ is a Hermitian matrix so it can be diagonalized with real eigenvalues

since $Z\u2194$ is an orthogonal matrix. The three eigenvalues are rates, with units of one over time, and ordered so *ν*_{+} ≥ *ν*_{0} ≥ *ν*_{−}. For a divergence-free flow they must sum to zero, *ν*_{+} + *ν*_{0} + *ν*_{−} = 0.

The unit vector $U\u0302$ rotates, Sec. III F, in such a way to come into alignment with $Z\u0302+$; at any instant

where the effective rate of growth of lnΛ_{u} satisfies *ν _{ef}* ≤

*ν*

_{+}. The infinite time Lyapunov exponent is defined as

Let $\nu ef=\lambda \u221e+\nu \u0303ef$, then when $\u222b0\u221e\nu \u0303efdt=ln\u2009A$ the largest singular value has the form

This form is equivalent to Eq. (1.11) of Goldhirsch *et al.*^{27} though other forms are possible in which the amplitude *A* has a sufficiently weak time dependences that (ln*A*)/*t* goes to zero as *t* goes to infinity. Reference 28 has a detailed discussion of these issues for two-dimensional divergence-free flows.

### F. Evolution of $U\u0302,\u2009M\u0302$, and $S\u0302$

Equation (35) for the time derivative of the Jacobian matrix can be used with Eq. (26), which gives the representation of the Jacobian matrix, to obtain the time derivatives of the $U\u0302,\u2009M\u0302$, and $S\u0302$ unit vectors. Taking the time derivatives with fixed Lagrangian coordinates

The time derivative of the Jacobian matrix can also be written as

Dotting on the left with $M\u0302$ and on the right with $u\u0302$ gives

Equation (48) for $m\u0302\xb7\u2202u\u0302/\u2202t=\u2212u\u0302\xb7\u2202m\u0302/\u2202t$, then implies

When $\Lambda u\u226b\Lambda m,\u2009M\u0302\xb7(\u2202U\u2192/\u2202t)=M\u0302\xb7G\u2194\xb7U\u0302$, which implies that $U\u0302$ has a continual rotational evolution.

## IV. EVOLUTION OF THE MAGNETIC FIELD

Equation (1) for the ideal evolution of a magnetic field can be solved using Lagrangian coordinates, Eq. (2)

*J* is the Jacobian of Lagrangian coordinates, which is the determinant of the matrix $\u2202x\u2192/\u2202x\u21920$. Equation (62) has a long history, which was reviewed by Stern^{33} in 1966, its importance was recognized in the 2017 review of magnetic reconnection by Zweibel and Yamada,^{13} and a derivation was given in Ref. 7.

A related proof of Eq. (62) is given here for completeness. The equation $(\u2202B\u2192/\u2202t)x\u2192=\u2207\u2192\xd7(u\u2192\xd7B\u2192)$ implies $(\u2202B\u2192/\u2202t)x\u2192=\u2212B\u2192\u2207\u2192\xb7u\u2192+B\u2192\xb7\u2207\u2192u\u2192\u2212u\u2192\xb7\u2207\u2192B\u2192$. Equation (31), which is $(\u2202J/\u2202t)x\u21920=J\u2207\u2192\xb7u\u2192$, and the equation $(\u2202B\u2192/\u2202t)x\u21920=(\u2202B\u2192/\u2202t)x\u2192+u\u2192\xb7\u2207\u2192B\u2192$, imply $(\u2202(JB\u2192)/\u2202t)x\u21920=JB\u2192\xb7\u2207\u2192u\u2192$. Equation (62) is valid if

which can be shown to hold using Eq. (32).

When the magnetic field evolves ideally, Eq. (1), the magnetic field $B\u2192(x\u2192,t)$ at any time *t* is given by Eq. (62). This expression for $B\u2192(x\u2192,t)$ depends on the initial magnetic field $B\u21920(x\u21920)$ and the Jacobian matrix $J\u2194=\u2202x\u2192/\u2202x\u21920$ of Lagrangian coordinates, Eq. (2).

The expression for $B\u2192$ of Eq. (65) can be dotted with itself to obtain the square of the magnetic field strength

The three terms in Eq. (66) for *B*^{2} have fundamentally different time dependencies as Λ_{u} goes to infinity and Λ_{s} goes to zero exponentially. The term in *B*^{2} proportional to $(u\u0302\xb7B\u21920)2$ goes to infinity, the term proportional to $(s\u0302\xb7B\u21920)2$ goes to zero exponentially, while the term proportional to $(m\u0302\xb7B\u21920)2$ changes only moderately. The term proportional to $(u\u0302\xb7B\u21920)2$ is important in dynamo theory,^{30,32} but as a system evolves toward a rapidly reconnecting state, an exponentially large increase in the magnetic field strength is not expected. During the period in which $u\u0302(x\u21920,t)$ relaxes to its steady-state value, Sec. III E, $u\u0302$ must rotate to a direction orthogonal to $B\u21920$, which implies $u\u0302\xb7B\u21920\u21920$.

## V. EVOLUTION OF THE CURRENT DENSITY AND THE LORENTZ FORCE

### A. The current density

The expanded form for the magnetic field, Eq. (65), allows the current density, $j\u2192=\u2207\u2192\xd7B\u2192/\mu 0$, to be determined using the mathematics of general coordinate systems. Unfortunately, these mathematical methods are not known by most plasma physicists, but a two-page derivation is given in the Appendix to Ref. 34.

While applying the method of general coordinates to Lagrangian coordinates, $x\u2192(x\u21920,t)$, the standard convention of using superscripts will be used to number the coordinates. As before, the three Lagrangian coordinates will be denoted using Latin superscripts $x0i$ and the three coordinates in ordinary Cartesian space will be denoted using Greek superscripts *x ^{α}*. The Jacobian matrix $J\u2194\u2261\u2202x\u2192/\u2202x\u21920$ has components $Ji\alpha $, the metric tensor of Lagrangian coordinates, $g\u2194\u2261J\u2194\u2020\xb7J\u2194$ has components

*g*, and the metric tensor of ordinary cartesian coordinates is the unit tensor $1\u2194$, which has the Kronecker delta function,

_{ij}*δ*, as components. With these conventions, dot products are always sums over one superscript and one subscript of the same type, Latin or Greek.

_{αβ}The curl of the magnetic field, when calculated using Lagrangian coordinates, is

where *ϵ ^{ijk}* is the fully antisymmetric tensor:

*ϵ*

^{123}= 1 as are

*ϵ*

^{312}= 1 and

*ϵ*

^{231}= 1, which are even permutations of the indices. The other three permutations are negative, such as

*ϵ*

^{213}= −1. When two indices are identical,

*ϵ*is zero.

^{ijk}$Bj$ is the coefficient of the covariant representation of the magnetic field, $B\u2192=\u2211jBj\u2207\u2192x0j$. When $B\u2192$ is known in covariant form, the curl can be simply calculated. The coefficients $Bj$ can be obtained using the orthogonality relations

which follows from the chain rule. The expanded form of $B\u2192$, Eq. (65), implies

Let $e\u0302$ be any of the three eigenvectors, then the contribution of the $e\u0302\xb7B\u21920$ component of the magnetic field to the current density can be calculated using the representation of the Jacobian matrix expanded in the left and right eigenvectors of the singular value decomposition, Eq. (26). The current density is

### B. The Lorentz force

The force exerted on the plasma, the Lorentz force, can be calculated by crossing the expression for the current density, Eq. (73), with Eq. (65) for the magnetic field. When the magnetic field strength does not increase exponentially with time, only the $m\u0302\xb7B\u21920$ term remains nonzero, and the Lorentz force is

where $U\u0302\xd7M\u0302=S\u0302$ and $S\u0302\xd7M\u0302=\u2212U\u0302$. The term

since $u\u0302\xd7m\u0302=s\u0302$. The definition of $u\u0302\xb7(\u2207\u21920\xd7m\u0302)$ is

The term

The largest term in the Lorentz force is

which predicts an exponentially large force as *t* → *∞* unless $u\u0302\xb7\u2207\u21920\xd7(hmm\u0302)$ relaxes exponentially rapidly to zero.

### C. The parallel current density

The parallel current is given by

As shown in Sec. III E, the quantity $m\u0302\xb7(\u2207\u21920\xd7m\u0302)$ quickly becomes time independent unless $|\u2207\u2192u\u2192|$ scales as Λ_{u}/Λ_{m} or Λ_{m}/Λ_{s}. The middle singular value Λ_{m} changes slowly with respect to time in comparison to the exponential dependences of Λ_{u} and 1/Λ_{s}. As shown in Sec. VI B, the separation between neighboring magnetic field lines is proportional to the variation in $\Lambda u(x\u21920,t)$ along a given magnetic field line, $dx\u21920/d\u2113=m\u0302(x\u21920)$ at a fixed time *t*.

### D. The gradient of $j||/B$

Equation (81) gives an expression for $j||/B$ as a function of $x\u21920$ at fixed time. Before calculating the gradient of $j||/B$, the gradient of an arbitrary function $f(x\u21920)$ will be obtained

where $\u2207\u21920f\u2261\u2202f/\u2202x\u21920$.

The implication is that the gradient of $j||/B$ is exponentially large in the $S\u0302$ direction and exponentially larger than the gradients in the other two directions. Regions of enhanced $j||/B$ are very narrow in the $S\u0302$ direction. The gradient in the $U\u0302$ direction is exponentially small, so regions of enhanced $j||/B$ are very extended in that direction. Both the $S\u0302$ and the $U\u0302$ directions are orthogonal to the magnetic field, which is in the $M\u0302$ direction.

## VI. MAGNETIC FIELD LINES

Appendix A is a direct calculation of the change in the magnetic field produced by nonideal effects when this change is small. What is found is that this change is to lowest order a change in the direction of the magnetic field, $b\u0302\u2248b\u0302I+b\u2192ni$, where $b\u0302I$ is the direction of the magnetic field would have had if the evolution had been ideal. Equation (A20) for $b\u2192ni$ shows the change in direction is proportional to nonideal part of the electric field times the exponentially large coefficient Λ_{u}.

This section shows that an ideal evolution generally leads to an exponentially increasing separation of neighboring field lines with the distance along the lines. This exponentiation in the separation leads to exponential sensitivity to nonideal effects.

The separation of neighboring magnetic field lines can be determined using the Clebsch potentials *α* and *β*, Sec. II B, which have constant values along a magnetic field line, $B\u2192\xb7\u2207\u2192\alpha =B\u2202\alpha /\u2202\u2113=0$ and $B\u2192\xb7\u2207\u2192\beta =B\u2202\beta /\u2202\u2113=0$. The Clebsch potentials can be determined from the starting points of magnetic field lines calculated at a fixed time. The starting points should be on a surface that is nowhere tangential to $B\u2192$.

As discussed in Sec. II B, the equation for the ideal evolution of a magnetic field is satisfied when the Clebsch potentials are functions of the Lagrangian coordinates $x\u21920,\u2009\alpha (x\u21920)$ and $\beta (x\u21920)$.

### A. The gradients of the Clebsch potentials

The gradients of the Clebsch potentials can be calculated using Eq. (84) for the gradient in ordinary space of a function known in Lagrangian coordinates

The initial magnetic field is $B\u21920=\u2207\u21920\alpha \xd7\u2207\u21920\beta $. Since

which implies that $B\u2192=\u2207\u2192\alpha \xd7\u2207\u2192\beta $ reproduces Eq. (66) for the expanded form for $B\u2192$.

### B. Distance between magnetic field lines

Integrations along the field lines of the initial magnetic field, $B\u21920(x\u21920)=\u2207\u21920\alpha \xd7\u2207\u21920\beta $, can be used to determine the functions $\alpha (x\u21920)$ and $\beta (x\u21920)$.

Equation (84) implies that the spatial derivatives in ordinary space of *α* are

When the magnetic field is in the direction $M\u0302$, the derivative $M\u0302\xb7\u2207\u2192\alpha =0$. Analogous equations hold for the derivatives of *β*.

The implication of these expressions is that *α* changes weakly in the $U\u0302$ direction, by an amount proportional to 1/Λ_{u}. Magnetic field lines that are infinitesimally separated in *α* have a separation in the $U\u2192$ direction proportional to $\Lambda u(x\u21920,t)$.

The dependence of the separation along a magnetic field line depends on the variation in Λ_{u} along the line at fixed time *t*. This variation is given by Eq. (56) and is comparable to Λ_{u} itself.

An analogous argument implies that the separation of magnetic field lines in the $S\u0302$ direction is small, proportional to $\Lambda s(x\u21920,t)$.

## VII. LAGRANGIAN COORDINATES IN OTHER FIELDS

### A. Mixing in fluids

A closely related situation to fast magnetic reconnection is mixing in stirred fluids—the advection-diffusion problem—in which it is found that the time required for complete mixing depends only logarithmically on the diffusion coefficient. Aref *et al.*^{35} have written an informative review—particularly the first section. The advection-diffusion equation in Cartesian coordinates, $(\u2202n/\u2202t)x\u2192+u\u2192\xb7\u2207\u2192n=D\u22072n$, when written in Lagrangian coordinates becomes

*J* is the coordinate Jacobian, which is the determinant of the Jacobian matrix $\u2202x\u2192/\u2202x\u21920$, and $\u2202x\u21920/\u2202x\u2192$ is the matrix inverse of the Jacobian matrix. As has been shown, the Jacobian matrix and its inverse generally have elements that become exponentially large as time advances, $\u221d\u2009exp\u2009(\lambda t)$, where *λ* is called a Lyapunov exponent. In a stirred fluid, the effective diffusion coefficient increases exponentially in magnitude until it becomes sufficiently large to flatten spatial variations in *n*. As anyone who has observantly stirred coffee or paint is aware, a certain time is required before the mixing rather suddenly occurs.

Methods of enhancing the effectiveness of stirring in fluids are of practical importance and have received both mathematical and experimental attention. For example, Boyland *et al.*^{36} have used topological concepts to study various protocols for fluid stirring.

The parallel current $j||$ in a reduced-MHD model^{7} of an evolving magnetic field obeys an equation mathematically identical to Eq. (93) with *D* replaced by *η*/*μ*_{0}. The time delay, which is intrinsic to solutions of the advection-diffusion equation, explains the trigger for reconnection. Both the trigger and the speed of fast magnetic reconnection are easily explained in a three-dimensional evolving magnetic field.

### B. Lagrangian coherent structures

The enhancement of mixing is undoubtedly the property described by Lagrangian coordinates with which we are mostly familiar from daily life. Nevertheless, the Lagrangian description is used not only to explain exponentially enhanced mixing but also to describe spatial regions in moving fluids that form barriers to enhanced mixing.

Research on describing barriers to enhanced mixing in fluids provides important insights into incomplete magnetic reconnection when the reconnection proceeds at an Alfvénic rate. A review of the situation in fluids has been written by Haller,^{37} who notes “Lagrangian fluid motion is inherently unstable owing to its sensitivity with respect to initial conditions.” Although the basic result of this instability is rapid mixing, coherent structures are observed that “describe the most repelling, attracting, and shearing material surfaces that form the skeletons of Lagrangian particle dynamics. Uncovering such surfaces from experimental and numerical flow data promises a simplified understanding of the overall flow geometry, an exact quantification of material transport, and a powerful opportunity to forecast, or even influence, large-scale flow features and mixing events.” These are called Lagrangian coherent structures.

Lagrangian coherent structures are given by the ridges and trenches of $\Lambda u(x\u21920,t)$ over some time period. There is large literature on the behavior of $\Lambda u(x\u21920,t)$. The part of the literature that is associated with magnetic field problems was discussed in Sec. III E. Work in the fluid mechanics community is discussed by Haller.^{37}

### C. The standard map

A direct simulation of fast magnetic reconnection or of solutions to the advection-diffusion equation will always be beyond the capability of any computer as the nonideal effects become extremely small, Sec. 3.8 of Ref. 7. Nevertheless, a laptop computer can be used to study the detailed mathematical properties of Lagrangian coordinates and magnetic field lines by iterating the standard map.^{38} The standard map in the *n*th iteration is

where *δφ* is an arbitrary constant. A physical interpretation of the variables of the standard map is that *θ* is a poloidal angle and *φ* is a toroidal angle with *ψ* proportional to the magnetic flux enclosed by a *ψ* surface. The rotational transform $\iota =\psi /\delta \phi $. The parameter *k* is the strength of a perturbation. Since

this map represents the divergence-free character of the magnetic field and can be iterated to plot possible trajectories of magnetic field lines—or of a divergence-free flow. It has the oddity that it is not only periodic *θ* but also in *ψ*. That is *ψ* + 2*π* obeys the same equation as *ψ*. For simplicity of discussion, *k *≥* *0 is assumed.

Greene^{39} found that for $k<0.971635\u2026$ the field lines cover only a limited range of *ψ*. That is, *ψ* remains within the 2*π* periodicity of its initial value. For larger *k*, some lines, though not all, would cover an unbound range of *ψ*. Even for *k *≪* *1, small regions exist in (*ψ*, *θ*) space in which neighboring magnetic field lines separate exponentially as the map is iterated. But, a barrier with a complicated shape exists, crudely at *ψ* = *π*, which field lines cannot cross for $k<0.971635\u2026$. For *k* slightly greater than this critical value, the remnants of the barrier exist in the form of a cantorus^{40} or a Lagrangian coherent structure, which greatly slows the mixing of trajectories from its two sides. Even when *k *=* *2, substantial but isolated regions exist in (*ψ*, *θ*) space in which field lines do not exponentially separate and enhanced mixing would not occur. Indeed, each of these bound regions, which are called islands, is surrounded by its own cantorus. All of these complicated structures have analogues in the behavior of magnetic field lines or of divergence-free flows.

## VIII. DISCUSSION

Much of the existing literature on magnetic reconnection is two dimensional, which in three-dimensional space implies a continuous perfect symmetry in the third direction. Two-dimensional models are inadequate when nonideal effects are weak:

Where is the locus of reconnection?

In two-dimensional theories, reconnection occurs where the two-dimensional part of the magnetic field vanishes independent of the magnitude of the field in the direction of symmetry. In a true three-dimensional system, such as a tokamak subjected to weak three dimensional perturbations, the reconnection occurs at rational surfaces on which magnetic field lines close on themselves. On irrational magnetic surfaces or in stochastic regions, the response to perturbations is fundamentally different because the perturbation is spread over the full volume covered by a magnetic field line—or more precisely the volume covered by an Alfvén wave moving along a field line during the time taken to produce the perturbation. The third coordinate is ignorable only when the symmetry of the original magnetic field and the perturbation holds over the spatial scale sampled by an Alfvén wave during the time required to turn on the perturbation.

Are exponentially large terms ignorable?

For all but exceptional magnetic evolutions, nonideal effects, which are represented by $Eni$, are multiplied by an exponentially large factor, $\Lambda u=e\lambda ut$, to determine the deviation $b\u2192ni$ in the magnetic field direction from the direction it would have had in an ideal evolution. Evolutions that preserve a continuous perfect symmetry are examples of the exceptional evolutions for which the Lyapunov exponent

*λ*= 0._{u}

Three-dimensionality is irrelevant when nonideal effects are sufficiently strong that reconnection occurs before Alfvén waves can propagate a significant distance or before Λ_{u} can become large. But these conditions are very restrictive on the applicability of two-dimensional reconnection theories.

The obvious argument for the use of two-dimensional theories is mathematical simplicity. But in conversations, those who publish work based on two-dimensional analysis claim agreement with observations shows correctness. Nevertheless, no doubt is expressed on the validity of Maxwell's equations or of mathematics. The implication is that it should be mathematically possible to show a more general applicability of two-dimensional models to three dimensional space.

The Lagrangian solutions to the magnetic evolution equations that were developed in this paper give important constraints on the validity of two-coordinate models—the weaker the nonideal effects the more restrictive are the constraints. The smallness of $Eni$, the nonideal part of the electric field, when fast magnetic reconnection occurs makes the conservation of magnetic helicity obvious but the transfer of energy from the magnetic field to the plasma subtle. These issues are discussed in Ref. 41.

Numerical simulation is a way forward for understanding magnetic reconnection in three-dimensional space. But no computer will ever be able to carry out a direct numerical simulation of fast magnetic reconnection in the limit as nonideal effects approach zero. See Sec. 3.8 of Ref. 7. Many reconnection problems of importance will remain beyond the power of direct computations. What can be done is to couple numerical results that are sufficiently restricted to be consistent with existing computational power to extrapolate results to relevant regimes using constraints such as those obtained through methods based on Lagrangian coordinates.

What is clear is that computational and theoretical work on magnetic reconnection must move beyond the two-dimensional models that have dominated the field for more than sixty years and still form the basis of most papers being published.

## ACKNOWLEDGMENTS

This material is based upon the work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences under Award Nos. DE-FG02-95ER54333, DE-FG02-03ER54696, DE-SC0018424, and DE-SC0019479.

### APPENDIX A: NEAR IDEAL CORRECTION TO $B\u2192$

When the electric field has the nonideal form $E\u2192+u\u2192\xd7B\u2192=\u2212\u2207\u2192\Phi +Eni(\alpha ,\beta ,t)\u2207\u2192\u2113$, the equations derived in Sec. II B imply Clebsch potentials have the form

where $Ani$ is the nonideal part of the vector potential, which will be assumed to produce only a small change to the magnetic field from its ideal form $B\u2192I(x\u2192,t)=\u2207\u2192\alpha I\xd7\u2207\u2192\beta I$.

##### 1. Ideal Clebsch potentials as Lagrangian coordinates

The perturbative calculation of the deviation of the magnetic field from that of an ideal evolution uses the Clebsch potentials of an ideal evolution *α _{I}* and

*β*as Lagrangian coordinates. That this is possible is implied by Sec. II B. A third Lagrangian coordinate is required. What will be shown is that

_{I}*ℓ*, the distance along field lines, can be used. To show this, the fundamental definition of Clebsch potentials for a given ideally evolving magnetic field $B\u2192I(x\u2192,t)$ must be considered. This definition uses a two-dimensional Clebsch surface in three-dimensional space.

Let $x\u2192c(\rho ,\beta I,t)$ be a Clebsch surface in three space, which means that it is nowhere tangential to the magnetic field in the region of interest. The two coordinates that give positions on the Clebsch surface are *ρ* and *β _{I}*. Ordinary Cartesian coordinates $x\u2192$ are functions of where a particular magnetic field line penetrates the Clebsch surface and the distance

*ℓ*that the field line must be followed to reach that surface. At each instant of time, the Cartesian coordinates can be given as $x\u2192(\rho ,\beta I,\u2113,t)$ with $\u2202x\u2192/\u2202\u2113=b\u0302I$ with $b\u0302I$ the unit vector along the ideal magnetic field. Since

*ρ*and

*β*are constant along the magnetic field, $B\u2192I=f\u2207\u2192\rho \xd7\u2207\u2192\beta I$. The constraint $\u2207\u2192\xb7B\u2192I=0$ requires

*f*satisfy

*∂f*/

*∂ℓ*= 0. The implication is that $B\u2192I=\u2207\u2192\alpha I\xd7\u2207\u2192\beta I$, where

*∂α*(

_{I}*ρ*,

*β*,

_{I}*t*)/

*∂ρ*=

*f*(

*ρ*,

*β*,

_{I}*t*).

Although *α _{I}* and

*β*can be taken to be constant in Lagrangian coordinates, it remains to be shown that

_{I}*ℓ*can be taken to be a Lagrangian coordinate. As with any function of position and time, the time derivative of

*ℓ*in Lagrangian coordinates is $(\u2202\u2113/\u2202t)x\u21920=(\u2202\u2113/\u2202t)x\u2192+u\u2192\xb7\u2207\u2192\u2113$. Since

*ℓ*is to be a Lagrangian coordinate, $(\u2202\u2113/\u2202t)x\u21920=0$, and the Clebsch surface must be given the velocity $v\u2192c$ required for $(\u2202\u2113/\u2202t)x\u2192=\u2212u\u2192\u22a5\xb7\u2207\u2192\u2113$ with

*ℓ*= 0 on the Clebsch surface. In order not to change

*α*and

_{I}*β*while moving the Clebsch surface, the motion must be along the magnetic field, so the required velocity of the Clebsch surface is

_{I}The quantities that define the velocity of the Clebsch surface are evaluated on the instantaneous Clebsch surface.

##### 2. Calculation of the nonideal field

Using (*α _{I}*,

*β*,

_{I}*ℓ*) as Lagrangian coordinates, Eqs. (A1)–(A3) imply that the magnetic field is

using the dual relations, $Jc(\u2207\u2192\u2113\xd7\u2207\u2192\alpha I)=\u2202x\u2192/\u2202\beta I$ with $1/Jc=(\u2207\u2192\alpha I\xd7\u2207\u2192\beta I)\xb7\u2207\u2192\u2113=BI$ since $B\u2192I\xb7(\u2202x\u2192/\u2202\u2113)=BI$. The dual relations are derived and explained in the Appendix to Ref. 34.

The freedom of canonical transformations contained in *g _{a}*(

*α*,

*β*,

*t*) implies that the two ideal Clebsch potentials can be chosen as

Consistency with Eq. (65) implies

the initial magnetic field strength.

The orthogonality relations, which are derived in the Appendix to Ref. 34, can then be shown to imply

*B _{α}* and

*B*are coefficients in the expansion of the ideal magnetic field in the gradients of the Clebsch coordinates,

_{β}where the expression for $(\u2202x\u2192/\u2202\alpha I)\xb7M\u0302$ follows from Eq. (A12) since $(\u2202x\u2192/\u2202\alpha I)\xb7\u2207\u2192\u2113=0$.

As Λ_{u} → *∞*, the expression for the magnetic field $B\u2192$ has the form

Since the nonideal correction to the magnetic field is an orthogonal direction to the ideal field, the correction to the field strength is negligible in a first order calculation.

In the asymptotic limit as Λ_{u} → *∞*, the deviation in the magnetic field line direction due to nonideal effects is

### APPENDIX B: CURRENT DENSITY IN AN IDEAL-EVOLUTION MODEL

Consider a model in (*x*, *y*, *z*) Cartesian coordinates^{7} in which initially straight magnetic field lines go from a perfectly conducting stationary wall at *z *=* *0 to another perfectly conducting but moving wall at *z *=* L*. The plasma between the walls is assumed to be ideal, $E\u2192+u\u2192\xd7B\u2192=0$. The flow velocity $u\u2192w$ of the *z *=* L* wall is assumed to be on a scale $a\u226aL$ and very slow compared to $(a/L)VA$, so Alfvén waves can keep the magnetic field in a force-free, $j\u2192\xd7B\u2192=0$, state. The evolution of the magnetic field in the plasma, 0 < *z *<* L*, is determined by the vorticity of the slow, divergence-free, flow in the wall $\Omega w(x,y,t)=z\u0302\xb7\u2207\u2192\xd7u\u2192w$. In this model, the distance along the field lines *ℓ* equals *z*

where $\u2207\u22a52\u2261\u22022/\u2202x2+\u22022/\u2202y2$ and *H*(*x*, *y*, *z*, *t*) gives the magnetic field lines *dx*/*dz* = *∂H*/*∂y* and *dy*/*dz* = −*∂H*/*∂x*. A solution to the equations is

where the distribution of parallel current, $K\u2261\mu 0j||/B$, is independent of *ℓ*. When *ϕ _{w}* is independent of time $K$ increases linearly with time.

The magnetic field line labels *α*(*x*, *y*, *z*, *t*) and *β*(*x*, *y*, *z*, *t*) are time independent functions of the (*x*, *y*) coordinates in the *z *=* *0 stationary wall but are generally complicated functions of (*x*, *y*, *z*, *t*) for *z* in the range 0 ≤ *z *≤* L*. On the perfectly conducting wall at *z *=* L*, the magnetic field line labels are equivalent to two-dimensional Lagrangian coordinates for the flow of the wall. $\u2202\alpha (x,y,L,t)/\u2202t=\u2212u\u2192w\xb7\u2207\u2192\alpha $ and $\u2202\beta (x,y,L,t)/\u2202t=\u2212u\u2192w\xb7\u2207\u2192\beta $.

The solution given by Eqs. (B5) and (B6) may become unstable; when this occurs the relevant solution is more complicated.

When the flow velocity in the *z *=* L* wall is written in terms of its stream function, $u\u2192w=\u2207\u2192\xd7(\varphi w(x,y,t)z\u0302)$, the streamlines of the wall flow are *dx*/*dt* = *∂ϕ _{w}*/

*∂y*and

*dy*/

*dt*= −

*∂ϕ*/

_{w}*∂x*, which are of identical form to the equations of classical mechanics of one and a half degrees of freedom,

*H*(

*p*,

*q*,

*t*). When

*∂ϕ*/

_{w}*∂t*=

*0, the stream lines remain on constant*

*ϕ*surfaces, and the regions in the

_{w}*z*=

*L*surface in which the streamlines increase the separation by

*e*over time occupy an

^{σ}*e*

^{−2σ}fraction of the surface area. When

*ϕ*is time dependent, then generally a large fraction of the area is occupied by streamlines that increase their separation exponentially in time, as $e\lambda st$, where

_{w}*λ*> 0 is called a Lyapunov exponent. Even simple stream functions can produce complicated spatial distributions of $j||/B$. An example is a Fourier representation with a periodicity distance

_{s}*a*that has at least two terms

where the *ϕ _{mn}* and the

*ω*are distinct constants.

_{mn}