Abstract
This paper shows the analysis of the thin film flow of fourth-grade fluid on the outer side of a vertical cylinder. Solution of the governing nonlinear equation is obtained by Rational Homotopy Perturbation Method (RHPM); comparison with exact solution reflects the reliability of the method. Analysis shows that this method is reliable for even high nonlinearity. Graphs and tables strengthen the idea.
1. Introduction
Recently, a number of new and modified techniques have been introduced by various scientists which subsequently proved extremely useful to tackle various nonlinear problems of diversified physical nature. It is an established fact that most of the physical problems in nature are nonlinear and hence solutions of such problems are of utmost importance. In a similar context, Rational Homotopy Perturbation Method (RHPM) [1] is a newly developed modified form of Homotopy Perturbation Method (HPM). HPM is developed by the coupling of Homotopy and perturbation is very useful for solving nonlinear problems [2–16]. The literature reveals that Jalaal et al. applied Homotopy Perturbation Method to find the velocity profile of a spherical solid particle in plane Couette fluid flow [2]. Liao compared Homotopy Analysis and Homotopy Perturbation Method [3], and Ghorbani and Nadjafi made an elegant comparison of He’s polynomials which are obtained from Homotopy Perturbation Method (HPM) Adomian’s Polynomials [4]. They worked on new developments of the HPM [5]. Mohyud-Din applied HPM on various problems like partial differential equations, higher dimensional initial boundary value problem, and fourth-order, nonlinear higher-order, and sixth-order boundary value problems [6–10]. Xu solved boundary layer equation in unbounded domain by HPM [11]. Vázquez-Leal [1] modified the idea of Homotopy Perturbation Method. He used two power series of Homotopy parameters in a quotient, resulting in a series of linear differential equations.
Siddiqui applied Homotopy Perturbation Method on thin film flow of fourth-grade fluid on the outer wall of vertical cylinder. Hayat applied Homotopy Analysis Method on the same problem and proved convergence. Afterwards, Sajid introduced the slip effect for this problem and presented solution using Homotopy Analysis Method [12–16].
It is worth mentioning that Hayat et al. [17, 18] made an appropriate use of Homotopy technique and applied the same on thin film flow of an Oldroyd 6-constant fluid over a moving belt [17] and also on Couette and Poiseuille flows for fourth-grade fluid [18]. It is to be highlighted that the proposed method (RHPM) is highly suitable for a wide range of physical problems including heat transfer and thermodynamics. The detailed study of the literature reconfirms the reliability of the suggested Rational Homotopy Perturbation Method (see [19–23] and the references therein). In a similar context, Ganji et al. [20, 21] solved nonlinear equations arising in heat transfer by applying the coupling of Homotopy and perturbation. Moreover, Chowdhury and Hashim [22] and Islam et al. [23] extended modified version of Homotopy Perturbation Method for problems related to heat transfer. Inspired and motivated by the ongoing research in this area, we apply a relatively new modified version of Homotopy Perturbation Method which is called Rational Homotopy Perturbation Method (RHPM) [24] on the thin film flow of fourth-grade fluid coupled with slip effect. A convergent and reliable solution is presented with the aid of figures and tables. Analysis shows that Rational Homotopy Perturbation Method is adequate to present the flow behavior of this problem. It is observed that RHPM is very efficient for such problems. Moreover, RHPM is equally applicable on the mathematical models derived from nature like Stiff system of equations [24], transient of nonlinear circuits [19], and heat transfer problems (see [19–24] and the references therein).
2. Basic Idea of RHPM
Foundation of Rational Homotopy Perturbation Method and Homotopy Perturbation Method is the same. Consider the nonlinear differential equation to explain both methods: and consider the boundary condition where and are nonlinear and linear operators, respectively, and is a boundary operator. is an analytic function and domain has the boundary . Homotopy can be written as where is an initial approximation for (1) which satisfies the boundary condition and Homotopy parameter is . For and , (3) can be written as For HPM, let us take the solution for (3) asconsider to get approximate solution Let us take the solution of (3) for RHPM as where known analytic functions have independent variables and unknown functions are to be determined by the Rational Homotopy Perturbation Method. An approximate solution of (1) is obtained when limiting case of (7) is taken as . Hence, The limiting case in the above equation is associated with existence of limits
3. Solution by RHPM
Let us consider the dimensionless equation representing thin film flow of fourth-grade fluid which is falling on the outer side of a vertical cylinder having an infinite length of radius [13]:where is nondimensional fluid parameter, is dimensionless constant which corresponds to gravity, and corresponds to partial slip effect. Boundary condition is Exact solution is where Integral in (12) can be calculated numerically.
Consider the following Homotopy equation:where is initial approximation. Let us consider the solution of (14) with order of approximation as where and are the adjustment parameters. By considering (14) and (15) and regrouping according to the -powers, we get the following.
Zeroth-Order Problem. The zeroth-order problem is as follows:subject to the boundary conditions
First-Order Problem. The first-order problem is as follows:and relevant boundary conditions are
Second-Order Problem. The second-order problem is as follows:and relevant boundary conditions are Solving these problems with the corresponding boundary conditions, we have Substituting (22) into (15) and assuming result in Now, let us discuss the solution for different parameters.
Case 1 (, , , and ). Choosing the adjustment parameters for (23) as and , we have
Case 2 (, , , and . Choosing the adjustment parameters for (23) as and , we have
4. Conclusion
In this section, we will give analysis of the present attempt. Figure 1 shows that solution found by Rational Homotopy Perturbation Method overlaps with exact solution. Figure 2 represents the notion that there exists a gradual shift of velocity with varying slip effect. By increasing the slip effect, velocity also increases. Figure 3 shows that, by increasing the nonlinearity and slip effect, Rational Homotopy Perturbation Method and exact solutions are in good agreement. Figure 4 represents the effect of on the velocity. Tables 1-2 show that error is negligible. Table 3 shows that, for different values of slip effects, Rational Homotopy Perturbation Method efficiently represents the flow behavior. Table 4 reflects the reliability of the method for different values of . Tables 3-4 demand that and be determined individually by the method described in Cases 1 and 2. Mathematica is used in all calculations.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This project was funded by the National Plan for Science, Technology and Innovation (MAARIFA), King Abdulaziz City for Science and Technology, Kingdom of Saudi Arabia, Award no. 15-MAT4688-02.