We employ the multiplier approach (variational derivative method) to derive the conservation laws for the Degasperis Procesi equation and a coupled variable-coefficient modified Korteweg-de Vries system in a two-layer fluid model. Firstly, the multipliers are computed and then conserved vectors are obtained for each multiplier.
1. Introduction
The conservation laws are important in the solution and reductions of partial differential equations. Conservation law, also called law of conservation, in physics, is several principles that state that certain physical properties (i.e., measurable quantities) do not change in the course of time within an isolated physical system. In classical physics, laws of this type govern energy, momentum, angular momentum, mass, and electric charge. In particle physics, other conservation laws apply to properties of subatomic particles that are invariant during interactions. An important function of conservation laws is that they make it possible to predict the macroscopic behaviour of a system without having to consider the microscopic details of the course of a physical process or chemical reaction. Many powerful methods have been developed for the construction of conservation laws, such as The Laplace Direct method [1], multiplier approach [2, 3], Kara and Mahomed symmetry condition [4], Wolf [5, 6], Göktaş and Hereman [7], Hereman et al. [8–10], and Cheviakov [11] who developed powerful software packages to compute conservation laws for partial differential equations. Infinitely many conservation laws are obtained based on the Lax pair via the Hirota method and symbolic computation, bilinear forms, bilinear Backlund transformations, and one- and two-soliton-like solutions are also derived. With different coefficients, bell-shaped, periodic-changing, quadratic-varying, exponential-decreasing, and exponential-increasing soliton-like profiles are obtained in [12]. Also by the spectral analysis the Hamiltonian and periodicity of the qZK equation are investigated by usig the Hirota method [13]. The nonautonomous matter waves with time-dependent modulation in a one-dimensional trapped spin-1 Bose-Einstein condensate and the generalized three-coupled Gross-Pitaevskii equations by means of the Hirota bilinear method are studied in [14]. The multiplier approach (also known as variational derivative method) was proposed by Steudel [15] who wrote the conservation law in characteristic form as DiTi=AαEα. Later, Olver [16] modified the method of determining the characteristics (multipliers) by taking the variational derivative of DiTi=QαEα not only for the arbitrary functions, but also for solutions of system of partial differential equations. The outline of the paper is as follows. In Section 2, some definitions related to the multiplier approach are given. In Section 3, conservation laws for the Degasperis Procesi equation are derived by first computing the multipliers. The conservation laws for a coupled variable-coefficient modified Korteweg-de Vries system in a two-layer fluid model are derived in Section 4. Finally, conclusions are summarized in Section 5.
2. Necessary Preliminaries
Let xi, i=1,2,…,n, be n independent variables and let u be the dependent variable.
The total derivative operator with respect to xi is
(1)Di=∂∂xi+ui∂∂u+uij∂∂uj+⋯i=1,2,…,n,
where ui denotes the derivative of u with respect to xi. Similarly uij denotes the derivative of u with respect to xi and xj.
The Euler operator is defined by
(2)δδu=∂∂u-Di∂∂ui+Dij∂∂uij-Dijk∂∂uijk+⋯.
Consider a kth-order partial differential equation of n independent and one dependent variable
(3)E(x,u,u(1),u(2),…,u(k))=0.
An n-tuple T=(T1,T2,…,Tn), i=1,2,…,n, such that
(4)DiTi=0
holds for all solutions of (3) is known as the conserved vector of (3).
The multiplier A of system (3) has the property
(5)DiTi=AE
for arbitrary function u(x1,x2,…,xn) [16].
The determining equations for multipliers are obtained by taking the variational derivative of (5) (see [16]):
(6)δδu(AE)=0.
Equation (6) holds for arbitrary function u(x1,x2,…,xn) not only for solutions of system (3).
Once the multipliers are computed from (6), the conserved vectors can be derived systematically using (5) as the determining equation. But in some problems it is not difficult to construct the conserved vectors by elementary manipulations after the determination of the multipliers.
3. Conservation Laws for the Degasperis Procesi Equation
The Degasperis Procesi equation [17] takes the form
(7)ut-uxxt+4uux-3uxuxx-uuxxx=0.
The Degasperis Procesi equation (7) is very interesting as it is an integrable shallow water equation and presents a quite rich structure. Also it can be used to model wave perturbations in relaxing media.
We will derive the conservation laws for (7) by the multiplier approach. The determining equation for multiplier A(t,x,u), from (6), is
(8)δδuAut-uxxt+4uux-3uxuxx-uuxxx=0.
The standard Euler operator δ/δu from (2) can be defined as
(9)δδu=∂∂u-Dt∂∂ut-Dx∂∂ux+Dt2∂∂utt+Dx2∂∂uxx+DxDt∂∂utx-⋯,
and total derivative operators Dt and Dx using (1) are
(10)Dt=∂∂t+ut∂∂u+utt∂∂ut+utx∂∂ux+⋯,Dx=∂∂t+ux∂∂u+uxx∂∂ux+utx∂∂ut+⋯.
Equation (8) after expansion and simplification takes the following form:
(11)uxx(Aut+3uAxu+3uuxAuu-3uxAu+utAuu)WW+ux2(uAuux+Auut+utAuuu)+uux3AuuuWW+uxt2Axu+2uxAuu+ux2Axtu+3uAxxuWW+ut(Axxu+2uxAxuu)+u(Axxx-4Ax)WW+Axxt-At=0,
which yields
(12)A=c1+c2e-2x+c3e2x.
From (5) and (12), we have
(13)(c1+c2e-2x+c3e2x)×(ut-uxxt+4uux-3uxuxx-uuxxx)=Dtc1(u-uxx)+c2(ue-2x-uxxe-2x)WWWW+c3ue2x-uxxe2x+Dxc1(2u2-ux2-uuxx)WWWWiii+c2(-2uuxe-2x-ux2e-2x-uuxxe-2x)WWWWiii+c3(2uuxe2x-ux2e2x-uuxxe2x)
for arbitrary functions u(t,x). When u(t,x) is solutions of (7) then left hand side of (13) vanishes and we obtain
(14)Dtc1(u-uxx)+c2(ue-2x-uxxe-2x)iiiiiii+c3(ue2x-uxxe2x)+Dxc1(2u2-ux2-uuxx)WWWiii+c2(-2uuxe-2x-ux2e-2x-uuxxe-2x)WWWiii+c3(2uuxe2x-ux2e2x-uuxxe2x)=0.
Therefore the conserved vectors for the Degasperis Procesi equation (7) are
(15)T11=u-uxx,T12=2u2-ux2-uuxx,T21=ue-2x-uxxe-2x,T22=-2uuxe-2x-ux2e-2x-uuxxe-2x,T31=ue2x-uxxe2x,T32=2uuxe2x-ux2e2x-uuxxe2x.
The variational derivative approach for the Degasperis Procesi equation gives three multipliers of the form A(t,x,u) and hence three conserved vectors are obtained.
4. Conservation Laws for a Coupled Variable-Coefficient Modified Korteweg-de Vries System in a Two-Layer Fluid Model
In this section we recall some basic definitions related to the multiplier approach.
Let (x,t) be two independent variables and let (u,v) be dependent variables.
The total derivative operators Dt and Dx are
(16)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 standard Euler operators δ/δu and δ/δv are
(17)δδu=∂∂u-Dt∂∂ut-Dx∂∂ux+Dt2∂∂utt+Dx2∂∂uxx+DxDt∂∂utx-⋯,(18)δδv=∂∂v-Dt∂∂vt-Dx∂∂vx+Dt2∂∂vtt+Dx2∂∂vxx+DxDt∂∂vtx-⋯.
Consider a kth-order system of two partial differential equations of two independent and two dependent variables
(19)E1t,x,u,v,ut,vt,…,E2t,x,u,v,ut,vt,….
A vector T=(T1,T2) satisfying
(20)DtT1+DxT2=0
for all solutions of (19) is known as the conserved vector of (19).
The multipliers A1,A2 of system (19) have the property
(21)DtT1+DxT2=A1E1+A2E2,
for the arbitrary functions u(x,t),v(x,t).
The determining equations for the multipliers are obtained by taking variational derivative of (21):
(22)δδuA1E1+A2E2=0,δδvA1E1+A2E2=0.
Equation (22) holds for the arbitrary functions u(x,t),v(x,t) not only for the solutions of system (19). Equation (22) yields multipliers for all local conservation laws. Then conserved vectors can be derived systematically using (21) as the determining equation. But in some problems it is not difficult to construct the conserved vectors by elementary manipulations once the multiplier has been determined.
Example 1.
Consider a coupled variable-coefficient modified Korteweg-de Vries system in a two-layer fluid model
(23)ut-αtuxxx+6u2-v2ux-12uvvx-4β(t)ux=0,vt-α(t)[vxxx+6(u2-v2)vx+12uvux]-4β(t)vx=0.
System (23) was proposed in [18] as an important particular case of the formidable generalized coupled variable-coefficient modified Korteweg-de Vries (CVmKdV) system. The (CVmKdV) system was derived by Gao and Tang [19] as a two-layer model describing atmospheric and oceanic phenomena like interactions between the atmosphere and ocean, atmospheric blocking, oceanic circulations, hurricanes, typhoons, and so forth.
The determining equations for multipliers of the forms A1(x,t,u,v) and A2(x,t,u,v) from (22) are
(24)δδu(u2-v2)A1ut-α(t)[uxxx+6(u2-v2)ux-12uvvx]iiiiiiiiiiiii-4β(t)ux(u2-v2)iiiiii+A2vt-α(t)[vxxx+6(u2-v2)vx+12uvux]iiiiiiiiiiiiiiii-4β(t)vx(u2-v2)(u2-v2)=0,δδv(u2-v2)A1ut-α(t)[uxxx+6(u2-v2)ux-12uvvx]iiiiiiiiiiiiii-4β(t)ux(u2-v2)iiiiiii+A2vt-α(t)[vxxx+6(u2-v2)vx+12uvux]iiiiiiiiiiiiiiiii-4β(t)vx(u2-v2)(u2-v2)=0,
where the standard Euler operators δ/δu and δ/δv are defined in (17) and (18), respectively. Expansion of (24) yields
(25)A1uut-α(t)[uxxx+6(u2-v2)ux-12uvvx]iiiiiii-4β(t)ux(u2-v2)+A2uvt-α(t)[vxxx+6(u2-v2)vx+12uvux]iiiiiiiiiiiiiii-4β(t)vx(u2-v2)+12αvvx-12αuuxA1-12αuvx+12αvuxA2-Dt(A1)+Dx6αu2-v2+4βA1+12αuvA2+Dx3αA1,(26)A1vut-α(t)[uxxx+6(u2-v2)ux-12uvvx]iiiiii-4β(t)ux(u2-v2)+A2vvt-α(t)[vxxx+6(u2-v2)vx+12uvux]iiiiiiiiiiiiii-4β(t)vx(u2-v2)+12αvvx-12αuuxA2+12αuvx+12αvuxA1-Dt(A2)+Dx[(6α(u2-v2)+4β)A2-12αuvA1]+Dx3αA2.
Equations (25) and (26) are separated according to different combinations of derivatives of u and v and, after some simplification following system of equations for A1,A2 is obtained:
(27)A1uu=0,A1vv=0,A1ux=0,A1vx=0,A1uv=0,A1v-A2u=0,A1u+A2v=0,6αu2-v2+4βA1x+12αuvA2x-A1t+αA1xxx=0,A2uu=0,A2vv=0,A2ux=0,A2vx=0,A2uv=0,6αu2-v2+4βA2x-12αuvA1x-A2t+αA2xxx=0.
Solution of system (27) yields
(28)A1=c1+c3u+c4v,A2=c2-c3v+c4u,
where c1,c2,c3, and c4 are constants.
Equations (21) and (28) give the following conserved vectors satisfying (20):(29)T11=u,T12=-2αu3-αuxx+6αuv2-4βu,T21=v,T22=2αv3-αvxx-6αu2v-4βv,T31=u22-v22,T32=αvvxx-uuxx+9αu2v2+αux22-vx22+v22β-32αv2-u232αu2+2β,T41=uv,T42=-6αu3v+αuxvx-vuxx+u6αv3-4vβ-αvxx.
5. Conclusions
The conservation laws for the Degasperis Procesi equation and a coupled variable-coefficient modified Korteweg-de Vries system in a two-layer fluid model were established with the help of the multiplier approach. The multiplier approach on the Degasperis Procesi equation yielded three multipliers and thus three local conserved vectors were obtained in each case. The multiplier approach when applied to a coupled variable-coefficient modified Korteweg-de Vries system in a two-layer fluid model gave four multipliers of form A(x,t,u,v). Each multiplier corresponds to a conserved vector and thus four local conserved vectors were obtained.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
LaplaceP. S.1966New York, NY, USACelestrial MechanicsOlverP. J.1993New York, NY, USASpringerMR124005610.1007/978-1-4612-4350-2AncoS. C.BlumanG.Direct construction method for conservation laws of partial differential equations. Part I: examples of conservation law classifications200213554556610.1017/S0956792501004661MR19391602-s2.0-0036815244KaraA. H.MahomedF. M.Relationship between symmetries and conservation laws2000391234010.1023/A:1003686831523MR17550152-s2.0-0034378632WolfT.A comparison of four approaches to the calculation of conservation laws2002132129152MR18962252-s2.0-003601371110.1017/S0956792501004715WolfT.BrandA.MohammadzadehM.Computer algebra algorithms and routines for the computation of conservation laws and fixing of gauge in differential expressions199927222123810.1006/jsco.1998.0250MR1672132ZBL0919.650762-s2.0-0000070068GöktaşÜ.HeremanW.Symbolic computation of conserved densities for systems of nonlinear evolution equations199724559162110.1006/jsco.1997.0154MR14840712-s2.0-0031281728HeremanW.AdamsP. J.EklundH. L.HickmanM. S.HerbstB. M.YanZ.Direct methods and symbolic software for conservation laws of nonlinear equations2009chapter 2New York, NY, USANova Science Publishers1979HeremanW.ColagrossoM.SayersR.RinglerA.DeconinckB.NivalaM.HickmanM.WangD.ZhengZ.Continuous and discrete homotopy operators and the computation of conservation laws2005Birkhäuser249285MR218738410.1007/3-7643-7429-2_15HeremanW.Symbolic computation of conservation laws of nonlinear partial differential equations in multi-dimensions2006106127829910.1002/qua.20727ZBL1188.683642-s2.0-30544436592CheviakovA. F.GeM software package for computation of symmetries and conservation laws of differential equations20071761486110.1016/j.cpc.2006.08.001MR22794592-s2.0-33751539441WangY.-F.TianB.LiM.WangP.WangM.Integrability and soliton-like solutions for the coupled higher-order nonlinear Schrödinger equations with variable coefficients in inhomogeneous optical fibers201419617831791MR314475910.1016/j.cnsns.2013.09.0032-s2.0-84890789083Hui-LingZ.BoT.Yu-FengW.HuiZ. H.Wen-RongS.Dynamic behavior of the quantum Zakharov-Kuznetsov equations in dense quantum magnetoplasmas201421012304ShenY.-J.GaoY.-T.ZuoD.-W.SunY.-H.FengY.-J.XueL.Nonautonomous matter waves in a spin-1 Bose-Einstein condensate20148961206291510.1103/PhysRevE.89.062915SteudelH.Über die Zuordnung zwischen Invarianzeigenschaften und Erhaltungssatzen196217A129132MR0156591OlverP. J.1993New York, NY, USASpringer10.1007/978-1-4612-4350-2MR1240056DegasperisA.ProcesiM.DegasperisA.GaetaG.Asymptotic integrability1999SingaporeWorld Scientific Publishers2327MR1844104ZhuS.-H.GaoY.-T.YuX.SunZ.-Y.GaiX.-L.MengD.-X.Painlevé property, soliton-like solutions and complexitons for a coupled variable-coefficient modified Korteweg-de Vries system in a two-layer fluid model2010217129530710.1016/j.amc.2010.05.061MR26725892-s2.0-77955416661GaoY.TangX.-Y.A coupled variable coefficient modified KdV equation arising from a two-layer fluid system200748696197010.1088/0253-6102/48/6/001MR2445636ZBL1267.351942-s2.0-38549092717