Next Article in Journal
IRS-Assisted Hybrid Secret Key Generation
Next Article in Special Issue
On Geometric Interpretations of Euler’s Substitutions
Previous Article in Journal
Audio-Visual Effects of a Collaborative Robot on Worker Efficiency
Previous Article in Special Issue
Remarks on Approximate Solutions to Difference Equations in Various Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solutions of Magnetohydrodynamics Equation through Symmetries

by
Rangasamy Sinuvasan
1,
Amlan K. Halder
2,
Rajeswari Seshadri
3,
Andronikos Paliathanasis
4,5,* and
Peter G. L. Leach
4
1
Department of Mathematics, VIT-AP University, Amaravati 522237, Andhra Pradesh, India
2
School of Technology, Woxsen University, Hyderabad 502345, Telangana, India
3
Department of Mathematics, Pondicherry University, Puducherry 605014, Kalapet, India
4
Institute of Systems Science, Durban University of Technology, Durban 4000, South Africa
5
Departamento de Matemáticas, Universidad Católica del Norte, Avda. Angamos 0610, Casilla, Antofagasta 1280, Chile
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(10), 1908; https://doi.org/10.3390/sym15101908
Submission received: 9 August 2023 / Revised: 22 September 2023 / Accepted: 10 October 2023 / Published: 12 October 2023

Abstract

:
The magnetohydrodynamics (1 + 1) dimension equation, with a force and force-free term, is analysed with respect to its point symmetries. Interestingly, it reduces to an Abel’s Equation of the second kind and, under certain conditions, to equations specified in Gambier’s family. The symmetry analysis for the force-free term leads to Euler’s Equation and to a system of reduced second-order odes for which singularity analysis is performed to determine their integrability.

1. Introduction

In this paper, we focus on the magnetohydrodynamics (MHD) ( 1 + 1 ) -dimension equation, as proposed in reference [1]. Fleischer et al. presents the effect of intermittency in the turbulence of plasmas, by specifying a replica of the Burger’s equation for the magnetohydrodynamics (MHD) equation. This paper discusses in detail the physical significance of each parameter and the impact on the properties of the ( 1 + 1 ) dimensional MHD under a nonzero external force. It is to be noted here that the model proposed by Fleischer and Diamond is a plainer version of the presence of magnetic field to ( 1 + 1 ) -dimensional MHD equation. As suggested in reference [1], it is to be considered as an enhanced model in comparison to the one cited in reference [2]. Later on, Basu et al. [3] discussed in detail the physical significance of the equation and its corresponding similarity to certain well-known equations. Moreover, a comprehensive list of historical details of the equation are provided in references [3,4,5,6,7,8,9]. The equation in ( 1 + 1 ) dimension is of the following form:
u t + B 0 v x + u u x + v v x = ν u x x + f ( t , x ) , v t + B 0 u x + ( u v ) x = μ v x x + g ( t , x ) ,
where f ( t , x ) and g ( t , x ) are external forces, ν , μ are fluid kinematic viscosity and magnetic diffusivity, respectively. u and v are velocity field and magnetic field, respectively. B 0 is a mean magnetic field (independent of x and t).
In the literature, a comprehensive work has been carried out on the various dimensions of MHD equations. In reference [10], Fuchs discussed the symmetries and the corresponding reductions at a stretch for the three-dimensional MHD equation. The references of the paper also cited some interesting work done by Nucci [11] on the group analysis of ( 1 + 3 ) dimensional MHD equations along with the equations with respect to its magnetic component. Further, we suggest that readers consult references [12,13] for an in-depth study. We also recommend references [14,15] which discuss the impact of symmetries on the stabilitiy of the MHD equations. Also, in references [16,17] the constitution of the MHD equation is specified elaborately.
Our main objective in this work is to compute certain new similarity reductions and solutions for Equation (1) using Lie group analysis. The Lie point symmetry analysis does provide new results for Equation (1). The equation reduces to Abel’s equation of the second kind and, under certain conditions, leads to Gambier’s equation [18], the classification number of which is presented against each equation. Also, we performed a symmetry analysis of the force-free term. The subsequent reductions leads to an Euler’s equation and to a system of second-order odes, with respect to travelling-wave, for which we conduct a singularity analysis.
The application of the symmetry analysis can also be extended to problems consisting of boundary conditions. The main criteria is that the equation under study remains invariant under the transformations and all the properties of the domain are preserved. It is to be noted that, for certain equations, a composition of solutions can be constructed from the derived invariant solutions and that particular composition could possibly satisfy the boundary conditions. Here, it is important to emphasize that the invaraint solutions obtained by the symmetries can be considered to be more general, which remains valid for any bounded or unbounded domain and, hence, it describes the system properly and efficiently. Therefore, our main objective of the computation of solutions through the similarity transformations is more generic, which could possibly satisfy the boundary conditions.
This paper is organised as follows. In Section 2, the preliminaries of the computation of the point symmetries are mentioned. In Section 3, the point symmetries and subsequent computations are done for the external forces. In Section 4, the reduction to the corresponding odes and its symmetries, which are used for further reductions, are mentioned. In Section 5 and Section 6, the particular cases of a second-order equation obtained in Section 4 is discussed in detail. Finally, in Section 7, a symmetry analysis of the Equation (1) without an external force is presented and a singularity analysis test for certain reduced equations is elaborated. The conclusion, acknowledgements, and proper references follow henceforth.

2. Preliminaries

Any differential equation is understood to possess a symmetry, provided the equation possesses its basic properties and characteristics, even after undergoing a transformation with respect to its coordinates. Mathematically, it can be explained as follows:
Consider a PDE of the form F ( x ^ , u , u ^ ) = 0 , where x ^ represents the array of independent variables, u represents the dependent variable, and u ^ represents the array of derivatives of the dependent variable. Our objective is to compute the point symmetries of the PDE. First of all, certain transformations need to be defined as follows:
x ^ ̲ = x ^ + ϵ ξ ( x ^ , u ) + O ( ϵ 2 ) u ̲ = u + ϵ η ( x ^ , u ) + O ( ϵ 2 ) .
The transformations with respect to the each derivative term in u ^ can be derived easily using (2) with the aid of the Chain rule. Now, the PDE, to be invariant, need to satisfy the condition
F ( x ^ , u , u ^ ) = F ( x ^ ̲ , u ̲ , u ^ ̲ ) .
The set of transformations (2) forms a one-parameter Lie group, the generator of which is represented as
Γ = ξ ( x ^ , u ) + η ( x ^ , u ) .
The operator (3) is the generator of the infinitesimal transformations ξ ( x ^ , u ) and η ( x ^ , u ) . This opertor is being used to define the characteristics using the method of the associated Lagarange’s system and it is to be considered a point symmetry of the PDE under study. Further, these generators are mainly used to reduce the order of the PDE, to compute conservation laws, and, also, it has been used recently to derive nonlocal symmetries. The reader can follow references [19,20,21,22,23,24,25] to understand the basics of point symmetries and the procedure of computation.

3. Lie Symmetries of the Magnetohydrodynamics Equation

Corresponding to the differential Equation (1), we write a general element of the symmetry vector field of Equation (1), in the form
V = ξ 1 t + ξ 2 x + η 1 u + η 2 v .
where ξ 1 , ξ 2 , η 1 , and η 2 are functions of t , x , u , and v.
Requiring that the second prolongation P r ( 1 ) V of Equation (4) should annihilate (1) on the solution set, we obtain the determining equations for the coefficients ξ 1 , ξ 2 , η 1 , and η 2 . The solution of the determining equations are given by
ξ 1 = a 1 ( t )
ξ 2 = a 2 ( t ) + x a 1 ˙ 2
η 1 = a 2 ¨ + x a 1 ¨ 2 u a 1 ˙ 2
η 2 = 1 2 ( v + B 0 ) a 1 ˙
Provided f ( t , x ) and g ( t , x ) satisfies the following linear partial differential equations,
2 a 1 f t + ( 2 a 2 + x a ˙ 1 ) f x + 3 f a ˙ 1 2 a ¨ 2 x a 1 = 0
2 a 1 g t + ( 2 a 2 + x a ˙ 1 ) g x + 3 g a ˙ 1 = 0
where a 1 and a 2 are arbitrary functions of t. Then, the Lie symmetry algebra of Equation (1) is
Γ = a 1 ( t ) t + a 2 ( t ) + x a 1 ˙ 2 x + a 2 ˙ + x a 1 ¨ 2 u a 1 ˙ 2 u + 1 2 ( v + B 0 ) a 1 ˙ v
Since (11) is a first-order partial differential equation, the equation of the characteristic curve may be expressed invariantly using the Lagrange–Charpit method. Equation (11) can be rewritten in the Lagrange form.
d t a 1 = d x 1 2 a ˙ 1 x + a 2 = d f 3 2 a ˙ 1 f + 1 2 x a 1 + a ¨ 2 .
The characteristic systems are
d x d t a ˙ 1 2 a 1 x a 2 a 1 = 0
d f d t + 3 a ˙ 1 2 a 1 f a 1 2 a 1 x a ¨ 2 a 1 = 0 .
The first equation gives a characteristic z = ( x α ) ρ = c o n s t a n t , where a 1 = ρ 2 and α = ρ a 2 ρ 3 d t .
Then, the Equation (14) becomes
d f d t + 3 a ˙ 1 2 a 1 f 1 2 ( a 1 1 2 z + α ) a 1 a 1 a ¨ 2 a 1 = 0
and gives the characteristic
H 1 ( z ) = f ρ 3 ρ 2 ρ ¨ x ( α ¨ ρ α ρ ¨ ) ρ 2 = c o n s t a n t .
Then, the solution of Equation (9) is
f ( x , t ) = ρ ¨ x α ρ + α ¨ + 1 ρ 3 H 1 x α ρ ,
where H is an arbitrary function of its argument.
Similarly,
g ( x , t ) = 1 ρ 3 H 2 x α ρ ,
is the solution of Equation (10), where H 2 is an arbitrary function of its argument.
It is important to note, here, that the analysis of the linear counterpart of Equation (1) leads to certain important observations. The linear counterpart of Equation (1) is obtained by considering the terms u u x , v v x , and ( u v ) x equal to zero. The subsequent search for symmetries lists certain trivial and nontrivial symmetries and certain equations similar to (9) and (10). The explicit symmetries are as follows:
Γ 1 = t Γ 2 = x Γ 3 = u u + v v .
Now, along with these symmetries, certain arbitrary functions representing infinite-dimensional symmetries are also obtained, which can be used productively to obtain the reduced ordinary differential equations of the subsequent linear counterpart of Equation (1). It is also important to note here that the symmetries mentioned as Γ 1 , Γ 2 , and Γ 3 are the trivial symmetries inherent to any linear equations and, hence, forms the main difference between the linear and nonlinear equations. The reduced equations, using these symmetries, may lead to certain trivial equations which are easily solvable and, hence, the plot of these solutions can be easily understood.

4. Reduction to an Ordinary Differential Equation

The characteristics corresponding to the generator (4) are found by solving the equation
d t a 1 = d x 1 2 a 1 ˙ x + a 2 = d u 1 2 a 1 ¨ x 1 2 a 1 ˙ u + a 2 ˙ = d v 1 2 a ˙ 1 ( v + B 0 ) .
The first characteristic for the Equation (18) is z = ( x α ) ρ , which is the new similarity variable, where α and ρ are defined as a 1 = ρ 2 and α = ρ a 2 ρ 3 d t . The second characteristic can be obtained by solving the equation
d u d t + a ˙ 1 2 a 1 u a ¨ 1 2 a 1 ( ρ z + α ) a ˙ 2 a 1 = 0 .
The solution is
G ( z ) = ρ u ρ ρ ˙ z ρ α ˙ .
From the above solution, u can be written as
u = ρ ˙ z + α ˙ + G ( z ) ρ .
The third characteristic can be obtained by solving the equation
d v d t + a ˙ 1 2 a 1 v B 0 a ¨ 1 2 a 1 = 0 .
The solution is
F ( z ) = ρ v + ρ B 0 .
From the above solution, v can be written as
v = F ( z ) ρ B 0 .
In the above, G ( z ) and F ( z ) are similarity functions, when one substitutes the values of u , v , f , and g into Equation (1), we have
H 1 ( z ) + G ( z ) G ( z ) + F ( z ) F ( z ) ν G ( z ) = 0 H 2 ( z ) + F ( z ) G ( z ) + G ( z ) F ( z ) μ F ( z ) = 0
Letting H 1 ( z ) = H 1 ( z ) and H 2 ( z ) = H 2 ( z ) and integrating (25) with respect to z, we have a system of first-order ordinary differential equations:
G ( z ) 2 2 + F ( z ) 2 2 + H 1 ( z ) ν G ( z ) = 0
F ( z ) G ( z ) + H 2 ( z ) + μ F ( z ) = 0
We can find the value of G ( z ) from Equation (27) and, substituting into Equation (26), we get the second-order ordinary differential equation:
F = μ 2 ν + 1 F 2 F + 1 ν + 1 μ H 2 F F + 1 2 ν μ F 3 + 1 ν μ H 1 F 1 μ H 2 + H 2 2 2 μ ν F
Now, apply the transformation F ( z ) = 1 w ( z ) and the above equation becomes
w = 1 μ 2 ν w 2 w + 1 ν + 1 μ H 2 w w H 2 2 2 ν μ w 3 + 1 μ H 2 2 H 1 μ ν w 1 2 μ ν w
The general case does not fit into the known equation. However, we perform another symmetry reduction and identify it to an Abel’s Equation of the second kind.
F = 1 + μ 2 ν F 2 F + 1 μ + 1 ν H 2 F F + F 3 2 μ ν + H 1 F μ ν H 2 μ + H 2 2 2 μ ν F
The infinitesimal generator is
Γ = p 2 p 1 z μ μ ν z + F p 1 μ μ ν F ,
provided H 1 ( z ) = k H 2 ( z ) and H 2 ( z ) = h 1 ( p 1 z μ + p 2 ( ν μ ) ) 2 .
The canonical coordinates corresponding to the vector field (31) are
r = F ( p 1 z μ + p 2 ( ν μ ) ) , s = ( μ ν ) log ( p 1 z μ + p 2 ( ν μ ) ) p 1 μ .
Then, it is easy to see that
d s d r = ν μ ( p 1 z μ + p 2 ( ν μ ) ) ( p 1 μ F ( p 1 z μ + p 2 ( ν μ ) ) F ) .
If we choose v to be d s d r 1
d v d r = p 1 2 μ 2 F + ( p 1 z μ + p 2 ( ν μ ) ) ( 3 p 1 μ F + ( p 1 z μ + p 2 ( ν μ ) ) F ) ( ν μ ) ( p 1 μ F + ( p 1 z μ + p 2 ( ν μ ) ) F )
r 2 4 μ 3 ( μ + 6 ν ) + r 2 + 2 r ( 2 μ ( μ ν ) + 5 r ) + 4 μ ( μ ν ) v ( r )             μ + ν + 2 μ 2 r ( μ + 5 ν ) + μ ν r ( μ ν ) d v d r + 4 μ 3 ( μ ν ) 2 ( μ + 2 ν ) v ( r ) 2 + 1 = 0

5. Particular Cases of Equation (28)

5.1. Case 1 (Gambier.B 19)

Let μ = ν = i 2 2 , H 2 = 0 , H 1 = 2 ν 2 and apply the transformation z x , F y , we obtain
y = y 2 y + 4 y 3 + 2 y ,
with the analytic solution
d y 4 y 2 + 4 ln y + y 0 y = x x 0 .
The numerical evolution of the latter equation is presented in Figure 1.

5.2. Case 2 (Gambier.B 29)

Let μ = ν = i 3 , H 1 = H 2 = 0 , and apply the transformation z x , F y , we obtain
y = y 2 y + 3 2 y 3 ,
with the analytic solution
y x = 8 y 0 e 2 2 y 1 x y 0 2 e 2 y 1 x 48 and y x = 8 y 0 e 2 2 y 1 x y 0 2 e 2 y 1 x 48 .
The numerical evolution of the latter equation is presented in Figure 2.

5.3. Case 3 (Gambier.B 30)

Let μ = ν = i 3 , H 1 = 2 β μ ν , H 2 = γ ν , and apply the transformation z x , F y , we obtain
y = y 2 2 y + 3 2 y 3 + 2 β y γ 2 2 y
where the solution is given in terms of the elliptic integral
2 d y 4 y 2 y 0 + 6 y 4 + 16 y 2 β ln a + 2 γ 2 = x x 0 .
This can be considered a variant of Gambier.B 30 for α = 0 ,
y = y 2 2 y + 3 2 y 3 + 4 α y 2 + 2 β y γ 2 2 y .
The numerical evolution of the latter equation is presented in Figure 3.

5.4. Case 4 (Second Painlevé Transcendent)

Let μ = 2 ν = i 2 , H 2 = 0 , H 1 = 2 ν 2 z and apply the transformation F w , z x . We obtain the second Painlevé transcendent:
w 3 + x w = 0 .
The numerical evolution of the latter equation is presented in Figure 4.

5.5. Case 5 (Kummer–Schwarz Equation)

If μ = ν and H 2 = 0 , and apply the transformation z x , F y , we obtain the second-order Kummer–Schwarz equation
F = 3 2 F 2 F + 1 2 μ 2 F 3 + H 1 μ 2 F .
The numerical evolution of the latter equation is presented in Figure 5.

5.6. Case 6 (Duffing Equation)

If H 1 is constant H 2 = 0 and μ = 2 ν , and apply the transformation z x , F y , we obtain the undamped and unforced Duffing equation:
F + 1 4 ν 2 F 3 + 1 2 ν 2 H 1 F = 0
The numerical evolution of the latter equation is presented in Figure 6.

6. General Cases of Equation (29)

Equation (29) can be considered as a particular case of the following two Gambier’s equations:

6.1. Case 7 Gambier.B 28

y = y 2 2 y y y + q y + y 3 2 2 q y 2 + 3 ( q + q 2 2 ) y 72 r 2 y

6.2. Case 8 Gambier.B 27

y = ( 1 1 n ) y 2 2 y + f n ( q , r ) y y + ϕ n ( q , r ) y n 2 n y y + n f n 2 ( n + 2 ) 2 y 3 + n ( f n + f n ϕ n ( n + 2 ) y 2 + ψ n ( q , r ) y + ϕ n + 1 n y

7. The Case f = 0, g = 0

The equation has the Lie point symmetries
Γ 1 = t Γ 2 = x Γ 3 = t x + u
In the subsequent subsections, we study the reductions with respect to Γ 1 , Γ 2 , Γ 3 , and the travelling-wave. In Section 7.2, the singularity analysis of the reduced ode is also mentioned explicitly.

7.1. Case 6a

The reduction, with respect to Γ 1 , leads to the reduced ode
F 1 ( x ) F 1 ( x ) + B 0 G 1 ( x ) + G 1 ( x ) G 1 ( x ) = ν F 1 ( x ) , B 0 F 1 ( x ) + G 1 ( x ) F 1 ( x ) + F 1 ( x ) G 1 ( x ) = μ G 1 ( x ) ,
where u ( t , x ) = F 1 ( x ) and v ( t , x ) = G 1 ( x ) . The point symmetries of Equation (35) are
Γ 1 a = x , Γ 2 a = x x F 1 F 1 B 0 + G 1 G 1 .
The reductions with respect to Γ 1 a and Γ 2 a does not yield any positive result. The singularity analysis of the Equation (35), too, does not provide any fruitful result.
Next, we look for the reduction with respect to Γ 2 , which leads to the simpler first-order ode F 2 ( t ) = 0 , and G 2 ( t ) = 0 , where u ( t , x ) = F 2 ( t ) and v ( t , x ) = G 2 ( t ) .
Further, the reduction with respect to Γ 3 , leads to an Euler’s equation, which is
F 3 ( t ) + t F 3 ( t ) = 0 , B 0 + G 3 ( t ) + t G 3 ( t ) = 0 ,
where u ( t , x ) = x t + F 3 ( t ) and v ( t , x ) = G 3 ( t ) .

7.2. Case 6b

In this subsection, we perform the reductions with respect to the travelling-wave. The reduction leads to a system of two second-order odes, which are
B 0 F 1 ( z ) + F 1 ( z ) F 1 ( z ) c G 1 ( z ) + G 1 ( z ) G 1 ( z ) = ν G 1 ( z ) , c F 1 ( z ) + G 1 ( z ) F 1 ( z ) + B 0 G 1 ( z ) + F 1 ( z ) G 1 ( z ) = μ F 1 ( z ) ,
where the similarity variables are
x c t = z , u ( t , x ) = G 1 ( z ) , v ( t , x ) = F 1 ( z ) .
The Lie point symmetries of equation (38) are
Γ 1 b = z , Γ 2 b = z z B 0 + F 1 F 1 + c G 1 G 1 .
The reduction, with respect to Γ 1 b and Γ 2 b , does not provide any fruitful result, similar to Equation (35). Therefore, we dwell on the singularity analysis to comprehend its integrability.
Singularity Analysis of Equation (38).
We substitute F 1 α z p and G 1 β z q , in Equation (38). The substitution leads to
B 0 p z 1 + p α + p z 1 + 2 p α 2 c q z 1 + q β + q z 1 + 2 q β 2 + q z 2 + q β ν q 2 z 2 + q β ν = 0 , c p z 1 + p α + B 0 q z 1 + q β + p z 1 + p + q α β + q z 1 + p + q α β + p z 2 + p α μ p 2 z 2 + p α μ = 0 ,
One of the possible values for the exponents p and q is 1 . After substituting the values of p and q in Equation (39), we use the dominant terms to compute the leading-order coefficients α and β . The terms are
α 2 z 3 β 2 z 3 2 β ν z 3 , 2 α β z 3 2 α μ z 3
Solving Equation (40) for α and β leads to the following possible values
α 0 , β 0 , α 2 μ ν μ 2 , β μ , α 2 μ ν μ 2 , β μ , α 0 , β 2 ν .
Next, to compute the resonance, we substitute F 1 α z p + m 1 z p + s , G 1 β z q + m 2 z q + s in Equation (38). The substitution leads to
0 = B 0 m 1 z 2 + s + c m 2 z 2 + s + B 0 m 1 s z 2 + s c m 2 s z 2 + s m 1 2 z 3 + 2 s m 2 2 z 3 + 2 s + m 1 2 s z 3 + 2 s + m 2 2 s z 3 + 2 s B 0 α z 2 2 m 1 z 3 + s α + m 1 s z 3 + s α α 2 z 3 + c β z 2 2 m 2 z 3 + s β + m 2 s z 3 + s β β 2 z 3 2 m 2 z 3 + s ν + 3 m 2 s z 3 + s ν m 2 s 2 z 3 + s ν 2 β ν z 3 , 0 = c m 1 z 2 + s B 0 m 2 z 2 + s c m 1 s z 2 + s + B 0 m 2 s z 2 + s 2 m 1 m 2 z 3 + 2 s + 2 m 1 m 2 s z 3 + 2 s + c α z 2 2 m 2 z 3 + s α + m 2 s z 3 + s α B 0 β z 2 2 m 1 z 3 + s β + m 1 s z 3 + s β 2 α β z 3 2 m 1 z 3 + s μ + 3 m 1 s z 3 + s μ m 1 s 2 z 3 + s μ 2 α μ z 3 .
We select the linear terms with respect to m 1 and m 2 , corresponding to various powers of z . For z 2 + s , we have
B 0 + B 0 s m 1 + c c s m 2 = 0 , c c s m 1 + B 0 s B 0 m 2 = 0 .
For system (43) to possess a non-trivial solution, the necessary requirement is
B 0 ( s 1 ) c ( 1 s ) c ( 1 s ) B 0 ( s 1 ) = 0 .
The possible values of s = 1 ( 2 ) . Similarly, for z 3 + s , we have
2 m 1 α + m 1 s α 2 m 2 β + m 2 s β 2 m 2 ν + 3 m 2 s ν m 2 s 2 ν = 0 , 2 m 2 α + m 2 s α 2 m 1 β + m 1 s β 2 m 1 μ + 3 m 1 s μ m 1 s 2 μ = 0
Similarly to Equation (43), for the existence of a non-trivial solution, the system
( s 2 ) α ( β + ν s ν ) α ( β μ s μ ) = 0 ,
leads to the following possible values of s.
s 2 , s β μ + β ν + 2 μ ν β 2 μ 2 + 4 α 2 μ ν 2 β 2 μ ν + β 2 ν 2 2 μ ν , s β μ + β ν + 2 μ ν + β 2 μ 2 + 4 α 2 μ ν 2 β 2 μ ν + β 2 ν 2 2 μ ν .
Now, for α , β = 0 , 2 ν , the possible values of s are
s ν μ 4 μ 2 ν 2 8 μ ν 3 + 4 ν 4 2 μ ν , s ν μ + 4 μ 2 ν 2 8 μ ν 3 + 4 ν 4 2 μ ν .
For α 2 μ ν μ 2 , β μ , the possible values of s are
s 1 2 μ 2 ν μ 4 2 μ 3 ν + μ 2 ν 2 + 4 μ ν ( μ 2 + 2 μ ν ) 2 μ ν , s 1 2 μ 2 ν + μ 4 2 μ 3 ν + μ 2 ν 2 + 4 μ ν ( μ 2 + 2 μ ν ) 2 μ ν .
Similar results are obtained for α 2 μ ν μ 2 , β μ . The generic value of the resonance s being 1 is not obtained for both the computations, and, hence, as cited in references [26,27,28,29], the integrability of Equation (38) cannot be determined.
Still, a possibility exists when ν , fluid kinematic viscosity, and μ magnetic diffusivity are equal. This leads to s = 1 . Therefore, under a special circumstance of ν = μ , we have s = 1 , and from the analysis of Equation (44) s also takes the value 2. The readers can easily verify the consistency test and, hence, we conclude that the Laurent expansion leads to a right Painlevé series.

8. Conclusions

In this work, reductions concerned with that of the arbitrary functions are deduced for the magnetohydrodynamics(MHD) ( 1 + 1 ) -dimension equation. The reduction leads to a certain well-known family of equations which are connected to Gambier’s, Abel’s, and Kummer–Schwarz equations. The reductions of the force-free equation are generally trivial but the traveling-wave reduction leads to equations which are devoid of point symmetries. Hence, the computation of a series solution using singularity analysis leads to a right Painlevé series, provided the kinematic viscosity and magnetic diffusivity parameters are equal. The physical significance of such a situation is yet to be ascertained and forms the basis of our future work.

Author Contributions

Methodology, R.S. (Rangasamy Sinuvasan), A.K.H. and P.G.L.L.; Investigation, A.K.H. and R.S. (R. Seshadri); Writing—original draft, A.K.H. and A.P.; Writing #x2014;review & editing, A.P. and P.G.L.L.; Supervision, A.P. and P.G.L.L. All authors have read and agreed to the published version of the manuscript.

Funding

PGLL acknowledges the support of the National Research Foundation of South Africa, the University of KwaZulu-Natal, and the Durban University of Technology.

Data Availability Statement

The authors declare that no datasets were used in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fleischer, J.; Diamond, P.H. Compressible Alfven turbulence in one dimension. Phys. Rev. E 1998, 58, R2709. [Google Scholar] [CrossRef]
  2. Politano, H.; Pouquet, A. Model of intermittency in magnetohydrodynamic turbulence. Phys. Rev. E 1995, 52, 636. [Google Scholar] [CrossRef] [PubMed]
  3. Basu, A.; Bhattacharjee, J.K.; Ramaswamy, S. Mean magnetic field and noise cross-correlation in magnetohydrodynamic turbulence: Results from a one-dimensional model. Eur. J. B Condens. Matter Complex Syst. 1999, 9, 725–730. [Google Scholar] [CrossRef]
  4. Basu, A.; Sain, A.; Dhar, S.K.; Pandit, R. Multiscaling in models of magnetohydrodynamic turbulence. Phys. Rev. Lett. 1998, 81, 2687. [Google Scholar] [CrossRef]
  5. Bhattacharjee, J.K. Randomly stirred fluids, mode coupling theories and the turbulent Prandtl number. J. Phys. A Math. Gen. 1988, 21, L551. [Google Scholar] [CrossRef]
  6. Camargo, S.J.; Tasso, H. Renormalization group in magnetohydrodynamic turbulence. Phys. Fluids B Plasma Phys. 1992, 4, 1199–1212. [Google Scholar] [CrossRef]
  7. Forster, D.; Nelson, D.R.; Stephen, M.J. Large-distance and long-time properties of a randomly stirred fluid. Phys. Rev. A 1977, 16, 732. [Google Scholar] [CrossRef]
  8. Lahiri, R.; Ramaswamy, S. Are steadily moving crystals unstable? Phys. Rev. Lett. 1997, 79, 1150. [Google Scholar] [CrossRef]
  9. Yakhot, V.; Orszag, S.A. Renormalization group analysis of turbulence. I. Basic theory. J. Sci. Comput. 1986, 1, 3–51. [Google Scholar] [CrossRef]
  10. Fuchs, J.C. Symmetry groups and similarity solutions of MHD equations. J. Math. Phys. 1991, 32, 1703–1708. [Google Scholar] [CrossRef]
  11. Nucci, M.C. Group analysis for MHD equations. Atti Sem. Mat. Fis. Univ. Modena 1984, 33, 21–34. [Google Scholar]
  12. Gross, J. Invariante L Sungen der Eindimensionalen Nichtstationa ren Realen MHD-Gleichungen; Technische Universita t Carolo-Wilhelmina zu Braunschweig: Braunschweig, Germany, 1983. [Google Scholar]
  13. Grundland, A.M.; Lalague, L. Lie subgroups of symmetry groups of fluid dynamics and magnetohydro-dynamics equations. Can. Phys. 1995, 73, 463–477. [Google Scholar] [CrossRef]
  14. Goedbloed, J.; Lifschitz, A. Stationary symmetric magnetohydrodynamic flows. Phys. Plasmas 1997, 4, 3544–3564. [Google Scholar] [CrossRef]
  15. Knobloch, E. Symmetry and instability in rotating hydrodynamic and magnetohydrodynamic flows. Phys. Fluids 1996, 8, 1446–1454. [Google Scholar] [CrossRef]
  16. Kraichnan, R.H. Inertial-range spectrum of hydromagnetic turbulence. Phys. Fluids 1965, 8, 1385–1387. [Google Scholar] [CrossRef]
  17. She, Z.S.; Leveque, E. Universal scaling laws in fully developed turbulence. Phys. Rev. Lett. 1994, 72, 336. [Google Scholar] [CrossRef]
  18. Gambier, B. Sur les equations differentielles du second ordre et du premier degre dont l’integrale generale est a points critiques fixes. Acta Math. 1910, 33, 1–55. [Google Scholar] [CrossRef]
  19. Meleshko, S.V. Methods for Constructing Exact Solutions of Partial Differential Equations; Springer Science: New York, NY, USA, 2005. [Google Scholar]
  20. Dimas, S.; Tsoubelis, D. A new Mathematica-based program for solving overdetermined systems of PDEs. In Proceedings of the 8th International Mathematica Symposium, Avignon, France, 19–23 June 2006. [Google Scholar]
  21. Polyanin, A.D.; Zaitsev, V.F. Handbook of Nonlinear Equations of Mathematical Physics. Exact Solutions; Fizmatlit: Moscow, Russia, 2002. (In Russian) [Google Scholar]
  22. Dimas, S.; Tsoubelis, D. SYM: A new symmetry-finding package for Mathematica. In Proceedings of the 10th International Conference in Modern Group Analysis; University of Cyprus: Nicosia, Cyprus, 2004; pp. 64–70. [Google Scholar]
  23. Ince, E.L. Ordinary Differential Equations; Longmans, Green & Co.: London, UK, 1927. [Google Scholar]
  24. Dimas, S. Partial differential equations, algebraic computing and nonlinear systems. Ph.D. Thesis, University of Patras, Patras, Greece, 2008. [Google Scholar]
  25. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
  26. Andriopoulos, K.; Leach, P.G.L. An interpretation of the presence of both positive and negative nongeneric resonances in the singularity analysis. Phys. Lett. A 2006, 359, 199–203. [Google Scholar] [CrossRef]
  27. Andriopoulos, K.; Leach, P.G.L. Singularity analysis for autonomous and nonautonomous differential equations. Appl. Discret. Math. 2011, 5, 230–239. [Google Scholar] [CrossRef]
  28. Andriopoulos, K.; Leach, P.G.L. The occurrence of a triple-1 resonance in the standard singularity. Nuovo C. Della Societa Ital. Fis. B Gen. Phys. Relativ. Astron. Math. Phys. Methods 2009, 124, 1–11. [Google Scholar]
  29. Andriopoulos, K.; Leach, P.G.L. Symmetry and singularity properties of second-order ordinary differential equation of Lie’s Type III. J. Math. Anal. Appl. 2007, 328, 860–875. [Google Scholar] [CrossRef]
Figure 1. Qualitative evolution of the Gambier.B 19 equation.
Figure 1. Qualitative evolution of the Gambier.B 19 equation.
Symmetry 15 01908 g001
Figure 2. Qualitative evolution of the Gambier.B 29 equation.
Figure 2. Qualitative evolution of the Gambier.B 29 equation.
Symmetry 15 01908 g002
Figure 3. Qualitative evolution of the Gambier.B 30 equation.
Figure 3. Qualitative evolution of the Gambier.B 30 equation.
Symmetry 15 01908 g003
Figure 4. Qualitative evolution of the second Painlevé transcendent equation.
Figure 4. Qualitative evolution of the second Painlevé transcendent equation.
Symmetry 15 01908 g004
Figure 5. Qualitative evolution of the Kummer–Schwarz equation.
Figure 5. Qualitative evolution of the Kummer–Schwarz equation.
Symmetry 15 01908 g005
Figure 6. Qualitative evolution of the Duffing equation.
Figure 6. Qualitative evolution of the Duffing equation.
Symmetry 15 01908 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sinuvasan, R.; Halder, A.K.; Seshadri, R.; Paliathanasis, A.; Leach, P.G.L. Solutions of Magnetohydrodynamics Equation through Symmetries. Symmetry 2023, 15, 1908. https://doi.org/10.3390/sym15101908

AMA Style

Sinuvasan R, Halder AK, Seshadri R, Paliathanasis A, Leach PGL. Solutions of Magnetohydrodynamics Equation through Symmetries. Symmetry. 2023; 15(10):1908. https://doi.org/10.3390/sym15101908

Chicago/Turabian Style

Sinuvasan, Rangasamy, Amlan K. Halder, Rajeswari Seshadri, Andronikos Paliathanasis, and Peter G. L. Leach. 2023. "Solutions of Magnetohydrodynamics Equation through Symmetries" Symmetry 15, no. 10: 1908. https://doi.org/10.3390/sym15101908

APA Style

Sinuvasan, R., Halder, A. K., Seshadri, R., Paliathanasis, A., & Leach, P. G. L. (2023). Solutions of Magnetohydrodynamics Equation through Symmetries. Symmetry, 15(10), 1908. https://doi.org/10.3390/sym15101908

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop