Skip to main content

Lie group analysis and similarity solutions for hydro-magnetic Maxwell fluid through a porous medium

Abstract

The equations of two dimensional incompressible fluid flow for hydro-magnetic Maxwell fluid through a porous medium have been studied. Lie group analysis has been employed and the group invariant solutions are obtained. Solutions corresponding to translational and rotational symmetries are obtained. A boundary value problem for the translational symmetry is investigated and the results are also sketched graphically. The effects of physical parameters have been noticed.

MSC 2011: 53C11; 76S05.

1 Introduction

Non-Newtonian fluid behavior, which is characterized by a nonlinear viscosity dependence on the strain, can be observed in many complex fluids, for example, polymers, dense colloidal dispersions, surfactant solutions, micellar solutions chemical, and petroleum industries [1]. In addition to shear-thinning and shear-thickening behavior, a dynamic or even chaotic response can be found in some fluids subjected to a steady shear flow. Because of the difficulty to suggest a single model which exhibits all properties of non-Newtonian fluids, they cannot be described as simply as Newtonian fluids. Due to this fact many models of constitutive equations have been proposed and most of them are empirical or semi empirical [2]. Amongst these the differential type fluid model gained considerable attention of many researchers. The flows of non-Newtonian fluids are not only important because of their technological significance but also in the interesting mathematical features presented by the equations governing the flow. However on the other hand there are much controversies on these models as well. Such fluids are also inadequate to describe the relaxation phenomena. For a complete and detailed discussion of the relevant issues for differential type fluids, we refer the readers to Dunn and Rajagopal [3] and Aksel [4].

The non-Newtonian fluids are mainly classified into three types namely differential, rate and integral. The simplest subclass of the rate type fluids is the Maxwell model [5]. This fluid model can very well describe the relaxation time effects. Specifically the Maxwell fluid model has been used for the viscoelastic flows where the dimensionless relaxation time is small. However in some more concentrated polymeric fluids the Maxwell model is also useful for large dimensionless relaxation time. Some recent investigations dealing with the flows of Maxwell fluids are given in the references [69].

Modified Darcy's law for a Maxwell fluid including the Hall current has been used for the modeling. In fact, the Hall effect is important when the Hall parameter, which is the ratio between the electron-cyclotron frequency and the electronatom-collision frequency, is high. This happens, when the magnetic field is high or when the collision frequency is low. In most cases, the Hall term has been ignored in applying Ohm's law as it has no marked effects for small and moderate values of the magnetic field. However, the current trend in the application of magnetohydrodynamics is towards a strong magnetic field, so that the influence of electromagnetic force is noticeable. Under these conditions, the Hall current is important and it has marked effects on the magnitude and direction of the current density and consequently on the magnetic-force term. Therefore, it is of interest to study the influence of the Hall current on the flow.

In the Earth there are a large number of problems that can be described by the interaction of a low viscosity fluid (water, oil, gas, magma) in a permeable (possibly deformable) matrix. Darcy's Law is the classic, empirically derived equation for the flux of a low viscosity fluid in a permeable matrix. This equation assumes that flow in the pores or cracks of the medium is essentially laminar and provides the average flux through a representative area that is larger than the pore scale and smaller than the scale of significant permeability variation (if such a scale exists). Various approaches have been used to justify this rule from first principles (see e.g., Dagan [10]) but it generally seems to work.

In this article, we apply the so-called symmetry methods for a particular problem of fluid mechanics. The main advantage of such methods is that they can successfully be applied to nonlinear differential equations [1113]. The similarity solutions are quite popular because they result in the reduction of the independent variables of the problem. The symmetry transformations method transform the given family of equations of n independent variables, say, to another family of equations of n - 1 independent variables, which can further be solved [14, 15]. The fundamental concepts of this approach can be found in [1619]. In our case, the problem under investigation is (2 + 1)-nonlinear partial differential equations (PDEs). Hence, any similarity solution will transform the system of (n + 1)-nonlinear PDEs into a system of (n)-nonlinear PDEs and any similarity solution will transform the system of (2)-nonlinear PDEs into a system of ordinary differential equations (ODEs).

Many authors used Lie group analysis method to obtain the exact solutions for some problems in fluid mechanics. Yurusoy and Pakdemirli [20] investigated the boundary layer equations of a non-Newtonian fluid model in which the shear stress is an arbitrary function of the velocity gradient. Yurusoy et al. [21] have obtained the solution for the creeping flow of the second grade fluid. Also the two-dimensional equations of motions for the slowly flowing and heat transfer in second grade fluid in cartesian coordinates neglecting the inertial terms are considered by Yürüsoy [22]. Shahzad et al. [23] found the analytical solution of a micropolar fluid by using the Lie group analysis. Recently, Mekheimer et al. studied the Lie group analysis and similarity solutions for a couple stress fluid with heat transfer [24], Lie point symmetries and similarity solutions for an electrically conducting Jeffrey fluid [25] and Lie point symmetries and similarity solution for a micro-polar fluid through a porous medium [26].

From discussion above, we attend to find the analytical (similarity) solutions for the flow problem of an incompressible hydro-magnetic Maxwell fluid through a porous medium using Lie group analysis. The problem is presented as follows, in Section 2, the equations governing two-dimensional motion of an incompressible, MHD Maxwell fluid are introduced. In Section 3, the basic idea of the Lie group analysis method are given and used to find the isovector field of our equations. The similarity solutions corresponding to translational and rotational symmetry are obtained in Sections 3.1 and 3.2. Also a boundary value problem for the similarity solutions corresponding to translational symmetry are obtained in Section 4. The graphs for a boundary value problem (magma flow) are plotted and discussed in Section 5. Finally a concluding remarks are pointed in Section 6.

2 Equations of motion

The continuity and momentum equations governing the two-dimensional motion of an incompressible hydro-magnetic Maxwell fluid through a porous medium can be written as:

E 1 = u ˜ x ˜ + v ˜ y ˜ = 0 , E 2 = ( 1 + λ ˜ t ) [ u ˜ t ˜ + u ˜ u ˜ x ˜ + v ˜ u ˜ y ˜ ] + ( 1 + λ ˜ t ˜ ) 1 ρ p ˜ x ˜ μ ρ 2 u ˜ + ( 1 + λ ˜ t ˜ ) [ σ B 0 2 θ ρ ( u ˜ m v ˜ ) ] + μ φ ρ k ˜ u ˜ = 0 , E 3 = ( 1 + λ ˜ t ˜ ) [ v ˜ t ˜ + u ˜ v ˜ x ˜ + v ˜ v ˜ y ˜ ] + ( 1 + λ ˜ t ˜ ) 1 ρ p ˜ y ˜ μ ρ 2 v ˜ + ( 1 + λ ˜ t ˜ ) [ B 0 2 θ ρ ( v ˜ + m u ˜ ) ] + μ φ ρ k ˜ v ˜ = 0 ,
(1)

where 2 = 2 x ˜ 2 + 2 y ˜ 2 , θ = 1 1 + m 2 , m = B 0 e n e are the fluid velocities in the x ̃ , directions, p ̃ is the pressure, and t ̃ is the time. Here λ ̃ , ρ, μ, φ, k, e, n e , σ, B0, and m are the relaxation time, density, coefficient of viscosity, porosity of the porous medium, permeability, electric charge, the number density of electrons, electrical conductivity of the fluid, magnetic field and Hall parameter respectively.

Using the following dimensionless parameters

u = u ˜ U , v = v ˜ U , x = x ˜ L , y = y ˜ L , t = U L t ˜ , p = p ˜ ρ U 2 , λ = U L λ ˜ , k = k ˜ φ L 2 ,
(2)

the system (1) becomes

E 1 = u x + v y = 0 , E 2 = 1 + λ t u t + u u x + v u y + 1 + λ t p x - 1 R 2 u + 1 + λ t M θ ( u - m v ) + 1 R k u = 0 , E 3 = 1 + λ t v t + u v x + v v y + 1 + λ t p y - 1 R 2 v + 1 + λ t M θ ( v + m u ) + 1 R k v = 0 ,
(3)

where R= ρ L U μ is the Reynolds number, M= σ B 0 2 L ρ U is the Hartmann number and L, U are the dimensionless length and velocity, respectively.

3 Lie group analysis and isovector fields

In order to obtain the analytical solution, we apply the Lie group analysis theory to system (3). For this we write

x i * = x i + ε ξ i ( x j , u α ) + o ( ε 2 ) , u α * = u α + ε η α ( x j , u β ) + o ( ε 2 ) , i , j = 1 , 2 , 3 , α , β = 1 , 2 , 3 ,
(4)

as the infinitesimal Lie point transformations. We have assumed that system (3) is invariant under the transformations given in Eq. (4). The corresponding infinitesimal generator of Lie groups (symmetries) is given by

X = ξ i x i + η α u α ,
(5)

with summation convention over the repeated index and x1 = x, x2 = y, x3 = t, u1 = u, u2 = v, u3 = p. The coefficients ξ1, ξ2, ξ3, η1, η2, and η3 are the functions of all independent and dependent variables. There coefficients are the components of the infinitesimals symmetries corresponding to x, y, t, u, v, and p, respectively to be determined from the invariance conditions:

Pr ( 2 ) X ( E a ) E α = 0 = 0 , a = 1 , 2 , 3 ,
(6)

where E a = 0, i = 1, 2, 3 represent the system of Eq. (3) and Pr(2) is the second prolongation of the isovector field X. Since the system (3) is of order two, then our prolongation will be in the form

P r ( 1 ) X = X + η α i u α , i , P r (2) X = P r ( 1 ) X + η α i j u α , i j .
(7)

where

η α i = D i η α - ξ j u α , j + ξ j u α , j i , η α i j = D i j η α - ξ k u α , k + ξ k u α , k i j .
(8)

and the operator D i 1 i 2 . . . i s is called the total derivative (Hash operator) and has the following form:

D i = x i + u α , i u α + u α , i j u α , j + u α , i j k u α , j k ,
(9)

where D ij = D i (D j ) = D j (D i ) = D ji and u α , i = u α x i .

Expanding the system of Eq. (6) with the aid of Mathematica programm, along with the original system of Eq. (3) to eliminate u x , p xt , p yt and setting the coefficients involving u y , u yy , v x , v y , v xx , v xy , v yy and various products to zero give rise the essential set of over-determined equations. Solving these set of determining equations we obtain the required components of isovector field as follows:

ξ 1 = a 2 - a 1 y , ξ 2 = a 3 + a 1 x , ξ 3 = a 4 , η 1 = - α 1 v , η 2 = a 1 u , η 3 = δ ( t ) ,
(10)

where a i , i = 1,..., 5 are arbitrary constants, δ(t) is arbitrary function of the variable t only.

3.1 Translational symmetry

In this case we take a1 = 0. The characteristic equations corresponding to the translational symmetry are:

d x a 2 = d y a 3 = d t a 4 = d u 0 = d v 0 = d p δ ( t ) .
(11)

By solving the ODEs (11), we can obtain the similarity variables and similarity functions as follows:

ϕ = x - m 1 t , ψ = y - m 2 t ,
(12)
u ( x , y , t ) = u ^ ( ϕ , ψ ) , v ( x , y , t ) = v ^ ( ϕ , ψ ) , p ( x , y , t ) = p ^ ( ϕ , ψ ) + δ 1 ( t ) ,
(13)

where m 1 = a 2 a 4 , m 2 = a 3 a 4 are arbitrary constants and δ1(t) = ∫ δ(t) dt is an arbitrary function. Substituting the transformations (12), (13) in the Eq. (3) lead to the following system of PDEs:

{ E 1 = u ^ ϕ + v ^ ψ = 0 , E 2 = u ^ [ 1 R k + M θ + u ^ ϕ λ ( m 2 u ^ ϕ ψ + m 1 u ^ ϕ ϕ ) ] 1 R ( u ^ ϕ ϕ + u ^ ψ ψ ) + λ m 2 2 u ^ ϕ ϕ + p ^ ϕ v ^ [ m M θ u ^ ψ + λ ( m 2 u ^ ψ ψ + m 1 u ^ ϕ ψ ) ] m 2 [ u ^ ψ ( 1 + M λ θ + λ v ^ ψ + λ u ^ ϕ ) + λ ( m M θ v ^ ψ 2 m 1 u ^ ϕ ψ + p ^ ϕ ψ ) ] m 1 [ u ^ ϕ ( 1 + M λ θ + λ u ^ ϕ ) + λ ( u ^ ψ v ^ ϕ m M θ v ^ ϕ m 1 u ^ ϕ ϕ + p ^ ϕ ϕ ) ] = 0 , E 3 = v ^ [ 1 R k + M θ + v ^ ψ λ ( m 2 v ^ ψ ψ + m 1 v ^ ϕ ψ ) ] 1 R ( v ^ ϕ ϕ + v ^ ψ ψ ) λ m 2 2 v ^ ψ ψ p ^ ψ + u ^ [ m M θ v ^ ϕ λ ( m 2 v ^ ϕ ψ + m 1 v ^ ϕ ϕ ) ] + m 2 [ v ^ ψ ( 1 + M λ θ + λ v ^ ψ ) + λ ( u ^ ψ ( m M θ + v ^ ϕ ) 2 m 1 v ^ ϕ ψ + p ^ ψ ψ ) ] + m 1 [ v ^ ϕ ( 1 + M λ θ + λ v ^ ψ ) + λ ( u ^ ϕ v ^ ϕ + m M θ u ^ ϕ + m 1 v ^ ϕ ϕ + p ^ ϕ ψ ) ] = 0 ,
(14)

To transform Eq. (14) to an (ODEs), we use the Lie group analysis again and obtain the infinitesimal generator corresponding to system of equation (14) in the following form:

ξ 11 = b 1 , ξ 21 = b 2 , η 11 = 0 , η 21 = 0 , η 31 = β ( ϕ , ψ ) ,
(15)

where b i , i = 1, 2 are arbitrary constants and β(ϕ, ψ) an arbitrary function that satisfy two conditions:

β ϕ - λ m 2 β ϕ ψ + m 1 β ϕ ϕ = 0 , β ψ - λ m 2 β ψ ψ + m 1 β ϕ ψ = 0 ,
(16)

By solving the characteristic equations,

d ϕ b 1 = d ψ b 2 = d u ^ 0 = d v ^ 0 = d p ^ β ( ϕ , ψ ) ,
(17)

if we take β(ϕ, ψ) = 0 then we can obtain the similarity variable and similarity functions as the following:

χ = ψ - m 0 ϕ , u ^ ( ϕ , ψ ) = f ( χ ) , v ^ ( ϕ , ψ ) = g ( χ ) , p ^ ( ϕ , ψ ) = h ( χ ) ,
(18)

where m 0 = b 2 b 1 an arbitrary constant. Substituting the transformations (18) in Eq. (14) lead to the following system of ODEs:

g - m 0 f = 0 f 1 R k + M θ - m 0 f - λ m 0 m 0 m 1 - m 2 f - 1 R 1 + m 0 2 f - g m M θ - f - λ m 0 m 1 - m 2 f - m 2 1 + M λ θ + λ g f - λ m M θ g + m 2 f + m 0 m 1 1 + M λ θ + λ g f - λ m M θ g + 2 m 2 f + λ m 0 - m 0 m 1 - m 2 ( f ) 2 + m 0 m 1 2 f + λ m 0 m 0 m 1 - m 2 h - m 0 h = 0 , g 1 R k + M θ + g + λ m 0 m 1 - m 2 g - 1 R ( 1 + m 0 2 ) g - f m M θ - m 0 g - m 0 λ m 0 m 1 - m 2 g - m 2 1 + M λ θ + λ g g + λ m M θ f - m 2 g + m 0 m 1 1 + M λ θ + λ g g + λ m M θ f - 2 m 2 g + λ m 0 - m 0 m 1 - m 2 f g + m 0 m 1 2 g + λ m 0 m 1 - m 2 h + h = 0 .
(19)

Integrating the first equation in (19) yields

g = m 0 f + c 1 ,
(20)

where c1 is an arbitrary constant. Eliminating h(ξ) from the second and third equations in (19) along with Eq. (20) we get the following equation:

A f + B f + C f + D = 0 .
(21)

where

A = - 1 R 1 + m 0 2 + λ m 0 2 m 1 2 + m 2 2 - 2 m 0 m 1 m 2 + c 1 m 0 m 1 - m 2 , B = c 1 + 1 + M λ θ m 0 m 1 - m 2 , C = 1 R k + M θ , D = c 1 1 + m 0 2 m 0 1 R k + M θ - m M θ .
(22)

By solving equation (21) we get

f ( χ ) = c 2 exp [ α 1 χ ] + c 3 exp [ α 2 χ ] - D ̃ ,
(23)

where D ̃ = D C and c2, c3 are arbitrary constants and α1, α2 are roots of the following equation:

A γ 2 + B γ + C = 0 ,
(24)

From equations (20) and (23) the expression of the function g(χ) becomes

g ( χ ) = m 0 c 2 exp [ α 1 χ ] + c 3 exp [ α 2 χ ] - D ̃ + c 1 ,
(25)

and from the second and third equations in (19) we get

h ( χ ) = D ̃ m M θ - c 1 ( 1 + m 0 2 ) 1 R k + M θ + m M θ m 0 χ - m M θ c 2 α 1 exp [ α 1 χ ] + c 3 α 2 exp [ α 2 χ ] - λ c 4 m 0 m 1 - m 2 exp - χ λ ( m 0 m 1 - m 2 ) + c 5 ,
(26)

where c4 and c5 are arbitrary constants.

In the form of the original variables, our exact solutions can be written as follows:

u ( x , y , t ) = c 2 exp α 1 y - m 0 x + m 0 m 1 - m 2 t + c 3 exp α 2 y - m 0 x + m 0 m 1 - m 2 t - D ̃ , v ( x , y , t ) = m 0 c 2 exp α 1 y - m 0 x + m 0 m 1 - m 2 t + c 3 exp α 2 y - m 0 x + m 0 m 1 - m 2 t - D ̃ + c 1 , p ( x , y , t ) = D ̃ m M θ - c 1 ( 1 + m 0 2 ) 1 R k + M θ + m M θ m 0 y - m 0 x + m 0 m 1 - m 2 - m M θ c 2 α 1 exp α 1 y - m 0 x + m 0 m 1 - m 2 t + c 3 α 2 exp α 2 y - m 0 x + m 0 m 1 - m 2 t λ - c 4 m 0 m 1 - m 2 exp - y - m 0 x + m 0 m 1 - m 2 t λ ( m 0 m 1 - m 2 ) + c 5 + δ 1 ( t ) .
(27)

3.2 Rotational symmetry

In this section, the parameter a1 is taken to be an arbitrary non-zero constant. The characteristic equations corresponding to the rotational symmetry are:

d x - a 1 y + a 2 = d y a 1 x + a 3 = d t a 4 = d u - a 1 v = d v a 1 u = d p δ ( t ) .
(28)

Integrating equations (28) using the Lie group analysis method, we get the rotationally invariant solutions for our problem in the following form:

ψ = ( x + β 1 ) 2 + ( y + β 2 ) 2 , ϕ = tan - 1 y + β 2 x + β 1 , u = G 1 ( ψ , t ) cos [ ϕ ] + G 2 ( ψ , t ) sin [ ϕ ] , v = G 1 ( ψ , t ) sin [ ϕ ] - G 2 ( ψ , t ) cos [ ϕ ] , p = G 3 ( ψ , t ) + δ 1 ( t ) ,
(29)

where β 1 = a 3 a 1 , β 2 = - a 2 a 1 , and G3 are functions of ψ and t.

Substituting the new (similarity) variables (ϕ, ψ) and functions (G1, G2, G3) into the original system (3) yields the following set of equations:

G 1 ψ + G 1 ψ = 0 , G 1 k + ψ 2 1 + R k M θ G 1 ψ + λ 2 G 1 ψ t + k ψ R G 2 m M θ ψ - G 2 - 2 λ G 2 t - G 1 ψ - k ψ 2 2 G 1 ψ 2 + R k ψ 2 1 + M λ θ + λ G 1 ψ G 1 t + G 3 ψ + λ m M θ G 2 t + 2 G 1 t 2 + 2 G 3 ψ t = 0 , G 2 k + ψ 2 1 + R k M θ + R k λ ψ G 1 t + k ψ R G 1 G 2 + λ G 2 t - G 2 ψ - ψ 2 G 2 ψ 2 + R k ψ 2 ( 1 + M λ θ ) G 2 t + λ 2 G 2 t 2 + G 1 t G 2 ψ + G 1 2 G 2 ψ t - m M θ G 1 - m M λ θ G 1 t + G 1 G 2 ψ = 0
(30)

To transform Eq. (30) to (ODEs), we use the Lie group analysis again and obtain the infinitesimal generator corresponding to system of equations (30) in the following form:

ξ 12 = 0 , ξ 22 = d 1 , η 11 = 0 , η 21 = 0 , η 32 = η ( ψ , t , G 1 ) ,
(31)

where d1 is an arbitrary constant and η(ψ,t,G1) an arbitrary function that satisfy the following condition:

ψ η ψ - G 1 η G 1 = 0 .
(32)

If we take η(ψ,t,G1) = 0 is a simple solution of Eq. (32) then the characteristic equations are:

d ψ 0 = d t d 1 = d G 1 0 = d G 2 0 = d G 3 0 .
(33)

The similarity variables and resulting functions are

ζ = ψ , G 1 = F 1 ( ζ ) , G 2 = F 2 ( ζ ) , G 3 = F 3 ( ζ ) .
(34)

One now substitutes the similarity variable and the functions into the equations (30) and obtains

F 1 + F 1 ζ = 0 , F 1 k + ζ 2 1 + R k M θ + R k F 1 + k ζ R F 2 m M θ ζ - F 2 - F 1 + R ζ F 3 - ζ F 1 = 0 , F 2 k + ζ 2 1 + R k M θ + R k ζ F 1 F 2 + ζ F 2 - m M θ - k ζ F 2 + ζ F 2 = 0 .
(35)

Integrating the first equation in the system (35) we get:

F 1 = g 1 ζ ,
(36)

where g1 is an arbitrary constant. From the second and the third equations in (35) along with (36) we get:

F 3 - F 2 2 ζ + m M θ F 2 - g 1 2 ζ 3 + 1 + R k M θ g 1 R k ζ = 0 ,
(37)
F 2 k + ζ 2 1 + R k M θ + R k g 1 F 2 + ζ F 2 - m M θ - k ζ F 2 + ζ F 2 = 0 ,
(38)

by integrating equation (38) we obtain

F 2 = 1 k 2 ( 1 + R g 1 2 ) ζ R g 1 2 ( i ζ 1 k + R M θ ) ( 1 + R g 1 2 ) [ 2 2 + R g 1 k g 3 csc [ 1 2 π R g 1 ] 0 F 1 R e g u l a r i z e d ( R g 1 2 , ζ 2 ( 1 + R k M θ ) 4 k ) ζ 2 ( i ζ 1 k + R M θ ) ( 1 + R k M θ ) [ g 2 + g 3 cot [ 1 2 π R g 1 ] ] 0 F 1 R e g u l a r i z e d ( 2 + R g 1 2 , ζ 2 ( 1 + R k M θ ) 4 k ) ] R k m M θ g 1 ζ ( 1 + R k M θ ) [ 1 + 0 F 1 ( R g 1 2 , ζ 2 ( 1 + R k M θ ) 4 k ) ] ,
(39)

and

F 3 = F 2 2 ζ - m M θ F 2 + g 1 2 ζ 3 - 1 + R k M θ g 1 R k ζ d ζ ,
(40)

where g2 and g3 are arbitrary constants, and the function p F q and p F q Regularized are defined as:

p F q a ; b ; z = HypergeometricPFQ { a 1 , . . . , a p } ; { b 1 , . . . , b q } ; z = k = 0 ( a 1 ) k . . . ( a p ) k ( b 1 ) k . . . ( b q ) k z k k ! ,
(41)
p F q Regularized = p F q Γ ( b 1 ) ... Γ ( b q ) .
(42)

Then the solution of our problem in the original variables using symbolic computations is:

u ( x , y , t ) = g 1 ψ cos [ ϕ ] + { 1 k 2 ( 1 + R g 1 2 ) ψ R g 1 2 ( i ψ 1 k + R M θ ) ( 1 + R g 1 2 ) [ 2 2 + R g 1 k g 3 csc [ 1 2 π R g 1 ] 0 F 1 R e g u l a r i z e d ( R g 1 2 , ψ 2 ( 1 + R k M θ ) 4 k ) ψ 2 ( i ψ 1 k + R M θ ) ( 1 + R k M θ ) [ g 2 + g 3 cot [ 1 2 π R g 1 ] ] 0 F 1 R e g u l a r i z e d ( 2 + R g 1 2 , ψ 2 ( 1 + R k M θ ) 4 k ) ] R k m M θ g 1 ψ ( 1 + R k M θ ) [ 1 + 0 F 1 ( R g 1 2 , ψ 2 ( 1 + R k M θ ) 4 k ) ] } sin [ ϕ ] v ( x , y , t ) = g 1 ψ sin [ ϕ ] { 1 k 2 ( 1 + R g 1 2 ) ψ R g 1 2 ( i ψ 1 k + R M θ ) ( 1 + R g 1 2 ) [ 2 2 + R g 1 k g 3 csc [ 1 2 π R g 1 ] 0 F 1 R e g u l a r i z e d ( R g 1 2 , ψ 2 ( 1 + R k M θ ) 4 k ) ψ 2 ( i ψ 1 k + R M θ ) ( 1 + R k M θ ) [ g 2 + g 3 cot [ 1 2 π R g 1 ] ] 0 F 1 R e g u l a r i z e d ( 2 + R g 1 2 , ψ 2 ( 1 + R k M θ ) 4 k ) ] R k m M θ g 1 ψ ( 1 + R k M θ ) [ 1 + 0 F 1 ( R g 1 2 , ψ 2 ( 1 + R k M θ ) 4 k ) ] } cos [ ϕ ] , p ( x , y , t ) = ( F 2 2 ψ m M θ F 2 + g 1 2 ψ 3 ( 1 + R k M θ ) g 1 R k ψ ) d ψ + δ 1 ( t ) ,
(43)

where ψ and ϕ are the same in (29).

4 Solutions for hydro-magnetic Maxwell fluid through a porous medium: (magmatic fluid) problem

One of the important applications in geology is the magmatic fluid. Consider a magmatic fluid as an incompressible hydro-magnetic Maxwell fluid through a porous medium and a plate over it. The plate occupies the position y = 0, where the positive y goes deep into the fluid beneath the plate. The relevant boundary conditions are of the form:

u ( x , 0 , 0 ) = U 0 , u ( x , , t ) = 0 , u y ( x , 0 , 0 ) = 0 , v ( x , 0 , 0 ) = - V 0 , p ( x , , 0 ) = P 0 , p ( x , 0 , 0 ) = P a ,
(44)

where U0 is the velocity of the plate, V0 is the magmatic fluid velocity penetrating into the plate, P0 is the pressure deep in the magmatic fluid and P a is the atmosphere pressure. The expressions (27) for the translational symmetry case solution after using conditions (44) give

u ( x , y , t ) = U 0 α 2 - α 1 α 2 exp α 1 ( y + W t ) - α 1 exp α 2 ( y + W t ) , v ( x , y , t ) = V 0 α 2 - α 1 α 2 exp α 1 ( y + W t ) - α 1 exp α 2 ( y + W t ) , p ( x , y , t ) = m M θ U 0 α 1 α 2 ( α 2 - α 1 ) - α 1 2 exp α 2 ( y + W t ) + α 2 2 exp α 1 ( y + W t ) - p 0 - p a + m M θ U 0 α 1 α 2 ( α 1 + α 2 ) exp - ( y + W t ) C λ + P 0 ,
(45)

where W = m0m1 - m2, α1 and α2 are the negative roots of Eq. (24).

5 Discussion of the magmatic fluid problem

This section deals with the graphics on the magmatic fluid. So, the interpretation of the relaxation time λ, Reynolds number R, Hartmann number M, Hall parameter m, the time parameter t, and the permeability parameter k have been studied on the pressure p, and the x and y components of the velocity distributions u and v.

Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, and 14 describe the variations of the velocity components u and v with the time t at y = 0 for different values of permeability parameter k, Hall parameter m, the relaxation time λ, Reynolds number R, and Hartmann number M. For all of these figures at y = 0 we note that as the time t increases the variation of each component of the velocity decreases and vanishes for large values of t. This is expected, where for small values of t and at the magma plate (y = 0), the variation of the velocity components is obvious. Also, the gap between the curves for small values of t at the magma plate increases than those as t increases.

Figure 1
figure 1

Variation of the dimensionless velocity distribution along the x -axis with t for different values of permeability parameter k ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; y = 0; m = 0.5; M = 0.5; R = 0.5; λ = 50).

Figure 2
figure 2

Variation of the dimensionless velocity distribution along the y -axis with t for different values of permeability parameter k ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; y = 0; m = 0.5; M = 0.5; R = 0.5; λ = 50).

Figure 3
figure 3

Variation of the dimensionless velocity distribution along the x -axis with t for different values of Hall parameter m ( U 0 = V 0 = 2; y = 0; m 0 = m 1 = m 2 = 2; k = 0.05; M = 0.5; R = 0.5; λ = 50).

Figure 4
figure 4

Variation of the dimensionless velocity distribution along the y -axis with t for different values of Hall parameter m ( U 0 = V 0 = 2; y = 0; m 0 = m 1 = m 2 = 2; k = 0.05; M = 0.5; R = 0.5; λ = 50).

Figure 5
figure 5

Variation of the dimensionless velocity distribution along the x -axis with t for different values of relaxation time λ ( U 0 = V 0 = 2; y = 0; m 0 = m 1 = m 2 = 2; k = 0.05; M = 0.5; R = 0.5; m = 0.5).

Figure 6
figure 6

Variation of the dimensionless velocity distribution along the y -axis with t for different values of relaxation time λ( U 0 = V 0 = 2; y = 0; m 0 = m 1 = m 2 = 2; k = 0.05; M = 0.5; R = 0.5; m = 0.5).

Figure 7
figure 7

Variation of the dimensionless velocity distribution along the x -axis with t for different values of Reynolds number R ( U 0 = V 0 = 2; y = 0; m = 0.5; m 0 = m 1 = m 2 = 2; k = 0.05; M = 0.5; λ = 50).

Figure 8
figure 8

Variation of the dimensionless velocity distribution along the y -axis with t for different values of Reynolds number R ( U 0 = V 0 = 2; y = 0; m = 0.5; m 0 = m 1 = m 2 = 2; k = 0.05; M = 0.5; λ = 50).

Figure 9
figure 9

Variation of the dimensionless velocity distribution along the x -axis with t for different values of Hartmann number M ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; y = 0; k = 0.05; R = 0.5; λ = 50; m = 0.5).

Figure 10
figure 10

Variation of the dimensionless velocity distribution along the y -axis with t for different values of Hartmann number M ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; y = 0; k = 0.05; R = 0.5; λ = 50; m = 0.5).

Figure 11
figure 11

Variation of the dimensionless velocity distribution along the x -axis with t for different values y ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; M = 0.5; k = 0.05; R = 0.5; λ = 50; m = 0.5).

Figure 12
figure 12

Variation of the dimensionless velocity distribution along the y -axis with t for different values y ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; M = 2; k = 0.05; R = 50; λ = 50; m = 0.5).

Figure 13
figure 13

Variation of the dimensionless velocity distribution along the x -axis with t for different values of the velocity U 0 ( V 0 = 2; y = 0; m 0 = m 1 = m 2 = 2; k = 0.05; R = 0.5; λ = 50; m = 0.5; M = 0.5).

Figure 14
figure 14

Variation of the dimensionless velocity distribution along the y -axis with t for different values of magmatic velocity V 0 ( U 0 = 2; y = 0; m 0 = m 1 = m 2 = 2; k = 0.05; R = 0.5; λ = 50; m = 0.5; M = 0.5).

Figures 1 and 2 show that as the permeability parameter k increases the horizontal velocity component u increases, while the vertical velocity component v decreases. Figures 3 and 4 illustrate the variation of the velocity components u and v with the Hall parameter m, which indicate that for small values of t (or at initial values of t) the curves are the same with no obvious different which for t > 2, the gap between the curves appears. Also, we can see that curves with small values of m (m = 0, 0.5) are vanishing rabidly than those for (m = 1,1.5) i.e., as the Hall parameter m increases the disturbance of the velocity components increase. (decreasing the number of density electrons or the electronic charges).

Figures 5 and 6 illustrate the variations of u and v with t for different values of the relaxation time λ, which show that for small values λ the disturbance in u and v will vanish rapidly than those as λ increases. Also, the figures show that the disturbance in u and v for a Newtonian fluid less than those for a Non-Newtonian fluid in the case of magma flow.

Figures 7 and 8 show that the variation with the Reynolds number R. As R increases the velocity components u and v increase. Figures 9 and 10 show that as the Hartmann number M increases the velocity components u and v decrease, i.e., the fluid moves as a block and takes a constant value for large values of M. Figures 11 and 12 illustrate the variation of u and v with t for different values of the y axis, which show that the velocity components take the initial values of the magma plate at y = 0 and the velocity components decreases as y increases. Figures 13 and 14 describe the variations of u and v with t for different values of U0 and V0 (velocities of the magma plate), the figures show that the gab between the curves decreases with time and finally vanishes and for certain values of t the velocity components u and v equal to zero. Also, the magnitudes of u and v increase with increasing U0 and V0.

Figures 15, 16, 17, 18, 19, and 20 illustrate the variation of the pressure p with y for different values of k, m, λ, R, M, and t. We can see that the pressure decreases as the permeability parameter k increases and take a constant value for large values of k. However an inverse effective behavior for p with the Hartmann number M is shown in Figure 19. Figures 16 and 17 describe the variation of p with the Hall parameter m and the relaxation time λ, which shows that as y increases the pressure increases and the gabs between the curves are more obvious near to the magma plate. Also, the pressure takes the same values of the pressure deep in the magmatic fluid P0 for the large values of y (as we move deep into the fluid) and the same effect is shown with λ. The pressure decreases as the Reynolds number R increases as shown in Figure 18. Figure 20 shows the variation of the pressure with y for different values of t. We can see that the pressure increases as t increases for small values of y.

Figure 15
figure 15

Variation of the dimensionless pressure distribution with y for different values of permeability parameter k ( U 0 = V 0 = 2; p 0 = 5; pa = 1; m 0 = m 1 = m 2 = 2; t = 0; m = 0.5; M = 2; R = 50; λ = 50).

Figure 16
figure 16

Variation of the dimensionless pressure distribution with y for different values of Hall parameter m ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; t = 0; k = 0.5; M = 2; R = 50; λ = 50; p 0 = 5;p a = 1).

Figure 17
figure 17

Variation of the dimensionless pressure distribution with y for different values of relaxation time λ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; t = 0; k = 0.5; M = 2; R = 50; m = 0.5; p 0 = 5; p a = 1).

Figure 18
figure 18

Variation of the dimensionless pressure distribution with y for different values of Reynolds number R ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; t = 0; k = 0.5; M = 2; λ = 50; m = 0.5; p 0 = 5; p a = 1).

Figure 19
figure 19

Variation of the dimensionless pressure distribution with y for different values of Hartmann number M ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; t = 0; k = 0.5; R = 50; λ = 50; m = 0.5; p 0 = 5; p a = 1).

Figure 20
figure 20

Variation of the dimensionless pressure distribution with y for different values of time t ( U 0 = V 0 = 2; m 0 = m 1 = m 2 = 2; M = 2; k = 0.5; R = 50; λ = 50; m = 0.5; p 0 = 5; p a = 1; y = 0).

Other cases of symmetry will be considered for other boundary value problems else where for other applications.

6 Concluding remarks

The significant features of Lie group analysis for hydro-magnetic Maxwell fluid through a porous medium have been presented. Similarities solutions are obtained and applied to an important phenomena in geology, which is the magmatic fluid. The main points have been summarized as follows:

  • As the Hall parameter m increases the disturbances of the velocity components are increase.

  • The disturbances in the fluid velocity components for a Newtonian fluid are less than those for a non-Newtonian fluid (magmatic fluid)

  • The magmatic fluid moves as a block for large values of the Hartmann number M.

  • The pressure near to the magma plate is higher for a magneto-magma flow than that for a magma flow without a magnetic field. Also, this pressure for a porous medium is less than that for a medium with high permeability.

  • The pressure increase near to the magma plate and take the constant value of the pressure deep in the magmatic fluid for large values of y.

References

  1. Zhaosheng Y, Jianzhong L: Numerical research on the coherent structure in the viscoelastic second-order mixing layers. Appl Math Mech 1998, 19: 671–677.

    Article  Google Scholar 

  2. Shifang H: Constitutive Equation and Computational Analytical Theory of Non-Newtonian Fluids. Science Press, Beijing; 2000.

    Google Scholar 

  3. Zhaosheng JE, Rajagopal KR: Fluids of differential type: critical review and thermodynamic analysis. Int J Eng Sci 1995, 33: 689–729. 10.1016/0020-7225(94)00078-X

    Article  Google Scholar 

  4. Aksel N: A brief note from the editor on the second-order fluid. Acta. Mech 2002, 157: 235–236. 10.1007/BF01182167

    Article  MATH  Google Scholar 

  5. Maxwell JC: On the dynamical theory of gases. Philos Trans R Soc Lond A 1866, 157: 26–78.

    Google Scholar 

  6. Fetecau C, Fetecau C: Decay of a potential vortex in a Maxwell fluid. Int J Nonlinear Mech 2003, 38: 985–990. 10.1016/S0020-7462(02)00042-2

    Article  MATH  MathSciNet  Google Scholar 

  7. Fetecau C, Zierep J, Angew Z: The RayleighStokes-problem for a Maxwell fluid. Math Phys 2003, 54: 1086–1093.

    MATH  Google Scholar 

  8. Hayat T, Nadeem S, Asghar S: Periodic unidirectional flows of a vis-coelastic fluid with fractional Maxwell model. Appl Math Comput 2004, 151: 153–161. 10.1016/S0096-3003(03)00329-1

    Article  MATH  MathSciNet  Google Scholar 

  9. Tan WC, Pan WX, Xu MY: A note on unsteady flows of a viscoelastic fluid with the fractional Maxwell model between two parallel plates. Int J Nonlinear Mech 2003, 38: 645–650. 10.1016/S0020-7462(01)00121-4

    Article  MATH  Google Scholar 

  10. Dagan G: Flow and Transport in Porous Formations. Springer-Verlag, Berlin; 1989.

    Chapter  Google Scholar 

  11. Ali AT: A note on the Exp-function method and its application to nonlinear equations. Phys Scr 2009., 79: 025006

    Google Scholar 

  12. El-Sabbagh MF, Ali AT: New generalized Jacobi elliptic function expansion method. Commun Nonl Sci Numer Simul 2008, 13: 1758–1766. 10.1016/j.cnsns.2007.04.014

    Article  MATH  MathSciNet  Google Scholar 

  13. Ali AT: New generalized Jacobi elliptic function rational expansion method. J Comput Appl Math 2011, 235(14):4117–4127. 10.1016/j.cam.2011.03.002

    Article  MATH  MathSciNet  Google Scholar 

  14. Ali AT: New exact solutions of Einstein vacuum equations for rotating axially symmetric fields. Phys Scr 2009., 79: 035006

    Google Scholar 

  15. Attallah SK, El-Sabbagh MF, Ali AT: Isovector fields and similarity solutions of Einstein vacuum equations for rotating fields. Commun Nonl Sci Numer Simul 2007, 12: 1153–1161. 10.1016/j.cnsns.2006.02.004

    Article  MATH  MathSciNet  Google Scholar 

  16. Bluman GW, Kumei S: Symmetries and differential equations. In Applied Mathematical sciences. Volume 81. Springer-Verlag, New York; 1989.

    Google Scholar 

  17. Olver PJ: Equivalence, Invariance and Symmetry. Cambridge University, Cambridge; 1995.

    Chapter  Google Scholar 

  18. Ovsiannikov LV: Group Analysis of Differential Equations. Cambridge ca-demic Press, New York; 1982.

    Google Scholar 

  19. Stephani H, MacCallum M: Differential equations: Their solutions using symmetries. Cambridge University Press, Cambridge; 1989.

    Google Scholar 

  20. Yürüsoy M, Pakdemirli M: Group classification of a non-Newtonian fluid model using classical approach and equivalence transformations. Int J Nonlinear Mech 1999, 34: 341–346. 10.1016/S0020-7462(98)00037-7

    Article  MATH  Google Scholar 

  21. Yürüsoy M, Pakdermirli M, Noyan OF: Lie group analysis of creeping flow of a second grade. Int J Nonlinear Mech 2001, 36: 955–960. 10.1016/S0020-7462(00)00060-3

    Article  MATH  Google Scholar 

  22. Yürüsoy M: Similarity solutions for creeping flow and heat transfer in second grade fluids. Int J Nonlinear Mech 2004, 39: 665–672. 10.1016/S0020-7462(03)00020-9

    Article  MATH  Google Scholar 

  23. Shahzad F, Sajid M, Hayat T, Ayub M: Analytic solution for flow of a micropolar fluid. Acta Mech 2007, 188: 93–102. 10.1007/s00707-006-0398-4

    Article  MATH  Google Scholar 

  24. Mekheimer KhS, Husseny SZA, Ali AT, Abo-Elkhair RE: Lie Group Analysis and Similarity Solutions for a Couple Stress Fluid with Heat Transfer. J Adv R Appl Math 2010, 2: 1–17.

    Article  MathSciNet  Google Scholar 

  25. Mekheimer KhS, Husseny SZA, Ali AT, Abo-Elkhair RE: Lie point symmetries and similarity solutions for an electrically conducting Jeffrey fluid. Phys Scr 2011., 83: 015017

    Google Scholar 

  26. Mekheimer KhS, Husseny SZA, Ali AT, Abo-Elkhair RE: Similarity Solution for Flow of a Micro-Polar Fluid Through a Porous Medium. Applications and Applied Mathematics 2011, 6: 2082–2093.

    MathSciNet  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Khaled Saad Mekheimer.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Khaled S. Mekheimer, Mostafa F El-Sabbagh and and Rabea E. Abo-Elkhair contributed to each part of this work equally.

All the authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Authors’ original file for figure 11

Authors’ original file for figure 12

Authors’ original file for figure 13

Authors’ original file for figure 14

Authors’ original file for figure 15

Authors’ original file for figure 16

Authors’ original file for figure 17

Authors’ original file for figure 18

Authors’ original file for figure 19

Authors’ original file for figure 20

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mekheimer, K.S., El-Sabbagh, M.F. & Abo-Elkhair, R.E. Lie group analysis and similarity solutions for hydro-magnetic Maxwell fluid through a porous medium. Bound Value Probl 2012, 15 (2012). https://doi.org/10.1186/1687-2770-2012-15

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2012-15

Keywords