﻿ Stokes First Problem for an Unsteady MHD Third-Grade Fluid in a Non-Porous Half Space with Hall Currents

Open Journal of Applied Sciences
Vol.4 No.3(2014), Article ID:43476,11 pages DOI:10.4236/ojapps.2014.43010

Stokes First Problem for an Unsteady MHD Third-Grade Fluid in a Non-Porous Half Space with Hall Currents

Haider Zaman, Arif Sohail, Ubaidullah

Faculty of Numerical Sciences, Islamia College University, Peshawar, Pakistan

Email: haiderzaman67@yahoo.com

Received 6 January 2014; revised 11 February 2014; accepted 20 February 2014

ABSTRACT

The well-known problem of unidirectional plane flow of a fluid in a non-porous half-space due to the impulsive motion of the rigid plane wall it rests upon is discussed in the context of an unsteady MHD third-grade fluid in presence of Hall currents. The governing non-linear partial differential equations describing the problem are converted to a system of non-linear ordinary differential equations by using the similarity transformations. The complex analytical solution is found by using the homotopy analysis method (HAM). The existing literature on the topic shows that it is the first study regarding the effects of Hall current on flow of an unsteady MHD third-grade fluid over an impulsively moving plane wall. The convergence of the obtained complex series solutions is carefully analyzed. The effects of dimensionless parameters on the velocity are illustrated through plots and the effects of the pertinent parameters on the local skin friction coefficient at the surface of the wall are presented numerically in tabular form.

1. Introduction

The present investigation is to analyze the Stokes first problem for an unsteady MHD third-grade fluid in a non-porous half space with Hall currents. The arising non-linear problem is solved by the homotopy analysis method (HAM), which is a novel technique and has been used by many researchers [37] -[45] . The method gives more realistic complex analytic solution which is uniformly valid for all values of the dimensionless time. The convergence region for the complex series solution is found with the help of. The effects of the material parameters of the third-grade fluid, Hall parameter, Hartmann number and homotopy parameter on the velocity and its time series are investigated for the impulsive motion of the wall. When we introduce Hall currents in the problem then velocity field becomes complex and it has two parts, real and imaginary. It is not possible to plot complex velocity as a whole on the real plane, but to analyze the flow we plot real and imaginary parts separately. Sometimes our graphs for real and imaginary parts of velocity look likes that they are not satisfying boundary conditions, but it is not true because we are not plotting the complete solution, we are taking either real part or imaginary part of the solution.

2. Mathematical Formulation of the Problem

We consider the unsteady flow of an electrically conducting incompressible third-grade fluid past a rigid plane wall coinciding with the plane. The fluid over the plane wall is initially at rest and it sets in motion due to the sudden jerked of the wall subjected to a uniform transverse magnetic field. It is assumed that the flow takes place in the upper half plane with the wall on the. The Cauchy stress tensor for a third-grade fluid is given as [46]

, (1)

where is the scalar pressure, is the identity tensor, is the coefficient of viscosity, , are the material parameters of third-grade fluid, and are the first three Rivlin-Ericksen tensors defined by [46]

, (2)

. (3)

The equations governing the magnetohydrodynamic flow with Hall effect are:

Velocity field:

, (4)

Continuity equation:

, (5)

Equation of motion:

, (6)

Maxwell equations:

(7)

Generalized Ohm’s law:

, (8)

where is the velocity component in the, is time, is the total magnetic field, is the applied magnetic field, is the induced magnetic field, is the current density, is the electrical conductivity of the fluid, is the electric field, is the magnetic permeability, is the fluid density, and are the cyclotron frequency and collision time of the electrons respectively. We assume that, the quantities, and are constants throughout the flow field, the magnetic field is normal to the velocity vector and the induced magnetic field is neglected compared with the imposed magnetic field, so that the magnetic Reynolds number is small [47] . We also suppose that and (where and are the cyclotron frequency and collision time for ions respectively). The velocity is zero far from the wall and the pressure is uniform, so we neglect the pressure gradient term as the flow is the only kept in motion by the motion of the wall. Continuity Equation (5) is identically satisfied by the velocity taken in Equation (4). Under the aforementioned assumptions the unsteady Stokes first problem with Hall currents become

, (9)

The boundary and initial conditions are

, (10)

where is the kinematic viscosity, is the Hall parameter and is the velocity of the wall. In order to non-dimensionalize the problem let us introduce the similarity transformations

, (11)

where is the dimensionless velocity function, is the dimensionless distance from the wall and is the dimensionless time. Equations (9)-(11) become

, (12)

(13)

where prime denotes differentiation with respect to, is dimensionless second-grade parameter, , are dimensionless third-grade fluid material parameters and is the dimensionless modified Hartmann number [47] . The local skin friction coefficient or fractional drag coefficient on the surface of the moving wall is

, (14)

Now using Equations (1)-(4) and (11) the Equation (14) can be written in dimensionless variables as

, (15)

3. Analytic Solution with Convergence Check

The boundary conditions (13) leads us to take base functions for the velocity as

(16)

The velocity can be expressed in terms of base functions as

, (17)

To start with the homotopy analysis method, due to the boundary conditions (13) it is reasonable to choose the initial guess approximation

, (18)

and the auxiliary linear operator

, (19)

with the property:

, (20)

where and are arbitrary constants. Following the HAM and trying higher iterations with the unique and proper assignment of the results converge to the exact solution:

, (21)

At, , , , , , using the symbolic computation software such as MATLAB, MAPLE, MATHEMATICA to successively obtain

, (22)

similarly, , and so on are calculated. The obtained values of, , , …, leads us to write

. (23)

By avoiding the detailed recurrence relations which we calculated but not presenting here due to space, the total complex analytic solution in compact form is

(24)

where from initial guess in Equation (18) we obtain

. (25)

Clearly Equation (24) contains the auxiliary parameter, which gives the convergence region and rate of approximation for the homotopy analysis method. For this purpose, the is plotted for the dimensionless velocity for the sudden jerked motion of the plane wall at. It is obvious from Figure 1 that the range for the admissible value for is. Our calculations depict that the series of the dimensionless velocity in Equation (24) converges in the whole region of and for.

4. Graphs, Table and Discussion

In order to discuss the effects of various parameters on the flow field the graphs are drawn for the variation of the dimensionless velocity with distance from the surface. These graphs are plotted for different values of the Hartmann number, Hall parameter, homotopy parameter, second-grade material parameter, third-grade material parameters and. Figures 2 to 7 are plotted in absence of Hall currents.

Figure 2 indicates the variation of the velocity with time at different distances from the wall. From Figure 2 it is observed that the velocity decays with the increase of the distance from the wall and also the fluid near the wall moves together with the wall with the same velocity as of the wall. Figure 3 shows that for fixed values of, , , and after the sudden jerked of the wall with increase in dimensionless time the velocity decreases at all points and decays with increase of time, boundary layer thickness decreases, which results in shear thinning effects, which are used in recognition of the fluid. Clearly the viscosity induces in the fluid with the passage of time and after some time the steady state is achieved which proves the uniform validity of the solution for all time in the whole spatial region. Figure 4 illustrates that with increase in Hartmann number the velocity decreases and boundary layer thickness also decreases at all points. When magnetic field is applied transverse to the fluid velocity then it gives rise to a drag-like or resistive force which slow down or suppress the motion of the fluid on the wall. This leads to a reduction in the velocity of the fluid as seen in Figure 4. Figure 5 demonstrate the velocity changes with distance from the wall

Figure 1. for.

Figure 2. Influence of on.

Figure 3. Influence of on.

Figure 4. Influence of on.

at given points and. It is clearly seen that the smaller the, the more rapidly the velocity changes very close to the wall and the more like a solid the third grade fluid behaves. After velocity decreases, indicating larger apparent viscosity for third-grade fluid. Figure 6 displays that for fixed values of,

Figure 5. Influence of on.

Figure 6. Influence of on.

Figure 7. Influence of on.

, , and with increase in third-grade material parameter the velocity and boundary layer thickness decreases at all points. Figure 7 illustrates that for fixed values of, , , and with increase in third-grade material parameter the velocity and boundary layer thickness decreases. Inclusion of Hall current makes the velocity field complex, so we plot real and imaginary parts of the velocity to see the effects of Hall current on the flow field. Figure 8 describes that for fixed values of, , , , and with increase in Hall parameter real part of the velocity decreases in magnitude, boundary layer thickness increases, which results in shear thickening effects. Figure 9 shows that for fixed values of, , , , and with increase in Hall parameter imaginary part of the velocity decreases in magnitude and boundary layer thickness increases.

Here the table is prepared for the variation of the absolute values of the skin friction coefficient given by Equation (15) for sudden jerked motion of the wall. The table shows the impact of the Hartmann number and Hall current with time on. It is observed from table that for fixed values of, , , and with increase in Hartmann number absolute value of the skin friction coefficient increases for all values of the time. Also for fixed values of, , , and with increase in Hall current absolute value of the skin friction coefficient decreases for all values of the time. For fixed values of, , , , and with increase in dimensionless time the absolute value of the skin friction coefficient increases at all points. When magnetic field is applied normal to the fluid velocity then it gives rise to a drag-like or resistive force which slow down or suppress the motion of the fluid on the moving surface. This leads to a reduction in the velocity of the fluid and flow rates. With the increase in the strength of the magnetic field the motion of the particulate suspension on the surface reduces due to which shear stress at the wall reduces with increase in Hall current, as observed in the Table 1.

5. Conclusion

The Stokes’ first problem of an unsteady MHD third grade fluid in a non-porous half space is discussed by tak-

Figure 8. Influence of on.

Figure 9. Influence of on.

Table 1. Absolute values of the skin friction coefficient with, , ,.

ing Hall currents into account. Analytical solution for the non-linear problem is given and convergence of the solution is appropriately discussed. The non-linear effect on the velocity is shown and discussed. It reveals that plots for shear thickening/shear thinning behavior of a fluid are dependent upon the rheological properties of the fluid. The numerical results for Hall parameter, Hartmann number, second-grade parameter, thirdgrade material parameters, and dimensionless time reveal that Hall parameter, Hartmann number and dimensionless time have significant influence on the local skin-friction coefficient on the surface of the impulsively moving wall.

References

1. Stokes, G.G. (1851) On the Effect of the Internal Friction of Fluids on the Motion of Pendulums. Transactions of the Cambridge Philosophical Society, Part II, 9, 8-106.
2. Tanner, R.I. (1962) Note on the Rayleigh Problem for a Visco-Elastic Fluid. Zeitschrift für Angewandte Mathematik und Physik, 13, 573-580. http://dx.doi.org/10.1007/BF01595580
3. Soundalgekar, V.M. (1974) Stokes Problem for Elastico-Viscous Fluid. Rheologica Acta, 13, 177-179. http://dx.doi.org/10.1007/BF01520872
4. Teipel, I. (1981) The Impulsive Motion of a Flat Plate in a Viscoelastic Fluid. Acta Mechanica, 39, 277-279. http://dx.doi.org/10.1007/BF01170349
5. Puri, P. (1984) Impulsive Motion of a Flat Plate in a Rivlin-Ericksen Fluid. Rheologica Acta, 23, 451-453. http://dx.doi.org/10.1007/BF01329198
6. Erdogan, M.E. (1995) Plane Surface Suddenly Set in Motion in a Non-Newtonian Fluid. Acta Mechanica, 108, 179- 187. http://dx.doi.org/10.1007/BF01177337
7. Zeng, Y. and Weinbaum, S. (1995) Stokes Problems for Moving Half-Planes. Journal of Fluid Mechanics, 287, 59-74. http://dx.doi.org/10.1017/S0022112095000851
8. Tan, W.C. and Xu, M.Y. (2002) The Impulsive Motion of Flat Plate in a Generalized Second Grade Fluid. Mechanics Research Communications, 29, 3-9. http://dx.doi.org/10.1016/S0093-6413(02)00223-9
9. Tan, W.C. and Xu, M.Y. (2002) Plane Surface Suddenly Set in Motion in a Viscoelastic Fluid with Fractional Maxwell Model. Acta Mechanica Sinica, 18, 342-349. http://dx.doi.org/10.1007/BF02487786
10. Erdogan, M.E. (2002) On the Unsteady Unidirectional Flows Generated by Impulsive Motion of a Boundary or Sudden Application of a Pressure Gradient. International Journal of Non-Linear Mechanics, 37, 1091-1106. http://dx.doi.org/10.1016/S0020-7462(01)00035-X
11. Fetecau, C. and Fetecau, C. (2003) The First Problem of Stokes for an Oldroyd-B Fluid. International Journal of NonLinear Mechanics, 38, 1539-1544. http://dx.doi.org/10.1016/S0020-7462(02)00117-8
12. Tan, W.C. and Masuoka, T. (2005) Stokes’ First Problem for an Oldroyd-B Fluid in a Porous Half Space. Physics of Fluids, 17, 023101. http://dx.doi.org/10.1063/1.1850409
13. Tan, W.C. and Masuoka, T. (2005) Stokes’ First Problem for a Second Grade Fluid in a Porous Half-Space with Heated Boundary. International Journal of Non-Linear Mechanics, 40, 515-522. http://dx.doi.org/10.1016/j.ijnonlinmec.2004.07.016
14. Zierep, J. and Fetecau, C. (2007) Energetic Balance for the Rayleigh-Stokes Problem of a Second Grade Fluid. International Journal of Engineering Science, 45, 155-162. http://dx.doi.org/10.1016/j.ijengsci.2006.09.001
15. Zierep, J., Bohning, R. and Fetecau, C. (2007) Rayleigh-Stokes Problem for Non-Newtonian Medium with Memory. Zeitschrift für Angewandte Mathematik und Mechanik, 87, 462-467. http://dx.doi.org/10.1002/zamm.200710328
16. Vieru, D., Nazar, M., Fetecau, C. and Fetecau, C. (2008) New Exact Solutions Corresponding to the First Problem of Stokes for Oldroyd-B Fluids. Computers & Mathematics with Applications, 55, 1644-1652. http://dx.doi.org/10.1016/j.camwa.2007.04.040
17. Shahzad, F., Hayat, T. and Ayub, M. (2008) Stokes’ First Problem for the Rotating Flow of a Third Grade Fluid. Nonlinear Analysis: Real World Applications, 9, 1794-1799. http://dx.doi.org/10.1016/j.nonrwa.2007.05.008
18. Hayat, T., Shahzad, F., Ayub, M. and Asghar, S. (2008) Stokes’ First Problem for a Third Grade Fluid in a Porous Half Space. Communications in Nonlinear Science and Numerical Simulation, 13, 1801-1807. http://dx.doi.org/10.1016/j.cnsns.2007.04.015
19. Fakhari, K., Zainal, A.A. and Kara, A.H. (2011) A Note on the Interplay between Symmetries, Reduction and Conservation Laws of Stokes’ First Problem for Third-Grade Rotating Fluids. Pramana Journal of Physics, 77, 439-445. http://dx.doi.org/10.1007/s12043-011-0164-6
20. Sajid, M., Ali, N., Javed T. and Abbas, Z. (2010) Stokes’ First Problem for a MHD Third Grade Fluid in a Porous Half Space. Journal of Porous Media, 1, 279-284.
21. Sato, H. (1961) The Hall Effects in the Viscous Flow of Ionized Gas between Parallel Plates under Transverse Magnetic Field. Journal of the Physical Society of Japan, 16, 1427-1433. http://dx.doi.org/10.1143/JPSJ.16.1427
22. Cramer, K. and Pai, S. (1973) Magnetofluid Dynamics for Engineers and Applied Physicists. McGraw-Hill, New York.
23. Ayub, M., Zaman, H. and Ahmad, M. (2010) Series Solution of Hydromagnetic Flow and Heat Transfer with Hall Effect in a Second Grade Fluid over a Stretching Sheet. Central European Journal of Physics, 8, 135-149. http://dx.doi.org/10.2478/s11534-009-0110-0
24. Ahmad, M., Zaman, H. and Rehman, N. (2010) Effects of Hall Current on Unsteady MHD Flows of a Second Grade Fluid. Central European Journal of Physics, 8, 422-431. http://dx.doi.org/10.2478/s11534-009-0083-z
25. Zaman, H. (2013) Hall Effects on the Unsteady Incompressible MHD Fluid Flow with Slip Conditions and Porous Walls. Applied Mathematics and Physics, 1, 31-38.
26. Hayat, T., Zaman, H. and Ayub, M. (2010) Analytic Solution of Hydromagnetic Flow with Hall Effect over a Surface Stretching with a Power Law Velocity. Numerical Methods for Partial Differential Equations, 27, 937-959.
27. Hayat, T., Naz, R. and Asghar, S. (2004) Hall Effects on Unsteady Duct Flow of a Non-Newtonian Fluid in a Porous Medium. Applied Mathematics and Computation, 157, 103-114. http://dx.doi.org/10.1002/num.20562
28. Hayat, T. and Nawaz, M. (2011) Hall and Ion-Slip Effects on Three-Dimensional Flow of a Second Grade Fluid. International Journal for Numerical Methods in Fluids, 66, 183-193. http://dx.doi.org/10.1016/j.amc.2003.08.069
29. Pop I. and Soundalgekar, V.M. (1974) Effects of Hall Current on Hydromagnetic Flow near a Porous Plate. Acta Mechanica, 20, 315-318. http://dx.doi.org/10.1007/BF01175933
30. Gupta, A.S. (1975) Hydromagnetic Flow Past a Porous Flate Plate with Hall Effects. Acta Mechanica, 22, 281-287. http://dx.doi.org/10.1007/BF01170681
31. Debnath, L., Ray, S.C. and Chatterjee, A.K. (1979) Effects of Hall Current on Unsteady Hydromagnetic Flow Past a Porous Plate in a Rotating Fluid System. Zeitschrift für Angewandte Mathematik und Mechanik, 59, 469-471. http://dx.doi.org/10.1002/zamm.19790590910
32. Katagiri, M. (1969) The Effect of Hall Currents on the Magnetohydrodynamic Boundary Layer Flow Past a Semi-Infinite Flate Plate. Journal of the Physical Society of Japan, 27, 1051-1059. http://dx.doi.org/10.1143/JPSJ.27.1051
33. Abo-Eldahab, E.M. and Elbarbary, M.E. (2001) Hall Current Effect on Magnetohydrodynamic Free Convection Flow Past a Semi-Infinite Vertical Plate with Mass Transfer. International Journal of Engineering Science, 39, 1641-1652. http://dx.doi.org/10.1016/S0020-7225(01)00020-9
34. Abo-Eldahab, E.M. and Salem, A.M. (2004) Hall Effects on MHD Free Convection Flow of a Non-Newtonian Power Law Fluid at a Stretching Surface. International Communications in Heat and Mass Transfer, 31, 343-354. http://dx.doi.org/10.1016/j.icheatmasstransfer.2004.02.005
35. Attia, H.A. (2006) Hall Effects on the Flow of a Dusty Bingham Fluid in a Circular Pipe. Turkish Journal of Engineering and Environmental Sciences, 30, 14-21.
36. Attia, H.A. (1998) Hall Current Effects on the Velocity and Temperature Fields of an Unsteady Hartman Flow. Canadian Journal of Physics, 76, 739-746.
37. Liao, S.J. (2003) Beyond Perturbation: Introduction to Homotopy Analysis Method. Chapman and Hall/CRC Press, Florida. http://dx.doi.org/10.1201/9780203491164
38. Liao, S.J. (1992) The Proposed Homotopy Analysis Technique for the Solution of Nonlinear Problem. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai.
39. Liao, S.J. (2012) Homotopy Analysis Method in Nonlinear Differential Equations. Springer-Verlag, Berlin. http://dx.doi.org/10.1007/978-3-642-25132-0
40. Liao, S.J. (2013) Advances in the Homotopy Analysis Method. World Scientific Publishing Company, Singapore.
41. Liao, S.J. (2009) Notes on the Homotopy Analysis Method: Some Definitions and Theorems. Communications in Nonlinear Science and Numerical Simulation, 14, 983-997.
42. Liao, S.J. (2006) An Analytic Solution of Unsteady Boundary-Layer Flows Caused by an Impulsively Stretching Plate. Communications in Nonlinear Science and Numerical Simulation, 11, 326-339.
43. Ayub, M., Zaman, H., Sajid, M. and Hayat, T. (2008) Analytical Solution of Stagnation-Point Flow of a Viscoelastic Fluid towards a Stretching Surface. Communications in Nonlinear Science and Numerical Simulation, 13, 1822-1835. http://dx.doi.org/10.1016/j.cnsns.2007.04.021
44. Zaman, H. and Ayub, M. (2010) Series Solution of Unsteady Free Convection Flow with Mass Transfer along an Accelerated Vertical Porous Plate with Suction. Central European Journal of Physics, 8, 931-939. http://dx.doi.org/10.2478/s11534-010-0007-y
45. Zaman, H., Hayat, T., Ayub, M. and Gorla, R.S.R. (2011) Series Solution for Heat Transfer from a Continuous Surface in a Parallel Free Stream of Viscoelastic Fluid. Numerical Methods for Partial Differential Equations, 27, 1511-1524. http://dx.doi.org/10.1002/num.20593
46. Rivlin, R.S. and Ericksen, J.L. (1955) Stress Deformation Relations for Isotropic Materials. Journal of Rational Mechanics and Analysis, 4, 323-425.
47. Sutton, G.W. and Sherman, A. (1965) Engineering Magnetohydrodynamics. McGraw-Hill, New York.