MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 845843 10.1155/2013/845843 845843 Research Article Lie Group Analysis of a Forced KdV Equation Molati Motlatsi 1, 2 0000-0002-1986-4859 Khalique Chaudry Masood 1 Mahomed Fazal M. 1 Department of Mathematical Sciences International Institute for Symmetry Analysis and Mathematical Modelling North-West University Mafikeng Campus Private Bag X 2046, Mmabatho 2735 South Africa nwu.ac.za 2 Department of Mathematics and Computer Science National University of Lesotho P.O. Roma 180 Lesotho nul.ls 2013 21 3 2013 2013 26 01 2013 26 02 2013 2013 Copyright © 2013 Motlatsi Molati and Chaudry Masood Khalique. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The Korteweg-de Vries (KdV) equation considered in this work contains a forcing term and is referred to as forced KdV equation in the sequel. This equation has been investigated recently as a mathematical model for waves on shallow water surfaces under the influence of external forcing. We employ the Lie group analysis approach to specify the time-dependent forcing term.

1. 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  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 (1)ut+cux+αuux+βuxxx=F(t), where u represents the elevation of free water surface and α and β are arbitrary constants which depend upon the long wave speed, c. The arbitrary smooth function F(t) is the forcing term. The variables t and x represent time and space, respectively. This equation was proposed in  and the Hirota's bilinear approach was used to obtain the multiple soliton solutions. However, it is worth mentioning that F(t) remained unspecified. In  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 ; 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.

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.

2. Lie Group Analysis

In Lie's algorithm (see  for more details), a vector field (2)X=ξ1(t,x,u)t+ξ2(t,x,u)x+η(t,x,u)u is a generator of Lie point symmetries of (1) if and only if (3)X(ut+cux+αuux+βuxxx-F(t))(1)=0, where (4)X=X+ζtut+ζxux+ζxxxuxxx is the third prolongation of the vector field X. The variables ζ's are given by the formulae (5)ζt=Dt(η)-utDt(ξ1)-uxDt(ξ2),ζx=Dx(η)-utDx(ξ1)-uxDx(ξ2),ζxxx=Dx(ζxx)-uxxtDx(ξ1)-uxxxDx(ξ2), where (6)Dt=t+utu+,Dx=x+uxu+ are the total derivative operators .

The invariance conditions (3) are separated with respect to the powers of the derivatives of u, 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 (7)ξu1=0,ξu2=0,ηuu=0,ξx1=0,ηxu-ξxx2=0,3ξx2-ξt1=0,ξ1F(t)+(ξt1-ηu)F(t)-(c+αu)ηx-βηxxx+ηt=0,(c+αu)(ξx2-ξt1)-ξt2+βξxxx2-αη-3βηxxu=0, where the subscripts denote partial differentiation with respect to the indicated variables, and a “prime” represents total derivative with respect to t.

The manipulation of (7) leads to the coefficients of Lie point symmetry generator given by (8)ξ1=k1t+k2,    ξ2=13k1x+a(t),      η=13α[3a(t)-2k1(c+αu)], where k1, k2 are arbitrary constants, and a(t) is an arbitrary function which satisfies the classifying relation (9)5αk1F(t)+3α(k1t+k2)F(t)=3a′′(t). Assume that F(t) is an arbitrary smooth function of t; then from the classifying (9) we obtain (10)k1=0,k2=0,a(t)=k3t+k4, where k3 and k4 are arbitrary constants. Thus, we have a two-dimensional principal Lie algebra, which is spanned by the operators (11)X1=x,X2=αtx+u. Our goal is to obtain the functional forms of the forcing term, F(t), for which the principal Lie algebra is extended. Therefore, the analysis of the classifying relation (9) considering the cases a′′(t)=0 and a′′(t)0 yields the various forms of F(t), 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  a more general case is considered.

Classification results.

No. F Condition on consts. Extension of principal Lie algebra
1. F 0 ( n + t ) - 5 / 3 + p 5 α α , F0, p0 X 3 = 6 α ( n + t ) t + α ( p t 2 + 2 x ) x + 2 [ p t - 2 ( c + α u ) ] u
2. F 0 + q t 3 α α , q0 X 3 = 6 α t + α q t 2 x + 2 q t u
3. F 0 ( n + t ) - 5 / 3 F 0 0 X 3 = 3 α ( n + t ) t + α x x - 2 ( c + α u ) u
4. F 0 F 0 0 X 3 = t ,
X 4 = 6 α t + α ( 5 α F 0 t 2 + 2 x ) x + 2 [ 5 α F 0 t - 2 ( c + α u ) ] u

Here F0, n, p and q are arbitrary constants.

3. 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 four-dimensional 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 (12)X1=x,X2=αtx+u,      X3=3α(n+t)t+αxx-2(c+αu)u, which are the symmetries of the equation (13)ut+cux+αuux+βuxxx=F0(n+t)-5/3. 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  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.).

Table of commutators.

[ X i , X j ]    X 1    X 2    X 3
X 1 0 0 α X 1
X 2 0 0 - 3 n α 2 X 1 - 2 α X 2
X 3 - α X 1 3 n α 2 X 1 + 2 α X 2 0

Here [Xi,Xj]=Xi(Xj)-Xj(Xi); i,j=1,2,3 is the commutator operation.

Ad ( e ϵ X i ) X j X 1    X 2    X 3
X 1 X 1 X 2 X 3 - α ϵ X 1
X 2 X 1 X 2 X 3 + ϵ ( 3 n α 2 X 1 + 2 α X 2 )
X 3 e ϵ α X 1 e - 2 α ϵ X 2 + 3 n α ( e - ϵ α - 1 ) X 1 X 3

Here Ad(eϵXi)Xj=Xj-ϵ[Xi,Xj]+(1/2!)ϵ2[Xi,[Xi,Xj]]- is the adjoint representation where ϵ is a real number.

We use Table 3 to simplify the linear combination of operators (12) given by (14)Γ=a1X1+a2X2+a3X3 for some constants a1, a2, and a3.

Firstly, we let a30 (take a3=1). The operator (14) becomes (15)Γ=a1X1+a2X2+X3. We eliminate a2X2 by acting Ad(ea2X2) on (15) and obtain (16)Γ=a1X1+X3 for some constant a1. Likewise, in order to eliminate a1X1, we act on (16) by Ad(e(a1/α)X1) to get Γ′′=X3.

Next we let a3=0 (a20). We take a2=1, and from (14) we have (17)Γ=a1X1+X2. If a10 (i.e., a1>0 or a1<0 ), then we obtain Γ=λX1+X2, where λ=±1.

Finally, we let a2=a3=0. The operator (14) reduces to Γ=X1 for a1=1.

Therefore, an optimal system of one-dimensional subalgebras is given by {X1,X2,λX1+X2,X3}.

Next we utilize the optimal system to construct the invariant solutions of (13). However, the invariance under space translation, that is, X1=x, is trivial hence, it is not considered. The other cases are presented as follows.

Case 1.

Invariance under X2: the corresponding characteristic system is given by (18)dt0=dxαt=du1, the solution of which leads to the invariants (19)C1=t,C2=u-xαt. Therefore, the invariant solution takes the form (20)u(t,x)=f(t)+xαt, where f(t) satisfies the reduced equation (21)cαt-F0(n+t)5/3+ft+dfdt=0. Now solving (21) for f(t) and substituting into (20), we obtain the exact solution (22)u(t,x)=1t[f0+F0{3+3n2(n+t)}(n+t)1/3+x-cα], where f0 is an arbitrary constant.

Case 2.

The invariance under λX1+X2(λ0) yields the invariant solution (23)u(t,x)=f(t)-xαt+λ, where f(t) is an arbitrary smooth function of its argument. Upon substitution of (23) into (13) and solving the resulting ordinary differential equation (ODE), we have (24)f(t)=1αt+λ{f0-ct+3F0(λ-3nα-2tα)2(n+t)2/3}, for an arbitrary contant f0. Therefore, we obtain the exact solution (25)u(t,x)=1αt+λ{f0-ct-x+3F0(λ-3nα-2tα)2(n+t)2/3}.

Case 3.

Invariance under X3: in this case, the invariant solution assumes the form (26)u(t,x)=f(z)(n+t)2/3-cα, where z=x(n+t)-1/3 is the similarity variable. The function f(z) is an arbitrary function which satisfies the third-order ODE (27)3βd3fdz3+(3αdfdz+2)f+zdfdz-3F0=0.

4. 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 three-and 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.

Acknowledgment

M. Molati thanks the North-West University, Mafikeng Campus, for the Postdoctoral Fellowship.

Korteweg D. J. de Vries G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves Philosophical Magazine 1895 39 240 422 443 10.1080/14786449508620739 Salas A. H. Computing solutions to a forced KdV equation Nonlinear Analysis: Real World Applications 2011 12 2 1314 1320 10.1016/j.nonrwa.2010.09.028 MR2736312 ZBL1206.35222 Wazwaz A.-M. A study on KdV and Gardner equations with time-dependent coefficients and forcing terms Applied Mathematics and Computation 2010 217 5 2277 2281 10.1016/j.amc.2010.06.038 MR2727974 ZBL1202.35266 Jefferson G. F. On the second-order approximate symmetry classification and optimal systems of subalgebras for a forced Korteweg–de Vries equation Communications in Nonlinear Science and Numerical Simulation 2013 10.1016/j.cnsns.2012.12.022 Gandarias M. L. Bruzón M. S. Some conservation laws for a forced KdV equation Nonlinear Analysis: Real World Applications 2012 13 6 2692 2700 10.1016/j.nonrwa.2012.03.013 MR2927217 Listopadova V. Magda O. Pobyzh V. How to find solutions, Lie symmetries, and conservation laws of forced Korteweg-de Vries equations in optimal way Nonlinear Analysis: Real World Applications 2013 14 1 202 205 10.1016/j.nonrwa.2012.05.013 MR2969829 ZBL1253.35151 Johnpillai A. G. Khalique C. M. Group analysis of KdV equation with time dependent coefficients Applied Mathematics and Computation 2010 216 12 3761 3771 10.1016/j.amc.2010.05.043 MR2661741 ZBL1197.35231 Molati M. Ramollo M. P. Symmetry classification of the Gardner equation with time-dependent coefficients arising in stratified fluids Communications in Nonlinear Science and Numerical Simulation 2012 17 4 1542 1548 10.1016/j.cnsns.2011.09.002 MR2855446 ZBL1245.35110 Bluman G. W. Kumei S. Symmetries and Differential Equations 1989 81 New York, NY, USA Springer MR1006433 ZBL0698.35001 Olver P. J. Applications of Lie Groups to Differential Equations 1986 107 New York, NY, USA Springer 10.1007/978-1-4684-0274-2 MR836734 Ovsiannikov L. V. Group Analysis of Differential Equations 1982 New York, NY, USA Academic Press MR668703