A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction
A 2+1-dimensional anisentropic magnetogasdynamic system with a polytropic gas law is shown to admit an integrable elliptic vortex reduction when γ=2 to a nonlinear dynamical subsystem with underlying integrable Hamiltonian-Ermakov structure. Exact solutions of the magnetogasdynamic system are thereb...
Gespeichert in:
| Veröffentlicht in: | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Datum: | 2012 |
| Hauptverfasser: | , |
| Format: | Artikel |
| Sprache: | English |
| Veröffentlicht: |
Інститут математики НАН України
2012
|
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/148449 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction / H. An, C. Rogers // Symmetry, Integrability and Geometry: Methods and Applications. — 2012. — Т. 8. — Бібліогр.: 22 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
nasplib_isofts_kiev_ua-123456789-148449 |
|---|---|
| record_format |
dspace |
| spelling |
An, H. Rogers, C. 2019-02-18T12:44:36Z 2019-02-18T12:44:36Z 2012 A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction / H. An, C. Rogers // Symmetry, Integrability and Geometry: Methods and Applications. — 2012. — Т. 8. — Бібліогр.: 22 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 34A34; 35A25 DOI: http://dx.doi.org/10.3842/SIGMA.2012.057 https://nasplib.isofts.kiev.ua/handle/123456789/148449 A 2+1-dimensional anisentropic magnetogasdynamic system with a polytropic gas law is shown to admit an integrable elliptic vortex reduction when γ=2 to a nonlinear dynamical subsystem with underlying integrable Hamiltonian-Ermakov structure. Exact solutions of the magnetogasdynamic system are thereby obtained which describe a rotating elliptic plasma cylinder. The semi-axes of the elliptical cross-section, remarkably, satisfy a Ermakov-Ray-Reid system. This paper is a contribution to the Special Issue “Geometrical Methods in Mathematical Physics”. The full collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction Article published earlier |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| title |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction |
| spellingShingle |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction An, H. Rogers, C. |
| title_short |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction |
| title_full |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction |
| title_fullStr |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction |
| title_full_unstemmed |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction |
| title_sort |
2+1-dimensional non-isothermal magnetogasdynamic system. hamiltonian-ermakov integrable reduction |
| author |
An, H. Rogers, C. |
| author_facet |
An, H. Rogers, C. |
| publishDate |
2012 |
| language |
English |
| container_title |
Symmetry, Integrability and Geometry: Methods and Applications |
| publisher |
Інститут математики НАН України |
| format |
Article |
| description |
A 2+1-dimensional anisentropic magnetogasdynamic system with a polytropic gas law is shown to admit an integrable elliptic vortex reduction when γ=2 to a nonlinear dynamical subsystem with underlying integrable Hamiltonian-Ermakov structure. Exact solutions of the magnetogasdynamic system are thereby obtained which describe a rotating elliptic plasma cylinder. The semi-axes of the elliptical cross-section, remarkably, satisfy a Ermakov-Ray-Reid system.
|
| issn |
1815-0659 |
| url |
https://nasplib.isofts.kiev.ua/handle/123456789/148449 |
| citation_txt |
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System. Hamiltonian-Ermakov Integrable Reduction / H. An, C. Rogers // Symmetry, Integrability and Geometry: Methods and Applications. — 2012. — Т. 8. — Бібліогр.: 22 назв. — англ. |
| work_keys_str_mv |
AT anh a21dimensionalnonisothermalmagnetogasdynamicsystemhamiltonianermakovintegrablereduction AT rogersc a21dimensionalnonisothermalmagnetogasdynamicsystemhamiltonianermakovintegrablereduction AT anh 21dimensionalnonisothermalmagnetogasdynamicsystemhamiltonianermakovintegrablereduction AT rogersc 21dimensionalnonisothermalmagnetogasdynamicsystemhamiltonianermakovintegrablereduction |
| first_indexed |
2025-11-24T11:44:40Z |
| last_indexed |
2025-11-24T11:44:40Z |
| _version_ |
1850463976560787456 |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 8 (2012), 057, 15 pages
A 2+1-Dimensional Non-Isothermal
Magnetogasdynamic System.
Hamiltonian–Ermakov Integrable Reduction?
Hongli AN † and Colin ROGERS ‡§
† College of Science, Nanjing Agricultural University, Nanjing 210095, P.R. China
E-mail: kaixinguoan@163.com
‡ School of Mathematics and Statistics, The University of New South Wales,
Sydney, NSW 2052, Australia
§ Australian Research Council Centre of Excellence for Mathematics & Statistics
of Complex Systems, School of Mathematics, The University of New South Wales,
Sydney, NSW2052, Australia
E-mail: c.rogers@unsw.edu.au
URL: http://web.maths.unsw.edu.au/~colinr/
Received May 27, 2012, in final form August 02, 2012; Published online August 23, 2012
http://dx.doi.org/10.3842/SIGMA.2012.057
Abstract. A 2+1-dimensional anisentropic magnetogasdynamic system with a polytropic
gas law is shown to admit an integrable elliptic vortex reduction when γ = 2 to a nonlinear
dynamical subsystem with underlying integrable Hamiltonian–Ermakov structure. Exact
solutions of the magnetogasdynamic system are thereby obtained which describe a rotating
elliptic plasma cylinder. The semi-axes of the elliptical cross-section, remarkably, satisfy
a Ermakov–Ray–Reid system.
Key words: magnetogasdynamic system; elliptic vortex; Hamiltonian–Ermakov structure;
Lax pair
2010 Mathematics Subject Classification: 34A34; 35A25
1 Introduction
Neukirch et al. [6, 7, 8] have investigated 2+1-dimensional magnetogasdynamic systems via
a solution approach in which the nonlinear acceleration terms in the Lundquist momentum
equation either vanish or are conservative. By contrast, in recent work [12, 20] an elliptic
vortex ansatz was adopted in 2+1-dimensional isothermal magnetogasdynamics and underlying
integrable Ermakov–Ray–Reid structure was isolated. In [20], magnetogasdynamic pulsrodon-
type solutions were constructed analogous to those originally derived in elliptic warm-core theory
in [13]. The pulsrodons describe an elliptical plasma cylinder bounded by a vacuum. The time-
dependent semi-axes of the elliptical cross-section of the cylinder were shown to be governed by
an integrable Hamiltonian Ermakov system.
The present work concerns an extension of that of [12] to a non-isothermal rotating mag-
netogasdynamic version of a spinning non-conducting gas cloud system with origin in work of
Ovsiannikov [9] and Dyson [2]. A nonlinear dynamical subsystem is derived which is again
remarkably, shown to have integrable Hamiltonian Ermakov–Ray–Reid structure. Moreover,
a Lax pair for the dynamical system is constructed.
?This paper is a contribution to the Special Issue “Geometrical Methods in Mathematical Physics”. The full
collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html
mailto:kaixinguoan@163.com
mailto:c.rogers@unsw.edu.au
http://web.maths.unsw.edu.au/~colinr/
http://dx.doi.org/10.3842/SIGMA.2012.057
http://www.emis.de/journals/SIGMA/GMMP2012.html
2 H.L. An and C. Rogers
2 The magnetogasdynamic system
Here, we consider a 2+1-dimensional anisentropic magnetogasdynamic system incorporating
rotation, namely,
∂ρ
∂t
+ div(ρq) = 0, (2.1)
ρ
[
∂q
∂t
+ (q · ∇)q + fk× q
]
− µ curlH×H +∇p = 0, (2.2)
divH = 0, (2.3)
∂H
∂t
= curl(q×H), (2.4)
∂S
∂t
+ q · ∇S = 0, (2.5)
where the velocity q and magnetic field H are given by
q = ui + vj, H = ∇A× k + hk (2.6)
respectively, while the gas law adopts the polytropic form
S = − ln ρ+
1
γ − 1
lnT, γ 6= 1 (2.7)
with
p = ρT. (2.8)
In the above, the magneto-gas density ρ(x, t), pressure p(x, t), entropy S(x, t), temperatu-
re T (x, t) and magnetic flux A(x, t) are all assumed to be dependent only on x = xi + yj
and time t. In addition, f is the Coriolis constant, µ the magnetic permeability and h(x, t) the
transverse component of the magnetic field.
Insertion of the representation (2.6) into Faraday’s law (2.4) produces the convective con-
straint
∂A
∂t
+ q · ∇A = 0 (2.9)
together with
∂h
∂t
+ div(hq) = 0,
which holds automatically if we set
h = λρ, λ ∈ R.
Here, a novel two-parameter (m,n) pressure-density ansatz
p = ε0(t)ρ
2 + ε1(t)ρ
n + ε2(t)ρ
m. (2.10)
is introduced. In the magnetogasdynamic study of [7], a relation p ∼ ρ was adopted, while
pressure-density relations of the type p ∼ ρ2 arise in astrophysical contexts [21]. A parabolic
pressure-density law was recently employed in 2+1-dimensional isothermal magnetogasdynamics
in [20] and pulsrodon-type solutions were isolated.
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 3
In the present non-isothermal context, substitution of (2.10) in (2.8) produces the temperature
distribution
T = ε0(t)ρ+ ε1(t)ρ
n−1 + ε2(t)ρ
m−1, (2.11)
while the entropy distribution adopts the form
S = − ln ρ+
1
γ − 1
ln
(
ε0(t)ρ+ ε1(t)ρ
n−1 + ε2(t)ρ
m−1). (2.12)
The energy equation now requires that
(ρt + q · ∇ρ)
[
ε0 + (n− 1)ε1ρ
n−2 + (m− 1)ε2ρ
m−2
(γ − 1)(ε0ρ+ ε1ρn−1 + ε2ρm−1)
− 1
ρ
]
+
ε̇0ρ+ ε̇1ρ
n−1 + ε̇2ρ
m−1
(γ − 1)(ε0ρ+ ε1ρn−1 + ε2ρm−1)
= 0,
whence, on use of the continuity equation (2.1)
−divq
[
(2− γ)ε0ρ+ (n− γ)ε1ρ
n−1+ (m− γ)ε2ρ
m−1]+ ε̇0ρ+ ε̇1ρ
n−1+ ε̇2ρ
m−1 = 0 (2.13)
that is
1
γ − 1
Ṫ
T
[
(2− γ)ε0ρ+ (n− γ)ε1ρ
n−1 + (m− γ)ε2ρ
m−1]+ ε̇0ρ+ ε̇1ρ
n−1 + ε̇2ρ
m−1 = 0.
On substitution of (2.6) and (2.10) into the momentum equation (2.2), it is seen that
∂q
∂t
+ q · ∇q + fk× q +
1
ρ
[
µ(∇2A)(∇A) + ε1∇ρn
]
+ (µλ2 + 2ε0)∇ρ+
ε2
ρ
∇ρm = 0
together with
Ayρx −Axρy = 0,
so that
A = A(ρ, t).
Attention is here restricted to the separable case
A = Φ(ρ)Ψ(t)
whence, substitution into (2.9) and use of the continuity equation yields
Ψ̇(t)
Ψ(t)
= ρ
Φ′(ρ)
Φ(ρ)
divq.
Here, we proceed with
Φ = ρn,
where n is the parameter involving in the relation (2.11), so that
divq =
(
1
n
)
Ψ̇
Ψ
(2.14)
4 H.L. An and C. Rogers
and
A = ρnΨ(t).
Hence, as in the case of the spinning non-conducting gas cloud analysis of Ovsiannikov [9] and
Dyson [2], the divergence of the velocity is dependent only on time. Moreover, the relation (2.13)
shows that
− 1
n
Ψ̇
Ψ
[
(2− γ)ε0ρ+ (n− γ)ε1ρ
n−1 + (m− γ)ε2ρ
m−1]+ ε̇0ρ+ ε̇1ρ
n−1 + ε̇2ρ
m−1 = 0
and it is observed that this condition holds identically with
ε0 = α0Ψ
2−γ
n , (2.15)
ε1 = α1Ψ
n−γ
n , (2.16)
ε2 = α2Ψ
m−γ
n , (2.17)
where αi (i = 0, 1, 2) are arbitrary constants of integration. In addition, the isentropic con-
dition (2.5) together with the polytropic gas law (2.7) and the continuity equation (2.1) show
that
divq =
1
1− γ
Ṫ
T
whence, on use of (2.11),
(2− γ) divq =
ε̇0
ε0
, (n− γ) divq =
ε̇1
ε1
, (m− γ) divq =
ε̇2
ε2
.
It is seen that in view of (2.14), these relations are indeed consistent with (2.15)–(2.17).
In summary, the magnetogasdynamic system now reduces to consideration of the nonlinear
coupled system
∂ρ
∂t
+ div(ρq) = 0,
Ψ̇
Ψ
= n divq, (2.18)
∂q
∂t
+ q · ∇q + fk× q +
1
ρ
(
µΨ2∇2ρn + ε1
)
∇ρn + (µλ2 + 2ε0)∇ρ+
ε2m
m− 1
∇ρm−1 = 0,
where m 6= 1, together with the additional conditions (2.15)–(2.17). The inherent nonlinearity
of the system (2.18) remains a major impediment to analytic progress. It is noted also that the
system (2.18)1,3 is overdetermined since it is implicitly constrained by the requirement (2.18)2
that divq be a function of t only.
3 An elliptic vortex ansatz. A dynamical system reduction
Here, integrable nonlinear dynamical subsystems of the magnetogasdynamic system (2.18) are
sought via an elliptic vortex ansatz of the type
q = L(t)x + M(t), ρ =
(
xTE(t)x + ρ0
)m−1
, m 6= 1, x =
(
x− q̄(t)
y − p̄(t)
)
(3.1)
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 5
with
L(t) =
(
u1(t) u2(t)
v1(t) v2(t)
)
, E(t) =
(
a(t) b(t)
b(t) c(t)
)
, M(t) =
(
˙̄q(t)
˙̄p(t)
)
. (3.2)
Insertion of (3.1) into the continuity equation yieldsȧḃ
ċ
+
2u1 + (m− 1)(u1 + v2) 2v1 0
u2 m(u1 + v2) v1
0 2u2 2v2 + (m− 1)(u1 + v2)
ab
c
= 0 (3.3)
together with
ρ̇0 + ρ0(m− 1)(u1 + v2) = 0, (3.4)
whence
ρ0 = const Ψ(1−m)/n. (3.5)
If we now proceed with
n = m− 1 (3.6)
together with
2ε0 + µλ2 = 0 (3.7)
and
ε1 + 2µΨ2(a+ c) = 0 (3.8)
then it is seen that (2.18)3 reduces to
∂q
∂t
+ q · ∇q + fk× q +
m
m− 1
ε2∇ρm−1 = 0. (3.9)
The relation (3.7) implies that ε̇0 = 0 whence (2.15) shows that the adiabatic index γ = 2,
while (3.8) and (2.16) together require
a+ c = −α1
2µ
Ψ−(n+2)/n = −α1
2µ
Ψ(1+m)/(1−m), n 6= 1.
Substitution of (3.1) into (3.9) now gives
u̇1
u̇2
v̇1
v̇2
+
(
LT −fI
fI LT
)
u1
u2
v1
v2
+ 2ε2
m
m− 1
a
b
b
c
= 0 (3.10)
augmented by the auxiliary linear equations
¨̄p+ f ˙̄q = 0, ¨̄q − f ˙̄p = 0. (3.11)
It is noted that the relation (2.18)2 together with (3.6) shows that
Ψ̇ = (m− 1)(u1 + v2)Ψ.
6 H.L. An and C. Rogers
While ρ0 is given in terms of Ψ via (3.5). The constraints (3.7) and (3.8) are to be adjoined and
their admissibility will be examined subsequently.
In what follows, it proves convenient to proceed in terms of new variables, namely
G = u1 + v2, GR =
1
2
(v1 − u2), GS =
1
2
(v1 + u2), GN =
1
2
(u1 − v2),
B = a+ c, BS = b, BN =
1
2
(a− c).
These quantities were originally introduced in a hydrodynamic context (see e.g. [13]). Therein,
G and GR correspond, in turn, to the divergence and spin of the velocity field, while GS and GN
represent shear and normal deformation rates. The system (3.3) and (3.4) together with (3.10)
now reduces to the nonlinear dynamical system
ρ̇0 + (m− 1)ρ0G = 0,
Ḃ +mBG+ 4(BNGN +BSGS) = 0,
ḂS +mBSG+BGS − 2BNGR = 0,
ḂN +mBNG+BGN + 2BSGR = 0,
Ġ+
1
2
G2 + 2
(
G2
N +G2
S −G2
R
)
− 2fGR + 2
ε2m
m− 1
B = 0,
ĠN +GGN − fGS + 2
ε2m
m− 1
BN = 0,
ĠS +GGS + fGN + 2
ε2m
m− 1
BS = 0,
ĠR +GGR +
1
2
fG = 0
(3.12)
together with
Ψ̇ = (m− 1)ΨG. (3.13)
It is observed that the introduction of the pressure-density parameters (m,n) and ε2 leads to
a generalisation of the nonlinear dynamical systems obtained in [12, 13, 14, 20].
If we now introduce the quantity Ω via
G =
2Ω̇
Ω
then (3.12)1 and (3.12)8 show, in turn, that
ρ0 =
cI
Ω2(m−1) (3.14)
and
GR =
c0
Ω2
− 1
2
f. (3.15)
While the relation (3.13) yields
Ψ = νΩ2(m−1), (3.16)
where c0, cI and ν denote arbitrary constants of integration.
Two conditions which are key to the subsequent development and which may be established
by appeal to the original system (3.12) are now recorded. These represent extensions of results
obtained in a hydrodynamic context [13, 14].
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 7
Theorem 1.
Ṁ + (m+ 1)GM = 0, Q̇+ (m+ 1)GQ = 0,
where
M = 2(BNGS −BSGN )−B
(
GR +
1
2
f
)
, 4 =
1
4
B2 −B2
S −B2
N ,
Q = −B
(
G2
S +G2
N +G2
R +
1
4
G2
)
+ 4GR(BNGS −BSGN )
+ 2G(BSGS +BNGN ) + 4
ε2m
m− 1
4− 4
m
m− 1
4Ωm−1
∫
ε̇2Ω
1−mdt.
New Ω-modulated variables involving the pressure-density parameter m are now introduced
according to
B̄ = Ω2mB, B̄S = Ω2mBS , B̄N = Ω2mBN , ḠS = Ω2GS , ḠN = Ω2GN , (3.17)
whence the dynamical system (3.12) reduces to
˙̄B +
4(B̄N ḠN + B̄SḠS)
Ω2
= 0,
˙̄BS + fB̄N +
B̄ḠS − 2c0B̄N
Ω2
= 0,
˙̄BN − fB̄S +
B̄ḠN + 2c0B̄S
Ω2
= 0,
˙̄GS + fḠN +
2ε2m
m− 1
B̄S
Ω2(m−1) = 0,
˙̄GN − fḠS +
2ε2m
m− 1
B̄N
Ω2(m−1) = 0
(3.18)
augmented by the relations (3.14) and (3.15) together with a nonlinear equation for Ω, namely
Ω3Ω̈ +
1
4
f2Ω4 + Ḡ2
N + Ḡ2
S − c20 +
ε2m
m− 1
B̄
Ω2(m−2) = 0. (3.19)
The reduced dynamical system (3.18) together with (3.19) and the constraints given by (2.15)–
(2.17) and (3.7), (3.8) will now be examined in detail.
Thus, if we turn to the expressions for ε0, ε1 and ε2 as given by (2.15)–(2.17), it is seen
immediately that consistency of (2.15) and (3.7) requires that the adiabatic index γ = 2. Further,
comparison of the expressions for ε1 in (2.16) and (3.8) now yields
α1Ψ
n−2
n + 2µΨ2(a+ c) = 0,
where the relations (3.6), (3.16) and (3.17) combine to show that
α1ν
m−3
m−1 + 2µν2Ω2B̄ = 0,
whence
ν = 0 or Ω2B̄ = −α1
2µ
ν
1+m
1−m := δ.
In the former case, by virtue of (3.16), the magnetic flux A vanishes so that the magnetic field H
is purely transverse and the dynamical system (3.18) and (3.19) is not thereby constrained. Here,
8 H.L. An and C. Rogers
we proceed with the latter case, so that the system (3.18) and (3.19) is additionally constrained
by the requirement Ω2B̄ = const and (3.18)1 yields
B̄N ḠN + B̄SḠS −
δΩ̇
2Ω
= 0.
Finally, for ε2, the relations (2.17), (3.6) and (3.16) combine to show that
ε2 = α2ν
m−2
m−1 Ω2(m−2),
which it subsequently proves convenient to re-write in this form
ε2 = α
(
m− 1
m
)
Ω2(m−2). (3.20)
4 Integrals of motion and parametrisation
Under the constraint (3.20), the nonlinear dynamical system (3.18) is readily shown to admit
the key integrals of motion
B̄2
S + B̄2
N −
1
4
B̄2 = cII, (4.1)
Ḡ2
S + Ḡ2
N − αB̄ = cIII, (4.2)
2(B̄N ḠS − B̄SḠN )− c0B̄ = cIV, (4.3)
2(GR + c0B̄) + 2G(B̄SḠS + B̄N ḠN ) + 4αcIIΩ
−2m− 1
m− 3
− B̄Ω2
(
Ḡ2
S + Ḡ2
N
Ω4
+
G2
4
+G2
R
)
= cV, (4.4)
where cII, cIII, cIV and cV are constants of integration.
The relations (4.1) and (4.2) may be conveniently parametrised according to
B̄S = −
√
cII +
1
4
B̄2 cosφ(t), B̄N = −
√
cII +
1
4
B̄2 sinφ(t),
ḠS = −
√
cIII + B̄ sin θ(t), ḠN = +
√
cIII + B̄ cos θ(t). (4.5)
Substitution of this parametrisation into (3.18)1 yields
˙̄B +
4
Ω2
√
(cII + B̄2/4)(cIII + αB̄) sin(θ − φ) = 0, (4.6)
while conditions (3.18)2,3 reduce to a single relation, namely(
φ̇− f +
2c0
Ω2
)√
cII +
1
4
B̄2 − B̄
Ω2
√
cIII + αB̄ cos(θ − φ) = 0. (4.7)
Similarly, (3.18)4,5 produce another single requirement
(θ̇ − f)
√
cIII + αB̄ +
2α
Ω2
√
cII +
1
4
B̄2 cos(θ − φ) = 0. (4.8)
Moreover, substitution of the representations (4.5) into (4.3) yields
c0B̄ = −cIV + 2
√(
cII + B̄2/4
) (
cIII + αB̄
)
cos(θ − φ), (4.9)
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 9
while elimination of θ − φ in (4.7) and (4.8) respectively shows that
φ̇ = f +
2
Ω2
[
δ(c0δ + cIVΩ2)
δ2 + 4cIIΩ4
− c0
]
(4.10)
and
θ̇ = f − α
Ω2
[
c0δ + cIVΩ2
αδ2 + cIIIΩ2
]
. (4.11)
It remains to consider the nonlinear equation (3.19) for Ω, namely
Ω3Ω̈ +
f2
4
Ω4 + cIII − c20 +
ε2m
m− 1
B̄
Ω2(m−2) = 0
which, by virtue of (3.20), reduces to a generalisation of the classical Steen–Ermakov equa-
tion [3, 22], namely
Ω̈ +
1
4
f2Ω =
c20 − cIII
Ω3
− 2αδ
Ω5
. (4.12)
Further, on use of Theorem 1, it may be readily shown that there is the necessary requirement
(cf. [13])
¨(Ω2B) + f2Ω2B̄ = −2(Q∗ + fM∗)Ω2(m+1) = −2(cV + fcIV).
This holds automatically here with
Ω2B̄ = const = δ = −2(cV + fcIV)/f2, f 6= 0. (4.13)
Elimination of θ − φ between (4.6) and (4.9) now yields
˙̄B2 +
4
Ω4
(c0B̄ + cIV)2 =
4
Ω4
(
4cII + B̄2
)
(cIII + αB̄),
whence on use of (4.13)
δ2Ω̇2 +
(c20 − cIII)δ2
Ω2
+
(
c2IV − 4cIIcIII
)
Ω2 − αδ3
Ω4
+ 2c0cIV − 4αcIIδ = 0. (4.14)
The latter equation is required to be compatible with the 1st integral of (4.12), namely
Ω̇2 +
1
4
f2Ω2 +
(c20 − cIII)
Ω2
− αδ
Ω4
+ k = 0
and these are indeed seen to be consistent subject to the relations
c2IV − 4cIIcIII =
δ2f2
4
,
and
k =
2c0cIV − 4αcIIδ
δ2
.
In summary, a multi-parameter class of exact vortex solutions of the original 2+1-dimensional
magnetogasdynamic system has been generated with the velocity components u1, u2, v1, v2 and
the quantities a, b, c, ρ0 in the density representation given, in turn, by
u1 =
Ω̇
Ω
+
1
Ω3
√
αδ + cIIIΩ2 cos θ(t), v1 =
c0
Ω2
− f
2
− 1
Ω3
√
αδ + cIIIΩ2 sin θ(t),
u2 = − c0
Ω2
+
f
2
− 1
Ω3
√
αδ + cIIIΩ2 sin θ(t), v2 =
Ω̇
Ω
− 1
Ω2
√
αδ + cIIIΩ2 cos θ(t) (4.15)
10 H.L. An and C. Rogers
together with
a =
1
2Ω2(m+1)
[
δ −
√
4cIIΩ4 + δ2 sinφ(t)
]
, b =
1
2Ω2(m+1)
√
4cIIΩ4 + δ2 cosφ(t),
c =
1
Ω2(m+1)
[
δ +
√
4cIIΩ4 + δ2 sinφ(t)
]
, ρ0 =
cI
Ω2(m−1) ,
where the angles φ and θ are obtained by integration of (4.10) and (4.11), respectively while Ω
is given by an elliptic integral resulting from (4.14).
The magnetic flux A is given by
A = νρm−1Ω2(m−1) = ν
[
a(x− q̄)2 + 2b(x− q̄)(y − p̄) + c(y − p̄)2 + ρ0
]
Ω2(m−1),
while the temperature T and entropy distribution S are determined by (2.11) and (2.12), re-
spectively.
5 Hamiltonian Ermakov structure
The nonlinear dynamical system (3.12) may be shown to have remarkable underlying structure
in that it will be seen to reduce to consideration of an integrable Ermakov–Ray–Reid type system
α̈+ ω2(t)α =
1
α2β
F (β/α), β̈ + ω2(t)β =
1
αβ2
G(α/β).
Such systems have their origin in the work of Ermakov [3] and were introduced by Ray and Reid
in [10, 11]. Extension to 2+1-dimensions were presented in [15] and to multi-component systems
in [18]. The main theoretical interest in the system resides in its admittance of a distinctive
integral of motion, namely, the Ray–Reid invariant
I =
1
2
(αβ̇ − βα̇)2 +
∫ β/α
F (z)dz +
∫ α/β
G(w)dw.
Ermakov–Ray–Reid systems arise, in particular, in a variety of contexts in nonlinear optics (see
e.g. [16, 17] and references cited therein).
Here, we proceed with p̄(t) = q̄(t) = 0 in the ansatz (3.1), since the translation terms p̄(t), q̄(t)
are readily re-introduced by use of a Lie group invariance of the magnetogasdynamic system.
The semi-axes of the time-modulated ellipse
a(t)x2 + 2b(t)xy + cy2 + h0(t) = 0, ac− b2 > 0,
are now given by
Φ =
√
2ρ0√
(a− c)2 + 4b2 − (a+ c)
=
√
− cI
2cII
√
−δ −
√
4cIIΩ4 + δ2
and
Ψ =
√
2ρ0
−
√
(a− c)2 + 4b2 − (a+ c)
=
√
− cI
2cII
√
−δ +
√
4cIIΩ4 + δ2,
where it is required that
cI > 0, cII < 0, δ < 0, δ2 + 4cIIΩ
4 > 0.
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 11
It is readily established that the semi-axes Φ, Ψ are governed by a Ermakov–Ray–Reid system,
namely
Φ̈ +
1
4
f2Φ =
1
Φ2Ψ
[
ZZ ′
1 + (Ψ/Φ)2
−
(
Ψ
Φ
)
(Z2 + k/4)
[1 + (Ψ/Φ)2]2
]
,
Ψ̈ +
1
4
f2Ψ =
1
ΦΨ2
[
− ZZ ′
1 + (Ψ/Φ)2
−
(
Φ
Ψ
)
(Z2 + k/4)
[1 + (Φ/Ψ)2]2
]
, (5.1)
where
Z = Z(Φ/Ψ) = ΨΦ̇− Ψ̇Φ =
2cI
Ω
√
−cII
√
(δ2 + 4cIIΩ4)(αδ + cIIIΩ2)− Ω2(c0δ + cIVΩ2)2
δ2 + 4cIIΩ4
and Ω is given in terms of the ratio of the semi-axes via the relation
Ω =
(
− δ
2
cII
)1/4(
Ψ
Φ
+
Φ
Ψ
)−1/2
.
In addition, the Ermakov–Ray–Reid system (5.1) is seen to be Hamiltonian with invariant
H =
1
2
(
Φ̇2 + Ψ̇2
)
− 1
2(Φ2 + Ψ2)
[
Z2 − f2
4
(
Φ2 + Ψ2
)2
+
k
4
]
= −1
4
f2
cIcIv
cII
, cII 6= 0,
and accordingly, is amenable to the general procedure described in detail in [14].
It is remarkable indeed that the semi-axes Φ and Ψ of the time modulated ellipse associated
with the density representation in (3.1), are governed by an integrable Ermakov–Ray–Reid
system, albeit of some complexity. In fact, a Ermakov–Ray–Reid system may also be associated
with the velocity components, at least, in a particular reduction. Attention is here restricted,
as in the work of Dyson [2] on non-conducting gas clouds, to irrotational motions in the absence
of a Coriolis term.
Thus, here we set
L =
(
α̇(t)/α(t) 0
0 β̇(t)/β(t)
)
, E =
(
a(t) 0
0 c(t)
)
in (3.2) corresponding to the subclass of exact solutions in (4.15) with θ = 0, φ = π/2 and
α̇
α
=
Ω̇
Ω
+
1
Ω2
√
cIII +
αδ
Ω2
,
β̇
β
=
Ω̇
Ω
− 1
Ω2
√
cIII +
αδ
Ω2
,
a =
1
2Ω2(m+1)
[
δ −
√
4cIIΩ4 + δ2
]
, c =
1
2Ω2(m+1)
[
δ +
√
4cIIΩ4 + δ2
]
.
The continuity equation, via (3.3), yields
ȧ
a
+
α̇
α
(m+ 1) +
β̇
β
(m− 1) = 0,
ċ
c
+
α̇
α
(m− 1) +
β̇
β
(m+ 1) = 0,
whence
a = cIα
−(m+1)β1−m, c = cIIα
1−mβ−(m+1). (5.2)
Moreover, (3.4) shows that
ρ0 = cIII(αβ)1−m = c∗IIIΩ
2(1−m).
12 H.L. An and C. Rogers
In the above, cI, cII, cIII and c∗III are arbitrary non-zero constants of integration. The momentum
equation gives
α̈+ 2ε2(t)
m
m− 1
aα = 0, β̈ + 2ε2(t)
m
m− 1
cβ = 0 (5.3)
together with
¨̄p = 0, ¨̄q = 0.
Insertion of the expressions (5.2) into (5.3) gives
α̈+ 2ε2(t)
m
m− 1
cI
α2β
(αβ)2−m = 0, β̈ + 2ε2(t)
m
m− 1
cII
αβ2
(αβ)2−m = 0,
whence, in view of the relation (3.20), we again obtain a Ermakov–Ray–Reid system, namely
α̈ =
c∗I
α2β
, β̈ =
c∗II
αβ2
, (5.4)
with the Ray–Reid invariant
I =
1
2
(α̇β − αβ̇)2 + c∗I
β
α
+ c∗II
α
β
,
where
c∗I = −2αcI
(
c∗III
cIII
)m−2
m−1
, c∗II = −2αcII
(
c∗III
cIII
)m−2
m−1
.
It is observed moreover, that the system (5.4) is also Hamiltonian with additional integral of
motion
H =
1
2
(
cIβ̇
2 + cIIα̇
2
)
+
c∗I c
∗
II
αβ
.
6 A Lax pair formulation
It is now shown, following a procedure analogous to that set down in the spinning gas cloud
analysis of [16], that the nonlinear dynamical system (3.3) and (3.10) admits an associated Lax
pair representation. In this connection, it is seen that the nonlinear dynamical system given
by (3.3) together with (3.10) arising from the ansatz (3.1) and (3.2) may be written in the
compact matrix form as
Ė + EL + LTE + (m− 1)E trL = 0,
L̇ + L2 + fPL + 2ε2(t)
m
m− 1
E = 0, m 6= 1,
where L, E are given by (3.2) and
P =
(
0 −1
1 0
)
.
Moreover, the relations (3.4) and (3.11) yield
ρ̇0 + (m− 1)ρ0 trL = 0 and Ṁ + fPM = 0.
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 13
A gauge transformation is now introduced via
L̃ = DLD−1 +
1
2
fP, Ẽ = DED−1,
where
D = exp
(
1
2
Pft
)
to obtain
˙̃E + ẼL̃ + L̃T Ẽ + (m− 1)Ẽ tr L̃ = 0 (6.1)
and
˙̃L + L̃2 +
1
4
f2I + 2ε2(t)
m
m− 1
Ẽ = 0. (6.2)
On use of the relation
PHP = HT − (trH)I
together with the Cayley–Hamilton identity
L̃2 − (tr L̃)L̃ + (det L̃)I = 0
it is seen that (6.2) yields
˙̃L + (tr L̃)L̃− (det L̃)I +
1
4
f2I + 2ε2(t)
m
m− 1
Ẽ = 0. (6.3)
Moreover, on introduction of a new trace-free matrix Q̃ according to
Q̃ = PẼ
the matrix equation (6.1) becomes
˙̃Q + [Q̃, L̃] +mQ̃(tr L̃) = 0. (6.4)
Since trL = tr L̃ = 2Ω̇/Ω, it is natural to introduce the scaling
L̄ = L̃Ω2, Ē = ẼΩ2m, Q̄ = Q̃Ω2m,
whence (6.3) and (6.4) reduce, in turn, to
˙̄L− Ω−2(det L̄)I +
f2
4
Ω2I + 2ε2(t)
m
m− 1
Ω2(1−m)Ē = 0 (6.5)
and
˙̄Q + Ω−2[Q̄, L̄] = 0, (6.6)
where
L̄∗ = L̄− 1
2
(tr L̄)I
14 H.L. An and C. Rogers
denotes the trace-free part of L̄. Moreover, the trace-free part of (6.5) yields
˙̄L∗ + ε2(t)
m
m− 1
Ω2(1−m)[Q̄,P] = 0, (6.7)
while its trace gives
(tr L̄)· − 2Ω−2(det L̄∗)− 1
2
Ω−2(tr L̄)2 +
1
2
f2Ω2 + 2ε2(t)
m
m− 1
Ω2(1−m)(tr Ē) = 0. (6.8)
Insertion of the expression (3.20) for ε2 (with α = 1) into (6.7) yields
˙̄L∗ + Ω−2[Q̄,P] = 0 (6.9)
and on introduction of the new time measure τ according to
dτ = Ω−2dt
the systems (6.6) and (6.9) become in turn
Q̄′ + [Q̄, L̄∗] = 0 and L̄∗′ + [Q̄,P] = 0, (6.10)
where the prime denotes d/dτ . The matrix system (6.10) constitutes the compatibility condition
M′(λ) + [M(λ),L(λ)] = 0
for the Lax pair
Ψ′ = L(λ)Ψ, µΨ =M(λ)Ψ, (6.11)
where
L(λ) = L̄∗ + λP, M(λ) = Q̄ + λL̄∗ + λ2P.
An analogous result has been obtained in the case of non-conducting rotating gas clouds in [19].
As in that work, there is an interesting Steen–Ermakov connection. Thus, on setting
Σ = Ω−1
then the relation (6.8) is readily shown to reduce to a Steen–Ermakov equation, namely
Σ′′ + (det L̄∗ − tr Ē)Σ =
f2
4Σ3
.
Results of [19] related to the Lax pair for a spinning gas cloud system carry over mutatis
mutandis to the Lax pair (6.11) obtained in the present magnetogasdynamic study. Thus, the
linear system (6.11) is gauge equivalent to the standard Lax pair for the stationary reduction
of the integrable cubic nonlinear Schrödinger equation. The connection may be made in the
manner set down in [19].
7 Conclusion
It has been shown via an elliptic vortex ansatz that there is hidden integrable structure of
Ermakov–Ray–Reid type underlying a 2+1-dimensional non-isothermal magnetogasdynamic sys-
tem. The Ermakov variables turn out to have a natural physical interpretation as the semi-axes
of the time-modulated density representation. Moreover, a Lax pair for the original nonlinear
dynamical subsystem has been constructed. The preceeding and previous studies such as that
in [4] suggest that a general investigation of the occurrence of integrable Ermakov–Ray–Reid
structure in 2+1-dimensional hydrodynamic systems would be of interest. It is noted that
Hamiltonian–Ermakov type systems have been additionally investigated in [1, 5].
A 2+1-Dimensional Non-Isothermal Magnetogasdynamic System 15
References
[1] Cerveró J.M., Lejarreta J.D., Ermakov Hamiltonians, Phys. Lett. A 156 (1991), 201–205.
[2] Dyson F.J., Dynamics of a spinning gas cloud, J. Math. Mech. 18 (1969), 91–101.
[3] Ermakov V.P., Second-order differential equations: conditions for complete integrability, Univ. Izv. Kiev 20
(1880), no. 9, 1–25.
[4] Ferapontov E.V., Khusnutdinova K.R., The characterization of two-component (2+1)-dimensional integrable
systems of hydrodynamic type, J. Phys. A: Math. Gen. 37 (2004), 2949–2963, nlin.SI/0310021.
[5] Haas F., Goedert J., On the Hamiltonian structure of Ermakov systems, J. Phys. A: Math. Gen. 29 (1996),
4083–4092, math-ph/0211032.
[6] Neukirch T., Quasi-equilibria: a special class of time-dependent solutions of the two-dimensional magneto-
hydrodynamic equations, Phys. Plasmas 2 (1995), 4389–4399.
[7] Neukirch T., Cheung D.L.G., A class of accelerated solutions of the two-dimensional ideal magnetohydro-
dynamic equations, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 457 (2001), 2547–2566.
[8] Neukirch T., Priest E.R., Generalization of a special class of time-dependent solutions of the two-dimensional
magnetohydrodynamic equations to arbitrary pressure profiles, Phys. Plasmas 7 (2000), 3105–3107.
[9] Ovsiannikov L.V., A new solution of the hydrodynamic equations, Dokl. Akad. Nauk SSSR 111 (1956),
47–49.
[10] Ray J.R., Nonlinear superposition law for generalized Ermakov systems, Phys. Lett. A 78 (1980), 4–6.
[11] Reid J.L., Ray J.R., Ermakov systems, nonlinear superposition, and solutions of nonlinear equations of
motion, J. Math. Phys. 21 (1980), 1583–1587.
[12] Rogers C., A Ermakov–Ray–Reid reduction in 2+1-dimensional magnetogasdynamics, in Group Analysis
of Differential Equations and Integrable Systems, Editors N.M. Ivanova, P.G.L. Leach, R.O. Popovych,
C. Sophocleous, P.A. Damianou, Department of Mathematics and Statistics, University of Cyprus, Nicosia,
2011, 164–177.
[13] Rogers C., Elliptic warm-core theory: the pulsrodon, Phys. Lett. A 138 (1989), 267–273.
[14] Rogers C., An H., Ermakov–Ray–Reid systems in (2+1)-dimensional rotating shallow water theory, Stud.
Appl. Math. 125 (2010), 275–299.
[15] Rogers C., Hoenselaers C., Ray J.R., On (2+1)-dimensional Ermakov systems, J. Phys. A: Math. Gen. 26
(1993), 2625–2633.
[16] Rogers C., Malomed B., An H., Ermakov–Ray–Reid reductions of variational approximations in nonlinear
optics, Stud. Appl. Math., to appear.
[17] Rogers C., Malomed B., Chow K., An H., Ermakov–Ray–Reid systems in nonlinear optics, J. Phys. A:
Math. Theor. 43 (2010), 455214, 15 pages.
[18] Rogers C., Schief W.K., Multi-component Ermakov systems: structure and linearization, J. Math. Anal.
Appl. 198 (1996), 194–220.
[19] Rogers C., Schief W.K., On the integrability of a Hamiltonian reduction of a 2+1-dimensional non-isothermal
rotating gas cloud system, Nonlinearity 24 (2011), 3165–3178.
[20] Rogers C., Schief W.K., The pulsrodon in 2+1-dimensional magneto-gasdynamics: Hamiltonian structure
and integrability, J. Math. Phys. 52 (2011), 083701, 20 pages.
[21] Schäfer G., Gravity-wave astrophysics, in Relativistic Gravity Research with Emphasis on Experiments and
Observations, Lecture Notes in Physics, Vol. 410, Springer-Verlag, Berlin, 1992, 163–183.
[22] Steen A., Om Formen for Integralet af den lineare Differentialligning af anden Orden, Forhandl. af Kjobn-
haven (1874), 1–12.
http://dx.doi.org/10.1016/0375-9601(91)90139-Y
http://dx.doi.org/10.1088/0305-4470/37/8/007
http://arxiv.org/abs/nlin.SI/0310021
http://dx.doi.org/10.1088/0305-4470/29/14/029
http://arxiv.org/abs/math-ph/0211032
http://dx.doi.org/10.1063/1.870995
http://dx.doi.org/10.1098/rspa.2001.0842
http://dx.doi.org/10.1063/1.874167
http://dx.doi.org/10.1016/0375-9601(80)90789-6
http://dx.doi.org/10.1063/1.524625
http://dx.doi.org/10.1016/0375-9601(89)90275-2
http://dx.doi.org/10.1111/j.1467-9590.2010.00488.x
http://dx.doi.org/10.1111/j.1467-9590.2010.00488.x
http://dx.doi.org/10.1088/0305-4470/26/11/012
http://dx.doi.org/10.1111/j.1467-9590.2012.00557.x
http://dx.doi.org/10.1088/1751-8113/43/45/455214
http://dx.doi.org/10.1088/1751-8113/43/45/455214
http://dx.doi.org/10.1006/jmaa.1996.0076
http://dx.doi.org/10.1006/jmaa.1996.0076
http://dx.doi.org/10.1088/0951-7715/24/11/009
http://dx.doi.org/10.1063/1.3622595
http://dx.doi.org/10.1007/3-540-56180-3_8
1 Introduction
2 The magnetogasdynamic system
3 An elliptic vortex ansatz. A dynamical system reduction
4 Integrals of motion and parametrisation
5 Hamiltonian Ermakov structure
6 A Lax pair formulation
7 Conclusion
References
|