# Heat Generation and Thermal Radiation Effects over a Stretching Sheet in a Micropolar Fluid

Authors: M. Gnaneswara Reddy
Journal: ISRN Thermodynamics
 (2012)

---

## Abstract

The effects of radiation and heat generation on steady thermal boundary layer flow induced by a linearly stretching sheet immersed in an incompressible micropolar fluid with constant surface temperature are investigated. Similarity transformation is employed to transform the governing partial differential equations into ordinary ones, which are then solved numerically using the Runge-Kutta fourth order along shooting method. Results for the local Nusselt number as well as the temperature profiles are presented for different values of the governing parameters. It is observed that the velocity increases with an increase in the material parameter. It is seen that the temperature profile is influenced considerably and increases when the value of heat generation parameter increases along the boundary layer. Also, the temperature distribution of the fluid increases with an increase in the radiation parameter. Comparisons with previously published work are performed and the results are found to be in very good agreement.


---

## Body

## 1. Introduction

Flow of a viscous fluid past a stretching sheet is a classical problem in fluid dynamics. The development of boundary layer flow induced solely by a stretching sheet was first studied by Crane [1] who first obtained an elegant analytical solution to the boundary layer equations for the problem of steady two dimensional flow due to a stretching surface in a quiescent incompressible fluid. Flow and heat transfer characteristics due to a stretching sheet in a stationary fluid occur in a number of industrial manufacturing processes and include both metal and polymer sheets, for example, the cooling of an infinite metallic plate in a cooling bath, the boundary layer along material handling conveyers, the aerodynamic extrusion of plastic sheets, paper production, metal spinning, and drawing plastic films. The quality of the final product depends on the rate of heat transfer at the stretching surface. This problem was then extended by P. S. Gupta and A. S. Gupta [2] to a permeable surface. The flow problem due to a linearly stretching sheet belongs to a class of exact solutions of the Navier-Stokes equations. Thus, the exact solutions reported by Crane [1] and P. S. Gupta and A. S. Gupta [2] are also the exact solutions to the Navier-Stokes equations. The heat transfer aspects of similar problems were studied by Grubka and Bobba [3], Chen and Char [4], Dutta et al. [5], Ali [6, 7], Afzal and Varshney [8], Afzal [9], and many others. On the other hand, the effects of buoyancy force on the development of velocity and thermal boundary layer flows over a stretching sheet have been investigated by Chen [10], Ali and Al-Yousef [11], Daskalakis [12], Partha et al. [13], Abd El-Aziz [14], Mahapatra et al. [15], and Ishak et al. [16–20], among others.The study of flow and heat transfer past a stretching sheet has gained tremendous interest among researchers due to its industrial and engineering applications. This includes extrusion of plastic sheets, annealing and tinning of copper wire, paper production, crystal growing, and glass blowing. The final products depend mainly on the stretching and cooling rates at the surface. Their studies are not only restricted to the Newtonian fluids but also include the non-Newtonian fluids such as micropolar fluids. Such studies have been carried out by Chiam [21], Heruska et al. [22], Agarwal et al. [23], Hassanian and Gorla [24], Kelson and Desseaux [25], Kelson and Farrell [26], Nazar et al. [27], and very recently by Hayat et al. [28].In the present paper we study the development of thermal boundary layer flow induced by a stretching sheet immersed in a micropolar fluid with the effect of heat generation and thermal radiation is taken into consideration. The governing partial differential equations are transformed into ordinary ones using similarity transformation, before being solved numerically by the Runge-Kutta fourth order along shooting method.


## 2. Mathematical Model

A steady laminar boundary layer flow over a stretching sheet immersed in a quiescent and incompressible micropolar fluid with uniform surface temperatureTw is considered. It is assumed that the sheet is stretched with a linear velocity Uw=ax, where a is a positive constant and x is the distance from the slit where the sheet is issued. The simplified two-dimensional equations governing the flow may be written asContinuity equation(1)∂u∂x+∂v∂y=0.Momentum equation(2)u∂u∂x+v∂u∂y=(ν+κρ)∂2u∂y2+κρ∂N∂y.Angular Momentum equation,(3)u∂N∂x+v∂N∂y=γρj∂2N∂y2-κρj(2N+∂u∂y).Energy equation,(4)u∂T∂x+v∂T∂y=kρcp∂2T∂y2-1ρcp∂qr∂y+Q0ρcp(T-T∞).The spin-gradient viscosityγ (Ahmadi [29]) can be defined as
(5)γ=(μ+κ2)j=μ(1+K2)j,
where μ is the dynamic viscosity, K=κ/μ is the dimensionless viscosity ratio and is called the material parameter, and we take j=ν/a as a reference length. The Equation (5) is invoked to allow the field of equations predicts the correct behavior in the limiting case when the microstructure effects become negligible and the total spin N reduces to the angular velocity.The boundary conditions for the velocity, angular velocity, and fields are(6)u=Uw,v=0,N=-m∂u∂y,T=Twaty=0u⟶0,N⟶0,T⟶T∞asy⟶∞,
where T∞ is the ambient fluid temperature and m is the boundary parameter with 0≤m≤1.By using the Rosseland approximation, Brewster [30], the radiative heat flux qr is given by
(7)qr=-4σs3ke∂T4∂y,
where σs is the Stefan-Boltzmann constant and ke is the mean absorption coefficient. It should be noted that by using the Rosseland approximation, the present analysis is limited to optically thick fluids. If the temperature differences within the flow are sufficiently small, then (7) can be linearized by expanding T4 into the Taylor series about T∞, which after neglecting higher order terms takes the form
(8)T4≅4T∞3T-3T∞4.
In view of (7) and (8), (4) reduces to
(9)u∂T∂x+v∂T∂y=α(1+R)∂2T∂y2+Q0ρcp(T-T∞),
where α=k/ρcp is the thermal diffusivity and R=16σ*T∞3/3kk* is the radiation parameter.We introduce now the following similarity transformation:(10)η=(Uwνx)1/2y,ψ=(νxUw)1/2f(η),N=Uw(Uwνx)1/2h(η),θ(η)=T-T∞Tw-T∞,Q=Q0νρcpUw2,
where η is the similarity variable and ψ is the stream function defined as u=∂ψ/∂y,v=-∂ψ/∂x which identically satisfies the mass conservation equation (1). Substituting (10) into (2), (3), and (9) we obtain the following ordinary differential equations:
(11)(1+K)f′′′+ff′′-f′2+Kh′=0,(12)(1+K2)h′′+fh′-f′h-K(2h+f′′)=0,(13)(1+R)θ′′+Prfθ′+PrQθ=0,
where primes denote differentiation with respect to η and Pr=μ/α is the Prandtl number.The corresponding dimensionless boundary conditions are(14)f=0,f′=1,h=-mf′′(0),θ=1atη=0f′⟶0,h⟶0,θ⟶0asη⟶∞.It is note thatK= 0 corresponds to viscous fluid.


## 3. Results and Discussion

The set of nonlinear ordinary differential equations (11)–(13) with boundary conditions (14) have been solved by using the Runge-Kutta fourth order along with Shooting method. First of all, higher order nonlinear differential Equations (11)–(13) are converted into simultaneous linear differential equations of first order and they are further transformed into initial value problem by applying the shooting technique (Jain et al. [31]). The resultant initial value problem is solved by employing Runge-Kutta fourth order technique. The step size Δη=0.01 is used to obtain the numerical solution with five decimal place accuracy as the criterion of convergence. To analyze the results, numerical computation has been carried out using the method described in the previous section for variations in the governing parameters, namely, material parameter K, radiation parameter R, heat generation parameter Q, the Prandtl number Pr, and boundary parameter m. In the present study following default parameter values are adopted for computations: K=1.0, R=1.0, Q=1.0, Pr=0.71, and m=0.5. All graphs therefore correspond to these values unless specifically indicated on the appropriate graph.In order to assess the accuracy of our computed results, the present result has been compared with Ishak [32] for different values of K as shown in Figure 1 with Q=0.0. It is observed that the agreements with the solution of velocity profiles are excellent.Figure 1
Comparison of temperature profiles.Figure2 presents the velocity profiles for various values of K when m=0.5. We note that the parameters R and Pr have no influence on the flow field, which is clear from (11)–(13). It is evident from this figure that the boundary layer thickness increases with K. The velocity gradient at the surface f′′(0) decreases (in absolute sense) as K increases. Thus, micropolar fluids show drag reduction compared to viscous fluids. The negative velocity gradient at the surface f′′(0)<0 as shown in Figure 2 means the stretching sheet exerts a drag force on the fluid. This is not surprising since the development of the boundary layer is solely induced by it. The effect of material parameter K on the temperature is shown in Figure 3. The temperature decreases with an increase in the material parameter.Figure 2
Velocity profiles for different values ofK.Figure 3
Temperature profiles for different values ofK.Figure4 shows the temperature profiles for different values of R. The radiation parameter R defines the relative contribution of conduction heat transfer to thermal radiation transfer. It is obvious that an increase in the radiation parameter results in increasing temperature within the boundary layer. The effect of heat generation parameter Q on the temperature is shown in Figure 5. From this figure, we observe that when the value of heat generation parameter increases, the temperature distribution also increases along the boundary layer.Figure 4
Temperature profiles for different values ofR.Figure 5
Temperature profiles for different values ofQ.The influence of the Prandtl numberPr on temperature field is shown in Figure 6. The numerical results show that the effect of increasing values of the Prandtl number results in a decreasing temperature. It is observed that an increase in the Prandtl number results in a decrease of the thermal boundary layer thickness and in general lower average temperature within the boundary layer. The reason is that smaller values of Pr are equivalent to increasing the thermal conductivities, and therefore heat is able to diffuse away from the heated surface more rapidly than for higher values of Pr. Hence in the case of the smaller Prandtl numbers as the boundary layer is thicker and the rate of heat transfer is reduced.Figure 6
Temperature profiles for different values ofPr.The effects ofm on velocity, angular velocity, and temperature profiles are depicted in Figures 7, 8, and 9, respectively Figure 7 shows that the velocity gradient at the surface is larger for larger values of m. Different behaviors are observed for the effect of m on the heat transfer rate at the surface as presented in Figure 9. As expected, the couple stress h(0) is more dominant for larger values of m, as shown in Figure 8.Figure 7
Velocity profiles for different values ofm.Figure 8
Angular velocity profiles for different values ofm.Figure 9
Temperature profiles for different values ofm.


## 4. Conclusions

A steady two-dimensional laminar and heat transfer due to a stretching sheet immersed in an incompressible micropolar fluid has been investigated. The effect of thermal radiation and heat generation on the development of the thermal boundary layer flow has been taken into consideration. The effects of the governing parametersK,R,Q,Pr, and m on the fluid flow and heat transfer characteristics are discussed. It is found that the temperature distribution increases as the radiation parameter R increases. It is observed that the velocity gradient at the surface is larger for larger values of m.



---

*Source: 795814-2012-08-23.xml*