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.
1. Introduction
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 [1], Soundalgekar and Bhat [2], Vajravelu [3], and Attia and Kotb [4]. A survey of MHD studies in the technological fields can be found in Moreau [5]. 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. [6] have analysed the hydromagnetic free convection flow through a porous medium between two parallel plates. Aldoss et al. [7] have studied mixed convection flow from a vertical plate embedded in a porous medium in the presence of a magnetic field. Makinde and Mhone [8] 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 [8] to the case of viscoelastic fluid characterised by second-order fluid.
The constitutive equation for the incompressible second-order fluid is of the form
σ=-pI+μ1A1+μ2A2+μ3(A1)2,
where σis the stress tensor, pis the hydrostatic pressure, Iis the unit tensor, An(n=1,2)are the kinematic Rivlin-Ericksen tensors, μ1,μ2,andμ3are the material coefficients describing viscosity, elasticity, and cross-viscosity, respectively. The material coefficients μ1,μ2,andμ3have taken constants with μ1and μ3as positive and μ2as negative (Markovitz and Coleman [9]). Equation (1) was derived by Coleman and Noll [10] 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 x-axis is taken along the centre of the channel, and the y-axis is taken normal to it. Then, assuming a Boussinesq incompressible fluid model, the equations governing the motion are given
∂u∂t=-1ρ∂P∂x+υ1∂2u∂y2+υ2∂3u∂y2∂t-υ1uk-σeB02uρ+gβ(T-T0),∂T∂t=kρCp∂2T∂y2-1ρCp∂q∂y,
subject to boundary conditions
u=0,T=Twony=1,u=0,T=T0ony=0,
where u is the axial velocity, t is the time, T is the fluid temperature, P is the pressure, g is the gravitational force, q is the radiative heat flux, β is the co-efficient of volume expansion due to temperature, Cpis the specific heat at constant pressure, kis the thermal conductivity, Kis the porous medium permeability co-efficient, B0(=μeH0) is the electromagnetic induction, μeis the magnetic permeability, H0 is the intensity of the magnetic field, σe is the conductivity of the fluid, ρ is the fluid density, and υi=μi/ρ,(i=1,2). It is assumed that both walls of temperature T0,Tw are high enough to induce radiative heat transfer. Following Cogley et al. [11], it is assumed that the fluid is optically thin with a relatively low density and the radiative heat flux is given by
∂q∂y=4α2(T0-T),
where αis the mean radiation absorption co-efficient.
Geometry of the problem.
The following nondimensional quantities are introduced:
Re=Uaυ1,x̅=xa,y̅=ya,u̅=uU,θ=T-T0Tw-T0,H2=a2σeB02ρυ1,t̅=tUa,P̅=aPρυ1U,Da=Ka2,Gr=gβ(Tw-T0)a2υ1U,Pe=UaρCpk,N2=4α2a2k,
where U is the flow mean velocity.
The dimensionless governing equations together with the appropriate boundary conditions (neglecting the bars for clarity) can be written as
Re∂u∂t=-∂P∂x+∂2u∂y2-(s2+H2)u+GrT+γ∂3u∂y2∂t,Pe∂θ∂t=∂2θ∂y2+N2θ,
with
u=0,θ=1ony=1,u=0,θ=0ony=0,
where Gr,H,N,Pe,Re,Da,S(=1/Da), andγ=(υ2Re)/a2 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
In order to solve (7) and (8) for purely oscillatory flow, let
-∂P∂x=λeiωt,u(y,t)=u0(y)eiωt,θ(y,t)=θ0(y)eiωt,
where λis a constant and ω is the frequency of oscillation.
Substituting the above expressions into (7) and (8) and using (9), we get
(1+iγω)d2u0dy2-m22u0=-λ-Grθ0,d2θ0dy2+m12θ0=0,
subject to boundary conditions
u0=0,θ0=1ony=1,u0=0,θ0=0ony=0,
where m1=N2-iωPe and m2=S2+H2+iωRe.
Equations (11) and (12) are solved, and the solution, for the fluid velocity and temperature are given as follows:
u(y,t)={Grsin(m1y)(m12L+m22)sin(m1)M1e(m2y)/L+(M1-λm22)e-(m2y)/L+λm22+Grsin(m1y)(m12L+m22)sin(m1)}eiωt,θ(y,t)=sin(m1y)sinm1eiωt,
where M1=((λ/m22)e-m2/L-(λ/m22)-Gr)/(m12L+m22))/(em2/L-e-m2/L) and L=1+iωγ.
The nondimensional shear stress σat the wall y=0 is given by
σ=σ̅(μ1U/a)=[∂u∂y+γ∂2u∂y∂t]y=0.
The rate of heat transfer across the channel’s wall is given as
Nu=-∂θ∂y=-m1cos(m1)sin(m1)eiωt.
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 γ=0,and it is worth mentioning here that these results coincide with that of Makinde and Mhone [8]. We have considered the real parts of the results throughout for numerical validation. The velocity profile u against y is plotted in Figures 2–4 to observe the viscoelastic effects for various sets of values of Hartmann number Hand radiation parameter N(H=0.5,N=1.5;H=0.5,N=2.5;H=1.5,N=2.5) with fixed values of other flow parameters, namely, Pe=2,Re=2,s=1,t=0,Gr=2,λ=1,and ω=1. It is evident from Figures 2–4 that the velocity profile is parabolic in nature, and the values of velocity u increase with the increasing values of the viscoelastic parameter |γ|(γ=0,-0.10,-0.20) 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 Nincrease with the fixed values of the magnetic field parameter H(Figures 2 and 3), (ii) the values of the magnetic field parameter Hincrease with fixed values of the radiation parameter N(Figures 3 and 4), and (iii) both the values of Hand Nincrease (Figures 2 and 4).
Variation of u against y when H = 0.5, N = 1.5.
Variation of u against y when H = 0.5, N = 2.5.
Variation of u against y when H = 1.5, N = 2.5.
Figures 5, 6, 7, and 8 exhibit the effects of the viscoelastic parameter |γ|on skin friction σ against magnetic field parameter Hand radiation parameter N,respectively with Pe=2,Re=2,s=1,t=0,Gr=2,λ=1,andω=1. Figures 5 and 6 show that for radiation parameter N=1.5and N=2.5, the shear stress decreases with the increasing values of Hfor both Newtonian and non-Newtonian fluids, while shear stress increases for increasing values of |γ|(γ=0,-0.10,-0.20) 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 Nfor 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 Hincrease.
Variation of σ against H when N = 1.5.
Variation of σ against H when N = 2.5.
Variation of s against N when H = 1.5.
Variation of s against N when H = 2.5.
It has also been observed that the temperature field is not significantly affected by the viscoelastic parameter.
NigamS. D.SinghS. N.Heat transfer by laminar flow between parallel plates under the action of transverse magnetic field196013185972-s2.0-000246001210.1093/qjmam/13.1.85SoundalgekarV. M.BhatJ. P.Oscillatory MHD channel flow and heat transfer197115819828VajraveluK.Exact periodic solution of a hydromagnetic flow in a horizontal channel19885549819832-s2.0-0024144237AttiaH. A.KotbN. A.MHD flow between two parallel plates with heat transfer19961172152202-s2.0-0029748464MoreauR.1990Dordrecht, The NetherlandsKluwar AcademicRaptisA.MassalasC.TzivanidisG.Hydromagnetic free convection flow through a porous medium between two parallel plates19829062882892-s2.0-49049136188AldossT. K.Al-NimrM. A.JarrahM. A.Al-ShaerB.Magnetohydrodynamic mixed convection from a vertical plate embedded in a porous medium19952856356452-s2.0-0029403492MakindeO. D.MhoneP. Y.Heat transfer to MHD oscillatory flow in a channel filled with porous medium200550931938MarkovitzH.ColemanB. D.Incompressible second-order fluids19648691012-s2.0-000297502610.1016/S0065-2156(08)70353-3ColemanB. D.NollW.An approximation theorem for functionals, with applications in continuum mechanics1960613553702-s2.0-3425096306110.1007/BF00276168CogleyA. C. L.VinventW. G.GilesE. S.Differential approximation for radiative heat transfer in non-linear equations-grey gas near equilibrium19686551553