`Abstract and Applied AnalysisVolume 2012 (2012), Article ID 415431, 22 pageshttp://dx.doi.org/10.1155/2012/415431`
Research Article

## On the Exact Analytical and Numerical Solutions of Nano Boundary-Layer Fluid Flows

1Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Cairo 11757, Egypt
3Department of Mathematics, Faculty of Science, Tabuk University, Tabuk 71491, Saudi Arabia

Received 31 May 2012; Revised 24 July 2012; Accepted 25 July 2012

Copyright © 2012 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

The nonlinear boundary value problem describing the nanoboundary-layer flow with linear Navier boundary condition is investigated theoretically and numerically in this paper. The -expansion method is applied to search for the all possible exact solutions, and its results are then validated by the Chebyshev pseudospectral differentiation matrix (ChPDM) approach which has been recently introduced and successfully used. This numerical technique is firstly applied and, on comparing with the other recent work, it is found that the results are very accurate and effective to deal with the current problem. It is then used to examine and validate the present analytical analysis. Although the -expansion method has been used widely to solve nonlinear wave equations, its application for nonlinear boundary value problems has not been discussed yet, and the present paper may be the first to address this point. It is clarified that the exact solutions obtained via the -expansion method cannot be obtained by using some of the other methods. In addition, the domain of the physical parameters involved in the current boundary value problem is also discussed. Furthermore, the convex, vicinity of zero, and asymptotic solutions are deduced.

#### 1. Introduction

At nanoscale diameters, it was shown in an interesting paper by He et al. [1] that a fascinating phenomena arise when the diameter of the electrospun nanofibers is less than 100 nm. The nanoeffect has been demonstrated for unusual strength, high surface energy, surface reactivity, and high thermal and electric conductivity.

The notion of a boundary layer was first introduced by Prandtl [2] over a hundred years ago to explain the discrepancies between the theory of inviscid fluid flow and experiment. In the classical boundary-layer theory, the condition of no-slip near the solid walls is usually applied, where the fluid velocity component is assumed to be zero relative to the solid boundary. However, this is not true for fluid flows at the micro- and nanoscale. Investigations show that the no-slip condition is no longer valid, and instead, a certain degree of tangential slip must be allowed [3]. In words, nanoboundary-layer fluid flows mean nanoscale flows which have many applications in microelectromechanical systems. Because of the microscale dimensions of these devices, the fluid flow behavior deviates significantly from the traditional no-slip flow. In recent years, some interest has been given to the study of this type of flow, and some useful results have been introduced by many authors, see, for example, [412].

In this paper we consider the model proposed by Wang [8] to describe the viscous flow due to a stretching surface with both surface slip and suction (or injection). He considered two geometries situations, namely, the two-dimensional and axisymmetric of a stretching surface. Wang [8] applied a similarity transform to convert the Navier-Stokes equations into a 3rd-order nonlinear ordinary differential equation given by where is a parameter describing the type of stretching, where describes the two-dimensional stretching, while is for axisymmetric stretching. The flow is subjected to the following boundary conditions: where is a nondimensional slip parameter, and when injection from the surface occurs and for suction.

During the past two decades, much effort has been spent on searching for exact solutions of nonlinear equations due to their importance in understanding its phenomena. In order to achieve this goal, various direct methods have been proposed, such as tanh-function [13], Jacobi-elliptic function [14], -expansion [15], exp-function ([1620]), the generalized exp-function [21], -expansion ([22, 23]), generalized -expansion [24], and simplified -expansion [25]. However, a little attention was devoted for their applications to solving nonlinear boundary value problems (BVPs). It should be noted that the current work may be the first to indicate the way of applying the -expansion method to solve BVPs, where the advantages of the current method over some of the other ones mentioned above are clarified later in this paper.

In order to solve the BVP given by (1.1) and (1.2), Van Gorder et al. [11] applied the homotopy analysis method. They have also discussed the effects of the slip parameter and the suction parameter on the fluid velocity and on the tangential stress. As expected, they found that for such fluid flows at nanoscales, the shear stress at the wall decreases (in an absolute sense) with an increase in the slip parameter . The existence and uniqueness results for each of the two problems were discussed in [8, 9] along with some numerical results.

It is well known that the exact solutions of nonlinear differential equations are not available in most cases, the reason we sometimes resort to implement accurate numerical methods instead. However, the numerical methods as declared by Rashidi and Erfani [10] gave discontinuous points of a curve, and thus, they are often costly and time consuming to get a complete curve of results. In addition, the stability and convergence of these methods should be considered to avoid divergence or inappropriate results; the numerical method should be therefore chosen carefully.

Chebyshev pseudospectral differentiation matrix (ChPDM) approach has been introduced successfully and applied by Aly et al. [26] 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 have been examined by this technique. Recently, Guedda et al. [27] 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 a water at 4°C (maximum density) and an 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. Very recently, Aly et al. [28] have investigated the effect of magnetic field on viscoelastic fluid flow in boundary-layer through porous media by applying the ChPDM, where the resulting equations for the similar stream function, velocity, and skin friction coefficient were discussed for various parameters. It is found that the results were more accurate than those in the literature.

The motivation of this paper is therefore to extend the applications of both the -expansion method and the ChPDM to solve nonlinear BVPs with nonclassical boundary conditions, where the ChPDM is modified to treat both the infinity and the mixed boundary conditions in a direct manner. With this modification, we are able to obtain very accurate numerical results as will be shown later. The procedure we followed in this paper is found effective in studying the current BVP and may be useful for similar nonlinear problems. The suggested procedure is based first on obtaining all the possible exact solutions together with prescribing the domains of the physical parameters. 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. Moreover, various types of solutions such as the convex, vicinity to zero, and asymptotic are also obtained.

#### 2. Previous Results

In this section, we report some previous results obtained for (1.1) and (1.2). At and , Crane [29] gave the exact solution At arbitrary values of and , Wang [8] obtained a solution in the following form: where is the positive root of the cubic equation When there is no suction, (2.2) reduces to that of Andersson [30]. Moreover, when there is no slip, it reduces to that of P. S. Gupta and A. S. Gupta [31]. Finally, Crane's solution [29] was recovered when both suction and slip are absent.

#### 3. Validation of ChPDM Technique

A new numerical technique, namely, Chebyshev pseudospectral differentiation matrix (ChPDM), introduced by Aly et al. [26, 28] and Guedda et al. [27], has been briefly introduced in this section. 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 kth 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 compared to the other elements. It is shown that the error in is of order , where is the machine precision.

As shown in [2628], on applying the new ChPDM approach, the derivatives of the function at the points are given by Therefore, (1.1) and (1.2) become respectively.

Before starting the current analysis, ChPDM approach is therefore applied by using the system (3.7). Figures 1(a) and 1(b) show the comparison between ChPDM solutions and homotopy analysis method [11] over the current problem for (a) and (b) for various values of the investigated parameters, , , and . As shown, these figures are exactly the same as Figures  1(a) and  1(b) given in [11]. In addition, Figures 2(a) and 2(b) show the results of applying ChPDM approach at and , respectively, in injection () and suction () cases. These figures present exactly the same as Figures  2(c),  3(c),  2(d), and  3(d), respectively, in [11] and as Figures  2(b) and  4(b), respectively, for , in [10]. Hence, without any hesitation, ChPDM technique may be applied with highly trust in the next sections.

Figure 1: ChPDM solutions (solid line) and homotopy analysis method (circles) [11] for (a) and (b) for various values of the investigated parameters.
Figure 2: ChPDM solutions when at (a) (exactly same as Figures  2(c) and  3(c) in [11]) and (b) (exactly same as Figures  2(d) and  3(d) in [11]) for various values of , in injection and suction cases.

#### 4. The Generalized-Expansion Method

In the next subsections, we introduce the basic concept of the generalized -function method. It is then applied to solve the BVP given by (1.1) and (1.2).

##### 4.1. Description of the Method

Consider a given nonlinear ordinary differential equation The generalized -expansion method is then based on the assumption that the exact solution can be expressed in the following form: where and satisfy the following second-order linear ODE: where and are constants to be determined. Degree of the polynomial can be determined by balancing the highest-order derivative with the highest nonlinear terms. Substituting (4.2) into (4.1), using the second-order linear ODE (4.3), and then equating each coefficient of the resulted polynomial to zero yield a set of algebraic equations with respect to , , and . On solving this algebraic system, we may evaluate the values of unknowns. In addition, the solutions of (4.3) depend on whether : where , , and are constants.

##### 4.2. Application to the Problem

As mentioned before, we discuss here the applicability of the generalized -expansion method to solve (1.1) and (1.2). On using the ansatz (4.2) and considering the homogeneous balance between and or in (1.1), we get , so the solution can be supposed in the form By substituting (4.5) and (4.3) into (1.1) and collecting all terms with the same power of and together, the left-hand side of (4.3) can be converted into another polynomial in and . Equating each coefficient of this polynomial to zero, we obtain the following equations: Solving the above algebraic equations by using MATHEMATICA 6 yields the following cases:

(1)

(2)

(3) This is valid in the case of , where the expression under the root of should be a negative to avoid the singularity of the trigonometric function :

(4)

(5)

Remark 4.1. It should be noted that Fang et al. [32] studied recently the viscous flow over a shrinking sheet with a second-order slip flow model. It is observed in their paper that the full Navier-Stocks equations reduce to the same nonlinear ordinary differential equation given by (1.1) in the present paper at . However, boundary conditions are slightly different at the second one and given in [32] as In fact, we can easily deduce the exact solution of the model in [32] by using some of our results as follows. In particular, the exact solution of equation (7) in Fang et al. [32] is already obtained in the current paper and given by (4.7) as The third boundary condition at infinity is satisfied if is replaced by , where . Therefore, the exact solution takes the form where and are unknown parameters determined by applying the boundary conditions (4.12). So, and are obtained by solving two resulted algebraic equations. Therefore, by substituting their values into (4.14), we get the same results given by equations (11) and (12) in [32]. Here, it may be important to refer to the importance of the present result in (4.7) which can be used to obtain many exact solutions for different models with different boundary conditions. We just apply the boundary conditions to (4.7), and the exact solution is immediately resulted.
The advantages of the generalized -expansion method over some of the other methods are clarified at the end of the next section.

#### 5. Exact and Numerical Solutions

It is important here to mention that the exact solutions in (4.8) and (4.9) are both in the final form as they satisfied all the boundary conditions. For the other cases (4.7), (4.10), and (4.11), the boundary conditions are applied to construct more possible exact analytical and numerical solutions of the current boundary value problem. On the spirit of [9], where was examined for any , and to study the possibility of a general solution for the BVP system (1.1) and (1.2), it should be noted that the values of , , and have been considered as any parameters in the following subsections.

Case 1 (at ). Applying the first two boundary conditions to the solution given by (4.7), it then follows that By solving (5.1) for and then substituting it into (4.7) and (5.2), respectively, we obtain the closed form solution where satisfies the cubic equation From (5.3), it should be noted that the third boundary condition at infinity is satisfied providing that is negative. So, has to be a negative root of (5.4). Alternatively, the solution can be written as where is a positive root of the cubic equation The solution given by (5.5) and (5.6) is the same solution reported in [8]. Figures 3(a) and 3(b) show comparison between the solution of ChPDM and -expansion techniques for plane at for (a) , and (b) , . It can be seen from this figure that the solutions are identical and matching with those results in [11].

Figure 3: Profiles of the similar stream function as a function of at for (a) and (b) , solid and dashed lines are for ChPDM and -expansion solutions, respectively.

Convex Solution. In the current value of and at , we set where , then is a solution of (1.1) and (1.2) if and only if satisfies We are now looking for functions such that both hand sides of (5.7) vanish. For , we get that the functions given by [33] with are convex solutions of (1.1) and (1.2).

Case 2 (at ). In this case, we obtain the following closed-form solution: where is a real root of the following cubic equation: Figures 4(a) and 4(b) show comparison between the solution of ChPDM and -expansion techniques for plane at for (a) , and (b) , . In this figure, the solutions are very closed for in the range and slightly closed in the range . The solution of the present case is physically acceptable at a certain range for , which is discussed in the following lemma.

Figure 4: Profiles of the similar stream function as a function of at for (a) and (b) , solid and dashed lines are for ChPDM and -expansion solutions, respectively.

Lemma 5.1. The given solution by (5.10) and (5.11) is valid and unique only if .

Proof. To get the range of at which the solution (5.10) is unique, should satisfy the condition of obtaining only one real root for the cubic equation , which is well known as Inserting , , and into (5.12) yields On noting that , it follows from (5.13) that , and this completes the proof.

Case 3 (at ). (1) When .
As mentioned before, the solution of this case is given in the final form by (4.9). However, it should be noted that and are to be chosen so that is complex to avoid the singularities of . For more specification, we can write the solution in the following simplest form: Therefore, and have to be specified such that Figure 5 presents then existence domain of (5.14) and (5.15) in the parameter plane in the ranges and for and , respectively.
It is observed from the solution given above that it posses singular points. These singularities can be obtained by For given values of and , we can find the singular points of the solution. However, must be continuous in the domain of definition of the current BVP. This means that the solution at some given values of the parameters and is physically acceptable if the singularities lie outside the physical domain, that is, should be negative.
Figures 6(a)6(e) and 7(a)7(d) present profiles of the similar stream function as a function of at the case under consideration for some values of in the regions and , respectively. These figures show that the ChPDM and -expansion solutions are identical inside the domains with slightly difference as near the boundaries of and . However, it is found that for any chosen values of and in , the singular points lie within the domain of interest , and hence, the solution is discontinuous in the physical domain. Table 1 shows singularities in for some values of and . It should be also noted that the behaviour of the solutions in , as becomes negative, inverses those in , as shown in Figure 8. However, more investigations regarding these regions may be taken into account as a future work. (2) When .
Using the first two boundary conditions to the solution given by (4.10) yields By solving (5.17) for , we get By substituting from (5.19) into (5.18) and solving the resulted equation for , we obtain Therefore, the exact solution in this case, , is given by in (4.10), where and are defined by (5.19) and (5.20), respectively. This solution can be easily checked by a direct substitution.

Table 1: Some values of and with their singularities in .
Figure 5: Existence four domains of (5.14) and (5.15) in the parameter plane at and .
Figure 6: Profiles of the similar stream function as a function of at in for (a) and , (b) and , (c) and , (d) and , and (e) and , where solid and dashed lines are for ChPDM and -expansion solutions, respectively.
Figure 7: Profiles of the similar stream function as a function of at in for (a) and , (b) and , (c) and , and (d) and , where solid and dashed lines are for ChPDM and -expansion solutions, respectively.
Figure 8: Profiles of the similar stream function as a function of at in , for and , and , for and , where solid and dashed lines are for ChPDM and -expansion solutions, respectively.

Case 4 (a fractional ). Proceeding as above, we obtain the following equations from applying the first two boundary conditions to the solution in (4.11): By solving (5.21) for and substituting the result into in (4.11) and (5.22), then the following exact solution can be derived: provided that , , and are related by the equation The solution given by (5.23) and (5.24) can be also checked by a direct substitution.

Lemma 5.2. From (5.24), it is observed that at the following cases:(i)at when ,(ii)at when ,(iii)at when .

Proof. On solving (5.24) for , we then get (i)At , this leads to . Now setting , we obtain , that is, . However, and .(ii)When , this means that . In this case, the denominator in (5.25) should be negative so that . Hence, must be negative real number, that is, .(iii)Here, we may rewrite in (5.25) as At , we note that . Therefore, is positive if the numerator and denominator have the same sign, that is, the following conditions have to be held: On solving the inequalities in (5.27), we obtain the range of as

##### 5.1. Important Remark:-Expansion Method Over the Others

Here, we aim to show that the exact solutions obtained in the previous subsections cannot be achieved by using some other methods. Ebaid [18] pointed out that the Jacobi-elliptic function and -expansion methods cannot be used to search the exact solutions for nonlinear differential equations that include both odd and even-order derivative terms. However, He indicated the applicability of the standard exp-function method to solve such kind of equations. Despite this ability, on applying the exp-function method [17] to solve (1.1) and (1.2), it does not provide any of the exact solutions obtained by the current method. To clarify this point, this method is applied to the present BVP when , and the results are found as where and are governed by the following equation: In view of (5.29), the solution is trivial at . In the case of , this exact solution is physically acceptable in a very short range for the parameter (). The solution given by (5.29) and (5.30) is therefore just a special case of that obtained by generalized -expansion method and given by (5.5) and (5.6). The same conclusion can be also deduced when applying the standard exp-function method at . Moreover, the rest of solutions given by (5.10) and (5.11) (at ) and (5.23) and (5.24) cannot be recovered by using the exp-function method. In view of this discussion, it may be concluded that the generalized -expansion method has many advantages over the exp-function method for the present BVP. Further, in order to make this point as clear as possible, an appendix containing the mathematical details of applying the exp-function method to the current problem is added.

#### 6. Solution for

In the view of Aly et al. [34], the main attention is paid to the construction of solutions for (1.1) and (1.2), when . We look for a solution which has the form Substituting (6.1) into (1.1) and (1.2), we obtain the equations We obtain at where at , we get Equation (6.3) appears in the study of similarity solutions to problems of boundary-layer theory in some contexts of fluid mechanics, see, for example, [26, 27, 33, 34]. Exact solutions for these equations can be easily found, for different values of , by applying the technique in Section 5. We substitute then this solution in the system (6.4) to get the construction of . This means that the form of can be therefore evaluated by (6.1).

#### 7. Asymptotic Solution ()

As in [26] and by means of a shooting method, the boundary condition at infinity is replaced by the condition where is the shooting parameter which has to be determined. Regarding the difficulties of obtaining the numerical solution of the system (1.1) and (1.2) in the case of and as in [26, 27, 34], we now seek a new set of full equations which do not contain on using the following transformation: where ; and are constants to be determined. On substituting expressions (7.2) into (1.1), we obtain In order to ensure that the highest derivative remains present in the resulting equation, so avoiding the need to disregard any of the boundary conditions, we look for a balance within the equation of this term. Hence, we obtain . Therefore, when (i.e., ), we obtain the following new set of full equations: The boundary conditions (7.5) do not contain . As in the last section, (7.4) and (7.5) can be solved by meaning of Section 5.

#### 8. Conclusion

Third-order nonlinear differential equations describing the nano boundary-layer flow have been investigated theoretically, using -expansion method, and numerically, applying ChPDM approach. The present results are itemized as follows.(i)It is found that the ChPDM results are very accurate in an excellent manner on comparing to those published in the literature using the homotopy analysis method and the modified differential transform-Padé method.(ii)For the first time, we have showed the way of applying -expansion method to solve nonlinear BVPs. In addition, it has been proven that it has many advantages over some of the other methods on solving the present BVP, where four certain domains for the physical parameters have been discussed.(iii)ChPDM technique has been successfully applied to validate and evidence the resulted exact solutions, for different positive and negative values of the investigated parameters, , , and .(iv)Convex solutions have been obtained at in the special case of .(v)Vicinity of zero and asymptotic solutions when and , respectively, are also deduced. It should be noted that, in both cases, does not exist. Therefore, the resulting equations are easy to deal with analytically and numerically.

#### Appendix

In this section, details of applying the exp-function method to the investigated problem (1.1) and (1.2) are introduced. The aim is to confirm that the solution obtained through this method given by (5.26) and (5.27) is in fact a special case of the exact one obtained by using the -expansion method. Very recently, Ebaid [35] proved that on searching for exact solution by using the exp-function method, one can go directly by assuming the solution in the form: where the tedious calculations of the balancing procedure are not required. On substituting (A.1) into (1.1), multiplying by , and then equating the coefficients of each exp-function to zero, we obtain the following system of algebraic equations: On solving this system, we obtain a nontrivial solution at as On applying the first boundary condition, we obtain . The exact solution hence becomes which is equivalent to (5.29). Now on applying the second boundary condition, this leads directly to (5.30).

In conclusion, the above discussion shows that application of the exp-function method to the present problem gives the same solution in (5.29) and (5.30) which is already shown in Section 5.1. Exp-function method came therefore as a special case from one of the solutions provided by the -expansion method.

#### Acknowledgment

This paper was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. (105-130-D1432). The authors, therefore, acknowledge with thanks DSR technical and financial support.

#### References

1. J.-H. He, Y.-Q. Wan, and L. Xu, “Nano-effects, quantum-like properties in electrospun nanofibers,” Chaos, Solitons and Fractals, vol. 33, no. 1, pp. 26–37, 2007.
2. L. Prandtl, “Über Flüssigkeitsbewegung bei sehr kleiner Reibung,” in Proceedings of the 3rd International Mathematical Congress, 1904.
3. M. Gad-el-Hak, “The fluid mechanics of macrodevices-the Freeman scholar lecture,” Journal of Fluids Engineering, vol. 121, pp. 5–33, 1999.
4. M. T. Matthews and J. M. Hill, “Nano boundary layer equation with nonlinear Navier boundary condition,” Journal of Mathematical Analysis and Applications, vol. 333, no. 1, pp. 381–400, 2007.
5. M. T. Matthews and J. M. Hill, “Micro/nano thermal boundary layer equations with slip-creep-jump boundary conditions,” IMA Journal of Applied Mathematics, vol. 72, no. 6, pp. 894–911, 2007.
6. M. T. Matthews and J. M. Hill, “A note on the boundary layer equations with linear slip boundary condition,” Applied Mathematics Letters, vol. 21, no. 8, pp. 810–813, 2008.
7. J. Cheng, S. Liao, R. N. Mohapatra, and K. Vajravelu, “Series solutions of nano boundary layer flows by means of the homotopy analysis method,” Journal of Mathematical Analysis and Applications, vol. 343, no. 1, pp. 233–245, 2008.
8. C. Y. Wang, “Analysis of viscous flow due to a stretching sheet with surface slip and suction,” Nonlinear Analysis. Real World Applications, vol. 10, no. 1, pp. 375–380, 2009.
9. F. Talay Akyildiz, H. Bellout, K. Vajravelu, and R. A. Van Gorder, “Existence results for third order nonlinear boundary value problems arising in nano boundary layer fluid flows over stretching surfaces,” Nonlinear Analysis. Real World Applications, vol. 12, no. 6, pp. 2919–2930, 2011.
10. M. M. Rashidi and E. Erfani, “The modified differential transform method for investigating nano boundary-layers over stretching surfaces,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 21, pp. 864–883, 2011.
11. 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.
12. C. Y. Wang, “Flow due to a stretching boundary with partial slip-an exact solution of the Navier-Stokes equations,” Chemical Engineering Science, vol. 57, pp. 3745–3747, 2002.
13. E. J. Parkes and B. R. Duffy, “An automated tanh-function method for finding solitary wave solutions to nonlinear evolution equations,” Computer Physics Communications, vol. 98, pp. 288–300, 1996.
14. S. Liu, Z. Fu, S. Liu, and Q. Zhao, “Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations,” Physics Letters A, vol. 289, no. 1-2, pp. 69–74, 2001.
15. A. Ebaid and E. H. Aly, “Exact solutions for the transformed reduced Ostrovsky equation via the F-expansion method in terms of Weierstrass-elliptic and Jacobian-elliptic functions,” Wave Motion, vol. 49, pp. 296–308, 2012.
16. J.-H. He and X.-H. Wu, “Exp-function method for nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 30, no. 3, pp. 700–708, 2006.
17. X.-H. Wu and J.-H. He, “EXP-function method and its application to nonlinear equations,” Chaos, Solitons and Fractals, vol. 38, no. 3, pp. 903–910, 2008.
18. A. Ebaid, “Exact solitary wave solutions for some nonlinear evolution equations via Exp-function method,” Physics Letters A, vol. 365, no. 3, pp. 213–219, 2007.
19. A. Ebaid, “Application of the exp-function method for solving some evolution equations with nonlinear terms of any orders,” Zeitschrift für Naturforschung, vol. 65, pp. 1039–1044, 2010.
20. C. Chun, “New solitary wave solutions to nonlinear evolution equations by the Exp-function method,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 2107–2110, 2011.
21. A. Ebaid, “Generalization of He’s exp-function method and new exact solutions for Burgers equation,” Zeitschrift für Naturforschung, vol. 64, pp. 604–608, 2009.
22. M. Wang, X. Li, and J. Zhang, “The ${G}^{\prime }/G$-expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics,” Physics Letters A, vol. 372, no. 4, pp. 417–423, 2008.
23. R. Abazari, “Application of ${G}^{\prime }/G$-expansion method to travelling wave solutions of three nonlinear evolution equation,” Computers & Fluids, vol. 39, no. 10, pp. 1957–1963, 2010.
24. S. Yu, “A generalized ${G}^{\prime }/G$-expansion method and its application to the MKdV equation,” International Journal of Nonlinear Science, vol. 8, no. 3, pp. 374–378, 2009.
25. X. Fan, S. Yang, and D. Zhao, “Travelling wave solutions for the Gilson-Pickering equation by using the simplified ${G}^{\prime }/G$-expansion method,” International Journal of Nonlinear Science, vol. 8, no. 3, pp. 368–373, 2009.
26. 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.
27. 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°C,” Applied Mathematical Modelling, vol. 35, no. 10, pp. 5182–5197, 2011.
28. E. H. Aly, N. T. El-Dabe, and A. S. Al-Bareda, “ChPDM analysis for MHD flow of viscoelastic fluid through porous media,” Journal of Applied Sciences Research. In press.
29. L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik, vol. 21, pp. 645–647, 1970.
30. H. I. Andersson, “Slip flow past a stretching surface,” Acta Mechanica, vol. 158, pp. 121–125, 2002.
31. P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction or blowing,” The Canadian Journal of Chemical Engineering, vol. 55, pp. 744–746, 1977.
32. T. Fang, S. Yao, J. Zhang, and A. Aziz, “Viscous flow over a shrinking sheet with a second order slip flow model,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 7, pp. 1831–1842, 2010.
33. B. Brighi and J.-D. Hoernel, “Similarity solutions for high frequency excitation of liquid metal in an antisymmetric magnetic field,” in Self-Similar Solutions of Nonlinear PDE, vol. 74, pp. 41–57, Polish Academy of Science, Warsaw, Poland, 2006.
34. E. H. Aly, L. Elliott, and D. B. Ingham, “Mixed convection boundary-layer flow over a vertical surface embedded in a porous medium,” European Journal of Mechanics B, vol. 22, no. 6, pp. 529–543, 2003.
35. A. Ebaid, “An improvement on the exp-function method when balancing the highest order linear and nonlinear terms,” Journal of Mathematical Analysis and Applications, vol. 392, pp. 1–5, 2012.