The flow and heat transfer of a micropolar fluid past a nonlinearly stretching plate is studied numerically, by taking into account the viscous dissipation effect. It is assumed that the plate is stretched nonlinearly from the slot where it is issued. The governing system of partial differential equations is transformed into ordinary differential equations, which are then solved numerically using a finite-difference scheme known as the Keller-box method. The effects of the governing parameters, namely, the material parameter K, the Eckert number Ec, the Prandtl number Pr, and the nonlinear stretching parameter n, on the flow field and the heat transfer characteristics are obtained and discussed. The velocity and the temperature profiles are also illustrated to aid the validity of the numerical results obtained. It is found that both the local Nusselt number and the magnitude of the skin friction coefficient increase with the nonlinear stretching parameter n, and the opposite trend occurs as K increases for fixed n.
1. Introduction
The initial study by Crane [1] has attracted many researchers to investigate similar problems on the boundary layer flow caused by a stretching sheet, as it has a number of applications in industry such as the extrusion of polymer sheet from a dye, continuous casting, drawing of plastic films, and crystal growing. The desired characteristics of the final product strictly depend on the rate of cooling and the process of stretching. This stretching sheet may not necessarily be linear, as it can also take nonlinear fashion, even though the problem might not have obvious technological relevance [2]. In view of this, Vajravelu [3] studied the flow over a nonlinearly stretching sheet, and Cortell [4, 5] examined the flow and heat transfer on a nonlinear stretching sheet for two different types of thermal boundary conditions on the sheet, namely, constant surface temperature (CST) and prescribed surface temperature (PST). Ganji et al. [6] reported the analytical solution of the magnetohydrodynamic flow over a nonlinearly stretching sheet. Similar problems have been studied recently by Ishak et al. [7], Prasad et al. [8], Van Gorder et al. [2], Raftari et al. [9], Abbas and Hayat [10], and Abel et al. [11], among others.
The Navier-Stokes equations alone do not adequately describe the flow properties for fluids with diversity in the physical structure which occur in most non-Newtonian fluids. Hence, the theory of micropolar fluid which was proposed by Eringen [12, 13] perhaps gives an accurate model for fluids where the polymeric and rotating particles need to be taken into account, by including the microrotational momentum equation in addition to the classical momentum equation. Studies of micropolar fluids have recently received considerable attention due to their applications in a number of processes that occur in industry. Such applications include suspension solutions, solidification of liquid crystals, animal bloods, and exotic lubricants. Some of the researchers who investigated such fluid with nonlinear stretching sheet are Hayat et al. [14, 15] and Rahman et al. [16], who investigated the micropolar fluid flow past a nonlinear stretching sheet taking into account the effects of a temperature-dependent viscosity and variable surface temperature. El-Aziz [17] examined the micropolar boundary-layer flow and heat transfer characteristics associated with a heated exponential stretching continuous sheet being cooled by a mixed convection flow, and Hsiao [18] has done an analysis for the heat and mass transfer with radiation effect of a micropolar fluid past a nonlinearly stretching sheet. Ishak et al. [19] studied the behavior of micropolar fluid flow over a fixed or continuous moving surface. The flow over a stretching/shrinking sheet immersed in a micropolar fluid has been considered by Ishak et al. [20], Yacob and Ishak [21, 22], and Rosali et al. [23], among others. They showed the existence of dual solutions for the shrinking case.
Motivated by the above investigations, the present paper studies the problem of micropolar fluid flow and heat transfer over a nonlinearly stretching plate with viscous dissipation. It is worth mentioning that viscous dissipation plays an important role in the polymer processing where it is considered like an energy source for producing heat which in turn delays the process of solidification and as a result more time is required to cool the final product (Hassan et al. [24]). It will affect the heat transfer problem as the source term does appear in the energy equation when the viscosity of the fluid or the velocity gradient is high.
2. Problem Formulation
Consider a steady two-dimensional laminar boundary layer flow over a nonlinearly stretching plate immersed in an incompressible micropolar fluid of ambient temperature T∞ as shown in Figure 1. It is assumed that the plate is stretched with velocity Uw=axn, and the plate temperature varies like Tw=T∞+bx2n, where x is the distance from the slot where the plate is issued and a, b, and n are constants. The boundary layer equations are [25, 26]
(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=μ+κρ∂2u∂y2+κρ∂N∂y,(3)ρj(u∂N∂x+v∂N∂y)=γ∂2N∂y2-κ(2N+∂u∂y),(4)u∂T∂x+v∂T∂y=kρcp∂2T∂y2+μ+κρcp(∂u∂y)2,
where u and v are the velocity components in the x and y directions, respectively. Further more, μ is the dynamic viscosity, κ is the vortex viscosity (or the microrotation viscosity), ρ is the fluid density, k is the thermal conductivity, cp is the specific heat at constant pressure, T is the fluid temperature, j is the microinertia density, N is the microrotation (or angular velocity), and γ is the spin gradient viscosity. The boundary conditions are
(5)u=Uw,v=0,N=-m∂u∂y,T=Twaty=0,u⟶0,N⟶0,T⟶T∞asy⟶∞,
where m is a constant with 0≤m≤1. The case m=1/2 indicates the vanishing of antisymmetric part of the stress tensor and denotes weak concentrations (Ahmadi [27]), which is considered in the present paper. Furthermore, we follow the work of many recent authors by assuming that γ=(μ+κ/2)j=μ(1+K/2)j, where K=κ/μ is the micropolar or material parameter. This assumption is invoked to allow the field of equations to predict the correct behavior in the limiting case when the microstructure effects become negligible and the total spin N reduces to the angular velocity (see Ahmadi [27] or Yücel [28]).
Physical model and coordinate system.
In order to solve (1)–(4) subject to the boundary conditions (5), we introduce the following similarity transformation (see Ziabakhsh et al. [26] and Ishak et al. [29]):
(6)η=(Uwνx)1/2y,ψ=(νxUw)1/2f(η),N=Uw(Uwνx)1/2h(η),h(η)=-12f′′(η),θ(η)=T-T∞Tw-T∞,
where η is the similarity variable, primes denote differentiation with respect to η, ν=μ/ρ is the kinematic viscosity, f is the dimensionless stream function, h is the dimensionless microrotation, and ψ is the stream function defined as u=∂ψ/∂y and v=-∂ψ/∂x which identically satisfies (1). Using transformation (6) and utilising the boundary condition N=-(1/2)(∂u/∂y) from (5), (2) and (3) reduce to the single equation
(7)(1+K2)f′′′+n+12ff′′-nf′2=0,
and the energy equation (4) becomes
(8)1Prθ′′+n+12fθ′-2nf′θ+Ec(1+K)f′′2=0,
where Pr is the Prandtl number and Ec is the Eckert number defined as
(9)Pr=μcpk,Ec=Uw2cp(Tw-T∞).
The transformed boundary conditions are
(10)f(0)=0,f′(0)=1,θ(0)=1,f′(η)⟶0,θ(η)⟶0asη⟶∞.
The physical quantities of interest are the skin friction coefficient Cf and the local Nusselt number Nux which are defined as
(11)Cf=τwρUw2,Nux=xqwk(Tw-T∞),
where the surface shear stress τw and the surface heat flux qw are given by
(12)τw=[(μ+κ)∂u∂y+κN]y=0,qw=-k(∂T∂y)y=0.
Using the dimensionless variables in (6), we obtain
(13)CfRex1/2=(1+K2)f′′(0),NuxRex1/2=-θ′(0),
where Rex=Uwx/ν is the local Reynolds number.
3. Results and Discussion
The transformed system of (7) and (8) with the boundary conditions (10) was solved numerically using the Keller-box method as described in [30] for some values of parameters. In order to validate the numerical results obtained, we compare our results with those obtained by Grubka and Bobba [31] as shown in Table 1, which shows a favorable agreement.
Values of θ′(0) for various values of Pr when K=0,n=1 and Ec=0.
Pr
Grubka and Bobba [31]
Present results
0.72
−1.0885
−1.088524
1.0
−1.3333
−1.333333
3.0
−2.5097
−2.509725
10.0
−4.7969
−4.796873
Figures 2 and 3 show the effect of the nonlinear stretching parameter n on the skin friction coefficient CfRex1/2 and the local Nusselt number Nux/Rex1/2, respectively, as the material parameter K varies when Pr=0.7 and Ec=1. As discussed by Alam et al. [32], Pr=0.7 corresponds to air that may be considered as a micropolar fluid when it is heavily contaminated with suspended particles of dirt. It is found that all the values of the skin friction coefficients CfRex1/2 are negative and decrease slowly as K increases from K=0 (Newtonian fluid) to K=2 (micropolar fluid) for all values of the nonlinearity parameter n. The same trend is observed as the plate is being stretched (n>0); that is, the values of the skin friction coefficient CfRex1/2 decrease with respect to K. It is worth mentioning that the plate is being as more stretched, it will decrease the values of the skin friction coefficients. But opposite trend is observed for the local Nusselt number as depicted in Figure 3. As the plate is being stretched, the values of the local Nusselt number Nux/Rex1/2 are greater and for a particular value of n, the local Nusselt number decreases slowly as K increases from 0 (Newtonian fluid) to positive values of K (micropolar fluid).
Variation of the skin friction coefficient CfRex1/2 with K for various values of n.
Variation of the local Nusselt number Nux/Rex1/2 with K for various values of n when Pr=0.7 and Ec=1.
The velocity and temperature distributions obtained by the finite-difference method for various values of the governing parameters are displayed in Figures 4 to 7. Figures 4 and 5 depict the effect of the micropolar/material parameter K when Pr=1, Ec=1, and n=0.1 and n=1.0 on the velocity and temperature profiles. On the other hand, Figures 6 and 7 display the velocity and temperature distributions for different values of the Eckert number Ec when Pr=0.7, K=1, and n=0 and n=1. It is seen that an increase in the Eckert number Ec increases the thermal boundary layer thickness, while stretching the plate decreases the thermal boundary layer. There are no changes in the velocity boundary layer thickness as illustrated in Figure 6. This is expected as the viscous dissipation only affects the thermal field. We can see that the velocity and the temperature profiles for all the figures subside monotonously to zero as η increases which validate the boundary conditions (10), asymptotically.
Velocity profiles f′(η) for various values of K when n=0.1 and n=1.0.
Temperature profiles θ(η) for various values of K when Pr=1, Ec=1, and n=0.1 and n=1.
Velocity profiles f′(η) for n=0 and n=1.
Temperature profiles θ(η) for various values of Ec when Pr=0.7, K=1, and n=0 and n=1.
4. Conclusions
We have done a numerical investigation on the effects of the material parameter K, the viscous dissipation Ec, the nonlinear stretching parameter n, and the Prandtl number Pr on the fluid flow and heat transfer characteristics toward a nonlinear stretching sheet immersed in a micropolar fluid. It is found that both the magnitude of the skin friction coefficient and the local Nusselt number increase with the nonlinear stretching parameter n and decrease when K increases for fixed values of n.
Acknowledgments
The authors wish to express their thanks to the anonymous referees for their valuable comments and suggestions. The financial support received from the Universiti Kebangsaan Malaysia (Project code: DIP-2012-31) is gratefully acknowledged.
CraneL. J.Flow past a stretching plate19702146456472-s2.0-3364609921210.1007/BF01587695Van GorderR. A.VajraveluK.AkyildizF. T.Existence and uniqueness results for a nonlinear differential equation arising in viscous flow over a nonlinearly stretching sheet201124223824210.1016/j.aml.2010.09.011MR2735149ZBL1229.34041VajraveluK.Viscous flow over a nonlinearly stretching sheet2001124328128810.1016/S0096-3003(00)00062-XMR1855769ZBL1162.76335CortellR.Viscous flow and heat transfer over a nonlinearly stretching sheet2007184286487310.1016/j.amc.2006.06.077MR2294952ZBL1112.76022CortellR.Heat and fluid flow due to non-linearly stretching surfaces2011217197564757210.1016/j.amc.2011.02.029MR2799770GanjiD. D.BararniaH.SoleimaniS.GhasemiE.A nalytical solution of the magneto-hydrodynamic flow over a nonlinear stretching sheet20092320-21254125562-s2.0-69849105914IshakA.NazarR.PopI.Unsteady mixed convection boundary layer flow due to a stretching vertical surface2006312165182MR2310278PrasadK. V.VajraveluK.DattiP. S.Mixed convection heat transfer over a non-linear stretching surface with variable fluid properties20104533203302-s2.0-7554908949810.1016/j.ijnonlinmec.2009.12.003RaftariB.Mohyud-DinS. T.YildirimA.Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method20115423423452-s2.0-7995592804110.1007/s11433-010-4180-1AbbasZ.HayatT.Stagnation slip flow and heat transfer over a nonlinear stretching sheet201127230231410.1002/num.20523MR2752873AbelM. S.KumarK. A.RavikumarR.MHD flow and heat transfer with effects of buoyancy, viscous and Joules dissipation over a nonlinear vertical stretching porous sheet with partial slip20113285291EringenA. C.Theory of micropolar fluids196616118MR0204005EringenA. C.Theory of thermomicrofluids19723824804962-s2.0-33847495482HayatT.AbbasZ.JavedT.Mixed convection flow of a micropolar fluid over a non-linearly stretching sheet200837256376472-s2.0-3804903632910.1016/j.physleta.2007.08.006HayatT.JavedT.AbbasZ.MHD flow of a micropolar fluid near a stagnation-point towards a non-linear stretching surface20091031514152610.1016/j.nonrwa.2008.01.019MR2502962RahmanM. M.RahmanM. A.SamadM. A.AlamM. S.Heat transfer in a micropolar fluid along a non-linear stretching sheet with a temperature-dependent viscosity and variable surface temperature2009305164916702-s2.0-7044953077210.1007/s10765-009-0656-5El-AzizM. A.Viscous dissipation effect on mixed convection flow of a micropolar fluid over an exponentially stretching sheet20098743593682-s2.0-6804912947410.1139/P09-047HsiaoK. L.Heat and mass transfer for micropolar flow with radiation effect past a nonlinearly stretching sheet20104644134192-s2.0-7795201060610.1007/s00231-010-0580-zIahakA.NazarR.PopI.Boundary-layer flow of a micropolar fluid on a continuous moving or fixed surface20068453994102-s2.0-3375123763710.1139/P06-059IshakA.LokY. Y.PopI.Stagnation-point flow over a shrinking sheet in a micropolar fluid201019711141714272-s2.0-7795168667710.1080/00986441003626169YacobN. A.IshakA.Micropolar fluid flow over a shrinking sheet20124722932992-s2.0-8485565207110.1007/s11012-011-9439-8YacobN. A.IshakA.Stagnation point flow towards a stretching/shrinking sheet in a micropolar fluid with a convective surface boundary condition20129036216262-s2.0-8485963560010.1002/cjce.20517RosaliH.IshakA.PopI.Micropolar fluid flow towards a stretching/shrinking sheet in a porous medium with suction201239826829HassanH.RegnierN.PujosC.DefayeG.The effect of viscous dissipation on the polymer temperature during injection moldingProceedings of the 5th European Thermal-Sciences Conference2008Eindhoven, The NetherlandsIshakA.NazarR.PopI.Post-stagnation-point boundary layer flow and mixed convection heat transfer over a vertical, linearly stretching sheet2008604303322MR2479870ZBL1175.76134ZiabakhshZ.DomairryG.BararniaH.Analytical solution of non-Newtonian micropolar fluid flow with uniform suction/blowing and heat generation20094044434512-s2.0-6734909038310.1016/j.jtice.2008.12.005AhmadiG.Self-similar solution of imcompressible micropolar boundary layer flow over a semi-infinite plate19761476396462-s2.0-0017138904YücelA.Mixed convection in micropolar fluid flow over a horizontal plate with surface mass transfer19892712159316022-s2.0-0024929901IshakA.NazarR.PopI.Heat transfer over a stretching surface with variable heat flux in micropolar fluids200837255595612-s2.0-3804916185110.1016/j.physleta.2007.08.003CebeciT.BradshawP.1988New York, NY, USASpringerxii+487Springer Study Editions10.1007/978-1-4612-3918-5MR969366GrubkaL. J.BobbaK. M.Heat transfer characteristics of a continuous, stretching surface with variable temperature198510712482502-s2.0-0022013926AlamM. S.SattarM. A.RahmanM. M.PostelnicuA.Local similarity solution of unsteady hydromagnetic convection flow of a micropolar fluid along a continuously moving permeable plate2010282951052-s2.0-79954534974