Journal of Computational Engineering

Volume 2015, Article ID 293105, 6 pages

http://dx.doi.org/10.1155/2015/293105

## New Application for the Generalized Incomplete Gamma Function in the Heat Transfer of Nanofluids via Two Transformations

Department of Mathematics, Faculty of Science, University of Tabuk, P.O. Box 741, Tabuk 71491, Saudi Arabia

Received 20 September 2014; Accepted 20 February 2015

Academic Editor: Clement Kleinstreuer

Copyright © 2015 Abdelhalim Ebaid and Hibah S. Alhawiti. 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

The boundary layer flow of nanofluids is usually described by a system of nonlinear differential equations with infinity boundary conditions. These boundary conditions at infinity are transformed into classical boundary conditions via two different transformations. Accordingly, the original heat transfer equation is changed into a new one which is expressed in terms of the new variable. The exact solutions have been obtained in terms of the exponential function for the stream function and in terms of the incomplete Gamma function for the temperature distribution. Furthermore, it is found in this project that a certain transformation reduces the computational work required to obtain the exact solution of the heat transfer equation. Hence, such transformation is recommended for future analysis of similar physical problems. Besides, the other published exact solution was expressed in terms of the WhittakerM function which is more complicated than the generalized incomplete Gamma function of the current analysis. It is important to refer to the fact that the analytical procedure followed in our project is easier and more direct than the one considered in a previous published work.

#### 1. Introduction

Nanofluids are a relatively new area of research which attracted attention in recent years because of their applications in engineering and applied sciences. The flow and heat transfer of such nanofluids are usually described by a system of nonlinear differential equations. Such system can be solved using numerical methods [1] or series methods such as Adomian decomposition method (ADM) [2–7], differential transformation/Taylor method (DTM) [8, 9], and homotopy perturbation method (HPM) [1, 10, 11]. The solutions using any of the just mentioned methods are usually called approximate numerical or analytical solutions. These approximate solutions have been extensively used to investigate various physical models because of the nonavailability of the exact solutions of these models. However, the exact solution of any system is the best if it can be obtained. Such exact solution gives us a better understanding about the phenomena involved in the physical model than the approximate solutions. Due to the difficulty of obtaining the exact solution many authors implement either approximate numerical methods or approximate analytical methods. However, we are interested in this paper in the exact solution of an example from nanofluid mechanics. Hence, an analytical approach will be suggested with the help of two transformations to facilitate the task of the current study. The example considered in this paper is governed by the following system of nonlinear differential equations [12, 13]:which represent the steady laminar two-dimensional flow of an incompressible viscous nanofluid past a linearly semi-infinite stretching sheet under the influence of a constant magnetic field. The primes denote the differentiation with respect to a similarity variable . and are the dimensionless stream function and temperature, respectively, is the solid volume fraction, and are densities, and are the heat capacitances, is the magnetic parameter, is the Prandtl number, and is the thermal conductivity defined as follows [14, 15]:where and are the thermal conductivities, where and denote the basic fluid and solid fractions, respectively. The flow is subject to the boundary conditionsThe exact solution of the stream function will be discussed in the next section. In the subsequent section, we discuss the effectiveness of two transformations in obtaining the exact solution of the heat transfer equation to get for the present physical model.

#### 2. Exact Solution for the Stream Function

Equations (1) can be rewritten aswhereIn order to solve system (5)-(6) with the boundary conditions (3)-(4) we begin with solving the -differential equation. To do that, the following assumption is assumed:Here, it should be noted that the infinity boundary condition in (3) is already satisfied provided that . On using (8) into (5), we havewhich leads toApplying the first two boundary conditions given in (3), we obtainIn view of (10)-(11), we haveTherefore,Inserting (13) into (6) yieldsRegarding the heat transfer equation, that is, (14), we will discuss in the next section the use of two transformations for obtaining the exact solution in terms of the generalized incomplete Gamma function.

#### 3. Exact Solution for the Heat Transfer

##### 3.1. Transformation 1

Suppose the following transformation: . The unbounded domain of the independent variable can be changed into a bounded one by using a new independent variable (say) using the transformation [16]Accordingly, the governing system should be expressed in terms of the new variable . In order to do that, we introduce the following relations between the derivatives with respect to and the derivatives with respect to :The relations given by (16) are obtained by using the chain rule in the differential calculus. Therefore, (14) becomeswhich can be simplified tosubject to the following set of boundary conditions:From (18) and using the separation of variables, we getorIntegrating (21) once with respect to from 0 to , we haveThereforeIntegrating (23) once with respect to from 0 to , we obtainIn view of the first condition in (19), we havewhereThe integration in the right-hand side can be analytically solved in terms of a well-known special function as declared by the following procedure. We first suppose thatAccordinglyWe then assume thatAccordinglyUsing the definition of the generalized incomplete Gamma function, we haveOn inserting (31) into (25), we getApplying the boundary condition , we obtainFrom (32) and (33), we havewhich can be simplified toand in terms of , we have the following final form for the temperature distribution:This exact solution clearly satisfies the boundary conditions and can be easily verified by direct substitution into (14). It is also important here to mention that the exact solution of the heat transfer equation has been already obtained in [12] in terms of the WhittakerM function. More discussion about the advantages and effectiveness of the current approach will be clarified later.

##### 3.2. Transformation 2

Suppose the following transformation: ; hence we haveTherefore, (14) becomessubject to the following set of boundary conditions:We rewrite (38) asIntegrating (40) once with respect to from 0 to , we haveorIntegrating (42) once with respect to from 0 to and using the boundary condition , we obtainWe then assume thatAccordinglyIn terms of the generalized incomplete Gamma function, (45) can be rewritten asApplying the boundary condition yields

Inserting this value of into (46), we have

which in terms of gives the following exact solution:Here, we can easily observe that less computational work is needed to get the exact solution of the heat transfer equation. Therefore, the transformation may be recommended for any future analysis to analyze similar physical models.

#### 4. Comparison with Other Analytical Published Results

In the previous sections, we have obtained the exact solution of the heat transfer differential equation described by (14) via two different transformations. It has been also shown that the two transformations lead to the same exact solution. However, one of them, the latter one, was found more direct than the first and hence few steps were sufficient to achieve the exact solution. It has been also shown that the current exact solution was expressed in terms of the generalized incomplete Gamma function. In order to compare our analytical method with the one published by Hamad [12], let us recall the analytical solution obtained in [12] (equation (28)) in the formwhere are parameters having the following forms:where the WhittakerM function and the Hypergeometric function are defined as follows:In view of our exact solution (49) and the exact solution given by (50) and (51) obtained in [12], we may summarize the advantages of our approach as follows. (i)The current exact solution is expressed in terms of the generalized incomplete Gamma function, while the other published solution was expressed in terms of two special functions, namely, the WhittakerM and the Hypergeometric functions. The properties of these two latest functions are more difficult than the properties of the generalized incomplete Gamma function. Furthermore, our exact solution agrees with the exact solution obtained very recently by Ebaid et al. [13] (equation (23)) in the absence of the slip parameter.(ii)Our exact solution is expressed in the form of a very simple expression if compared with the complex published solution obtained in [12] and given here by (50) and (51).(iii)Our simple exact solution can be easily verified by inserting only one expression into the heat transfer equation (14), while the other published solution needs more effort for achieving this task.(iv)The verification of the boundary conditions can be checked in a very simple way through our exact solution.(v)Our exact solution can be easily plotted when compared with the other published solution. With a very simple programme using Mathematica, we have plotted two figures for the purpose of comparisons between the current results and those obtained in [12]. In Figure 1, the effect of the magnetic parameter on the temperature distribution for the Cu-water nanofluid is depicted at and . Figure 2 shows the effect of the solid volume fraction of the Cu-nanoparticles on the variation of the temperature distribution at and . It is observed from these two figures that there are full agreements between the current results and those obtained by Hamad [12] (see Figure 2 and Figure 3 in [12], resp.). However, our approach is analytically and numerically much easier than the one followed in [12].