American Journal of Computational Mathematics
Vol.04 No.05(2014), Article ID:51291,9 pages
10.4236/ajcm.2014.45034

Hall Effects on Unsteady MHD Three Dimensional Flow through a Porous Medium in a Rotating Parallel Plate Channel with Effect of Inclined Magnetic Field

P. Sulochana

Intell Engineering College, Anantapuramu, Andhra Pradesh, India   Received 16 July 2014; revised 18 August 2014; accepted 3 September 2014

ABSTRACT

In this paper, we make an initial value investigation of the unsteady flow of incompressible viscous fluid between two rigid non-conducting rotating parallel plates bounded by a porous medium under the influence of a uniform magnetic field of strength H0 inclined at an angle of inclination α with normal to the boundaries taking hall current into account. The perturbations are created by a constant pressure gradient along the plates in addition to the non-torsional oscillations of the upper plate while the lower plate is at rest. The flow in the porous medium is governed by the Brinkman’s equations. The exact solution of the velocity in the porous medium consists of steady state and transient state. The time required for the transient state to decay is evaluated in detail and the ultimate quasi-steady state solution has been derived analytically. Its behaviour is computationally discussed with reference to the various governing parameters. The shear stresses on the boundaries are also obtained analytically and their behaviour is computationally discussed.

Keywords:

Hall Effects, Unsteady Rotating Flows, Three-Dimensional Flows, Parallel Plate Channels, Incompressible Viscous Fluids, Brinkman’s Model 1. Introduction

2. Formulation and Solution of the Problem

We consider the unsteady flow of an incompressible electrically conducting viscous fluid bounded by porous medium with two non-conducting rotating parallel plates. A uniform transverse magnetic field is applied to z-axis. In the presence of strong magnetic field a current is inclined in a direction normal to the both electric and magnetic field viz. Magnetic field of strength H0 inclined at angle of inclination to the normal to the boundaries in the transverse xz-plane. The inclined magnetic field gives rise to a secondary flow transverse to the channel. The hydro magnetic flow is generated in a fluid system by non-torsional oscillations of the upper plate. The lower plate is at rest. The origin is taken on the lower plate and the x-axis parallel to the direction of the upper plate. Since the plates are infinite in extent, all the physical quantities except the pressure depend on z and t only. In the equation of motion along x-direction, the x-component current density― and the z- component current density . We choose a Cartesian system 0 (x, y, z) such that the boundary walls are at z = 0 and z = l. The flow through porous medium governed by the Brinkman equations. The unsteady hydro magnetic equations governing flow through porous medium under the influence of a transverse magnetic field with reference to a rotating frame are (1) (2)

where, (u, w) is the velocity components along O (x, z) directions respectively. is the density of the fluid, is the magnetic permeability, is the coefficient of kinematic viscosity, k is the permeability of the medium, is the applied magnetic field. When the strength of the magnetic field is very large, the generalized Ohm’s law is modified to include the Hall current, so that (3)

where, q is the velocity vector, H is the magnetic field intensity vector, E is the electric field, J is the current density vector, is the cyclotron frequency, is the electron collision time, is the fluid conductivity and, is the magnetic permeability. In Equation (3) the electron pressure gradient, the ion-slip and thermo- electric effects are neglected. We also assume that the electric field E = 0 under assumptions reduces to (4) (5)

where is the Hall parameter.

On solving Equations (4) and (5) we obtain

(6)

(7)

Using the Equations (6) and (7), the equations of the motion with reference to rotating frame are given by

(8)

(9)

By combining the Equations (8) and (9), we get.

Let

(10)

The boundary and initial conditions are

(11)

(12)

We introduce the following non dimensional variables are

Using non-dimensional variables, the governing equations are (dropping asterisks)

(13)

where,

is the Hartmann number;

is the rotation parameter;

is inverse Darcy parameter and;

is the Hall parameter.

We choose is the prescribed of pressure gradient, then the Equation (13) reduces to

(14)

Corresponding initial and boundary conditions are

(15)

(16)

Taking Laplace transform of Equation (14) using initial condition (15) the governing equations in terms of the transformed variable reduces to

(17)

The relevant transformed boundary conditions are

(18)

(19)

Solving the Equation (17) and making use of the boundary conditions (18) and (19), we obtain

(20)

where

,

Taking inverse Laplace transform to the Equation (20), we obtain

(21)

The shear stresses on the upper plate and the lower plate are given by

(22)

3. Results and Discussion

The flow is governed by the non-dimensional parameters M the Hartman number, D−1 the inverse Darcy parameter, K is the rotation parameter and m is the Hall parameter. The velocity field in the porous region is evaluated analytically its behaviour with reference to variations in the governing parameters has been computationally analyzed. The profiles for u and w have been plotted in the entire flow field in the porous medium. The solution for the velocity consists of three kinds of terms 1) steady state, 2) the quasi-steady state terms associated with non-torsional oscillations in the boundary, 3) the transient term involving exponentially varying time dependence. From the expression (21), it follows that the transient component in the velocity in the fluid region

decays in dimensionless time. When the transient terms decay the steady oscillatory

solution in the fluid region is given by

(23)

(24)

We now discuss the quasi steady solution for the velocity for different sets of governing parameters namely viz. M the Hartman number and D−1 the inverse Darcy parameter, K the rotation parameter, m is the Hall parameter, P0 & P1 the non dimensional pressure gradients, the frequency oscillations ω, a and b the constants related to non torsional oscillations of the boundary, for computational analysis purpose we are fixing the axial pressure gradient as well as a and b, and, , ,. Figures 1-8 corresponding to the velocity components u and w along the prescribed pressure gradient for different sets of governing parameters when the upper boundary plate executes non-torsional oscillations. The magnitude of the velocity u and w increases for the sets of values 0.1 ≤ z ≤ 0.3 as well as which reduces for all values of z with increase in the intensity of the magnetic field (Figure 1 and Figure 5). The resultant velocity q decreases with increasing the Hartmann number M. The magnitude of the velocity u decreases in the upper part of the fluid region 0.1 ≤ z ≤ 0.2 while it experiences enhancement lower part 0.3 ≤ z ≤ 0.9 with increasing the inverse Darcy parameter D−1 (Figure 2). The magnitude of the velocity w increases in the upper part of the fluid region 0.1 ≤ z ≤ 0.3, while it reduces in lower part 0.4 ≤ z ≤ 0.9 with increasing the inverse Darcy parameter D−1 (Figure 6). The resultant velocity q reduces with increasing the inverse Darcy parameter D−1. The magnitude of velocity u decreases in the upper part of the fluid region while it experiences enhancement lower part 0.3 ≤ z ≤ 0.9 and also the magnitude of velocity w increases throughout the fluid region (Figure 3 and Figure 7). However the resultant velocity q enhances with increasing the Hall parameter m. Finally we notice that, from (Figure 4 and Figure 8) the magnitude of the velocity component enhances for 0.1 ≤ z ≤ 0.3 and z = 0.7, and reduces within the region 0.4 ≤ z ≤ 0.6 and 0.8 ≤ z ≤ 0.9 with increase in rotation parameter K. while the velocity component w enhances for 0.3 ≤ z ≤ 0.4 and z = 0.9, and reduces for 0.1 ≤ z ≤ 0.2, with increase in rotation parameter K.

The shear stresses and on the upper plate have been calculated for the different variations in the governing parameters and are tabulated in the Table 1, Table 2. On the upper plate we notice that the magnitudes of enhances the inverse Darcy parameter D−1, the hall parameter m, rotation parameter K decreases with increase in the Hartmann number M (Table 1). The magnitude of decreases with increase in the Hartmann number M, the inverse Darcy parameter D−1 rotation parameter K and the Hall parameter m fixing the other parameters (Table 2). The similar behaviour is observed on the lower plate (Table 3, Table 4). We also notice that the magnitude of the shear stresses on the lower plate is very small compare to its values of the upper plate.

Figure 1. The velocity profile for u with M.

Figure 2. The velocity profile for u with D1.

Figure 3. The velocity profile for u with m.

Figure 4. The velocity profile for u with K.

Figure 5. The velocity profile for w with M.

Figure 6. The velocity profile for w with D1.

Figure 7. The velocity profile for w with m.

Figure 8. The velocity profile for w with K.

Table 1. The shear stress on the upper plate.

Table 2. The shear stress on the upper plate.

Table 3. The shear stress on the lower plate.

Table 4. The shear stress on the lower plate.

4. Conclusions

1) The resultant velocity q enhances with increasing hall parameter m and rotation parameter K, and decreases with increasing inverse Darcy parameter D−1 as well as the Hartmann number M.

2) On the upper plate the magnitude of enhances when increasing the hall parameter m; rotation parameter K and the inverse Darcy parameter D−1 decrease with increase in the Hartmann number M.

3) On the upper plate the magnitude of shear stress enhances when increasing the hall parameter M; rotation parameter K and the inverse Darcy parameter D−1 decrease with increase in the Hartmann number M.

4) The similar behaviour is observed on the lower plate.

5) The magnitude of the shear stresses on the lower plate is very small than the values of the upper plate.

References

1. Raptis, A., Perdikis, C. and Tzivanidis, G. (1981) Free Convection Flow through a Porous Medium Bounded by a Vertical Surface. Journal of Physics D: Applied Physics, 14, 99-102. http://dx.doi.org/10.1088/0022-3727/14/7/001
2. Raptis, A., Tzivanidis, G. and Kafousias, N. (1981) Free Convection and Mass Transfer Flow through a Porous Medium Bounded by an Infinite Vertical Limiting Surface with Constant Suction. Letters Heat Mass Transfer, 8, 417-424. http://dx.doi.org/10.1016/0094-4548(81)90029-1
3. Raptis, A., Kafousias, N. and Massalas, C. (1982) 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, 62, 489-491. http://dx.doi.org/10.1002/zamm.19820620911
4. Singh, K.D. and Verma, G.N. (1995) Three-Dimensional Oscillatory Flow through a Porous Medium with Periodic Permeability. Zeitschrift für Angewandte Mathematik und Mechanik, 75, 599-604. http://dx.doi.org/10.1002/zamm.19950750811
5. Singh, K.D., Sharma, R. and Chand, K. (2000) Three-Dimensional Fluctuating Flow and Heat Transfer through a Porous Medium with Variable Permeability. Zeitschrift für Angewandte Mathematik und Mechanik, 80, 473-480. http://dx.doi.org/10.1002/1521-4001(200007)80:7<473::AID-ZAMM473>3.0.CO;2-1
6. Shercliff, T.A. (1965) A Text-Book of Magneto Hydro Dynamics. Pergamon Press, London.
7. Ferraro, V.C.A. and Plumpton, C. (1966) An Introduction to Magneto Fluid Mechanics. Clarandon Press, Oxford.
8. Crammer, K.P. and Pai, S.I. (1973) Magneto Fluid Dynamics for Engineers and Applied Physicist. McGraw Hill Book Co., New York.
9. Yan, J.T. and Chang, C.C. (1964) Magneto Hydrodynamic Couette Flow as Affected by Wall Electrical Conductance. Zeitschrift für Angewandte Mathematik und Physik, 15, 400-407. http://dx.doi.org/10.1007/BF01601291
10. Chang, C.C. and Lundgren, T.S. (1961) Duct Flow in Magneto Hydrodynamics. Zeitschrift für Angewandte Mathematik und Physik, 12, 100-114. http://dx.doi.org/10.1007/BF01601011
11. Attia, H.A. and Kotb, N.A. (1996) MHD Flow between Two Parallel Plates with Heat Transfer. Acta Mechanica, 117, 215-220. http://dx.doi.org/10.1007/BF01181049
12. Singh, K.D. and Alphonsa, M. (2008) Injection/Suction Effects on an Oscillatory Hydromagnetic Flow in a Rotating Horizontal Porous Channel. Indian Journal of Physics, 82, 435-445.
13. Attia, H.A. (2006) Time Varying Hydro Magnetic Couette Flow with Heat Transfer of a Dusty Fluid in the Presence of Uniform Suction and Injection Considering the Hall Effect. Turkish Journal of Engineering and Environmental Sciences, 30, 285-297.
14. Singh, K.D. and Rakesh, S. (2001) MHD Three-Dimensional Couette Flow with Transpiration Cooling. Zeitschrift für Angewandte Mathematik und Mechanik, 81, 715-720. http://dx.doi.org/10.1002/1521-4001(200110)81:10<715::AID-ZAMM715>3.0.CO;2-A
15. Singh, K.D. (2004) Influence of Moving Magnetic Field on Three-Dimensional Couette Flow. Zeitschrift für Angewandte Mathematik und Physik, 55, 894-902. http://dx.doi.org/10.1007/s00033-004-2065-8
16. Mazumder, B.S. (1991) An Exact Solution of Oscillatory Couette Flow in a Rotating System. Journal of Applied Mechanics, 58, 1104-1107. http://dx.doi.org/10.1115/1.2897694
17. Ganapathy, R.A. (1994) A Note on Oscillatory Couette Flow in a Rotating System. Journal of Applied Mechanics, 61, 208-209. http://dx.doi.org/10.1115/1.2901403
18. Mazumder, B.S., Gupta, A.S. and Datta, N. (1976) Hall Effects on Combined Free and Forced Convective Hydromagnetic Flow through a Channel. International Journal of Engineering Science, 14, 285-292. http://dx.doi.org/10.1016/0020-7225(76)90045-8
19. Singh, K.D. (2000) An Oscillatory Hydromagnetic Couette Flow in a Rotating System. Zeitschrift für Angewandte Mathematik und Mechanik, 80, 429-432. http://dx.doi.org/10.1002/1521-4001(200006)80:6<429::AID-ZAMM429>3.0.CO;2-1
20. Hartmann, J. and Lazarus, F. (1937) Kongelige danske videnskabernes selskab. Matematisk-Fysiske Meddelelser, 15, 6-7.
21. 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
22. Prasada Rao, D.R.V. and Krishna, D.V. (1981) Hall Effect on Unsteady Hydromagnetic Flow. Indian Journal of Pure and Applied Mathematics, 12, 270-276.
23. Prasada Rao, D.R.V., Krishna, D.V. and Debnath. L. (1982) Hall Effects on Free and Forced Convective Flow in a Rotating Channel. Acta Mechanica, 43, 49-59. http://dx.doi.org/10.1007/BF01175815
24. Veera Krishna, M. and Suneetha, S.V. (2009) Hall Effects on Unsteady MHD Rotating Flow of an Incompressible Viscous Fluid through a Porous Medium. Journal of Pure and Applied Physics, 21, 143-156.
25. Suneetha, S.V., Veera Krishna, M. and Siva Pradad, R. (2010) Hall Effect on Unsteady Rotating Hydro Dynamic Flow of an Incompressible Second Grade Fluid in a Porous Half Space. Journal of Pure and Applied Physics, 22, 143-156.
26. Das, S., Mandal, H.K. and Jana, R.N. (2013) Hall Effect on Unsteady MHD Flow through a Porous Channel in a Rotating System with Variable Pressure Gradient. Ph.D. Thesis, Vidya Sagar University, Midnapur, 1-23.