Original article

On Lie symmetries and invariant solutions of Broer-Kaup-Kupershmidt equation in shallow water of uniform depth

  • Dig Vijay Tanwar , a, * ,
  • Mukesh Kumar , b
Expand
  • a Department of Mathematics, Graphic Era Deemed to be University, Dehradun 248002, India
  • b Department of Mathematics, Motilal Nehru National Institute of Technology Allahabad, Prayagraj 211004, India
*E-mail addresses: (D.V. Tanwar);
(M. Kumar)

Received date: 2022-02-28

  Revised date: 2022-04-12

  Accepted date: 2022-04-24

  Online published: 2022-04-29

Abstract

The dynamics of atmosphere and ocean can be examined under different circumstances of shallow water waves like shallow water gravity waves, Kelvin waves, Rossby waves and inertio-gravity waves. The influences of these waves describe the climate change adaptation on marine environment and planet. Therefore, the present work aims to derive symmetry reductions of Broer-Kaup-Kupershmidt equation in shallow water of uniform depth and then a variety of exact solutions are constructed. It represents the propagation of nonlinear and dispersive long gravity waves in two horizontal directions in shallow water. The invariance of test equations under one parameter transformation leads to reduction of independent variable. Therefore, twice implementations of symmetry method result into equivalent system of ordinary differential equations. Eventually, the exact solutions of these ODEs are computed under parametric constraints. The derive results entail several arbitrary constants and functions, which make the findings more admirable. Based on the appropriate choice of existing parameters, these solutions are supplemented numerically and show parabolic nature, intensive and non-intensive behavior of solitons.

Highlights

● Lie symmetry classification of Broer-Kaup-Kupershmidt equation in shallow water.

● The system describes the propagation of nonlinear and dispersive long gravity waves in two horizontal directions in shallow water.

● Symmetry reductions and invariant solutions are carried out.

● The obtained solutions are analyzed graphically.

● The solutions show parabolic nature, multi soliton, soliton fission phenomena.

Cite this article

Dig Vijay Tanwar , Mukesh Kumar . On Lie symmetries and invariant solutions of Broer-Kaup-Kupershmidt equation in shallow water of uniform depth[J]. Journal of Ocean Engineering and Science, 2024 , 9(3) : 199 -206 . DOI: 10.1016/j.joes.2022.04.027

1. Introduction

The study of solitary and shallow water waves has gained a great significance to describe characteristics of nonlinear wave phenomena. This nonlinear feature arises due to interaction of propagating high and low amplitude waves. Due to highly nonlinear behavior, such wave phenomena are expressed by nonlinear partial differential equations (NPDEs). These NPDEs most often illustrate the dynamics of shallow water waves, interaction of propagating solitary waves, flow of heat and mass, ion acoustic waves in plasma physics and electromagnetic waves in optical fibers etc. [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41]. Now a days, numerous great efforts have been spent to investigate exact solutions of NPDEs and to understand the dynamics of associated nonlinear physical models. Thus, several adequate tools notably variable separation approach [4], extended mapping approach [6], rational expansion method [7], extended Fan sub-equation method [8], [9], Exp-function method [10], [18], (GG) -expansion method [11], Hirota bilinear method [12], [16], Bäcklund transformation method [14], [26], modified Kudryashov method [21], [22], [23], modified Sardar sub-equation [24], generalized auxiliary equation method [25] and Lie symmetry approach [5], [15], [27], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41] have been evolved for deep observations of these NPDEs.
To reveal the significant features, we are interested to enlighten some new invariant solutions of Broer-Kaup-Kupershmidt (BKK) equation
utyuxxy+2uxuy+2uuxy+2vxx=0,
vt+vxx+2uvx+2vux=0,
where u and v represent the height and horizontal velocity of water waves respectively. It represents the propagation of nonlinear and dispersive long gravity waves in horizontal directions x and y in shallow water of uniform depth. The BKK equation can be deduced from Kadomtsev-Petviashvili (KP) equation using inner parameter symmetry constraints [1]. While x=y in Eq. (1.1) results in one dimensional Broer-Kaup equations, which describes dispersive long waves traveling in shallow water. The shallow water waves are the progressive waves, which classify marine environment, examine the ocean dynamics and model the equatorial tsunami waves. These waves arise due to action of wind on water and describe the long wavelength phenomena where horizontal length is greater than fluid depth. During propagation, these waves are affected by ocean bottom and thus orbital motion of water is disturbed. Therefore, it may cause underwater earthquakes, landslides, unimaginable destruction to the coastal ecology etc. Moreover, the internal waves arising from temperature variation can destruct marine vessels and can affect ecological as well as socio-economic environment.
A rich literature of BKK equation is reviewed here as Yu and Xin [2] analyzed PainlevȨ property and derived infinitely many symmetries of Eq. (1.1). Zhang et al. [3] used truncated PainlevȨ expansion method and obtained some soliton solutions of generalized BK equations. With the aid of variable separation approach, H.M. LI [4] found localized coherent structures of Eq. (1.1) and introduced its soliton solutions. H.C. Ma [5] derived Lie symmetries and then found exact dromion solutions. Thereafter, Fang and Zheng [6] studied BKK equation and derived nonelastic soliton solutions using extended mapping approach. Some soliton and periodic solutions were constructed by wang and Chen [7]. Zhang and Xia [8] proposed Fan sub-equation method and dummyTXdummy- used it to generated soliton and periodic solutions of Eq. (1.1). Continuing with same method, E. Yomba [9] investigated some general non-traveling wave solutions. Using Exp-function method via Riccati equation, some exact solutions of BKK equation were derived by S. Zhang [10]. Later on, Zhang et al. [11] employed (GG) -expansion method to Eq. (1.1) and listed some exact solutions. By means of symbolic computation and Hirota bilinear method, some soliton solutions were constructed by X.Y. Wen [12]. Moreover, Hu and Chen [13] derived nonlocal symmetries for BKK equation and studied its soliton and periodic solutions. By virtue of symbolic computation and Bell polynomial approach, Lan et al. [14] investigated some soliton solutions of Eq. (1.1). Furthermore, Kassem and Rashed [15] generated optimal system of Lie symmetries and derived some particular exact solutions. Some Wronskian solutions were introduced by Tang et al. [16] using Hirota method and bell polynomial approach. Recently, Kumar et al. used Lie symmetry method and found some exact solutions [17].
Being motivated with rich literature, we derive Lie symmetries and soliton solutions of BKK Eq. (1.1). Lie symmetry analysis under one parameter group of transformations ensures the invariance of system, which is used make the PDEs integrable and to reduce its order. It is also used to derive point symmetries of PDEs. Therefore, the given system of PDEs produces an equivalent system of ODEs and results in exact solutions. Based on characteristic of system and existing parameters, these solutions exhibit parabolic nature, soliton nature, intensive and non-intensive behavior of solitons. The solitons have wide applications in diverse physical contexts like plasma physics, quantum fluid dynamics and stratified fluid flows. Thus, Lie symmetry method serves as most applicable mathematical tool to find exact solutions having applications in various physical contexts.
The remaining texts are structured as: Section 2 presents the basic terms of symmetry method and invariance of BKK equation. The symmetry reductions and exact solutions are listed in Section 3. The derived results are analyzed in Section 4. The Section 5 is furnished with concluding remarks of findings in the end.

2. Lie symmetry analysis

The section briefly defines the termilogy to exatract infinitesimal generators. The detailed methodology can be reviewed from literature [27], [28] and existing references. Let us consider the one-parameter (ϵ ) transformations
xα˜=xα+ϵηα(xα,uβ)+O(ϵ2),uβ˜=uβ+ϵψβ(xα,uβ)+O(ϵ2),
where x1=t , x2=x , x3=y , u1=u , u2=v are variables, and η1=τ , η2=ξ , η3=η , ψ1=ϕ , ψ2=ψ are infinitesimals designated to t , x , y , u , v for maintaining the invariance of Eqs. Δ1=0 and Δ2=0 . The infinitesimal vector is v=τt+ξx+ηy+ϕu+ψv. Employing prolongation formulas Pr(3)v(Δ1)=0 and Pr(1)v(Δ2)=0 to Eq. (1.1), we have invariance condition as
ϕytϕxxy+2(uxϕy+uyϕx)+2(uxyϕ+uϕxy)+2ψxx=0,
ψt+ψxx+2(vxϕ+uψx)+2(uxψ+vϕx)=0,
where the extension coefficients ϕx , ϕy , ϕyt , ϕxxy , ψt , ψx , ψxx , ψxy are defined as ϕx=DxϕuxDxξuyDxηutDxτ,ϕy=DyϕuxDyξuyDyηutDyτ,ψt=DtψuxDtξuyDtηutDtτ,ψx=DxψuxDxξuyDxηutDxτ,ψxx=Dx2ψuxDx2ξuyDx2ηutDx2τ2uxxDxξ2uxyDxη2uxtDxτ,ψxy=DxDyψuxDxDyξuyDxDyηutDxDyτuxyDxξuxxDyξuyyDxηuxyDyηuytDxτuxtDyτ,ϕyt=DyDtϕuxDyDtξuyDyDtηutDyDtτuxtDyξuytDyηuttDyτuxyDtξuyyDtηuytDtτ,ϕxxy=Dx2DyϕuxxxDyξ2uxxyDxξ2uxyyDxη2uxytDxτuxxyDyηuxxtDyτ2uxxDxDyξ2uxyDxDyη2uxtDxDyτuxDx2DyξuxyDx2ξuyyDx2ηuytDx2τuyDx2DyηutDx2Dyτ, with total derivative operators Dx , Dy , Dt .
The Eqs. (2.1) and (2.2) along with given expressions yield determining equations as ξu=0,ξv=0,ξy=0,ηx=0,ηt=0,ηu=0,ηv=0,τx=0,τy=0,τu=0,τv=0,2ξxτt=0,2ϕ+uτtξt=0,2ψ+vτt+2vηy=0. This system generates the desired infinitesimals ξ=12xf2¯(t)+f3(t),η=f1(y),τ=f2(t),ϕ=14(2uf2¯(t)+xf2¯¯(t)+2f3¯(t)),ψ=12v(2f1¯(y)+f2¯(t)), with arbitrary functions f1(y) , f2(t) , f3(t) while the bar denotes derivative.
Hence, algebra of infinitesimal symmetries are generated by following vectors
v1(f1)=f1yvf1¯v,v2(f2)=x2f2¯x+f2tu2f2¯u+x4f2¯¯u12vf2¯v,v3(f3)=f3x+12f3¯u. The commutative relations of above vectors are determined as
[v1,v2]=0,[v1,v3]=0,[v2,v3]=v3(f2f312f2f3),[v2,v1]=0,[v3,v1]=0,[v3,v2]=v3(f2f312f2f3),[v1(g1),v1(h1)]=v1(g1h1h1g1),[v2(g2),v2(h2)]=v2(g2h2h2g2),[v3(g3),v3(h3)]=0, which constitutes infinite dimensional Lie algebra of symmetries.

3. Invariant solutions

To derive invariant solutions of BKK equation, the characteristic equation is furnished by
2dxxf2¯(t)+2f3(t)=dyf1(y)=dtf2(t)=4du2uf2¯(t)+xf2¯¯(t)+2f3¯(t)=2dvv(2f1¯(y)+f2¯(t)).
To continue integration further, the following two cases are taken under consideration as:
Case (I): If f2(t)=1a , where a0 , then Eq. (3.1) reforms to
dxf3(t)=dyf1(y)=adt1=2duf3¯(t)=dvvf1¯(y). Then, group invariants of BKK equation yield
u=a2f3(t)+G(r,s)andv=H(r,s)f1(y)
with symmetry variables r=xaf3(t)dtands=dyf1(y)at. Thereupon, first symmetry reduction of BKK equation results
Grrs+aGss2GrGs2GGrs2Hrr=0,
Hrr+2(GHr+HGr)aHs=0.
For further reduction, infinitesimal generators ξ˜ , η˜ , ϕ˜ and ψ˜ are
ξ˜=a1r+a2,η˜=2a1s+a3,ϕ˜=a1G,ψ˜=3a1H with constants a1 , a2 and a3 . So, one gets
dra1r+a2=ds2a1s+a3=dGa1G=dH3a1H.
Here, two intermediary cases are illustrated as:
Case (Ia ): For a10 , Eq. (3.5) turns to drr+A1=ds2(s+A2)=dGG=dH3H.(A1=a2a1andA2=a32a1) Eventually, the group invariants are G=G1(X)(s+A2)12andH=H1(X)(r+A1)3, with symmetry variable X=(r+A1)2s+A2 .
Substitution of above invariants into Eqs. (3.3) and (3.4) yields
4X2G1¯¯¯+12XG1¯¯+3G1¯a4(4X2G1¯¯+12XG1¯+3G1)4X12(XG1G1¯¯+2G1G1¯+XG1¯2)+4X52(2X2H1¯¯5XH1¯+6H1)=0,
2(2X2H1¯¯5XH1¯+6H1)6X12G1H1+4X32(G1H1¯+H1G1¯)+aX2H1¯=0.
The desired primitives are G1=a2X12andH1=c1X2
where c1 is constant.
Consequently, one derives the solutions
u=a2f3(t)a2(xaf3(t)dt+A1)(1f1(y)dyat+A2),
v=c1(xaf3(t)dt+A1)f1(y)(1f1(y)dyat+A2)2.
Case (Ib): For a1=0 and a30 , Eq. (3.5) reforms to
drA3=ds1=dG0=dH0.(A3=a2a3) It provides
G=G2(Y)andH=H2(Y) with variable Y=rA3s .
Availing above invariants into Eq. (3.6) and (3.7), we have
A3G2¯¯¯aA32G2¯¯2A3G2¯22A3G2G2¯¯+2H2¯¯=0,
H2¯¯+2(G2H2¯+H2G2¯)+aA3H2¯=0.
The primitives are read as G2=aA321Y+c3andH2=c2(Y+c3) which rises to
u=a2f3(t)aA321(xaf3(t)dtA31f1(y)dy+aA3t+c3),
v=c2f1(y)(xaf3(t)dtA31f1(y)dy+aA3t+c3).
Case (II): For f1(y)=b and f3(t)=cf2¯(t) , Eq. (3.1) reforms to
2dx(x+2c)f2¯(t)=dyb=dtf2(t)=4du(x+2c)f2¯¯(t)2uf2¯(t)=2dvvf2¯(t). Then, group invariants yields
u=(x+2c)f2¯4f2+G(r,s)f212(t)andv=H(r,s)f212(t)
with symmetry variables r=x+2cf212(t)ands=yb1f2(t)dt. Thus, first symmetry reduction produces
Grrs2GGrs2GrGs+bGss2Hrr=0,
Hrr+2(GHr+HGr)bHs=0.
To proceed further, infinitesimal generators ξ˜ , η˜ , ϕ˜ and ψ˜ are ξ˜=a4r+a5,η˜=2a4s+a6,ϕ˜=a4G,ψ˜=3a4H,
where a4 , a5 and a6 are constants.
Therefore, characteristic equation is
dra4r+a5=ds2a4s+a6=dGa4G=dH3a4H.
Again, two intermediary cases are illustrated as:
Case (IIa): For a40 , Eq. (3.17) generates drr+A4=ds2(s+A5)=dGG=dH3H
where A4=a5a4 and A5=a62a4 . Then, group invariants are obtained as
G=G1(X)(s+A5)12andH=H1(X)(r+A4)3 with symmetry variable X=(r+A4)2s+A5 .
Using above values into system (3.15) and (3.16), one gets
4X2G1¯¯¯+12XG1¯¯+3G1¯b4(4X2G1¯¯+12XG1¯+3G1)4X12(XG1G1¯¯+2G1G1¯+XG1¯2)+4X52(2X2H1¯¯5XH1¯+6H1)=0,
2(2X2H1¯¯5XH1¯+6H1)6X12G1H1+4X32(G1H1¯+H1G1¯)+bX2H1¯=0.
The primitives of these are G1=b2X12andH1=c4X2.(c4beingarbitraryconstant.)
Consequently, one derives the solution
u=(x+2c)f2¯(t)4f2(t)b[x+2c+A4f212(t)]2(ybdtf2(t)+A5),
v=c4[x+2c+A4f212(t)][ybdtf2(t)+A5]2f2(t).
Case (IIb): For a4=0 and a60 , Eq. (3.17) transforms to drA6=ds1=dG0=dH0.(A6=a5a6) It leads to
G=G2(Y)andH=H2(Y) with symmetry variable Y=rA6s .
Availing these invariants into Eqs. (3.18) and (3.19), we have
A6G2¯¯¯2A6(G2G2¯¯+G2¯2)bA62G2¯¯+2H2¯¯=0,
H2¯¯+2(G2H2¯+H2G2¯)+bA6H2¯=0.
It has primitives
G2=bA621Y+c6andH2=c5(Y+c6) and lead to
u=(x+2c)f2¯(t)4f2(t)bA62f212(t)1[x+2cA6(ybdtf2(t))f212(t)+c6f212(t)],
v=c5f212(t)[x+2cf212(t)A6(ybdtf2(t))+c6]
where c5 and c6 are constants.
A comparison of derived solutions (3.8), (3.9), (3.12), (3.13), (3.20), (3.21), (3.24) and (3.25) with previous results is presented here. The derived solutions involve the arbitrary functions f1(y) , f2(t) , f3(t) and arbitrary constants as well. The derived solutions are entirely new and general than previous finding [15], [17], [19] as
● Kassem and Rashed [15] derived their solutions by restricting the arbitrary functions and taking their values as constants while the solutions listed in this work contain several arbitrary functions.
● Kumar et al. [17] derived symmetry reductions and exact solutions for the single vector or combination of two vector fields only, while the present work deals with the general vector field v=v1(f1)+v2(f2)+v3(f3) . Therefore, the derived solutions involves more arbitrary functions existed in infinitesimals, which may be helpful to derives several physical phenomena and enhance the significance of results.
● Li et al. [19] derived the results involving arbitrary function of y only, while the solutions listed in the present work contain arbitrary function of y and t . Moreover, the present results are entirely different from Li et al. [19].
Thus, the comparison with previous results ensures the novelty of results and existence of arbitrary functions shows enrichment of present work.

4. Analysis and discussions

The graphical interpretation of invariant solutions encourages to reveal significant features physical phenomena. It demonstrates the precise informations to predict phenomena dynamics. This segment reveals the physical behavior of the phenomena associated with invariant solutions (3.8), (3.9), (3.12), (3.13), (3.20), (3.21), (3.24) and (3.25) using numerical simulation. The solutions involve arbitrary functions f1(y) , f2(t) , f3(t) and several arbitrary constants so simulation is performed for their appropriate values f1(y)=12(α+y) , α=0.8351 , a=0.5825 , b=0.8929 , c=0.7032 , A1=0.5827 , A2=0.8549 , A3=0.587 , A4=0.2176 , A5=0.251 and rest values are listed in adjacent figures. Due to existence of these functions and constants, the dynamical system has rich physical structures. Consequently, doubly soliton, multi-soliton, parabolic nature and soliton fission behavior have been analyzed. Solitons are nonlinear solitary waves, which generally preserve their identities after mutual collision and show elastic behavior. But for instance, interactions of solitons may not be elastic and result into soliton fusion or soliton fission. The detailed analysis of solution profiles is performed as follows:
Fig. 1: The surface plots of solution profiles u and v for Eqs. (3.8) and (3.9) are shown in this figure. The components u and v show doubly and multisolition profiles for f3(t)=2(β+t) , β=0.8954 , c1=0.0348 at t=0.7551 . The interaction of solitons is observed completely elastic with high peaks and deep tails under the plane wave.
Fig.1 Doubly soliton and multisoliton profile of u and v for expressions (3.8) and (3.9) at t=0.7551 .
Fig. 2: The spatio-temporal profiles for wave components u and v listed in Eqs. (3.12) and (3.13) are traced with adjacent figure. The values of functions and parameters are considered as f3(t)=2(β+t) , β=0.8954 , c2=0.9421 , c3=0.9564 at y=0.1319 . The component u shows soliton profile while v shows parabolic nature. Thus during propagation, the waves interact and result into nonlinear parabolic nature.
Fig.2 Soliton and parabolic nature for Eqs. (3.12) and (3.13).
Figs. 3 and 4: The Fig. 3 reveals the propagation of multisoliton and doubly soliton for Eqs. (3.20) and (3.21). The appropriate values of functions and constants are taken f2(t)=(γ+t) , γ=0.5557 , c4=0.1844 . The interaction of multisoliton and doubly soliton is observed at different values of t and shown at t=0.3174 in this figure, which results into multisoliton after mutual collision. The yu and yv views show that waves attain highest amplitude around origin. Moreover, treating a periodic function f2(t)=sec(δ+t) with δ=0.0319 , solution profile preserves the similar soliton nature shown in the Fig. 4. However, the amplitude of waves vary with respect to time, different values of functions and parameters as well.
Fig.3 Multisoliton and doubly soliton profiles for Eqs.(3.20) and (3.21).
Fig.4 Soliton profile for Eqs. (3.20) and (3.21).
Fig.5 Spatio-temporal profile for Eqs. (3.20) and (3.21).
Fig. 5: The spatio-temporal profile for Eqs. (3.20) and (3.21) is traced via this figure. It reveals soliton fission phenomena for appropriate values f2(t)=(γ+t) , γ=0.5557 , c4=0.1844 at y=0.9597 . The fission of one soliton into two is exhibited in tu view and shows new phenomena. Thus, it is pointed out that physical phenomena preserve its behavior with respect to spatial variables while the behavior is changed for space-time variables.
Fig. 6: The wave components expressed in Eqs. (3.24) and (3.25) are analyzed through this spatio-temporal surface plot, which reflects the elastic behavior multisoliton at y=0.7922 . The values to the constants are provided as f2(t)=sec(δ+t) , δ=0.0319 , c5=0.9594 , c6=0.6557 . The waves with different amplitudes interact during propagation and ultimately re-merge to its own shape, which is represented in tu views.
Fig.6 Multisoliton profile for the expressions (3.24) and (3.25).

5. Concluding remarks

In this work, we study the BKK equation and analyze its physical significance in the shallow water of uniform depth. The BKK equation illustrates the nonlinear and dispersive long gravity waves propagating in two horizontal directions. The invariance of system under one parameter transformations leads to infinitesimal generators. The commutative relations of infinitesimal vectors are derived, which constitute infinite dimensional Lie algebra of symmetries. Thereafter, Lie symmetry method is employed to derive symmetry reductions and generate exact solutions listed in Eqs. (3.8), (3.9), (3.12), (3.13), (3.20), (3.21), (3.24) and (3.25). These solutions retain arbitrary functions involved in infinitesimals as well as some free parameters. Due the existence of arbitrary functions and parameters, these finding have rich physical structures and can be used to analyze the dynamics of real phenomena. The physical structures of solutions are analyzed using numerical simulation, which show the nonlinear parabolic nature and characteristics of soliton interactions. Thus, Lie symmetry method can be treated as an effective tool to handle system of nonlinear PDEs.

Declaration of Competing Interest

Authors declare that they have no conflict of interest.

[1]
S.Y. Lou, X.B. Hu, J. Math. Phys., 38, (1997), pp. 6401-6427

[2]
H.Y. Ruan, Y.X. Chen, Acta Phys. Sin., 7 (1998), pp. 241-248

[3]
S.L. Zhang, B. Wu, S.Y. Lou, Phys. Lett. A, 300 (2002), pp. 40-48

[4]
H.M. Li, Commun. Theor. Phys., 39 (2003), pp. 513-518

[5]
H.C. Ma, Commun. Theor. Phys., 43 (2005), pp. 1047-1052

[6]
J.P. Fang, C.L. Zheng, Chin. Phys., 14 (2005), pp. 669-675

[7]
Q. Wang, Y. Chen. Chaos, Solitons Fractals, 30 (2006), pp. 197-203

[8]
S. Zhang, T. C. Xia. Appl. Math. Comput., 182 (2006), pp. 1651-1660

[9]
E. Yomba. Chaos, Solitons Fractals, 27 (2006), pp. 187-196

[10]
S. Zhang. Phys. Lett. A, 372 (2008), pp. 1873-1880

[11]
S. Zhang, W. Wang, J.L. Tong, Appl. Math. Comput., 209 (2009), pp. 399-404

[12]
X.Y. Wen, Nonlinear Anal., 12 (2011), pp. 3346-3355

[13]
H.R. Hu,.Y. Chen. Chin. Phys. B, 24 (2015), p. 090203 DOI: 10.1088/1674-1056/24/9/090203

[14]
Z.Z. Lan, Y.T. Gao, J.W. Yang, C.Q. Su, B.Q. Mao. Commun, Nonlinear Sci. Numer. Simul., 44 (2017), pp. 360-372

[15]
M.M. Kassem, A.S. Rashed. Chin. J, Phys., 57 (2019), pp. 90-104

[16]
Y. Tang, M. Yuen, L. Zhang. Appl. Math. Lett., 105 (2020), p. 106285

[17]
S. Kumar, S.K. Dhiman, A. Chauhan. Math. Comput. Simul., 196 (2022), pp. 319-335

[18]
L.K. Ravi, S.S. Ray, S. Sahoo. J. Ocean Eng. Sci., 2 (2017), pp. 34-46

[19]
L.Q. Li, Y.T. Gao, X. Yu, G.F. Deng, C.C. Ding. Int. J, Numer. Methods Heat Fluid Flow (2021), DOI:

[20]
M. Åženol, L. Akinyemi, H. Nkansah,W. Adel. J. Ocean Eng. Sci. (2022), DOI:

[21]
K.U. Tariq, M. Inc, A. Pashrashid, M. Zubair,L. Akinyemi. J. Ocean Eng. Sci. (2022), DOI:

[22]
A.H. Arnous, M. Mirzazadeh, L. Akinyemi,A. Akbulut. J. Ocean Eng. Sci. (2022), DOI:

[23]
K. Debin, H. Rezazadeh, N. Ullah, J. Vahidi, K.U. Tariq,L. Akinyemi. J. Ocean Eng. Sci. (2022), DOI:

[24]
L. Akinyemi, U. Akpan, P. Veeresha, H. Rezazadeh,M. Inc. J. Ocean Eng. Sci. (2022), DOI:

[25]
L. Akinyemi, M.M.A. khater, H. Rezazadeh. Opt. Quantum Electron., 54 (191) (2022)

[26]
S. Singh, S.S. Ray. J. Ocean Eng. Sci. (2022), DOI:

[27]
G.W. Bluman, J.D. Cole, Similarity Methods for Differential Equations. Springer-Verlag, New York (1974)

[28]
P.J. Olver, Applications of Lie Groups to Differential Equations. Springer-Verlag, New York (1993)

[29]
S.Y. Lou, H.C. Ma. Chaos, Solitons Fractals, 30 (2006), pp. 804-821

[30]
M. Kumar, D.V. Tanwar, R. Kumar. Comput. Math. Appl., 75 (2018), pp. 218-234

[31]
M. Kumar, D.V. Tanwar. Int. J, Geom. Methods Mod. Phys., 16 (1950110) (2019), pp. 1-14

[32]
M. Kumar, D.V. Tanwar, Pramana, 94 (23) (2020), pp. 1-10

[33]
D.V. Tanwar, A.M. Wazwaz. Eur. Phys, J. Plus, 135 (520) (2020)

[34]
D.V. Tanwar,M. Kumar. Nonlinear Dyn., 106 (2021), pp. 3453-3468 DOI: 10.1007/s11071-021-06954-7

[35]
D.V. Tanwar, Phys.Scr. 96 (2021), p. 065215

[36]
S.S. Ray, Math. Methods Appl. Sci., 44 (2021), pp. 8379-8396

[37]
D.V. Tanwar, A.M. Wazwaz. Phys. Scr.,95 (2020), p. 065220

[38]
R. Kumar, R.S. Verma, A.K. Tiwari, J. Ocean Eng. Sci. (2021), DOI:10.1016/j.joes.2021.12.005.

[39]
S. Kumar, S. Rani, J. Ocean Eng. Sci. (2021), DOI:10.1016/j.joes.2021.10.002.

[40]
S. Rani, S. Kumar, R. Kumar, J. Ocean Eng. Sci. (2021), DOI:10.1016/j.joes.2021.12.007.

[41]
D.V. Tanwar, A.K. Ray, A. Chauhan. Qual. Theory Dyn. Syst., 21, (2022), pp. 1-14

Outlines

/