Lie Group Analysis of a Forced KdV Equation

Many studies in mathematical physics, engineering, life sciences, and in all other sciences use mathematical models to describe certain phenomenon. Such models are represented by the nonlinear evolution equations. The existence of analytical solutions to these equations becomes a basis for a better understanding of the dynamics involved. The KdV equation [1] being a nonlinear evolution equationmodels the nonlinear wave phenomena on the shallow water surfaces, and its various forms have been proposed depending upon the applications in diverse fields of science and engineering. To date, a lot of solution procedures, both analytical and numerical, have been developed to solve these types of differential equations. However, some of these approaches may involve approximation of solutions. The current work is based upon the Lie group theory approach as a solution procedure. The Lie point or higher-order symmetries of a differential equation enable one to obtain the solutions invariant under a particular symmetry or a linear combination of symmetries. The invariant solutions are a basis to finding exact solutions or numerical solutions. In most real-life applications, the differential equations, which are used to model a certain situation, contain arbitrary functions of dependent variable or its derivatives and independent variables. Instead of assuming the forms of these model parameters, the method of group classification can be employed to systematically specify their forms. There are various approaches to solving a group classification problem, namely, the direct analysis or the approach based upon the equivalence transformations. We consider the forced KdV equation


Introduction
Many studies in mathematical physics, engineering, life sciences, and in all other sciences use mathematical models to describe certain phenomenon.Such models are represented by the nonlinear evolution equations.The existence of analytical solutions to these equations becomes a basis for a better understanding of the dynamics involved.The KdV equation [1] being a nonlinear evolution equation models the nonlinear wave phenomena on the shallow water surfaces, and its various forms have been proposed depending upon the applications in diverse fields of science and engineering.To date, a lot of solution procedures, both analytical and numerical, have been developed to solve these types of differential equations.However, some of these approaches may involve approximation of solutions.The current work is based upon the Lie group theory approach as a solution procedure.
The Lie point or higher-order symmetries of a differential equation enable one to obtain the solutions invariant under a particular symmetry or a linear combination of symmetries.The invariant solutions are a basis to finding exact solutions or numerical solutions.In most real-life applications, the differential equations, which are used to model a certain situation, contain arbitrary functions of dependent variable or its derivatives and independent variables.Instead of assuming the forms of these model parameters, the method of group classification can be employed to systematically specify their forms.There are various approaches to solving a group classification problem, namely, the direct analysis or the approach based upon the equivalence transformations.
We consider the forced KdV equation where  represents the elevation of free water surface and  and  are arbitrary constants which depend upon the long wave speed, .The arbitrary smooth function () is the forcing term.The variables  and  represent time and space, respectively.This equation was proposed in [2] and the Hirota's bilinear approach was used to obtain the multiple soliton solutions.However, it is worth mentioning that () remained unspecified.In [3] the generalized wave transformation was used to transform a forced KdV equation with time-dependent coefficients into a homogeneous equation, and the soliton solutions were obtained by making use of the solitary wave ansatz.Likewise, the time-dependent coefficients remained arbitrary.The symmetry-based approach is given in [4][5][6]; the investigations in these references include approximate symmetry classification, derivation of some conservation laws, and the construction of the solutions.As mentioned earlier, the arbitrary functions assume their forms in a systematic fashion via the method of group classification.This approach has been utilized on similar equations [7,8] to the underlying equation.
Here  0 , ,  and  are arbitrary constants.
This work is organized as follows.The next section deals with the generation of the determining equations for Lie point symmetries and includes the classifying relation for the forcing term.The functional forms of the forcing term are specified via the direct analysis of the classifying relation.In Section 3, some results of the Lie group analysis are utilized for symmetry reductions and exact solutions.Finally, we summarize our investigations in Section 4.

Lie Group Analysis
In Lie's algorithm (see [9][10][11] for more details), a vector field is a generator of Lie point symmetries of (1) if and only if  [3] ( where is the third prolongation of the vector field .The variables 's are given by the formulae where are the total derivative operators [10].The invariance conditions (3) are separated with respect to the powers of the derivatives of , and this yields the determining equations, which are a system of linear partial differential equations of homogeneous type in  1 ,  2 , and .It is easy but tedious to generate the determining equations manually.Nowadays there are many computer software packages for symbolic computation, which have been developed to find symmetries interactively or automatically.However, we are yet to develop the software package which solves the group classification problem whether complete or partial group classification.
The coefficients of the generator of Lie point symmetries (2), namely,  1 ,  2 , and , satisfy the determining equations where the subscripts denote partial differentiation with respect to the indicated variables, and a "prime" represents total derivative with respect to .
The manipulation of ( 7) leads to the coefficients of Lie point symmetry generator given by where  1 ,  2 are arbitrary constants, and () is an arbitrary function which satisfies the classifying relation Assume that () is an arbitrary smooth function of ; then from the classifying (9) we obtain where  3 and  4 are arbitrary constants.Thus, we have a twodimensional principal Lie algebra, which is spanned by the operators Our goal is to obtain the functional forms of the forcing term, (), for which the principal Lie algebra is extended.Therefore, the analysis of the classifying relation ( 9) considering the cases   () = 0 and   () ̸ = 0 yields the various forms of (), and their corresponding extensions of the principal Lie algebra are given in Table 1.It is worth mentioning that some of the obtained symmetry classification results are comparable with those found in [4,5].It is, however, noted that in [6] a more general case is considered.

Symmetry Reductions and Exact Solutions
It can be seen from Table 1 that the symmetry Lie algebra is three-dimensional for the first three cases and fourdimensional in the last case.We consider Case 3 to illustrate the procedure involved in performing similarity reductions.Since the symmetry Lie algebra is three-dimensional, we look for solutions invariant under the linear combination of the operators which are the symmetries of the equation In order to obtain all the possible invariant solutions, the most systematic procedure is to determine the optimal system of one-dimensional subalgebras [10,11] for (13).We follow the approach given in [10] by firstly computing the commutators of the symmetry Lie algebra (12) and thereafter obtaining the adjoint representations (the calculations are summarized in Tables 2 and 3, resp.).We use Table 3 to simplify the linear combination of operators (12) given by for some constants  1 ,  2 , and  3 .Firstly, we let  3 ̸ = 0 (take  3 = 1).The operator (14) becomes We eliminate  2  2 by acting Ad(  2  2 ) on (15) and obtain for some constant   1 .Likewise, in order to eliminate   1  1 , we act on (16) by Ad( (  1 /) 1 ) to get Γ  =  3 .Next we let  3 = 0 ( 2 ̸ = 0).We take  2 = 1, and from (14) we have If  1 ̸ = 0 (i.e.,  1 > 0 or  1 < 0 ), then we obtain Γ =  1 +  2 , where  = ±1.
Next we utilize the optimal system to construct the invariant solutions of (13).However, the invariance under space translation, that is,  1 =   , is trivial hence, it is not considered.The other cases are presented as follows.
Case 1. Invariance under  2 : the corresponding characteristic system is given by the solution of which leads to the invariants Therefore, the invariant solution takes the form where () satisfies the reduced equation Now solving (21) for () and substituting into (20), we obtain the exact solution where  0 is an arbitrary constant.
where () is an arbitrary smooth function of its argument.Upon substitution of (23) into (13) and solving the resulting ordinary differential equation (ODE), we have for an arbitrary contant  0 .Therefore, we obtain the exact solution Case 3. Invariance under  3 : in this case, the invariant solution assumes the form where  = ( + )

Conclusion
In this work, the KdV equation with a forcing term was investigated using the Lie symmetry approach.The direct analysis of the classifying equation was employed to obtain the functional forms of the forcing term, which include power law and linear and constant time dependence.The threeand four-dimensional symmetry Lie algebras were obtained, respectively, for these forms of the forcing term.The optimal system of one-dimensional subalgebras of the Lie algebra of the invariant equation with the forcing term having the power law nonlinearity was obtained.As a result, for the same invariant equation, exact solutions were derived and the symmetry reduction was performed in the case where exact solutions were not obtained.