Abstract
The homogeneous balance of undetermined coefficients method is firstly proposed to solve such nonlinear partial differential equations (PDEs), the balance numbers of which are not positive integers. The proposed method can also be used to derive more general bilinear equation of nonlinear PDEs. The Eckhaus equation, the KdV equation and the generalized Boussinesq equation are chosen to illustrate the validity of our method. The proposed method is also a standard and computable method, which can be generalized to deal with some types of nonlinear PDEs.
1 Introduction
Nonlinear PDEs are known to describe a wide variety of phenomena not only in physics, but also in biology, chemistry and several other fields. In recent years, many powerful methods to construct the exact solutions of nonlinear PDEs have been established and developed, which are one of the most excited advances of nonlinear science and theoretical physics. The exact solutions of nonlinear PDEs play an important role since they can provide much physical information and more insight into the physica1 aspects of the problems and thus lead to further applications. Seeking the exact solutions of nonlinear PDEs has long been an interesting topic in the nonlinear mathematical physics. With the development of soliton theory, various methods for obtaining the exact solutions of nonlinear PDEs have been presented, such as the inverse scattering method [1], the Bäcklund and Darboux transformation method [2], the homotopy perturbation method [3], the first integral method [4], the variational iteration method [5], the Riccati-Bernoulli sub-ODE method [6], the Jacobi elliptic function method [7], the tanhsech method [8], the (G'/G)-expansion method [9, 10], the Hirota’s method [11], the homogeneous balance method (HBM) [12, 13], the differential transform method (DTM) [14–17] and so on.
In these traditional methods, the HBM and the DTM are straightforward and popular tools for handling many types of functional equations. Recently, the Adomian’s decomposition method (ADM) for solving differential and integral equations, linear or nonlinear, has been the subject of extensive analytical and numerical studies because the ADM provides the solution in a rapid convergent series with elegantly computable components [18–23]. The Jacobi pseudo spectral approximation method [24], the fully spectral collocation approximation method [25, 26], the Jacobi tau approximation method [27], and the homotopy perturbation Sumudu transform method [28,29] are powerful and effective tools for solving nonlinear PDEs.
As a direct method, the HBM provides a convenient analytical technique to construct the exact solutions of nonlinear PDEs and has been generalized to obtain multiple soliton (or multiple solitary-wave) solutions [12, 13]. Fan improved this method to investigate the Bäcklund transformation, Lax pairs, symmetries and exact solutions for some nonlinear PDEs [30]. He also showed that there are many links among the HBM, Weiss-Tabor-Carnevale method and Clarkson-Kruskal method [31].
However, the HBM usually encounters complicated and tedious algebraic calculation. The exact solutions are fixed and single types when we deal with nonlinear PDEs by using the HBM. Moreover, the balance numbers of all nonlinear PDEs dealt with by this method usually are limited to positive integers.
Based on these problems, the homogeneous balance of undetermined coefficients method is proposed to improve the HBM, to derive more general bilinear equation of the KdV equation and the generalized Boussinesq equation, and to the exact solution of the Eckhaus equation, the balance numbers of which is not a positive integer.
The remainder of this paper is organized as follows: the homogeneous balance of undetermined coefficients method is described in Section 2. In Section 3, the exact solutions of the Eckhaus equation are obtained by using the homogeneous balance of undetermined coefficients method. In Section 4, the bilinear equation of the KdV equation and the generalized Boussinesq equation are derived respectively. A brief conclusion is given in Section 5.
2 Description of the homogeneous balance of undetermined coefficients method
Let us consider a general nonlinear PDE, say, in two variables,
where P is a polynomial function of its arguments, the subscripts denote the partial derivatives. The homogeneous balance of undetermined coefficients method consists of three steps.
Step 1 Suppose that the solution of Eq. (1) is of the form
where u = (u(x,t), w = w(x,t),
By balancing the highest nonlinear terms and the highest order partial derivative terms, m, n and b can be determined.
Step 2 Substituting Eq. (2) into Eq. (1) and arranging it at each order of w yield an equation as follows:
where fl (l = 0, 1, ...) are differential-algebraic expressions of w and akj. Setting fl = 0 and using compatible condition (wxt = wxt) yield a set of differential-algebraic equations (DAEs).
Step 3 Solving the set of DAEs, w and akj (k = 0, 1, ...,m, j = 0, 1, ...,n) can be determined. By substituting w, m, n, b and akj into Eq. (2), the exact solutions of Eq. (1) can be obtained.
In the following section, the exact solutions of Eckhaus equation can be obtained by using the homogeneous balance of undetermined coefficients method.
3 Application to the Eckhaus equation
Let us consider the Eckhaus equation in the form
where
The Eckhaus equation was found as an asymptotic multiscale reduction of certain classes of nonlinear Schrödinger type equations [32]. In [33], a lot of the properties of the Eckhaus equation were obtained. In [34], the Eckhaus equation is linearized by making some transformations of dependent and independent variables. Exact traveling wave solutions of the Eckhaus equation can be obtained by the (G'/G)-expansion method [10] and the first integral method [4].
In this section, by applying the homogeneous balance of undetermined coefficients method to the Eckhaus equation, new exact solutions of the Eckhaus equation can be obtained.
Suppose that the solution of Eq. (4) is of the form
where m, n, b, c, d and akj (k = 0, 1, ...,m, j = 0, 1, ...,n) are constants to be determined later.
Balancing |Ψ|4Ψ and Ψxx in Eq. (4), it is required that 3mb + 1 = 5mb, 3nb = 5nb. Solving this set of equations, we get m = 1, n = 0,
where a, b, c(a10 = a, a00 = b) and d are real constants to be determined later.
Substituting Eq. (6) into Eq. (4), we get
where
Obviously, Ψ is a solution of Eq. (4) provided that fk = 0 (k = 0, 1, 2, 3, 4) and gj = 0 (j = 1, 2, 3). Firstly, suppose that wx = 0, from Eqs. (6) and (7), we get an exact solution of Eq. (4) as follows:
where C0 and c are arbitrary real constants.
Secondly, suppose that wx ≠ 0, setting f0 = 0 and g1 = 0, we get
Substituting Eq. (9) into Eqs. (7g) and (7h), we get g2 = g3 = 0.
Case A When
Case A-1 When b = 0, substituting b = 0 into Eqs. (10), we get
Setting f1 = f2 = 0 and noting that wx ≠ 0, we get
From Eqs. (9) and (12), we get
where C0, C1 and c are arbitrary real constants.
Substituting
where C0, C1 and c are arbitrary real constants.
Case A-2 When b ≠ 0, setting f1 = f4 = 0 and noting that wx ≠ 0, we get
Substituting Eqs. (15) into Eqs. (10b) and (10c), we get
Setting f2 = f3 = 0, we get bwx = 0 from Eqs. (16), which contradicts with b ≠ 0 and wx ≠ 0. Therefore, Eq. (4) has no solution in this case.
Case B When
Case B-1 When b = 0, substituting b = 0 into Eqs. (17), we get
Setting f2 = 0 yields an equation as follows:
where α = 4(c2 + d).
Using transformations Z = wx and
Solving Eq. (20), we get
namely
where β(t) is an arbitrary function of t.
(1) When α = β(t) = 0, from Eq. (21) we get
where C1(t) and C2(t) are arbitrary functions of t.
Substituting the above equation into Eq. (9), we get
Setting the coefficients of xj (j = 0, 1) to zero in the above equation, we get
Solving the above equations, we get
where Ck (k = 0, 1) are arbitrary real constants.
Then we get
Substituting
where Ck (k = 0, 1) and c are arbitrary real constants.
(2) When β(t) = 0, α > 0, similar to (1), we get
and an exact solution of Eq. (4) as follows:
where ε = ±1, Ck (k = 0, 1), c and d are arbitrary real constants.
Especially, if C0C1 > 0, then Eq. (24) can be reduced to
where Ck (k = 0, 1), c, d and ξ0 are arbitrary real constants.
If C0C1 < 0, then Eq. (24) can be reduced to
where Ck (k = 0, 1), c, d and ξ0 are arbitrary real constants.
(3) When β(t) ≠ 0, α = 0, similar to (1), we get
and an exact solution of Eq. (4) as follows:
where Ck (k = 0, 1) and c are arbitrary constants.
(4) When β(t) ≠ 0, α > 0, similar to (1), we get
and an exact solution of Eq. (4) as follows:
where
(5) When β(t) ≠ 0, α < 0, similar to (1), we get
and an exact solution of Eq. (4) as follows:
where
Case B-2 When b ≠ 0, setting f2 = f3 = f4 = 0 in Eqs. (17b), (17c) and (17d), we get
Solving Eqs. (32) and noting that Eq. (9), we get
Substituting
In this section, by applying the homogeneous balance of undetermined coefficients method to the Eckhaus equation, new exact solutions of the Eckhaus equation are obtained. Among the solutions of the Eckhaus equation, Ψj (j = 3, 4, 5, 6, 7, 8) are the same as the results of [4, 10]. To our knowledge, other solutions of Eckhaus equation, Ψj (j = 2, 9, 10, 11) have not been reported in any literature.
4 Application to derive the bilinear equation of nonlinear PDEs
In this section, firstly, we modify the homogeneous balance of undetermined coefficients method to derive the bilinear equation of nonlinear PDEs. Then, more general bilinear equation of the KdV equation and the generalized Boussinesq equation are obtained by using the proposed method.
Let us still consider Eq. (1). Suppose that the solution of Eq. (1) is Eq. (2). By balancing the highest nonlinear terms and the highest order partial derivative terms, m, n and b can be determined. Substituting Eq. (2) into Eq. (1) and balancing the terms with
Next, the KdV equation and the generalized Boussinesq equation are chosen as examples to illustrate our method.
Let us consider the celebrated KdV equation [9] in the form
Suppose that the solution of Eq. (34) is Eq. (2). Balancing uxxx and uux in Eq. (34), it is required that mb + 3 = 2mb + 1, nb = 2nb. Choosing m = 2, n = 0 and b = 1, Eq. (2) can be written as
where aj0(j = 0, 1, 2) are constants to be determined later.
From Eq. (35), one can calculate the following derivatives:
Equating the coefficients of
Solving the above algebraic equations, we get a20 = 12δ, a10 = 0. Substituting a20 and a10 back into Eq. (35), we get
where a00 is an arbitrary constant.
Substituting Eq. (37) into Eq. (34) and simplifying it, we get
where
Simplifying Eq. (38) and integrating once with respect to x, we get
Eq. (39) is identical to
where C(t) is an arbitrary function of t, and a00 is an arbitrary constant.
Especially, taking C(t) as zero in Eq. (40), we get the bilinear equation of Eq. (34) as follows:
Eq. (41) can be written concisely in terms of D-operator as
where
Applying the Hirota’s method [11], the bilinear equation of Eq. (34) can be written as
Eq. (43) is obtained by settinga00 = 0 in Eq. (42). Obviously, Eq. (43) is a special case of Eq. (42). Therefore, more general bilinear equation of the KdV equation is obtained by using our method.
Applying the perturbation method [11] to Eq. (42), we can get 1-soliton solution and 2-soliton solution of Eq. (42) as follows:
Substituting Eqs. (44) into Eq. (37), 1-soliton solution and 2-soliton solution of Eq. (34) can be obtained. Similarly, N-soliton solution of Eq. (34) can be obtained.
The generalized Boussinesq equation reads
Suppose that the solution of Eq. (45) is Eq. (2). In order to balance uxxxx and uux in Eq. (45), it is required that mb + 4 = 2mb + 2,nb = 2nb. Choosing m = 2, n = 0 and b = 1, Eq. (2) can be written as
where aj0 (j = 0, 1, 2) are constants to be determined later.
Substituting Eq. (46) into Eq. (45) and equating the coefficients of
where a00 is an arbitrary constant.
Substituting Eq. (47) into Eq. (45), we get
where
Simplifying Eq. (48) and integrating twice with respect to x, we get
Eq. (49) is identical to
where C1 (t) and C2 (t) are arbitrary functions of t, and a00 is an arbitrary constant.
Especially, letting C1 (t) = C2 (t) = 0 in Eq. (50), we get the bilinear equation of Eq. (45) as follows:
Eq. (51) can be written concisely in terms of D-operator as
where a00 is an arbitrary constant.
So far, the proposed method is successfully used to establish the bilinear equation of the KdV equation and the generalized Boussinesq equation. Our method can be used to derive the bilinear equations of some nonlinear PDEs.
5 Conclusions
The homogeneous balance of undetermined coefficients method is successfully used to solve such nonlinear PDEs, the balance numbers of which are not positive integers. By applying our method to the Eckhaus equation, more exact solutions are obtained. The proposed method can also be used to derive the bilinear equation of nonlinear PDEs. More general bilinear equation of the KdV equation and the generalized Boussinesq equation can be obtained by applying our method.
Once the bilinear equation of nonlinear PDE is obtained, the perturbation method can be employed to obtain the N -soliton solution for the nonlinear PDE. Many exact solutions of the nonlinear PDE are obtained by using the three-wave method and the homoclinic test approach. Generally, some nonlinear PDE can be linearized or homogenized by using the homogeneous balance of undetermined coefficients method.
Many well-known nonlinear PDEs can be handled by our method. The performance of our method is found to be simple and efficient. The availability of computer systems like Maple facilitates the tedious algebraic calculations. Our method is also a standard and computable method, which allows us to solve complicated and tedious algebraic calculations.
Acknowledgement
The research is supported by the National Natural Science Foundation of China (61402375), the Fundamental Research Funds for the Central Universities (G2015KY0103), the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University (CX201517).
References
[1] Ablowitz, M.J. and Clarkson, P.A., Solitons, nonlinear evolution equations and inverse scattering, Cambridge University, 199110.1017/CBO9780511623998Search in Google Scholar
[2] Rogers, C. and Schief, W.K., Bäcklund and Darboux transformation geometry and modern applications in solitons theory, Cambridge University, 200210.1017/CBO9780511606359Search in Google Scholar
[3] He, J.H., An approximate solution technique depending on an artificial parameter: A special example, Commun. Nonlinear Sci. Numer. Simul., 1998, 3, 92-9710.1016/S1007-5704(98)90070-3Search in Google Scholar
[4] Taghizadeh, N.; Mirzazadeh, M.; Filiz, T., The first-integral method applied to the Eckhaus equation, Appl. Math. Lett., 2012, 25, 798-80210.1016/j.aml.2011.10.021Search in Google Scholar
[5] Yusufoglu, E., The variational iteration method for studying the Klein-Gordon equation, Appl. Math. Lett., 2008, 21, 669-67410.1016/j.aml.2007.07.023Search in Google Scholar
[6] Yang,X.F.; Deng,Z.C.; Wei,Y., A Riccati-Bernoulli sub-ODE method for nonlinear partial differential equations and its application, Adv. Diff. Equa., 2015, 1, 117-13310.1186/s13662-015-0452-4Search in Google Scholar
[7] Yan, Z.L., Abunbant families of Jacobi elliptic function solutions of the (2+1)-dimensional integrable Davey-Stewartson-type equation via a new method, Chaos, Solitons and Fract., 2003, 18, 299-30910.1016/S0960-0779(02)00653-7Search in Google Scholar
[8] Wazwaz, A.M., The tanh method: exact solutions of the sine-Gordon and the sinh-Gordon equations, Appl. Math. Comput., 2005, 167, 1196-121010.1016/j.amc.2004.08.005Search in Google Scholar
[9] Wang, M.L.; Li, X.Z.; Zhang, J.L., The (G0=G)-expansion method and traveling wave solutions of nonlinear evolution equations in mathematical physics, Phys. Lett. A, 2008, 372, 417-42310.1016/j.physleta.2007.07.051Search in Google Scholar
[10] Zhang, H.Q., New application of the (G0=G)-expansion method, Commun. Nonlinear Sci. Numer. Simul., 2009, 14, 3220-322510.1016/j.cnsns.2009.01.006Search in Google Scholar
[11] Hirota, R., Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons, Phys. Rev. Lett., 1971, 27, 1192-119410.1103/PhysRevLett.27.1192Search in Google Scholar
[12] Wang, M.L., Solitary wave solutions for variant Boussinesq equations, Phys. Lett. A, 199 (1995) 169-17210.1016/0375-9601(95)00092-HSearch in Google Scholar
[13] Wang, M.L. and Zhou, Y.B., Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics, Phys. Lett. A, 1996, 216, 67-7510.1016/0375-9601(96)00283-6Search in Google Scholar
[14] Chen, C.K. and Ho, S.H., Solving partial differential equations by two-dimensional differential transform method, Appl. Math. Comput., 1999, 106, 171-17910.1016/S0096-3003(98)10115-7Search in Google Scholar
[15] Fatoorehchia, H. and Abolghasemi, H., An explicit analytic solution to the Thomas-Fermi equation by the improved differential transform method, Acta. Phys. Polonica A, 2014, 125, 1083-108710.12693/APhysPolA.125.1083Search in Google Scholar
[16] Fatoorehchi, H. and Abolghasemi, H., Improving the differential transform method: A novel technique to obtain the differential transforms of nonlinearities by the Adomian polynomials, Appl. Math. Model., 2013, 37, 6008-601710.1016/j.apm.2012.12.007Search in Google Scholar
[17] Kurnaz, A.; Oturanc, G.; Kiris, M.E., n-Dimensional differential transformation method for solving PDEs, Int. J. Comput. Math., 2005, 82, 369-38010.1080/0020716042000301725Search in Google Scholar
[18] Adomian, G., Solutions of nonlinear PDE, Appl. Math. Lett., 1998, 11, 121-12310.1016/S0893-9659(98)00043-3Search in Google Scholar
[19] El-Sayed, S.M., The decomposition method for studying the Klein-Gordon equation, Chaos, Soliton and Fract., 2003, 18, 1025-103010.1016/S0960-0779(02)00647-1Search in Google Scholar
[20] Fatoorehchia, H. and Abolghasemi, H., Adomian decomposition method to study mass transfer from a horizontal flat plate subject to laminar fluid flow, Adv. Nat. Appl.Sci., 2011, 5, 26-33Search in Google Scholar
[21] Fatoorehchi, H. and Abolghasemi, H., Approximating the minimum reflux ratio of multicomponent distillation columns based on the Adomian decomposition method, J. Taiwan Inst. Chem. Eng., 2014, 45, 880-88610.1016/j.jtice.2013.09.032Search in Google Scholar
[22] Fatoorehchi, H.; Abolghasemi, H.; Rach, R., A new parametric algorithm for isothermal flash calculations by the Adomian decomposition of Michaelis-Menten type nonlinearities, Fluid Phase Equilibria, 2015, 395, 44-5010.1016/j.fluid.2015.03.024Search in Google Scholar
[23] Fatoorehchi, H.; Abolghasemi, H.; Zarghami, R., Analytical approximate solutions for a general nonlinear resistor-nonlinear capacitor circuit model, Appl. Math. Model., 2015, 1, 38-4810.1016/j.apm.2015.01.038Search in Google Scholar
[24] Bhrawy, A.H., An efficient Jacobi pseudospectral approximation for nonlinear complex generalized Zakharov system, Appl. Math. Comput., 2014, 247, 30-4610.1016/j.amc.2014.08.062Search in Google Scholar
[25] Bhrawy, A.H. and Abdelkawy, M.A., A fully spectral collocation approximation formulti-dimensional fractional Schrödinger equations, J. Comput. Phys., 2015, 294, 462-48310.1016/j.jcp.2015.03.063Search in Google Scholar
[26] Bhrawy, A.H. and Zaky, M.A., Numerical simulation for two-dimensional variable-order fractional nonlinear cable equation, Nonlinear Dyn., 2015, 80, 101-11610.1007/s11071-014-1854-7Search in Google Scholar
[27] Bhrawy, A.H. and Zaky, M.A., A method based on the Jacobi tau approximation for solving multi-term time-space fractional partial differential equations, J. Comput. Phys., 2015, 281, 876-89510.1016/j.jcp.2014.10.060Search in Google Scholar
[28] Dubey, R.S.; Alkahtani, B.S.T.; Atangana, A., Analytical solution of space-time fractional Fokker-Planck equation by homotopy perturbation Sumudu transform method, 2014, Article ID 78092910.1155/2015/780929Search in Google Scholar
[29] Kumar, D.; Purohit, S.D.; Secer, A.; Atangana, A., On generalized fractional kinetic equations involving generalized Bessel function of the first kind, 2015, Article ID 28938710.1155/2015/289387Search in Google Scholar
[30] Fan, E.G., Bäcklund transformation and similarity reductions of nonlinear partial differential equations using extended homogeneous balance method, Commun. Theor. Phys., 2001, 35, 523-52610.1088/0253-6102/35/5/523Search in Google Scholar
[31] Fan, E.G., Connections among homogeneous balance method, Weiss-Tabor-Carnevale method and Clarkson-Kruskal method, Acta. Math. Phys. Sin., 2000, 49, 1409-141210.7498/aps.49.1409Search in Google Scholar
[32] Calogero, F. and Eckhaus, W., Nonlinear evolution equations, rescalings, model PDEs and their integrability: I, Inv. Probl., 1987, 3, 229-26210.7498/aps.49.1409Search in Google Scholar
[33] Calogero, F. and Lillo, S.D., The Eckhaus PDE i ψt + ψxx + 2(|ψ|2)xψ+|ψ|4ψ = 0, Inv. Probl., 1987, 4, 633-68210.1088/0266-5611/3/4/012Search in Google Scholar
[34] Calogero, F., The evolution partial differential equation ut = uxxx + 3(uxxu2 + 3
© 2016 Wei et al., published by De Gruyter Open
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.