Research Article | Open Access

Emad H. Aly, Abdelhalim 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. https://doi.org/10.1155/2013/219486

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

**Academic Editor:**Mohamed Fathy El-Amin

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

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.

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

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. View at: Google Scholar - 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. View at: Publisher Site | Google Scholar - 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
_{2}O_{3}SiO_{2}, and TiO_{2}ultra-fine particles),ā*Netsu Bussei*, vol. 7, pp. 227ā233, 1993. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - J.-H. He, āHomotopy perturbation technique,ā
*Computer Methods in Applied Mechanics and Engineering*, vol. 178, no. 3-4, pp. 257ā262, 1999. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | MathSciNet - J.-H. He, āHomotopy perturbation method: a new nonlinear analytical technique,ā
*Applied Mathematics and Computation*, vol. 135, no. 1, pp. 73ā79, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-H. He, āComparison of homotopy perturbation method and homotopy analysis method,ā
*Applied Mathematics and Computation*, vol. 156, no. 2, pp. 527ā539, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-H. He, āAsymptotology by homotopy perturbation method,ā
*Applied Mathematics and Computation*, vol. 156, no. 3, pp. 591ā596, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-H. He, āHomotopy perturbation method for solving boundary value problems,ā
*Physics Letters A*, vol. 350, no. 1-2, pp. 87ā88, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-H. He, āApplication of homotopy perturbation method to nonlinear wave equations,ā
*Chaos, Solitons and Fractals*, vol. 26, no. 3, pp. 695ā700, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. H. Nayfeh,
*Perturbation Methods*, John Wiley & Sons, New York, NY, USA, 1973. View at: MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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 ${4}^{\circ}$C,ā
*Applied Mathematical Modelling*, vol. 35, no. 10, pp. 5182ā5197, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

#### Copyright

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.