Journal of Computational Engineering
Volume 2015 (2015), Article ID 293105, 6 pages
http://dx.doi.org/10.1155/2015/293105

Research Article

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

Abdelhalim Ebaid and Hibah S. Alhawiti

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 conditions

The 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 as

where

In 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 have

which leads to

Applying the first two boundary conditions given in (3), we obtain

In view of (10)-(11), we have

Therefore,

Inserting (13) into (6) yields

Regarding 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) becomes

which can be simplified to

subject to the following set of boundary conditions:

From (18) and using the separation of variables, we get

or

Integrating (21) once with respect to from 0 to , we have

Therefore

Integrating (23) once with respect to from 0 to , we obtain

In view of the first condition in (19), we have

where

The 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 that

Accordingly

We then assume that

Accordingly

Using the definition of the generalized incomplete Gamma function, we have

On inserting (31) into (25), we get

Applying the boundary condition , we obtain

From (32) and (33), we have

which can be simplified to

and 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 have

Therefore, (14) becomes

subject to the following set of boundary conditions:

We rewrite (38) as

Integrating (40) once with respect to from 0 to , we have

or

Integrating (42) once with respect to from 0 to and using the boundary condition , we obtain

We then assume that

Accordingly

In terms of the generalized incomplete Gamma function, (45) can be rewritten as

Applying 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 form

where 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].

Figure 1: Effect of on temperature distribution for and .

Figure 2: Effect of on temperature distribution for and .

4.1. Conclusion

The nonlinear differential equations governing the flow and heat transfer of nanofluids in the presence of a magnetic field have been solved exactly. Two transformations have been suggested to change the domain from unbounded domain into a bounded one. It was observed that one of the two transformations is easier than the other where few steps were required in getting the exact solution. Moreover, the obtained exact solution for the heat transfer equation agreed with the results in literature at a special case. Furthermore, the other previous exact solution published in [12] (equation (23)) was expressed in terms of the WhittakerM function which is more complex special function than the generalized incomplete Gamma function of the current analysis. It may be reasonable to refer here to that the analytical procedure followed in our paper is easier and more direct than the one considered in [12].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

E. H. Aly and A. Ebaid, “New analytical and numerical solutions for mixed convection boundary-layer nanofluid flow along an inclined plate embedded in a porous medium,” Journal of Applied Mathematics, vol. 2013, Article ID 219486, 7 pages, 2013.

G. Adomian, Solving frontier problems of physics: the decomposition method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic, Boston, Mass, USA, 1994.

S. S. Siddiqi and M. Iftikhar, “Comparison of the Adomian decomposition method with homotopy perturbation method for the solutions of seventh order boundary value problems,” Applied Mathematical Modelling, vol. 38, no. 24, pp. 6066–6074, 2014.

A. Dib, A. Haiahem, and B. Bou-said, “An analytical solution of the MHD Jeffery-Hamel flow by the modified adomian decomposition method,” Computers & Fluids, vol. 102, pp. 111–115, 2014.

E. H. Aly, A. Ebaid, and R. Rach, “Advances in the Adomian decomposition method for solving two-point nonlinear boundary value problems with Neumann boundary conditions,” Computers & Mathematics with Applications, vol. 63, no. 6, pp. 1056–1065, 2012.

A.-M. Wazwaz, R. Rach, and J.-S. Duan, “Adomian decomposition method for solving the Volterra integral form of the Lane-Emden equations with initial values and boundary conditions,” Applied Mathematics and Computation, vol. 219, no. 10, pp. 5004–5019, 2013.

F. Shakeri Aski, S. J. Nasirkhani, E. Mohammadian, and A. Asgari, “Application of Adomian decomposition method for micropolar flow in a porous channel,” Original Research Article Propulsion and Power Research, vol. 3, no. 1, pp. 15–21, 2014.

G. Domairry and M. Hatami, “Squeezing Cu-water nanofluid flow analysis between parallel plates by DTM-Padé Method,” Journal of Molecular Liquids, vol. 193, pp. 37–44, 2014.

M. M. Rashidi, N. Laraqi, and S. M. Sadri, “A novel analytical solution of mixed convection about an inclined flat plate embedded in a porous medium using the DTM-Padé,” International Journal of Thermal Sciences, vol. 49, no. 12, pp. 2405–2412, 2010.

J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

R. A. van Gorder, “The variational iteration method is a special case of the homotopy analysis method,” Applied Mathematics Letters, vol. 45, pp. 81–85, 2015.

M. A. A. Hamad, “Analytical solution of natural convection flow of a nanofluid over a linearly stretching sheet in the presence of magnetic field,” International Communications in Heat and Mass Transfer, vol. 38, no. 4, pp. 487–492, 2011.

A. Ebaid, F. Al Mutairi, and S. M. Khaled, “Effect of velocity slip boundary condition on the flow and heat transfer of Cu-water and TiO₂-water nanofluids in the presence of a magnetic field,” Advances in Mathematical Physics, vol. 2014, Article ID 538950, 9 pages, 2014.

K. Khanafer, K. Vafai, and M. Lightstone, “Buoyancy-driven heat transfer enhancement in a two-dimensional enclosure utilizing nanofluids,” International Journal of Heat and Mass Transfer, vol. 46, no. 19, pp. 3639–3653, 2003.

K. Khanafer and K. Vafai, “A critical synthesis of thermophysical characteristics of nanofluids,” International Journal of Heat and Mass Transfer, vol. 54, no. 19-20, pp. 4410–4428, 2011.

E. H. Aly and A. Ebaid, “Exact solutions for boundary-layer flow of a nanofluid past a stretching sheet,” Journal of Computational and Theoretical Nanoscience, vol. 10, pp. 2591–2594, 2013.