1. Introduction
The generalized coupled KdV system given by [1]
(1)ut+auxxx-buux+cvvx=0,vt+dvxxx-euvx+fuxv=0,
where a, b, c, d, e, and f are real constants, describes the interaction of two long waves, whose dispersion relations are different. For the case when f=0, soliton solutions have been obtained in [2, 3]. Many other special cases of (1) have been considered in the literature, and various methods have been used to find its exact solutions. See, for example, [4–11].
In this study, we consider a special case of the generalized coupled KdV system given by
(2)ut+auxxx+buux+cvvx=0,vt+dvxxx+cuvx+cuxv=0
and construct conservation laws for (2). Recently, the conservation laws of system (2) for special values of the constants a=d=-1 and b=c=-6 were derived in [12] using the multiplier approach.
Many nonlinear partial differential equations (PDEs) of mathematical physics and engineering are continuity equations, which express conservation of mass, momentum, energy, or electric charge. It is well known that conservation laws play a crucial role in the solution and reduction of PDEs. For variational problems the conservation laws can be constructed by means of the Noether theorem [13]. The application of the Noether theorem depends upon the existence of a Lagrangian. However, there are nonlinear differential equations that do not have a Lagrangian. In such instances, researchers have developed several methods to derive conserved quantities for such equations. See, for example, [14–20].
The organization of this paper is as follows. In Section 2 we briefly recall some notations and fundamental relations concerning the Noether symmetries approach, which we utilize in the same section to obtain the Noether symmetries and the corresponding conserved vectors. The concluding remarks are summarized in Section 3.
2. Conservation Laws of Coupled KdV Equations
In this section we derive the conservation laws for the generalized coupled KdV system (2). This system does not have a Lagrangian. In order to apply the Noether theorem we transform our system (2) to a fourth-order system, using the transformations u=Ux and v=Vx. Then system (2) becomes
(3)Utx+aUxxxx+bUxUxx+cVxVxx=0,Vtx+dVxxxx+cUxVxx+cVxUxx=0.
It can readily be verified that the second-order Lagrangian for system (3) is given by
(4)L=12(aUxx2+dVxx2-13bUx3-cUxVx2-UxUt-VtVx)
because
(5)δLδU=0, δLδV=0,
where δ/δU and δ/δV are the standard Euler operators defined by
(6)δδU=∂∂U-Dt∂∂Ut-Dx∂∂Ux+Dt2∂∂Utt+Dx2∂∂Uxx+DxDt∂∂Utx-⋯,δδV=∂∂V-Dt∂∂Vt-Dx∂∂Vx+Dt2∂∂Vtt+Dx2∂∂Vxx+DxDt∂∂Vtx-⋯.
Consider the vector field
(7)X=ξ1(t,x,U,V)∂∂t+ξ2(t,x,U,V)∂∂x+η1(t,x,U,V)∂∂U+η2(t,x,U,V)∂∂V,
which has the second-order prolongation defined by
(8)X[2]=ξ1(t,x,U,V)∂∂t+ξ2(t,x,U,V)∂∂x+η1(t,x,U,V)∂∂U+η2(t,x,U,V)∂∂V+ζt1∂∂Ut+ζt2∂∂Vt+ζx1∂∂Ux+ζx2∂∂Vx+⋯.
Here
(9)ζt1=Dt(η1)-UtDt(ξ1)-UxDt(ξ2), ζx1=Dx(η1)-UtDx(ξ1)-UxDx(ξ2), ζt2=Dt(η2)-VtDt(ξ1)-VxDt(ξ2), ζx2=Dx(η2)-VtDx(ξ1)-VxDx(ξ2),Dt=∂∂t+Ut∂∂U+Vt∂∂V+Utt∂∂Ut+Vtt∂∂Vt+Utx∂∂Ux+Vtx∂∂Vx+⋯,Dx=∂∂x+Ux∂∂U+Vx∂∂V+Uxx∂∂Ux+Vxx∂∂Vx+Utx∂∂Ut+Vtx∂∂Vt+⋯.
The Lie-Bäcklund operator X defined in (7) is a Noether operator corresponding to the Lagrangian (4) if it satisfies
(10)X[2](L)+L[Dt(ξ1)+Dx(ξ2)]=Dt(B1)+Dx(B2),
where B1(t,x,U,V), B2(t,x,U,V) are the gauge terms. Expansion of (10) yields
(11)-12Ux[ηt1+UtηU1+VtηV1-Utξt1-Ut2ξU1-12Ux -UtVtξV1-Uxξt2-UtUtξU2-UxVtξV2] -12Vx[ηt2+UtηU2+VtηV2-Vtξt1-UtVtξU1 -12Vx -Vt2ξV1-Vxξt2-UtVxξU2-VtVxξV2] -12(bUx2+cVx2+Ut) ×[ηx1+UxηU1+VxηV1-Utξx1-UtUxξU1 -UtVxξV1-Uxξx2-Ux2ξU2-UxVxξV2] -12(cUxVx+Vt) ×[ηx2+UxηU2+VxηV2-Vtξx1-UxVtξU1 -VtVxξV1-Vxξx2-UxVxξU2-Vx2ξV2] +dVxx[Dx2(η1)-UtDx2(ξ1)-UxDx2(ξ2) +dVxx -2Utx(ξx1+UxξU1+VxξV1) +dVxx -2Uxx(ξx2+UxξU2+VxξV2)] +aUxx[Dx2(η2)-VtDx2(ξ1)-VxDx2(ξ2) +aUxx -2Vtx(ξx1+UxξU1+VxξV1) +aUxx -2Vxx(ξx2+UxξU2+VxξV2)] +12(aUxx2+dVxx2-13bUx3-cUxVx2-UxUt-VtVx) ×[ξt1+UtξU1+VtξV1+ξx2+UxξU2+VxξV2] =Bt1+UtBU1+VtBV1+Bx2+UxBU2+VxBV2.
The splitting of (11) with respect to different combinations of derivatives of U and V results in an overdetermined system of PDEs for ξ1,ξ2,η1,η2,B1, and B2. Solving this system of PDEs we arrive at the following two cases for which Noether symmetries exist.
Case 1.
b
≠
c
.
In this case we obtain the following Noether symmetries and gauge terms:
(12)ξ1=A1,ξ2=A2,η1=E(t),η2=F(t),B1=P(t,x),B2=-12UE′(t)-12VF′(t)+S(t,x),Pt+Sx=0.
The above results will now be used to find the components of the conserved vectors for the second-order Lagrangian. Here we can choose P=0, S=0 as they contribute to the trivial part of the conserved vector. We recall that the conserved vectors for the second-order Lagrangian are given by [13, 21]
(13)T1=-B1+ξ1L+W1[∂L∂Ut-Dt∂L∂Utt-Dx∂L∂Utx⋯]+W2[∂L∂Vt-Dt∂L∂Vxt-Dx∂L∂Vtt⋯]+Dt(W1)∂L∂Utt+Dt(W2)∂L∂Vtt,T2=-B2+ξ2L+W1[∂L∂Ux-Dt∂L∂Uxt-Dx∂L∂Uxx⋯]+W2[∂L∂Vx-Dt∂L∂Vxt-Dx∂L∂Vxx⋯]+Dx(W1)∂L∂Uxx+Dx(W2)∂L∂Vxx.
Here W1 and W2 are the Lie characteristic functions, given by W1=η1-Utξ1-Uxξ2 and W2=η2-Vtξ1-Vxξ2. Using (13) together with (12) and u=Ux, v=Vx we obtain the following independent conserved vectors for system (2):
(14)T11=12(aux2+dvx2-13bu3-cuv2),T12=12∫utdx∫utdx+12(bu2+cv2)×∫utdx+auxx∫utdx+12∫vtdx∫vtdx+dvxx∫vtdx+cuv∫vtdx-autux-dvtux,(15)T21=12(u2+v2),T22=auuxx+dvvxx-12aux2-12dvx2+13bu3+cuv2,
and for the arbitrary functions E(t) and F(t),
(16)T(E,F)1=-12uE(t)-12vF(t),T(E,F)2=12E′(t)∫udx+12F′(t)∫vdx-12E(t)∫utdx-12F(t)∫vtdx-12(bu2+cv2)E(t)-auxxE(t)-dvxxF(t)-cuvF(t).
Conserved vector (14) is a nonlocal conserved vector, and (15) is a local conserved vector for system (2). We now derive two particular cases from conserved vector (16) by letting E(t)=1 and F(t)=0, which gives a nonlocal conserved vector
(17)T31=-12u,T(3)2=-12(bu2+cv2)-auxx-12∫utdx,
and by choosing E(t)=0 and F(t)=1, we get the nonlocal conserved vector
(18)T41=-12v,T42=-cuv-dvxx-12∫vtdx.
Case 2.
b
=
c
.
The second case gives the following Noether symmetries and gauge terms:
(19)ξ1=A1,ξ2=cA2t+A3,η1=A2x+F(t),η2=G(t),B1=-12A2U+P(t,x),B2=-12UF′(t)-12VG′(t)+R(t,x),Pt+Rx=0.
Again we can set P=0 and R=0 as they contribute to the trivial part of the conserved vector. The independent conserved vectors for system (2), in this case, are
(20)T11=12(aux2+dvx2-13bu3-cuv2),T12=12∫utdx∫utdx+12(bu2+cv2)∫utdx+auxx∫utdx+12∫vtdx∫vtdx+dvxx∫vtdx+cuv∫vtdx-autux-dvtux,T21=12(-xu+ctu2+ctv2+∫udx),T22=aux+actuuxx+cdtvvxx+c2tuv2-axuxx+13cbtu3-12(actux2+cdtvx2+bxu2+cxv2+x∫utdx),(21)T31=12(u2+v2),T32=auuxx+dvvxx-12aux2-12dvx2+13bu3+cuv2,
and for the arbitrary functions E(t) and F(t), we obtain
(22)T(E,F)1=-12uE(t)-12vF(t),T(E,F)2=12E′(t)∫udx+12F′(t)∫vdx-12E(t)∫utdx-12F(t)∫vtdx-12(bu2+cv2)E(t)-auxxE(t)-dvxxF(t)-cuvF(t).
Conserved vectors (20) are nonlocal, whereas (21) is a local conserved vector for system (2). Conserved vector (22) for E(t)=1 and F(t)=0 gives a nonlocal conserved vector
(23)T31=-12u,T32=-12(bu2+cv2)-auxx-12∫utdx,
and for E(t)=0 and F(t)=1 it gives a nonlocal conserved vector
(24)T41=-12v,T42=-cuv-dvxx-12∫vtdx.
We note that for arbitrary values of E(t) and F(t) infinitely many nonlocal conservation laws exist for system (2).