arxiv: v2 [nlin.ps] 5 Apr 2018

Transcription

1 Nonlinear Dynamics manuscript No. will be inserted by the editor) Superposition solutions to the extended KdV equation for water surface waves Piotr Rozmej Anna Karczewska Eryk Infeld arxiv: v [nlin.ps] 5 Apr 08 Received: July 6, 08/ Accepted: Abstract The KdV equation can be derived in the shallow water limit of the Euler equations. Over the last few decades, this equation has been extended to include higher order effects. Although this equation has only one conservation law, exact periodic and solitonic solutions exist. Khare and Saxena [5, 6, 7] demonstrated the possibility of generating new exact solutions by combining known ones for several fundamental equations e.g., Korteweg - de Vries, Nonlinear Schrödinger). Here we find that this construction can be repeated for higher order, non-integrable extensions of these equations. Contrary to many statements in the literature, there seems to be no correlation between integrability and the number of nonlinear one variable wave solutions. Keywords Shallow water waves extended KdV equation analytic solutions nonlinear equations PACS 0.30.Jr a Bb Introduction A long time ago, Stokes opened the field of nonlinear hydrodynamics by showing that waves described by nonlinear P. Rozmej Faculty of Physics and Astronomy, University of Zielona Góra, Szafrana 4a, Zielona Góra, Poland Tel.: Fax: P.Rozmej@if.uz.zgora.pl A. Karczewska Faculty of Mathematics, Computer Science and Econometrics, University of Zielona Góra, Szafrana 4a, Zielona Góra, Poland A.Karczewska@wmie.uz.zgora.pl E. Infeld National Centre for Nuclear Research, Hoża 69, Warszawa, Poland Eryk.Infeld@ncbj.gov.pl models can be periodic []. Although several related results followed, it took half a century before the Korteweg - de Vries equation became widely known []. A more accurate equation system, Boussinesq, was formulated in 87. It is also the theme of several recent papers [3, 4]. Another direction research has gone in is including perpendicular dynamics in KdV, e.g., [5]. The KdV equation is one of the most succesful physical equations. It consists of the mathematically simplest possible terms representing the interplay of nonlinearity and dispersion. This simplicity may be one of the reasons for success. Here we investigate this equation, improved as derived from the Euler inviscid and irrotational water equations. Just as for conventional KdV, two small parameters are assumed: wave amplitude/deptha/h) and depth/wavelength squared H/l). These dimensionless expansion constants are called α and β. We take the expansion one order higher. The new terms will then be of second order. This procedure limits considerations to waves for which the two parameters are comparable. Unfortunately some authors tend to be careless about this limitation. The next approximation to Euler s equations for long waves over a shallow riverbed is η is the elevation above a flat surface divided by H) η t + η x + 3 α ηη x+ 6 β η 3x 3 8 α η η x ) 3 + αβ 4 η xη x + 5 ) ηη 3x β η 5x = 0. In ) and subsequently we use low indexes for derivatives η nx n η x ). n This second order equation was called by Marchant and Smyth [6,7] the extended KdV. It was also derived in a different way in [8] and [9,0]. We call it KdV. It is not integrable. However, by keeping the same terms but changing one numerical coefficient specifically, replacing 3 4 by 5 6 ) we can obtain an integrable equation [,].

2 Piotr Rozmej et al. Not only is KdV non-integrable, it only seems to have one conservation law volume or mass) [3]. However a simple derivation of adiabatically conserved quantities can be found in [4]. Recently, Khare and Saxena [5, 6, 7] demonstrated that for several nonlinear equations which admit solutions in terms of elliptic functions cnx, m), dnx, m) there exist solutions in terms of superpositions cnx,m)± mdnx,m). They also showed that KdV which admits solutions in terms of dn x,m) also admits solutions in terms of superpositions dn x,m)± mcnx,m)dnx,m). Since then we found analytic solutions to KdV in terms of cn x,m) [8,9] the results of Khare and Saxena [5,6,7] inspired us to look for solutions to KdV in similar form. Exact periodic solutions for KdV First, we repeat shortly the results obtained by Khare and Saxena [5], but formulating them for KdV in a fixed frame, that is, for the equation η t + η x + 3 α ηη x+ 6 β η 3x = 0. ) Assuming solution in the form ηx,t)=adn [Bx vt),m] 3) one finds A= 4 3 B β α and v=+ 3 β B m)=+ α A m). Next, the authors [5] showed that superpositions η ± x,t)= A dn [Bx vt),m] 4) ± mcn[bx vt),m]dn[bx vt),m] ) are solutions to ) with the same relation between A and B, but for a different velocity, v ± = + 6 β B 5 m). Now, we look for periodic nonlinear wave solutions of KdV ). Introduce y := x vt. Then ηx, t) = ηy), η t = vη y and equation ) takes the form of an ODE v)η y + 3 α ηη y+ 6 β η 3y 3 8 α η η y 5) 3 + αβ 4 η yη y + 5 ) ηη 3y β η 5y = 0.. Single periodic function dn First, we recall some properties of the Jacobi elliptic functions arguments are omitted) sn + cn =, dn + msn =. 6) Their derivatives are d sn dy = cndn, d cn dy = sndn, d dn = msncn. 7) dy Assume a solution of ) in the same form as KdV solution 3). Insertion of 3) into 5) yields ABm 80 cndnsn F 0 + F cn + F 4 cn 4) = 0. 8) Equation 8) holds for arbitrary arguments when F 0,F,F 4 vanish simultaneously. The explicit form of this set of equations is following F 0 = 35α A m ) +30αAm ) β B 63m 0)+8 ) 8 9β B 4 7m 7m+ ) + 30β B m )+45 ) + 360v= 0, 9) F = 30m [ 9α A m )+6αA β B 3m )+3 ) 8β B 9β B m )+3 )] = 0, 0) F 4 = 45m 3α A + 86αAβ B 5β B 4) = 0. ) Equation ) is equivalent to the [0, Eq. 6)] obtained for solitonic solutions to KdV. Denoting z := B β Aα one obtains from ) two possible solutions z = < 0 and z = > 0. ) The case z=z leads to B < 0 and has to be rejected as in previous papers [0,8]. Then for z=z the amplitude A is A= B β α Inserting this into 0) yields 305 ) 3 4 B = 703β m) > 0. 3) 4) and then using 3) one has 305 ) 3 5 A=. 5) 37αm ) Finally from 9) one obtains { v= ) 305 m ) } m ). / { 500m ) }. 6) Despite the same form of solutions to KdV and KdV there is a fundamental difference. KdV only imposes two conditions on coefficients A,B,v for given m, so there is one

3 Superposition solutions 3 parameter freedom. This is no longer the case for KdV..0 Comparison to KdV solutions Is a solution of KdV much different from the KdV solution for the same m? In order to compare solutions of both equations, remember that the set of three equations 9)-) fixes all A,B,v coefficients for KdV for given m. In the case of KdV the equation analogous to 8) only imposes two conditions on three parameters. Therefore one parameter, say amplitude A, can be chosen arbitrary. Then we compare coefficients of solutions to KdV and KdV chosing the same value of A, that is, A KdV. Such comparison is displayed in Fig. for α = β = As Bs Vs B V Fig. The amplitude Aα, m) given by 5) as function of α, m).. Superposition "dn + mcndn" Now assume the periodic solution to be in the same form as the corresponding superposition solution of KdV [5] functio Fig. Coefficients A,B,v, z=z ) as functions of m for the periodic solution of KdV in the form of a single function 3). Lines represent KdV coefficients, points KdV coefficients. A KdV = A KdV ). η + y)= A[ dn By,m)+ mcnby,m)dnby,m) ], 7) where A,B,v are yet unknown constants m is the elliptic parameter). We will need It is clear that v KdV and v KdV are very similar. We have the following relations: for KdV B A = 3α 4β, whereas for KdV B A = α β z. Since z, B KdV /B KdV = 3 4z.. The same relations hold between KdV and KdV coefficients for superposition solutions shown in Fig. 3. The above examples for the case α = β = 0 show that for somewhat small values of α the coefficients of KdV dn solutions are not much different from those of KdV. However, physically relevant exact solutions of dn to KdV can be found for much larger values of the parameter α. In figure the amplitude Aα,m) given by 5) is shown as a contour plot for the region α [0.0,], m [0,]. It is clear that reasonable amplitudes occur in wide regions of α and m values. The amplitudes become too big only for α 0, but in such cases KdV works very well. η y = AB m mcn+dn ) sn, 8) η y = AB m mcn+dn ) cndn+ msn ), 9) η 3y = AB3 m mcn+dn ) sn mcn + 6 mcndn+dn 4msn ), 0) η 5y = AB5 m mcn+dn ) ) [ sn m cn m 3/ cn 3 dn+ dn 4 44dn sn Denote 5) as + 6m sn 4 30 mcndn dn + 4msn ) +cn 74mdn 44m sn )]. E + E + E 3 + E 4 + E 5 + E 6 + E 7 = 0, )

4 4 Piotr Rozmej et al. where E = v)η y = AB v) m mcn+dn ) sn, 3) E = 3 α ηη y = 3 8 α A B m mcn+dn ) 3 sndn, 4) E 3 = 6 β η3y= β AB3 m mcn+dn ) sn 5) mcn + 6 mcndn+dn 4msn ), E 4 = 3 8 α η η y = 3 64 α A 3 B mdn mcn+dn ) 4 sn, 6) E 5 = 3 4 αβ η yη y = 3 96 αβ A B 3 m mcn+ dn ) 4 sn cndn+ msn ), 7) E 6 = 5 αβ ηη 3y = 5 48 αβ A B 3 mdn mcn+dn ) 3 sn mcn + 6 mcndn+dn 4msn ), 8) E 7 = β η 5y = 9 70 β AB 5 m mcn+dn ) sn [ m cn m 3/ cn 3 dn+dn 4 44dn sn + 6m sn 4 30 mcndn dn + 4msn ) 9) +cn 74mdn 44m sn )]. Then ) becomes AB m mcn+ dn ) sn 30) F 0 + F cd cndn+ F c cn + F c 3 d cn3 dn+ F c 4cn 4) = 0. Equation 30) is valid for arbitrary arguments when all coefficients F 0,F cd,f c,f c 3 d,f c4 vanish simultaneously. This gives us a set of equations for the coefficients A, B, v F 0 = 440v 35α A m ) 60αAm ) [ β B 48m 5)+8 ] 3) + 4 [ 9β B 4 6m 46m+ ) +60β B 5m )+360 ] = 0, F cd = 30 m [ 9α A m )+3αA β B 75m 3)+ ) 4β B 9β B 5m )+ )] = 0, 3) F c = 5m 7α A m )+αa β B 59m 37)+6 ) 3β B 9β B m )+3 )) = 0, 33) F c 3 d = 90m3/ 3α A +86αAβ B 5β B 4) = 0, 34) F c 4 = 90m 3α A + 86αAβ B 5β B 4) = 0. 35) Equations 34) and 35) are equivalent and give the same condition as ). Solving 34) with respect to B, we obtain the same relations as in [0, Eq. 8)] B / ) = Aα β ) 36) Denote z = and z = ) 5 It is clear that z < 0 and z > 0. B has to be real-valued. This is possible for the case z=z if A<0, and for z=z if A>0. The value of z is the same as that found for the exact soliton solution in [0, Eq. 8)]. In general B = Aα z. 38) β Now, we insert 38) into 3),3) and 33). Besides a trivial solution with A=0 we obtain 440 v)+aα m)[080 35Aα m)] 39) 40Aα 5m) 30Aα) 0 09m+96m )z +4Aα) 9 847m+59m )z = 0, 9[Aαm )+4]+3[Aα75m 3) 6]z 40) 76Aα5m )z = 0, 93Aαm )+8)+Aα59m 37) 8)z 4) From 40) we find 608Aαm )z = 0. 4z 3) A= α[76z 5m ) z5m 93) 9m )] but from 4) it follows that 4) 44z 3) A= α[608z m ) 4z77m ) 7m )]. 43) This looks like a contradiction, but substitution z = z = )/5 in both 4) and 43) gives the same result A = ) )αm 5). 44) For z=z = )/5 the common result is A = ) )αm 5). 45) This means that not only are the equations 34) and 35) equivalent, but also 3) and 33), as well. Therefore the equations 3)-35) supply only three independent conditions for the coefficients of KdV solutions in the form 7). Now, using z=z and A given by 44) we obtain from 39) v = vnum m) vden m) 46)

5 Superposition solutions 5 and with z=z and A given by 45) v = vnum +m) vden + m), 47) where { vnum m)= ) 305 m ) 305 m } 305, and vden m)= ) m 5). Coefficients A,B,v of superposition solutions 7) to KdV as functions of m are presented in Fig. 3 for α = β = 0 and compared to corresponding solutions to KdV. Here, similarly as in Fig., we assume that A KdV = A KdV. Physically relevant exact superposition solutions to KdV can be found for greater values of the parameter α than 0. In figure 4 the amplitude Aα,m) given by 48) is shown as a contour plot for the region α [0.0,], m [0,]. It is clear that reasonable amplitudes occur in wide regions of α and m values, similarly like in the case dn. The amplitudes become too big only for α 0, but in such cases KdV works very well. Discussion of mathematical solutions.0. From a strictly mathematical point of view we found two families of solutions determined by coefficients A, B, v as functions of the elliptic parameter m. There are two cases. Case. z=z = < 0. This 5 case leads to B < 0 and has to be rejected as in previous papers [0, 8] As Bs Vs B V m Fig. 3 The same as in Fig. but for superposition solutions η + 7) and η 5). Case. z=z = > 0. Then ) 5 A= > 0, 48) 37αm 5) 305 ) 4 B= 49) 7035 m)β and v is given by 47). Since m [0,], m 5) < 0 then B is real. The solution in this case is η x v t,m)= A [ dn B x v t),m) 50) + mcnb x v t),m)dnb x v t),m) ]. Fig. 4 The same as in figure. but for superposition solutions η +..3 Superposition "dn mcndn" Now we check the alternative superposition "dn mcndn" η y)= A[ dn By,m) mcnby,m)dnby,m) ]. 5) In this case the derivatives are given by formulas similar to 8)-) with some signes altered. Analogous changes occure in formulas 3)-9). Then ) has a similar form like 30) AB m mcn+dn ) sn 5) F 0 + F cd cndn+ F c cn + F c 3 d cn3 dn+f c 4cn 4) = 0. Equation 5) is valid for arbitrary arguments when all coefficients F 0,F cd,f c,f c 3 d,f c 4 vanish simultaneously. This

6 6 Piotr Rozmej et al. gives us a set of equations for the coefficients v,a,b. Despite some changes in signs on the way to 5) this set is the same as for "dn + mcndn" superposition 3) 35). Then the coefficients A,B,v for superposition "dn mcndn" are the same as for superposition "dn + mcndn" given above. This property for KdV is the same as for KdV, see [5]. It follows from periodicity of the Jacobi elliptic functions. From cny+km), m)= cny, m), dny+km), m)= dny, m) it follows that dn y+km),m)+ m)cnx+km),m)dnx+km),m) = dn x,m) m)cnx,m)dnx,m). 53) So both superpositions η + 7) and η 5) represent the same solution, but shifted by the period of the Jacobi elliptic functions. This property is well seen in figures x single KdV sup+ KdV sup- KdV sup- KdV Fig. 5 Profiles of KdV and KdV waves for m= x single KdV sup+ KdV sup- KdV sup- KdV Fig. 6 Profiles of KdV and KdV waves for m= x single KdV sup+ KdV sup- KdV sup- KdV Fig. 7 Profiles of KdV and KdV waves for m= Examples Below, some examples of wave profiles for both KdV and KdV are presented. We know from section that for a given m, the coefficients A,B,v of KdV solutions are fixed. As we have already written, this is not the case for A,B,v of KdV solutions. So, there is one free parameter. In order to compare KdV solutions to those of KdV for identical m, we set A KdV = A KdV. In figures 5-7 below, KdV solutions of the forms 3), 7) and 5) are drawn with solid red, green and blue lines, respectively. For KdV solutions the same color convention is used, but with dashed lines. In all the presented cases the parameters α = β = 0. were used. Comparison of wave profiles for different m suggests several observations. For small m, solutions given by the single formula 3) differ substantially from those given by superpositions 7) and 5). Note that 3) is equal to the sum of both superpositions and when m the distance between crests of η + and η increases to infinity in the m= limit). All three solutions converge to the same soliton. In order to check whether the obtained analytic solutions are really true solutions to KdV several numerical simulations were performed. In each of them the numerical FDM code used with success in previous studies [9, 0,,4,8,9] was applied. Since the calculations concerned periodic solutions the periodic boundary conditions were used with an x interval equal to the particular wavelength. In figures 5-7 dashed lines display profiles of single dn 3) and superposition η + and η 4) solutions for three values of m = 0.,0.9 and Below in figures 6-8 six examples of time evolution for these solutions obtained in numerics are presented. Profiles of solutions at time instants t = 0,T/4,T/,3T/4 and T, where T = λ/v are displayed. Open symbols represent the profiles at t = T which overlap with those at t = 0 with numerical deviations less than 0. In all the presented examples, as well as all others not shown here, numerics confirmed a uniform motion and a fixed shape for the considered solutions.

7 Superposition solutions Fig. 8 Time evolution of solution red lines) and η + KdV solution for m= Fig. 9 The same as in figure 6 but for m= Conclusions The most important results of the paper can be summarized as follows. It is shown that several kinds of analytic solutions of KdV have the same forms as corresponding solutions to KdV but with different coefficients. This statement is true for our single solitonic solutions [0], periodic solutions in the form of single Jacobi elliptic functions cn [8] or dn, and for periodic solutions in the form of superpositions dn ± mcndn this paper). Coefficients A,B,v of these solutions to KdV are fixed by coefficients of the equation, that is by values of α,β parameters. This is in contradiction to the KdV case where one coefficient usually A) is arbitrary. References. Stokes, G.G.: On the theory of oscillatory waves. Camb. Trans. 8, ). Papers, 97-9.) Fig. 0 The same as in figure 6 but for m= Korteweg, D.J., de Vries, G.: On the change of form of the long waves advancing in a rectangular canal, and on a new type of stationary waves. Phil. Mag. 5), 39, 4 985). 3. Benjamin, T.B., Bona, J.L., Mahony, J.J.: Model equations for long waves in nonlinear dispersive media. Phil. Trans. Roy. Soc. London, Series A 7, ). 4. Bona, J.L., Pritchard, W.G., Scott L.R.: An evaluation of a model equation for water waves. Phil. Trans. Roy. Soc. London, Series A ). 5. Infeld, E., Rowlands, G., Senatorski, A.: Instabilities and oscillations in one and two dimensional Kadomtsev-Petviashvili waves and solitons. Proc. Roy. Soc. London A 455, ). 6. Marchant, T.R., Smyth, N.F.: The extended Korteweg de Vries equation and the resonant flow of a fluid over topography. J. Fluid Mech., ). 7. Marchant, T.R., Smyth, N.F.: Soliton Interaction for the Kortewegde Vries equation. IMA J. Appl. Math. 56, ). 8. Burde, G.I., Sergyeyev, A.: Ordering of two small parameters in the shallow water wave problem. J. Phys. A: Math. Theor ). 9. Karczewska, A., Rozmej, P., Rutkowski, Ł.: A new nonlinear equation in the shallow water wave problem., Physica Scripta, ). 0. Karczewska, A., Rozmej, P., Infeld, E.: Shallow-water soliton dynamics beyond the Korteweg - de Vries equation. Phys. Rev. E, 90, ).. Karczewska, A., Rozmej, P., Infeld, E.: Energy invariant for shallow water waves and the Korteweg de Vries equation: Doubts about the invariance of energy. Phys. Rev. E, 9, ).. Ali, A., Kalisch, H.: On the formulation of mass, momentum and energy conservation in the KdV equation. Acta Appl. Math., 33, ). 3. Sergyeyev, A., Vitolo, R.: Symmetries and conservation laws for the Karczewska-Rozmej-Rutkowski-Infeld equation. Nonl. Analys.: Real World Appl. 3, -9 06). 4. Karczewska, A., Rozmej, P., Infeld, E., Rowlands, G.: Adiabatic invariants of the extended KdV equation. Phys. Lett. A, 38, ). arxiv: Khare, A., Saxena, A.: Linear superposition for a class of nonlinear equations, Phys. Lett. A 377, ). 6. Khare, A., Saxena, A. Superposition of elliptic functions as solutions for a large number of nonlinear equations, J. Math. Phys. 55, ). 7. Khare, A., Saxena, A.: Periodic and hyperbolic soliton solutions of a number of nonlocal nonlinear equations, J. Math. Phys. 56, ).

8 8 Piotr Rozmej et al. 8. Infeld, E., Karczewska, A., Rowlands, G., Rozmej, P.: Solutions to the extended KdV equation for water surface waves. arxiv: Infeld, E., Karczewska, A., Rowlands, G., Rozmej, P.: Possible existence of a new class of shallow water waves. 07) Submitted. 0. Gardner, C.S., Greene, J.M., Kruskal, M.D., Miura, R.M.: Method for solving the Korteweg de Vries equation. Phys. Rev. Lett. 9, ).. Ablowitz, M.J, Clarkson, P.A.: Solitons, Nonlinear Evolution Equations and Inverse Scattering, Cambridge University Press, Cambridge 99).. Zabusky, N.J., Kruskal, M.D.: Interaction of "Solitons" in a Collisionless Plasma and the Recurrence of Initial States. Phys. Rev. Lett. 5, ).