On 2+1-dimensional fibre-reinforced fluid motions: a magnetohydrodynamics nexus

An intrinsic geometric decomposition is applied to provide via a canonical third-order nonlinear equation a connection between 2+1-dimensional fibre-reinforced motions of a fluid and a magnetohydrodynamic system. A Lagrange-type parametrisation is introduced, whereby both geometric and algebraic properties of certain non-steady magnetohydrodynamic motions may be established.


Introduction
In [1], a novel intrinsic geometric formulation was introduced in the context of the kinematic analysis of certain hydrodynamic motions. The representations contained therein have subsequently proved to have diverse applications in both nonlinear continuum mechanics and soliton theory. In magnetohydrodynamics, application has been made to isolate conducting motions in which the Maxwellian surfaces contain the streamlines as geodesics and the magnetic lines as parallels [2]. In soliton theory, the anholonomic representations of [1] have been adopted in [3] to derive, in a purely geometric manner, both the auto-Bäcklund transformation and a Lax pair for the nonlinear Schrödinger (NLS) equation. NLS soliton surfaces of both 'smoke ring' and spatial breather-type may be thereby constructed [4]. A connection between a previously well-studied class of spatial hydrodynamic motions and the integrable solitonic Heisenberg spin equation subject to a natural geometric constraint was uncovered in [5] by means of an intrinsic geometric representation. This novel link was subsequently elaborated upon in [6]. Later a connection was thereby made in [7] between this Heisenberg spin reduction and a classical problem in hydrodynamics posed by Gilbarg [8]. In recent work, such Heisenberg spin connection has been extended to spatial relativistic gasdynamics [9]. Motions subject to a complex-lamellar geometric constraint have been investigated via intrinsic anholonomic representations in both spatial gasdynamics [10] and magnetohydrodynamics [11]. In magnetohydrostatics [12], an integrable Heisenberg spin reduction obtains mutatis mutandis to that originally derived in hydrodynamics and was applied to determine geometric configurations in which the constant total magnetic pressure surfaces comprise nested tori foliated in accordance with the admitted constraint.
The connection in geometric terms between privileged motions of inextensible curves and modern soliton theory may be said to have its origin in the analysis by Da Rios [13] of the spatial motion of an isolated vortex filament embedded in an unbounded viscous liquid. Therein, what has come to be termed the localised induction approximation was invoked to derive a pair of coupled evolution equations for the curvature and torsion of the filament. These results were later to be incorporated in a survey by Levi-Civita [14]. Hasimoto [15], motivated by a geometric study in [16], was to establish that the classical Da Rios system is equivalent, without approximation, to the canonical nonlinear Schrödinger equation of modern soliton theory. Indeed, the latter NLS equation may be generated in a purely geometric manner via privileged binormal motions of an inextensible curve [5].
A kinematic study of the motion of fibre-reinforced fluids has its origin in the work of Spencer [17]. The model therein comprised an incompressible fluid with embedded inextensible fibre lines convected by the motions. That the motion of inextensible curves appears in both the context of the motion of fibre-reinforced fluids and soliton theory indicated the potential existence of integrable structure of solitonic type admitted by the kinematic system of [17]. This indeed has proved to be the case. Thus, in [18] a class of deformations inherent in the kinematic system of Spencer [17] has been linked to the AKNS inverse scattering scheme via the integrable modified Korteweg-de-Vries (mkdV) solitonic hierarchy.
In [19,20], the geometric formalism previously adopted in hydrodynamics, magnetohydrostatics, magnetohydrodynamics, and soliton theory was used to investigate the kinematic conditions attendant upon the motion of a fibre-reinforced fluid. Kinematically admissible spatial motions were thereby obtained in which the fibres are geodesic windings on nested toroidal surfaces. In the case of purely planar motion, it was shown that the kinematic relations may be encapsulated in a novel third-order nonlinear equation. Remarkably, this admits a particular reduction to a solitonic system which is related to the classical sine-Gordon equation as originally obtained in the context of pseudospherical surfaces by Bour [21]. The kinematic conditions in this reduction admit both a Bäcklund transformation [22] and a duality property. In [23], the reduction of the kinematic system of [17] to a third-order nonlinear equation has been extended to 2+1-dimensional motions and underlying integrable solitonic structure again isolated. The geometry and integrability of steady motions of an ideal fibre-reinforced fluid in a curved stratum have been the subject of [24].
In [25], it has recently been established that the canonical third-order nonlinear equation derived via intrinsic geometric representations in [19,20] may likewise be shown to arise in steady two-dimensional magnetohydrodynamics subject to analogous kinematic conditions. This connection was used to isolate novel classes of time-independent conducting motions. Here, this nexus between the theory of the motion of fibrereinforced liquids and magnetohydrodynamics is extended to 2+1-dimensional motions. In that connection, a geometric Lagrangian-type parameterisation is introduced. Geometric and algebraic properties of the timedependent magnetohydrodynamic motions are described in detail. In overall terms, the preceding attests to the utility of the imposition of natural geometric constraints to isolate a hidden integrable structure in nonlinear continuum mechanics [26].

Planar fibre-reinforced fluid motions
The ideal fibre-reinforced fluid model considered here consists of an incompressible liquid which is inextensible along 'fibre' lines which occupy the volume of the fluid by which they are convected [17]. Thus, if a generic fibre direction is characterised by a unit vector t and q denotes the fluid velocity, then the kinematic condition which encodes both the convection requirement and inextensibility of the fibres is given by This evolution equation is augmented by the continuity equation which follows from the incompressibility of the fluid.

Geometric decomposition and reformulation
The equation of motion (1) may be regarded as the compatibility condition which guarantees the existence of a function ρ obeying the pair of equations Geometrically, this means that ρ constitutes the arc length along the fibres which is convected by the fluid flow. It is observed that this is consistent with the assumption that the fibres are inextensible. It is therefore 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 [27] (5) and, analogously, with θ = div t being the negative of the curvature of the trajectories which are orthogonal to the fibres. Furthermore, t and n may be parametrised by an angle ϕ according to The existence of the fibre lines gives rise to a natural time-dependent orthogonal coordinate system (x(s, n, τ ), y(s, n, τ )) on the plane, where τ = t and the fibres and orthogonal trajectories are parametrised by s and n, respectively. Remarkably, it has been shown [23] that, in terms of the independent variables s, n, and τ , the kinematic equations (1) and (2) may be reformulated as a third-order partial differential equation in which time only enters parametrically. This is the content of the following theorem.

Theorem 1 Given a solution (ρ, ψ(ρ, n)) of the partial differential equation
wherein ρ may depend parametrically on τ and a solution of the governing equations (1), (2) is given by where the angle ϕ parametrising t and n is determined by the compatible system and the velocity components v and w read with the parametrisation of the Eulerian coordinates being obtained via integration of the compatible system and t = τ . At any instant t, the one-parameter family of fibres labelled by n is parametrised by s according to It is observed that the relations (12) imply that the metric of the plane and the gradient are given by and Hence, insertion of the parametrisation (7) into the Serret-Frenet-type equations (5) and (6) reveals, on use of (10), that

Lagrangian coordinates
It is readily shown [23] that the convective derivative adopts the compact form in terms of the adapted coordinates s, n, and τ . If one applies the latter to the Eulerian coordinate vector r, then one obtains which is seen to be consistent with the components (11) of the velocity q. This observation suggests introducing the independent variables leading to or, equivalently, In particular, the equation of motion (1) is then simply represented by the commutativity and the particle lines of the fluid are given by r(s = const, n = const, t) since In summary, the above confirms that not only does ρ = s represent the arc length along the fibres which are convected by the fluid, it also labels the particles on the fibres n = const as they deform as a function of time.
Thus, the independent variables s and n constitute material (Lagrangian) coordinates. The above change of independent variables is naturally accompanied by the introduction of a new dependent variable given by The pair of equations (10) determining ϕ then adopts the form so that the associated compatibility condition (ϕ s ) n = (ϕ n ) s produces the second-order equation which constitutes the third-order equation (8) formulated in terms of the new (in)dependent variables. Once R is known, the function ρ is obtained by solving the first-order ordinary differential equation

The magnetohydrodynamic equations
Here, the concern is with the 2+1-dimensional reduction of the non-dissipative magnetohydrodynamic system div q = 0, (27.4) where H and p denote, in turn, the magnetic field and pressure. The density and magnetic permeability are assumed to be constant and have been scaled to unity. The solution of the magnetic induction and Faraday equations (27.3, 27.4) is readily seen to be where the magnetic flux A is convected by the fluid, that is, The (constant) unit vector b is orthogonal to the plane in which the motion takes place. If we reformulate the momentum equation (27.2) as then its compatibility condition guaranteeing the existence of the pressure p may be formulated as [25] with ω q and ω H being the vorticity and current density, respectively. In [25], it has been demonstrated in the steady case that it is natural to assume that there exist notional inextensible 'fibres' embedded in the fluid which are convected by the fluid flow. It turns out that this is also the case in the current setting so that we now require q to obey the kinematic equations (1) and (2) for a unit vector field t to be determined. An immediate consequence of this assumption is that the magnetic flux A is a function of ρ and n since both quantities are convected by the fluid. Here, for simplicity, we confine ourselves to the case of A being a function of ρ only, that is, On use of the decomposition (14), we conclude that so that The compatibility condition (32) may therefore be written as In order to find explicit expressions for the vorticity ω q and the current density ω H , we note that the curl of any vector field F = by virtue of the Serret-Frenet-type relations (5) and (6). Here, we have made the choice b := t × n. It is recalled that the curvature κ and divergence θ = div t of the fibre lines are given by (15.1,2). In the case of the fluid velocity q, the expression may be evaluated on use of the relations (10), (11), and (12). Remarkably, a straightforward but somewhat lengthy calculation reveals that the terms involving x τ and y τ cancel each other and the vorticity admits a compact representation, namely On the other hand, it may also be shown that, as in the steady case [25], In summary, the motions considered here are governed by common solutions (ρ, ψ(ρ, n), B(ρ)) of the thirdorder equation (8) together with the compatible pair (10) and the compatibility condition (36) with ω q and ω H given by (39) and (40), respectively.

Lagrangian coordinates
The transformation laws (20) show that, in terms of the Lagrangian coordinates (s, n, t), the compatibility condition (36) admits the conservative form In fact, by virtue of (39), we obtain Here, it is emphasised that ω H may be expressed entirely in terms of the new (in)dependent variables since does not depend on the quantity ρ s . Accordingly, if then the magnetic flux A(s) is unconstrained, and the compatibility condition (42) reduces to Furthermore, since the second-order equation (25) is a consequence of the pair (24), all solutions of the form (44) of the hybrid fibre-reinforced fluid-magnetohydrodynamic system considered here are encapsulated in the overdetermined system Prior to an analysis of the solutions of the above system, it is enlightening to derive some general properties of the motions encoded in the ansatz (44).

Geometric and algebraic properties
Integration of (26) leads to a relation of the form F (ρ(s, n, τ ), τ ) = n +s(s, τ ) wheres is a function of integration. The non-uniqueness of the coordinate s may then be exploited to sets = s without loss of generality so that, locally, In this case, the third-order equation (8) may be integrated once to obtain the nonlinear oscillator equation which, in turn, admits the first integral Here, the variable τ enters parametrically.
Remarkably, as pointed out in the steady case [19], the four equations (12) for the Eulerian coordinates x and y may be integrated explicitly to obtain for μ = 0, where the additive function of integration depending on τ has been neglected. It may be retrieved by applying a standard symmetry of the original (unconstrained) governing equations (1), (2), and (27.1), corresponding to the transition to an, in general, non-inertial frame. Hence, we conclude that the fibre-distribution is rotationally symmetric since However, in contrast to the steady case [19], in general, this does not apply to the fluid flow so that the motions considered here are genuinely 2+1-dimensional. Interestingly, the magnetic field is rotationally symmetric with the field lines constituting concentric circles. This may be concluded from the relations Furthermore, if the right-hand side of is bounded above and/or below at any fixed time t = τ , then the fluid is confined to a region which is bounded by one or two circles centred at the origin. In fact, it is shown in the following Section that the motion takes place within an annulus, the size of which changes as a function of time.

The solution (ϕ, R, ψ)
The compatibility conditions associated with the system (46.1-3) produce an overdetermined system for the function R(s, t) and ψ(s). In the case ψ s = 0, one obtains a system for R only, namely Now, the conditions (55.1,2) and (55.1,2) supplemented by their tand s-derivative, respectively, allow one to eliminate the quantities R ss , R tt and R sst and derive the necessary condition Hence, since R is real, we are led to the strong constraint R s R t = 0, corresponding to 'elementary' flows which we set aside.
In the case ψ s = 0, the compatibility of (46.1-3) leads to The overdetermined system of compatibility conditions for (46.1-3) may be closed in the sense of Riquier-Janet theory [28] by adding all integrability conditions. Specifically, the rif algorithm recorded in [29] and implemented in the computer algebra software Maple [30] may be successfully employed to show that the corresponding involutive system includes the constraint 2 + 4 = 0 in the case R tt = 0. Accordingly, the existence of real solutions requires that Since R n = 0, it is easy to solve the system (46.1-3) subject to (58). Indeed, after the dependence on n and t has been separated according to one is left with the ordinary differential equations If we take into account the non-uniqueness in the definition of the variables s, n, and ψ (subject to the preservation of (48)), the arbitrariness of the 'origin' of time t and that, without loss of generality, we may scale x, y, and t and rotate the (x, y)-plane in which the motion takes place, then the corresponding constants of integration in the solution of the system (60) may be neglected, and it may be assumed that α = 1 and β = 0. Hence, we obtain so that the solution of the system (46.1-3) is given by as may be directly verified. Integration of the differential equation (26) results in It is noted that even though this quantity changes its nature for τ < 1 and in the limit t → 1, it remains real. By construction, ρ is a solution of the differential equation (49) with μ = τ , that is, the 'stationary' double sine-Gordon equation [31] ρ zz + sin ρ cos ρ + τ sin ρ = 0.
It is observed that the above differential equation is invariant under τ → −τ , ρ → ρ + π which is consistent with (63) if the function of integration z 0 is taken into account. It is seen below that this invariance is reflected by a time reversal symmetry of the fluid flow. The curvature (15.1,2) of the fibre lines is given by so that the fibres constitute circles of unit radius. In fact, it is natural to consider the semicircular fibres r(s, n, t), s ∈ [0, π] which foliate the region occupied by the fluid at any time t = t. Indeed, evaluation of (51) delivers so that Accordingly, the fluid flow is confined to an annulus, namely the region between two concentric circles of radii given by r(s = 0, n, t) and r(s = π, n, t) for any fixed time t. It is evident that these two circles are convected by the fluid as required. As t → ∞, the two circles approach each other and coalesce into a unit circle centred at the origin.

The fluid flow
Since the mapping (x, y) → (s, n) is multi-valued if n is not restricted, the fluid velocity q = r t given by is multi-valued. In fact, on taking into account that s ∈ [0, π], inversion of (66) leads to so that the multiplicity of the coordinate transformation is encoded in the angle arg(x + iy). If we choose its principal value contained in the interval (−π, π], then we obtain a vector field q which is single-valued but discontinuous on the negative x-axis. This vector field is displayed in Figure 1 at time t = 2. The latter reveals the existence of a 'hyperbolic' stagnation point. Indeed, it is readily shown that q = 0 if and only if Accordingly, a (moving) stagnation point exists if |t| ≥ 1 and is located at Since the fibre lines are convected by the fluid, it is admissible to restrict the time-dependent region in which the motion takes place to a region which is bounded by two fibre lines corresponding to n = n min and n = n max and the concentric circles enclosing the annulus. For instance, if one chooses n min = − 1 4 and n max = 3 4 , then, at the initial time t = 2, the region indicated in Figure 1 is obtained in which the fluid velocity (69) is single-valued. Since this region contains the stagnation point, the fluid particles which are close to the fibre line n = n min move towards the unit circle in a clockwise manner, while the motion of the fluid particles close to the fibre line n = n max is anticlockwise. Figure 2 depicts a sample of the fibre lines at times t = 2, 3, 4, 5, 6.
The blue dotted lines are the curves s = const. Accordingly, the distance between (the particles on) any such curves along the fibre lines is constant. At time the fibre lines n = n min and n = n max meet, and the fluid particles begin to mix. It is noted that the same phenomenon occurs if the stagnation point is initially not part of the region in which the motion takes place. Eventually, the faster part of the fluid catches up with the slower part, and the particles begin to mix. Figure 3 depicts the fibre line n = 3 4 (red) at various times t ∈ [2, 6] and the associated particle lines s = const (blue).

A regular solution
As t decreases, the outer radius of the annulus increases and becomes infinite as t → 0. The inner radius decreases until it vanishes at t = 1, at which point the region in which the motion takes place becomes a disk. If t is reduced further, then the inner radius increases again and becomes infinite as t → 0. At t = 0, the fluid velocity also becomes singular. In order to obtain a better understanding of the properties of the flow as t → 0, it is convenient to introduce a non-vanishing function of integration in (51), corresponding to the transition to a non-inertial frame, which ensures that the annulus in which the motion takes place does not move to infinity as t → 0. As a consequence, the centre of the annulus becomes time dependent and approaches infinity. Similarly, at any fixed time, the magnetic field lines remain concentric circles but their centre moves since it coincides with the centre of the annulus as in the original solution. Specifically, we consider the transition y → y − t −1 , q (y) → q (y) + t −2 so that the solution of the governing equations is now given by and which is regular at t = 0. Indeed, at t = 0, we obtain x = sin s − n, y = cos s, q (x) = −n cos s, q (y) = n sin s − 1 2 n 2 (76) so that the fluid flow is confined to a horizontal strip bounded by the lines y = ±1. The multi-valued connection between the Eulerian coordinates x, y and the Lagrangian coordinates s, n is now given by (70) with y replaced by y + t −1 , that is, The fluid velocity q at t = 0 is depicted in Figure 4 where, once again, the principal value of arg has been taken. Figure 4 reveals two features of the fluid flow. On the one hand, there exists a fibre line on which the fluid velocity vanishes. In fact, inspection of the normalised vector field shows that the normalised velocity becomes tangent but changes its direction across this 'stagnation line'. This privileged fibre line is stationary and may be interpreted as a fixed boundary past which the fluid flows. Indeed, it is evident that q = 0 at n = 0 for all t and s ∈ [0, π]. Figure 5 displays a sample of the fibre lines at times t = 0, 1 2 , 1, 3 2 , 2, 5 2 . The orange semicircle x 2 + y 2 = 1, x ≥ 0 represents the stationary fibre line n = 0. The motion is chosen to take place in the region bounded by the fibres n = n min = −1 and n = n max = 1. As time progresses, the strip becomes an annulus of finite radius until it degenerates at t = 1 to a disk of radius 2. The inner radius then starts increasing and the outer radius continuous to decrease. Figure 6 depicts the fibre line n = 1 (red) at various times t ∈ [0, 5 2 ] and the associated particle lines s = const (blue). Secondly, Figure 4 indicates that the points on the fibre line n = 0 are not the only stagnation points. In fact, the 'elliptic' stagnation point at t = 0 is located at (−1, 0) since q(s = π 2 , n = 2) = 0 by virtue of (76). In general, excluding n = 0, the condition q = 0 leads to the pair which immediately implies that we require that n > 0 by virtue of s ∈ [0, π] and any additional stagnation point is located on the other half of the unit circle x 2 + y 2 = 1, x ≤ 0 since it is readily verified that x = − sin s, y = cos s.
If we focus on non-negative times, then its solution may be parametrised according to so that In Figure 7, the quantities t, n, and s are plotted against m. It is observed that, at time t = 0, there exists only one stagnation point (apart from those on the fibre line n = 0). As t increases, only one stagnation point remains present until t reaches a certain value t 2 at which another stagnation point appears. The latter subsequently splits into two stagnation points. As t approaches 1, more and more stagnation points emerge, and there exists an infinite number of coalescing stagnation points at t = 1 located at (0, −1) where all fibre lines meet. As t increases further, the number of stagnation points becomes finite again and decreases until only one stagnation point is left at some time t = t max . At any later time, no stagnation points exist. The fibre line n = 1 (red) at various times t ∈ [0, 5 2 ] and the associated particle lines s = const (blue).

Fig. 7
The values of t, n, s corresponding to additional stagnation points parametrised by m which do not lie on the fibre line n = 0 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funding Open Access funding enabled and organized by CAUL and its Member Institutions