## Abstract

This paper demonstrates an equivalence between rotating magnetized shear flows and a stressed elastic beam. This results from finding the same form of dynamical equations after an asymptotic reduction of the axis-symmetric magnetorotational instability (MRI) under the assumption of almost-critical driving. The analysis considers the MRI dynamics in a non-dissipative near-equilibrium regime. Both the magnetic and elastic systems reduce to a simple one-dimensional wave equation with a non-local nonlinear feedback. Under transformation, the equation comprises a large number of mean-field interacting Duffing oscillators. This system was the first proven example of a strange attractor in a partial differential equation. Finding the same reduced equation in two natural applications suggests the model might result from other applications and could fall into a universal class based on symmetry.

## 1. Introduction

In a non-dissipative, purely hydrodynamical context, the famous Rayleigh criterion [1,2] of outwardly increasing angular momentum governs the stability of rotating shear flows. That is,
*r* represents the outward directed radial coordinate, and *Ω*(*r*) represents the local angular rotation rate. During the 1950s and 1960s, Chandrasekhar [3,4] in The West, and Velikov [5] in the former Soviet Union considered the effect of an axial magnetic field on the stability of a cylindrically swirling conducting fluid, thereby modifying the Rayleigh criterion. In these works, both authors found independently that (in the ideal non-dissipative regime), the presence of a magnetic field can catalyse an instability for a merely inwardly increasing angular frequency profile; rather than angular momentum. That is, a magnetohydrodynamic (MHD) instability only requires

In the early 1990s, Balbus & Hawley [7] auspiciously applied the destabilizing of nature of a weak magnetic field to solve a previous astrophysical paradox. Specifically, in a fluid system orbiting a central mass, *M*, a balance between gravitational and centrifugal acceleration chiefly determines the rotation profile, *rΩ*(*r*)^{2}=*GM*/*r*^{2}; *G* represents Newton's gravitational constant. This *Keplerian* profile implies

The advent of the magnetorotational instability (MRI) in astrophysics produced many notable results over the last three decades [8]. Little doubt exists regarding the operation and efficiency of the MRI in hot astrophysical disc systems. But questions do remain concerning the instability's cessation and saturation [11,12]. This paper explores the dynamics of the MRI in a controlled setting. This simple approach helps explain some mysterious aspects regarding how the instability transports momentum and magnetic flux, which proves useful in understanding possible saturation mechanisms.

In contrast to a magnetized fluid, an elastic material is a continuum substance that returns to its original form with the removal of applied loads. In spite of clear differences, elastodynamics and MHD contain many deep physical and mathematical similarities. And just with forced fluids, elastic materials possess an extremely rich range of behaviours, specifically including dynamical instabilities.

When a solid column *buckles*, it loses the ability to support a load while retaining its elastic integrity. Because of its obvious importance to the stability of structures, the scientific investigation of buckling dates well back into antiquity, and quantitative investigations started at the very beginning of the modern scientific era. As early as the 1480s, Leonardo da Vinci produced empirical criteria addressing the stability and lateral deflection of columns under compression [13]. In the eighteenth century, Euler and Bernoulli began to consider time-dependent elastic deformations and buckling [14]. In particular, Euler derived a practical formula for calculating the critical load on a slender beam. Engineers still use this formula (and its extensions) in modern construction design [15].

Elastic substances exhibit many more types of instabilities than simple buckling [16]. After much success with linear dynamics of elastic solids, work in the 1970s began to focus heavily on nonlinear finite deformations. Aeronautical engineering and manufacturing particularly requires understanding the time-dependent effects of solid materials under a wide range of situations. The introduction of methods from mathematical dynamical systems theory allowed rigorous analysis for a wide range of scenarios [17]. In a series of work, Holmes and Marsden considered a class of nonlinear, non-local models for buckling and fluttering [18,19]. They proved the first example in a partial differential equation of chaos and a strange attractor; i.e. an infinity of arbitrarily long periodic orbits.

This paper finds an almost identical equation arising naturally from the weakly nonlinear theory of the MRI in a simple geometry. This link gives qualitative insight into the astrophysical setting, and a more quantitative understanding of magneto-Taylor–Couette flow in laboratory experiments [20]. In contrast to astrophysical discs, the model may also apply to the inside of some stars with stronger ambient magnetic field and more modest differential rotation. However, even if no astrophysical object exists near-equilibrium, understanding nonlinear behaviour near a phase transition can uncover fundamental interactions. In this case, the current model helps elucidate transport of momentum and magnetic flux in an important MHD instability. To use a zoological analogy: underneath many obvious differences, both the buckling beam and the wild-type MRI contain a morphologically equivalent primitive skeleton.

The remainder of this paper is as follows: §2 analyses both the linear and nonlinear theory of the MRI. The main result of the paper follows from (2.63) and (2.64) and its analogy to (3.24); §3 derives the buckling instability of a slender elastic beam; §4 provides a detailed discussion on several aspects of the MRI and the buckling beam; §5 shows computational results of the nonlinear solutions and §6 provides conclusion.

## 2. Instability of a rotating magnetohydrodynamic shear flow

The first part of this paper considers a two-dimensional model for the MRI of an incompressible magnetized fluid with linear shear flow. The system remains invariant along the shear direction, but contains all three components of flow and magnetic field. The full primitive equations describe the dynamics in a rotating coordinate system,
_{y}=0.

The stream-function and magnetic-scalar-potential formulations enforce the solenoidal character of flow and magnetic field (2.2) and (2.4),
*v*(*x*,*z*), *b*(*x*,*z*) represent the flow and magnetic field in the stream-wise direction, respectively, and *ψ*(*y*,*z*), *a*(*y*,*z*) represent the poloidal stream function and scalar potential in the plane perpendicular to the streaming flow. Both *v*(*x*,*z*) and *a*(*x*,*z*) represent perturbations relative to the linear background profiles
*x*≤*d*. This represents an arbitrary gauge for the magnetic potential, but complies with the definition of the global rotations rate for the shear flow. The linear background magnetic potential implies a uniform background magnetic field in the vertical direction, *S* represents the local background shear rate. Figure 1*a* depicts the basic geometry and background configuration.

Obtaining equations for the stream function and magnetic scalar potential requires taking the two-dimensional curl of the *x* and *z* components of the flow and magnetic field equations. The incompressible MHD equations for each scalar variable are
*μ*_{0}*ρ*_{0}=1. We restore these parameters at the conclusion of the derivations. For the remaining parameters, *f*=2*Ω* represents the background vorticity owing to the frame rotation, *Ω*. The dissipation parameters, *ν* and *η* represents viscosity and magnetic diffusivity, respectively. The following analysis considers the cases of both small and/or completely vanishing diffusion coefficients.

### (a) Linear analysis

Temporarily neglecting the nonlinear and dissipative terms in (2.8)–(2.11) helps determine the relevant spatio-temporal scales. After linearizing, and with some simplifications,
*ψ*(*t*,*x*,*z*), (2.13), (2.14) and (2.15) determine *b*(*t*,*x*,*z*), *a*(*t*,*x*,*z*) and *v*(*t*,*x*,*z*), respectively. (2.16) determines the stream function, subject to the impenetrable boundary conditions
*ω*, versus wavenumber *k*, and other parameters,
*ω*^{2}>0) to purely imaginary frequencies (*ω*^{2}<0) characterizes any transition of stability. For any real *k*,
*ω*^{2}=0) between these two regimes. Instability (*ω*^{2}<0) exists for

The difference between the critical shear and actual shear helps simplify the analysis in the vicinity of the instability,
*kd*|≪1, and near-critical shear |*σ*|≪|*S*_{crit.}|, we can expand the dispersion relation to lowest order and find
*f*|*σ*|*d*^{2}/*B*^{2}≪1. Figure 2*a* shows the behaviour of (2.25).

Reintroducing space and time variables via *ω*→−*i*∂_{t} and *k*→−*i*∂_{z} implies (heuristically) that (2.25) represents the simplified (unstable) wave equation

### (b) Nonlinear asymptotic analysis

Rendering the dynamical equations in dimensionless form often helps to simplify derivations. In this case, working with scaled quantities partially obscures the physical meaning of different dynamical ingredients. Therefore, the following analysis introduces a small non-dimensional parameter (*ε*≪1) only to keep track of the relative magnitudes of the different terms. We rescale all relevant dynamical variables in terms of particular powers of *ε* and conduct an asymptotic analysis accordingly. From (2.25), we deduce the consistent rescaling
*b*, *a* and *v* in terms of *ψ*. The fully nonlinear dissipative system, (2.8)–(2.11), constrains the amplitude of *ψ*, and also gives the magnitude of *ν* and *η* that allow dynamically significant dissipation. That is,

Upon rescaling, (2.8)–(2.11) become

Considering this system order-by-order, the analysis halts when the first time derivative from the right-hand side enters the balance. Therefore, we may automatically neglect terms that are formally smaller than the highest order time-derivative. Other than not writing higher order terms, (2.30)–(2.33) are completely equivalent to (2.8)–(2.11).

We expand all dynamical variables in the following series:
*ε*. The leading-order balance produces

Solving (2.36) and (2.37) in general form requires introducing a potential function *φ*(*t*,*z*)

At this point, *φ*(*t*,*z*) represents all the dynamical degrees of freedom left undetermined by (2.36) and (2.37). This minimalist approach solves (2.36) and (2.37) without adding any more information about the form of the solution until nonlinearity dictates it. This single unknown function represents a scalar-valued ‘order parameter’ in the language of critical phenomena. The typical value of *φ*(*t*,*z*) grows in amplitude and becomes increasingly disordered as the system becomes more unstable.

(2.42) matches with the critical value resulting from linear theory, (2.23). We choose boundary conditions for the magnetic field that match the natural non-dissipative profiles. Other choice of boundary conditions would bring boundary layers into the analysis. We filter these effects for simplicity.

The next-order balance produces
*z*-dependent components of *v*_{1} and *a*_{1}, which exactly mimic the form of *v*_{0} and *a*_{0}. (2.43) and (2.44) do not determine the *z*-mean components of these variables.

Averaging (2.45) and (2.46) in the vertical direction eliminates all pure *z* derivatives and produces the balances
_{z}*φ* drives purely *x*-dependent corrections to the mean background field and shear momentum.

Closing the system in terms of *φ* requires generating a second-order equation for *a*_{2}. After combining second-order versions of (2.36),
*a*_{2} satisfying the boundary conditions [21]

Grouping all evolution equations together,
*q*=1 corresponds to the Rayleigh instability threshold. In the case of fixed *q*, we may think of the background magnetic field taking the role of the critical parameter.

(2.56) contains the mean-field shear parameter

In terms of *q*,
*q*<1, both *γ*_{3}>0 and *γ*_{4}>0 are well-defined averages of *γ*_{1} and *γ*_{2}. The dissipation coefficients take only one of the two possible orderings

Assuming *η*=*ν*=0 allows reducing the system even further. In this case, (2.56) and (2.57) integrate explicitly so that

Equation (2.63) differs significantly from previous attempts at a weakly nonlinear model for the MRI [22,23]. Those cases consider geometry more similar to accretion discs, and as a result required significantly higher (leading-order) dissipation. Past results produce a dissipative Ginzburg–Landau equation, not the high Reynolds number model produced here. However, for large dissipation (2.56) reduces to a first-order equation in ∂_{t} and hence more like Ginzburg–Landau.

Liverts *et al*. [24] derived an ordinary differential Duffing equation for the thin-disc MRI. We see in §5 that (2.63) represents a possible infinity of coupled Duffing equations. This correspondence shows strong evidence that the ideal Duffing-type dynamics underlies the MRI at a fundamental level.

## 3. Buckling instability of an elastic beam

This section derives a dynamical nonlinear model for an elastic beam in the vicinity of a buckling instability. Holmes gave a very short equivocal version of this same derivation [18]. The current context needs enough details to see the similarities and differences to the magnetic case. We use an asymptotic expansion in terms of the beam aspect ratio. In this section,
*z*≤*L* represents the length of the beam along the direction of applied load, and −*d*/2≤*x*≤*d*/2 represents cross-sectional thickness on the beam. Figure 1*b* depicts the basic geometry and background configuration. A fully complete treatment of fully nonlinear elastodynamics would lead to a very complicated asymptotic analysis in powers of *ε*. Even though we apply some intuitive reasoning, a completely systematic analysis gives the same eventual answer.

The nonlinear strain tensor characterizes a general Lagrangian deformation of a continuum solid [25]. Therefore,

where |d** x**|

^{2}=d

*x*

^{2}+d

*z*

^{2}represents the

*original*distance between nearby points, and |d

**+d**

*x***|**

*ξ*^{2}represents the distance between nearby points after a Lagrangian displacement,

**(**

*ξ**t*,

**). The explicit components of Green's strain tensor are**

*x*We consider a homogeneous isotropic Hookean solid with linear stress–strain relationship
** E**)=

*E*

_{x,x}+

*E*

_{z,z}, and

**represents the identity matrix. The parameters λ and**

*I**μ*represent Lamb's constants, but using the alternative definitions proves advantageous,

*ν*denotes Poisson's ratio (note: this is

*not*the same as the viscosity parameter in the MRI analysis), and

*Y*represents Young's modulus. Some texts use

*Y*/(1−

*ν*

^{2}) to denote Young's modulus e.g. [25]. Our current definition coincides with [18] and allows for a simpler derivation.

The thin aspect ratio leads to small displacements and variation in the longitudinal direction, while the perpendicular displacements and variation remain order unity. Kinetic energy also must balance stresses. Therefore, we replace
*u*(*z*) and *w*(*z*) represent unknown mean displacements in the *x* and *z* directions, respectively, and primes, e.g. *u*′(*z*), denote *z* derivatives. (3.8) and (3.9) imply
*E*_{z,z} and *S*_{z,z} components contribute to leading order. To leading order, only the horizontal displacement contributes to the kinetic energy
*x* dependence in the potential energy allows integrating out this dimension, which produces the effective one-dimensional Lagrangian density
*w*(*t*,*z*) implies
*z*, but depending on *t*). Varying the action with respect to *u*(*t*,*z*) implies that
*L*≡*w*(*z*=*L*)−*w*(*z*=0) gives the total change in the length of the rod.

From the definition of Young's modulus applied to the entire solid,
*Γ* denote the total applied compression (*Γ*>0), or tensile (*Γ*<0) load. Therefore,

## 4. Discussion

### (a) Dynamical correspondences

Inspecting each term in (2.63) and (3.24) highlights the analogies between the various physical parameters in the two cases. To simplify the derivations in §2, (2.8)–(2.11) use Alfvén units such that *μ*_{0}*ρ*_{0}=1, where *μ*_{0}=4*π* in cgs units. This section restores the general parameters in order to make comparisons between MHD and elastic parameters.

The shear criticality in the MHD system and the applied compression/tension loading correspond such that,
*B*^{2}/*μ*_{0} lends the correct units to the right-hand side of the relationship.

The magnetic tension in the background magnetic field corresponds to Young's modulus in the elastic setting

Together, (4.1) and (4.2) relate the shear criticality in the fluid to the total strain (relative compression) in the solid

Lastly, the comparing the nonlinear terms in both systems gives interpretation for the scalar potential in MHD in terms of the horizontal beam displacement.
*u*(*t*,*z*)∝*φ*(*t*,*z*)/*d*.

The potential *φ* does not relate naturally to the horizontal Lagrangian displacement in the fluid system. The definition of the potential function is somewhat arbitrary, but (2.38)–(2.41) produce the most straightforward solution to (2.36) and (2.37) in terms of derivatives; *b*,*ψ*∝∂_{t}*φ*; *a*,*v*∝∂_{z}*φ*. The constants of proportionality are more chosen according to style. In terms of (fluid) Lagrangian displacements, *ξ*_{x}∝∂_{z}*φ*; *ξ*_{z}∝*φ*, whereas for the elastic solid *ξ*_{x}∝*u*,∂_{z}*w*; *ξ*_{z}∝*w*,∂_{z}*u*. The correspondence between the parameter in the MRI and the buckling beam would allow the demonstration of some aspects of MHD with a more manageable slender elastic rod.

### (b) Saturation

How the MRI saturates remains an open question in the context of accretion disc modelling [12,24]. Intuitive understanding in the case of the elastic beam, helps clarify the saturation mechanism of the MRI in this simple model. An elastic beam saturates a buckling instability by setting one side of the beam under tension and the other side under compression. In a sense, the instability transports stress (equivalently strain and density) dynamically from the tension side to the compression side. In the MHD case, (2.51) and (2.62) imply that the instability puts the left-half (0<*x*<*d*/2) of the domain under a slight excess of magnetic flux and slight deficit of linear momentum. The right-half (*d*/2<*x*<*d*) receives the opposite feedback such that the total of both quantities remains conserved. In the linear phase of the instability, momentum and magnetic flux compete to both drive and suppress growing perturbations. The nonlinear transport rearranges the background so that both halves of the domain experience more stringent stability criteria. In their early MRI studies, Balbus & Hawley [10] demonstrated positive outward angular momentum flux. From (2.62), 〈*V* 〉<0 also corresponds to positive outward transport. Even in the local, incompressible, Cartesian geometry, the system can correctly discern inward and outward. Along with angular momentum transport, the exchange of magnetic potential reconciles with past accretion disc models [26].

### (c) Symmetry

The buckling beam analogy is not the first between MHD and elastic systems [27,28]. There are deep mathematical reasons for this relationship. Chandrasekhar speculated that the reason for not recovering Rayleigh's criteria (1.1) in the limit of zero magnetic field ‘must lie in the circumstances that … the lines of magnetic force are permanently attached to the fluid’ [6], section 81, p. 389. This is another way of saying that specifying the location of the magnetic field lines is the same as specifying the location of the fluid parcels.

Like an elastic medium, MHD supports shear stresses and shear waves in the form of Alfvén displacements. Symmetry breaking provides the mechanism underlying this simple physical fact. When considering a collection of *N* particles, phase space generically comprises 3*N* momenta and 3*N* coordinates for a total of 6*N* dimensions. For a fluid, the range scales of fluctuations determines the effective *N*. Independence of the system on one or more of the coordinates implies conservation of the corresponding momenta. Ordinary fully compressible hydrodynamics is only 5*N* dimensional, 3*N* for the velocities, 1*N* for the density, and 1*N* for entropy (or equivalently the pressure). The missing third set of coordinates implies the conservation of potential vorticity (effectively a component of momenta). There is no restoring force for one component of the momentum (more accurately one subset of size *N*). For incompressible hydrodynamics a much smaller set of unique coordinate labels implies Kelvin's Circulation Theorem which implies a much larger set of conserved momenta. Magnetic field changes this picture.

The frozen-in condition [29] provides the magnetic fluid with enough Lagrangian coordinate labels to break potential vorticity conservation. In fact, MHD seems to provide more coordinate labels than necessary to specify a location in phase space completely. The divergence condition, and other constraints, remedy the over-counting problem; see [30] for more details. With a full accounting for fluid parcels, there remain no point-wise conserved momenta, and almost all degrees of freedom must experience a non-zero force of some kind. Lack of Lagrangian-particle relabelling symmetry is an imperative feature MHD shares with elastodynamics. In the latter case, the dynamics depends explicitly on all components of the displacement field through the strain tensor.

Along with their primary discovery, Balbus & Hawley [31] produced an analogy between magnetic field lines between Keplerian fluid parcels and a simple elastic spring connecting two orbiting masses. In this example, the connecting spring breaks the conservation of the individual masses angular momentum and Laplace–Runge–Lenz vectors. This example shows that the tension in the magnetic field and/or spring provides both stabilization via tension, and acts to liberate rotational energy. The buckling beam differs in this respect. The product of magnetic field and shear compares to the compression in (4.1). The magnetic field alone (without destabilizing shear) compares to Young's modulus in (4.2).

### (d) Non-locality

The non-local nature of the nonlinear feedback requires some consideration. The nonlinear feedback in (2.63) appears surprising at first. Especially given that the original MHD model does not seem to contain non-local terms. How does the non-locality arise?

Answering this requires pointing out a number of small facts. (i) In spite of appearances the original equations (2.8)–(2.11) are actually non-local. This results from evolving *ω*_{y}=∇^{2}*ψ*, rather than *ψ* alone in (2.11). Computing the inverse Laplacian amounts to convolving against the appropriate Green's function. This is a non-local operation. (ii) Unlike the incompressible MHD equations, the full elastodynamic equations are completely local. (iii) Therefore, the reason for non-locality in (2.63) cannot result from the non-locality (2.11). But similar physics underlies both situations. (iv) Non-locality in incompressible MHD results physically from fast acoustic waves equalizing the pressure field on timescales much faster than dynamical timescales. Rather than experiencing a small delay, distant points respond instantaneously to fluid motions, but in a fashion that attenuates with distance away from the disturbance. (v) The buckling beam achieves non-locality through of the separation of timescales between extremely fast elastic waves, and unstable flexural waves. (vi) In addition to the unstable MRI modes, the full linear dispersion relation (2.19) contains two additional fast modes for each wavenumber. (vii) These mixed Alfvén–Coriolis modes propagate with the phase speed

Incompressibility implies that non-locality is more common than often realized. The surprising aspect about (2.63) is that feedback is entirely non-local. (2.64) shows no spatial variation or attenuation. We could express (2.64) in local form by stating

### (e) Dissipative stability

Also, in contrast to the ideal stability theory in §2a the nonlinear analysis elucidates stability with a small amount of diffusion. In that case there exist a critical shear curve as a function of wavenumber, i.e.
*k*, gives
*b* shows the difference in stability for dissipative versus non-dissipative dynamics. This adds further insight into breaking the apparent degeneracies associated with the small-dissipation limit of the MRI; also see [33,34]. In the limit as the background magnetic field vanishes, we find a necessary condition in terms of the magnetic Reynolds number,

## 5. Nonlinear solutions

Rescaling the time, space and amplitude variables simplifies the discussion of nonlinear solutions to (2.63) and (3.24). Assuming *η*=*ν*=0, and defining a non-dimensional length, time and amplitude leads to
*μ*=±1. It is straight forward in both the MHD and elastic cases to find the relevant re-scalings that yield (5.1).

The completely non-local character of the nonlinearity in (5.1) implies a particularly interesting class of solutions. The spatial Fourier transform
*φ*(*t*,*y*) is real implies the Hermitian symmetry

The original nonlinear system (with *η*=*ν*=0) is conservative, and the reduced system remains Hamiltonian with the total conserved energy

The completely mean-field character of (5.3) implies interesting symmetries when considering the evolution of the complex-valued amplitude in terms of the real amplitude and phase
*for each* individual value of *k*, each ‘angular momentum’ remains constant in time
*L*_{k} results from the symmetry for each *k*, *χ*_{k}.

The system further reduces to the following real equations:
*L*_{k} is an arbitrary function of *k* and is independent of time. The total energy now becomes

In terms of the complex initial data,

Furthermore, the inverse Fourier transform of *L*_{k} gives a point-wise conserved quantity in the spatial domain,
*z*, ∂_{t}*Λ*(*t*,*z*)=0. The conservation of *Λ*(*t*,*z*) is equivalent to the conservation of
*s*. The case *s*=1,3, respectively, correspond to magnetic-helicity and cross-helicity conservation in the original MHD system. Even though the system contains an infinite number of conserved quantities, it is not apparently integrable.

(5.8) shows that the dynamics is equivalent to a collections of particles all moving in an averaged potential, and each individually conserving its angular momentum. The non-local character of the nonlinearity implies that *every* Fourier truncation of the system represent an *exact* solution. In other words, the system does not cascade energy to smaller scales than are present in the initial conditions. This unusual form of memory implies an extremely diverse class of solutions. The model displays one of the simplest types of mean-field dynamics, but with non-trivial results.

Figure 3 shows a typical solution for a modest number of initial modes. The solution lies in a periodic domain 0≤*z*≤2*π*, with *μ*=+1. The initial conditions consist of *K*=20 Gaussian-random unit-amplitude complex-valued Fourier modes for both *φ* and

Figure 3 shows a pattern of large-scale travelling waves with faster small-scale dynamics superimposed. The small-scale fast dynamics show up clearly in the velocity pattern, ∂_{t}*φ*. The qualitative aspects of this pattern seem typical for an array of initial conditions. A more comprehensive numerical study of the space of solutions both with and without dissipation lies beyond the scope of this paper. For an example of the possible range-rich behaviours, Eugeni *et al*. [35] recently studied the forced nonlinear beam dynamics with multiple interacting degrees of freedom.

## 6. Conclusion

Heuristically, one can use the concept of entropy to interpret two different physical systems producing the same reduced governing equations. Near an instability, only a small number of degrees of freedom grow to significance. There are many fewer ways to organize a small number of modes than the large number active in a strongly driven system. Or, to borrow from Tolstoy, ‘Weakly nonlinear systems are all alike; every strongly nonlinear system is nonlinear in its own way’.^{1} This is especially true in the presence of symmetry. This fact is apparent in that so many systems reduced to a small set of well-known canonical equations. For example, the complex Ginzburg–Landau equation (and its further simplifications) derives systematically from an extremely wide range of physical starting points [36].

This paper points out a new class of *non-local* models deriving from at least two different natural systems. More physical examples likely exist. In a particularly fundamental way, MHD is more similar to elastodynamics than pure hydrodynamics. The magnetic field provides a final material coordinate label, and therefore eliminates fluid particle relabelling symmetry, and conservation of potential vorticity. In elastodynamics, material particles contain explicit labels, and the issue of potential vorticity never really arises. This implies that fluid systems without conserved potential vorticity, and a spring-like restoring mechanism could produce similar dynamics as the MRI or buckling beam. Over-stable double-diffusive convection satisfies these conditions, and (speculatively) may provide an additional example. It is quite possible that this mean-field network of Duffing oscillators forms a new class of universal near-equilibrium dynamics.

More generally, there are infinitely more ways to produce a non-local interaction than a local one. The multiple scales assumption allows for information to travel asymptotically fast and produce the instantaneous interaction we see in our model. In both the buckling beam, and MRI, the reason for the particular type of nonlinearity arises from conservation laws. In the elastic case, the integral terms result from fixing the total mass of the beam. In the MHD case, the integral terms arise from conservation of total linear momentum, and magnetic flux. In both cases, these quantities remain conserved with the inclusion of dissipative effects. Rotating and/or magnetized free-slip Rayleigh–Bénard convection, displays non-local nonlinear terms arising from conservation of total momentum and/or magnetic flux; just as in the MRI problem [37,38]. These correspondences perhaps imply that global conservation principles and non-locality relate more deeply. Intriguingly, the non-local terms in the MRI model allow for a very large class of symmetries and this hints at the deeper link to conservation principles.

Elastodynamics and MHD probably part ways with the introduction of more physical ingredients. The systems are likely very different from each other in three dimensions. Dissipative boundary layer dynamics (filtered in this paper with judicious choice of boundary conditions), interactions with other fluid instabilities, and more complex geometric effects all likely pull the correspondence further apart. Nevertheless, the link between elastic buckling and the MRI produces interesting insight regarding the growth and saturation of the MRI, albeit in a regime far removed from traditional astrophysical applications. The type of derivation presented in this paper would probably work in a thin disc-like geometry, which is more suitable to accretion discs. This interesting case will produce several additional complications resulting from boundary conditions alone. Attempting to understand more complex systems requires adding dynamical elements back into the minimalist model and considering the consequences. Hopefully adding richer dynamical ingredients will still allow simple mathematical progress. But even when interactions become too numerous to consider analytically, understanding the relevant spatio-temporal scales and possible couplings will greatly help streamline numerical simulation of more complex models.

## Data accessibility

This work does not have any experimental data.

## Funding statement

The Australian Research Council supports the author through a Discovery Early Career Researcher Award, no. DE140101960. Part of this work was conducted with funding from the University of California Berkeley Theoretical Astrophysics Center.

## Conflict of interests

I have no competing interests.

## Acknowledgements

The author would like to thank Edgar Knobloch for originally pointing out Holmes and Marsden's work on elastic buckling. The computations used the Dedalus code: dedalus-project.org. The author thanks the remainder of the Dedalus collaboration: Keaton Burns, Daniel Lecoanet, Jeff Oishi and Ben Brown for their significant contributions to the code project. The author thanks two anonymous referees for suggesting improvements to the manuscript.

## Footnotes

↵1 Tolstoy's original quote: ‘Happy families are all alike; every unhappy family is unhappy in its own way.’, concisely describes how some macro states can contain vastly different numbers of possible micro states. Anna Karenina was published within a couple years of Boltzmann's original mathematical formulation of entropy.

- Received September 16, 2014.
- Accepted March 12, 2015.

- © 2015 The Author(s) Published by the Royal Society. All rights reserved.