Heat Transfer to MHD Oscillatory Viscoelastic Flow in a Channel Filled with Porous Medium
The combined effect of a transverse magnetic field and radiative heat transfer on unsteady flow of a conducting optically thin viscoelastic fluid through a channel filled with saturated porous medium and nonuniform walls temperature has been discussed. It is assumed that the fluid has small electrical conductivity and the electromagnetic force produced is very small. Closed-form analytical solutions are constructed for the problem. The effects of the radiation and the magnetic field parameters on velocity profile and shear stress for different values of the viscoelastic parameter with the combination of the other flow parameters are illustrated graphically, and physical aspects of the problem are discussed.
The flow of an electrically conducting fluid has important applications in many branches of engineering science such as magnetohydrodynamics (MHD) generators, plasma studies, nuclear reactor, geothermal energy extraction, electromagnetic propulsion, and the boundary layer control in the field of aerodynamics. In the light of these applications, MHD flow in a channel has been studied by many authors; some of them are Nigam and Singh , Soundalgekar and Bhat , Vajravelu , and Attia and Kotb . A survey of MHD studies in the technological fields can be found in Moreau . The flow of fluids through porous media is an important topic because of the recovery of crude oil from the pores of the reservoir rocks; in this case, Darcy’s law represents the gross effect. Raptis et al.  have analysed the hydromagnetic free convection flow through a porous medium between two parallel plates. Aldoss et al.  have studied mixed convection flow from a vertical plate embedded in a porous medium in the presence of a magnetic field. Makinde and Mhone  have considered heat transfer to MHD oscillatory flow in a channel filled with porous medium.
In this study, an attempt has been made to extend the problem studied by Makinde and Mhone  to the case of viscoelastic fluid characterised by second-order fluid.
The constitutive equation for the incompressible second-order fluid is of the form where is the stress tensor, is the hydrostatic pressure, is the unit tensor, are the kinematic Rivlin-Ericksen tensors, are the material coefficients describing viscosity, elasticity, and cross-viscosity, respectively. The material coefficients have taken constants with and as positive and as negative (Markovitz and Coleman ). Equation (1) was derived by Coleman and Noll  from that of the simple fluids by assuming that stress is more sensitive to the recent deformation than to the deformation that occurred in the distant past.
2. Mathematical Formulation of the Problem
Consider the flow of a conducting optically thin fluid in a channel filled with saturated porous medium under the influence of an externally applied homogeneous magnetic field and radiative heat transfer as shown in Figure 1. It is assumed that the fluid has small electrical conductivity and the electromagnetic force produced is very small. The -axis is taken along the centre of the channel, and the -axis is taken normal to it. Then, assuming a Boussinesq incompressible fluid model, the equations governing the motion are given subject to boundary conditions where is the axial velocity, is the time, is the fluid temperature, is the pressure, is the gravitational force, is the radiative heat flux, is the co-efficient of volume expansion due to temperature, is the specific heat at constant pressure, is the thermal conductivity, is the porous medium permeability co-efficient, is the electromagnetic induction, is the magnetic permeability, is the intensity of the magnetic field, is the conductivity of the fluid, is the fluid density, and . It is assumed that both walls of temperature are high enough to induce radiative heat transfer. Following Cogley et al. , it is assumed that the fluid is optically thin with a relatively low density and the radiative heat flux is given by where is the mean radiation absorption co-efficient.
The following nondimensional quantities are introduced: where is the flow mean velocity.
The dimensionless governing equations together with the appropriate boundary conditions (neglecting the bars for clarity) can be written as with where and are Grashoff number, Hartmann number, Radiation parameter, Péclet number, Reynolds number, Darcy number, porous medium shape factor parameter, and viscoelastic parameter, respectively.
3. Method of Solution
The nondimensional shear stress at the wall is given by The rate of heat transfer across the channel’s wall is given as
4. Discussions and Conclusion
The purpose of this study is to bring out the effects of the viscoelastic parameter on the governing flow with the combination of the other flow parameters. The corresponding results for Newtonian fluid can be deduced from the above results by setting and it is worth mentioning here that these results coincide with that of Makinde and Mhone . We have considered the real parts of the results throughout for numerical validation. The velocity profile against is plotted in Figures 2–4 to observe the viscoelastic effects for various sets of values of Hartmann number and radiation parameter with fixed values of other flow parameters, namely, and . It is evident from Figures 2–4 that the velocity profile is parabolic in nature, and the values of velocity increase with the increasing values of the viscoelastic parameter in comparison with the Newtonian fluid. It is also noted from the figures that the behaviours of the velocity profiles remain the same with the increasing values of the viscoelastic parameter when (i) the values of the radiation parameter increase with the fixed values of the magnetic field parameter (Figures 2 and 3), (ii) the values of the magnetic field parameter increase with fixed values of the radiation parameter (Figures 3 and 4), and (iii) both the values of and increase (Figures 2 and 4).
Figures 5, 6, 7, and 8 exhibit the effects of the viscoelastic parameter on skin friction against magnetic field parameter and radiation parameter respectively with . Figures 5 and 6 show that for radiation parameter and , the shear stress decreases with the increasing values of for both Newtonian and non-Newtonian fluids, while shear stress increases for increasing values of in comparison with Newtonian fluid. From Figures 7 and 8, it is evident that the shear stress increases firstly and then decreases with the increasing values of for both Newtonian and non-Newtonian cases. Also, Figures 7 and 8 depict that the shear stress increases with the increasing values of the viscoelastic parameter in comparison with Newtonian fluid when the values of the magnetic field parameter increase.
It has also been observed that the temperature field is not significantly affected by the viscoelastic parameter.
V. M. Soundalgekar and J. P. Bhat, “Oscillatory MHD channel flow and heat transfer,” International Journal of Pure & Applied Mathematical, vol. 15, pp. 819–828, 1971.View at: Google Scholar
K. Vajravelu, “Exact periodic solution of a hydromagnetic flow in a horizontal channel,” Journal of Applied Mechanics, vol. 55, no. 4, pp. 981–983, 1988.View at: Google Scholar
H. A. Attia and N. A. Kotb, “MHD flow between two parallel plates with heat transfer,” Acta Mechanica, vol. 117, pp. 215–220, 1996.View at: Google Scholar
R. Moreau, Magnetohydrodynamics, Kluwar Academic, Dordrecht, The Netherlands, 1990.
A. Raptis, C. Massalas, and G. Tzivanidis, “Hydromagnetic free convection flow through a porous medium between two parallel plates,” Physics Letters A, vol. 90, no. 6, pp. 288–289, 1982.View at: Google Scholar
T. K. Aldoss, M. A. Al-Nimr, M. A. Jarrah, and B. Al-Shaer, “Magnetohydrodynamic mixed convection from a vertical plate embedded in a porous medium,” Numerical Heat Transfer A, vol. 28, no. 5, pp. 635–645, 1995.View at: Google Scholar
O. D. Makinde and P. Y. Mhone, “Heat transfer to MHD oscillatory flow in a channel filled with porous medium,” Romanian Journal of Physics, vol. 50, pp. 931–938, 2005.View at: Google Scholar
A. C. L. Cogley, W. G. Vinvent, and E. S. Giles, “Differential approximation for radiative heat transfer in non-linear equations-grey gas near equilibrium,” American Institute of Aeronautics and Astronautics, vol. 6, pp. 551–553, 1968.View at: Google Scholar