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

Authors: Emad H. Aly; Abdelhalim Ebaid
Journal: Abstract and Applied Analysis
 (2012)

---

## Abstract

The nonlinear boundary value problem describing the nanoboundary-layer flow with linear
Navier boundary condition is investigated theoretically and numerically in this paper. TheG′/G-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 G′/G-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 G′/G-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.


---

## Body

## 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, [4–12].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
(1.1)f′′′(η)+mf(η)f′′(η)-[f′(η)]2=0,
where m is a parameter describing the type of stretching, where m=1 describes the two-dimensional stretching, while m=2 is for axisymmetric stretching. The flow is subjected to the following boundary conditions:
(1.2)f(0)=s,f'(0)=1+λf′′(0),f′(∞)=0,
where λ>0 is a nondimensional slip parameter, and s<0 when injection from the surface occurs and s>0 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], F-expansion [15], exp-function ([16–20]), the generalized exp-function [21], G'/G-expansion ([22, 23]), generalized G'/G-expansion [24], and simplified G'/G-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 G'/G-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 s 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 theG′/G-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 m=1 and s=λ=0, Crane [29] gave the exact solution
(2.1)f(η)=1-e-η.
At arbitrary values of λ and s, Wang [8] obtained a solution in the following form:
(2.2)f(η)=γ-(γ-s)e-γη,
where γ is the positive root of the cubic equation
(2.3)λγ3+(1-λs)γ2-sγ-1=0.
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 [0,η∞], then the following algebraic mapping:
(3.1)z=2ηη∞-1
transfers the domain to the Chebyshev one, that is, [-1,1]. It is known that the Chebyshev polynomials are usually taken with their associated collocation points in the interval [-1,1] given by
(3.2)zj=cos(πNj),j=0,1,…,N.
Therefore, the kth derivative of any function, say F(z), at these collocation points can be approximated by the equation
(3.3)F(k)=D(k)F,
where D(k)F is the Chebyshev pseudospectral approximation of F(k) where F=[F(z0),F(z1),…,F(zN)]T and F(k)=[F(k)(z0),F(k)(z1),…,F(k)(zN)]T. The entries of the matrix D(k) are given by
(3.4)di,j(k)=2θjN∑r=kN‍∑n=0(n+r-k)evenr-k‍θrbn,rk(-1)[(rj+ni)/N]zrj-N[rj/N]zni-N[ni/N],
where θj=1, except for θ0=θN=1/2 and
(3.5)bn,rk=2kr(k-1)!cn(v-n+k-1)!(v+k-1)!(v)!(v-n)!,
where 2v=r+n-k and c0=2, cj=1, j≥1. The elements d0,1(k) 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 d0,1(1) is of order O(N2εr), where εr is the machine precision.As shown in [26–28], on applying the new ChPDM approach, the derivatives of the function f(η) at the points zi are given by
(3.6)f(k)(zi)=∑j=0N‍di,j(k)f(zj),k=1,2,3,i=1,2,…,N.
Therefore, (1.1) and (1.2) become
(3.7)∑j=0N‍di,j(3)f(zj)+mf(zi)(η∞2)∑j=0N‍di,j(2)f(zj)-(η∞2)(∑j=0N‍di,j(1)f(zj))2=0,f(zN)=s,(η∞2)∑j=0N‍dN,j(1)f(zj)=(η∞2)2+λ∑j=0N‍dN,j(2)f(zj),∑j=0N‍d0,j(1)f(zj)=0,
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) f(η) and (b) f'(η) for various values of the investigated parameters, m, λ, and s. 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 λ=0.5 and λ=5, respectively, in injection (s<0) and suction (s>0) 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 λ=0.5, in [10]. Hence, without any hesitation, ChPDM technique may be applied with highly trust in the next sections.ChPDM solutions (solid line) and homotopy analysis method (circles) [11] for (a) f(η) and (b) f'(η) for various values of the investigated parameters.
(a)
(b)ChPDM solutions whenm=2 at (a) λ=0.5 (exactly same as Figures  2(c) and  3(c) in [11]) and (b) λ=5 (exactly same as Figures  2(d) and  3(d) in [11]) for various values of s, in injection and suction cases.
(a)
(b)
(c)
(d)


## 4.The GeneralizedG′/G-Expansion Method

In the next subsections, we introduce the basic concept of the generalizedG′/G-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(4.1)N(f,dfdη,d2fdη2,d3fdη3,…)=0.
The generalized G′/G-expansion method is then based on the assumption that the exact solution can be expressed in the following form:
(4.2)f(η)=∑i=-nn‍ai(G′G)i,
where ai≠0 and G=G(η) satisfy the following second-order linear ODE:
(4.3)G′′+σG′+μG=0,
where σ and μ are constants to be determined. Degree of the polynomial n 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 ai, σ, and μ. On solving this algebraic system, we may evaluate the values of unknowns. In addition, the solutions of (4.3) depend on whether σ2-4μ(>,<or=)0:
(4.4)G'G={σ2-4μ2tanh(σ2-4μ2η+η0)-σ2,σ2-4μ>0,|tanh(η0)|>1σ2-4μ2coth(σ2-4μ2η+η0)-σ2,σ2-4μ>0,|coth(η0)|<14μ-σ22cot(4μ-σ22η+η0)-σ2,σ2-4μ<0C2C1+C2η-σ2,σ2-4μ=01-1/σ+beση,μ=0μ[C3cos(ημ)-sin(ημ)]cos(ημ)+C3sin(ημ),σ=0,
where C1, C2, and C3 are constants.


### 4.2. Application to the Problem

As mentioned before, we discuss here the applicability of the generalizedG′/G-expansion method to solve (1.1) and (1.2). On using the ansatz (4.2) and considering the homogeneous balance between f′′′(η) and f(η)f′′(η) or [f′(η)]2 in (1.1), we get n=1, so the solution can be supposed in the form
(4.5)f(η)=a0+a1(G′G)+a-1(G′G)-1.
By substituting (4.5) and (4.3) into (1.1) and collecting all terms with the same power of G′/G and G/G′ together, the left-hand side of (4.3) can be converted into another polynomial in G'/G and G/G'. Equating each coefficient of this polynomial to zero, we obtain the following equations:
(4.6)a-1[(2m-1)a-1+6μ]=0,μa-1[(3m-2)σa1+2mμa0+12μσ]=0,a-1[(m-1)(σ2+2μ)a-1+μ(2(m+1)μa1+3mσa0+8μ+7σ2)]=0,a-1[8σμ+σ3+(m-2)σa-1+m(σ2+2μ)a0+4(m+1)σμa1]=0,a-1[mσa0+(σ2+2μ)(1+2(m+1)a1)]-μa1(σ2+2μ-mσa0+μa1)-a-12=0,a1[4(m+1)σa-1-(σ2+8μ)σ+m(σ2+2μ)a0+(m-2)σμa1]=0,a1[2(m+1)a-1+(m-1)(σ2+2μ)a1+3mσa0-7σ2-8μ]=0,a1[2ma0+σ((3m-2)a1-12)]=0,a1[(2m-1)a1-6]=0.
Solving the above algebraic equations by using MATHEMATICA 6 yields the following cases:(1)(4.7)m=1,μ=0,a1=0,a-1=σ(σ+a0),f(η)=-σ+σ(σ+a0)beση,(2)(4.8)m=0,a-1=0,a1=-6,σ2=4μ,f(η)=s+3η8γ2+2γη,where16γ2-6γ-3λ=0,(3)(4.9)m=-1,a-1=2μ,a0=0,a1=-2,Γ=±s2-s3λ2+s2(-1+sλ),f(η)=sΓ(Γ+tan((s/2Γ)η)1-Γtan((s/2Γ)η)).
This is valid in the case of σ=0, where the expression under the root of Γ should be a negative to avoid the singularity of the trigonometric function tan((s/2Γ)η):(4)(4.10)m=-1,a0=σ,a1=0,a-1=2μ,f(η)=σ-4μσ-σ2-4μtanh[((σ2-4μ)/2)η+η0],(5)(4.11)σ2=4μ,a0=3σ1-2m,a1=0,a-1=6μ1-2m,f(η)=6σC2(2m-1)(C1σ-2C2+σC2η).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 m=1. However, boundary conditions are slightly different at the second one and given in [32] as
(4.12)f(0)=s,f'(0)=-1+λf′′(0)+δf′′′(0),f'(∞)=0.
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
(4.13)f(η)=-σ+σ(σ+a0)beση.
The third boundary condition at infinity is satisfied if σ is replaced by -β, where β>0. Therefore, the exact solution takes the form
(4.14)f(η)=-β+β(β-a0)be-βη,
where a0 and b are unknown parameters determined by applying the boundary conditions (4.12). So, a0 and b 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 generalizedG'/G-expansion method over some of the other methods are clarified at the end of the next section.


## 4.1. Description of the Method

Consider a given nonlinear ordinary differential equation(4.1)N(f,dfdη,d2fdη2,d3fdη3,…)=0.
The generalized G′/G-expansion method is then based on the assumption that the exact solution can be expressed in the following form:
(4.2)f(η)=∑i=-nn‍ai(G′G)i,
where ai≠0 and G=G(η) satisfy the following second-order linear ODE:
(4.3)G′′+σG′+μG=0,
where σ and μ are constants to be determined. Degree of the polynomial n 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 ai, σ, and μ. On solving this algebraic system, we may evaluate the values of unknowns. In addition, the solutions of (4.3) depend on whether σ2-4μ(>,<or=)0:
(4.4)G'G={σ2-4μ2tanh(σ2-4μ2η+η0)-σ2,σ2-4μ>0,|tanh(η0)|>1σ2-4μ2coth(σ2-4μ2η+η0)-σ2,σ2-4μ>0,|coth(η0)|<14μ-σ22cot(4μ-σ22η+η0)-σ2,σ2-4μ<0C2C1+C2η-σ2,σ2-4μ=01-1/σ+beση,μ=0μ[C3cos(ημ)-sin(ημ)]cos(ημ)+C3sin(ημ),σ=0,
where C1, C2, and C3 are constants.


## 4.2. Application to the Problem

As mentioned before, we discuss here the applicability of the generalizedG′/G-expansion method to solve (1.1) and (1.2). On using the ansatz (4.2) and considering the homogeneous balance between f′′′(η) and f(η)f′′(η) or [f′(η)]2 in (1.1), we get n=1, so the solution can be supposed in the form
(4.5)f(η)=a0+a1(G′G)+a-1(G′G)-1.
By substituting (4.5) and (4.3) into (1.1) and collecting all terms with the same power of G′/G and G/G′ together, the left-hand side of (4.3) can be converted into another polynomial in G'/G and G/G'. Equating each coefficient of this polynomial to zero, we obtain the following equations:
(4.6)a-1[(2m-1)a-1+6μ]=0,μa-1[(3m-2)σa1+2mμa0+12μσ]=0,a-1[(m-1)(σ2+2μ)a-1+μ(2(m+1)μa1+3mσa0+8μ+7σ2)]=0,a-1[8σμ+σ3+(m-2)σa-1+m(σ2+2μ)a0+4(m+1)σμa1]=0,a-1[mσa0+(σ2+2μ)(1+2(m+1)a1)]-μa1(σ2+2μ-mσa0+μa1)-a-12=0,a1[4(m+1)σa-1-(σ2+8μ)σ+m(σ2+2μ)a0+(m-2)σμa1]=0,a1[2(m+1)a-1+(m-1)(σ2+2μ)a1+3mσa0-7σ2-8μ]=0,a1[2ma0+σ((3m-2)a1-12)]=0,a1[(2m-1)a1-6]=0.
Solving the above algebraic equations by using MATHEMATICA 6 yields the following cases:(1)(4.7)m=1,μ=0,a1=0,a-1=σ(σ+a0),f(η)=-σ+σ(σ+a0)beση,(2)(4.8)m=0,a-1=0,a1=-6,σ2=4μ,f(η)=s+3η8γ2+2γη,where16γ2-6γ-3λ=0,(3)(4.9)m=-1,a-1=2μ,a0=0,a1=-2,Γ=±s2-s3λ2+s2(-1+sλ),f(η)=sΓ(Γ+tan((s/2Γ)η)1-Γtan((s/2Γ)η)).
This is valid in the case of σ=0, where the expression under the root of Γ should be a negative to avoid the singularity of the trigonometric function tan((s/2Γ)η):(4)(4.10)m=-1,a0=σ,a1=0,a-1=2μ,f(η)=σ-4μσ-σ2-4μtanh[((σ2-4μ)/2)η+η0],(5)(4.11)σ2=4μ,a0=3σ1-2m,a1=0,a-1=6μ1-2m,f(η)=6σC2(2m-1)(C1σ-2C2+σC2η).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 m=1. However, boundary conditions are slightly different at the second one and given in [32] as
(4.12)f(0)=s,f'(0)=-1+λf′′(0)+δf′′′(0),f'(∞)=0.
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
(4.13)f(η)=-σ+σ(σ+a0)beση.
The third boundary condition at infinity is satisfied if σ is replaced by -β, where β>0. Therefore, the exact solution takes the form
(4.14)f(η)=-β+β(β-a0)be-βη,
where a0 and b are unknown parameters determined by applying the boundary conditions (4.12). So, a0 and b 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 generalizedG'/G-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 m was examined for any m>0, 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 m, s, and λ have been considered as any parameters in the following subsections.Case 1 (atm=1).
Applying the first two boundary conditions to the solution given by (4.7), it then follows that
(5.1)-σ+σ(σ+a0)b=s,(5.2)λσ3(σ+a0)b-σ2(σ+a0)b+1=0.
By solving (5.1) for a0 and then substituting it into (4.7) and (5.2), respectively, we obtain the closed form solution
(5.3)f(η)=-σ+(σ+s)eση,
where σ satisfies the cubic equation
(5.4)λσ3-σ2(1-λs)-σs+1=0.
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
(5.5)f(η)=σ-(σ-s)e-ση,
where σ is a positive root of the cubic equation
(5.6)λσ3+σ2(1-λs)-σs-1=0.
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 G'/G-expansion techniques for (η,f(η)) plane at m=1 for (a) s=1, λ=0.5 and (b) s=0, λ=1. It can be seen from this figure that the solutions are identical and matching with those results in [11].Profiles of the similar stream functionf(η) as a function of η at m=1 for (a) s=1,λ=0.5 and (b) s=0,λ=1, solid and dashed lines are for ChPDM and G'/G-expansion solutions, respectively.
(a)
(b)Convex Solution. In the current value of m and at λf′′(0)=-2, we set f=g+k where k∈ℝ, then f is a solution of (1.1) and (1.2) if and only if g satisfies
(5.7)g′′′+kg′′=(g′)2-gg′′,(5.8)g(0)=s-k,g'(0)=-1,g(∞)=0.
We are now looking for functions g such that both hand sides of (5.7) vanish. For s≥2, we get that the functions f1,f2:[0,∞)→ℝ given by [33]
(5.9)fi(x)=Ki+1Kie-Kixfori=1,2
with K1,2=(1/2)(s∓s2-4) are convex solutions of (1.1) and (1.2).Case 2 (atm=0).
In this case, we obtain the following closed-form solution:(5.10)f(η)=s+3η8γ2+2γη,
where γ is a real root of the following cubic equation:
(5.11)16γ3-6γ-3λ=0.
Figures 4(a) and 4(b) show comparison between the solution of ChPDM and G'/G-expansion techniques for (η,f(η)) plane at m=0 for (a) s=1, λ=0.5 and (b) s=0, λ=1. In this figure, the solutions are very closed for η in the range [0,20] and slightly closed in the range [20,η∞]. The solution of the present case is physically acceptable at a certain range for λ, which is discussed in the following lemma.Profiles of the similar stream functionf(η) as a function of η at m=0 for (a) s=1,λ=0.5 and (b) s=0,λ=1, solid and dashed lines are for ChPDM and G'/G-expansion solutions, respectively.
(a)
(b)Lemma 5.1.
The given solution by (5.10) and (5.11) is valid and unique only if λ∈]2/3,∞[.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 aγ3+bγ+c=0, which is well known as
(5.12)4b3-27ac2<0.
Inserting a=16, b=-6, and c=-λ into (5.12) yields
(5.13)2-9λ2<0.
On noting that λ>0, it follows from (5.13) that λ>2/3, and this completes the proof.Case 3 (atm=-1).
(1)
Whenσ=0.
As mentioned before, the solution of this case is given in the final form by (4.9). However, it should be noted that λ and s are to be chosen so that Γ is complex to avoid the singularities of tan((s/2Γ)η). For more specification, we can write the solution in the following simplest form:
(5.14)f(η)=sΛ[Λ-tanh((s/2Λ)η)1-Λtanh((s/2Λ)η)],whereΛ=λs3-s2λs3-s2+2,Γ=iΛ.
Therefore, λ and s have to be specified such that
(5.15)s3-s2+2<0,λs3-s2<0,s3-s2+2>0,λs3-s2>0.
Figure 5 presents then existence domain of (5.14) and (5.15) in the parameter plane (λ,s) in the ranges [-5,5] and [-3,3] for λ and s, respectively.
It is observed from the solution given above that it posses singular points. These singularities can be obtained by(5.16)ηsingular=2Λsarctanh(1Λ).
For given values of λ and s, we can find the singular points of the solution. However, f(η) must be continuous in the domain of definition [0,∞[ of the current BVP. This means that the solution at some given values of the parameters λ and s is physically acceptable if the singularities lie outside the physical domain, that is, ηsingular should be negative.
Figures6(a)–6(e) and 7(a)–7(d) present profiles of the similar stream function f(η) as a function of η at the case under consideration for some values of (λ,s) in the regions R1 and R2, respectively. These figures show that the ChPDM and G'/G-expansion solutions are identical inside the domains with slightly difference as η→∞ near the boundaries of R1 and R2. However, it is found that for any chosen values of λ and s in R3, the singular points lie within the domain of interest [0,∞[, and hence, the solution is discontinuous in the physical domain. Table 1 shows singularities in R3 for some values of λ and s. It should be also noted that the behaviour of the solutions in R4, as λ becomes negative, inverses those in R1, as shown in Figure 8. However, more investigations regarding these regions may be taken into account as a future work. (2)
Whenσ≠0.
Using the first two boundary conditions to the solution given by (4.10) yields
(5.17)σ-4μσ-σ2-4μtanh(η0)=s,(5.18)(σcosh(η0)-σ2-4μsinh(η0))3=2μ(4μ-σ2)[(σ+4λμ-λσ2)cosh(η0)+(λσ-1)σ2-4μsinh(η0)].
By solving (5.17) for η0, we get
(5.19)η0=tanh-1(σ(s-σ)+4μ(s-σ)σ2-4μ).
By substituting η0 from (5.19) into (5.18) and solving the resulted equation for σ, we obtain
(5.20)σ=±λs3-s2+2+4μ(λs-1)λs-1.
Therefore, the exact solution in this case, m=-1, is given by f(η) in (4.10), where η0 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 s with their singularities in R3.
λ
s
η
singular
0.2
2
1.16
0.1
3
0.95
−0.1
2
1.32
−0.1
3
1.12
−1
2
1.69
−3
2
2.07
−4
1
2.66Figure 5
Existence four domains of (5.14) and (5.15) in the parameter plane (λ,s) at m=-1 and σ=0.Profiles of the similar stream functionf(η) as a function of η at m=-1,σ=0 in R1 for (a) λ=-0.1 and s=-2.5, (b) λ=0.5 and s=-2, (c) λ=1 and s=-3, (d) λ=2 and s=-2, and (e) λ=5 and s=-1, where solid and dashed lines are for ChPDM and G'/G-expansion solutions, respectively.
(a)
(b)
(c)
(d)
(e)Profiles of the similar stream functionf(η) as a function of η at m=-1,σ=0 in R2 for (a) λ=2 and s=1, (b) λ=3 and s=2, (c) λ=4 and s=1, and (d) λ=5 and s=3, where solid and dashed lines are for ChPDM and G'/G-expansion solutions, respectively.
(a)
(b)
(c)
(d)Figure 8
Profiles of the similar stream functionf(η) as a function of η at m=-1,σ=0 in R1, for λ=2 and s=-2, and R4, for λ=-2 and s=-2, where solid and dashed lines are for ChPDM and G′/G-expansion solutions, respectively.Case 4 (a fractionalm).
Proceeding as above, we obtain the following equations from applying the first two boundary conditions to the solutionf(η) in (4.11):
(5.21)6σC2(2m-1)(C1σ-2C2)=s,(5.22)(2m-1)(C1σ-2C2)3=6σ2C22[σC1+2(λσ-1)C2].
By solving (5.21) for C1 and substituting the result into f(η) in (4.11) and (5.22), then the following exact solution can be derived:
(5.23)f(η)=6s6+(2m-1)sη,
provided that m, λ, and s are related by the equation
(5.24)λ(2m-1)2s3+3(2m-1)s2+18=0.
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 λ>0 at the following cases:(i)
at0<m<1/2 when s∈]-6/(2m-1),∞[,(ii)
atm>1/2 when s<0,(iii)
atm<0 when s∈]6/(1-2m),∞[∪]-6/(1-2m),0[.Proof.
On solving (5.24) for λ, we then get
(5.25)λ=-3(2m-1)s2+18(2m-1)2s3.(i)
At0<m<1/2, this leads to 3(2m-1)s2<0. Now setting 3(2m-1)s2+18<0, we obtain s2>-6/(2m-1), that is, s∈]-∞,--6/(2m-1)[⋃‍]-6/(2m-1),∞[. However, s∈]-∞,--6/(2m-1)[⇒λ<0 and s∈]-6/(2m-1),∞[⇒λ>0.(ii)
Whenm>1/2, this means that 3(2m-1)s2+18>,foralls∈ℜ. In this case, the denominator in (5.25) should be negative so that λ>0. Hence, s must be negative real number, that is, s<0.(iii)
Here, we may rewriteλ in (5.25) as
(5.26)λ=3(1-2m)s2-18(2m-1)2s3. At m=0, we note that 3(1-2m)s2>0. Therefore, λ is positive if the numerator and denominator have the same sign, that is, the following conditions have to be held:
(5.27)3(1-2m)s2-18>0,s>0,3(1-2m)s2-18<0,s<0.
On solving the inequalities in (5.27), we obtain the range of s as
(5.28)s∈]61-2m,∞[⋃‍]-61-2m,0[.


### 5.1. Important Remark:G′/G-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 F-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 m=1, and the results are found as
(5.29)f(η)=1-(1-s)e-η,
where λ and s are governed by the following equation:
(5.30)λ(1-s)-s=0orλ=s1-s,s≠1.
In view of (5.29), the solution is trivial at s=1. In the case of λ>0, this exact solution is physically acceptable in a very short range for the parameter (s∈[0,1[). The solution given by (5.29) and (5.30) is therefore just a special case of that obtained by generalized G'/G-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 m=-1. Moreover, the rest of solutions given by (5.10) and (5.11) (at m=0) 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 G'/G-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.


## 5.1. Important Remark:G′/G-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 F-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 m=1, and the results are found as
(5.29)f(η)=1-(1-s)e-η,
where λ and s are governed by the following equation:
(5.30)λ(1-s)-s=0orλ=s1-s,s≠1.
In view of (5.29), the solution is trivial at s=1. In the case of λ>0, this exact solution is physically acceptable in a very short range for the parameter (s∈[0,1[). The solution given by (5.29) and (5.30) is therefore just a special case of that obtained by generalized G'/G-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 m=-1. Moreover, the rest of solutions given by (5.10) and (5.11) (at m=0) 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 G'/G-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 for0<λ≪1

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 λ→0. We look for a solution which has the form
(6.1)f(η)=f0+λf1+⋯.
Substituting (6.1) into (1.1) and (1.2), we obtain the equations
(6.2)f0′′′+λf1′′′+m(f0+λf1)(f0′′+λf1′′)-(f0′+λf1′)2=0,f0(0)+λf1(0)+𝒪(λ2)=s,f0′(0)+λf1′(0)+𝒪(λ2)=1+λ[f0′′(0)+λf1′′(0)+𝒪(λ2)],f0′(∞)+λf1′(∞)+𝒪(λ2)=0.
We obtain at λ0(6.3)f0′′′+mf0f0′′-(f0′)2=0,f0(0)=s,f0′(0)=1,f0′(∞)=0,
where at λ1, we get
(6.4)f1′′′+m(f0f1′′+f1f0′′)-2f0′f1′=0,f1(0)=0,f1′(0)=f0′′(0),f1′(∞)=0.
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 m, by applying the technique in Section 5. We substitute then this solution in the system (6.4) to get the construction of f1. This means that the form of f(η) can be therefore evaluated by (6.1).


## 7. Asymptotic Solution (λ≫1)

As in [26] and by means of a shooting method, the boundary condition at infinity is replaced by the condition
(7.1)f′′(0)=α,
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 λ≫1 and as in [26, 27, 34], we now seek a new set of full equations which do not contain λ on using the following transformation:
(7.2)f(η)=η+ϵϑH(ζ)whereζ=ϵγη,
where ϵ=λα; ϑ and γ are constants to be determined. On substituting expressions (7.2) into (1.1), we obtain
(7.3)ϵ2γH′′′+ϵγmηH′′+ϵϑ+γ[mHH′′-H'2]-2H′-ϵ-ϑ=0.
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 ϑ=γ=1/2. Therefore, when ϵ→∞ (i.e., λ→∞), we obtain the following new set of full equations:
(7.4)H′′′+mHH′′-H'2=0,(7.5)H(0)=0,H′(0)=1,H′′(0)=0.
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, usingG'/G-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 applyingG'/G-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,m, s, and λ.(iv)
Convex solutions have been obtained atm=1 in the special case of λf′′(0)=-2.(v)
Vicinity of zero and asymptotic solutions when0<λ≪1 and λ≫1, 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.



---

*Source: 415431-2012-09-25.xml*