Journal of Applied Mathematics

Volume 2013 (2013), Article ID 219486, 7 pages

http://dx.doi.org/10.1155/2013/219486

Research Article

New Analytical and Numerical Solutions for Mixed Convection Boundary-Layer Nanofluid Flow along an Inclined Plate Embedded in a Porous Medium

Emad H. Aly1,2 and Abdelhalim Ebaid³

¹Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia

²Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Cairo 11757, Egypt

³Department of Mathematics, Faculty of Science, University of Tabuk, Tabuk 71491, Saudi Arabia

Received 17 July 2013; Accepted 23 August 2013

Academic Editor: Mohamed Fathy El-Amin

Copyright © 2013 Emad H. Aly and Abdelhalim Ebaid. 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.

Two different analytical and numerical methods have been applied to solve the system describing the mixed convection boundary-layer nanofluids flow along an inclined plate embedded in a porous medium, namely, homotopy perturbation method (HPM) and Chebyshev pseudospectral differentiation matrix (ChPDM), respectively. Further, ChPDM is used as a control method to check the accuracy of the results obtained by HPM. The analytical method is applied using a new way for the deformed equations, and the resulted solution was expressed in terms of a well-known entire error function. In addition, using only two terms of the homotopy series, the approximate analytical solution is compared with the numerical solution obtained by the accurate ChPDM approach. The results reveal that good agreements have been achieved between the two approaches for various values of the investigated physical parameters.

1. Introduction

The concept of nanofluids is put into practice particularly after the tremendous development of nanotechnologies in the last decade, as they are of great importance in many industrial and engineering applications, especially in heat transfer enhancement [1–4]. Nanofluids are produced by dispersing nanometer-scale solid particles into base liquids such as water, ethylene glycol, and oils. Normally, if the particle sizes are in the 1–100 nm ranges, they are generally called nanoparticles.

Recently, the boundary-layer flow of nanofluids has attracted much attention [5–8]. The flow field of such problems is governed at first by a system of linear and nonlinear partial differential equations which are then transformed into a new system of linear and non-linear ordinary differential equations with applying the appropriate transformations. Due the difficulties that arise in obtaining the analytical solutions for the resulted system, most of the authors resort to numerical methods.

Although the numerical methods have the upper hand over the analytical methods in solving the complex systems of differential equations, the analytical solutions can be viewed as an optimal choice when available. In addition, the obtained solutions can be directly used for getting the required plots without needing subprograms as in the numerical methods. In conclusion, one can say that each of the two approaches has its own advantage over the other. Therefore, comparing the results obtained by two different analytical and numerical approaches is generally of great importance in the applied sciences. Two different analytical methods (homotopy perturbation method, HPM) [9–18] and numerical method (Chebyshev pseudospectral differentiation matrix, ChPDM) [19–21] are proposed in this paper for the mixed convection boundary-layer nanofluids flow along an inclined plate embedded in a porous medium. More details about these techniques and their advantages are introduced in the next two paragraphs.

The homotopy perturbation method (HPM) was proposed by He [9–15] as an analytical technique to solve nonlinear differential equations. This method has been widely used by many authors to investigate various models [16–18]. Unlike, the requirement for the regular perturbation techniques [19], the homotopy perturbation method is always valid no matter whether there exists small physical parameter or not. It combines the traditional perturbation method and the homotopy technique to deform a nonlinear problem into a simple solving one. The solution using this method is expressed as the summation of an infinite series, which usually converges rapidly to the exact solution in most cases. Several terms of the series solution can be used for approximating the exact solution with a high degree of accuracy.

On the other hand, the ChPDM is a known highly accurate numerical approach which was introduced and successfully applied by Aly et al. [20] to analyze the two-dimensional MHD boundary-layer flow over a permeable surface with a power law stretching velocity in the presence of a magnetic field applied normally to the surface. Under certain circumstances, it is shown that the problem has an infinite number of solutions which were examined by this technique. Further, Guedda et al. [21] have applied this method to validate and evidence the analysis of two-dimensional mixed convection boundary-layer flow over a vertical flat plate embedded in a porous medium saturated with water at 4°C (maximum density) and applied magnetic field. Both cases of the assisting and opposing flows are considered. Multiple similarity solutions are obtained and investigated by ChPDM under the power law variable wall temperature, or variable heat flux, or variable heat transfer coefficient. Recently, Aly and Ebaid [22] have applied this technique to solve the current problem in the absence of the magnetic and permeability of parameters. They found that the ChPDM results are very accurate in an excellent manner compared to those published in the literature using the homotopy analysis [23]. It is then successfully applied to validate and evidence the resulted exact solutions for different positive and negative values of the investigated parameters.

The motivation of presenting this paper is therefore to extend the applicability of both HPM and ChPDM approaches. The suggested procedure is based first on obtaining all the possible exact solutions. Then, the second step of the suggested procedure is to validate these results numerically to explore the effectiveness and efficiency of the proposed numerical approach. Besides, comparisons with other published results are also presented, where a full agreement is observed. In addition, it is shown that the solution up to the first order is sufficient to obtain accurate results.

2. The Physical Model

Rana et al. [24] considered the steady, incompressible, laminar, boundary-layer flow of a nanofluid along a semi-infinite inclined flat plate in a nanofluid-saturated porous medium with an acute angle to the vertical. The mathematical model is governed by the following system of nonlinear ordinary differential equations:

which has to be solved, subject to the boundary conditions:

where primes denoting differentiation with respect to and , , , , , and are the Lewis number, buoyancy ratio parameter, Brownian motion parameter, thermophoresis parameter, local Darcy-Rayleigh number, and local Péclet number, respectively. On integrating (1) with respect to from to and using the boundary conditions in (3)-(4), we have

3. Application of the HPM

In this section, the homotopy perturbation method (HPM) is applied to search for the approximate analytical solutions for the system of (6) and (2) with the boundary conditions (3)–(5). In order to apply this method, we rewrite the investigated system in terms of an embedding parameter as

According to the HPM, the following assumptions are assumed:

Substituting (8) into (7) and by comparing the like powers of the embedding parameter , we get the following systems:

and for , we have

The last systems are subjected to the boundary conditions:

3.1. Solution of the Zeroth-Order System

The exact solutions of the zeroth-order system are given as

where is the error function; see Abramowitz and Stegun [25]. These exact solutions can be verified by the direct substitution.

3.2. Solution of the First-Order System

The system of first-order is obtained by considering in (10) and given by

where the appropriated boundary conditions are

On using the initial condition (14), the exact solution of the -equation in system (13) is given as follows:

Regarding the view of -equation in (13) and boundary conditions (15), we obtain in the following closed form:

where is given by

where

The double integral in (19) can be expressed in terms of single integrations by using the integral identity:

and hence,

In a similar way, we obtain the solution of the -equation in (13) with the boundary conditions (16) in closed form as follows:

where is given by

where

Following the above analysis of in (22), can be given as

Therefore, at , the approximate solutions for the temperature distribution and the nanoparticle concentration up to the first order are given by

The current approximate solution obtained in this section is checked via comparing with those obtained in the literature.

3.3. Special Case: Horizontal Plate

At , we have the following simple expressions for the temperature distribution and the nanoparticle concentration distribution :

where

4. The ChPDM Technique

Chebyshev pseudospectral differentiation matrix (ChPDM) technique is briefly introduced; see Guedda et al. [21], Aly and Ebaid [22], and Van Gorder et al. [23]. On supposing that the domain of the problem is , then the following algebraic mapping transfers the domain to the Chebyshev one, that is . It is known that the Chebyshev polynomials are usually taken with their associated collocation points in the interval given by

Therefore, the th derivative of any function, say , at these collocation points can be approximated by the equation:

where is the Chebyshev pseudospectral approximation of , where and . The entries of the matrix are given by

where , except for , and

where and , , . The elements are the major elements concerning its values. Accordingly, they bear the major error responsibility when compared to the other elements. It is shown that the error in is of order , where is the machine precision [26].

On the view of [21–23], derivatives of the function , , and at the points are given by

Therefore, (1)–(5) become, respectively,

5. Results and Discussion

In Sections 3 and 4, mixed convection boundary-layer nanofluid flow along an inclined plate embedded in a porous medium has been solved analytically using two terms of the HPM and numerically by the ChPDM technique, respectively. The solutions were obtained for the temperature distribution and nano-particle concentration distribution . In the current section, we aim to compare the results derived from HPM by those obtained by the ChPDM, which may be considered as a reference approach because it was successfully compared before with homotopy analysis method in [22]. In each forthcoming figure, it should be noted that solid and dotted curves refer to the results obtained by ChPDM and HPA, respectively.

Figure 1 shows the comparison between the HPM and ChPDM solutions for the nanoparticle concentration distribution at different values of Brownian motion parameter , where , , , , and . From this figure, it is clear that the profiles of are very close to each other and to those obtained in Figure 3 by Rana et al. [24].

Figure 1: Comparison of ChPDM (solid) and HPM (dotted) for the effect of the Brownian motion parameter on the nanoparticle concentration distribution at , , , , and .

The temperature distribution is depicted against at different values of thermophoresis parameter for , , , , and in Figure 2. It is also observed from it that excellent agreement is achieved on comparing with Figure 4 in [24]. These results may refer to that two terms approximate solutions of the HPM are sufficient to give accurate numerical results.

Figure 2: Comparison of ChPDM (solid) and HPM (dotted) for the effect of the thermophoretic parameter on the temperature distribution at , , , , and .

The last conclusion is also detected in Figure 3 for the effect of the buoyancy ratio parameter and the mixed convection parameter on the temperature distribution at , , , and . From this figure, we can notice excellent agreements between the two approaches and further of those plotted in Figure 6 by Rana et al. [24].

Figure 3: Comparison of ChPDM (solid) and HPM (dotted) for the effect of the buoyancy ratio parameter and the mixed convection parameter on the temperature distribution at , , , and .

Figure 4 indicates the effect of the inclination angle on the nanoparticle concentration distribution at , , , , and . It can be easily observed from this figure that the slight difference between the curves obtained via the two approaches increases with increasing , especially at . A similar observation is detected from Figure 5 for the effect of Lewis number on at , , , , and , where the slight difference between the curves increases as Lewis number decreases.

Figure 4: Comparison of ChPDM (solid) and HPM (dotted) for the effect of the angle of inclination on the nanoparticle concentration distribution at , , , , and .

Figure 5: Comparison of ChPDM (solid) and HPM (dotted) for the effect of the Lewis number on the nanoparticle concentration distribution at , , , , and .

From the above comparisons, one can conclude that although slight differences are observed between the two approaches in some cases, two terms approximate solutions for and are still effective and of course can be enhanced via adding more terms for the homotopy series solutions. Accordingly, the HPM can be used as a method of solution with highly trust for investigating several similar models by following the proposed analysis of the current paper.

6. Conclusion

In the present paper, the system of differential equations, describing the mixed convection boundary-layer nanofluids flow along an inclined plate embedded in a porous medium, has been successfully solved by applying two methods, analytically by the homotopy perturbation and numerically by Chebyshev pseudospectral differentiation matrix, which was used as a control technique. In order to obtain meaningful physical solution, a new way was deduced to formulate the deformed system of differential equations. Accordingly, the approximate solution was expressed in terms of an entire error function. This proposed way of applying HPM may be of interest for investigating similar models in the future. On applying HPM, it was found that two terms of the homotopy series are enough to achieve very good solutions when compared with the accurate results obtained numerically via ChPDM. This was investigated for various values of the appropriated physical parameters: temperature distribution and nanoparticle concentration distribution at variation of the thermophoretic parameter , Brownian motion parameter , buoyancy ratio parameter , angle of inclination , Lewis number , and mixed convection parameter .

Acknowledgment

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant no. (233/965/1433). The authors, therefore, acknowledge with thanks DSR technical and financial support.

References

S. U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition, pp. 99–105, ASME, San Francisco, Calif, USA, 1995.

S. U. S. Choi, Z. G. Zhang, W. Yu, F. E. Lockwood, and E. A. Grulke, “Anomalous thermal conductivity enhancement in nanotube suspensions,” Applied Physics Letters, vol. 79, no. 14, pp. 2252–2254, 2001.

H. Masuda, A. Ebata, K. Teramae, and N. Hishinuma, “Alterlation of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles (Dispersion of g-Al₂O₃SiO₂, and TiO₂ ultra-fine particles),” Netsu Bussei, vol. 7, pp. 227–233, 1993.

S. Lee, S. U.-S. Choi, S. Li, and J. A. Eastman, “Measuring thermal conductivity of fluids containing oxide nanoparticles,” Journal of Heat Transfer, vol. 121, no. 2, pp. 280–288, 1999.

A. V. Kuznetsov and D. A. Nield, “Natural convective boundary-layer flow of a nanofluid past a vertical plate,” International Journal of Thermal Sciences, vol. 49, no. 2, pp. 243–247, 2010.

W. A. Khan and I. Pop, “Boundary-layer flow of a nanofluid past a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 53, no. 11-12, pp. 2477–2483, 2010.

N. Bachok, A. Ishak, and I. Pop, “Boundary-layer flow of nanofluids over a moving surface in a flowing fluid,” International Journal of Thermal Sciences, vol. 49, no. 9, pp. 1663–1668, 2010.

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

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

J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

J.-H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003.

J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004.

J.-H. He, “Asymptotology by homotopy perturbation method,” Applied Mathematics and Computation, vol. 156, no. 3, pp. 591–596, 2004.

J.-H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87–88, 2006.

J.-H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005.

P. D. Ariel, “The three-dimensional flow past a stretching sheet and the homotopy perturbation method,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 920–925, 2007.

S. Pamuk and N. Pamuk, “He's homotopy perturbation method for continuous population models for single and interacting species,” Computers & Mathematics with Applications, vol. 59, no. 2, pp. 612–621, 2010.

H. Aminikhah, “The combined Laplace transform and new homotopy perturbation methods for stiff systems of ODEs,” Applied Mathematical Modelling, vol. 36, no. 8, pp. 3638–3644, 2012.

A. H. Nayfeh, Perturbation Methods, John Wiley & Sons, New York, NY, USA, 1973.

E. H. Aly, M. Benlahsen, and M. Guedda, “Similarity solutions of a MHD boundary-layer flow past a continuous moving surface,” International Journal of Engineering Science, vol. 45, no. 2-8, pp. 486–503, 2007.

M. Guedda, E. H. Aly, and A. Ouahsine, “Analytical and ChPDM analysis of MHD mixed convection over a vertical flat plate embedded in a porous medium filled with water at C,” Applied Mathematical Modelling, vol. 35, no. 10, pp. 5182–5197, 2011.

E. H. Aly and A. Ebaid, “On the exact analytical and numerical solutions of nano boundary-layer fluid flows,” Abstract and Applied Analysis, vol. 2012, Article ID 415431, 22 pages, 2012.

R. A. van Gorder, E. Sweet, and K. Vajravelu, “Nano boundary layers over stretching surfaces,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 6, pp. 1494–1500, 2010.

P. Rana, R. Bhargava, and O. A. Bég, “Numerical solution for mixed convection boundary layer flow of a nanofluid along an inclined plate embedded in a porous medium,” Computers & Mathematics with Applications, vol. 64, no. 9, pp. 2816–2832, 2012.

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, NY, USA, 4th edition, 1972.

E. M. E. Elbarbary and S. M. El-Sayed, “Higher order pseudospectral differentiation matrices,” Applied Numerical Mathematics, vol. 55, no. 4, pp. 425–438, 2005.