﻿ Heat and Mass Transfer in Visco-Elastic Fluid through Rotating Porous Channel with Hall Effect

Open Journal of Fluid Dynamics
Vol.06 No.01(2016), Article ID:64426,19 pages
10.4236/ojfd.2016.61002

Heat and Mass Transfer in Visco-Elastic Fluid through Rotating Porous Channel with Hall Effect

Pradip Kumar Gaur, Abhay Kumar Jha

Department of Mathematics, JECRC University, Jaipur, India    Received 28 March 2015; accepted 8 March 2016; published 11 March 2016

ABSTRACT

This paper examined the hydromagnetic boundary layer flow of viscoelastic fluid with heat and mass transfer in a vertical channel with rotation and Hall current. A constant suction and injection is applied to the plates. A strong magnetic field is applied in the direction normal to the plates. The entire system rotates with uniform angular velocity (W), about the axis perpendicular to the plates. The governing equations are solved by perturbation technique to obtain an analytical result for velocity, temperature, concentration distributions and results are presented graphically for various values of viscoelastic parameter (K2), Prandtl number (Pr), Schmidt number (Sc), radiation parameter (R), heat generation parameter (Qh) and Hall parameter (m).

Keywords:

Visco-Elastic Fluid, MHD, Hall Effect, Porous Medium 1. Introduction

2. Mathematical Formulation

The constitutive equations for the rheological equation of state for an elastico-viscous fluid (Walter’s liquid B') are . (1) . (2)

in which . (3)

N(t) is the distribution function of relaxation times. In the above equations pik is the stress tensor, p an arbitrary isotropic pressure, gik is the metric tensor of a fixed co-ordinate system xi and is the rate of strain tensor. It was shown by Walter’s  that Equation (2) can be put in the following generalized form which is valid for all types of motion and stress (4)

where is the position at time t' of the element that is instantaneously at the print xi at time “t”. The fluid with equation of state (1) and (4) has been designated as liquid B'. In the case of liquids with short memories, i.e. short relaxation times, the above equation of state can be written in the following simplified form . (5)

In which is the limiting viscosity at small rates of shear, and de-

notes the convected time derivative. We consider Oscillatory free convective flow of a viscous incompressible and electrically conducting fluid between two insulating infinite vertical permeable plates. A strong transverse magnetic field of uniform strength is applied along the axis of rotation by neglecting induced electric and magnetic fields. The fluid is assumed to be a gray, emitting, and absorbing, but non scattering medium. The radiative heat flux term can be simplified by using the Rosseland approximation. It is also assumed that the chemically reactive species undergo first-order irreversible chemical reaction.

The equations governing the flow of fluid together with Maxwell’s electromagnetic equations are as follows.

Equation of Continuity . (6)

Momentum Equation . (7)

Energy Equation . (8)

Concentration Equation . (9)

The generalized Ohm’s law, in the absence of the electric field  , is of the form.

(10)

where and pe are velocity, the electrical conductivity, the magnetic permeability, the cyclotron frequency, the electron collision time, the electric charge, the number density of the electron and the electron pressure, respectively. Under the usual assumption, the electron pressure (for a weakly ionized gas), the thermoelectric pressure, and ion slip are negligible, so we have from the Ohm’s law.

. (11)

From which we obtain that

. (12)

The solenoidal relation for the magnetic field where gives (con-

stant) everywhere in the flow, which gives. If are the component of electric current density, then the equation of conservation of electric charge gives.

Since the plates are infinite in extent, all the physical quantities except the pressure depend only on and. A cartesian coordinate system is assumed and -axis is taken normal to the plates, while -axis and -axis are in the upward and perpendicular directions on the plate (origin), respectively .The velocity components, , are in the -, -, -directions respectively. The governing equations in the rotating system in presence of Hall current, thermal radiation and chemical reaction are given by the following equations.

(13)

(14)

(15)

(16)

(17)

where is the Hall parameter, and are coefficients of thermal and solutal expansion, cp is the specific heat at constant pressure, is the density of the fluid, v is the kinematics viscosity k is the fluid thermal conductivity, is the acceleration of gravity, is the additional heat source, is the radiative heat flux, Dm is the molecular diffusivity, K1 is the chemical reaction rate constant .The radiative heat flux is

given by, in which and are the Stefan-Boltzmann constant and the mean absorp-

tion coefficient, respectively.

The initial and boundary conditions as suggested by the physics of the problem are

(18)

where e is a small constant.

We now introduce the dimensionless variables and parameter as follows:

. (19)

After combining (14) and (15) and taking, then (14)-(17) reduce to

(20)

where is the modified Grashof number, is the Prandtl number,

is the modified solutal Grashof number, is the Hartmann

number, is the radiation parameter, is the Schmidt number, and is the reaction parameter.

The boundary conditions (9) can be expressed in complex form as:

. (21)

3. Method of Solution

The set of partial differential Equations (20) cannot be solved in closed form. So it is solved analytically after these equations are reduced to a set of ordinary differential equations in dimensionless form. We assume that

(22)

where R stands for q or or and which is applicable for small perturbation.

Substituting (22) into (20) and comparing the harmonic and non harmonic terms, we obtain the following ordinary differential equations:

(23)

where and dashes denote the derivatives w.r.t.

The Transformed boundary conditions are

. (24)

The solutions of (23) under the boundary conditions are

(25)

(26)

(27)

(28)

(29)

(30)

(31)

(32)

. (33)

Equation (25) corresponds to the steady part, which gives 0 as the primary and V0 as secondary velocity components. The amplitude (resultant velocity) and phase difference due to these primary and secondary velocities for the steady flow are given by

(34)

where

Equation (26) and (27) together give the unsteady part of the flow. Thus unsteady primary and secondary velocity components and, respectively, for the fluctuating flow can be obtained from the following.

. (35)

The amplitude (resultant velocity) and the phase difference of the unsteady flow are given by

(36)

where

The amplitude (resultant velocity) and the phase difference are given by

(37)

where u = Real part of q and v = Imaginary part of q.

3.2. Amplitude and Phase Difference of Shear Stresses Due to Steady and Unsteady Flow at the Plate

The amplitude and phase difference of shear stresses at the stationary plate (η = 0), the steady flow can be obtained as

. (38)

For the unsteady part of flow, the amplitude and phase difference of shear stresses at the stationary plate () can be obtained as.

where

(39)

(40)

where

. (41)

The amplitude and phase difference of shear stresses at the stationary plate (η = 0) can be as

(42)

where Real part of and Where Imaginary part of.

The Nusselt number

. (43)

The rate of heat transfer (i.e. heat flux) at the plate in terms of amplitude and phase difference is given by

. (44)

The Sherwood number

. (45)

The rate of mass transfer (i.e. mass flux) at the plate in terms of amplitude and phase difference is given by.

. (46)

4. Results and Discussion

The system of ordinary differential Equation (23) with boundary condition (24) is solved analytically using perturbation technique. The solutions are obtained for the steady and unsteady velocity fields from (25)-(27), temperature fields from (28)-(30) and concentration fields are given by (31)-(33). The influences of each of the parameters on the thermal mass and hydrodynamic behaviors of buoyancy-induced flow in a rotating vertical channel are studied. The results are presented graphically. Temperature of the heated wall (left wall) is a function of time as given in the boundary conditions and the cooled wall at is maintained at a constant temperature. Further it is assumed that the temperature difference is small enough so that the density changes of the fluid in the system will be small. When the injection/suction parameter l is positive, fluid is injected through the hot wall into the channel and sucked out through the cold wall. The effect of various physical parameters on flow, heat, concentration fields, skin-friction Nusselt number, and Sherwood number are presented graphically in Figures 1-14. The profiles for resultant velocity Rn for the flow are in Figures 1-4 for suction/injection parameter l, rotation parameter W, viscoelastic parameter K2, and e respectively. Figure 1 shows that increase in the suction parameter l leads to an increase of Rn within the stationary plates. Similarly the

Figure 1. Resultant velocity Rn due to u and v versus l for different values of h at t = π/4.

Figure 2. Resultant velocity Rn due to u and v versus h for different values of W at t = π/4.

Figure 3. Resultant velocity Rn due to u and v versus h for different values of e at t = π/4.

Figure 4. Resultant velocity Rn due to u and v versus h for different values of K2 at t = π/4.

Figure 5. Phase difference a due to u and v versus h for different values of l at t = π/4.

Figure 6. Phase difference a due to u and v versus h for different values of W at t = π/4.

Figure 7. Phase difference a due to u and v versus h for different values of K2 at t = π/4.

Figure 8. Phase difference a due to u and v versus h for different values of Gr, Gm, M and m at π/4.

Figure 9. Concentration profiles against f for different values of x at t = π/4.

Figure 10. Concentration profiles against f for different values of Sc at t = π/4.

Figure 11. Temperature profiles q against h for different values of QH at t = π/4.

Figure 12. Temperature profiles q against h for different values of R at t = π/4

Figure 13. Nusselt number against for different values of QH with l = 0.5, Pr = 0.71, R = 2, e = 0.01 at t = π/4.

Figure 14. Sherwood number against w for different values of x with l = 0.5, e = 0.01, Sc = 0.3 at t = π/4.

resultant velocity increases with increasing values of rotation parameter W. This is due to the fact that the rotation effects being more dominant near the walls, so when W reaches high values secondary velocity component v decreases with increases in W as shown in Figure 2. From Figure 3, it is observed that the increase in the e leads to an increase of Rn within the stationary plates. From Figure 4, it inferred that resultant velocity Rn goes on increasing with increasing value of viscoelastic parameter K2.

The phase difference a for the flow is shown graphically in Figures 5-8. Figure 5 shows phase angle for various positive values of suction/injection parameter l. The figure shows that the phase angle a decreases with the increases of suction parameter. Figure 6 is the phase angle for various values of rotation parameter W. From this figure, it is observed that the phase angle a decreases with an increase in rotation parameter. From Figure 7, it is observed that the phase angle a increases with an increase in visco elastic parameter. Figure 8 shows the variation of a against h for different values of thermal Grashof number Gr, Solutal Grashof number Gm, and Hartmaan number. From this figure, it is found that the values of phase a decreases with increasing values of Gr, Gm and M while reversed effect is observed for the Hall parameter m.

The concentration profile f for the flow is shown graphically in Figure 9, Figure 10. From Figure 9, it is observed that with the increasing the value of the chemical reaction parameter decreases the concentration of species in the boundary layer; this due to the fact that destructive chemical reduces the solutal boundary thickness and increases the mass transfer. Opposite trend is seen in the case when Schmidt number is increased as noted in Figure 10. It may also be observed from this figure that the effect of Schmit number Sc is to be increases the concentration distribution in the solutal boundary layer.

The temperature profiles q are shown graphically in Figure 11 and Figure 12. Figure 11 has been plotted to depict the variation of temperature profiles against h for different values of heat absorption parameter QH by fixing other physical parameters. From this figure, we observe that temperature q decreases with increase in the heat absorption parameter QH because when heat is absorbed, the buoyancy force decreases the temperature profile. Figure 12 represents graph of temperature distribution with h for different values of radiation parameter. From this figure, we note that initial temperature q = 1 decreases zero satisfying boundary condition at h = 1.0 Further, it is observed from this figure that increase in the radiation parameter decreases the temperature distribution in the thermal boundary layer due to decreases in the thickness of the thermal boundary layer with thermal radiation parameter R. This is because large values of radiation parameter correspond to an increase in dominance of conduction over radiation, thereby decreasing the buoyancy force and temperature in the thermal boundary layer.

Figure 13 and Figure 14 show the amplitude of skin-friction, Nusselt number, and Sherwood number against frequency parameter w for different values of QH, and, respectively. The amplitude of Nusselt number decreases with increasing the value of heat source parameter QH which is shown in Figure 13. Figure 14 shows the variation of Sherwood number with and w. From this figure, it is observed that the Sherwood number decreases with increasing the value of chemical reaction parameter, and opposite trend is seen with increasing the values of w.

Cite this paper

Pradip KumarGaur,Abhay KumarJha, (2016) Heat and Mass Transfer in Visco-Elastic Fluid through Rotating Porous Channel with Hall Effect. Open Journal of Fluid Dynamics,06,11-29. doi: 10.4236/ojfd.2016.61002

References

1. 1. Cowling, T. (1957) Magnetohydrodynamics. Interscience Publications Inc., New York.

2. 2. Soundalgekar, V.M. (1979) Hal and Ion-Slip Effects in MHD Coutte Flow with Heat Transfer. IEEE Transactions on Plasma Science, PS-7, 178-186.
http://dx.doi.org/10.1109/TPS.1979.4317226

3. 3. Soundalgekar, V.M. and Uplekar, A.G. (1986) Hal Effects in MHD Coutte Flow with Heat Transfer. IEEE Transactions on Plasma Science, PS-7, 579-587.
http://dx.doi.org/10.1109/TPS.1986.4316600

4. 4. Hossain, M.A. and Rashid, R.I.M.I. (1987) Hall Effect on Hydromagnetic Free Convection Flow along a Porous Flat Plate with Mass Transfer. Journal of the Physical Society of Japan, 56, 97-104.
http://dx.doi.org/10.1143/JPSJ.56.97

5. 5. Attia, H.A. (1998) Hall Current Effects on Velocity and Temperature Fields of an Unsteady Hartmann Flow. Canadian Journal of Physics, 76, 73-78.
http://dx.doi.org/10.1139/cjp-76-9-739

6. 6. Singh, K.D. and Garg, B.P. (2010) Radiative Heat Transfer in MHD Oscillatory Flow through Porous Medium Bounded by Two Vertical Porous Plates. Bulletin of Calcutta Mathematical Society, 102, 129-139.

7. 7. Saha, L.K., Siddqia, S. and Hossain, M. A. (2011) Effect of Hall Current on MHD Natural Convection Flow from Vertical Permeable Flat Plate with Uniform Surface Heat Flux. Applied Mathematics and Mechanics, 32, 1127-1146.
http://dx.doi.org/10.1007/s10483-011-1487-9

8. 8. Aboeldahab, E.M. and Elbarby, E.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

9. 9. Abo-Eldahab, M. and ElAziz, M.A. (2005) Viscous Dissipation and Joule Heating Effects on MHD Free Convection from a Vertical Plate with Power-Law Variation in Surface Temperature in the Presence of Hall and Ion-Slip Currents. Apllied Mathematical Modeling, 29, 579-595.
http://dx.doi.org/10.1016/j.apm.2004.10.005

10. 10. Cogley, A.C., Vinceti, W.C. and Gilles, S.E. (1968) Differential Approximation for Radiation Transfer in a Nongray Gas near Equilibrium. American Institute of Aeronautics and Astronautics Journal, 6, 551-555.
http://dx.doi.org/10.2514/3.4538

11. 11. Mansour, M.A. (1990) Radiative and Free Convection Effects on the Oscillatory Flow past a Vertical Plate. Astrophysics and Space Science, 166, 269-275.
http://dx.doi.org/10.1007/BF01094898

12. 12. Hossain, M.A. and Thakar, H.S. (1996) Radiation Effect on Mixed Convection along a Vertical Plate with Uniform Surface Temperature. Heat and Mass Transfer, 314, 243-248.
http://dx.doi.org/10.1007/BF02328616

13. 13. Hossain, M.A., Alim, M.A. and Rees, S. (1999) The Effect of Radiation on Free Convection from a Porous Vertical Plate. International Journal of Heat and Mass Transfer, 42, 181-191.
http://dx.doi.org/10.1016/S0017-9310(98)00097-0

14. 14. Seddek, M.A. (2002) Effect of Radiation and Variable Viscosity on a MHD Free Convection Flow past a Semi Infinite Flat Plate with an Aligned Magnetic Field in the Case of Unsteady Flow. International Journal of Heat and Mass Transfer, 45, 931-935.
http://dx.doi.org/10.1016/S0017-9310(01)00189-2

15. 15. Muthucumarswamy, R. and Senthil, G.K. (2004) Studied the Effect of Heat and Mass Transfer on Moving Vertical Plate in the Presence of Thermal Radiation. Journal of Theoretical and Applied Mechanics, 31, 35-46.
http://dx.doi.org/10.2298/TAM0401035M

16. 16. Pal, D. (2009) Heat and Mass Transfer in Stagnation-Point Flow towards a Stretching Surface in the Presence of Buoyancy Force and Thermal Radiation. Meccanica, 44, 145-158.
http://dx.doi.org/10.1007/s11012-008-9155-1

17. 17. Aydin, A. and Kaya, A. (2008) Radiation Effect on MHD Mixed Convection Flow about a Permeable Vertical Plate. Heat and Mass Transfer, 45, 239-246.
http://dx.doi.org/10.1007/s00231-008-0428-y

18. 18. Mohamed, R.A. (2009) Double-Diffusive Convection Radiation Interaction on Unsteady MHD Flow over a Vertical Moving Porous Plate with Heat Generation and Soret Effects. Applied Mathematical Sciences, 3, 629-651.

19. 19. Chauhan, D.S. and Rastog, P. (2010) Radiation Effects on Natural Convection MHD Flow in a Rotating Vertical Porous Channel Partially Filled with a Porous Medium. Applied Mathematical Sciences, 4, 643-655.

20. 20. Ibrahim, S.Y. and Makinde, O.D. (2011) Radiation Effect on Chemically Reacting MHD Boundary Layer Flow of Heat and Mass Transfer through a Porous Vertical Flat Plate. International Journal of Physical Sciences, 6, 1508-1516.

21. 21. Pal, D. and Mondal, H. (2011) The Influence of Thermal Radiation on Hydromagnetic Darcy-Forchheimer Mixed Covection Flow past a Stretching Sheet Embedded in a Porous Medium. Mecccanica, 46, 739-753.
http://dx.doi.org/10.1007/s11012-010-9334-8

22. 22. Kandasamy, R., Periasamy, K. and Sivagnana Prabhu, K.K. (2005) Chemical Reaction, Heat and Mass Transfer on MHD Flow over a Vertical Stretching Surface with Heat Source and Thermal Stratification Effects. International Journal of Heat and Mass Transfer, 48, 4457-4561.
http://dx.doi.org/10.1016/j.ijheatmasstransfer.2005.05.006

23. 23. Muthucumaraswamy, R. and Janakiraman, B. (2008) Mass Transfer Effects on Isothermal Vertical Oscillating Plate in the Presence of Chemical Reaction. International Journal of Applied Mathematics and Mechanics, 4, 66-74.

24. 24. Sharma, P.R. and Singh, K.D. (2009) Unsteady MHD Free Convective Flow and Heat Transfer along a Vertical Porous Plate with Variable Suction and Internal Heat Generation. International Journal of Applied Mathematics and Mechanics, 4, 1-8.

25. 25. Sudheer Babu, M. and Satya Narayan, P.V. (2009) Effects of the Chemical Reaction and Radiation Absorption on Free Convection Flow through Porous Medium with Variable Suction in the Presence of Uniform Magnetic Field. Journal of Heat and Mass Transfer, 3, 219-234.

26. 26. Makinde, D. and Chinyoka, T. (2011) Numerical Study of Unsteady Hydromagnetic Generalized Coutte Flow of a Reactive Third-Grade Fluid with Asymmetric Convective Cooling. Computers and Mathematics with Applications, 61, 1167-1179.
http://dx.doi.org/10.1016/j.camwa.2010.12.066

27. 27. Pal, D. and Talukdar, B. (2011) Combined Effects of Joule Heating and Chemical Reaction on Unsteady Magnetohydrodynamic Mixed Convection of a Viscous Dissipating Fluid over a Vertical Plate in Porous Media with Thermal Radiation. Mathematical and Computer Modeling, 54, 3016-3036.
http://dx.doi.org/10.1016/j.mcm.2011.07.030

28. 28. Aldoss, T.K., Nimar, A.L., Jarrah, M.A. and Shaer, A.L. (1995) Magnetohydrodynamic Mixed Convection from a Vertical Plate Embedded in a Porous Medium. Numerical Heat Transfer, 28, 635-642.
http://dx.doi.org/10.1080/10407789508913766

29. 29. Rajgopal, K., Veena, P.H. and Parvin, V.K. (2006) Oscillatory Motion of an Electrically Conducting Viscoelastic Fluid over a Stretching Sheet in Saturated Porous Medium with Suction/Blowing. Mathematical Problems in Engineering, 1, 1-14.
http://dx.doi.org/10.1155/MPE/2006/60560

30. 30. Attia, H.A. (2004) Unsteady MHD Coutte Flow of a Viscoelastic Fluid Considering the Hall Effect. Canadian Journal of Physics, 82, 127-135.
http://dx.doi.org/10.1139/p03-117

31. 31. Chaudhary, R.C. and Jha, A.K. (2012) Heat and Mass Transfer in Elastic-Viscous Fluid past an Impulsively Started Infinite Vertical Plate with Hall Effect. Latin American Applied Research, 38, 17-26.

32. 32. Singh, K.D. (2012) Viscoelastic Mixed Convection MHD Oscillatory Flow through a Porous Medium Filled in a Vertical Channel. International Journal of Physics and Mathematical Sciences, 3, 194-205.

33. 33. Walters, K. (1962) The Solution of Flow Problems in the Case of Materials with Memories. Journal de Mécanique, 1, 474-478.

34. 34. Meyer, R.C. (1958) On Reducing Aerodynamics Heat Transfer Rates by Magnetohydrodynamic Techniques. Journal of the Aerospace Science, 25, 561-566.
http://dx.doi.org/10.2514/8.7781

Appendix

, ,

,

,

Nomenclature

Dimensional concentration

x-component of current density

Concentration at the left plate

Mean absorption coefficient

Specific heat at constant pressure

Chemical reaction rate constant

Specific heat at constant pressure

Hall parameter

Distance of the plate

Hartmann number

Chemical molecular diffusivity

Nusselt number

Electric charge

Number density of the electron

Acceleration due to gravity

Dimensional pressure

Modified Grashof number for mass transfer

Electron pressure

Modified Grashof number for heat transfer

Prandtl number

Magnetic field

Magnetic field of uniform strength

Dimensional heat source

x-Component of magnetic field

Heat source parameter

Current density

Temperature at the right wall

Resultant velocity

Dimensional time

Schmidt number

Non zero constant mean velocity

Sherwood number

Dimensional temperature

Electron velocity

Temperature at the left wall

Visco-elastic parameter

Dimensional injection /suction velocity

Primary velocity component for unsteady flow

Primary velocity component for steady flow

Secondary velocity component for steady flow

Secondary velocity component for unsteady flow

Velocity components are in the directions respectively.

Greek Symbols

Injection/suction parameter

Dynamic viscosity

Phase difference for the flow

Magnetic permeability

Phase difference of mass flux

Kinematic viscosity

Small positive constant

Oscillation parameter

Dimensionless distance

Dimensional parameter

Phase difference of heat flux

Angular velocity

Fluid thermal conductivity

Cyclotron frequency

Electric conductivity

Amplitude of mass flux

Stefan-Boltzmann constant

Amplitude of heat flux

Non dimensional concentration

Density

Non-dimensional temperature

Coefficient of thermal expansion

Coefficient of thermal expansion

Coefficient of solutal expansion

Phase difference of shear stresses for the steady flow

Phase difference of shear stresses for the unsteady flow

Phase difference of shear stresses for the flow

Amplitude of shear stresses for the flow

Amplitude of shear stresses for the steady flow

Amplitude of shear stresses for the unsteady flow