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 [14]. 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 [58]. 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) [918] and numerical method (Chebyshev pseudospectral differentiation matrix, ChPDM) [1921] 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 [915] as an analytical technique to solve nonlinear differential equations. This method has been widely used by many authors to investigate various models [1618]. 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 [2123], 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.