Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

1 Issue per year

IMPACT FACTOR 2017: 0.755
5-year IMPACT FACTOR: 0.820

CiteScore 2017: 0.83

SCImago Journal Rank (SJR) 2017: 0.241
Source Normalized Impact per Paper (SNIP) 2017: 0.537

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

Optimal system, nonlinear self-adjointness and conservation laws for generalized shallow water wave equation

Dumitru Baleanu
  • Corresponding author
  • Cankaya University, Department of Mathematics, Ögretmenler Cad. 1406530, Ankara, Turkey
  • Institute of Space Sciences, Magurele, Bucharest, Romania
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Mustafa Inc / Abdullahi Yusuf / Aliyu Isa Aliyu
Published Online: 2018-06-26 | DOI: https://doi.org/10.1515/phys-2018-0049


In this article, the generalized shallow water wave (GSWW) equation is studied from the perspective of one dimensional optimal systems and their conservation laws (Cls). Some reduction and a new exact solution are obtained from known solutions to one dimensional optimal systems. Some of the solutions obtained involve expressions with Bessel function and Airy function [1,2,3]. The GSWW is a nonlinear self-adjoint (NSA) with the suitable differential substitution, Cls are constructed using the new conservation theorem.

Keywords: GSWW; optimal system; Cls; infinitesimal generators; NSA

PACS: 02.20.Qs; 02.20.Sv; 02.30.Hq; 02.30.Jr

1 Introduction

Lie symmetry methods play a vital role in the study and finding solutions for nonlinear partial differential equations (NLPDE) [4,5, 6,7,8,9,10,11,12,13,14,15,16]. Different techniques are used in the literature for the construction of Cls for different system of equation and these Cls are important for the investigation of a physical system [4,5, 6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]. Moreover, Authors made rigorous attempts for construction of the construction of one-dimensional and higher-dimensional optimal systems optimal system of Lie algebra [22,23,24,25,26]. Cls and symmetries have many application in science, physics and engineering [32,33,34,35,36,37].

In this work, we obtain one-dimensional optimal system, exact solutions and Cls for the GSWW equation given by


where a ≠ 0, b ≠ 0 are arbitrary constants. Eq. (1) have been studied by different authors using a variety of techniques. For example [27] introduced exact solutions for Eq. (1) by the general projective Ricatti equations method. Periodic wave solution for Eq. (1) by the improved Jacobi elliptic function method was investigated by [28] . Homogeneous balance method [29] was applied to investigate some solutions for Eq. (1) and some new solution of Eq. (1) with extended elliptic function method was proposed in [30] and many more.

2 One-dimensional optimal system of subalgebras of GSWW

In this section, we establish the optimal system of one-dimensional subalgebras of L4 and their corresponding exact solutions. Consider one parameter Lie group of the infinitesimal transformation below




where ϵ is the group parameter. The corresponding Lie algebra of the infinitesimal symmetries is the set of vector field of the form


Considering the fourth order prolongation Pr(4) of the vector field X such that


where Δ = (1), whenever Δ = 0. Using SYM package introduced in [31], the determining equations for Eq. (1) are obtained. Solving for η(x, t, u), ξ2(x, t, u), and ξ1(x, t, u) from the obtained determining equations, we get




where c1, c2, c3 and c4 are arbitrary constants and F(t) is an arbitrary function of t. The Lie symmetry algebra admitted by Eq. (1) is spanned by four infinitesimals generators below





The corresponding commutator table of the infinitesimal generators is given by

Table 1

Commutator table of the Lie algebra for GSWWE

2.1 Construction of one-dimensional optimal system of subalgebras

The Lie algebra L4 spanned by the given generators X1, X2, X3, and X4 can guarantee a possibility to obtain invariant solutions of Eq. (1). This will be based mainly on one-dimensional subalgebra of L4. There may be an infinite number of one-dimensional subalgebras of L4. Therefore, one can write an arbitrary generators from L4 as


which depend on the four arbitrary constants l1, l2, l3, and l4. We construct the one-dimensional optimal system of subalgebras using the method introduced in [22,23,24,25,26]. After the transformation of L4, we can get a 4-parameter group of linear transformations of the generators as


where l1, l2, l3, and l4 are the coefficients in Eq. (13).

2.2 Linear transformation

Here, we investigate the linear transformations by using their generators which is given as


and the structure constants of the Lie algebra L4 defined by cμvλ is given as


Consider the following cases:

  • Case 1: For μ = 1, v = 4, and λ = 1, 2 in Table 1. [X1, X4] = c141 X1 c142 X2 and the non vanishing structure constants are ( cμvλ ) are c141 = a, c142 = 2.

  • Case 2: For μ = 2, v = 4, and λ = 2 in Table 1. [X2, X4] = c242 X2 and the non vanishing structure constants ( cμvλ ) are c242 = −a.

  • Case 3: For μ = 3, v = 4, and λ = 3 in Table 1. [X3, X4] = c343 X3 and the non vanishing structure constants ( cμvλ ) are c343 = ∗a.

  • Case 4: For μ = 4, v = 1, and λ = 1, 2 in Table 1. [X4, X1] = c411 X1 c412 X2 and the non vanishing structure constants ( cμvλ ) are c411 = −a, c412 = −2. Setting v = 2, λ = 2 row four column two, we get [X4, X2] = c522 X2 and c422 = a, Setting v = 3, λ = 3 row four column three, we get [X4, X3] = c433 X3 and c433 = a.

Now, Substituting the values of the non-vanishing structure constants in Eq. (15) for μ = 1, 2, …, 4, we obtain


2.3 Lie equation

To obtain the Lie equation, we integrate the generators E1, E2, E3, E4 in Eq. (17) using the initial condition l|ϵ=0 = l.

  • For the generator E1, the Lie equation with the parameter ϵ are given by l1¯ϵ=al4¯,l2¯ϵ=2l4¯l3¯ϵ=0,l4¯ϵ=0. Integrating and using the initial condition we obtain l1 = 1 l4 + l1, l2 = 2ϵ1 l4 + l2, l3 = l3, l4 = l4.

    similarly for the other generators by following the same approach we get

  • For E2, we obtain l1 = l1, l2 = −2 l4 + l2, l3 = l3, l4 = l4.

  • For E3, we have l1 = l1, l2 = l2, l3 = −3 l4 + l3 and l4 = l4.

  • For E4, we have l1¯=l11+aϵ4,l2¯=2ϵ1+l21aϵ4,l3¯=l31aϵ4 and l4 = l4.

using SYM package [31], we can get the optimal system raw data in matrix form as:


and the number of the functionally invariants, which is found to be l4. The corresponding one-dimensional optimal system of subalgebras are found to be the following:

  1. X3,

  2. αX2 + X3,

  3. X4,

  4. X1 + X4,

  5. αX3 + X4,

  6. αX1 + βX2 + X3,

where α, β ∈ ℝ. In the following, we list the corresponding similarity variables, similarity solutions as well as the reduced PDEs obtained from the generators of optimal system and their exact solutions.

  1. Similarity variable related to X3 is u(x, t) = F(x) and F(x) satisfies Fxx = 0 two times integration implies that F(x) = c1 + xc2 and we have the exact solution


  2. Similarity variable related to αX2 + X3 is u(x, t) = αt + F(x) and F(x) satisfies Fxx = 0 which after integrating twice gives F(x) = c1 + xc2 and we have the exact solution


  3. Similarity variable related to X4 is u(x, t) = x2+aF(ζ)ax, ζ = tx and F(ζ) satisfies


    thrice integration of Eq. (20) and letting c1 = 0 yields


    solving for F(ζ) in Eq. (21) we obtain




    Hence by back substituting the similarity variables we get the exact solution as




    L1, L2 are as stated above and ζ = tx.

  4. Similarity variable related to X1 + X4 is u(x, t) = x2+F(ζ)1+ax, ζ = t(1 + ax) and F(ζ) satisfies the following


    three times integration and letting c1 = c2 = c3 = 0 in the above equation, gives


    solving for F(ζ) we have the following equation in Bessel function form




    and the exact solution is


    where ζ = t(1 + ax).

  5. Similarity variable related to βX3 + X4 is u(x, t) = x2+aF(ζ)ax,ζ=(atα)xa and F(ζ) satisfies the following


    here the reduced PDE is the same as that in reduction 3, the only different is the variable ζ. Therefore, we get




    and L1, L2, Q0, Q1, Q2, and Q3 are as stated in Eq. (23) and ζ = (atα)xa.

  6. Similarity variable related to αX1 + βX2 + X3 is u(x, t) = βx+αF(ζ)α,ζ=αtxα and F(ζ) satisfies the following


    three times integration of Eq. (30) with c1 = 0 leads


    solving for Fζ, we obtain




    Thus, by back substituting the similarity variables we get


    where ζ=αtxα.

3 Physical interpretation of the solutions (23) and (33)

In order to have clear and proper understanding of the physical properties of the power series solution, the 3-D, 2-D and contour plots for the solution Eqs. (23) and (33) are plotted in Figures 1-4 by using suitable parameter values.

3D plot of (33) α = 4, a = 1, b = 2.5, c4 = 13
Figure 1

3D plot of (33) α = 4, a = 1, b = 2.5, c4 = 13

contour plot of (33) α = 80, c1 = 5, c2 = 0.5
Figure 2

contour plot of (33) α = 80, c1 = 5, c2 = 0.5

3D plot of (23) α = 80, β = 10, a = 3, b = 2.5, c3 = 0.5, c2 = c4 = 10
Figure 3

3D plot of (23) α = 80, β = 10, a = 3, b = 2.5, c3 = 0.5, c2 = c4 = 10

contour plot of (23) α = 80, β = 10, a = 3, b = 2.5, c3 = 0.5, c2 = c4 = 10
Figure 4

contour plot of (23) α = 80, β = 10, a = 3, b = 2.5, c3 = 0.5, c2 = c4 = 10

4 Nonlinear self-adjointness

The system of m [20, 21] differential equations


α = 1, …, m, with m dependent variables u = (u1,…,um) is said to be NSA if the adjoint equations


α = 1, …,m, are satisfied for all solutions u of the original system Eq. (34) upon a substitution


α = 1, …, m, such that


On the other hand, the equation below holds:


α = 1, …, m, where λαβ¯ are undetermined coefficients, and φ(σ) are derivatives of Eq. (36),


,σ = 1, …, s. Here v and φ are the m-dimensional vectors v = (v1, …, v(m)), φ = (φ1, …, φm), and also, it is worth noting that not all components φα(x, u) of φ vanish simultaneously from Eq. (37).

4.1 Test for self-adjointness for GSWW

Here, we want to test the self-adjointness of Eq. (1). The adjoint equation for Eq. (1) is given by


Let v = ϕ(x, t, u), after some calculations and equating the coefficients of the derivatives ut, ux, uxt, uxx, uxxt and uxxx to zero, we have


The solution for ϕ(t, x, u) from the above equation is simply found to be the following


where C1 is an arbitrary constant. Therefore, Eq. (1) is NSA with the substitution in Eq. (40).

5 Conservation laws for GSWW

In this section, we establish Cls for Eq. (1) [19,20,21].

The reality that Eq. (1) is NSA with the obtained differential substitution in Eq. (40), we can use the Noether operator 𝓝 to obtain its conserved vectors (C1, C2) [17,18,19,20]. The obtained conserved vectors will satisfy the conservation equation DxC1 + DtC2 = 0. Moreover, the non local variables appearing in that formula must be substituted according to equation (40).

The following are the conserved vectors obtained from the four infinitesimals say X1, X2, X3, and X4 when C1 = 1 respectively.



Similarly, one can verify and see that DxC1 + DtC2 = 0 which is a trivial conservation laws.








6 Conclusion

In this study, Lie symmetry analysis, one dimensional optimal system and Cls for GSWW equation were studied. Some reductions and their solutions were reported from the obtained one dimensional optimal system. We presented some figures for some of the obtained exact solutions. The exact solutions include an expression with a Bessel function and an Airy function. We verified the authenticity of the solution by substitution into the original equation. The GSWW equation is a NSA with the obtained differential substitution, we obtained Cls using the new conservation theorem presented by Ibragmov.


  • [1]

    Temme N.M., Special Functions: An Introduction to the Classical Functions of Mathematical Physics, 2011, John Wiley and Sons,  CrossrefGoogle Scholar

  • [2]

    Vallee O., Soares M., Airy Functions and Applications to Physics, 2004, Imperial College Press Google Scholar

  • [3]

    Temme N.M., Special Functions: An Introduction to the Classical Functions ofMathematical Physics, 1996, John Wiley and Sons Google Scholar

  • [4]

    Emrullah Y., San S., Özkan Y.S. “Nonlinear self ad-jointness, conservation laws and exact solutions of ill-posed Boussinesq equation, Open Physics, 2016, 14(1), 37-43. Google Scholar

  • [5]

    San S., Akbulut A., Ünsal Ö, Tascan F. Conservation laws and double reduction of (2+1) dimensional Calogero- Bogoyavlenskii-Schiff equation, Math. Meth. Appl. Sci., 2017, 40(5), 1703-1710. CrossrefWeb of ScienceGoogle Scholar

  • [6]

    Emrullah Y., San S., A procedure to construct conservation laws of nonlinear evolution equations, Zeitschrift für Naturforschung A, 2016, 71(5), 475-480. Web of ScienceGoogle Scholar

  • [7]

    Inc M., Aliyu A.I., Yusuf A., Dark optical, singular solitons and conservation laws to the nonlinear Schrödinger’s equation with spatio-temporal dispersion, Mod. Phys. Let. B, 2017, 31(14), 1750163 CrossrefGoogle Scholar

  • [8]

    Tchier F., Yusuf A., Aliyu A.I., Inc M., Soliton solutions and Conservation laws for Lossy Nonlinear Transmission line equation, Superlat. Microstr., 2017, 107, 320-336. CrossrefGoogle Scholar

  • [9]

    Inc M., Aliyu A.I., Yusuf A., Traveling wave solutions and conservation laws of some fifth-order nonlinear equations, Eur. Phys. J. Plus, 2017, 132, 224. Web of ScienceCrossrefGoogle Scholar

  • [10]

    Bokhari A.H., Al Dweik A.Y., Kara A.H., Zaman F.D. A symmetry analysis of some classes of evolutionary nonlinear (2+1)- diffusion equations with variable diffusivity, Nonlinear Dyn., 2010, 62, 127-138. Web of ScienceCrossrefGoogle Scholar

  • [11]

    Krishnan E.V., Triki H., Labidi M., Biswas A., A study of shallow waterwaveswith Gardner’s equation, Nonlinear Dyn, 2011, 66, 497-507. CrossrefGoogle Scholar

  • [12]

    Biswas A., Kara A.H., 1-Soliton solution and conservation laws of the generalized Dullin-Gottwald-Holm equation, Appl. Math. Comp., 2010, 217, 929-932. CrossrefGoogle Scholar

  • [13]

    Biswas A., Kara A.H., Moraru L., Bokhari A.H., Zaman F.D., Conservation Laws of Coupled Klein-Gordon equation with cubic and power law nonlinearities, The publishing house proceeding of the Romanian Academy, Series A, 2014, 15(2), 123-129. Google Scholar

  • [14]

    Ebadi G., Kara A.H., Petkovic M.D., Biswas A., Soliton solutions and conservation laws of the Gilson-Pickering equation, Waves in Random and Complex Media, 2011, 21(2), 378-385. Web of ScienceCrossrefGoogle Scholar

  • [15]

    Morrisa R., Kara A.H., Biswas A., Soliton solution and conservation laws of the Zakharov equation in plasmas with power law nonlinearity, Nonlinear Analysis: Modelling and Control, 2013, 18(2), 153-159. Google Scholar

  • [16]

    Razborova P., Kara A.H., Biswas A., Additional conservation laws for Rosenau-KdV-RLW equation with power law nonlinearity by Lie symmetry, Nonlinear Dyn., 2015, 79, 743-748. CrossrefWeb of ScienceGoogle Scholar

  • [17]

    Noether E., Invariant variation problem, Mathematisch-Physikalische Klasse, 1918, 2, 235-257. Google Scholar

  • [18]

    Khamitova R., Symmetries and Conservation laws, PhD thesis, 2009, Växjö, Sweden Google Scholar

  • [19]

    Ibragimov N.H., A new Conservation laws theorem, J. Math. Anal., 2007, 333(1), 311-328. CrossrefGoogle Scholar

  • [20]

    Ibragimov N.H., Nonlinear self-adjointness in constructing conservation laws, Archives of ALGA, 2010, 7/4, 1-18. Google Scholar

  • [21]

    Ibragimov N.H., Nonlinear self-adjointness in constructing conservation laws, Archives of ALGA, 2011, 44(43), 2011 Google Scholar

  • [22]

    Ovsyannikov L.V., Group analysis of differential equations, Nuaka, Moscow, 1978, English transl., Ames W.F. (Ed.), 1982, Academic Press, New York. Google Scholar

  • [23]

    Ovsyannikov L.V., Group properties of differential equations, 1962 Siberian Branch, USSR Academy of Sciences, Novosibirirsk, (Russian) Google Scholar

  • [24]

    Patera J., Winternitz P., Zassenhaus H., Continuous subgroups of the fundamental groups of physics, General method and the Poincare group, J. Math. Phys., 1975, 16, 1597-1614. CrossrefGoogle Scholar

  • [25]

    Patera J., Sharp R.T., Winternitz P., Zassenhaus H., Invariants of real low dimension Lie algebras, J. Math. Phys., 1976,17, 986-994. CrossrefGoogle Scholar

  • [26]

    Galas F., Richter E.W., Exact similarity solutions of ideal MHD equations for plane motions, Phys. D, 1991, 50, 297-307. CrossrefGoogle Scholar

  • [27]

    Gomez C.A., Salas A., Exact solutions for the generalized shallow water wave equation by the general perspective Ricatti equations method, Math. Phys., 2006, 2006, 50-56. Google Scholar

  • [28]

    Inc M., Ergut M., Periodic wave solutions for the generalized shallow water wave equation by the improved Jacobi elliptic function method, Appl. Math. E-Notes, 2005, 5, 89-96. Google Scholar

  • [29]

    Elwakil S.A., El-labany S.K., Zahran M.A., Sabry R., Exact travelling wave solutions for the generalized shallow water wave equation, Chaos, Solitons & Fractals, 2003, 17(1), 121-126. CrossrefGoogle Scholar

  • [30]

    Bagchi B., Das S., Ganguly A., New exact solutions of a generalized shallow water wave equation, Phys. Scr., 2010, 82(2), 025003 CrossrefWeb of ScienceGoogle Scholar

  • [31]

    Dimas S., Tsoubelis D. SYM: A new symmetry finding pack- age for Mathematica, In: Ibragimov N.H., Sophocleous C., Damianou P.A., (Eds.), The 10th International Conference in Modern Group Analysis, 2005, 64-70, Nicosia, University of Cyprus Google Scholar

  • [32]

    Bilby B.A., Miller K.J., Willis J.R., Fundamentals of Deformation and Fracture, 1985, Cambridge Univ. Press, Cambridge Google Scholar

  • [33]

    Knops R.J., Stuart C.A., Quasiconvexity and uniqueness of equilibrium solutions in nonlinear elasticity, Arch. Rat. Mech. Anal., 1984, 86, 234-249. Google Scholar

  • [34]

    Harwitt M., Photon orbital angular momentum in astrophysics, Astrophys. J., 2003, 597, 1266-1270. CrossrefGoogle Scholar

  • [35]

    Elias N.M., Photon orbital angular momentum in astronomy, Astron. Astrophys, 2008, 492(3), 883-922. CrossrefWeb of ScienceGoogle Scholar

  • [36]

    Berkhout G.C.G., Beijersbergen M.W., Method for probing the orbital angular momentumof optical vortices in electromagnetic waves from astronomical objects. Phys. Rev. Lett., 2008, 101, 100801 CrossrefGoogle Scholar

  • [37]

    Thidé B., Then H., Sjöholm J., Palmer K., Bergman J., Carozzi T.D., Istomin Y.N., Ibragimov N.H., Khamitova R., Utilization of photon orbital angular momentum in the low-frequency radio domain, Phys. Rev. Lett., 2007, 99, 087701-1-087701-4 Web of ScienceGoogle Scholar

About the article

Received: 2017-08-19

Accepted: 2018-01-09

Published Online: 2018-06-26

Citation Information: Open Physics, Volume 16, Issue 1, Pages 364–370, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0049.

Export Citation

© 2018 Dumitru Baleanu et al.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in