A two-dimensional magnetohydrodynamic system: geometric decomposition and canonical reduction

A geometric decomposition previously adopted in areas of nonlinear continuum mechanics and soliton theory is applied to a two-dimensional time-independent magnetohydrodynamic system to obtain reduction to a canonical third order nonlinear equation subject to a constraint. Novel classes of exact solutions to the magnetohydrodynamic system are thereby constructed. Connection is made with kinematic conditions which attend certain steady two-dimensional motions of fibre-reinforced viscous fluids.


Introduction
In [1], an intrinsic geometric decomposition procedure was developed in a kinematic analysis of certain hydrodynamic motions. It has subsequently been applied in diverse areas such as, inter alia, nonlinear elasticity, gasdynamics, magnetohydrostatics, magnetohydrodynamics and modern soliton theory. Thus, in magnetohydrodynamics it has been applied to analyse motions with the admitted geometric characteristic that the Maxwellian surfaces contain the streamlines as geodesics and the magnetic lines as parallels on a normal congruence [2]. Viable helical geometries were thereby isolated. In soliton theory, this kind of anholonomic representation has been applied, in particular, in the derivation of both an auto-Bäcklund transformation and Lax pair in a purely geometric manner for the canonical nonlinear Schrödinger (NLS) equation [3]. NLS soliton surfaces may be thereby constructed systematically in a geometric way [4]. In [5], this intrinsic geometric approach uncovered a remarkable connection between a previously well-studied class of hydrodynamic motions [6][7][8][9] and the integrable solitonic Heisenberg spin equation subject to a geometric constraint. The nature of this Heisenberg spin link has subsequently been elaborated upon in [10][11][12][13] and most recently extended to spatial relativistic gasdynamics [14]. The corresponding solitonic connection in magnetohydrostatics to the S-integrable homogeneous Heisenberg spin equation has been employed in [15] to derive exact magnetic flux configurations which incorporate, in particular, nested toroidal geometries as originally described in [16] with relevance to tokamak design. An extension to magnetohydrodynamics leading to the solitonic Pohlmeyer-Lund-Regge model subject to a volumepreserving constraint was presented in [17]. This geometric procedure was again employed in soliton theory in [18] in an analysis of the spatial binormal motion of inextensible curves of either constant curvature or torsion. Novel integrable extensions with geometric origin of both the Dym equation [19] and classical sine-Gordon equation of Bour [20] together with associated Bäcklund transformations [21] and soliton surfaces were thereby constructed.
The important connection between certain motions of inextensible curves and the generation of associated soliton equations is well-established (see, e.g., [4,22] and literature cited therein). The study of the motion of fibre-reinforced fluids on the other hand, as described in the authoritative monograph of Spencer [23] likewise requires the analysis of viable motions of inextensible curves, in that physical context, embedded in a viscous liquid. This commonality with soliton theory involving the constrained motion of inextensible curves indicated potential integrable structure underlying certain admitted deformations in the motion of fibre-reinforced materials in formation processes. Such integrable connections have indeed been established in [24][25][26].
The basic mathematical theory of the deformation of fibre-reinforced materials is contained in [27]. However, it is remarked that certain aspects of the deformation of elastic materials reinforced by inextensible chords had been treated in the classic work by Green and Adkins [28]. In [24,25], it was established by geometric means that kinematic constraints as set down in [23] on the steady two-dimensional motion of an ideal fibre-reinforced fluid can be consolidated in a single 3rd order nonlinear equation. Importantly, the latter admits a particular solitonic reduction compatible with the classical sine-Gordon equation. The kinematic conditions in that case admit both a Bäcklund transformation and a novel duality property. In [24], kinematically admissible spatial motions were obtained in which the fibres constitute geodesic windings on nested toroidal surfaces. The analysis of the kinematics of 2+1-dimensional motions of a fibrereinforced fluid was subsequently undertaken in [26]. Under a constant divergence constraint, remarkably, the evolution was shown to be linked to the inverse scattering problem for the solitonic modified Korteweg-de Vries hierarchy which, in turn, may be generated in a purely geometric manner [4]. Motions were isolated in [26] wherein the fibres which are convected by the viscous fluid constitute generalized tractrices.
Two-dimensional steady magnetohydrodynamics has previously been the subject of analysis when the magnetic and velocity fields are everywhere orthogonal [29] or, more generally, inclined at a constant angle [30]. These cases were subsequently subsumed in non-aligned planar magnetohydrodynamics in [31] via a novel representation involving a quartet of Bernoulli integrals admitted by the governing nonlinear system. The present work serves to complement that of [31] and applies an intrinsic geometric formulation. It emerges here that, remarkably, a key 3rd order nonlinear equation of a geometric nature as originally derived in the fibre-reinforced viscous liquid context of [24,25] naturally arises in the present magnetohydrodynamic analysis. In most recent work [32] a time-dependent variant of this canonical equation has been derived in the description of the 2+1-dimensional motion of fibre-reinforced fluids. Importantly, this development indicates the potential to extend the present geometric analysis to 2+1-dimensional magnetodydrodynamics and magnetogasdynamics. It is remarked that, in the latter connection, novel integrable pulsrodon [33] and Ermakov-Ray-Reid structure [34,35] has been previously isolated in [36].

The magnetohydrodynamic system
Here, a geometric decomposition is undertaken of the two-dimensional reduction of the time-independent nonlinear magnetohydrodynamic system where H , q and p denote, in turn, the magnetic field, velocity and pressure respectively. The density and magnetic permeability are assumed to be constant and have been scaled to unity. It is remarked that the system (1) has been previously shown in [37] to be analogous mutatis mutandis to that descriptive of the plane elastostatic deformation of neo-Hookean materials [38]. In the current case of a translational symmetry in the direction of a unit vector b , the magnetic induction equation (1) 3 admits a representation where A is the magnetic flux. Faraday's law (1) 4 then shows that The general solution of the latter may be decomposed into a particular solution A 0 and a function of any non-constant quantity which is convected by q , that is, It is this observation which, in part, motivates the following alignment of magnetohydrodynamic motions with those of fibre-reinforced fluids in which quantities which are convected by the fluid flow arise naturally.

Fibre-reinforced fluid flows
In order to make a connection with the steady planar motion of ideal fibre-reinforced fluids, we now assume that there exist notional inextensible "fibres" embedded in the fluid which are convected by q . If t denotes the unit tangent to the fibre lines then the existence of these fibre lines is encoded in the commutativity of the vector fields q ⋅ ∇ and t ⋅ ∇ , that is, Accordingly, the compatibility condition which guarantees the existence of a function obeying the pair of equations is satisfied. Geometrically, this means that constitutes arc length along the fibres which is convected by the fluid flow. It is noted that this is consistent with the assumption that the fibres are inextensible.
In order to proceed, it is convenient to introduce the directional derivatives where n denotes the principal (unit) normal to the fibres. If is the curvature of the fibres then the orthonormal pair (t, n) obeys the Serret-Frenet equations [39] and, analogously, with = divt being the negative of the curvature of the trajectories which are orthogonal to the fibres. Furthermore, t and n are represented by and the fibres and orthogonal trajectories are parametrised by s and n respectively so that a canonical orthogonal coordinate system (x(s, n), y(s, n)) on the plane is obtained. In terms of this parametrisation, it has been shown that the kinematic equation (3) subject to the continuity equation (1) 1 may be brought into the following compact form [24,25].
sn n + �� s = 0, The one-parameter family of fibres is given by (x, y)(s, n = const) . Conversely, any generic motion is encoded in this manner.

Remark 1
The relations (10) show that the metric of the plane and the gradient are given by and so that combination of the Serret-Frenet-type equations (4), (5) and (9) reveals that It is also observed that by virtue of the representation (8) of the velocity q.

The momentum equation
In order to evaluate the compatibility condition for the momentum equation (1) 2 which guarantees the existence of the pressure p, we first derive an explicit expression for the magnetic field H . Inspection of the directional derivative (12) reveals that A 0 = nc∕c is a particular solution of (2) so that the magnetic flux is given by where g( ) = f � ( ) . Accordingly, we arrive at the decomposition which turns out to be convenient in the following discussion.
The momentum equation (1) The vorticity may now be calculated. Thus, since b is constant, we conclude that wherein the Serret-Frenet-type equations (4), (5) and the expressions (11) for the curvatures and have been used. Two of the above terms may be combined to obtain Moreover, since is linear in c , we may write it as = 0 +c 1 so that, on the one hand, On the other hand, Accordingly, the compatibility condition (15) is explicitly given by where We therefore conclude that any solution of the overdetermined system of the third-order partial differential equations (7) and gives rise to an essentially unique magnetohydrodynamic motion.

Classes of solutions
Even though the preceding analysis demonstrates that considering cognate magnetohydrodynamic and fibrereinforced fluid motions is natural, it is by no means evident that the overdetermined system (7), (16) admits non-trivial solutions.

Rotationally symmetric motions
A reduction to ordinary differential equations is obtained via the ansatz Here, it should be observed that even though it is evident that the ansatz (17) is compatible with the constraint (16), it is, a priori, not obvious that any other symmetry reduction of (7) is consistent with (16). In the case (17), the third-order equation (7) may be integrated once to obtain the nonlinear pendulumtype equation where ĉ is the constant of integration. Moreover, the left-hand side of (16) vanishes identically, while the right-hand constitutes a differential equation for ĝ(s + n) = g( (s + n)) if c ≠ 0 . If the velocity and the magnetic field are aligned ( c = 0 ) then the function g is arbitrary.

+ ( ) � ( ) =ĉ ( ),
Fibre-reinforced fluid motions generated by the ansatz (17) have been discussed in detail in [24]. In particular, it is readily verified that the Cartesian coordinate vector r = (x y) T obeying the system (10) is given by The vector-valued additive constant of integration which has been suppressed is irrelevant unless ĉ = 0 in which case the constant has to be chosen in such a manner that the limit ĉ → 0 exists. In the case ĉ ≠ 0 , we conclude that r 2 is a function of s + n and, hence, (8) and (13) reveal that the velocity and magnetic field are of the form If we set aside the genesis of the latter then direct evaluation of the original magnetohydrodynamic system (1) leads to and In the aligned case associated with H 10 = 0 , the function H 0 is arbitrary. This corresponds to the function g being arbitrary. In fact, one may show that the magnetohydrodynamic motions associated with the current ansatz = (s + n) and ĉ ≠ 0 completely represent the rotationally symmetric class (18). The latter is characterised by concentric circular streamlines and magnetic field lines of spiral type. In the case H 01 = 0 , one retrieves a sub-class of magnetohydrodnamic motions for which the magnetic field lines and streamlines meet at a constant angle. The entire class of this nature has been classified in [30].

Compatibility
We now present a route which may be taken to reveal to which extent the differential equations (7) and (16) H 00 x 2 + y 2 + H 01 for H 10 ≠ 0.
share non-trivial solutions. The classification of all common solutions will be presented elsewhere. However, here, we will demonstrate that there exist common solutions which may be expressed in terms of Whittaker, Bessel and Airy functions.
We first introduce a function R according to It is observed that if R = R( ) then, locally, = (S(s) + n) which is, in turn, equivalent to the ansatz (17) since a suitable reparametrisation of the fibre lines leads to S = s . If we take into account that then it is readily seen that the third-order equation (7) becomes the second-order equation Secondly, if, for simplicity, we confine ourselves to the aligned case then (16) reduces to U = l( ) so that Once again, the parametrisation (19) reduces the order of this differential equation to obtain where Elimination of R n between (20) and (22) then produces the ordinary differential equation The general solution of the second-order equation (23) is of the form In principle, for any given triple ( , h, k) , the functions W i are determined by (23). Alternatively, we may regard the functions W i as being given so that insertion of (24) into (23) leads to a system of ordinary differential equations for ( , h, k) . Now, "taking the square" of (22) leads to a functional equation of the form where the functions N k and R k depend on N i and h, k, W i respectively. A classification of the classes of solutions of this functional equation may be obtained by standard techniques [40], leading to a system of constraints on the functions N i , W i and h, k.

Separation of variables
As an illustration of the procedure outlined in the preceding, we focus on the case W 0 = W 2 = 0 so that R is a product of a function of n and a function of . In this case, the functional equation (25) admits only one class of solutions. Instead of deriving this class explicitly, we observe that R being separable is equivalent to the ansatz and, hence, we may directly deal with the pair of differential equations (7) and (21). Thus, insertion of (26) into (7) produces while (21) becomes The latter implies that where c 1 and c 2 are arbitrary constants, so that l and g may be determined by and (25)  (28) N �� + c 1 N �2 = c 2 , l = c 2 c 2 − g 2 2 F z once F is known. It remains to separate the variables in the differential equation (27). Comparison with (28) shows that Since z = � F z , (29) 1 may be formulated as and, hence, The second equation (29) 2 then becomes which may be integrated to obtain or, equivalently, It is noted that c 4 may be appropriately parametrised so that the above remains valid in the limit c 1 → 0 . We therefore conclude that the functions F and are obtained via the linear system with

Motions involving special functions
The solution of the differential equation (28) may be expressed in terms of elementary functions.