# Effects of Viscous Dissipation on the Thermal Boundary Layer of Pseudoplastic Power-Law Non-Newtonian Fluids Using Discretization Method and the Boubaker Polynomials Expansion Scheme

Authors: Karem Boubaker; Botong Li; Liancun Zheng; Ali H. Bhrawy; Xinxin Zhang
Journal: ISRN Thermodynamics
 (2012)

---

## Abstract

Heat transfer of pseudoplastic power-law non-Newtonian fluids aligned with a semi-infinite plate is studied. Unlike in most classical works, the effects of viscous dissipation which is coupled with the temperature-dependent thermal diffusivity are considered in the energy equation. The discretization method is used to convert the governing partial differential equations into a set of nonlinear ordinary differential equations. The solutions are presented numerically by using the shooting technique coupled with the Newtonian method and the Boubaker polynomials expansion scheme. The effects of power-law index and the Zheng number on the dynamics are analyzed. The associated heat-transfer characteristics are also tabulated in some domains of the said parameters.


---

## Body

## 1. Introduction

Since 1960, a considerable attention has been devoted to predict the drag force behavior and energy transport characteristics of the non-Newtonian fluid flows. The main reason for this is probably that fluids (such as molten plastics, pulps, slurries, emulsions), which do not obey the assumption of Newtonian fluids that the tress tensor is directly proportional to the deformation tensor, are found in various engineering applications. A variety of constitutive equations have been proposed to describe the flow and heat transfer non-Newtonian characteristics, among them the empirical Ostwaald-de Waele model, which is known as the so-called power-law model, gained much acceptance. For an incompressible power law non-Newtonian fluid, its shear stress is characterized asτxy=K|∂U/∂Y|n-1∂U/∂Y and the kinematics’ viscosity is ν=γ|∂U/∂Y|n-1 (where γ=K/ρ, ρ is the density and K is the consistency index). The case n=1 corresponds to a Newtonian fluid. 0<n<1 is descriptive of pseudoplastic fluids while n>1 describes dilatant’s fluids.Schowalter [1] and Acrivos et al. [2] successfully applied the power-law model to the boundary layer problems. The boundary-layer equations were formulated, and the conditions for the existence of similarity solutions were established. Following the pioneering works of Schowalter and Acrivos et al. [1, 2], Wang [3, 4], and Hady [5] analyzed the mixed convection heat transfer from a vertical or horizontal plate of non-Newtonian fluids with uniform surface heat flux.  Kumari et al. [6] considered the free-convection boundary-layer flow of a non-Newtonian fluid along a vertical wavy surface. Howell et al. [7] and Rao et al. [8] examined the momentum and heat transfer of the power-law fluids on a continuous moving surface. Hassanien et al. [9] investigated the heat transfer in power-law flow over a nonisothermal stretching sheet. Luna et al. [10] considered the case of conjugated heat transfer of a power law laminar fluid in circular ducts. Pinarbasi and Imal [11] extended their work to the viscous heat effects on the linear stability of Poiseuille flow of an inelastic fluid. Chen [12] analyzed the effect of viscous dissipation on heat transfer in non-Newtonian power law liquids over an unsteady stretching sheet. For the non-Newtonian power-law fluids, the hydrodynamic problem of the MHD boundary layer flow over a continuously moving surface has been dealt with by several authors (e.g., Andersson et al. [13], Cortell [14], and M. A. A. Mahmoud and M. A. E. Mahmoud [15]). The effect of magnetic field is found to decrease the velocity distribution and thus to increase the skin-friction coefficient.It seems that, in all the works cited above, the power law kinematic viscosity was applied only on the momentum equations and the thermal conductivity is still treated as constant. This is inconsistent with the known fact that the change of viscosity should affect both the momentum and heat equation of the power law flows. Pop et al. [16, 17] proposed a model which states that the thermal conductivity of non-Newtonian fluids was power-law-dependent on the velocity gradient. Based on this consumption, Gorla et al. and Ece and Buyuk [18, 19] performed a boundary-layer analysis for the free convection flow over a vertical flat plate embedded in a porous medium saturated by a power-law non-Newtonian fluid and gave the similarity solution to the classical boundary layer equations of the power-law wall plume problem.Motivated by works mentioned above, this paper takes the effects of power law kinematic viscosity into account in the boundary layer heat transfer equations. In terms of the analogy between viscosity diffusion and thermal diffusion, we assumed that the thermal conductivity will be chosen ask(T)=λ|∂T/∂Y|n-1 and the thermal diffusivity α be defined as α=ω|∂T/∂Y|n-1 with ω=λ/ρcP as a positive constant (cP is the specific heat at constant pressure) [20, 21]. We first formulate a thermal boundary layer equation for the power-law non-Newtonian fluids with a new variable thermal conductivity and provide similarity solutions.


## 2. Formalization of the Nonlinear Boundary Value Problem

Consider a semi-infinite plate at uniform wall temperature aligned with a power law fluid of constant speedU∞. In the presence of viscous dissipation, the laminar boundary layer equations are written as:(1)∂U∂X+∂V∂Y=0,(2)U∂U∂X+V∂U∂Y=1ρ∂τxy∂Y,(3)U∂T∂X+V∂T∂Y=∂∂Y(α∂T∂Y)+νcp(∂U∂Y)2,
where Uand Vare the velocity components along X and Ydirections, respectively. The boundary conditions are
(4)U|Y=0=0,V|Y=0=0,U|Y=+∞=U∞,(5)T|Y=0=Tw,T|Y=+∞=T∞.
The following dimensionless variables are introduced:(6)x=XL,y=[ρU∞2-nLnK]1/n+1YL,u=UU∞,v=[ρU∞2-nLnK]1/n+1VU∞,(7)θ=T-TwT∞-Tw,NZh=U∞n-1γ|T∞-Tw|n-1ω,Z*=U∞2cp(T∞-Tw),NZh is defined as the Zheng number.The boundary layer equations then become(8)∂u∂x+∂v∂y=0,(9)u∂u∂x+v∂v∂y=∂∂y(|∂u∂y|n-1∂u∂y),u∂θ∂x+v∂θ∂y=1NZh∂∂y(|∂θ∂y|n-1∂θ∂y)+Z*|∂u∂y|n-1(∂u∂y)2,
with the boundary conditions:
(10)u|y=0=0,v|y=0=0,u|y=+∞=1,(11)θ|y=0=0,θ|y=+∞=1.Equations (8)–(11) are super-nonlinear. Since the problem is very complex to solve analytically, a discretized model may be a considerable approach.The stream functionψ(x,y), similarity variable η, and dimensionless temperature function w(η) are defined as:(12)ψ=Axα̃f(η),θ(x,y)=w(η),η=Bxβ̃y,
where A,B, α̃, and β̃ are constant to be determined, and f(η) denotes the dimensionless stream function. The u  and v velocity component are(13)u=∂ψ∂y=ABxα̃+β̃f′(η),v=∂ψ∂x=-Aα̃xα̃-1[f(η)-ηf′(η)].Substitutingu and v into (8)–(11) combining with(14)β̃=-α̃,AB=1,α̃=1n+1,B=(n+1)(-1/(n+1)),
we arrive at the nonlinear boundary value problems of the form as:(15)(|f′′(η)|n-1f′′(η))′+f(η)f′′(η)=0,(16)f(0)=0,f′(0)=0,f′(η)|η=+∞=1,(17)f(η)w′(η)+1NZh(|w′(η)|n-1w′(η))′,+Z*|f′′(η)|n-1(f′′(η))2=0,(18)w(0)=0,w(η)|η=+∞=1.


## 3. Discretization Method Solution

In order to obtain numerical solutions, we transfer the problem (15)–(18) to a system of four first-order equations by denoting the f,f′,f′′,w,w′using variablesf,u,v,w,y, respectively:(19)f′=u,(20)u′=v,(21)v′=-1nfv2-n,(22)w′=y,(23)y′=-Z*⋅NZhny1-nvn+1-NZhnfy2-n.The corresponding boundary conditions are(24)f(0)=0,u(0)=0,w(0)=0(25)u(+∞)=1,w(+∞)=1.We introduce the parameterstand s as(26)v(0)=t,y(0)=s.Then, the problem is to find parameterst and s such that the solution to (19)–(24), (26) satisfies the boundary conditions (25).The solution to the initial-value problems (19)–(24), (26) can be denoted asf(η,t,s), u(η,t,s), v(η,t,s), w(η,t,s), and y(η,t,s). Thus the following equations are to be solved:(27)ϕ(t,s)=u(+∞,t,s)-1=0,(28)φ(t,s)=w(+∞,t,s)-1=0.The Newtonian method and Runge-Kutta scheme are convenient to solve the nonlinear problem above.


## 4. BPES Solution

The Boubaker polynomials expansion scheme (BPES) [22–37] solutions starts from assigning the following expression:(29)w(η)=12M∑j=1Mξj×B4j(ϖjε(η)),ε(η)=(η+1)e-η,
with B4j: 4j-order Boubaker polynomials, ϖj: B4j minimal positive roots, M: a prefixed integer, ϖj: B4j minimal positive roots, and ξj|j=1…M: unknown pondering real coefficients.Using expression (29), the boundary conditions:(30)w(η)|η=0=0,lim⁡η→+∞w(η)=1
are verified automatically with reference to the Boubaker polynomials properties:(31)∑q=1NB4q(x)|x=0=-2N≠0,∑q=1NB4q(x)|x=αq=0,(32)∑q=1NdB4q(x)dx|x=0=0,∑q=1NdB4q(x)dx|x=αq=∑q=1NHq,∑q=1Nd2B4q(x)dx2|x=0=83(N(N2-1)),∑q=1Nd2B4q(x)dx2|x=αq=∑q=1NGq,with(33)Hn=B4n′(αn)=(4αn[2-αn2]×∑q=1nB4q2(αn)B4(n+1)(αn)+4αn3),(34)Gq=d2B4q(x)dx2|x=αq=3αq(4qαq2+12q-2)Hq-8q(24q2αq2+8q2-3q+4)(αq2-1)(12qαq2+4q-2).
The final solution is derived by introducing the expression (29) in the system (15)–(18) and calculating the coefficients ξjoptim.|j=1…M which minimize the functional determinant Δ:(35)Δ=|(|∑j=1Mξjoptim.×Φj(2)|)n-1(∑j=1Mξjoptim.×Φj(2))-(∑j=1Mξjoptim.×Φj(0))(∑j=1Mξjoptim.×Φj(1))|,
where Φj(0), Φj(1), and Φj(2)are calculated using the properties expressed by (32)–(34).The solution is consequently:(36)wsol.(η)=12M∑j=1Mξjoptim.×B4j(ϖjε(η)),ε(η)=(η+1)e-η.


## 5. Numerical Results

The heat transfer of pseudoplastic power law fluids aligned with a semi-infinite plate is investigated in the presence of viscous dissipation. Numerical solutions are achieved using the discretization shooting method and the Boubaker polynomials expansion echeme (BPES). Figures1–4 show the dimensionless temperature profiles for different power law index n, the Zheng number NZh and Z*. It is observed that the length scale of thermal diffusion is significantly different for the values of parameters. On the whole, the thermal diffusion ratio increases with the increase of NZh and Z*.Figure 1
Dimensionless temperature forZ*=0.5,NZh=2.0.Figure 2
Dimensionless temperature forn=0.8,  NZh=2.0.Figure 3
Dimensionless temperature forn=0.8,Z*=1.0.Figure 4
Dimensionless temperature for=0.5,Z*=0.5.Figure1 depicts the effect of power-law index n on the temperature profiles using the two methods (Section 3, Section 4). It is a known fact that increasing values of n implies difference in temperature. The same is reiterated by Figure 1: the changing factor n arises some perturbation in temperature profiles and lead to the varying w for each η.Figure2 reveals the effect made by the viscous dissipation on the temperature. These profiles are well behaved and very little change occurs in the shapes of profiles with the varying parameter Z*. The different values of Z*contribute little to the thickness of the thermal boundary layer.Figures3-4 show temperature profiles for fixed Z*and n with different values of NZh. For each NZh(NZh>0), the thermal diffusion ratio increases with the increasing in NZh as recorded earlier by Zheng et al. [38] and Ece and Buyuk [19]. It is evident that small values of NZh result in thinning of the thermal boundary layer. This is similar to the effects of Prandtl number Pron heat transfer.


## 6. Conclusions

An investigation on heat transfer of pseudoplastic non-Newtonian fluids aligned with a semi-infinite plate is made. The boundary layer equations for both momentum and energy boundary layers have been developed. Numerical solutions are presented and the associated transfer behavior for different power-law index, and the Zheng numberNZh and Z* are discussed. Two significant physical phenomena are observed, that is, large values of NZh result in thicker thermal boundary layers, and, the different values of Z*contribute little to the thickness of the thermal boundary layer.



---

*Source: 181286-2012-02-29.xml*