ISRN Mathematical Physics

Volume 2013 (2013), Article ID 705296, 10 pages

http://dx.doi.org/10.1155/2013/705296

Research Article

A Note on the Unsteady Incompressible MHD Fluid Flow with Slip Conditions and Porous Walls

H. Zaman,¹ Z. Ahmad,¹ and M. Ayub²

¹Department of Mathematics, Islamia College Chartered University, Peshawar 25000, Pakistan

²Department of Mathematics, Quaid-i-Azam University, Islamabad 45320, Pakistan

Received 16 July 2013; Accepted 3 October 2013

Academic Editors: A. M. Gavrilik and A. Herrera-Aguilar

Copyright © 2013 H. Zaman et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This work is concerned with the influence of uniform suction or injection on flow and heat transfer analysis of unsteady incompressible magnetohydrodynamic (MHD) fluid with slip conditions. The resulting unsteady problem for velocity and heat transfer is solved by means of Laplace transform. The characteristics of the transient velocity, overall transient velocity, steady state velocity and heat transfer at the walls are analyzed and discussed. Graphical results reveal that the magnetic field, slip parameter, and suction (injection) have significant influences on the velocity, and temperature distributions, which also changes the heat transfer behaviors at the two plates. The results of Fang (2004) are also recovered by keeping magnetic field and slip parameter absent.

1. Introduction

Navier-Stokes equations are the basic equations of fluid mechanics. Exact solutions of Navier-Stokes equations are rare due to their inherent nonlinearity. Exact solutions are important because they serve as accuracy checks for numerical solutions. Complete integration of these equations is done by computer techniques, but the accuracy of the results can be established only by comparison with exact solutions. In the literature, there are a large number of Newtonian fluid flows for which exact solutions are possible [1–6]. The effects of transverse magnetic field on the flow of an electrically conducting viscous fluid have been studied extensively in view of numerous applications to astrophysical, geophysical, and engineering problems [7–15]. If the working fluid contains concentrated suspensions, then the wall slip can occur [16]. Khaled and Vafai [3] studied the effect of the slip on Stokes and Couette flows due to an oscillating wall. However, the literature lacks studies that take into account the possibility of fluid slippage at the walls. Applications of these problems appear in microchannels or nanochannels and in applications where a thin film of light oil is attached to the moving plates or when the surface is coated with special coating such as a thick monolayer of hydrophobic octadecyltrichlorosilane [17]. Yu and Ameel [18] imposed nonlinear slip boundary conditions on flow in rectangular microchannels. Erdogan [6] studied deeply the solution to the Stokes problem under nonslip conditions at the wall. Ayub and Zaman [19] studied the effects of suction and blowing for orthogonal flow impinging on a wall. Khan et al. [20] discussed the flow of Sisko fluid through a porous medium. Ariel et al. [21] considered the flow of elasticoviscous fluid with partial slip. Raptis et al. [22–25] studied steady and unsteady free convection and mass transfer flow through a porous medium. Penton [26] presented the transient solution for the flow due to the oscillating plate.

In this note, the flow of an incompressible, unsteady, viscous, MHD fluid with slip conditions is considered. Unsteady means time dependent flow, and we are looking for the effects of different parameters on flow with the variation of time. Unsteady and steady velocity profiles with mass transfer will be presented and solved exactly. There is mass injection from one plate and the same amount of suction on the other plate. The steady state temperature is also solved and discussed. When the fluid motion is set up from rest, the velocity field contains transients, determined by the initial conditions which gradually disappear in time. The effect of magnetic field and time on the transient velocity and on overall transient velocity has been seen graphically for both injection and suction. The effect of slip parameter on steady state velocity for injection/suction is shown graphically. Steady state temperature profiles and heat transfer rate at the walls (Nusselt number at the walls) are also discussed for injection/suction for different Prandtl and Reynolds numbers. The results of Fang [1] are recovered by taking magnetic field parameter and slip parameter to be zero.

2. Theoretical Derivation

2.1. Transient Velocity Profiles

Consider an incompressible, viscous, unsteady flow problem, in which there is slip between the bottom wall and fluid and also between top wall and fluid. There are mass injection velocityat the bottom wall and mass suction velocityat the top wall; corresponds to injection andcorresponds to suction. The governing equation for this problem can be obtained as [2]

whereis the velocity of the fluid in the-direction which is along the wall direction,is the distance alon -axis g,is the time,,,is the dynamic viscosity,is the kinematic viscosity,is the density of the fluid,is the electric conductivity of the fluid,is the applied magnetic field, andis the MHD factor or parameter. For the boundary conditions we consider the existence of slip between the velocity of the fluid at the walls and speed of the walls:

Initial condition is

whereis the slip parameter (gives the usual no slip condition at the wall) andis the velocity at the upper wall. The problem exactly reduces to the problem of Fang [1] if we take,Equations (1), (2), and (3) can be made dimensionless by defining

Then (1), (2), and (3) become

where,, and(Reynolds number). Decomposinginto two parts, say, transient part and a steady state part,

Then we have two separate problems and the steady state part will be

The solution of (10) with BCs (11) and (12) can be obtained as

whereand. When (no slip) and (no magnetic field), then (13) becomes

Equation (14) the result number (7) of Fang [1]. If there is no mass transfer at the walls, then, soand, and (13) becomes

If there is no mass transfer at the walls and magnetic field is absent, then,, and slip parameterthen (10), (11), and (12) will collapse into

Its solution is, which is conventional Couette flow. When,, and , then (13) gives, except at the bottom wall (at the bottom wallandWhen,, and , then (13) gives, except at the upper wall (at the upper wall and. The transient part problem becomes

The solution can be derived by using Laplace transformation techniques [27]. The Laplace transform pair is defined by

over that range of values offor which the integrals exist. Here,is a parameter, real or complex,is the operator that transformsinto, called Laplace transform operator, andis the inverse Laplace transform operator. The solution can be shown as

whereand. In inverse Laplace transform of the above equation we have simple pole atand infinite number of poles (located on the negative real axis) at,, where=and are given by

The transient part velocity will be

wherestands for residue andis given by

We have

where

The residue atgives steady velocity

Therefore,

By using (13), we get

The transient part velocity from (24) becomes

Therefore, the overall transient solution from (9) becomes

To recover Fang [1] result number (13) we substitute,,, andinto (32), and we get

Since, so (23) gives

Therefore; also we have, but, so

Using,, andinto (33) we obtain

Consider

So (36) becomes

Result number (13) of Fang [1] is exactly followed in (38).

2.2. Steady State Temperature

It is too difficult to exactly solve the transient energy equation with viscous dissipation for this problem by using. We solve the steady state energy equation, which in dimensionless form is given by [1]

whereis the Prandtl number,is the Eckert number, is temperature of the bottom fixed plate, is temperature of the top moving plate,is dimensionless temperature, and is the steady state velocity given by (13). The solution of (39) with BCs (35) for, is

where

When viscous dissipation term is negligible, then energy equation (39) is

Its solution is

So the first term in (40) is the solution of energy equation when viscous dissipation term is negligible and the second term in (40) is the temperature profile from viscous dissipation and MHD. It is found from (40) that the temperature profile is linearly dependent upon the Eckert number. The Nusselt number at the walls will be. The Nusselt number for the bottom wall, , is

The Nusselt number for the upper wall,, is

When, then the energy equation is

The solution of (46) is

where

The Nusselt number in this case is

It is worth recalling in the vicinity of (40) that when,, then result number (18) of Fang [1] can be derived.

3. Graphs and Discussion

In this part we discuss the variation of the transient part velocity, overall transient velocity, and steady state velocitywith distance from the wallfor different values of Reynolds number, magnetic field parameter, slip parameter , and time.

Figures 1 and 2 show the variation of transient part velocitywith distance from the wallfor several values of magnetic parameter, by keepingandfixed. Figure 1 shows that when there is mass suctionat the top wall, with increase in magnetic field, transient part velocity decreases in magnitude. Figure 2 shows that when there is mass injectionat the bottom wall, with increase in magnetic field, transient part velocitydecreases and will become weaker as compared to the case of suction. From Figure 3 it is observed that for suction at top wall and for fixed values of,transient part velocitydecreases in magnitude with increase in time. Figure 4 shows that for injection at bottom wall and for fixed values of,transient part velocitydecreases in magnitude with increase in time. From Figures 3 and 4 it is seen that the transient part velocity will decay with time, which is consistent with what we expected. From Figures 3 and 4 it is clear that, after a certain time, the transient part velocity will die away and velocity will become developed. Figures 5 and 6 indicate variation of overall transient velocitywithfor fixed values ofand. Figure 5 shows that, for, with increase in magnetic parameter, overall transient velocitydecreases. Figure 6 shows that, for, with increase in magnetic parameter, overall transient velocitydecreases but weaker. The overall transient velocities forat different times are depicted in Figure 7. Figure 7 shows that for mass injection at the bottom wall, overall transient velocityincreases with time.

Figure 1: Variation of  withfor several values ofwith,  , and .

Figure 2: Variation of  withfor several values ofwith,  , and .

Figure 3: Variation of  withfor several values of  with,  , and .

Figure 4: Variation of  withfor several values of  with,, and  .

Figure 5: Variation of  withfor several values of  with,, and  .

Figure 6: Variation of  withfor several values of  with,, and  .

Figure 7: Variation of  withfor several values ofwith,, and  .

Figures 8 and 9 illustrate the variation of steady state velocitywithfor several values of slip parameterand for fixed value of. Figureshows that, for suction at the top wall with increase in slip parameter, steady state velocityincreases. Figure 9 shows that, for injection at the bottom wall with increase in slip parameter, steady state velocitydecreases.

Figure 8: Variation of  withfor several values ofwith.

Figure 9: Variation of  withfor several values ofwith,.

Now we discuss the variation of the steady state temperature distributionwith distance from the surfacefor different values of five dimensionless parameters: Reynolds number, Prandtl number, Eckert number, slip parameter, and magnetic parameter. Variation of the Nusselt number at walls for different Prandtl numbers is also discussed. Figures 10 and 11 elucidate the variation of steady state temperature distributionwithfor several values ofand for fixed values of,, and. Figure 10 shows that for, (injection) temperatureat a point decreases with increase in. Figure 11 shows that, for, (suction) temperature with increase in increases at a point. Comparison of Figures 10 and 11 shows that, for mass suction at the bottom wall, the maximum temperature will exist in the fluid because of viscous dissipation.

Figure 10: Variation of  withfor several values ofwith,,, and .

Figure 11: Variation ofwithfor several values ofwith,,, and .

Figures 12 and 13 illustrate the variation of Nusselt numbers at the bottom wall with, for fixed values of,, andand for several values of. Figure 12 shows that when there is mass suction at the bottom wall, heat transfer rateincreases with increase in. Figure 13 shows that when there is mass injection at the bottom wall, heat transfer ratedecreases with increase in. Figures 14 and 15 depict the variation of Nusselt numbers at the top wall with, for fixed values of,, andand for several values of. Figure 14 shows that when there is mass suction at the top wall, heat transfer rateincreases in magnitude with increase in. Figure 15 shows that when there is mass injection at the top wall, heat transfer rateincreases with increase in.

Figure 12: Variation of  withfor several values ofwith,, and .

Figure 13: Variation of  withfor several values ofwith,, and  .

Figure 14: Variation of  withfor several values ofwith,, and .

Figure 15: Variation of  withfor several values ofwith,, and .

4. Final Remarks

In this study exact solutions for the velocity field and temperature field in the presence of magnetic field, porosity, and slip parameter are constructed. A uniform magnetic field is applied transversely to the flow. The expressions of the velocity field and temperature field for flow subjected to the slip conditions between the two parallel plates and fluid are obtained by means of Laplace transform. The so-obtained solutions, depending on the initial and the boundary conditions, are presented as sum of the steady state transient solutions. The results of Fang [1] are also recovered by taking ,. Graphical results for mass transfer reveal that it has significant influence on the velocity distribution, temperature distribution, and heat transfer rate at the walls. The current analysis will be useful in dealing with real engineering problems.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper. The authors do not have a direct financial relation with the commercial identity mentioned in the paper that might lead to a conflict of interests for any of them.

References

T. Fang, “A note on the incompressible couette flow with porous walls,” International Communications in Heat and Mass Transfer, vol. 31, no. 1, pp. 31–41, 2004.

H. Schlichting and K. Gersten, Boundary Layer Theory, Springer, Berlin, Germany, 8th edition, 2000.

A. R. A. Khaled and K. Vafai, “The effect of the slip condition on Stokes and Couette flows due to an oscillating wall: exact solutions,” International Journal of Non-Linear Mechanics, vol. 39, no. 5, pp. 795–809, 2004.

T. Hayat, M. Khan, A. M. Siddiqui, and S. Asghar, “Transient flows of a second grade fluid,” International Journal of Non-Linear Mechanics, vol. 39, no. 10, pp. 1621–1633, 2004.

M. Ahmad, H. Zaman, and N. Rehman, “Effects of Hall current on unsteady MHD flows of a second grade fluid,” Central European Journal of Physics, vol. 8, no. 3, pp. 422–431, 2010.

M. E. Erdogan, “Note on an unsteady flow of a viscous fluid due to an oscillating plane wall,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 1–6, 2000.

K. R. Crammer and S. I. Pai, Magnetofluid Dynamics for Engineers and Applied Physicists, McGraw-Hill, New York, NY, USA, 1973.

M. Ayub, H. Zaman, and M. Ahmad, “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, vol. 8, no. 1, pp. 135–149, 2010.

T. Hayat, T. Javed, and Z. Abbas, “MHD flow of a micropolar fluid near a stagnation-point towards a non-linear stretching surface,” Nonlinear Analysis: Real World Applications, vol. 10, no. 3, pp. 1514–1526, 2009.

T. Hayat, Z. Abbas, and T. Javed, “MHD stagnation point flow and heat transfer over a permeable surface through a porous space,” Journal of Porous Media, vol. 12, no. 2, pp. 183–195, 2009.

S. Asghar, M. Khan, and T. Hayat, “Magnetohydrodynamic transient flows of a non-newtonian fluid,” International Journal of Non-Linear Mechanics, vol. 40, no. 5, pp. 589–601, 2005.

M. Khan, S. Hyder Ali, T. Hayat, and C. Fetecau, “MHD flows of a second grade fluid between two side walls perpendicular to a plate through a porous medium,” International Journal of Non-Linear Mechanics, vol. 43, no. 4, pp. 302–319, 2008.

M. Khan, S. Asghar, and T. Hayat, “Hall effect on the pipe flow of a Burgers' fluid: an exact solution,” Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 974–979, 2009.

R. Cortell, “MHD flow and mass transfer of an electrically conducting fluid of second grade in a porous medium over a stretching sheet with chemically reactive species,” Chemical Engineering and Processing, vol. 46, no. 8, pp. 721–728, 2007.

T. Hayat, H. Zaman, and M. Ayub, “Analytical solution of hydromagnetic flow with Hall effect over a surface stretching with a power-law velocity,” Numerical Methods for Partial Differential Equations, vol. 27, no. 4, pp. 937–959, 2011.

F. Soltani and Ü. Yilmazer, “Slip velocity and slip layer thickness in flow of concentrated suspensions,” Journal of Applied Polymer Science, vol. 70, no. 3, pp. 515–522, 1998.

C. Derek, D. C. Tretheway, and C. D. Meinhart, “Apparent fluid slip at hydrophobic microchannel walls,” Physics of Fluids, vol. 14, no. 3, pp. L9–L12, 2002.

S. Yu and T. A. Ameel, “Slip-flow heat transfer in rectangular microchannels,” International Journal of Heat and Mass Transfer, vol. 44, no. 22, pp. 4225–4234, 2001.

M. Ayub and H. Zaman, “Analytical solution of orthogonal flow impinging on a wall with suction or blowing,” Journal of Basic & Applied Sciences, vol. 6, no. 2, pp. 93–97, 2010.

M. Khan, Z. Abbas, and T. Hayat, “Analytic solution for flow of Sisko fluid through a porous medium,” Transport in Porous Media, vol. 71, no. 1, pp. 23–37, 2008.

P. D. Ariel, T. Hayat, and S. Asghar, “The flow of an elastico-viscous fluid past a stretching sheet with partial slip,” Acta Mechanica, vol. 187, no. 1–4, pp. 29–35, 2006.

A. Raptis, C. Perdikis, and G. Tzivanidis, “Free convection flow through a porous medium bounded by a vertical surface,” Journal of Physics D, vol. 14, no. 7, article 001, pp. L99–L102, 1981.

A. Raptis, G. Tzivanidis, and N. Kafousias, “Free convection and mass transfer flow through a porous medium bounded by an infinite vertical limiting surface with constant suction,” Letters in Heat and Mass Transfer, vol. 8, no. 5, pp. 417–424, 1981.

A. Raptis, N. Kafousios, and C. Massalas, “Free convection and mass transfer flow through a porous medium bounded by an infinite vertical porous plate with constant heat flux,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 62, pp. 489–491, 1982.

A. A. Raptis, “Unsteady free convective flow through a porous medium,” International Journal of Engineering Science, vol. 21, no. 4, pp. 345–348, 1983.

R. Penton, “The transient for Stokes oscillating plane: a solution in terms of tabulated functions,” Journal of Fluid Mechanics, vol. 31, pp. 819–825, 1968.

G. E. Roberts and H. Kaufman, Tables of Laplace Transforms, WB Saunders, Philadelphia, Pa, USA, 1966.