Abstract
The equations of two dimensional incompressible fluid flow for hydromagnetic 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.
Keywords:
lie point symmetries; similarity solutions; Maxwell fluid; porous medium; MHD1 Introduction
NonNewtonian 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 shearthinning and shearthickening 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 nonNewtonian 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 nonNewtonian 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 nonNewtonian 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 electroncyclotron frequency and the electronatomcollision 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 magneticforce 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 socalled 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 nonNewtonian 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 twodimensional 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 micropolar fluid through a porous medium [26].
From discussion above, we attend to find the analytical (similarity) solutions for the flow problem of an incompressible hydromagnetic Maxwell fluid through a porous medium using Lie group analysis. The problem is presented as follows, in Section 2, the equations governing twodimensional 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 twodimensional motion of an incompressible hydromagnetic Maxwell fluid through a porous medium can be written as:
where
Using the following dimensionless parameters
the system (1) becomes
where
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
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
with summation convention over the repeated index and x_{1 }= x, x_{2 }= y, x_{3 }= t, u_{1 }= u, u_{2 }= v, u_{3 }= 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:
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
where
and the operator
where D_{ij }= D_{i}(D_{j}) = D_{j}(D_{i}) = D_{ji }and
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 overdetermined equations. Solving these set of determining equations we obtain the required components of isovector field as follows:
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 a_{1 }= 0. The characteristic equations corresponding to the translational symmetry are:
By solving the ODEs (11), we can obtain the similarity variables and similarity functions as follows:
where
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:
where b_{i}, i = 1, 2 are arbitrary constants and β(ϕ, ψ) an arbitrary function that satisfy two conditions:
By solving the characteristic equations,
if we take β(ϕ, ψ) = 0 then we can obtain the similarity variable and similarity functions as the following:
where
Integrating the first equation in (19) yields
where c_{1 }is an arbitrary constant. Eliminating h(ξ) from the second and third equations in (19) along with Eq. (20) we get the following equation:
where
By solving equation (21) we get
where
From equations (20) and (23) the expression of the function g(χ) becomes
and from the second and third equations in (19) we get
where c_{4 }and c_{5 }are arbitrary constants.
In the form of the original variables, our exact solutions can be written as follows:
3.2 Rotational symmetry
In this section, the parameter a_{1 }is taken to be an arbitrary nonzero constant. The characteristic equations corresponding to the rotational symmetry are:
Integrating equations (28) using the Lie group analysis method, we get the rotationally invariant solutions for our problem in the following form:
where
Substituting the new (similarity) variables (ϕ, ψ) and functions (G_{1}, G_{2}, G_{3}) into the original system (3) yields the following set of equations:
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:
where d_{1 }is an arbitrary constant and η(ψ,t,G_{1}) an arbitrary function that satisfy the following condition:
If we take η(ψ,t,G_{1}) = 0 is a simple solution of Eq. (32) then the characteristic equations are:
The similarity variables and resulting functions are
One now substitutes the similarity variable and the functions into the equations (30) and obtains
Integrating the first equation in the system (35) we get:
where g_{1 }is an arbitrary constant. From the second and the third equations in (35) along with (36) we get:
by integrating equation (38) we obtain
and
where g_{2 }and g_{3 }are arbitrary constants, and the function _{p}F_{q }and _{p}F_{q }Regularized are defined as:
Then the solution of our problem in the original variables using symbolic computations is:
where ψ and ϕ are the same in (29).
4 Solutions for hydromagnetic 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 hydromagnetic 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:
where U_{0 }is the velocity of the plate, V_{0 }is the magmatic fluid velocity penetrating into the plate, P_{0 }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
where W = m_{0}m_{1 } m_{2}, α_{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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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. Variation of the dimensionless velocity distribution along the xaxis 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. Variation of the dimensionless velocity distribution along the yaxis 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 NonNewtonian 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 U_{0 }and V_{0 }(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 U_{0 }and V_{0.}
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 P_{0 }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. 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. 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. 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. 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. 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. 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 hydromagnetic 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 nonNewtonian 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 magnetomagma 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.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
Khaled S. Mekheimer, Mostafa F ElSabbagh and and Rabea E. AboElkhair contributed to each part of this work equally.
All the authors read and approved the final manuscript.
References

Zhaosheng, Y, Jianzhong, L: Numerical research on the coherent structure in the viscoelastic secondorder mixing layers. Appl Math Mech. 19, 671–677 (1998)

Shifang, H: Constitutive Equation and Computational Analytical Theory of NonNewtonian Fluids. Science Press, Beijing (2000)

Zhaosheng, JE, Rajagopal, KR: Fluids of differential type: critical review and thermodynamic analysis. Int J Eng Sci. 33, 689–729 (1995). Publisher Full Text

Aksel, N: A brief note from the editor on the secondorder fluid. Acta. Mech. 157, 235–236 (2002). Publisher Full Text

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

Fetecau, C, Fetecau, C: Decay of a potential vortex in a Maxwell fluid. Int J Nonlinear Mech. 38, 985–990 (2003). Publisher Full Text

Fetecau, C, Zierep, J, Angew, Z: The RayleighStokesproblem for a Maxwell fluid. Math Phys. 54, 1086–1093 (2003)

Hayat, T, Nadeem, S, Asghar, S: Periodic unidirectional flows of a viscoelastic fluid with fractional Maxwell model. Appl Math Comput. 151, 153–161 (2004). Publisher Full Text

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. 38, 645–650 (2003). Publisher Full Text

Dagan, G: Flow and Transport in Porous Formations. SpringerVerlag, Berlin (1989)

Ali, AT: A note on the Expfunction method and its application to nonlinear equations. Phys Scr. 79, (2009) 025006

ElSabbagh, MF, Ali, AT: New generalized Jacobi elliptic function expansion method. Commun Nonl Sci Numer Simul. 13, 1758–1766 (2008). Publisher Full Text

Ali, AT: New generalized Jacobi elliptic function rational expansion method. J Comput Appl Math. 235(14), 4117–4127 (2011). Publisher Full Text

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

Attallah, SK, ElSabbagh, MF, Ali, AT: Isovector fields and similarity solutions of Einstein vacuum equations for rotating fields. Commun Nonl Sci Numer Simul. 12, 1153–1161 (2007). Publisher Full Text

Bluman, GW, Kumei, S: Symmetries and differential equations. Applied Mathematical sciences, SpringerVerlag, New York (1989)

Olver, PJ: Equivalence, Invariance and Symmetry. Cambridge University, Cambridge (1995)

Ovsiannikov, LV: Group Analysis of Differential Equations. Cambridge cademic Press, New York (1982)

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

Yürüsoy, M, Pakdemirli, M: Group classification of a nonNewtonian fluid model using classical approach and equivalence transformations. Int J Nonlinear Mech. 34, 341–346 (1999). Publisher Full Text

Yürüsoy, M, Pakdermirli, M, Noyan, OF: Lie group analysis of creeping flow of a second grade. Int J Nonlinear Mech. 36, 955–960 (2001). Publisher Full Text

Yürüsoy, M: Similarity solutions for creeping flow and heat transfer in second grade fluids. Int J Nonlinear Mech. 39, 665–672 (2004). Publisher Full Text

Shahzad, F, Sajid, M, Hayat, T, Ayub, M: Analytic solution for flow of a micropolar fluid. Acta Mech. 188, 93–102 (2007). Publisher Full Text

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

Mekheimer, KhS, Husseny, SZA, Ali, AT, AboElkhair, RE: Lie point symmetries and similarity solutions for an electrically conducting Jeffrey fluid. Phys Scr. 83, (2011) 015017

Mekheimer, KhS, Husseny, SZA, Ali, AT, AboElkhair, RE: Similarity Solution for Flow of a MicroPolar Fluid Through a Porous Medium. Applications and Applied Mathematics. 6, 2082–2093 (2011)