The (G′/G)-expansion method is used for the first time to find traveling wave solutions for thin film equations, where it is found that the related balance numbers are not the usual positive integers. The closed-form solution obtained via this method is in good agreement with the previously obtained solutions of other researchers. It is also noted that, for appropriate parameters, new solitary waves solutions are found.
1. Introduction
In this paper we are interested in the so-called standard thin film equation of the form
(1)ut=-(unuxxx)x,
which in general has important applications in geology, biophysics, physics, and engineering (see [1–3]). Also known as the lubrication equation [4], it models the spreading motion of the free surface of a thin film on a solid substrate [5]. In particular, the function u(x,t) is the thickness of the fluid film at position x and time t. Here the parameter n denotes the kind of flow. In the case n=1, the equation models the thickness of a thin film in a Hele-Shaw cell [6]. When n=2, it is the Navier slip thin film equation which arises in the study of wetting films with a free contact line between film and substrate [7]. Furthermore, when n=3, it corresponds to the surface-tension-driven spreading of a thin Newtonian fluid [8].
King [9] introduced a generalization of the above equation given by fourth-order nonlinear degenerate parabolic equations of the following form:
(2)ut=-(unuxxx+αu(n-1)uxuxx+βu(n-2)(ux)3)x,
where n, α, and β are constants, while the second is a doubly nonlinear equation:
(3)ut=-(un|uxxx|k-1uxxx)x,
where k>0 is a constant related to the flow. Note that when α=0, β=0, and k=1 each of the generalized thin film equations turned into the standard thin film equation (1). One of the most effective direct methods to build traveling wave solution of nonlinear PDEs is the (G′/G)-expansion method, which was first proposed by Wang et al. [10]. It is assumed that the traveling wave solutions can be expressed by a polynomial in (G′/G), where G=G(ζ) satisfies the following second-order linear ordinary differential equation G′′(ζ)+λG′(ζ)+μG(ζ)=0, where ζ=x-ct, λ, μ, and c are constants. Until now, (G′/G)-expansion method has been successfully applied to obtain exact solution for a variety of nonlinear PDEs [11–21].
Our main objective in this paper is to apply the (G′/G) method to provide closed-form travelling wave solutions of the generalized thin film equations and also standard thin film equation. To the best of our knowledge, this is the first time this method has been applied to such equations. In solving these equations, we found an instance where the related balance numbers are not the usual positive integers (see Zhang [22]). It is also noted that for appropriate parameters new solitary waves solutions are found. We compare our solutions with the solutions previously obtained by Bertozzi and Pugh [23] and King [9], where they proved the existence solution to thin film equation via separation of variables. The closed-form solution obtained via this method is in good agreement with the solutions reported in [9, 23].
Our paper is organized as follows: in Section 2, we present the summary of the (G′/G)-expansion method, in Section 3, we describe the applications of the (G′/G)-expansion method for two generalization thin film equations, standard thin film equation and a special case, and in Section 4, some conclusions are given.
2. Summary of the (G′/G)-Expansion Method
In this section, we describe the (G′/G)-expansion method for finding traveling wave solutions of nonlinear partial differential equations (PDEs). Suppose that a nonlinear partial differential equation, say in two independent variables x and t, is given by the following:
(4)P(u,ut,ux,uxt,utt,uxx,…)=0,
where u=u(x,t) is an unknown function and P is a polynomial in u=u(x,t) and its various partial derivatives, in which highest-order derivatives and nonlinear terms are involved. The procedure of the (G′/G)-expansion method can be presented in the following six steps.
Step 1.
To find the traveling wave solutions of (4), we introduce the wave variable
(5)u(x,t)=u(ζ),ζ=x-ct,
where the constant c is generally termed the wave velocity. Substituting (5) into (4), we obtain the following ordinary differential equations (ODE) in ζ (which illustrates a principal advantage of a traveling wave solution, i.e., a PDE is reduced to an ODE) as follows.
(6)P(u,cu′,u′,cu′′,c2u′′,u′′,…)=0.
Step 2.
If necessary we integrate (6) as many times as possible and set the constants of integration to be zero for simplicity. The solution process for (6) is based on the auxiliary conditions that the dependent variable and its first, second, and higher spatial derivatives tend to zero as ζ→∞, that is,
(7)u(ζ→±∞)=0,du(ζ→±∞)dζ=0,d2u(ζ→±∞)dζ2=0,….
From these conditions, we can take the constants of integration to be zero.
Step 3.
We suppose that the solution of nonlinear partial differential equation can be expressed by a polynomial in (G′/G) as follows:
(8)u(ζ)=∑i=0mai(G′G)i,
where G=G(ζ) satisfies the second-order linear ordinary differential equation as follows:
(9)G′′(ζ)+λG′(ζ)+μG(ζ)=0.
Here the prime denotes the derivative respective to ζ, and ai, λ, and μ are real constants with am≠0. Using the general solutions of (9), we have the following:
(10)G′G={-λ2+λ2-4μ2×[c1sinh((λ2-4μ/2)ζ)+c2cosh((λ2-4μ/2)ζ)c1cosh((λ2-4μ/2)ζ)+c2sinh((λ2-4μ/2)ζ)],λ2-4μ>0,-λ2+4μ-λ22×[-c1sin((4μ-λ2/2)ζ)+c2cos((4μ-λ2/2)ζ)c1cos((4μ-λ2/2)ζ)+c2sin((4μ-λ2/2)ζ)],λ2-4μ<0,(c2c1+c2ζ)-λ2,λ2-4μ=0.
Step 4.
The positive integer m can be accomplished by considering the homogeneous balance between the highest order derivatives and nonlinear terms appearing in (6) as follows: if we define the degree of u(ζ) as D[u(ζ)]=m, then the degree of other expressions is defined by
(11)D[ur(dqudζq)s]=mr+s(q+m),
where s is an integer. Therefore, we can get the value of m in (8).
Step 5.
Substitute (8) into (6) and use (9) and collect all terms with the same order of (G′/G) together, then set each coefficient of this polynomial to zero which yields a set of algebraic equations for ai, c, λ, and μ.
Step 6.
Substitute ai, c, λ, and μ obtained in Step 5 and the general solution of (9) into (8). Next, depending on the sign of the discriminant λ2-4μ, we get the solutions of (6). So, we can obtain exact solutions of the given (4).
The advantages of the approach taken in this paper are as follows.
It will be more important to seek solutions of higher-order nonlinear equations which can be reduced to ODEs of the order greater than 3.
In the (G′/G)-expansion method, there is no need to apply the initial and boundary conditions at the outset. The method yields a general solution with free parameters which can be identified by the above conditions.
The general solution obtained by (G′/G)-expansion method without approximation.
Finally, the solution procedure can be easily implemented in Mathematica or Maple.
3. Application of the (G′/G)-Expansion Method3.1. Exact Traveling Wave Solution of Standard Thin Film Equation (1)
Now we consider (1) which arises in the flow of a surface-tension dominated thin liquid film. Substituting (5) into (1) and integrating the result, and for simplicity equating the integration constant equal to zero, we get the following
(12)unu′′′-cu=0.
Suppose that the solution of (12) can be expressed by a polynomial in (G′/G) as follows:
(13)u(ζ)=E(G′G)m,
where E is a real constant to be determined later and G satisfies (9). Balancing between unu′′′ and u, we get m=-3/n. Now it is easy to deduce that
(14)u′=3E(G′/G)(-3+n)/n((G′/G)2+λ(G′/G)+μ)n,u′′=((-3E(G′G)2+λ(G′G)+μ)×((G′G)-3/n(n-3)-3(G′G)(-3+n)/n×λ-(G′G)(-3+2n)/n(3μ+nμ)))×(n2)-1.
With the aid of symbolic computation, substituting (13) along with (9) into (12), and setting the coefficients of all powers of (G′/G) to zero, we obtain the following system of nonlinear algebraic equations for E, c, λ, n, and μ:
(15)(G′G)-3/n:-27E1+nλμ+3E1+nn2λ2μ+27E1+nnλ2μ+81E1+nμ2+54E1+nλ2μ+81E1+nλ2+81E1+nμ+3E1+nn2λ2+27E1+nnμ2+6E1+nμn2+6E1+nμ2n2=0,(G′G)(-3+n)/n:9E1+nn2λμ2+54E1+nλnμ2-54E1+nλn+9E1+nλn2+81E1+nλ+81E1+nλμ2=0,(G′G)-3(1+n)/n:6E1+nn2λμ2+27E1+nλ3+162E1+nλμ=0,(G′G)3(n-1)/n:6E1+nn2-27E1+nn+27E1+n-cE1+nn3=0,(G′G)(-3+2n)/n:27E1+nμ3+6E1+nμ3n2+27E1+nμ3n2+27E1+nμ3n=0.
The solutions of this system are as follows:
(16)λ=0,μ=0,c=(3n-2)(3n-1)3nEn,
where E and n are arbitrary constants.
Consequently, we obtain the exact traveling wave solution of (1),
(17)u(x,t)=u(ζ)=E(c2c1+c2ζ)-3/n,
where ζ=x-ct=x-((3/n)-2)((3/n)-1)(3/n)Ent.
If we set c1=0 and c2=1 in (17), we obtain the solitary wave solution
(18)u(x,t)=u(ζ)=E(ζ)3/n,
where ζ is as above. This is exactly the same solution obtained by Bertozzi and Pugh [23]:
(19)u(x,t)={A(x-ct)3/n,x>ct,0,otherwise,
when c=((3/n)-2)((3/n)-1)(3/n)An.
We remark that if n≥3, the solution of the system using the (G′/G) cannot be solved due to the no-slip boundary condition on the liquid solid surface, in one space dimension, a similar result reached by Bertozzi and Pugh [24].
3.2. Exact Traveling Wave Solution of the Generalized Thin Film Equation (2)
We study nonnegative solutions of the generalized degenerate fourth-order parabolic equation of thin film equation (2). The solution of (2) as found by King [9] is as follows:
(20)u=[-n3c3[(3-n)(3-2n)+3(3-n)α+9β][-(x-ct)3]]1/n,
requiring n>0 and c>0 for 8β<(1-α)2 with 3(α-ν+3)/4<n<3(α+ν+3)/4.
We seek the traveling wave solution of (2) in the form (5). Now upon substituting of (5) into (2), one gets
(21)-cu′+[unu′′′+αu(n-1)u′u′′+βu(n-2)(u′)3]′=0,
and by integrating (21) and, for simplicity, equating the integration constant which is equal to zero, we get
(22)-cu+unu′′′+αu(n-1)u′u′′+βu(n-2)(u′)3=0.
Balancing between unu′′′ and u, we get m=-3/n. Then, suppose that (21) has the following formal solution:
(23)u(ζ)=E(G′G)-3/n,
where E is an unknown constant to be determined later.
Substituting (23), along with (9), into (22), and setting the coefficients of (G′/G)i,(i=0,1,…,5) to zero, we obtain a system of nonlinear algebraic equations as follows:
(24)(G′G)3(n-1)/n:27βE1+n+27E1+n+6E1+nn2-27E1+nn-cEn3-9αE1+n+27αE1+n=0,(G′G)(n-3)/n:3E1+nn2λ2-27E1+nλ2n+6E1+nμn2+6E1+nμ2n2+81E1+nλ2μ-27E1+nnμ+27E1+nnμ2+81E1+nμ+81E1+nμ2+81E1+nλ2+81E1+nλ2μ+3E1+nλ2nμ+9αE1+nμ2n-9αE1+nnμ+81αE1+nλ2μ-9αE1+nλ2nμ+81βE1+nμ+81βE1+nλ2+81βE1+nμ2+81αE1+nμ2+81αE1+nμ+81αE1+nλ2=0,(G′G)-3/n:126E1+nλμ+27E1+nλ3+162βE1+nλμ+6E1+nn2λμ+162αE1+nλμ+27βE1+nλ2+27αE1+nλ3=0,(G′G)(2n-3)/n:81E1+nλμ2-54E1+nλμ+9E1+nλn2+81E1+nλμ2+9E1+nn2λν2+54E1+nnλμ2+81αE1+nλμ2-18αE1+nnλμ2+81E1+nβλ+81αλE1+n=0,(G′G)-3(1+n)/n:27E1+nnμ3+27E1+nμ3+9αE1+nμ3n+6E1+nn2μ3λ3+27βE1+nμ3+27αE1+nμ3=0,
with the solutions
(25)λ=0,μ=0,c=3En[(3-n)(3-2n)+3(3-n)α+9β]n3,
where β, n, and α are arbitrary constants. Hence, we obtain the exact traveling wave solution of (2) as follows:
(26)u(ζ)=E[c2c1+c2ζ]-3/n.
For the comparison between our solution (26) with that of King’s as given in (20), first we assume c1=0 and c2=1 and then we get the same as that of King’s (20) if we take c as in (25) in (20).
3.3. Exact Traveling Wave Solution of the Second Generalization of the Thin Film Equation (3)
This part is primarily concerned with the Cauchy problem for the doubly degenerate equation (3). King [9] gave the solution of (3) in the form:
(27)u=-c|F|k-1F[-(x-ct)3k/(n+k-1)],
where F=(k+n-1)3/(3k(2k+1-n)(k+2-2n)) and here k=m as in King [9].
The traveling wave variable (5) permits us to convert (3) into an ordinary differential equation as follows:
(28)-cu+un|u′′′|k-1u′′′=0.
Considering the homogeneous balance between un|u′′′|k-1u′′′ and u in (28), we obtain m=-3k/(n+k-1). Therefore, we can write the following:
(29)u(ζ)=E(G′G)-3k/(n+k-1),
for the traveling wave solutions of (28). By substituting (29) together with (9) into (28), clearing the denominator, and setting the coefficients of (G′/G)i,(i=0,1,…,7) to zero, we have the following algebraic system for E, λ, μ, n, and k:
(30)(G′G)-3n/(n+k-1):6E1+nk+15E1+nk2+6E1+nk3-15E1+nk2n+6E1+nkn2-12E1+nkn-cE=0,(G′G)-k(5n+5k-2)/(n+k-1):-18E1+nkλnμ2+9E1+nkλμ2n2+72E1+nk2nμ2λ2+9E1+nkλμ2-72E1+nλk2μ2+144E1+nk3λμ2=0,(G′G)-k(4n+4k-1)/(n+k-1):-6E1+nλ2nμ+33E1+nμλ2k2+6E1+nkμ2+114E1+nk3μ2+3E1+nμλ2n2+6E1+nkμ2n2+111E1+nμλ2k3+39E1+nnμ2k2-39E1+nk2μ2-33E1+nμλ2k2+3E1+nkμλ2-12E1+nknμ2=0,(G′G)-3k(n+k)/(n+k-1):6E1+nkλμn2-12E1+nkλμn+12E1+nλμnk2+27E1+nk3λ3+6E1+nkλμ+168E1+nλμk3-12E1+nλμk2=0,(G′G)-k(n+k+2)/(n+k-1):36E1+nλk3+9E1+nkλ+36E1+nλk2-18E1+nknλ+E1+nkλn2-36E1+nλnk2=0,(G′G)-k(2n+2k+1)/(n+k-1):21E1+nλ2k2+57E1+nλ2k3+3E1+nkμ+3E1+nkλ2+15E1+nμk2+60E1+nμk3-6E1+nknλ2+6E1+nkμn2-15E1+nnμk2-21E1+nnk2+3E1+nkn2λ2-12E1+nknμ=0,(G′G)-3k(2n+2k-1)/(n+k-1):60E1+nk3μ3-234E1+nk3μ6+6E1+nkμ3n2-12E1+nknμ3+39E1+nnμ3k2=0.
Solving this algebraic system by the use of Maple, we get the solutions for (3) as follows:
(31)λ=0,μ=0,c=3Enk(2k+1-n)(k+2-2n),
where k, n, and E are arbitrary constants. From (31) and (29), we obtain the exact traveling wave solution as follows:
(32)u(x,t)=u(ζ)=E[c2c1+c2ζ]-3k/(n+k-1).
Choosing c1=0 and c2=1 in (32), we get
(33)u(x,t)=u(ζ)=E(ζ)-3k/(n+k-1).
Now this is exactly the same as (27) if we substitute c as given in (31) into (27).
4. Conclusion
In this paper, we provide another instance of the applications of the (G′/G)-expansion method to the still very limited case whereby the balance numbers are not positive integers; see Zhang [22]. We have obtained some new exact traveling wave solutions of the thin film equation and its two generalizations. The solitary wave solutions are derived from these functions when the parameters are taken as special values. The Zhang technique [22] used in this paper is more effective and more general than that originally proposed by Wang et al. [10]. In all the general solutions (17), (26), and (32), we have the additional arbitrary constants c1 and c2. We note that the special case c1=0 and c2=1 reproduced the results of Bertozzi and Pugh [23] and King [9] with an appropriate choice of c. The new type of exact traveling wave solutions obtained in this paper for thin film equation and its two generalizations could be of beneficial use in future studies.
O'BrienS. B. G.SchwartzL. W.Theory and modeling of thin film flows2002635283BuckinghamR.ShearerM.BertozziA.Thin film traveling waves and the Navier slip condition200263272274410.1137/S0036139902401409MR1951957KumarM.SinghN.Phase plane analysis and traveling wave solution of third order nonlinear singular problems arising in thin film evolution20126492886289510.1016/j.camwa.2012.05.003MR2980481MyersT. G.Thin films with high surface tension199840344146210.1137/S003614459529284XMR1642807ZBL0908.35057MomoniatE.An investigation of an Emden-Fowler equation from thin film flow201228230030710.1007/s10409-012-0007-9MR2927872ConstantinP.DupontT. F.GoldsteinR. E.KadanoffL. P.ShelleyM. J.ZhouS.-M.Droplet breakup in a model of the Hele-Shaw cell19934764169418110.1103/PhysRevE.47.4169MR1377907KnüpferH.Well-posedness for the Navier slip thin-film equation in the case of partial wetting20116491263129610.1002/cpa.20376MR2839301ZBL1227.35146GreenspanH. P.On the motion of a small viscous droplet that wets a surface19788411251432-s2.0-0017926990KingJ. R.Two generalisations of the thin film equation2001347-873775610.1016/S0895-7177(01)00095-4MR1858796WangM.LiX.ZhangJ.The (G′/G)-expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics2008372441742310.1016/j.physleta.2007.07.051MR2381823KilicmanA.AbazariR.Travelling wave solutions of the Schrdinger-Boussinesq system201220121110.1155/2012/198398198398JafariH.KadkhodaN.BiswasA.The (G′/G)-expansion method for solutions of evolution equations from isothermal magnetostatic atmospheres20121576210.1016/j.jksus.2012.02.002AyhanB.BekirA.The (G′/G)-expansion method for the nonlinear lattice equations2012173490349810.1016/j.cnsns.2012.01.009MalikA.ChandF.KumarH.MishraS.C.Exact solutions of the Bogoyavlenskii equation using the multiple (G'/G)-expansion method20126492850285910.1016/j.camwa.2012.04.018MR2980477ZayedE. M. E.AbdelazizM. A. M.The two-variable (G′/G, 1/G)-expansion method for solving the nonlinear KdV-mKdV equation201220121476857310.1155/2012/725061NaherH.AbdullahF. A.Some new traveling wave solutions of the nonlinear reaction diffusion equation by using the improved (G'/G)-expansion method201220121710.1155/2012/871724871724MR2922132ElboreeM. K.Hyperbolic and trigonometric solutions for some nonlinear evolution equations201217114085409610.1016/j.cnsns.2012.03.029NaherH.AbdullahF. A.The (G′/G)-expansion method for abundant traveling wave solutions of Caudrey-Dodd-Gibbon equation201120111121821610.1155/2011/218216KabirM. M.BorhanifarA.AbazariR.Application of (G′/G)-expansion method to regularized long wave (RLW) equation20116182044204710.1016/j.camwa.2010.08.064MR2785540FengJ.LiW.WanQ.Using (G′/G)-expansion method to seek the traveling wave solution of Kolmogorov-Petrovskii-Piskunov equation2011217125860586510.1016/j.amc.2010.12.071MR2770205MalikA.ChandF.MishraS. C.Exact travelling wave solutions of some nonlinear equations by (G′/G)-expansion method201021692596261210.1016/j.amc.2010.03.103MR2653072ZhangH.New application of the (G′/G)-expansion method20091432203225BertozziA. L.PughM.The lubrication approximation for thin viscous films: regularity and long-time behavior of weak solutions19964928512310.1002/(SICI)1097-0312(199602)49:2<85::AID-CPA1>3.3.CO;2-VMR1371925ZBL0863.76017BertozziA. L.PughM.The lubrication approximation for thin viscous films: the moving contact line with a “porous media” cut-off of van der Waals interactions19947615351564MR130443810.1088/0951-7715/7/6/002ZBL0811.35045