** Engineering ** Vol. 3 No. 1 (2011) , Article ID: 3722 , 4 pages DOI:10.4236/eng.2011.31005

Time Dependent Pressure Gradient Effect on Unsteady MHD Couette Flow and Heat Transfer of a Casson Fluid

Department of Engineering Mathematics and Physics, Faculty of Engineering,

Fayoum University, Fayoum, Egypt

E-mail: mes00@fayoum.edu.eg

Received November 21, 2009; revised July 28, 2010; accepted August 3, 2010

**Keywords:** MHD Flow, Heat Transfer, Non-Newtonian Fluids, Unsteady Pressure, Numerical Solution

ABSTRACT

The unsteady magnetohydrodynamic flow of an electrically conducting viscous incompre-ssible non-Newtonian Casson fluid bounded by two parallel non-conducting porous plates has been studied with heat transfer considering the Hall effect. The fluid is acted upon by a uniform and exponential decaying pressure gradient. An external uniform magnetic field is applied perpendicular to the plates and the fluid motion is subjected to a uniform suction and injection. The lower plate is stationary and the upper plate is suddenly set into motion and simultaneously suddenly isothermally heated to a temperature other than the lower plate temperature. Numerical solutions are obtained for the governing momentum and energy equations taking the Joule and viscous dissipations into consideration. The effect of unsteady pressure gradient, the Hall term, the parameter describing the non-Newtonian behavior on both the velocities and temperature distributions have been studied.

1. Introduction

The study of Couette flow in a rectangular channel of an electrically conducting viscous fluid under the action of a transversely applied magnetic field has immediate applications in many devices such as magnetohydrodynamic (MHD) power generators, MHD pumps, accelerators, aerodynamics heating, electrostatic precipitation, polymer technology, petroleum industry, purification of crude oil and fluid droplets sprays. Channel flows of a Newtonian fluid with heat transfer were studied with or without Hall currents by many authors [1-10]. These results are important for the design of the duct wall and the cooling arrangements. The most important nonNewtonian fluid possessing a yield value is the Casson fluid, which has significant applications in polymer processing industries and biomechanics. Casson fluid is a shear thinning liquid which has an infinite viscosity at a zero rate of shear. Casson’s constitutive equation represents a nonlinear relationship between stress and rate of strain and has been found to be accurately applicable to silicon suspensions, suspensions of bentonite in water and lithographic varnishes used for printing inks [11-13]. Many authors [14-20] studied the flow or/and heat transfer of a non-Newtonian fluids in different geometries. The effect of time dependent pressure gradient on unsteady dusty fluid was studied by Rukmangadachari [21] in a rectangular duct and Gireesha et al. [22] in a nholonomic coordinates system.

Attia [10] studied the influence of the Hall current on the velocity and temperature fields of an unsteady Hartmann flow of a conducting Newtonian fluid between two infinite non-conducting horizontal parallel and porous plates. Attia and Sayed-Ahmed [17] studied the influence of the Hall current on the velocity and temperature fields of an unsteady Couettee flow of a conducting Casson fluid between two infinite non-conducting horizontal parallel plates with constant pressure gradient. The extension of such problem to the case of Couettee flow of non-Newtonian Casson fluid has been done in the present study. The upper plate is moving with a uniform velocity while the lower plate is stationary. The fluid is acted upon by an exponentially decaying pressure gradient, uniform suction and injection from above and below, respectively. The fluid is also, subjected to a uniform magnetic field perpendicular to the plates. The Hall current is taken into consideration while the induced magnetic field is neglected by assuming a very small magnetic Reynolds number [5]. The two plates are kept at two different but constant temperatures. This configuration is a good approximation of some practical situations such as heat exchangers, flow meters, and pipes that connect system components. The Joule and viscous dissipations are taken into consideration in the energy equation. The governing momentum and energy equations are solved numerically using the finite difference approximations. The inclusion of unsteady pressure gradient, the Hall current, the suction and injection, and the non-Newtonian fluid characteristics leads to some interesting effects on both the velocity and temperature fields.

2. Formulation of the Problem

The geometry of the problem is shown in Figure 1. The fluid is assumed to be laminar, incompressible and obeying a Casson model and flows between two infinite horizontal plates located at the y = ± h planes and extend from x = -¥ to ¥ and from z = -¥ to ¥. The upper plate is suddenly set into motion and moves with a uniform velocity U_{o} while the lower plate is stationary. The upper plate is simultaneously subjected to a step change in temperature from T_{1 }to T_{2}. Then, the upper and lower plates are kept at two constant temperatures T_{2} and T_{1} respectively, with T_{2} > T_{1}. The fluid is acted upon by an exponentially decaying pressure gradient ∂p/∂x in the x-direction, and a uniform suction from above and injection from below which are applied at t = 0. A uniform magnetic field B_{o} is applied in the positive y-direction and is assumed undisturbed as the induced magnetic field is neglected by assuming a very small magnetic Reynolds number. The Hall effect is taken into consideration and consequently a z-component for the velocity is expected to arise. The uniform suction implies that the y-component of the velocity v_{0} is constant. Thus, the fluid velocity vector is given by,

The fluid motion starts from rest at t = 0, and the no-slip condition at the plates in z-direction implies that the fluid velocity has no z-component at y = ± h. The initial temperature of the fluid is assumed to be equal to T_{1}. Since the plates are infinite in the x and zdirections, the physical quantities do not change in these directions.

The flow of the fluid is governed by the momentum equation

(1)

Figure 1. Geometry of the problem.

where r is the density of the fluid and m is the apparent viscosity of the model and is given by

(2)

where K_{c}^{2} is the Casson’s coefficient of viscosity and t_{o} is the yield stress. If the Hall term is retained, the current density J is given by

(3)

where s is the electric conductivity of the fluid and b is the Hall factor [5]. Equation (3) may be solved in J to yield

(4)

where m is the Hall parameter and m = _{o}. Thus, the two components of the momentum Equation (1) read

(5)

(6)

where is the unsteady pressure gradient The energy equation with viscous and Joule dissipations is given by

(7)

where c_{p} and k are, respectively, the specific heat capacity and the thermal conductivity of the fluid. The second and third terms on the right-hand side represent the viscous and Joule dissipations respectively. We notice that each of these terms has two components. This is because the Hall effect brings about a velocity w in the z-direction. The initial and boundary conditions of the problem are given by

at t £ 0, and at y =-h and y =h for t > 0

u = 0 at y =-h for t > 0, u = U_{o} at y =h for t > 0, (8)

T = T_{1} at t £ 0, T = T_{2} at y =h and T = T_{1} at y = -h for t > 0 (9)

It is expedient to write the above equations in the non-dimensional form. To do this, we introduce the following non-dimensional quantities

is the constant pressure gradient.

is the decaying parameter in the unsteady pressure gradient

is the Casson number (dimensionless yield stress)

is the Reynolds number,

is the suction parameter,

is the Prandtl number,

is the Eckert number,

is the Hartmann number squared

In terms of the above non-dimensional variables and parameters Eqs.(5-9) and (2) are, respectively, written as (where the hats are dropped for convenience);

(10)

(11)

(12)

(13)

(14)

(15)

where is the constant pressure gradient and a is the decaying parameter.

3. Numerical Solution

Equations (10,11,15) represent coupled system of nonlinear partial differential equations which are solved numerically under the initial and boundary conditions (13) using the finite difference approximations. A linearization technique is first applied to replace the nonlinear terms at a linear stage, with the corrections incorporated in subsequent iterative steps until convergence is reached. Then the Crank-Nicolson implicit method is used at two successive time levels [23]. An iterative scheme is used to solve the linearized system of difference equations. The solution at a certain time step is chosen as an initial guess for next time step and the iterations are continued till convergence, within a prescribed accuracy. Finally, the resulting block tridiagonal system is solved using the generalized Thomas-algorithm [23]. The energy Equation (12) is a linear non-homogeneous second-order partial differential equation whose righthand side is known from the solutions of the flow Equations (10,11,15) subject to the conditions (13). The values of the velocity components are substituted in the right-hand side of Equation (12) which is solved numerically with the initial and boundary conditions (14) using central differences and Thomas-algorithm to obtain the temperature distribution. Finite difference equations relating the variables are obtained by writing the equations at the mid-point of the computational cell and then replacing the different terms by their second order central difference approximations in the y-direction. The diffusion terms are replaced by the average of the central differences at two successive time-levels. The computational domain is divided into meshes of dimension Dt and Dy in time and space respectively as shown in Figure 2. We define the variables and to reduce the second order differential Equations (10,11,12) to first order differential equations. The finite difference representations for the resulting first order differential Equations (10,11) take the form (see Equations (16,17)).

where

and

The variables with bars are given initial guesses from the previous time step and an iterative scheme is used at every time to solve the linearized system of difference equations. Then the finite difference form for the energy Equation (12) can be written as

(18)

where DISP represents the Joule and viscous dissipation terms which are known from the solution of the momentum equations and can be evaluated at the mid point (i+1/2,j+1/2) of the computational cell. Computations have been made for α = 5, Pr = 1, Re = 1, Ha = 3 and Ec = 0.2. Grid-independence studies show that the computational domain 0

Figure 2. Mesh layout.

tions of the flow and heat parameters. Smaller step sizes do not show any significant change in the results. Convergence of the scheme is assumed when all of the unknowns u, v, w, B, θ and H for the last two approximations differ from unity by less than 10^{-6} for all values of y in –1 < y < 1 at every time step. Less than 7 approximations are required to satisfy this convergence criteria for all ranges of the parameters studied here.

4. Results and Discussions

Figures 3-5 show the variation of the velocity components u and w and the temperature θ at the central

(a)(b)(c)

Figure 3. Effect of decaying parameter a on u at y = 0 for various values of τ_{D} (m = 3, S = 1).

plane of the channel (y = 0) with time. These figures show the results for various values of the decaying parameter a = 0, 1 and 2 and for yield stress t_{D }= 0.0, 0.05 and 0.1. In these figures S = 1 and m = 3. Figure 3 shows that u decreases with increasing a for all values of t_{D}. It is observed also that the time at which u reaches its steady state value decreases with increasing a for a > 0 but that occurs earlier for constant pressure gradient (a = 0). Increasing t_{D} increases u for all values of a but with

(a)(b)(c)

Figure 4. Effect of decaying parameter a on w at y = 0 for various values of τ_{D} (m = 3, S = 1).

small differences.

In Figure 4, the velocity component w decreases with increasing a. This figure indicates that the influence of t_{D} on w depends on t and becomes more clear when the decaying parameter a = 0 but this influence is small for large a. It is observed that increasing t_{D} more decreases w for a = 0. Figure 5 shows that the influence of a on θ depends on t. It is observed that increasing a decreases θ

(a)(b)(c)

Figure 5. Effect of decaying parameter a on the time development of q at y = 0 for various values of τ_{D}_{ }(m = 3, S = 1).

while it is not greatly affected by changing t_{D}. The figure shows also that the time at which θ reaches its steady state value decreases with increasing a while it is not greatly affected by changing t_{D}.

Figures 6-8 show the profiles of the velocity components u and w and the temperature θ, respectively, for various values of time a and for t = 0.2, 1, and 2. The figures are evaluated for m =3, t_{D} = 0.05 and S = 1. It is clear from Figures 6 and 7 that the effect of decaying

(a) t = 0.2(b) t = 1(c) t = 2

Figure 6. Effect of decaying parameter a on the time development of the velocity u for various values of t (m = 3, S = 1, τ_{D}_{ }= 0.05 ).

parameter a on u and w depends on t and y. Figure 6 shows that, for small t, increasing a decreases u (with small differences) for all values of y and a. It is also observed that increasing a (a > 0) decreases u for all values of y with significant differences at medium and large t. It is also observed that the constant pressure gradient (a = 0) is greatly different from unsteady pressure gradients (a > 0). This can be attributed to the fact that increasing a will decrease the pressure gradient which mainly generates

(a) t = 0.2(b) t = 1(c) t = 2

Figure 7. Effect of decaying parameter a on the distribution of w with y for various values of t (m = 3, S = 1, τ_{D }= 0.05 ).

the velocity u. Figure 7 shows that, for small t, increasing a increases w (with no significant differences ) for all y. For large t, increasing a decreases w (with no signifycant difference between a=1 and a=2) for all y. The figures show also that the velocity components u and w do not reach their steady state monotonically. The velocities u and w increase with time up till a maximum value and then decrease up to the steady state. Figure 8 shows that the temperature profile does not reach its steady state

(a) t = 0.2(b) t = 1(c) t = 2

Figure 8. Effect of decaying parameter a on the distribution of w q with y for various values of t (m = 3, S = 1, τ_{D }= 0.05 ).

monotonically. Increasing a decreases θ for all y and t. It is observed also that the velocity component u reaches the steady state faster than w which, in turn, reaches the steady state faster than θ. This is expected as u is the source of w, while both u and w act as sources for the temperature.

Figures 9-11 show the variation of the velocity components u and w and the temperature θ at the central

(a)(b)(c)

Figure 9. Effect of the Hall current m on u at y = 0 for various values τ_{D} (S = 0, a=1).

plane of the channel (y = 0) with time for various values of the Hall parameter m and for t_{D }= 0.0, 0.05, and 0.1. In these figures S = 0. Figure 9 shows that u increases with increasing m for all values of t_{D} as the effective conductivity (s /(1+m^{2})) decreases with increasing m which reduces the magnetic damping force on u. It is observed also from the figure that the time at which u reaches its steady state value increases with increasing m. Increasing

(a)(b)(c)

Figure 10. Effect of the Hall current m on w at y = 0 for various values τ_{D} (S = 0, a = 1).

t_{D} increases u for all m and its effect on u becomes more pronounced for higher values of m. In Figure 10 the velocity component w increases with increasing m as w is a result of the Hall effect. On the other hand, at small times, w decreases when m increases. This happens due to the fact that, at small times w is very small and then the source term of w is proportional to (m u/(1+m^{2}) ) which decreases with increasing m (m > 1 ). This accounts for

the crossing of the curves of w with t for all values of t_{D}. Figure 10 indicates that the influence of t_{D} on w depends on t and becomes more clear when m is large. It is observed that, increasing t_{D} decreases w and increasing m increases w. Figure 11 shows that the influence of m on θ depends on t. Increasing m decreases θ at small times,

(a)(b)(c)

Figure 12. The variation of u profiles with t for various values of τ_{D} (m = 3, S = 1, a = 1).

but this is reversed at large times. This is due to the fact that, for small times, u and w are small and an increase in m increases u but decreases w. Then, the Joule dissipation which is also proportional to (1/1+m^{2}) decreases. For large times, increasing m increases both u and w and, in turn, increases the Joule and viscous dissipations. This accounts for the crossing of the curves of θ with time for

(a)(b)(c)

Figure 13. The variation of w profiles with t for various values of τ_{D} (m = 3, S = 1, a = 1).

all values of t_{D}. It is also observed that increasing t_{D} decreases the temperature θ for all values of m. This is because increasing t_{D} decreases both u and w and their gradients which decreases the Joule and viscous dissipations. The figure shows also that the time at which θ reaches its steady state value increases with increasing m while it is not greatly affected by changing t_{D}.

Figures (12-14) present the profiles of the velocity components u and w and the temperature θ for variousvalues of time t and for t_{D }= 0.0, 0.05 and 0.1. The fig-

(a)(b)(c)

Figure 14. The variation of q profiles with t for various values of τ_{D} (m = 3, S = 1, a = 1).

ures are evaluated for m = 3 and S = 1. It is clear from Figures 12 and 13 that the effect of yield stress t_{D} on u and w depends on t and y. Figure 12 shows that, for small t, increasing the yield stress t_{D} decreases u for small y, but this is reversed for large y. As time develops, increasing t_{D} increases u for all y. Figure 13 shows that increasing t_{D} increases w for all values of y, but increasing t_{D} more decreases w for large y.

For large t, increasing t_{D} decreases w for small t and all values of y. This can be attributed to the fact that increasing t_{D} will delay the attainment of maxima of u and w. It is also observed, from Figures 12 and 13, that the velocity components u and w do not reach their steady state monotonically. Figure 14 shows that the temperature profile does not reach its steady state monotonically. Increasing t_{D} increases θ for all y and t as a result of increasing the dissipations. It is observed also that the velocity component u reaches the steady state faster than w which, in turn, reaches the steady state faster than θ. This is expected as u is the source of w, while both u and w act as sources for the temperature.

5. Conclusions

A finite difference method is used to solve the transient Couette flow and heat transfer of a Casson non-Newtonian fluid under the influence of unsteady pressure gradient and uniform magnetic field. In the present work, we study Hall effect. The effects of the decaying parameter a, Casson yield stress t_{D}, and the Hall parameter m on the velocity and temperature distributions are studied. The decaying parameter a affects the main velocity components u and w and the temperature θ. The Hall term affects the main velocity component u in the x-direction and gives rise to another velocity component w in the z-direction. The results show that the influence of the parameters a and t_{D} on u and w depend on time and the Hall parameter m. It is also found that the effect of m on w and θ depends on time for all values of t_{D} which accounts for a crossover in the w-t and θ-t graphs for various values of m. The effect of m on the magnitude of θ depends on t_{D} and becomes more pronounced in case of small t_{D}. It is also found that the effect of a on the magnitude of θ depends on t_{D} and becomes more pronounced in case of smallt_{D.}

6. REFERENCES

- I. N. Tao, “Magnetohydrodynamic Effects on the Formation of Couette Flow,” Journal of Aerospace Science, Vol. 27, 1960, pp. 334-347.
- S. D. Nigam and S. N. Singh, “Heat Transfer by Laminar Flow Between Parallel Plates under the Action of Transverse Magnetic Field,” The Quarterly Journal of Mechanics and Applied Mathematics, Vol. 13, 1960, pp. 85-97. doi:10.1093/qjmam/13.1.85
- R. A. Alpher, “Heat Transfer in Magnetohydrodynamic Flow between Parallel Plates,” International Journal of Heat and Mass Transfer, Vol. 3, No. 2, 1961, pp. 108- 112. doi:10.1016/0017-9310(61)90073-4
- I. Tani, “Steady Motion of Conducting Fluids in Channels under Transverse Magnetic Fields with Consideration of Hall Effect,” Journal of Aerospace Science, Vol. 29, 1962, pp. 287-296.
- G. W. Sutton and A.Sherman, “Engineering Magnetohydrodynamics,” McGraw-Hill, New York, 1965.
- V. M. Soundalgekar, N. V. Vighnesam and H. S. Takhar, “Hall and Ion-Slip Effects in MHD Couette Flow with Heat Transfer,” IEEE Transactions on Plasma Science, Vol. PS-7, No. 3, 1979, pp. 178-182. doi:10.1109/TPS. 1979.4317226
- V. M. Soundalgekar and A. G. Uplekar, “Hall Effects in MHD Couette Flow with Heat Transfer,” IEEE Transactions on Plasma Science, Vol. PS-14, No. 5, 1986, pp. 579-583. doi:10.1109/TPS.1986.4316600
- E. M. H. Abo-El-Dahab, “Effect of Hall Currents on Some Magnetohydrodynamic Flow Problems,” Master Thesis, Helwan University, Ain Helwan, 1993.
- H. A. Attia and N. A. Kotb, “MHD Flow between Two Parallel Plates with Heat Transfer,” Acta Mechanica, Vol. 117, No. 1-4, 1996, pp. 215-220. doi:10.1007/BF011810 49
- H. A. Attia, “Hall Current Effects on the Velocity and Temperature Fields of an Unsteady Hartmann Flow,” Canadian Journal of Physics, Vol. 76, No. 9, 1998, pp. 739-746. doi:10.1139/cjp-76-9-739
- N. Casson, “A Flow Equation for Pigment Oil-Suspensions of the Printing Ink Type,” In: C. C. Mill, Ed., Rheolgy of Disperse Systems, Pergamon Press, London, 1959, p. 84.
- B. Tamamashi, “Consideration of Certain Hemorheological Phenomena from the Stand-Point of Surface Chemistry,” In: A. L. Copley, Ed., Hemorheology, Pergamon Press, London, 1968, p. 89.
- W. P. Walawander, T. Y. Chen and D. F. Cala, “An Approximate Casson Fluid Model for Tube Flow of Blood,” Biorheology, Vol. 12, No. 2, 1975, pp. 111-124.
- R. L. Batra and B. Jena, “Flow of a Casson Fluid in a Slightly Curved Tube,” International Journal of Engineering Science, Vol. 29, No. 10, 1991, pp. 1245-1258. doi:10.1016/0020-7225(91)90028-2
- B. Das and R. L. Batra, “Secondary Flow of a Casson Fluid in a Slightly Curved Tube,” International Journal of Non-Linear Mechanics, Vol. 28, No. 5, 1993, pp. 567- 577. doi:10.1016/0020-7462(93)90048-P
- M. E. Sayed-Ahmed and H. A. Attia, “Magnetohydrodynamic Flow and Heat Transfer of a Non-Newtonian Fluid in an Eccentric Annulus,” Canadian Journal of Physics, Vol. 76, No.5, 1998, pp. 391-401. doi:10.1139/cjp-76- 5-391
- H. A. Attia and M. E. Sayed-Ahmed, “Hydrodynamic Impulsively Lid-Driven Flow and Heat Transfer of a Casson Fluid,” Tamkang Journal of Science and Engineering, Vol. 9, No. 3, 2006, pp. 195-204.
- I. A. Abdalla, H. A. Attia and M. E. Sayed-Ahmed, “Adomian’s Polynomial Solution of Unsteady NonNewtonian MHD Flow,” Applied Mathematics and Information Science, Vol. 1, No. 3, 2007, pp. 227-245.
- H. A. Attia, “Unsteady MHD Couette Flow and Heat Transfer of Dusty Fluid with Variable Physical Properties,” Applied Mathematics and Computation, Vol. 177, No. 1, 2006, 308-318. doi:10.1016/j.amc.2005.11.010
- A. Mehmood and A. Ali, “The Effect of Slip Condition on Unsteady MHD Flow Oscillatory Flow of a Viscous Fluid in a Planner Channel,” Romanian Journal of Physics, Vol. 52, No. 1-2, 2007, pp. 85-91.
- E. Rukmangadachari, “Dusty Viscous Flow through a Cylinder of Rectangular Cross-Section under Time Dependent Pressure Gradient,” Defense Science Journal, Vol. 31, No. 2, 1981, pp. 143-153.
- B. J. Gireesha, C. S. Bagewadi and B. C. Prasanna Kamara, “Flow of Unsteady Dusty Fluid under Varying Pulastile Pressure Gradient in Anholonomic Coordinate System,” Electronic Journal of Theoretical Physics, Vol. 4, No. 14, 2007, pp. 9-16.
- M. Antia, “Numerical Methods for Scientists and Engineers,” Tata McGraw-Hill, New Delhi, 1991.