ISRN ThermodynamicsVolume 2012 (2012), Article ID 181286, 6 pagesdoi:10.5402/2012/181286

Research Article

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

Karem Boubaker,¹ Botong Li,² Liancun Zheng,² Ali H. Bhrawy,³ and Xinxin Zhang⁴

¹ESSTT, Université de Tunis, 63 Rue Sidi Jabeur, 5100 Mahdia, Tunisia²Department of Mathematics and Mechanics, University of Science and Technology Beijing, Beijing 100083, China³Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia⁴School of Mechanical Engineering, University of Science and Technology Beijing, Beijing 100083, China

Received 6 January 2012; Accepted 31 January 2012

Academic Editors: S. Hashimoto and G. Polidori

Copyright © 2012 Karem Boubaker et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 and the kinematics’ viscosity is (where , is the density and is the consistency index). The case corresponds to a Newtonian fluid. is descriptive of pseudoplastic fluids while 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 as and the thermal diffusivity be defined as with as a positive constant ( 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 speed . In the presence of viscous dissipation, the laminar boundary layer equations are written as: where and are the velocity components along and directions, respectively. The boundary conditions are The following dimensionless variables are introduced: is defined as the Zheng number.

The boundary layer equations then become with the boundary conditions:

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, similarity variable , and dimensionless temperature function are defined as: where A, B, , and are constant to be determined, and denotes the dimensionless stream function. The and velocity component are

Substituting and into (8)–(11) combining with we arrive at the nonlinear boundary value problems of the form as:

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 using variables, respectively:

The corresponding boundary conditions are

We introduce the parameters and as

Then, the problem is to find parameters and 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 as, , , , and . Thus the following equations are to be solved:

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: with : 4j-order Boubaker polynomials, : minimal positive roots, : a prefixed integer, : minimal positive roots, and : unknown pondering real coefficients.

Using expression (29), the boundary conditions: are verified automatically with reference to the Boubaker polynomials properties:

with The final solution is derived by introducing the expression (29) in the system (15)–(18) and calculating the coefficients which minimize the functional determinant : where , , and are calculated using the properties expressed by (32)–(34).

The solution is consequently:

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). Figures 1–4 show the dimensionless temperature profiles for different power law index , the Zheng number and . 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 and .

Figure 1: Dimensionless temperature for .

Figure 2: Dimensionless temperature for,  .

Figure 3: Dimensionless temperature for .

Figure 4: Dimensionless temperature for .

Figure 1 depicts the effect of power-law index on the temperature profiles using the two methods (Section 3, Section 4). It is a known fact that increasing values of implies difference in temperature. The same is reiterated by Figure 1: the changing factor arises some perturbation in temperature profiles and lead to the varying for each .

Figure 2 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 . The different values of contribute little to the thickness of the thermal boundary layer.

Figures 3-4 show temperature profiles for fixed and with different values of . For each , the thermal diffusion ratio increases with the increasing in as recorded earlier by Zheng et al. [38] and Ece and Buyuk [19]. It is evident that small values of result in thinning of the thermal boundary layer. This is similar to the effects of Prandtl number on 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 number and are discussed. Two significant physical phenomena are observed, that is, large values of result in thicker thermal boundary layers, and, the different values of contribute little to the thickness of the thermal boundary layer.

Nomenclature

: | Fluid specific heat (J kg⁻¹ C⁻¹)
: | Stream function (dimensionless)
: | Consistency index
: | Thermal conductivity (W kg⁻¹ K⁻¹)
: | Zheng number (dimensionless)
: | Power-law index (dimensionless)
: | Temperature (K)
: | Velocity components along and directions (m s⁻¹)
: | Cartesian coordinates along the plate and normal to it (m).

Greek symbols

: | Fluid thermal diffusivity (m² s⁻¹)
: | Similarity variable (dimensionless)
: | Normalized temperature (dimensionless)
: | Temperature function (dimensionless)
: | Stream function (dimensionless)
: | Kinematic viscosity (m s⁻¹)
: | Density (Kg m⁻³)
: | Constant (dimensionless).

Subscripts

:
Conditions at the surface and in the free stream, respectively.

Acknowledgment

The work is supported by the National Natural Science Foundations of China (No. 50936003), the Open Project of State Key Lab. for Adv. Metals and Materials (2009Z-02), and Research Foundation of Engineering Research Institute of USTB.

References

W. R. Schowalter, “The application of boundary-layer theory to power-law pseudo plastic fluids: similar solutions,” AIChE Journal, vol. 6, pp. 24–28, 1960.

A. Acrivos, M. J. Shah, and E. E. Petersen, “Momentum and heat transfer in laminar Boundary-Layer flows of non-Newtonian fluids past external surfaces,” AIChE Journal, vol. 6, pp. 312–317, 1960.

T. Y. Wang, “Mixed convection from a vertical plate to non-Newtonian fluids with uniform surface heat flux,” International Communications in Heat and Mass Transfer, vol. 22, no. 3, pp. 369–380, 1995.

T. Y. Wang, “Mixed convection heat transfer from a vertical plate to non-Newtonian fluids,” International Journal of Heat and Fluid Flow, vol. 16, no. 1, pp. 56–61, 1995.

F. M. Hady, “Mixed convection Boundary-Layer flow of non-Newtonian fluids on a horizontal plate,” Applied Mathematics and Computation, vol. 68, no. 2-3, pp. 105–112, 1995.

K. Kumari, I. Pop, and H. S. Takhar, “Free convection in power-law fluids near a three-dimensional stagnation point,” International Journal of Heat and Fluid Flow, vol. 18, pp. 625–631, 1997.

T. G. Howell, D. R. Jeng, and K. J. De Witt, “Momentum and heat transfer on a continuous moving surface in a power law fluid,” International Journal of Heat and Mass Transfer, vol. 40, no. 8, pp. 1853–1861, 1997.

J. H. Rao, D. R. Jeng, and K. J. De Witt, “Momentum and heat transfer in a power-law fluid with arbitrary injection/suction at a moving wall,” International Journal of Heat and Mass Transfer, vol. 42, no. 15, pp. 2837–3847, 1999.

I. A. Hassanien, A. A. Abdullah, and R. S. R. Gorla, “Flow and heat transfer in a power-law fluid over a nonisothermal stretching sheet,” Mathematical and Computer Modelling, vol. 28, no. 9, pp. 105–116, 1998.

N. Luna, F. Méndez, and C. Treviño, “Conjugated heat transfer in circular ducts with a power-law laminar convection fluid flow,” International Journal of Heat and Mass Transfer, vol. 45, no. 3, pp. 655–666, 2002.

A. Pinarbasi and M. Imal, “Viscous heating effects on the linear stability of Poiseuille flow of an inelastic fluid,” Journal of Non-Newtonian Fluid Mechanics, vol. 127, no. 2-3, pp. 67–71, 2005.

C.-H. Chen, “Effect of viscous dissipation on heat transfer in a non-Newtonian liquid film over an unsteady stretching sheet,” Journal of Non-Newtonian Fluid Mechanics, vol. 135, no. 2-3, pp. 128–135, 2005.

H. I. Andersson, K. H. Bech, and B. S. Dandapat, “Magnetohydrodynamic flow of a power-law fluid over a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 27, no. 6, pp. 929–936, 1992.

R. Cortell, “A note on magnetohydrodynamic flow of a power-law fluid over a stretching sheet,” Applied Mathematics and Computation, vol. 168, no. 1, pp. 557–566, 2005.

M. A. A. Mahmoud and M. A. E. Mahmoud, “Analytical solutions of hydromagnetic Boundary-Layer flow of a non-Newtonian power-law fluid past a continuously moving surface,” Acta Mechanica, vol. 181, no. 1-2, pp. 83–89, 2006.

I. Pop, M. Rashidi, and R. S. R. Gorla, “Mixed convection to power-law type non-Newtonian fluids from a vertical wall,” Polymer-Plastics Technology and Engineering, vol. 30, no. 1, pp. 47–66, 1991.

I. Pop, “Boundary layer flow at a three-dimensional stagnation point in power-law non-Newtonian fluids,” International Journal of Heat and Fluid Flow, vol. 14, no. 4, pp. 408–412, 1993.

R. S. R. Gorla, I. Pop, and J. K. Lee, “Convective wall plume in power-law fluid: second-order correction for the adiabatic wall,” Warme-und Stoffubertragung, vol. 27, no. 8, pp. 473–479, 1992.

M. C. Ece and E. Buyuk, “Similarity solutions for free convection to power-law fluids from a heated vertical plate,” Applied Mathematics Letters, vol. 33, no. 1, pp. 119–125, 1990.

L. Zheng, X. Zhang, and C. Lu, “Heat transfer for power law non-Newtonian fluids,” Chinese Physics Letters, vol. 23, no. 12, pp. 3301–3304, 2006.

L. Zheng, X. Zhang, and L. Ma, “Fully developed convective heat transfer of power law fluids in a circular tube,” Chinese Physics Letters, vol. 25, no. 1, pp. 195–197, 2008.

A. Belhadj, O. Onyango, and N. Rozibaeva, “Boubaker polynomials expansion scheme-related heat transfer investigation inside keyhole model,” Journal of Thermophysics and Heat Transfer, vol. 23, no. 3, pp. 639–640, 2009.

A. S. Kumar, “An analytical solution to applied mathematics-related Love's equation using the Boubaker polynomials expansion scheme,” Journal of the Franklin Institute, vol. 347, no. 9, pp. 1755–1761, 2010.

A. Belhadj, J. Bessrour, M. Bouhafs, and L. Barrallier, “Experimental and theoretical cooling velocity profile inside laser welded metals using keyhole approximation and Boubaker polynomials expansion,” Journal of Thermal Analysis and Calorimetry, vol. 97, no. 3, pp. 911–920, 2009.

P. Barry and A. Hennessy, “Meixner-type results for Riordan arrays and associated integer sequences, section 6: the Boubaker polynomials,” Journal of Integer Sequences, vol. 13, pp. 1–34, 2010.

M. Agida and A. S. Kumar, “A Boubaker polynomials expansion scheme solution to random Love's equation in the case of a rational Kernel,” Electronic Journal of Theoretical Physics, vol. 7, no. 24, pp. 319–326, 2010.

A. Yildirim, S. T. Mohyud-Din, and D. H. Zhang, “Analytical solutions to the pulsed Klein-Gordon equation using modified variational iteration method (MVIM) and Boubaker polynomials expansion scheme (BPES),” Computers and Mathematics with Applications, vol. 59, no. 8, pp. 2473–2477, 2010.

J. Ghanouchi, H. Labiadh, and K. Boubaker, “An attempt to solve the heat transfer equation in a model of pyrolysis spray using 4q-order m-Boubaker polynomials,” International Journal of Heat and Technology, vol. 26, no. 1, pp. 49–53, 2008.

O. B. Awojoyogbe and K. Boubaker, “A solution to Bloch NMR flow equations for the analysis of hemodynamic functions of blood flow system using m-Boubaker polynomials,” Current Applied Physics, vol. 9, no. 1, pp. 278–288, 2009.

S. Slama, J. Bessrour, and K. B. Ben Mahmoud, “Study of temperature 3D profile during weld heating phase using Boubaker polynomials expansion,” Thermochimica Acta, vol. 55, pp. 401–404, 2009.

S. Slama, M. Bouhafs, and K. B. Ben Mahmoud, “A Boubaker polynomials solution to heat equation for monitoring A3 point evolution during resistance spot welding,” International Journal of Heat and Technology, vol. 26, no. 2, pp. 141–146, 2008.

S. Tabatabaei, T. Zhao, O. Awojoyogbe, and F. Moses, “Cut-off cooling velocity profiling inside a keyhole model using the Boubaker polynomials expansion scheme,” Heat and Mass Transfer/Waerme- und Stoffuebertragung, vol. 45, no. 10, pp. 1247–1251, 2009.

S. Fridjine and M. Amlouk, “A new parameter: an ABACUS for optimizing PVT hybrid solar device functional materials using the Boubaker polynomials expansion scheme,” Modern Physics Letters B, vol. 23, no. 17, pp. 2179–2182, 2009.

H. Rahmanov, “A solution to the non linear korteweg-de-vries equation in the particular case dispersion-adsorption problem in porous media using the spectral Boubaker polynomials expansion scheme (BPES),” Studies in Nonlinear Sciences, vol. 2, no. 1, pp. 46–49, 2011.

A. Milgram, “The stability of the Boubaker polynomials expansion scheme (BPES)-based solution to Lotka-Volterra problem,” Journal of Theoretical Biology, vol. 271, no. 1, pp. 157–158, 2010.

M. Benhaliliba, C. E. Benouis, K. Boubaker, et al., “A new guide to thermally optimized doped oxides monolayer spray-grown solar cells: the amlouk-Boubaker optothermal expansivity ψab,” in Solar Cells, L. A. Kosyachenko, Ed., pp. 27–41, InTech, 2011.

M. Dada, M. A. Aweda, O. B. Awojoyogbe, and K. Boubaker, “A polynomial expression to the magnetic phase-shift induced in leon-vigmond model of the human heart,” Physical Journal of Mechanics in Medicine an Biology, vol. 12, no. 1, pp. 1–7, 2012.

L. C. Zheng, X. X. Zhang, K. Boubaker, U. Yücel, E. Gargouri-Ellouze, and A. Yıldırım, “Similarity and Boubaker polynomials expansion scheme BPES comparative solutions to the heat transfer equation for incompressible non-newtonian fluids: case of laminar boundary energy equation,” The European Physical Journal, vol. 55, no. 21102, 5 pages, 2011.