# MHD Flow and Heat Transfer in a Channel Bounded by a Shrinking Sheet and a Porous Medium Bed: Homotopy Analysis Method

Authors: Dileep Singh Chauhan; Rashmi Agrawal
Journal: ISRN Thermodynamics
 (2013)

---

## Abstract

MHD flow of viscous conducting fluid is considered between a shrinking sheet and a porous medium bed. Suction is applied at the upper shrinking sheet and its surface temperature is always maintained higher than the temperature of the lower porous bed surface. Similarity transformations and HAM are used to solve the governing equations for velocity and temperature fields. The effects of various pertinent parameters on the results are discussed graphically.


---

## Body

## 1. Introduction

The study of boundary layer viscous flow and heat transfer due to stretching/shrinking surfaces is very significant because of its several applications in engineering and industrial processes, such as extrusion of polymer sheets from a die, drawing of plastic films, polyester thin wall heat shrink tubing, shrink film, wire drawing, glass fiber, and paper production. It is very important to control drag and the rate of heat transfer from the sheet for a better quality of the final product. Crane [1] first examined such flow over a linearly stretching sheet and obtained an exact similarity solution for the steady boundary layer flow of a viscous fluid in a closed form. His work was extended under various situations and physical conditions, by several researchers, such as Rajagopal et al. [2], Dandapat and Gupta [3], Chiam [4], Cortell [5], Liu [6], Ariel et al. [7, 8], Abbas et al. [9], Abd El-Aziz [10], Chauhan and Rastogi [11], Chauhan and Agrawal [12], and Chauhan and Olkha [13].Recently the study of boundary layer flow and heat transfer adjacent to a shrinking sheet is given attention because of its several industrial applications, particularly in packaging industry. Flow due to a shrinking sheet is relatively a new phenomenon, where the velocity on the sheet is directed towards a fixed point. Steady flow in such a case does not occur normally, since the vorticity generated due to the shrinking of the sheet is not confined in the boundary layer unless a proper suction is applied at the sheet. Wang [14] reported the flow development adjacent to a shrinking sheet. Miklavčič and Wang [15] studied the existence and uniqueness of steady flow due to shrinking sheet and reported that for specific suction values dual solutions exist. Wang [16] studied stagnation-point flow over a shrinking sheet. Nadeem and Awais [17] investigated fluid flow through porous medium due to unsteady shrinking sheet with variable viscosity. Fang et al. [18] obtained a new solution for the Blasius equation of a shrinking sheet problem. Muhaimin et al. [19] investigated effects of heat and mass transfer on MHD flow over a shrinking sheet with suction. Rahimpour et al. [20] obtained the analytical solution of axisymmetric stagnation flow towards a shrinking sheet. Using HAM, Hayat et al. [21] investigated MHD rotating flow over a shrinking surface. Sajid and Hayat [22] examined MHD viscous flow due to a shrinking sheet using homotopy analysis method. The boundary layer flow adjacent to a shrinking sheet with power-law surface velocity is investigated by Fang [23] with mass transfer. A closed form analytical solution was presented by Fang and Zhang [24] for MHD flow over a shrinking sheet. Fang et al. [25] investigated the unsteady flow over a shrinking sheet with mass suction. Yao and Chen [26] presented the new analytical solution of Blasius equation for a shrinking sheet problem. Mohd and Hashim [27] studied MHD flow and heat transfer over a shrinking sheet embedded in a porous medium with suction. Cortell [28] examined MHD flow due to shrinking sheet for axisymmetric and two-dimensional cases with suction. Noor et al. [29] obtained a series solution for MHD flow due to shrinking sheet by ADM (Adomian decomposition method). Hayat et al. [30] developed an analytical solution for a rotating flow due to a shrinking sheet. Fang and Zhang [31] presented an exact analytical solution for thermal boundary layer over a shrinking sheet with mass transfer effects. Fang et al. [32] investigated analytically flow adjacent to a shrinking-sheet with a second-order slip flow model. Ishak et al. [33] investigated stagnation-point flow of a micropolar fluid over a shrinking sheet. The channel flows with stretching/shrinking walls are also very important problems in fluid dynamics. Brady and Acrivos [34], Borkakoti and Bharali [35], Chauhan and Jakhar [36], Hayat et al. [37], and Chauhan and Agrawal [38] investigated flow and heat transfer problem in channels with one wall as stretching or shrinking sheet.The purpose of the present investigation is to discuss the analytically solutions of electrically conducting viscous fluid flow and heat transfer in a channel bounded by a porous shrinking sheet and a porous medium bed of very small permeability, in the presence of a magnetic field. The equations governing velocity and temperature fields are simplified by a suitable similarity transformation, and then the reduced nonliner boundary value problem is solved analytically using homotopy analysis method. The physical effects of pertinent parameters on the flow and heat transfer are depicted graphically and discussed. The convergence of the solutions has also discussed by the so-called h-curves.


## 2. Formulation of the Problem

An electrically conducting, incompressible, viscous fluid flow in a horizontal channel bounded by a shrinking sheet and a porous medium bed is considered. A constant suctionw0 is applied at the upper shrinking sheet of the channel. The lower porous bed of the channel is of very small permeability k0. We model flow in porous bed by Darcy’s law. So, the flow in the porous bed is assumed to be zero in the absence of any external pressure gradient. A slip boundary condition suggested by Saffman [39] is applied at the fluid-porous medium interface. The permeability of the porous bed affects the flow in the channel through this boundary condition. Thus the effect of the porous bed is to introduce slip at the lower boundary of the horizontal channel. A constant magnetic field of strength B0 is applied to the normal direction of the horizontal channel. Magnetic Reynolds number is taken to be small; therefore, induced magnetic field is neglected. No external electric field is applied and the polarization effect of the ionized fluid is negligible. So the electric field E is assumed to be zero. A Cartesian coordinate system (x, y, z) is assumed with origin at the porous bed surface. The axis of x and axis of y are taken on the surface of the porous bed and the axis of z is normal to it (see schematic diagram, Figure 1). (u,v,w) are considered velocity components in x,y,and z directions, respectively, and T is the temperature of the fluid in the channel. The width of the channel is assumed h, so that walls of the channel are at z=0 and z=h.Figure 1
Schematic diagram.The governing equations for the present problem of MHD flow and heat transfer are given by(1)∂u∂x+∂v∂y+∂w∂z=0,(2)u∂u∂x+v∂u∂y+w∂u∂z=-1ρ∂P∂x+υ(∂2u∂x2+∂2u∂y2+∂2u∂z2)-σB02ρu,(3)u∂v∂x+v∂v∂y+w∂v∂z=-1ρ∂P∂y+υ(∂2v∂x2+∂2v∂y2+∂2v∂z2)-σB02ρv,(4)u∂w∂x+v∂w∂y+w∂w∂z=-1ρ∂P∂z+υ(∂2w∂x2+∂2w∂y2+∂2w∂z2),(5)ρCp(u∂T∂x+v∂T∂y+w∂T∂z)=k(∂2T∂x2+∂2T∂y2+∂2T∂z2).The corresponding boundary conditions are given by(6)atz=h,u=-ax,v=-a(c-1)y,w=w0,T=T1,atz=0,τzx=μαk0u,τzy=μαk0v,w=0,T=T0,
where P is the pressure; ρ is the density of the fluid; υ is the kinematic viscosity; σ is the electrical conductivity; Cp is the specific heat at constant pressure; k is the thermal conductivity in the clear fluid medium; T0 is the porous medium wall temperature; T1 is the temperature of upper stretching wall; k0 is the permeability of the porous medium; μ is the viscosity of the clear fluid; α is a nondimensional constant that depends only on the structure of the porous matrix.Here,a>0 is the shrinking constant. c is a constant, where c=1 indicates one-direction shrinking and c=2 corresponds to axisymmetric shrinking.


## 3. Method of Solution

We introduce the following similarity transformations:(7)u=axf′(η),v=a(c-1)yf′(η),w=-ahcf(η),η=zh,θ(η)=T-T0T1-T0,
where T1-T0=Dx,D being a nonzero constant.By using similarity transformations defined in (7), we find that the continuity equation (1) is automatically satisfied. Further, we see that (4) reduces to
(8)w∂w∂z=-1ρ∂P∂z+υ(∂2w∂z2).On integration, we get(9)w22=-Pρ+υdwdz+constant.Equations (2) and (3), using (9) and similarity transformations (7), are reduced to
(10)fiv-M2f′′-R((2-c)f′f′′-cff′′′)=0.Similarly, using (7), we see that (5) is reduced to
(11)θ′′=RPr(f′θ-cfθ′).The corresponding boundary conditions for the present problem are reduced to(12)atη=1,f′=-1,f=-w0ach=-λ,θ=1,atη=0,f′′=αKf′,f=0,θ=0,
where M=(σ/ρυ)B0h, the Hartmann number; R=ah2/υ, the nondimensional shrinking parameter; Pr=μCp/k, the Prandtl number; K=k0/h, the permeability parameter.


## 4. Homotopy Analysis Method

Following Liao [40], we now solve the nonlinear system consisting of (10)–(12) by HAM. For HAM solutions, we choose the initial guesses f0(η) and θ0(η) and the auxiliary linear operators Lf and Lθ in the following forms:
(13)f0(η)=2K(1-3λ)4K+αη+α(1-3λ)4K+αη2+α(2λ-1)-2K(1-λ)4K+αη3,θ0(η)=η,Lf(f)=d4fdη4,Lθ(θ)=d2θdη2.Here, auxiliary linear operators have the property that(14)Lf(f0)=Lf(a0,0+a0,1η+a0,2η2+a0,3η3)=0,Lθ(θ0)=Lθ(b0,0+b0,1η)=0.Letp∈[0,1] denote the embedding parameter and h1, h2 indicate the nonzero auxiliary parameters. Then we construct the following homotopies:
(15)Hf[f-(η;p);p]=(1-p)Lf[f-(η;p)-f0(η)]-pℏ1Nf[f-(η;p)],Hθ[θ-(η;p);p]=(1-p)Lθ[θ-(η;p)-θ0(η)]-pℏ2Nθ[θ-(η;p)],
where
(16)Nf[f-(η;p)]=∂4f-(η;p)∂η4-M2∂2f-(η;p)∂η2-R[(2-c)∂f-(η;p)∂η∂2f-(η;p)∂η2hhffff-cf-(η;p)∂3f-(η;p)∂η3],Nθ[θ-(η;p)]=∂2θ-(η;p)∂η2-RPr×[θ-(η;p)∂f-(η;p)∂η-cf-(η;p)∂θ-(η;p)∂η].Whenp=0, we have
(17)Hf[f-(η;0);0]=Lf[f-(η;0)-f0(η)],Hθ[θ-(η;0);0]=Lθ[θ-(η;0)-θ0(η)].Ifp=1, we have
(18)Hf[f-(η;1);1]=-ℏ1Nf[f-(η;1)],Hθ[θ-(η;1);1]=-ℏ2Nθ[θ-(η;1)].


### 4.1. Zeroth Order Deformation Equations

Thus, if we take(19)Hf[f-(η;p);p]=0,Hθ[θ-(η;p);p]=0,we obtain the family of equations as follows:(20)(1-p)Lf[f-(η;p)-f0(η)]=pℏ1Nf[f-(η;p)],(1-p)Lθ[θ-(η;p)-θ0(η)]=pℏ2Nθ[θ-(η;p),f-(η;p)].These equations are subjected to the boundary conditions:(21)f′-(1;p)=-1,f-(1;p)=-λ,θ-(1;p)=1,f′′-(0;p)=αKf′-(0;p),f-(0;p)=0,θ-(0;p)=0.Obviously, whenp=0, we see that
(22)f-(η;0)=f0(η),θ-(η;0)=θ0(η),
and when p=1, we have
(23)f-(η;1)=f(η),θ-(η;1)=θ(η).Thus, as the embedding parameterp increases from 0 to 1, the mappings f-(η;p) and θ-(η;p) vary smoothly from their known initial guesses, f0(η) and θ0(η), to the required solutions, f(η) and θ(η), respectively.Expandingf-(η;p) and θ-(η;p) by Taylor’s series in the embedding parameter p and using (22), we can write
(24)f-(η;p)=f0(η)+∑k=1∞fk(η)pk,
where fk(η)=[(1/k!)(∂kf-(η;p)/∂pk)]p=0 and
(25)θ-(η;p)=θ0(η)+∑k=1∞θk(η)pk,
where θk(η)=[(1/k!)(∂kθ-(η;p)/∂pk)]p=0.Now assume that for properly chosen values ofℏ1 and ℏ2 the series (24) and the series (25) are convergent at p=1. Thus we obtain
(26)f(η)=f0(η)+∑k=1∞fk(η),θ(η)=θ0(η)+∑k=1∞θk(η).


### 4.2.mth Order Deformation Equations

The governing equations offm(η) and θm(η) can be deduced from the zeroth order deformation equations (20) by differentiating these m times with respect to the embedding parameter p, setting p=0 and finally dividing them by |m_. Thus we have
(27)Lf[fm(η)-χmfm-1(η)]=ℏ1Rmf(η),Lθ[θm(η)-χmθm-1(η)]=ℏ2Rmθ(η),
subject to corresponding boundary conditions:
(28)fm′(1)=0,fm(1)=0,θm(1)=0,fm′′(0)=αKfm′(0),fm(0)=0,θm(0)=0,
where
(29)Rmf(η)=1|m-1_∂m-1Nf[f-(η;p)]∂pm-1|p=0=fm-1iv-M2fm-1′′+Re∑r=0m-1((2-c)fm-1-r′fr′′-cfm-1-rfr′′′),Rmθ(η)=1|m-1_∂m-1Nθ[θ-(η;p),f-(η;p)]∂pm-1|p=0=θm-1′′-RPr∑r=0m-1(fm-1-r′θr-cfm-1-rθr′).Here,χm is defined as
(30)χm={0,form≤1,1,form>1.Themth order (m≥1) deformation equations (27) are linear, and using the initial guesses f0(η), θ0(η), we can obtain f1(η),f2(η),f3(η), and so forth andθ1(η),θ2(η),θ3(η), and so forth using boundary conditions (28). With the help of the symbolic computation software MATHEMATICA, the solution series for (27) and (28) can be expressed in the following form:
(31)f(η)=∑m=0∞fm(η)=limM→∞[∑n=14M+3∑m=n-14M+2am,nηn],θ(η)=∑m=0∞θm(η)=limM→∞[∑n=14M+1∑m=n-14Mbm,nηn].The coefficientam,n and bm,n can be obtained using the following recurrence relations:
(32)am,0=0,am,1=χmam-1,1+2K4K+α∑k=44m+3Γm,k-4k(k-1)(k-2),am,2=χmam-1,2+α4K+α∑k=44m+3Γm,k-4k(k-1)(k-2),am,3=χmam-1,3-∑k=44m+3[{α+2K4K+α(k-3)+1}Γm,k-4k(k-1)(k-2)(k-3)],am,k=χmam-1,k+Γm,k-4k(k-1)(k-2)(k-3),4≤k≤4m+3,
where
(33)Γm,1k=ℏ1(a4m-1,k-(M2)a2m-1,khhh-R{(2-c)δ1m,k-cδ2m,k}),0≤k≤4m-1,δ1m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}a1r,ka2m-r-1,n-k,0≤n≤4m-1,δ2m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}ar,ka3m-r-1,n-k,0≤n≤4m-1,a4m,k=(k+1)a3m,k+1,0≤k≤4m-1,a3m,k=(k+1)a2m,k+1,0≤k≤4m,a2m,k=(k+1)a1m,k+1,0≤k≤4m+1,a1m,k=(k+1)am,k+1,0≤k≤4m+2,bm,0=0,bm,1=χmbm-1,1-∑k=24m+1γm,k-2k(k-1),bm,k=χmbm-1,k+γm,k-2k(k-1),2≤k≤4m+1,
where
(34)γm,1k=ℏ2(b2m-1,k-RPr{Δ1m,k-cΔ2m,k}),0≤k≤4m-1,Δ1m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}a1r,kbm-r-1,n-k,0≤n≤4m-1,Δ2m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}ar,kb1m-r-1,n-k,0≤n≤4m-1,b2m,k=(k+1)b1m,k+1,0≤k≤4m-1,b1m,k=(k+1)bm,k+1,0≤k≤4m.


## 4.1. Zeroth Order Deformation Equations

Thus, if we take(19)Hf[f-(η;p);p]=0,Hθ[θ-(η;p);p]=0,we obtain the family of equations as follows:(20)(1-p)Lf[f-(η;p)-f0(η)]=pℏ1Nf[f-(η;p)],(1-p)Lθ[θ-(η;p)-θ0(η)]=pℏ2Nθ[θ-(η;p),f-(η;p)].These equations are subjected to the boundary conditions:(21)f′-(1;p)=-1,f-(1;p)=-λ,θ-(1;p)=1,f′′-(0;p)=αKf′-(0;p),f-(0;p)=0,θ-(0;p)=0.Obviously, whenp=0, we see that
(22)f-(η;0)=f0(η),θ-(η;0)=θ0(η),
and when p=1, we have
(23)f-(η;1)=f(η),θ-(η;1)=θ(η).Thus, as the embedding parameterp increases from 0 to 1, the mappings f-(η;p) and θ-(η;p) vary smoothly from their known initial guesses, f0(η) and θ0(η), to the required solutions, f(η) and θ(η), respectively.Expandingf-(η;p) and θ-(η;p) by Taylor’s series in the embedding parameter p and using (22), we can write
(24)f-(η;p)=f0(η)+∑k=1∞fk(η)pk,
where fk(η)=[(1/k!)(∂kf-(η;p)/∂pk)]p=0 and
(25)θ-(η;p)=θ0(η)+∑k=1∞θk(η)pk,
where θk(η)=[(1/k!)(∂kθ-(η;p)/∂pk)]p=0.Now assume that for properly chosen values ofℏ1 and ℏ2 the series (24) and the series (25) are convergent at p=1. Thus we obtain
(26)f(η)=f0(η)+∑k=1∞fk(η),θ(η)=θ0(η)+∑k=1∞θk(η).


## 4.2.mth Order Deformation Equations

The governing equations offm(η) and θm(η) can be deduced from the zeroth order deformation equations (20) by differentiating these m times with respect to the embedding parameter p, setting p=0 and finally dividing them by |m_. Thus we have
(27)Lf[fm(η)-χmfm-1(η)]=ℏ1Rmf(η),Lθ[θm(η)-χmθm-1(η)]=ℏ2Rmθ(η),
subject to corresponding boundary conditions:
(28)fm′(1)=0,fm(1)=0,θm(1)=0,fm′′(0)=αKfm′(0),fm(0)=0,θm(0)=0,
where
(29)Rmf(η)=1|m-1_∂m-1Nf[f-(η;p)]∂pm-1|p=0=fm-1iv-M2fm-1′′+Re∑r=0m-1((2-c)fm-1-r′fr′′-cfm-1-rfr′′′),Rmθ(η)=1|m-1_∂m-1Nθ[θ-(η;p),f-(η;p)]∂pm-1|p=0=θm-1′′-RPr∑r=0m-1(fm-1-r′θr-cfm-1-rθr′).Here,χm is defined as
(30)χm={0,form≤1,1,form>1.Themth order (m≥1) deformation equations (27) are linear, and using the initial guesses f0(η), θ0(η), we can obtain f1(η),f2(η),f3(η), and so forth andθ1(η),θ2(η),θ3(η), and so forth using boundary conditions (28). With the help of the symbolic computation software MATHEMATICA, the solution series for (27) and (28) can be expressed in the following form:
(31)f(η)=∑m=0∞fm(η)=limM→∞[∑n=14M+3∑m=n-14M+2am,nηn],θ(η)=∑m=0∞θm(η)=limM→∞[∑n=14M+1∑m=n-14Mbm,nηn].The coefficientam,n and bm,n can be obtained using the following recurrence relations:
(32)am,0=0,am,1=χmam-1,1+2K4K+α∑k=44m+3Γm,k-4k(k-1)(k-2),am,2=χmam-1,2+α4K+α∑k=44m+3Γm,k-4k(k-1)(k-2),am,3=χmam-1,3-∑k=44m+3[{α+2K4K+α(k-3)+1}Γm,k-4k(k-1)(k-2)(k-3)],am,k=χmam-1,k+Γm,k-4k(k-1)(k-2)(k-3),4≤k≤4m+3,
where
(33)Γm,1k=ℏ1(a4m-1,k-(M2)a2m-1,khhh-R{(2-c)δ1m,k-cδ2m,k}),0≤k≤4m-1,δ1m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}a1r,ka2m-r-1,n-k,0≤n≤4m-1,δ2m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}ar,ka3m-r-1,n-k,0≤n≤4m-1,a4m,k=(k+1)a3m,k+1,0≤k≤4m-1,a3m,k=(k+1)a2m,k+1,0≤k≤4m,a2m,k=(k+1)a1m,k+1,0≤k≤4m+1,a1m,k=(k+1)am,k+1,0≤k≤4m+2,bm,0=0,bm,1=χmbm-1,1-∑k=24m+1γm,k-2k(k-1),bm,k=χmbm-1,k+γm,k-2k(k-1),2≤k≤4m+1,
where
(34)γm,1k=ℏ2(b2m-1,k-RPr{Δ1m,k-cΔ2m,k}),0≤k≤4m-1,Δ1m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}a1r,kbm-r-1,n-k,0≤n≤4m-1,Δ2m,n=∑r=0m-1∑k=max{0,n-4m+4r+1}min{n,4r+3}ar,kb1m-r-1,n-k,0≤n≤4m-1,b2m,k=(k+1)b1m,k+1,0≤k≤4m-1,b1m,k=(k+1)bm,k+1,0≤k≤4m.


## 5. Convergence of the HAM Solutions

Equations (31) are the series solutions of the present problem of MHD flow and heat transfer in a horizontal channel bounded by a shrinking sheet and a porous bed of very small permeability. As Liao [40] pointed out that the convergence region and convergence rate of the HAM solution series can be adjusted by means of the value of the auxiliary parameters ℏ1and ℏ2 here we have takenℏ1=ℏ2=ℏ. In order to find the convergence region for the series (31), we have plotted the so-called h-curves. We have plotted the curves of f′′(0) versus ℏ and θ′(0) versus ℏin Figure 2, for 20th, 21st, and 22nd order approximation. It clearly indicates that the valid region of ℏ is about-1.8≤ℏ≤-0.4. In all numerical calculations for plotting graphs, the order of approximation is taken as m=22 andℏ=-1.1.Figure 2
h-curve forf′′(0) and θ′(0) for M=1, c=3, K=0.001, λ=0.5, α=0.1, Re=2, and Pr=2.


## 6. Discussion

In the present investigation, the effects of a porous bed and a shrinking wall on the MHD flow and heat transfer are examined in a horizontal channel. The governing nonlinear partial differential equations for flow and heat transfer are converted into ordinary differential equations using similarity transformations. Then series solutions are obtained for velocity and temperature fields using homotopy analysis method (HAM). The convergence region of the solution series by HAM is obtained by plotting the so-called h-curve graph in Figure2. To ensure the steady flow near the shrinking sheet by confining the generated vorticity in the boundary layer, the opposite forces such as the wall mass suction and the magnetic field are taken as sufficiently high by assigning suitable values of M or λ.Figures3, 4, 5, 6, and 7 illustrate the variations of axial velocity for various values of the pertinent parameters. Figure 3 shows that magnitude of the axial velocity increases near the porous bed in the channel with the increase of the permeability parameter (K), while it decreases near the upper shrinking sheet. Permeability is the measure of ease to flow and it enhances fluid flow near the porous bed. The effect of the fluid saturated porous medium bed is to introduce a slip at one wall of the channel, and therefore the flow velocity at this wall is no longer equal to zero. This slip depends on the parameters of structure and permeability of the porous matrix and affects the flow inside the channel through boundary condition suggested by Saffman [39]. In Figure 4, we observe that the magnitude of the axial velocity decreases with the increase of the magnetic parameter (M) near the upper shrinking sheet since Lorentz force opposes the flow and leads to the deceleration of fluid flow near upper wall, while it increases near porous bed. Figure 5 illustrates that with increase in the value of suction parameter (λ), magnitude of the axial flow increases in the channel. Figure 6 shows that the effect of the shrinking parameter (R) is to reduce the flow near the upper sheet while the reverse effect is observed near the lower porous bed. Figure 7 depicts that the axial flow near the upper wall is greater when the sheet shrinks axisymmetrically than that when it shrinks only in the x direction, while near the porous bed reverse effect is observed.Figure 3
f
′
(
η
) for M=1, c=3, α=0.1, λ=0.5, and R=2.Figure 4
f
′
(
η
) for c=3, K=0.001, λ=0.5, α=0.1, and R=2.Figure 5
f
′
(
η
) for M=1, c=3, K=0.001, α=0.1, and R=2.Figure 6
f
′
(
η
) for M=1, c=3, K=0.001, λ=0.5, and α=0.1.Figure 7
f
′
(
η
) for M=1, K=0.001, α=0.1, and R=2.Figures8, 9, 10, 11, 12, and 13 illustrate profiles of the nondimensional temperature against the coordinateη. The effect of Prandtl number (Pr) on the temperature profiles can be seen in Figure 8. Increase in the Pr values leads to the decrease of the temperature in the channel and the thickness of the thermal boundary layer becomes smaller. As Pr value increases, the thermal diffusivity decreases which reduces the energy transfer ability; hence the temperature decreases. Note that very small values of Pr correspond to liquid metals, Pr~1 corresponds to air and diatomic gases, and very large values of Pr correspond to highly viscous oils. At room temperature,Pr~7 corresponds to water. Figure 9 depicts temperature distribution in the channel for various values of the permeability parameter (K). It is found that the effect of the permeability of lower bed is to decrease the temperature throughout the channel. It is also observed in Figures 10, 11, and 12 that the temperature in the channel decreases with increase in the value of suction parameter (λ) or magnetic parameter (M) or shrinking parameter (R). Further, it is observed in Figure 13 that the temperature throughout the channel for the axisymmetric shrinking case (c=2) is less than that for the case of one-direction shrinking (c=1).Figure 8
θ
(
η
) for M=1, c=3, K=0.001, λ=0.5, α=0.1, and R=2.Figure 9
θ
(
η
) for M=1, c=3, λ=0.5, α=0.1, R=2, and Pr=2.Figure 10
θ
(
η
) for M=1, c=3, K=0.001, R=2, Pr=2, and α=0.1.Figure 11
θ
(
η
) for Pr=2, c=3, K=0.001, λ=0.5, α=0.1, and R=2.Figure 12
θ
(
η
) for M=1, c=3, λ=0.5, α=0.1, K=0.001, and Pr=2.Figure 13
θ
(
η
) for M=1, λ=0.5, K=0.001, R=2, Pr=2, and α=0.1.Many industrial processes, for example, in metallurgical and packaging units, involve the drawing and cooling of continuous strips, sheets, and so forth and there quality depends on the rate of cooling and shear stresses on them. Therefore the skin friction coefficient and the rate of heat transfer are the physical quantities of interest to engineers and these are plotted in Figures14 and 15 for various values of the pertinent parameters. It is found that the magnetic field increases the surface shear stress for axial flow, that is, the local skin friction at the shrinking sheetη=1, because an increase in the strength of magnetic field leads to a thinner boundary layer, thereby causing an increase in the velocity gradient f′′(1) at the sheet. Shrinking parameter (R) or permeability parameter (K) also leads to an increase in the skin friction at the shrinking sheet, whereas the effect of the suction parameter (λ) is to reduce skin friction at the shrinking sheet. It is also noticed that the skin friction at shrinking sheet for axisymmetric shrinking is greater than that for the case of one-direction shrinking. It is observed in Figure 15 that the temperature gradient at the shrinking sheet increases as Prandtl number (Pr) increases, since a higher Pr fluid is more viscous and has a thinner thermal boundary layer which results at the sheet in a higher rate of heat transfer. It is also noticed that the temperature gradient at the shrinking sheet increases with the increase in the value of the shrinking parameter (R), or permeability (K), or suction parameter (λ), or magnetic parameter (M). It is also observed that the temperature gradient at the shrinking sheet is greater for the axisymmetric shrinking case than to one-direction shrinking.Figure 14
-
f
′′
(
1
) versus R for c=3, K=0.001, λ=0.5, and α=0.1.Figure 15
θ
′
(
1
) versus R for M=1, c=3, K=0.001, λ=0.5, and α=0.1.


## 7. Conclusion

The present paper studies the flow and heat transfer in a channel bounded by a porous medium bed and a shrinking sheet in the presence of a magnetic field. Using similarity transformations and homotopy analysis method, analytical series solutions are obtained and presented graphically to discuss features of these solutions for various values of the pertinent parameters. Convergence of the HAM series solution is ensured by taking suitable value of the parameterh. It is observed that the axial flow in the channel adjacent to upper wall is greater when it shrinks axisymmetrically to that when it shrinks in one direction, while reverse effect is observed at the other wall. Further, the effect of suction parameter is to increase the axial flow in the channel. Temperature in the channel decreases as the value of Prandtl number increases, and the effect of permeability parameter of the lower porous bed is to decrease temperature throughout the channel. It is also found that the temperature in the channel for the axisymmetric shrinking case is less than that for the case of one-direction shrinking. It is seen that the physical parameters such as suction and magnetic field control the rate of cooling and skin friction at the shrinking surface.



---

*Source: 291270-2013-01-17.xml*