# Effects of Thermal Radiation on Mixed Convection Flow of a Micropolar Fluid from an Unsteady Stretching Surface with Viscous Dissipation and Heat Generation/Absorption

Authors: Khilap Singh; Manoj Kumar
Journal: International Journal of Chemical Engineering
 (2016)

---

## Abstract

A numerical model is developed to examine the effects of thermal radiation on unsteady mixed convection flow of a viscous dissipating incompressible micropolar fluid adjacent to a heated vertical stretching surface in the presence of the buoyancy force and heat generation/absorption. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. The model contains nonlinear coupled partial differential equations which have been converted into ordinary differential equation by using the similarity transformations. The dimensionless governing equations for this investigation are solved by Runge-Kutta-Fehlberg fourth fifth-order method with shooting technique. Numerical solutions are then obtained and investigated in detail for different interesting parameters such as the local skin-friction coefficient, wall couple stress, and Nusselt number as well as other parametric values such as the velocity, angular velocity, and temperature.


---

## Body

## 1. Introduction

The micro polar fluids are those which contain microconstituents that can undergo rotation, the presence of which can affect the hydrodynamics of the flow. The classical Navier-Stokes theory does not describe the flow properties of micropolar fluids, for example, colloidal suspension, polymeric fluids, liquid crystals, fluids with additives, suspension solutions, animal’s blood, human blood, body fluids, biofluids, and fluids containing certain additives. Eringen [1] describes the theory of micropolar fluids, which show microrotation effects as well as microinertia. The theory of thermomicropolar fluids was developed by Eringen [2] by extending his theory of micropolar fluids. A good list of references for micropolar fluids is available in Łukaszewicz [3]. Many researchers [4–13] have studied the micropolar fluid flow for different fluid properties over different geometries.Unsteady mixed convection flow plays an important role in chemical engineering, turbomachinery, aerospace technology, geophysics, and so forth; Zueco et al. [14] studied the unsteady free convection flow of an MHD micropolar fluid through two parallel infinite porous vertical plates. Unsteady mixed convection flow of a micropolar fluid adjacent to a heated vertical surface along with viscous dissipation and the buoyancy force is analyzed by Abd El-Aziz [15]. Hussainn et al. [16] reported the radiation effects on the unsteady boundary layer flow of a micropolar fluid over a stretching permeable sheet. Oahimire and Olajuwon [17] studied the heat and mass transfer effects on an unsteady flow of a chemically reacting micropolar fluid over an infinite vertical porous plate. Rashad [18] studied the unsteady MHD boundary-layer flow and heat transfer for an electrically conducting rotating fluid due to a stretching surface in porous medium in the presence of thermal radiation. Abd El-Aziz [19] investigated the effects of variable viscosity on mixed convection flow along a semi-infinite unsteady stretching sheet with viscous dissipation.The heat transfer in the fluid flow due to a stretching sheet has attracted considerable attention during the last few decades due to its various applications in many industrial and engineering processes such as hot rolling, wire drawing, glass-fiber and paper production, drawing of plastic films, metal and polymer extrusion, and metal spinning. The pioneering work in this area was first made by Crane [20]; he studied boundary layer flow from a linearly stretching plate. Bhargava et al. [21] obtained a finite element solution for the mixed convection micropolar fluid flow in porous stretching sheet with suction. Eldabe and Ouaf [22] examined the heat and mass transfer flow of a micropolar fluid past a stretching surface with Ohmic heating and viscous dissipation. Pal et al. [23] analyzed the effects of thermal radiation and viscous dissipation on mixed convection flow of nanofluid over a nonlinear stretching/shrinking sheet. Singh and Kumar [24] studied boundary layer stagnation point flow of micropolar fluid towards a stretching/shrinking sheet in the presence of melting heat. Turkyilmazoglu [25] analyzed the micropolar fluid flow and heat transfer due to a porous stretching sheet.Influence of thermal radiation on flow and heat transfer study is much more important in different industries. The heat transfer and temperature profile of a fluid flow over different geometries can be affected significantly at high temperature. Mohamed and Abo-Dahab [26] studied thermal radiation effect on hydromagnetic free convection heat and mass transfer flow of a micropolar fluid over vertical porous plate. The micropolar fluid flow and heat transfer from a porous shrinking sheet were investigated by Bhattacharyya et al. [27]. Prakash and Muthtamilselvan [28] analyzed the effect of thermal radiation on fully developed flow of MHD micropolar fluid through two infinite parallel porous vertical plates.The effect of heat generation on heat transfer is an important issue in view of various physical problems. Ziabakhsh et al. [29] proposed micropolar fluid flow with heat generation. Bakr [30] investigated the effects of heat source on heat and mass transfer flow of a micropolar fluid in a rotating frame of reference. The heat generation/absorption effects on MHD flow of micropolar fluid through a stretching surface have been studied by Mahmoud and Waheed [31]. Abbasi et al. [32] investigated the Maxwell nanofluid fluid flow and heat transfer in the presence of heat generation/absorption. Mliki et al. [33] examined the influence of Brownian motion and heat generation/absorption over linear/sinusoidally heated cavity. Elgazery [34] analyzed the chemical reaction effect on MHD flow in the presence of temperature dependent viscosity and thermal diffusivity.The purpose of the present work is to study the effects of thermal radiation on mixed convection flow of a micropolar fluid through an unsteady stretching surface with viscous dissipation and heat generation/absorption. This problem is important in the processing of chemical engineering fluids including polymeric suspensions, lubricant manufacture, and so forth. The nonlinearity of basic equations associated with their inherent mathematical difficulties has led us to use numerical method. Thus the transformed dimensionless pertinent equations are solved numerically by using the Runge-Kutta-Fehlberg fourth fifth-order method along with shooting technique. The velocity, angular velocity, and temperature profiles are shown and the influences of the micropolar parameter, the thermal radiation parameter, the unsteadiness parameter, and the buoyancy parameter on the flow and heat transfer characteristics are discussed in detail. To the best of author’s knowledge such a study does not appear in the scientific literature.


## 2. Mathematical Formulation

We consider a two-dimensional unsteady mixed convection boundary layer flow of a viscous incompressible micropolar fluid over an elastic, vertical, and impermeable stretching sheet which emerges vertically in the upward direction from a narrow slot with velocity:(1)Uwx,t=ax1-αt,where both a and α are positive constants with dimension per time. The problem is in the presence of thermal radiation and heat generation/absorption. The flow configuration of this problem is illustrated in Figure 1. The positive x coordinate is measured along the stretching sheet with the slot as the origin and the positive y coordinate is measured normal to the sheet in the outward direction towards the fluid. The surface temperature Tw of the stretching sheet varies with the distance x from the slot and time t as (2)Twx,t=T∞+bx1-αt2,where b is a constant with dimension temperature over length and T∞ is the temperature of the ambient fluid. The expressions for Uw(x,t) and Tw(x,t) in (1) and (2) are valid only for time t<α-1 unless α=0. Further, expression (1) for velocity of the sheet Uw(x,t) reveals that the elastic surface which is fixed at the origin is stretched by applying force in the positive x-axis and the effective stretching rate a/(1-αt) increases with time. With the same analogy expression (2) for the surface temperature Tw(x,t) depicts a situation in which the surface temperature increases (or decreases) if b is positive (or negative) from T∞ at the slot in proportion to x and such that the amount of temperature increase (or decrease) along the sheet increases with time. It is further assumed that the fluid properties are taken to be constant except for the density variation with the temperature in the buoyancy term. Under these assumptions, the governing equations of boundary layer are given in the following form:(3)∂u∂x+∂v∂y=0,(4)∂u∂t+u∂u∂x+v∂u∂y=μ+κρ∂2u∂y2+κρ∂N∂x+g∗βT-T∞,(5)ρj∂N∂t+u∂N∂x+v∂N∂y=γ∂2N∂y2-κ2N+∂u∂y,(6)∂T∂t+u∂T∂x+v∂T∂y=kρCp∂2u∂y2+μ+κρCp∂u∂y2-1ρCp∂qr∂y+θ0ρCpT-T∞,where u and v are the components of velocity along x and y directions, respectively. Further μ is dynamic viscosity, κ is vortex viscosity, ρ is fluid density, g∗ is the gravitational acceleration, β is the volumetric coefficient of the thermal expansion, υ is kinematic viscosity, j is microinertia density, γ is the spin-gradient viscosity, N is the component of microrotation whose direction of rotation lies in the x-y plane, T is the temperature, k is the thermal conductivity of the fluid, Cp is the heat capacity at constant pressure p, qr is the radiative heat flux, θ0 is the heat generation/absorption coefficient, and T∞ is the temperature of the ambient fluid.Figure 1
Physical model and coordinate system.The appropriate boundary conditions for the problem are(7)u=Uwx,t,v=0,N=0,T=Twx,taty=0,u⟶0,N⟶0,T⟶T∞asy⟶∞.Using Roseland’s approximation, the radiative heat fluxqr is modeled as (8)qr=-4σ3k1∂T4∂y,where σ is the Stefan–Boltzmann constant and k1 is the absorption coefficient. Assuming that the difference in temperature within the flow is such that T4 can be expressed as a linear combination of the temperature, we expand T4 in Taylor’s series about T∞ as follows: (9)T4=T∞4+4T∞3T-T∞+6T∞2T-T∞2+⋯and neglecting higher order terms beyond the first degree in (T-T∞), we have(10)T4≈3T∞4+4T∞3T.Differentiating (8) with respect to y and using (10) we get(11)∂qr∂y=-16T∞3σ3k1∂2T∂y2.Using (11) in (6) we obtain (12)∂T∂t+u∂T∂x+v∂T∂y=kρCp∂2u∂y2+16σT∞33k1ρCp∂2T∂y2+μ+κρCp∂u∂y2+θ0ρCpT-T∞.The continuity equation (3) is satisfied by introducing the stream function ψ such that (13)u=∂ψ∂y=Uwf′η,v=-∂ψ∂x=-2υ1-αtfη.Equations (4), (5), and (12) can be transformed into a set of nonlinear ordinary differential equations by using the following similarity transformations:(14)η=aυ1-αty,ψ=υa1-αtxfη,N=a3υ1-αt3xgη,T=T∞+bx1-αt2θη.The transformed ordinary differential equations are (15)1+Kf′′′+ff′′-f′2-A22f′+ηf′′+Kg′+ξθ=0,λ0g′′η+fg′-f′h-A23g+ηg′-KB2g+f′′=0,3R+4θ′′+3RPrfθ′-f′g-A24θ+ηθ′+Ec1+Kf′′2+Hθ=0.The boundary conditions (7) reduce to(16)fη=0,f′η=1,gη=0,θη=1atη=0,f′η=0,gη=0,θη=0asη=∞,where prime denotes ordinary differentiation with respect to η, K=κ/ρ is the micropolar parameter, A=α/a is the unsteadiness parameter, θ=T-T∞/Tw-T∞ is the nondimensional temperature, ξ=gβb/a2 is the mixed convection or buoyancy parameter, λ0=γ/μj is the spin-gradient viscosity parameter, B=υ(1-αt)/jb=υx/jUw is the microinertia density parameter, Pr=υ/α is the Prandtl number, R=kk1/4σT∞3 is the thermal radiation parameter, Ec=Uw2/Cp(Tw-T∞) is the Eckert number, and H=θ0x/ρCpUw is the local heat generation/absorption parameter.It is worth mentioning that, forξ>0(Tw>T∞), buoyancy forces act in the direction of the mainstream and fluid is accelerated in the manner of a favorable pressure gradient (assisting flow). When ξ<0(Tw<T∞), buoyancy forces oppose the stretching induced flow, retarding the fluid in the boundary layer, acting as an adverse pressure gradient (opposing flow). Further, according to the definition of the viscous dissipation parameter, Ec<0 in the case of assisting flow (ξ>0) and Ec<0 in the case of opposing flow (ξ<0).The most important physical quantities for the problem are the local skin-friction coefficientCfx, the local wall couple stress coefficient Mwx, and the local Nusselt number Nux which are defined by(17)Cfx=2ρUw2μ+κ∂u∂yy=0=21+Kf′′0,Mwx=γ∂N∂yy=0=γaUwυ1-αtg′0,Nux=-xTw-T∞∂T∂yy=0=-Rex1/2θ′0,where Rex=xUw/υ is the local Reynolds number.


## 3. Method of Solution

The nonlinear differential equations (15) subject to boundary conditions (16) have been solved by using Runge-Kutta-Fehlberg fourth fifth-order method along with the shooting technique. This method is based on the discretization of the problem domain and the calculation of unknown boundary conditions.The domain of the problem is discretized and the boundary conditions forη=∞ are replaced by f′(ηmax)=0, g(ηmax)=0, and θ(ηmax)=0, where η∞ is a sufficiently large value of η at which boundary conditions (16) are satisfied. The authors ran the computer code written in MATLAB for different values of ηmax and step size Δη. They have seen that that there is no or negligible change in the velocity, angular velocity, and temperature for values of η greater than 7. Therefore in the present paper we have set η∞=7 and step size Δη=0.001. To solve the problem the nonlinear equations (15) have been first converted into seven first-order linear ordinary differential equations. There are four initial conditions at η=0 and three boundary conditions at η=∞. To find the solution of the problem, one will need three more conditions at η=0, that is, the values of f′′(0),g′(0), and θ′(0). These conditions have been found by the shooting technique. Finally the problem has been solved by the Runge-Kutta-Fehlberg method along with calculated boundary conditions.


## 4. Results and Discussion

In order to validate the numerical results obtained, we compare our results with those reported by Abd El-Aziz [11] as shown in Table 1, and they are found to be in a favorable agreement. In the simulation the default values of the parameters are considered as λ0=0.3,B=0.1,Pr=0.72,Ec=0.05, and H=0 unless otherwise specified.Table 1
Comparison of-θ′(0) for different values of A and ξ when K=Ec=R=0 and Pr=1.
A ξ Abd El-Aziz [15] Present results 0 0 1.00000000 1.00000000 0 1 1.08727816 1.08727818 0 2 1.14233927 1.14233928 0 3 1.18529030 1.18529032 1 0 1.68199254 1.68199255 1 1 1.70391279 1.70391281Figures2–4 show the effects of micropolar parameter K on the velocityf′(η), microrotation g(η), and temperature θ(η) profiles in the absence of unsteadiness (A=0) when R=0.5 and ξ=0. Figure 2 depicts that, with an increase in values of K, the velocity increases slightly near the surface and after a certain distance it decreases, but at K=2 it increases rapidly. It is clear from Figure 3 that the microrotation g(η) increases with increase in K, while at K=0.2 it decreases quickly. Further we notice from this figure that the microrotation takes negative values for large values of K (0.1 and 0.3). Physically it means that there is a reversed flow for this case. The fluid temperature increases with an increase in the values of K, and however at K=0.2 it decreases rapidly as shown in Figure 4.Figure 2
Velocity profiles for various values ofK when A=0, R=0.5, and ξ=0.Figure 3
Microrotation profiles for various values ofK when A=0, R=0.5, and ξ=0.Figure 4
Temperature profiles for various values ofK when A=0, R=0.5, and ξ=0.Figures5–7 depict the variation of the velocity f′(η), microrotation g(η), and temperature θ(η) profiles for different values of the micropolar parameter K when A=0.4, R=0.5, and ξ=0. It is clear from Figure 5 that velocity f′(η) increases with increasing value of K. It is noticed from Figure 6 that an increase in the values of micropolar parameter K leads to rise in the values of the microrotation g(η), also it can be seen that peak point of this profile shifts to the surface by increasing values of K (0.1 and 0.2), while for K=0 microrotation g(η) has constant value zero. It is evident from Figure 7 that the temperature decreases with the increasing values of micropolar parameter K.Figure 5
Velocity profiles for various values ofK when A=0.4, R=0.5, and ξ=0.Figure 6
Microrotation profiles for various values ofK when A=0.4, R=0.5, and ξ=0.Figure 7
Temperature profiles for various values ofK when A=0, R=0.5, and ξ=0.Figure8 presents the effects of the thermal radiation parameter R on temperature profile θ(η) when K=0.2, A=0, and ξ=0. It is clear from this figure that the temperature profile reduces with increasing values of the thermal radiation parameter R.Figure 8
Temperature profiles for various values ofR when K=0.2, A=0, and ξ=0.Figures9–11 exhibit the velocity, microrotation, and temperature profiles for various values of unsteadiness parameter A when K=0.2, R=0.5, and ξ=0. It is noticed from Figure 9 that, with an increase in unsteadiness parameter, the velocity f′(η) decreases and after a certain distance away from the surface it increases. It is seen from Figure 10 that the microrotation g(η) is initially smaller, in particular when η is less than approximately 4 and then g(η) increases slightly for larger A. Furthermore, the microrotation increases to a maximum point near the surface from which it starts to decrease and reaches zero near the free stream with the increasing values of η. It is clear from Figure 11 that the temperature reduces monotonically with an increase in the values of A.Figure 9
Velocity profiles for various values ofA when K=0.2, R=0.5, and ξ=0.Figure 10
Microrotation profiles for various values ofA when K=0.2, R=0.5, and ξ=0.Figure 11
Temperature profiles for various values ofA when K=0.2, R=0.5, and ξ=0.Figures12–14 depict the effects of buoyancy parameter ξ on the velocity f′(η), microrotation g(η), and temperature θ(η) profiles when A=0.2, R=0.5, and K=0.2. For assisting flow ξ>0, it is seen from Figure 12 that velocity increases with ξ; that is, increasing of ξ has tendency to induce more flow in the boundary layer. For the opposing flows (ξ<0), on the other hand the effect of buoyancy is to reduce the velocity compared to those for pure forced convection (ξ=0). This is because a positive ξ induces favorable pressure gradient that enhances the flow in boundary layer, while a negative ξ produces an adverse pressure gradient that slows down the fluid flow. Also for ξ=-0.2 the dimensionless fluid velocity f′(η) decreases rapidly with η first, arriving at a negative minimum value, and then increases to its free surface value. It is noted from Figure 13 that the angular velocity g(η) first decreases with ξ near the sheet surface where 0≤η<η0≅2.4 and increases with ξ for η>η0. Also it is seen from this figure that for ξ=-2 the angular velocity g(η) increases greatly near the sheet surface where 0≤η<η0≅1.25, whereas for η>η0 it decreases. Figure 14 shows that the effects of buoyancy parameter ξ are to decrease the temperature θ(η) in the case of assisting flow ξ>0 and increase it in the case of opposing flow ξ<0.Figure 12
Velocity profiles for various values ofξ when A=0.2, R=0.5, and K=0.2.Figure 13
Microrotation profiles for various values ofξ when A=0.2, R=0.5, and K=0.2.Figure 14
Temperature profiles for various values ofξ when A=0.2, R=0.5, and K=0.2.The effects of the unsteadiness parameterA, the buoyancy parameter ξ, the micropolar parameter K, and the thermal radiation parameter R on local skin-friction coefficient f′′(0), local wall couple stress g′(0), and local Nusselt number -θ′(0) are given in Table 2. It is clear from Table 2 that the local skin-friction coefficient increases with an increase in the values of K and in the case of the assisting flow ξ>0, while it decreases with increasing values of A and in the case of opposing flow ξ<0. Further, it is seen from Table 2 that the effect of increasing values of K is to increase local couple stress coefficient, whereas increasing A and ξ decreases local couple stress coefficient and there is no variation in local couple stress coefficient with R. In addition, the heat transfer rate increases with an increase in A,R and in the case of the assisting flow ξ>0. However, the heat transfer rate decreases with increasing the value of K and in the case of opposing flow ξ<0.Table 2
Computed values off′′(0),g′(0), and -θ′(0) for various values of A,ξ,K, and R.
A ξ K R f ′ ′ ( 0 ) g ′ ( 0 ) - θ ′ ( 0 ) 0 0 0.2 0.5 −0.91248 0.027431 0.36217 0.2 0 0.2 0.5 −0.97460 0.024750 0.45577 0.4 0 0.2 0.5 −1.03540 0.022970 0.53131 0.2 −0.2 0.2 0.5 −1.16166 0.034717 0.38816 0.2 −0.1 0.2 0.5 −1.04170 0.026600 0.44245 0.2 0.05 0.2 0.5 −0.94316 0.023959 0.46128 0.2 0.1 0.2 0.5 −0.91260 0.023225 0.46628 0 0 0 0.5 −1.08705 0 0.243760 0 0 0.1 0.5 −1.04570 0.026701 0.240340 0 0 0.3 0.5 −0.97566 0.072558 0.232590 0 0 0.2 0.2 −0.91248 0.027431 0.251050 0 0 0.2 1 −0.91248 0.027431 0.475040


## 5. Conclusions

The present work deals with the numerical analysis of thermal radiation effects of a mixed convection flow over an unsteady stretching surface. Fluid is a micropolar fluid in the presence of viscous dissipation and heat generation/absorption. The relevant nonlinear partial differential equations were transformed to a set of ordinary differential equations and then are solved numerically using the Runge-Kutta-Fehlberg fourth fifth-order method along with shooting technique. Conclusions drawn from the numerical results are as follows:(i)
Temperature reduces with increase in the values of the unsteadiness parameter, buoyancy parameter, and thermal radiation parameter.(ii)
The skin-friction is enhanced with an increase in the values of the micropolar parameter and buoyancy parameter, while it decreases with the increase of unsteadiness parameter.(iii)
The increasing value of the micropolar parameter is to increase the couple stress, whereas increasing unsteadiness and buoyancy parameter decreases the couple stress.(iv)
The rate of heat transfer increases with the unsteadiness parameter, buoyancy parameter, and thermal radiation parameter; however rate of heat transfer decreases with increasing micropolar parameter and thermal radiation parameter.



---

*Source: 8190234-2016-12-01.xml*