# Micropolar Fluid Flow and Heat Transfer over a Nonlinearly Stretching Plate with Viscous Dissipation

Authors: Kartini Ahmad; Anuar Ishak; Roslinda Nazar
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

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 parameterK, 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.


---

## Body

## 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 temperatureT∞ 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]).Figure 1
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 coefficientCf 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.Table 1
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.796873Figures2 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).Figure 2
Variation of the skin friction coefficientCfRex1/2 with K for various values of n.Figure 3
Variation of the local Nusselt numberNux/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 Figures4 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.Figure 4
Velocity profilesf′(η) for various values of K when n=0.1 and n=1.0.Figure 5
Temperature profilesθ(η) for various values of K when Pr=1, Ec=1, and n=0.1 and n=1.Figure 6
Velocity profilesf′(η) for n=0 and n=1.Figure 7
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 parameterK, 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.



---

*Source: 257161-2013-06-24.xml*