arxiv: v3 [physics.flu-dyn] 11 Sep 2015

Transcription

1 Energy invariant for shallow water waves and the Korteweg de Vries equation. Is energy always an invariant? Anna Karczewska Faculty of Mathematics, Computer Science and Econometrics University of Zielona Góra, Szafrana 4a, Zielona Góra, Poland arxiv:53.989v3 physics.flu-dyn] Sep 5 Piotr Rozmej Institute of Physics, Faculty of Physics and Astronomy University of Zielona Góra, Szafrana 4a, Zielona Góra, Poland Eryk Infeld National Centre for Nuclear Research, Hoża 69, -68 Warszawa, Poland Dated: September 4, 5 It is well known that the KdV equation has an infinite set of conserved quantities. The first three are often considered to represent mass, momentum and energy. Here we try to answer the question of how this comes about, and also how these KdV quantities relate to those of the Euler shallow water equations. Here Luke s Lagrangian is helpful. We also consider higher order extensions of KdV. Though in general not integrable, in some sense they are almost so, these with the accuracy of the expansion. PACS numbers:.3.jr, a, Bb, Fg Keywords: Soliton, shallow water waves, nonlinear equations, invariants of KdV I. INTRODUCTION There exists a vast number of papers dealing with the shallow water problem. Aspects of the propagation of weakly nonlinear, dispersive waves are still beeing studied. Last year we published two articles, ] in which Korteveg de Vries type equations were derived in weakly nonlinear, dispersive and long wavelength limit. The second order KdV type equation was derived. The second order KdV equation 3, 4], sometimes called "extended KdV equation", was obtained for the case with a flat bottom. In derivation of the new equation we adapted the method described in 4]. In ], an analytic solution of this equation in the form of a particular soliton was found, as well. It is well known, see, e.g. 5 8], that for the KdV equation there exists an infinite number of invariants, that is, integrals over space of functions of the wave profile and its derivatives, which are constants in time. Looking for analogous invariants for the second order KdV equation we met with some problems even for the standard KdV equation which is first order in small parameters. This problem appears when energy conservation is considered. In this paper we reconsider invariants of the KdV equation and formulas for the total energy in several different approaches and different frames of reference fixed and moving ones. We find that the invariant I 3, sometimes called the energy invariant, does not always have A.Karczewska@wmie.uz.zgora.pl P.Rozmej@if.uz.zgora.pl Eryk.Infeld@ncbj.gov.pl that interpretation. We also give a proof that for the second order KdV equation, obtained in 4], η dx is not an invariant of motion. There are many papers considering higher-order KdV type equations. Among them we would like to point out works of Byatt-Smith 9], Kichenassamy and Olver ], Marchant 3, 4], Zou and Su 5], Tzirtzilakis et.al. 6] and Burde 7]. It was shown that if some coefficients of the second order equation for shallow water problem are diferent or zero then there exists a hierarchy of solition solutions. Kichenassamy and Olver ] even claimed that for second order KdV equation solitary solutions of appropriate form can not exist. This claim was falsified in our paper ] where the analytic solution of the second order KdV equation was found. Concerning the energy conservation there are indications that collisions of solitons 8, 9] which are solutions of higher order equations of KdV type can be inelastic 5, 6]. The paper is organized as follows. In Section II several frequently used forms of KdV equations are recalled with particular attention to transformations between fixed and moving reference frames. In Section III the form of the three lowest invariants of KdV equations is derived for different forms of the equations. In Section IV we show that the energy calculated from the definition H = T+V has no invariant form. Section V describes the variational approach in a potential formulation which gives a proper KdV equation but fails in obtaining second order KdV equations. In the next section the proper invariants are obtained from Luke s Lagrangian density. Section VII summarizes conclusions on the energy for KdV equation. In section VIII we apply the same formalism to calculate energy for waves governed by the extended KdV equation second order. We found that energy is

2 not conserved neither in fixed coordinate system nor in the moving frame. II. THE EXTENDED KDV EQUATION The geometry of shallow water waves is presented in Fig.. In, ] we derived an equation, second order in small parameters, in the fixed reference system and with scaled nondimensional variables containing terms for bottom fluctuations. They will not be considered here. ηx,t α=a/h β=h/l a bottom undisturbed surface ηx,t FIG.. Schematic view of the geometry. For a flat bottom that equation reduces to the second order KdV type equation, identical with 4, Eq. ] for β = α, that is, η t +η x +α 3 ηη x +β 6 η 3x +α 3 8 η η x 3 +αβ 4 η xη x + 5 ηη 3x +β 9 36 η 5x =. Subscripts denote partial differentiation. Small parameters α,β are defined by ratios of the wave amplitude a, the average water depth h and mean wavelength l α = a h, β = h l. Equation was earlier derived in 3] and called "the extended KdV equation". Limitation to the first order in small parameters yields the KdV equation in a fixed coordinate system η t +η x +α 3 ηη x +β 6 η 3x =. Transformation to a moving frame in the form x = x t, t = t, η = η, 3 allows us to remove the term η x in the KdV equation in a frame moving withthe velocity of sound gh η t +α 3 η η x +β 6 η 3 x =. 4 h The explicit form of the scaling leading to equations 4 is given by 9. Problems with mass, momentum and energy conservation in the KdV equation were discussed in ] recently. In this paper the authors considered the KdV equations in the original dimensional variables. Then the KdV equatios are η t +cη x + 3 c h ηη x + ch 6 η xxx =, 5 in a fixed frame of reference and η t + 3 c h ηη x + ch 6 η xxx =, 6 in a moving frame. In both, c = gh, and 6 is obtained from 5 by setting x = x ct and dropping the prime sign. In our present paper we discuss the energy formulas obtained both in fixed and moving frames of reference for KdV, 4, 5, 6. There seem to be some contradictions in the literature because the form of some invariants and the energy formulas differ in different sources because of using different reference frames and/or different scalings. In this paper we address this problems. The second goal is to present some invariants for a KdV type equation of the second order. III. INVARIANTS OF KDV TYPE EQUATIONS What invariants can be attributed to equations and 5 6? It is well known, see, e.g. 7, Ch. 5], that an equation of the form T t + X =, 7 x where neither T an analog to density nor X an analog to flux contain partial derivatives with respect to t, corresponds to some conservation law. It can be applied, in particular, to KdV equations where there exist an infinite number of such conservation laws and to the equations of KdV type like. Functions T and X may depend on x,t,η,η x,η x,...,h,h x,..., but not η t. If both functions T and X x are integrable on, and lim X = const soliton solutions, then integration x ± of equation 7 yields d dt since T dx = or T dx = const., 8 X x dx = X,t X,t =. 9 The same conclusion applies for periodic solutions cnoidal waves, when in the integrals 8, 9 limits of integration, are replaced bya,b, whereb a = Λ is the space period of the cnoidal wave the wave length.

3 3 A. Invariants of the KdV equation For the KdV equation the two first invariants can be obtained easily. Writing in the form η t + η + 3 x 4 αη + 6 βη xx =. one immediately obtains the conservation of mass volume law I = ηdx = const. Similarly, multiplication of by η gives t η + x η + αη3 βη x + 6 βηη xx resulting in the invariant of the form I = =, η dx = const. 3 In the literature of the subject, see, e.g. 7, ], I is attributed to momentum conservation. Invariants I,I have the same form for all KdV equations, 4, A, 5, 6. Denote the left hand side of by KDVx,t and take 3η KDVx,t β 3α η x KDVx,t. 4 x The result, after simplifications is η 3 β t 3α η x + 9 x 8 αη4 + βη xη 5 βηx η +η3 + β 8 α η x β 9 α η xη 3x β 3α η x =. Then the next invariant for KdV in the fixed reference frame is I 3 fixed frame = η 3 β 3α η x dx = const. 6 The same invariant is obtained for the KdV in the moving frame 4. The same construction like 4 but for equation 4 results in t η 3 β 3 α η x + x 9 8 αη4 + βη xη 7 =. βηxη +η 3 + β 8 α η x β 9 α η xη 3x Then the next invariant for KdV equation in moving reference frame is I 3 moving frame = η 3 β 3α η x dx = const. 8 The procedure similar to those described above leads to the same invariants for both equations 5 and 6 where KdV equations are written in dimensional variables. To see this, it is enough to take 3η kdvx,t kdvx,t =, where kdvx,t is the lhs either of 3 h3 x 5 or 6. For equation 5 the conservation law is η 3 h3 t 3 η x + cη 3 9c x 8h η4 3 ch3 ηx 9 =, ch ηη x + ch η η xx + 8 ch5 η xx 9 ch5 η x η xxx whereas for equation 6 the flux term is different η 3 h3 t 3 η x + 9c x 8h η4 ch ηηx + ch η η xx + 8 ch5 η xx 9 ch5 η x η xxx =. But in both cases the same I 3 invariant is obtained as I 3 dimensional η = 3 h3 3 η x dx = const. Conclusion InvariantsI 3 have the same form for fixed and moving frames of reference when the transformation from fixed to moving frame scales x and t in the same way e.g. x = x t and t = t. When scaling factors are different, like in A, then the form of I 3 in the moving frame differs from the form in the fixed frame, see Appendix A. For those solutions of KdV which preserve their shapes during the motion, that is, for cnoidal solutions and single soliton solutions, integrals of any power of ηx,t and any power of arbitrary derivative of the solution with respect to x are invariants. That is, I a,n = η nx a dx = const, where n =,,,..., and a R is an arbitrary real number. Then an arbitrary linear combination of I a,n is an invariant, as well. B. Invariants of the second order equations Can we construct invariants for KdV type equations of the second order? Let us try to take T = η for equation. Then we find that all terms, except η t, can be written as X x, as η x +α 3 ηη x +β 6 η 3x +α 38 η η x +αβ 3 4 η xη x + 5 ηη 3x +β 9 ] 36 η 5x dx = η αη + 6 βη x 8 α η αβ 48 η x + 5 ηη x β η 4x.

4 4 As 3 depends on η and space derivatives and also since all those functions vanish when x ±, the conservation law for mass volume ηx,tdx = const., 4 holds for the second order equation. Conservation law 4 holds for the equation with an uneven bottom, as well. Until now we did not find any other invariants for the secod order equations. Moreover, we can show that the integral I 3 is no longer an invariant of the second order KdV equation. Upon multiplication of equation by η one obtains = t η + x 3 3 α η β η + αη3 + 6 β η x +ηη x η xx η x η 3x +ηη 4x αβη η x ]+ 8 αβηη xη x. The last term in 5 can not be expressed as x Xη,η x,... Therefore η dx is not a conserved quantity. IV. ENERGY The invariant I 3 is usually referred to as the energy invariant. Is this really the case? A. Energy in a fixed frame as calculated from the definition The hydrodynamic equations for an incompressible, inviscid fluid, in irrotational motion and under gravity in a fixed frame of reference, lead to a KdV equation of the form We will find the function η t + η x +α 3 η η x +β 6 η 3 x =. 6 f x = η 4 α η + 3 β η x x, 7 obtained as a byproduct in derivation of KdV, useful in what follows. For more details see Appendix B or 4, Chapter 5]. Tildas denote scaled dimensionless quantities. Now construct the total energy of the fluid from the definition. The total energy is the sum of potential and kinetic energy. In our two-dimensional system the energy in original dimensional coordinates is + h+η ρv E = T +V = dy dx 8 + h+η + ρgydy dx. Equations 6 and 7 are obtained after scaling,, 4]. We now have dimesionless variables, according to φ = h la gh φ, x = x l, η = η a, ỹ = y h, t = and T = ρgh l V = ρgh l +α η +α η t l/ gh, 9 ρỹdỹd x, 3 α φ x + α β φ ỹ dỹd x. 3 Note, that the factor in front of the integrals has the dimension of energy. In the following, we omit signs, having in mind that we are working in dimensionless variables. Integration in 3 with respect to y yields V = gh lρ α η +αη + dx 3 = gh lρ α η +αη ] dx+ dx. After renormalization substraction of constant term dx one obtains V = gh lρ α η +αη dx. 33 In kinetic energy we use the velocity potential expressed in the lowest first order φ x = f x βy f xxx and φ y = βyf xx, 34 where f x was defined in 7. Now the bracket in the integral 3 is α φ x + α β φ y = α fx +βy f x f xxx +fxx. 35 Inegration with respect to the vertical corrdinate y gives, up to the same order, T = ρgh l = ρgh l α f x +αη +β f x f xxx +f xx 3 +αη3 ]dx 36 α f x +αf xη + 3 β f xx f x f xxx ] dx.

5 In order to express energy through the elevation funcion η we use 7. We then substitute f x = η in terms of the third order and f x = η αη3 + 3 βηη xx in terms of the second order T = ρgh l = ρgh l α + α η αη3 + 3 βηη xx +αη β η x ηη xx ] dx η + αη3 dx 37 ] 3 β ηx +ηη xx dx. The last term vanishes as V. VARIATIONAL APPROACH + η x +ηη xx dx = ηx dx+ηη x + ηx dx =. A. Lagrangian approach, potential formulation 38 Therefore the total energy in the fixed frame is given by Some attempts at the variational approach to shallow E tot =T +V = ρgh l αη +αη + 4 water problems are summarized in G.B. Whitham s book αη3 dx 39, Sect 6.4]. =ρgh l αi +α I + 4 α I 3 + α β ηx dx The energy 39 in a fixed frame of reference has non invariant form. The last term in 39 results in small deviations from energy conservation only when η x changes in time in soliton s reference frame, what occurs only during soliton collision. This deviations are discussed and illustrated in Section VI E. The result 39 gives the energy in powers of η only. The same structure of powers in η was obtained by the authors of ], who work in dimensional KdV equations 5 and 6. On page they present a non-dimensional energy density E in a frame moving with the velocity U. Then, if U = is set, the energy density in a fixed frame is proportional to αη + α η as the formula is obtained up to second order in α. However, the third order term is 4 α3 η 3, so the formula up to the third order in α becomes E αη +α η + 4 α3 η This energy density contains the same terms like 39 and does not contain the term η x, as well. Energy calculated from the definition does not contain a proper invariant of motion. by Ali and Kalisch 8,Sect. 3] working in such system one has to replace φ x by the horizontal velocity in a moving frame, that is by φ x α = α η 4 α η +β 3 y η x x α. Then repeating the same steps as in the previous subsection yields the energy expressed by invariants E tot =ρgh l α η+4 α η + α3 η 3 ] β d x 3α η x = ρgh l αi + 4 α I + α3 I 3. 4 The crucial term 6 α β η x in 4 appears due to integration over vertical variable of the term α η x x β supplied by φ x α. For KdV as it stands, we can not write a variational principle directly. It is necessary to introduce a velocity potential. The simplest choice is to take η = ϕ x. Then equation in the fixed frame takes the form ϕ xt +ϕ xx + 3 αϕ xϕ xx + 6 βϕ xxxx =. 4 The appropriate Lagrangian density is L fixed frame := ϕ tϕ x ϕ x α 4 ϕ3 x + β ϕ xx. 4 Indeed, the Euler Lagrange equation obtained from Lagrangian 4 is just 4. For our moving reference frame the substitution η = ϕ x into 4 gives ϕ xt + 3 αϕ xϕ xx + 6 βϕ xxxx =. 43 So, the appropriate Lagrangian density is L moving frame := ϕ tϕ x α 4 ϕ3 x + β ϕ xx. 44 Again, the Euler Lagrange equation obtained from Lagrangian 44 is just B. Energy in a moving frame Now consider the total energy according to 8 calculated in a frame moving with the velocity of sound c = gh. Using the same scaling 9 to dimensionless variables we note that in these variables c =. As pointed B. Hamiltonians for KdV equations in the potential formulation The Hamiltonian for the KdV equation in a fixed frame can be obtained in the following way. Defining generalized momentum π = L ϕ t, where L is given by 4,

6 6 one obtains H = = π ϕ L]dx = η + 4 α η 3 β 3α η x ϕ x + α 4 ϕ3 x β ] ϕ xx dx ] dx. 45 The energy is expressed by invariants I,I 3 so it is a constant of motion. The same procedure for KdV in a moving frame 4 gives H = = 4 α π ϕ L]dx = η 3 β 3α η x α 4 ϕ3 x β ] ϕ xx dx dx. 46 The Hamiltonian 46 only consists I 3. The constant of motion in a moving frame is E = 4 I3 = const. 47 The potential formulation of the Lagrangian, described above, is succesful for deriving KdV equations both for fixed and moving reference frames. It fails, however, for the second order KdV equation. We proved that there exists a nonlinear expression of Lϕ t,ϕ x,ϕ xx,..., such that the resulting Euler Lagrange equation differs very little from equation. The difference lies only in the value of one of the coefficients in the second order term αβ 3 4 η xη x + 5 ηη 3x. Particular values of coefficients in this term also cause the lack of the I invariant for second order KdV equation, see 5. VI. LUKE S LAGRANGIAN AND KDV ENERGY The full set of Euler equations for the shallow water problem, as well as KdV equations, A, and second order KdV equation can be derived from Luke s Lagrangian 3], see, e.g. 3]. Luke points out, however, that his Lagrangian is not equal to the difference of kinetic and potential energy. Euler Lagrange equations obtained from L = T V do not have the proper form at the boundary. Instead, Luke s Lagrangian is the sum of kinetic and potential energy suplemented by the φ t term which is necessary in dynamical boundary condition. A. Derivation of KdV energy from the original Euler equations according to 4] In Chapter 5. of the Infeld and Rowlands book the authors present a derivation of the KdV equation from the Euler hydrodynamic equations using a single small parameter ε. Moreover, they show that the same method allows us to derive the Kadomtsev Petviashvili KP equation 9] for water waves 3, 3, 33, 34] and also nonlinear equations for ion acoustic waves in a plasma 3, 35]. The authors first derive equations of motion, then construct a Lagrangian and look for constants of motion. For the purpose of this paper and for comparison to results obtained in the next subsections it is convenient to present their results starting from Luke s Lagrangian density. That density can be written as here g = +η L = φ t + ] φx +φz +z dz. 48 In Chapter 5.. of 4] the authors introduce scaled variables in a movimg frame ε plays a role of small parameter and if ε = α = β, then KdV equation is obtained. Then for details, see Appendix B or 4, Chapter 5.] φ z = ε 3 zfξξ, φ t = εf ξ +ε f τ + z f ξξξ φ x = εf ξ ε z f ξξξ, ε 3z f ξξτ. 49 Substitution of the above formulas into the expression ] under the integral in 48 gives ] = z εf ξ +ε f τ + f ξ + z f ξξξ 5 +ε 3z fξξτ +fξξ f ξ f ξξξ ] +Oε 4. Remark The full Lagrangian is obtained by integration of the Lagrangian density 48 with respect to x. Integration limits are, for a soliton solutions, or a,b], where b a = X wave length space period for cnoidal solutions. Integration by parts and properties of the solutions at the limits, see 9, allow us to use the equivalence f ξξ f ξf ξξξ dξ = f ξξ dξ. Therefore ] = z εf ξ +ε f τ + f ξ + z f ξξξ 5 +ε 3z fξξτ +f ] ξξ +Oε 4. Integration over y gives note that +η = +εη L = +εη ++εη εf ξ +ε f τ + ] f ξ + 3 +εη3 ε f ξξξ ] ε3 f ξξτ +ε 3 fξξ. 5 Write 5 up to third order in ε L = L +εl +ε L +ε 3 L 3 +Oε 4. It is easy to show, that L =, L = η f ξ, L = f τ + η ηf ξ + f ξ + 6 f ξξξ, 53 L 3 = ηf τ + ηf ξ + ηf ξξξ 6 f ξξτ + 3 f ξξ.

7 7 The Hamiltonian density reads as L L H = f τ +f ξξτ L 54 f τ f ξξτ = +εη f ξ+ε η ηf ξ + f ξ + 6 f ξξξ +ε 3 ηf ξ + ηf ξξξ + ] 3 f ξξ. Dropping the constant term one obtains the total energy as E = εη f ξ +ε η ηf ξ + f ξ + 6 f ξξξ +ε 3 ηf ξ + ηf ξξξ + ] 3 f ξξ dξ. 55 Now, we need to express f ξ and its derivatives by η and its derivatives. We use 7 replacing α and β by ε, that is, f ξ = η 4 εη + 3 εη ξξ. 56 Then the total energy in a moving frame is expressed in terms of the second and the third invariants E = ε η dx+ε 3 η 3 ] 4 3 η ξ dx. 57 Note that the term 3 η ξ occuring in the third order invariant originates from three terms appearing in φ z, φ x and φ t see terms f ξξ and f ξξξ in 49. B. Luke s Lagrangian The original Lagrangian density in Luke s paper 3] is L = hx After scaling as in,, 4] ρ φ t + ] φ x +φ y+gy dy. 58 φ = h la gh φ, x = x l, η = η a, ỹ = y h, t = we obtain t l/ gh, 59 φ t = ghα φ t, φ x = ghα φ x, φ y = ghα β φ ỹ. 6 The Lagrangian density in scaled variables becomes dy = hdỹ +αη L = ρgha φ t + ] φ x + α β φ ỹ dỹ + ρgh +αη. 6 So, in dimensionless quatities L +αη ρgha = φ t + α φ x + α β φ ỹ ] dỹ + αη, 6 where the constant term and the term proportional to η in the expansion of +αη are omitted. The form 6 is identical with Eq..9 in Marchant & Smyth 3]. The full Lagrangian is obtained by integration over x. In dimensionless variables dx = ld x it gives +αη L=E φ t + α φ x + α β φ ỹ ]dỹ+ ] αη d x. 63 The factor in front of the integral, E = ρghal = ρgh lα, has the dimension of energy. Next, the signs will be omitted, but we have to remember that we are working in scaled dimensionless variables in a fixed reference frame. C. Energy in the fixed reference frame Express the Lagrangian density by η and f = φ. Now, up to first order in small parameters φ = f βy f xx, φ t = f t βy f xxt, φ x = f x βy f xxx, φ y = βyf xx. 64 Then the expression under the integral in 6 becomes ] = f t βy f xxt + αf x + αβy f x f xxx +f xx. 65 From properties of solutions at the limits fx f xxx +fxx f xx. Integration of 65 over y yields L ρgha = f t + αf x +αη βf xxt +αβf xx 3 +αη3 3 +αη3 + αη. 66 The dimensionless Hamiltonian density is L L f t f t +f xxt f xxt L H ρgh l = α αf x +αη+αβf xx 3 +αη3 + ] αη. 67 Again, we need to express the Hamiltonian by η and its derivatives, only. Inserting f x = η 4 αη + 3 βη xx 68 into 67 and leaving terms up to third order one obtains H ρgh l = α αη + 4 α η 3 + ] 3 αβη x +ηη xx. 69

8 8 The energy is αη + 4 α η 3 + ] 3 αβη x +ηη xx dx E ρgh l = α = α η dx+ 4 α3 ] η 3 dx 7 since the integral of the αβ term vanishes. Here, in the same way as in calculations of energy directly from the definition 39, the energy is expressed by integrals of η and η 3. The term proportional to αη is not present in 7, because it was dropped earlier 3]. D. Energy in a moving frame Transforming into the moving frame x = x t, t = αt, x = x, t = x +α t. 7 φ = f βy f x x, φ x = f x βy f x x x, φ y = βyf x x, 7 φ t = f x + βy f x x x +αf t βy f x x t. 73 Up to second order αφ x + α β φ y = αf x +αβy f x f x x x +f x x ] = αf x +αβy f x x. 74 Therefore the expression under the integral in 6 is ] = f x + βy f x x x +αf t βy f x x t+ αf x+αβy f x x. 75 Integration yields L ρgha = f x +αf t + αf x +αη αη3 βf x x x f x x t+αβf x x + αη. Like in 67 above, the Hamiltonian density is H ρgh l = α f x + αf x +αη αη3 βf x x x +αβf + x x ] αη. Expressing f x by 68 one obtains H ρgh l = α 4 αη + 3 βη xx α η αβ 4 η x 5 ηη xx ] 8 β η xxxx. Finally the energy is given by E ρgh l = α 4 η dx+α 3 η 3 β 3α η x dx 79 since integrals from terms with β,β vanish at integration limits, and 5 ηη xx 5 η x. The invariant term proportional to αη is not present in 79, because it was dropped in 6. If we include that term, the total energy is a linear combination of all three lowest invariants, I,I 3,I 3. Comment An almost identical formula for the energy in a moving frame, for KdV expressed in dimensional variables 6, was obtained in ]. That energy is expressed by all three lowest order invariants E = c + ηdx+ 4 c h c h η 3 h3 3 η x dx, η dx 8 as well. Translation of 8 to nondimensional variables yields E = gh l αi + 4 α I + α3 I 3. E. How strongly is energy conservation violated? E.8935 E FIG.. Precision of energy conservation for 3-soliton solution. Energies are plotted as open circles E and open squares E for 4 time instants. The total energy in the fixed frame is given by equation 39. Taking into account its non-dimensional part we may write E t = T +V gh l = = αi +α I + 4 t αη +αη + 4 αη3 ] dx αη 3 dx 8

9 9 In order to see how much the changes ofe violate energy conservation we will compare it to the same formula but expressed by invariants E t = αi +α I + 4 α3 I 3. 8 The time dependence of E and E is presented in Fig. for a 3-soliton solution of KdV. Presented is time evolution in the interval t,]. The shape of the 3-soliton solution is presented only for three times t =, 6, in order to show shapes changing during the collision. For presentation the example of a 3-soliton solution with amplitudes equal,5, and.5 was chosen. In Fig. 3 the positions of solutions at given times were artificially shifted to set them closer to each other. The plots in Figs. and 3 for t > are symmetric to those which are shown in the figures. For this example the relative discrepancy of the enregy E from the constant value, is very small δe = E t = E t = E t = However, the E energy is conserved with numerical precision of thirteen decimal digits in this example. In a similar example with a -soliton solution with apmlitudes and.5 the relative error 83 was even smaller, with the value δe.4. This suggests that the degree of nonconservation of energy increases with n, where n is the number of solitons in the solution. ηx,t x t=- t = -6 t = FIG. 3. Shape evolutiom of 3-soliton solution during collision. We confirmed some known facts. Firstly, that the usual form of the energy H = T +V is not always expressed by invariants only. The reason lies in the fact, as pointed out by Luke in 3], that the Euler Lagrange equations obtained from the Lagrangian L = T V do not supply the right boundary conditions. Secondly, the variational approach based on Luke s Lagrangian density provides the right Euler equations at the boundary and allows for a derivation higher order KdV equations. In the frame moving with the velocity of sound all energy components are expressed by invariants. Energy is conserved. Numerical calculations confirm that invariants I,I,I 3 in the forms, 3, 6, 8 are exact constants of motion for two- and threesoliton solutions, both for fixed and moving coordinate systems. In all performed tests the invariants were exact up to fourteen digits in double precision calculations. For the extended KdV equation we have only found one invariant of motion I 4. The total energy in the fixed coordinate system as calculated in 39 is not exactly conserved but only altered during collisions, even then by minute quantities an order of magnitude smaller than expected. Details in figure caption of figure. A summary of these conclusions can be found in Table I. VIII. EXTENDED KDV EQUATION In this section we calculate energy formula coresponding to a wave motion governed by second order equations in scaled variables, that is the equation for the fixed coordinate system and the correponding equation for a moving coordinate system. As previosly we compare energies calculated from the definition with those Luke s Lagrangian. A. Energy in a fixed frame calculated from definition VII. CONCLUSIONS FOR KDV EQUATION The main conclusions can be formulated as follows The invariants of KdV in fixed and moving frames have the same form. Of course when we have the same scaling factor for x and t in the transformation between frames. Now, instead of we consider the second order KdV equation, that is called by Marchant & Smyth 3] "extended KdV". In section IV A, total energy of the wave governed by KdV equation, that is the equation with terms only up to first order in small parameters was obtained in 39. In calculation according to eq. the potential energy is expressed by the same formula 3 as previously for KdV equation. In the expression for kinetic energy the

10 TABLE I. Comparison of different energy formulas. Here η 3 = η 3 dx. Formulas in this column are written in Euler Luke s Integrals Potential KdV equations Lagrangian T +V Lagrangian dimensional E gh l. Fixed frame αi +α I + 4 α3 η 3 α I + 4 α3 η 3 αi +α I + 4 α3 η 3 I + 4 αi3 αi +α I + 4 α3 η Moving frame 4 α I + α3 I 3 4 α I + α3 I 3 αi + 4 α I + α3 I 3 4 αi3 αi + 4 αi + α I velocity potential has to be expanded to second order in small parameters φ = f βy f xx + 4 β y 4 f xxxx, 84 with derivatives φ x = f x βy f xxx + 4 β y 4 f xxxxx, φ y = βyf xx + 6 β y 3 f xxxx. 85 Integrating over y and retaining terms up to fourth order yields T = ρgh l α fx +αηfx + 3 β fxx f x f xxx +αβηfxx ηf xf xxx 86 +β f xxx 5 f xxf xxxx + ] 6 f xf xxxxx dx. Expression 86 limited to first line gives kinetic energy for KdV equation, see 36. Now, we use the expression for f x and its derivatives up to second order, see e.g. 3, Eq..7],, Eq. 7] f x = η 4 αη + 3 βη xx + 8 α η αβ 6 η x + ηη xx + β η xxxx. Insertion 87 and its derivatives into 86 gives T = ρgh l α η + αη3 + 3 β ηx +ηη xx α η 4 +αβ 4 ηη x η η xx 88 +β η xx η xη xxx ηη xxxx ] dx. From properties of solutions at x ± terms with β and β in square bracket vanish and the term with αβ can be written form. Finally one obtains T = ρgh l α η + αη3 3 6 α η 4 7 ] 4 αβηη x dx. 89 Then total energy is the sum of 33 and 89 E tot = ρgh l αη +αη + 4 αη αη4 7 ] 48 α3 βηηx dx. The first three terms are identical as in KdV energy formula 39, the last two terms are new for extended KdV equation. B. Energy in a fixed frame calculated from Luke s Lagrangian Calculate energy in the same way as in Section VI, C, but in one order higher. In scaled coordinates Lagrangian density is expressed by 6 here we keep infinite constant term { +αη L = ρgh l α φ t + αφ x + α ] β φ y dy + } +αη. 9 From 84 we have φ t = f t βy f xxt + 4 β y 4 f xxxxt. 9 Inserting 9 and 85 into 9, integrating over y and retaining terms up to third order one obtains constant term id dropped { L ρgh l = α η +f t +α η +ηf t + f x βf xxt + α ηfx +αβ 6 f xx ηf xxt 6 f xf xxx + β f xxxxt + α β ηf xx η f xxt ηf x f xxx +αβ 4 f xxx 3 f xxf xxxx ηf xxxxt + f xf xxxxx β 3f xxxxxxt 54 }.

11 The the Hamiltonian density H = f t L f t +f xxt L f xxt +f 4xt L f 4xt +f 6xt L f 6xt L is H ρgh l = αη α η +fx α3 ηfx +α β 6 fxx + 6 f xf xxx +α 3 β ηfxx + ηf xf xxx +α β 4 f xxx + 3 f xxf xxxx f xf xxxxx Now, we use f x in the second order 87 and its derivatives. Insertion these expressions into 94 nd retention terms up to thired order yields H ρgh l = αη α η 4 α3 η α4 η 4 +α β 6 η x 6 ηη xx +α 3 β 9 48 ηη x 3 8 η η xx +α β 4 η xx 7 9 η xη xxx 9 36 ηη xxxx 95. The energy is obtained by integration of 95 over x using integration by parts and properties ofη and its derivatives at x ±. Then terms with αβ and αβ vanish. The final result is E = ρgh l αη +αη + 4 αη αη4 7 ] 48 α3 βηηx dx, the same as 9 but with the opposite sign. C. Energy in a moving frame from definition Let us follow arguments given by Ali and Kalisch, Sec. 3] and used already in Section IV B. Working in a moving frame one has to replace φ x by the horizontal velocity in a moving frame, that is, φ x α. Then in a frame moving with the sound velocity we have φ x = f x βy f xxx + 4 β y 4 f xxxxx α, φ y = βyf xx + 6 β y 3 f xxxx Then the expression under integral over y in 3 becomes in the following terms up to fourth order are kept α φ x + α β φ y = αf x +α fx +y α βfxx 98 +y αβf xxx y α βf x f xxx y4 αβ f xxxxx +y 4 α β 4 f xxx 3 f xxf xxxx + f xf xxxxx. After integration over y one obtains T = ρgh l +αη fx +α ηf x +f x + 3 αβf xxx +α 3 ηfx 6 αβ f xxxxx +α β 3 f xx +ηf xxx 3 f xf xxx +α 3 β ηfxx +η f xxx ηf x f xxx +α β f xxx 5 f xxf xxxx 99 ηf xxxxx + 6 f xf xxxxx ] dx. Then insertion f x 87 and its derivatives yields T = ρgh l αη α η 3 αβη xx α3 η 3 α β 4 η x + ηη xx 9 8 αβ η xxxx α4 η 4 +α 3 β ηη x η η xx + α β 3 η xx η xη xxx ηη xxxx + ] 36 αβ3 η xxxxxx dx, where constant term is dropped. Using properties of solutions at x ± this expression can be simplified to T = ρgh l αη α η α3 η α4 η α βηx + α3 βηηx + ] α β ηxx dx. Then total energy is E tot = ρgh l αη + 4 α η α3 η α4 η α βηx + 4 α3 βηηx + ] 4 α β ηxx dx. In special case α = β this formula simplifies to E tot = ρgh l αη α η +α 3 8 η η x +α η4 + 4 ηη x + ] 4 η xx dx. 3

12 D. Energy in a moving frame from Luke s Lagrangian Follow considerations in Section VI, but with KdV equation. Transforming into the moving frame through 7 we have now φ = f βy f x x + 4 β y 4 f x x x x, 4 φ x = f x βy f x x x + 4 β y 4 f x x x x x, 5 φ y = βyf x x + 6 β y 3 f x x x x, 6 φ t = f x + βy f x x x 4 β y 4 f x x x x x 7 +αf t βy f x x t + 4 β y 4 f x x x x t. Inserting 4 7 into 9 one obtains Lagrangian density in moving frame as constant term is dropped as previously L ρgh = αη f x+α l η +f t ηf x + f x + 6 αβf x x x +α ηf t 3 + ηf x αβ f x x x x x + α β 6 f x x 6 f x x t + ηf x x x 6 f xf x x x + α 3 β ηf x x ηf x x t + η f x x x ηf xf x x x + α β 4 f x x x 3 f x xf x x x x + 4 ηf x x x x x + f xf x x x x x Then Hamiltonian density L H = f t f t L L +f x x t +f x x x x t f x x t f x x x x t after insertion of 8 into 9 yields H ρgh = α η +f x+α l f x x x x t. 8 L 9 η +ηf x f x 6 αβf x x x α3 ηf x + αβ f x x x x x + α β 6 f x x ηf x x x + 6 f xf x x x + α 3 β ηf x x η f x x x + ηf xf x x x + α β 4 f x x x + 3 f x xf x x x x + 4 ηf x x x x x f xf x x x x x. In order to express by η only we use f x in the form 87 and its derivatives. It gives H ρgh l = 4 α η + 6 αβη x x 3 8 α3 η α β 48 η x + 4 ηη xx αβ η x x x x α4 η 4 α 3 β 4 ηη x η η x x α β 6 η x x η xη x x x ηη x x x x 7 αβ3 η x x x x x x Then energy is given by the integral E = gh l 4 α η + 6 αβη x x 3 8 α3 η 3 5 +α β 48 η x + 4 ηη xx αβ η x x x x α4 η 4 7 α 3 β 4 ηη x η η x x 7 αβ3 η x x x x x x α β 6 η x x η xη x x x + 9 ] 7 ηη x x x x dx. From properties of solution integrals of terms with αβ,αβ,αβ 3 vanish and terms with α β,α 3 β,α β can be simplified. Finally, energy is given by the follwing expression E = gh l 4 α η 3 8 α3 η α4 η α βη x 4 α3 βηη x 4 α β η x x ] dx. In special case when β = α the result is E = gh l 3 4 α η α 3 8 η η x 7 +α 4 3 η4 4 ηη x ] 4 η x x dx. 4 If the invariant term I αηdx is dropped in 9 or 96 then the energy calculated in the moving frame 3 have the same value but with oposite sign. E. Numerical tests. Fixed frame In order to check energy conservation for the extended KdV equation we performed several numerical tests. First, discuss energy conservation in a fixes frame. We calculated time evolution governed by the equation of waves which initial shape was given by -, - and 3- soliton solutions of the KdV first order equations. For

13 3 ηx,t t= t= t=4 t=35 range considered they do not extend.%,.4% and.5% for -, -, 3-soliton waves, respectively they increase with more complicated waves. For additional check of numerics the invariant I = αηx,tdx for the eaquation was plotted as Mass. In spite of approximate integration the value of I was obtained constant up to digits for all initial conditions FIG. 4. Example of time evolution of 3-soliton solution. presentation the following initial conditions were chosen. 3-soliton solution have amplitudes.5, and.5, - soliton solution have amplitudes and.5 and -soliton solution the amplitude. The changes of energy presented in Figs. 6 and 5 are qualitatively the same also for different amplitudes. An example of such time evolution for 3-soliton solution is presented in Fig. 4. Time range in Fig. 4 contains initial shape of 3-soliton solution with almost separated solitons at t =, intermediate shapes and almost ideal overlap of solitons at t = 35. In order to do not obscure details the subsequent shapes are shifted verticaly with respect to the previous ones. Note additional slower waves after the main one which are generated by second order terms of the equation, that were already discussed in ]. E/E En -sol En -sol En 3-sol Mass FIG. 5. Energy nonconservation for the extended KdV equation in the fixed frame. Symbols represent values of the total energy given by formulas 9 or 96. Full square symbols represent the invariant I. We see that the total energy for waves which move according to the extended KdV equation is not conserved. Although energy variations are generally small in time x t F. Moving frame Here we present variations of the energy calculated in a moving frame. The time evolution of the wave is given by the equation transformed with 7, that is the equation η t+ 3 ηη x + β 6α η 3 x αη η x +β 4 η xη x + 5 ηη 3 x + 9 β 36 α η 5 x =. E/E En-sol En-sol En-3sol Mass FIG. 6. Energy nonconservation for the extended KdV equation in the moving frame 5. Symbols represent values of the total energy given by the formula. Full square symbols represent the invariant I. The time range of the evolution was chosen for a convenient comparison with the numerical results obtained in fixed reference frame, that is - and 3-soliton waves move from separate solitons to fully colliding time instant. The convention of symbols is the same as in Fig. 5. the energy is calculated according to the formula. In moving coordinate system energy variations are even greater than in the fixed reference frame, because in the time period considered it approaches values of.%,.% and.% for -, - and 3-soliton waves, respectively. This increase of relative time variations of energy can not be atributed only to two times smaller leading term αη in with respect to 9. Again, in spite of approximate integration the value of I was obtained constant up to digits for all initial conditions. t

14 4. Conclusions for extended KdV equation We calculated energy of the fluid governed by the extended KdV equation in two cases:. In a fixed frame sections VIII A and VIII B.. In the frame moving with the sound velocity sections VIII C and VIII D. In both cases we calculated energy using two methods: from definition and from Luke s Lagrangian. Both methods give consistent results. For fixed frame energies 9 and 96 are the same. For moving frame the energy calculated from the definition contains one term more then energy calculated from Luke s Lagrangian, but this term αηdx is the invariant I. When this term is dropped both energies in moving coordinate system and 3 are the same and energies in both coordinate systems differ only by sign. The general conclusion concerning energy conservation for shallow water wave problem can be formulated as follows. Since there exists the Lagrangian of the system Luke s Lagrangian then exact solutions of Euler equations have to conserve energy. However, when approximate equations of different orders resulting from exact Euler equations are considered, energy conservation is not a priori determined. The KdV equations obtained in first order approximation has a miraculous property, an infinite number of invariants with energy among them. However, this astonishing property is lost in second order approximation to Euler equations and energy in this order may be conserved only approximately. Appendix A The simplest, mathematical form of the KdV equation is obtained from by passing to the moving frame with additional scaling x = 3 x t, t = 4 3 αt, u = η, A which gives a standard, mathematical form of the KdV equation u t +6uu x + β α u x x x =, or u t +6uu x +u x x x = for β = α. A Equations A, particularly with β = α are favored by mathematicians, see, e.g. 5]. This form of KdV is the most convenient for ISM the Inverse Scattering Method, see, e.g. 6 8]. For the moving reference frame, in which the KdV equation has a standard mathematical form A, the invariant I 3 is slightly different. To see this difference denote the lhs of A by KDVmx,t and construct 3η KDVmx,t β α η x KDVmx,t =. x Then after simplifications one obtains η 3 β t α η x + 9 x η4 6 β α ηη x A3 +3 β α η η xx β β α η xx + η x η xxx] =, α which implies the invariant I 3 in the following form I 3 moving frame = η 3 β α η x dx = const or A4 I 3 moving frame = η 3 η x dx = const for β = α. We see, however, that the difference between A4 and 8 is caused by additional scaling. In the Lagrangian approach as described in Sect. V, the substitution u = ϕ x into A gives ϕ xt +6ϕ x ϕ xx +ϕ xxxx =. A5 Then the appropriate Lagrangian density for equation A with α = β is L standard KdV := ϕ tϕ x ϕ 3 x + ϕ xx. A6 Indeed, the Euler Lagrange equation obtained from the Lagrangian A6 is just A5. The Hamiltonian for KdV A can be found e.g. in ]. Defining generalized momentum π = L, where ϕ t L is given by A6, one obtains ] L H = π ϕ L]dx = dx A7 = ϕ 3 x ] ϕ xx dx = ϕ t ϕ t L η 3 η x ] dx. This is the same invariant as I 3 moving frame in A4. Appendix B The set of Euler equations for irrotational motion of an incompresible and inviscid fluid can be written neglecting surface tension in dimensionless form: φ = φ z = on z = B B η t +φ x η x φ x = on z = +η B3 φ t + φ x +φ z +η = on z = +η. B4 We look for solutions to the Laplace equation B in the form φ = z n f n x,y,t B5 n=

15 yielding n= nn z n f n +z n f n] =. B6 They then applied a transformation to a frame moving with velocity of sound. The coordinates scales as 5 ξ = ε x t, τ = ε 3 t B In two dimensions x,z we obtain f n+ = f n n+n+ x. B7 The boundary condition at the bottom, φ z = at z= implies f = and then all odd f k+ =. Now φ = n= m zm m! m f x m, B8 where f := f. In the stretched coordinates x = ε ξ so φ = ε f + m zm ε m m! ξ f. B9 n= Now both B and B are satisfied. We must also satisfy the boundary conditions on z = +η. In the derivation of KdV and Kadomtsev-Petiashvili 9] from the Euler equations B B4 Infeld and Rowlands 4] applied scaling assuming the following relations vawelength : depth : amplitude as ε / : : ε. t = ε ξ +ε 3 τ, x = ε ξ. B For the wave amplitude and velocity potential the appropriate scaling was and η = εη +ε η +..., φ = ε φ +ε 3 φ +... Then the lowest order expression for φ is B B3 φ ε 3 z f ε f ξξ B4 Next, Infeld and Rowlands show that in order to simultaneously satisfy B3 and B4 the next order contributions to η and φ cancel. It is enough to keep +η = +εη and φ = ε φ and drop upper index in what follows. B5 ] A. Karczewska, P. Rozmej and Ł. Rutkowski, A new nonlinear equation in the shallow water wave problem, Physica Scripta 89, ] A. Karczewska, P. Rozmej and E. Infeld, Shallow-water soliton dynamics beyond the Korteweg de Vries equation, Phys. Rev. E 9, ] T.R. Marchant and N.F. Smyth, The extended Korteweg de Vries equation and the resonant flow of a fluid over topography, J. Fluid Mech., ] G.I. Burde, A. Sergyeyev, Ordering of two small parameters in the shallow water wave problem, J. Phys. A: Math. Theor. 46, ] R.M. Miura, KdV equation and generalizations I. A remarkable explicit nonlinear transformation, J. Math. Phys. 9, ] R.M. Miura, C.S. Gardner and M.D. Kruskal, KdV equation and generalizations II. Existence of conservation laws and constants of motion, J. Math. Phys. 9, ] P.G. Drazin and R.S. Johnson, Solitons: An Introduction, Cambridge University Press, Cambridge, ] A.C. Newell, Solitons in Mathematics and Physics, Philadelphia: Society for Industraial and Applied Mathematics, ] J.G.B. Byatt-Smith, On the change of amplitude of interacting solitary waves, J. Fluid Mech. 8, ] S. Kichenassamy and P. Olver, Existence and nonexistence of solitary wave solutions to higher-order model evolution equations, SIAM J. Math. Anal., 3, ] T.R. Marchant and N.F. Smyth, Soliton Interaction for the Korteweg-de Vries equation, IMA J. Appl. Math. 56, ] T.R. Marchant, Coupled Korteweg - de Vries equations describing, to higher-order, resonant flow of a fluid over topography, Phys. Fluids, No. 7, ] T.R. Marchant, High-order interaction of solitary waves on shallow water, Studies in Appl. Math. 9, -7. 4] T.R. Marchant, Asymptotic solitons for a higher-order modified Korteveg-de Vries Equations, Phys. Rev. E 66, ] Q. Zou and CH-H. Su, Overtaking collision between two solitary waves, Phys. Fluids 9, No. 7, ] E. Tzirtzilakis, V. Marinakis, C. Apokis and T. Bountis, Soliton-like solutions of higher order wave equations of Korteweg-de Vries type, J. Math. Phys. 43, No., ] G.I. Burde, Solitary wave solutions of the higher-order KdV models for bi-directional water waves, Commun. Nonlinear Sci. Numerical Simulat. 6, ] R. Hirota, Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons, Phys. Rev. Lett.

16 6 7, ] R. Hirota, The Direct Method in Soliton Theory, Cambridge University Press, Cambridge, 4, first published in Japanese 99. ] A. Ali and H. Kalisch, On the formulation of mass, momentum and energy conservation in the KdV equation, Acta Appl. Math. 33, ] G.B. Whitham, Linear and Nonlinear Waves, Wiley, New York, 974. ] F. Cooper, C. Lucheroni, H. Shepard and P. Sodano, Variational Method for Studying Solitons in the Korteweg-deVries Equation, Phys. Lett. A ] J.C. Luke, A variational principle for a fluid with a free surface, J. Fluid Mech. 967, 7, part, ] E. Infeld and G. Rowlands, Nonlinear Waves, Solitons and Chaos, nd edition, Cambridge University Press, Cambridge,. 5] P. Lax, Integrals of nonlinear equations of evolution and solitary waves, Comm. Pure Applied Math. 5, ] C.S. Gardner, J.M. Greene, M.D. Kruskal, and R.M. Miura, Method for Solving the Korteweg-deVries Equation, Phys. Rev. Lett. 9, ] M. Ablowitz, H. Segur, Solitons and the Inverse Scattering Transform, SIAM, Philadelphia, 98. 8] M. Ablowitz, P. Clarkson, Solitons, Nonlinear Evolution Equations and Inverse Scattering, Cambridge University Press, Cambridge, 99. 9] B.B. Kadomtsev and V.I. Petviashvili, On the stability of solitary waves in weakly dispersive media, Dokl. Akad. Nauk SSSR 9, ; Sov. Phys. Dok. 5, ] E. Infeld, G. Rowlands and M. Hen, Three dimenional stability of KdV waves and solitons, Acta Phys. Polon. A 54, 3-43, ] E. Infeld and G. Rowlands, Three-dimensional stability of Korteweg de Vries waves and solitons, II, Acta Phys. Polon. A 56, 39-33, ] E. Infeld, Three-dimensional stability of Korteweg de Vries waves and solitons. III. Lagrangian methods, KdV with positive dispersion, Acta Phys. Polon. A , ] T.B. Benjamin, The stability of solitary waves, Proc. R. Soc. Lond. A 38, ] M.S. Longuet-Higgins and J.D. Fenton, On the mass, momentum, energy and circulation of a solitary wave. II, Proc. R. Soc. Lond. A 34, , ] E. Infeld and G. Rowlands, Stability of nonlinear ion sound waves and solitons in plasma, Proc. R. Soc. Lond. A 366,