1Multisoliton and other invariant solutions of combined KdVnKdV equation by using Lie Symmetry Approach
Abstract.
Lie symmetry method is applied to investigate symmetries of the combined KdVnKdV equation, that is a new integrable equation by combining the KdV equation and negative order KdV equation. Symmetries which are obtained in this article, are further helpful for reducing the combined KdVnKdV equation into ordinary differential equation. Moreover, a set of eight invariant solutions for combined KdVnKdV equation is obtained by using proposed method. Out of the eight solutions so obtained in which two solutions generate progressive wave solutions, five are singular solutions and one multisoliton solutions which is in terms of WeierstrassZeta function.
Key words and phrases:
Combined KdVnKdV equation; Lie symmetry method; Invariant solutions; Similarity reduction.1. Introduction
Korteweg and Vries derived KdV equation [6] to model Russell’s phenomenon of solitons like shallow water waves with small but finite amplitudes [26]. Solitons are localized waves that propagate without change of it’s shape and velocity properties and stable against mutual collision [7, 22]. It has also been used to describe a number of important physical phenomena such as magneto hydrodynamics waves in warm plasma, acoustic waves in an inharmonic crystal and ionacoustic waves [17, 2, 8]. A special class of analytical solutions of KdV equation, the socalled traveling waves, for nonlinear evolution equations (NEEs) is of fundamental importance because a lot of mathematicalphysical models are often described by such a wave phenomena. Thus, the investigation of traveling wave solutions is becoming more and more attractive in nonlinear science nowadays. However, not all equations posed in these fields are solvable. As a result, many new techniques have been successfully developed by a diverse group of mathematicians and physicists, such as Rational function method [31, 32], Bäcklund transformation method [14], Hirotaâs bilinear method [5, 13], Lie symmetry method[1, 30, 18, 9, 10, 11], Jacobi elliptic function method [12], Sinecosine function method [27], Tanhcoth function method [15], Weierstrass function method [20], Homogeneous balance method [28], Expfunction method [4], expansion method [21], etc. But, it is extremely difficult and time consuming to solve nonlinear problems with the wellknown traditional methods. This work investigates the combined KdVnKdV equation
(1) 
where . We apply Lie symmetry analysis on combined KdVnKdV equation, first constructed by Wazwaz [25] using recursion operator [19]. In addition, the combined KdVnKdV equation (1) possesses the Painlevé property for complete integrability [3]. In this paper, Lie point symmetry generators of the combined KdVnKdV equation were derived. Similarity reductions and number of explicit invariant solutions for the equation using Lie symmetry method were obtained. All the new invariant solutions of combined KdVnKdV were analyzed graphically. Also, 1multisoliton solution obtained in terms of WeierstrassZeta function which appear in classical mechanics such as, motion in cubic and quartic potentials, description of the movement of a spherical pendulum, and in construction of minimal surfaces [29]. Some of the outcomes are interesting in physical sciences and are beautiful in mathematical sciences.
The organization of the paper is as follows. In Sec. 2, we discuss the methodology of Lie symmetry analysis of the general case. In Sec. 3, we obtain infinitesimal generators and the Lie point symmetries of the Eq. (1). In Sec. 4, symmetry reductions and exact group invariant solutions for the combined KdVnKdV eqation were obtained. In Sec. 5, we discussed all the invariant solutions graphically by Figures 1, 2, 3 and 4. Finally, concluding remarks are summarized in Section 6.
2. Method of Lie symmetries
Let us consider a system of partial differential equations as follows:
(2) 
where denotes all the derivatives of of all orders from to . The oneparameter Lie group of infinitesimal transformations for Eq. (2) is given by
where is the group parameter, and the Lie algebra of Eq. (1) is spanned by vector field of the form
(3) 
A symmetry of a partial differential equation is a transformation which keeps the solution invariant in the transformed space. The system of nonlinear PDEs leads to the following invariance condition under the infinitesimal transformations
In the above condition, is termed as order prolongation [16] of the infinitesimal generator which is given by
(4) 
the second summation being over all (unordered) multiindices , . The coefficient functions of are given by the following expression
(5) 
where , and denotes total derivative.
3. Lie symmetry analysis for the combined KdVnKdV equation
In this section, authors explained briefly all the steps of the STM method to keep this work selfconfined. The Lie symmetries for the Eq. (1) have generated and then its similarity solutions are found. Therefore, one can consider the following oneparameter () Lie group of infinitesimal transformations
(6) 
where and are infinitesimals for the variables and respectively, and is the solution of Eq. (1). Therefore, the associated vector field is
(7) 
Lie symmetry of Eq. (1) will be generated by Eq. (7). Use fourth prolongation gives rise to the symmetry condition for Eq. (1) as follows:
(8) 
where , and are the coefficient of , values are given in many references [16, 1]. Incorporating all the expressions into Eq. (3), and then equating the various differential coefficients of to zero, we derive following system of Eight determining equation
(9) 
Solving the above system of equations, we obtain following infinitesimals for (1) using software Maple,
(10) 
where and are arbitrary constants whereas is an arbitrary function.
The symmetries under which Eq. (1) is invariant can be spanned by the following four infinitesimal generators if we assume , a constant Then all of the infinitesimal generators of Eq. (1) can be expressed as
where
(11) 
*  

0  
0  0  0  
0  0  0  
0  0  0 
The vector field yield commutation relations through the Table 1. The th entry in Table 1 is the Lie bracket . Table 1 is skew  symmetric with zero diagonal elements. Table 1 shows that the generators and are linearly independent. Thus, to obtain the similarity solutions of Eq. (1), the corresponding associated Lagrange system is
(12) 
4. Invariant solutions of the combined KdVnKdV equation
To proceed further, selection of and by assigning the particular values to ’s , provide new physically meaningful solutions of Eq. (1). In order to obtain symmetry reductions and invariant solutions, one has to solve the associated Lagrange equations given by (12).
Now, let us discuss the following particular cases for various forms of :
Case 1: For , then Eqs. (10) and (12) gives
(13) 
The similarity form suggested by Eq. (13) is given by
(14) 
with similarity variable
(15) 
where
(16) 
Inserting the value of from Eq. (14) into Eq. (1), we get the following fourth order ordinary differential equation
(17) 
where is given by Eq. (15) and , etc.
Any how, we could not find the general solution of Eq. (17) still two particular solutions are found as below
(18) 
where and are arbitrary constants. Thus, from Eqs. (14) and (18), we get two invariant solutions of Eq. (1) given below
(19)  
(20) 
where , and are given by Eq. (16).
Case 2: For , Eq. (12) are of the form
(21) 
The similarity solution to Eq. (21) can be written as
(22) 
where is a similarity variable.
Inserting the value of from Eq. (22) in Eq. (1), we get the fourth order ordinary differential equation in
(23) 
Assume . Without loss of generality, we can assume , then Eq. (34) reduces to
(24) 
The general solution of Eq. (24) in terms of WeiestrassZeta function as
(25) 
where and and are arbitrary constants. From Eq. (25) with Eq. (22), we have another invariant solution of Eq. (1)
(26) 
Case 2A: , Eq. (13) becomes
(27) 
In this case, we get
(28) 
where . Substituting the value of in Eq. (1), again we obtain the same fourth order ordinary differential equation (17). Some particular solutions are given below
(29) 
where and are arbitrary constants. Therefore, using Eq. (29) in Eq. (28), we obtain the following two exact solutions for Eq. (1)
(30)  
(31) 
Case 3: For , Eq. (13) modified as
(32) 
The group invariant solution is given as
(33) 
with . Substituting the value of from Eq. (33) into Eq. (1), we get the following reduced ordinary differential equation
(34) 
Hence, from Eq.(34) we found two values of given as
(35)  
(36) 
where is an arbitrary constant of integration. Here, Eq. (35) gives same solution as Eq. (20) with . Also, using Eq. (36) in Eq. (33) gives new invariant solution of Eq. (1) as
(37) 
where and are given by Eq. (16).
Case 3A: For , from Eq. (32), we have
(38) 
Therefore, the similarity transformation method gives
(39) 
where and is the similarity function. Substituting the value of in Eq. (1) we get reduced ordinary differential equation given as
(40) 
Again it gives two values of as
(41)  
(42) 
where is constant of integration. Here, Eq. (41) gives same solution as Eq. (20). Also, using (42) with (39) gives another new invariant solution of (1) as
(43) 
Case 4: For , from Eq. (32) we get similarity form suggested by Eq.(13) is given by
(44) 
with similarity variable and is given by (16). In this case, we get same ordinary differential equation as Eq. (17), one particular solution is as follows
(45) 
where is arbitrary constant. Thus, from Eqs. (44) and (45) we get another new invariant solution of Eq. (1) given below
(46) 
5. Discussion
The results of the combined KdVnKdV equation presented in this paper have richer physical structure then earliar outcomes in the literature [25]. The recorded results are significant in the context of nonlinear dynamics, physical science, mathematical physics etc. The invariant solutions obtained can illustrate various dynamic behaviour due to existence of arbitrary constants. The nonlinear behaviour of the results are analyzed in the following manner:
Figure 1:
Fig 1(a) shows progressive wave soluion in Eq. (19) for .
Fig 1(b) shows presence of singularity in the plane in Eq. (20) with values .
Figure 2: The evolution profiles of 1Multisoliton solution is given by Eq. (26) as shown in this Figure. We have recorded the physical nature with variation in parameters.
Fig 2(a) For and a plot of WeiestrassZeta function is shown;
Fig 2(b) For and only few solitons are shown;
Fig 2(c) Front orthographic projection is shown for and ;
Fig 2(d) Front orthographic projection is shown for and .
Figure 3:
Fig 3(a) For and , Eq. (30) exhibits progressive wave;
Fig 3(b), shows singularity in planes for Eqs. (31) for and ;
Fig 3(c, d) shows singularity wave profile for Eq. (37) and Eq. (43) which explain the transition of nonlinear behaviour in the form of opposite rotatory folded sheets.
Figure 4:
exhibits singularity near but asymptotic structures is observed near for parameters .
6. Conclusion
In this paper, the similarity reductions and invariant solutions for the combined KdVnKdV are presented. This paper obtained the 1multisoliton and other invariant solution of the equation. The method that was used to obtain the exact group invariant solutions is the Lie symmetry analysis approach. All the solutions are different from earlier work which have been obtained by Wazwaz [25]. Eventually, the stucture of combined KdVnKdV equation is an nontrivial one, which can be clearly seen from the graphically results of invariant solutions. The Lie symmetry analysis method extracts the new forms of analytic solutions which are of physical importance such as condensed matter physics and plasma physics.
Acknowledgment
The second author sincerely and genuinely thanks Department of Mathematics, SGTB Khalsa College, University of Delhi for financial support.
References
 [1] G. W. Bluman and S. Kumei, Symmetries and Differential Equations, (SpringerVerlag, New York, 1989).
 [2] K. Brauer, The Kortewegde Vries equation: history, exact solutions, and graphical representation, University of Osnabruck, (2006).
 [3] A. S. Fokas, Symmetries and integrability, Stud. Appl. Math. 77(3) (1987) 253–299.
 [4] J. H. He and X. H. Wu, Expfunction method for nonlinear wave equations, Chaos Solitons Fractals 30 (2006) 700–708.
 [5] R. Hirota, The Direct Method in Soliton Theory, (Cambridge University Pres, Cambridge, 2004).
 [6] D. J. Korteweg and G. de Vries, On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary wave, Philos. Mag. 39 (1895) 422–443.
 [7] A. J. Khattak and SirajulIslam, A comparative study of numerical solutions of a class of KdV equation, App. Math. Comput. 199 (2008) 425–434.
 [8] S. A. Khuri, Soliton and periodic solutions for higher order wave equations of KdV type (I), Chaos Solitons Fractals 26(1) (2005) 25–32.
 [9] S. Kumar and Y. K. Gupta, Generalized Invariant Solutions for Spherical Symmetric Nonconformally Flat Fluid Distributions of Embedding Class One, Int. J. Theor. Phys. 53 (2014) 2041–2050.
 [10] S. Kumar, Pratibha and Y. K. Gupta, Invariant Solutions of Einstein Field equation for nonconformally flat fluid spheres of embedding class one, Int. J. Mod. Phys. A 25 (2010) 3993–4000.
 [11] M. Kumar and A. K. Tiwari, Soliton solutions of BLMP equation by Lie symmetry approach, Comput. Math. Appl. 75 (2018) 1434–1442.
 [12] G. T. Liu and T.Y. Fan, New applications of developed Jacobi elliptic function expansion methods. Phys. Lett. A 345 (2005) 161–166.
 [13] J.J. Mao et al, Optical solitons, complexitons, Gaussian soliton and power series solutions of a generalized Hirota equation, Mod. Phys. Lett. B, 32, 1850127 (2018) https://doi.org/10.1142/S0217984918501439
 [14] R. M. Miura, Backlund Transformation, (Springer Verlag, New York, 1973).
 [15] W. Malfliet and W. Hereman, The tanh method I: exact solutions of nonlinear evolution and wave equations. Phys. Scripta 54 (1996) 563–568.
 [16] P. J. Olver, Applications of Lie Groups to Differential Equations Graduate Texts in Mathematics, Vol. 107 (SpringerVerlag, New York, 1986).
 [17] T. Ozis and S. Ozer, A simple similarity transformation iterative scheme applied to Korteweg de Vries equation, Appl. Math. Comput. 173 (2006) 19–32.
 [18] S. S. Ray, Lie symmetry analysis and reduction for exact solution of (2+1) dimensional BogoyavlenskyâKonopelchenko equation by geometric approach, Mod. Phys. Lett. B, 32, (2018) 1850127
 [19] J. M. Verosky, Negative powers of Olver recursion operators, J. Math. Phys. 32(7) (1991) 1733–1736.
 [20] M. L. Wang, Exact solutions for a compound KdVBurgers equation, Phys. Lett. A, 213 (1996) 279–287.
 [21] M. Wang, X. Li and J. Zhang, The (G’/G)expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics, Phys. Lett. A 372 (2008) 417–423.
 [22] A. M. Wazwaz, Partial Differential Equations and Solitary Waves Theory, (Springer, New York 2009).
 [23] A. M. Wazwaz, A study on the (2 + 1)dimensional KdV4 equation derived by using the KdV recursion operator, Math. Methods Appl. Sci. 36 (2013) 1760–1767.
 [24] A. M. Wazwaz, A new integrable equation constructed via combining the recursion operator of the CalogeroBogoyavlenskiiSchiff (CBS) equation and its inverse operator, Appl. Math. Inf. Sci. 11 (2017) 1241–1246.
 [25] A. M. Wazwaz, A new integrable equation that combines the KdV equation with the negativeorder KdV equation, Math. Methods Appl. Sci. 41 (2018) 80–87.
 [26] L. Wazzan, A modified tanhcoth method for solving the KdV and the KdVBurgers equations, Commun. Nonlinear Sci. Numer. Simul. 14(2) (2009) 443–450.
 [27] C. Yan, A simple transformation for nonlinear waves, Phys. Lett. A 224 (1996) 77–84.
 [28] Z.Y. Yan, An improved algebra method and its applications in nonlinear wave equations, Chaos Solitons Fractals 21 (2004) 1013–1021.

[29]
Applications of Weierstrass functions and inverses, (19982018) Wolfram Research, Inc., http://functions.wolfram.com/EllipticFunctions/WeierstrassZeta/introductions/Weierstrass
/ShowAll.html.  [30] Z. Zhao and B. Han, Lie symmetry analysis, Bäcklund transformations, and exact solutions of a (2 + 1)dimensional BoitiLeonPempinelli system, J. Math. Phys. 58, (2017) 101514 .
 [31] Y. Zhang and W.X. Ma, Rational solutions to a KdVlike equation, Appl. Math. Comp. 256 (2015) 252–256.
 [32] Y. Zhang and W.X. Ma, A Study on Rational Solutions to a KPlike Equation, Z. Naturforsch. 70(4)a (2015) 263–268.