Research article

Dynamics of some new solutions to the coupled DSW equations traveling horizontally on the seabed

  • Raj Kumar , * ,
  • Ravi Shankar Verma
Expand
  • Department of Mathematics, Faculty of Engineering & Technology, Veer Bahadur Singh Purvanchal University, Jaunpur 222003, India
*E-mail addresses: (R. Kumar)
(R.S.Verma).

Received date: 2022-02-12

  Revised date: 2022-04-15

  Accepted date: 2022-04-18

  Online published: 2022-04-22

Abstract

The system of (1+1 )-coupled Drinfel'd-Sokolov-Wilson equations describes the surface gravity waves travelling horizontally on the seabed. The objective of the present research is to construct a new variety of analytical solutions for the system. The invariants are derived with the aid of Killing form by using the optimal algebra classification via Lie symmetry approach. The invariant solutions involve time, space variables, and arbitrary constants. Imposing adequate constraints on arbitrary constants, solutions are represented graphically to make them more applicable in designing sea models. The behavior of solutions shows asymptotic, bell-shaped, bright and dark soliton, bright soliton, parabolic, bright and kink, kink, and periodic nature. The constructed results are novel as the reported results [26,28,29,30,33,38,42,49] can be deduced from the results derived in this study. The remaining solutions derived in this study, are absolutely different from the earlier findings. In this study, the physical character of analytical solutions of the system could aid coastal engineers in creating models of beaches and ports.

Highlights

● System of coupled Drinfel'd-Sokolov-Wilson equations describes waves travelling horizontally on the seabed.

● The killing form is used to get the invariants.

● One dimensional optimal sub algebras are classified and analytical solutions are attained via one parameter Lie symmetry analysis.

● The solutions are asymptotic, bell-shaped, bright and dark soliton, bright soliton, parabolic, bright and kink, kink, and periodic in nature.

Cite this article

Raj Kumar , Ravi Shankar Verma . Dynamics of some new solutions to the coupled DSW equations traveling horizontally on the seabed[J]. Journal of Ocean Engineering and Science, 2024 , 9(2) : 154 -163 . DOI: 10.1016/j.joes.2022.04.015

1. Introduction

1.1. Scope

In the physical sciences, nonlinear partial differential equations (NPDEs) are frequently employed to describe complicated processes. Thus, the problems/models arising in ocean science [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], oceanography [24], mathematical physics [25], complex fluid flows [26], [27], plasma physics [28], electromagnetic theory, fluid dynamics [29], [30], nuclear physics, chemical physics, to name a few [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [58], [59], [60], [61], [62], [63], [64], [65], [66], [67], [68], [69] are not easily solvable and finding their analytical solutions is critical. Reviewing the existing literature [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50], the authors are motivated to solve the following (1+1 )-coupled Drinfel'd-Sokolov-Wilson equations (CDSWEs) system analytically. The system is governed by
ut+avvx=0,vt+bvxxx+cuvx+dvux=0,
where a , b , c and d are non-zero parameters, x , t are space and time variables, respectively, while u and v depict the components of nonlinear surface gravity waves travelling horizontally on the seabed.
Gravity waves are formed in a fluid medium or at the interface of two media when gravity or buoyancy trying to restore equilibrium. The contact between the atmosphere and the ocean, which causes wind waves, is an example of such an interaction. In deeper water, a long water wave moves faster. In deeper water, gravity waves have a faster phase speed than in shallow water. Shallow water waves (SWW) are helpful to classify the marine environment, investigate ocean dynamics, and model equatorial tsunami waves. During their propagation, SWW are influenced by the ocean floor, causing the orbital motion of water to be disrupted. As a result, they may cause underwater earthquakes, and unimaginable damage to the coastal ecology. On the other hands, the heights, wave lengths, and time durations of gravity waves frequently vary as they advance in different directions [1]. Gravity waves in shallow water are dispersion (water waves) and non-dispersive as the depth is substantially smaller than the wavelength. Nonlinear interactions between triads of wave components with frequencies (ω) and vector wave numbers (k) satisfying following conditions which affect ocean surface gravity waves [2].
|ω1±ω2±ω3|=δω,|k1±k2±k3|=δk.

1.2. Origin of the problem

It is essential to consider the following form Hirota et al. [3] of coupled Korteweg-de Vries (KdV):
2ut(uxxx+6uux)=4bvvx,vt+vxxx+3uvx=0,
which represents the interaction of two long waves with distinct dispersion relations, in which b is a non-zero parameter and u and v are similar as described in CDSWEs (1).
The system (3) is integrable and can be reformulated as six-reductions of the Kadomtsev-Petviashvili hierarchy. The general construction of the KdV system (3) can be deduced using affine Lie algebras [3], while Drinfel'd and Sokolov [4] established that the general construction of the KdV system (3) can be derived using affine Lie algebras (ALA). Wilson established a link between the Kac-Moody Lie algebra [5] and KdV (3). Wilson [6] used ALA to obtain a particular variant of Eq. (1) with a=3 , b=2 , c=2 , and d=1 as inputs. One of the universal models is the CDSWEs (1), which have an infinite number of conservation laws and an integrable system [3].

1.3. Literature review

The authors are inspired by the characteristics of CDSWEs (1). The research community [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50] solved the system (1) using different approaches and addressing different concerns. Singh et al. [25] employed a fractional time derivative and a homotopy perturbation technique to get a solitary wave solution in the Caputo sense. The travelling wave solutions were obtained by Khan et al. [26] using a modified simple equation approach. Lu et al. [27] have employed first integral technique and obtained exponential and rational functions types solutions. Homotopy analysis has been used by Arora and Kumar [28] to establish a convergent series with two initial conditions. To get a solitary wave solution, Arnous et al. [29] applied two approaches to obtain a solitary wave solution which were the Bäcklund transformation of the Riccati equation and the trial function approach. Using the F-expansion method, Akbar and Ali [30] derived hyperbolic, rational, single soliton, and periodic solutions.
For solving time-fractional DSW systems and obtaining a solution in the form of a truncated series, Chen et al. [31] compared the coupled fractional reduced differential transform method with the residual power series method. Tariq and Seadawy [32] tackled the equation using an auxiliary equation and derived traveling and solitons solutions. The Jacobi elliptic function method was employed by Yao [33] and attained kink, bell-shaped, singular, and periodic solutions. Lu et al. [34] exploited the extended auxiliary and equation mapping methods to obtain a solution in the form of a moving and solitary wave for the same system of equations with odd and even order partial derivatives. For the CDSWEs (1), Wen et al. [35] applied bifurcation and qualitative analysis of dynamical systems to obtain solitary wave, periodic, kink shaped solutions.
Javeed et al. [36] used the first integral approach to construct a nontrivial solution to a coupled space-time fractional DSW system with Riemann-Liouville and conformable types of derivatives. For the space-time fractional derivative of the classical DSW system, Shehata et al. [37] employed the Kudryashov method to obtain a travelling wave solution. In order to determine unique periodic and solitary wave solutions for CDSWEs (1), Zhang [38] employed a semi-inverse variational approach, whereas Bibi et al. [39] used the Tanh and expanded Tanh methods and generated kink and soliton solutions. Bhatter et al. [40] exploited fractional homotopy analysis and obtained bell-shaped solutions to the nonlinear CDSWEs using a fractional operator. The GG -expansion method was used by Matjila et al. [41], Shi et al. [42], and Khan et al. [43] who obtained different forms of the solutions like travelling waves, hyperbolic, and trigonometric, rational, bell-shaped soliton, kink, and solitary waves.
Hirota et al. [3] established a special type of solitons for CDSWEqs. (1) i.e. static solitons, if a=3 , b=c=2 , and d=1 . They [3] claimed that when using Painlevé to deform CDSWEs into Bilinear form, such solitons can not be deformed when interacting with other moving solitons, whereas Ali Akbar et al. [44] used the modified alternative GG -expansion method to extract travelling wave solutions for the same form of equations.
Using a=3 , b=c=2 , and d=1 in system (1), Zhao and Zhi [45] utilised the F-expansion approach to get Jacobi elliptic, rational, solitary, and periodic solutions; and Ullah et al. [46] used the optimal homotopy method to provide doubly periodic wave solutions. The Adomian decomposition method was used by Inc [47] to determine the system's Jacobi elliptic solution
A different form of CDSWEs is examined for specific values of a=2 , b=a , c=3k , and d=3b in Eq. (1), Morris and Kara [48] established a travelling wave solution using Lie symmetry and conservation laws. Using Ibragimov's technique, Zhang and Zhao [49] created a one-dimensional optimal system. They [49] solved the additional reductions using the simple equation approach and then extracted the polynomial solutions, whereas Zhao et al. [50] applied Lie symmetry analysis and conservation laws to the identical form of DSWEs but did not solve the reductions further. This paved the way for more solutions of DSWEs.

1.4. Motivation

Making use of the similarity transformations method (STM) via one-parameter Lie symmetry analysis, the authors derived several types of novel analytical solutions. Sophus Lie created the Lie symmetry analysis as an ad-hock integration technique for the system of PDEs. The system of differential equations is still invariant after each similarity reduction. In a system of PDEs, the STM decreases the number of independent variables by one. Using the STM repeatedly can obviously reduce the system to an equivalent system of ordinary differential equations (ODEs). Due to the invariance of STM, a system of PDEs might result in a linear system of new over-determining system of PDEs that enables one to find infinitesimal generators that are functions of the independent and dependent variables. Such adjustments, which employ Lagrange's equation as a tool, result in similarity variables and similarity functions. Substituting the value of the dependent variable in terms of a similarity function can reduce the given system of PDEs to a new system. This process is repeated until one gets a system of ODEs. On solving it, it provides solutions to the system. For a description of the STM and its uses, one might look through the extensive literature [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], [42], [43], [44], [45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [58], [59], [60], [61], [62], [63], [64], [65], [66], [67], [68], [69] and the references therein.

1.5. Outline

This article is structured as follows: In the next section, invariants are derived by killing form and then one-dimensional optimal sub algebras are generated. In Section 3, a new class of analytical solutions is constructed by using Lie symmetry analysis. Section 4 depicts the comparison with the reported results till now. Section 5 comprises a detailed discussion about the nature of solutions and their physical analysis. Finally, a summary of the whole article is given in the conclusion.

2. Invariants by using Killing form

The essential terminology for the methodology is presented in the following. The one-parameter (ϵ) Lie symmetry transformations can be considered as:
x*=x+ϵξ(χ)+O(ϵ2),t*=t+ϵτ(χ)+O(ϵ2),u*=u+ϵη(1)(χ)+O(ϵ2),andv*=v+ϵη(2)(χ)+O(ϵ2)
where ξ , τ , η(1) and η(2) are the infinitesimals of the variables x , t , u , and v respectively. The notation x* means the value of x in transformed space. The notation (χ) is equivalent to (x , t , u , v ). The vector field V generated by the infinitesimal transformations
V=ξx+τt+η(1)u+η(2)v.
The invariance conditions (termed in Bluman and Cole [68], Olver [69]) for the CDSWEs (1) can be written as:
Pr(1)V[ut+avvx]=0,andPr(3)V[vt+bvxxx+cuvx+dvux]=0,
where the first prolongation Pr(1) and the third prolongation Pr(3) are
Pr(1)V=V+[ηt(1)](ut)+[ηx(2)](vx),andPr(3)V=V+[ηt(2)](vt)+[ηx(1)](ux)+[ηx(2)](vx)+[ηxxx(2)](vxxx),
where
[ηt(1)]=η(1)t+η(1)uθt+[η(1)vτt]ϕtξtϕxξuθtϕxτuθtϕtξvϕxϕtτvϕt2and[ηx(2)]=η(2)x+η(2)uθx+[η(2)vξx]ϕxτxϕtξuθxϕxτuθxϕtξvϕx2τvϕxϕt,etc.
are extensions.
Then over-determining equations are
ξt=ξu=ξv=τx=τtt=τu=τv=0,ξx=13τt,η(1)=23uτt,andη(2)=23vτt.
By solving them, one can get
ξ=a1+a2x,τ=a3+3a2t,η(1)=2a2u,andη(2)=2a2v,
where a1 , a2 , and a3 are arbitrary constants.
The Lie symmetry algebra for CDSWEs (1) can be generated by
X1=x,X2=xx+3tt2uu2vv,andX3=t.
Lemma 2.1. Let X=i=13aiXi be an arbitrary element of L3 , where a1 , a2 , a3 R . Then Symmetry algebra admits an arbitrary invariant function of the form F(a2) , where F is an arbitrary function.
Proof. Let ψ be an invariant function in Lie algebra L3 then, ψ[Ad(g)X]=ψ(X) , XL3 and gG , where G stands for the symmetry group. Now X=i=13aiXi and W=i=13biXi , then
Ad(exp(ϵW))X=eϵWXeϵW=Xϵ[W,X]+ϵ22[W,[W,X]]=(a1X1+a2X2+a3X3)(Θ1X1+Θ2X2+Θ3X3)+0(ϵ2).
Now, with the assistance of a commutator Table 1, one can get
Θ1=b1a2b2a1,Θ2=0,andΘ3=3b2a3+3b3a2.
Then the corresponding system of linear DEs recasts as:
a2ψa1=0,a1ψa1+3a3ψa3=0,and3a2ψa3=0.
Table 1 Commutator table for (1).
*X1 X2 X3
X1 0X1 0
X2 X1 03X3
X3 03X3 0
Lemma 2.2. The killing form is K(X,X)=10a22 , where X=i=13aiXi for Lie algebra L3 .
Proof. The Lie algebra's killing form is provided by K(X,X)=Trace[adXadX] , where
AdX=(a2a1000003a33a2)
So, the killing form is K(X,X)=10a22 . □
The adjoint relation is calculated using a Lie series that follows
Ad(exp(ϵW))X=Xϵ[W,X]+12ϵ2[W,[W,X]]
Lemma 2.3. B=a2 is an invariant.
Proof. Table 3 is a direct consequence of the proof. □
Table 2 Adjoint table for (1).
Ad X1 X2 X3
X1 X1 X2 - ϵX1 X3
X2 eϵX1 X2 e3ϵX3
X3 X2 X23ϵX3 X3
Table 3 Construction of invariant function for (1).
Coeff of X1 X2 X3
Ad(expϵX1)X a1a2ϵ a2 a3
Ad(expϵX2)X eϵa1 a2 e3ϵa3
Ad(expϵX3)X a1 a2 a33ϵa2
Lemma 2.4. I is an invariant, where
I={1,ifa22+a3200otherwise.
Proof. Table 3 shows that the adjoint action Ad(exp(ϵX1)) keeps invariant to the coefficients of the X2 and X3 . Hence I is an invariant.
a¯22+a¯32=a22+(a3e3ϵ)2
a¯22+a¯32=a22+(a33ϵa2)2
From Eqs. (12), (13), it is clear that a¯2=a¯3=0 iff a2=a3=0 , where a¯i , i=2 , 3 is an adjoint action. □
Lemma 2.5. Let
S={sgn(a3),a30,a2=00otherwise
where S is a sgn function.
Then S(X)=S(Ad(exp(ϵX))), for any vector X=i=13aiXi with a2=0 and Ad(exp(ϵXi))=i=13aiXi with a¯2=0 for i=1,2,3 .
Proof. The adjoint action Ad(exp(ϵX1)) do not the change sgn of a3 , i.e. the coefficient of X3 . □
Lemma 2.6. Let
K={sgn(a1),a10,a2=00otherwise.
where K is a sgn function.
Then K(X)=K(Ad(exp(ϵX))), for X with a2=0 and Ad(exp(ϵXi))=i=13aiXi with a¯2=0 and i=1,2,3 .
Proof. By Lemma 2.5, that can be proved. □
Construction of adjoint matrix:
Apply the adjoint action of X1 to X,
Ad(exp(ϵX1))X=a1X1+a2X2+a3X3=(a1+a2+a3)A1(X1+X2+X3)T,
where
A1=(100ϵ110001).
The matrices A2 and A3 can be calculated as:
A2=(eϵ20001000e3ϵ2),A3=(100013ϵ3001)
and
A=A1A2A3=(eϵ200ϵ1eϵ213ϵ300e3ϵ2)
A is termed as global adjoint matrix.
Lemma 2.7. The one-dimensional optimal system of Lie algebra L3 for CDSWEs (1) is P1=X2 , and P2=X1+a3X3 where, a3{1,1,0} .
Proof. The adjoint transformation can be calculated as
(a¯1,a¯2,a¯3)=(a1,a2,a3)A
It recasts as
a¯1=a1eϵ2a2ϵ1eϵ2,a¯2=a2,anda¯3=3ϵ3a2+a3e3ϵ2.
Further, solutions to (16) for ϵ1 , ϵ2 , and ϵ3 can be obtained as:
Step (I): If a2=1 , taking ϵ2=0 , for a¯1=0 , a¯2=1 , and a¯3=0 , then the system has a solution ϵ1=a1 , and ϵ3=a33 , thus the representative element is P1=X2 .
Step (II): If a2=0 , and a10 , then the authors took a1=1 , and ϵ2=0 , then system (16) follows three representative as X1+a3X3 , where a3{1,1,0} . i.e. P2=X1+X3 , P3=X1X3 , and P4=X1 .

3. Similarity solutions

This section shows how to derive invariant solutions and similarity reductions for different classes of optimum sub algebra. All of the cases can be categorized as:
Case (1) For P1=X2=xx+3tt2uu2vv , the characteristics equations can be written as
dxx=dt3t=du2u=dv2v,

which yields the similarity variable X=xt13 and similarity functions as u=F(X)t23 and v=G(X)t23 . Thus CDSWEs (1) reduces into an equivalent form
XF3aGG+2F=0,and
3bG+3cFG+3dFGXG2G=0.

Prime shows the derivative of a function with respect to the indicated variable.
Multiply Eq. (18) by G(X) and by F(X) in Eq. (19), then adding them we have an ODE. In which taking F=G , and after twice integration, one can get
b2(F)2+d+ca6F3=C1F+C2,

where C1 and C2 are constants of integration.
If C2=0 , then solution for Eq. (1) can be explored as
u1(x,t)=v1(x,t)=Ft23,

where F can be expressed as:
EllipticF(1+Fα,12)=(xt13+α+C3)C1/b,
where α=6C1d+ca and C3 is a constant of integration.
Another particular solution of (18) and (19) can be obtained treating F and G as
F=Xc+d and G=±1a(c+d)X ,

then solution of the system (1) is
u2(x,t)=x(c+d)t,andv2(x,t)=1a(c+d)xt.
Case (2) For P2=X1±X3=x±t , the corresponding Lagrange's characteristic equation is
dx=±dt=du0=dv0.

It yields X1=xt and u=F(X1) , v=G(X1). Then, similarity reduction of the CDSWEs (1) is
aGGF=0,
bG+cFG+dFGG=0.

The authors have used double signs (± or ) wherever in the whole manuscript, which are considered upper sign with upper and lower sign with lower.
Taking Eq. (24) and integrating it yields
F=±a2G2+C4,

where C4 is a constant of integration.
From Eqs. (25), (26), getting
bG±α12G2G+α2G=0,

where α1=ac+2ad and α2=cC41 are constants and twice integrating yields
(G)2±α112bG4+α2bG2C5G=C6,

where C5 and C6 are constants of integration.
Case (2a) For α1=±3k3 , b=k , α2=14k3 , and C5=C6=0 , the solution for the CDSWEs (1) is
u3(x,t)=±a2sech2[k(xt)2+C7]+C4,andv3(x,t)=sech[k(xt)2+C7].
Case (2b) Assuming α1=3k5 , b=k , α2=12k3 , C5=0 , and C6=14 , the solution for the system (1) can be obtained as
u4(x,t)=±a2k2tan2[k(xt)2+C8k]+C4,andv4(x,t)=1ktan[k(xt)2+C8k].
Case (2c) For α1=3k5 , b=k , α2=14k5 , and C5=C6=0 , the solution for the system (1) are
u5(x,t)=±a2sec2[k2(xt)2+C9]+C4,andv5(x,t)=sec[k2(xt)2+C9];
u6(x,t)=±a2csc2[k2(xt)2+C10]+C4,andv6(x,t)=csc[k2(xt)2+C10].
Case (2d) For α1=3k3 , b=k , and α2=C5=C6=0 , the solution is
u7(x,t)=±a2[2±k(xt)+C11]2+C4,andv7(x,t)=2±k(xt)+C11.
Case (2e) Putting α1=3k5 , b=k , α2=12k3,C5=0 , and C6=14 , the solutions can be expressed as
u8(x,t)=±a2k2tanh2[±k(xt)2+C12]+C4,andv8(x,t)=±1ktanh[±k(xt)2+C12];
u9(x,t)=±a2[21exp(k(xt)+C13)1]2+C4,andv9(x,t)=1k[21exp(k(xt)+C13)1].
Case (2f) Assuming α1=3k3 , b=k , α2=±14k3 , and C5=C6=0 , the solution of CDSWEs is
u10(x,t)=±a2csch2[k(xt)2+C14]+C4,andv10(x,t)=csch[k(xt)2+C14].
Case (2g) For c=2k1 , d=k1 , α2=k2 , C5=0 , C6=k2b , and b0 , the solution of system (1) is
u11(x,t)=±a2sin2[k(xt)b+C15]+C4,andv11(x,t)=sin[k(xt)b+C15];
u12(x,t)=±a2cos2[k(xt)b+C16]+C4,andv12(x,t)=cos[k(xt)b+C16].
Case (2h) Assuming c=2k1 , d=k1 , α2=k2 , C5=0 , C6=k2b , and b0 , the solution of Eq. (1) is
u13(x,t)=±a2sinh2[k(xt)b+C17]+C4,andv13(x,t)=sinh[k(xt)b+C17].
Case (2i) Taking c=2k1 , d=k1 , α2=k2 , C5=0 , C6=k2b , and b0 , the solution of CDSWEs is
u14(x,t)=±a2cosh2[k(xt)b+C18]+C4,andv14(x,t)=cosh[k(xt)b+C18].
Case (3) : For P4=X1=x , the characteristic equation is
dx=dt0=du0=dv0.
The trivial solution for system (1) is
u15(x,t)=v15(x,t)=C19.

4. Comparison with results reported earlier

The novelty of the results can be confirmed as some of the reported results [26], [28], [29], [30], [33], [38], [42], [49] are derived from the solutions attained in this article. Other solutions are absolutely different from earlier findings as far as authors are aware.
1. If C4=C7=0 is substituted in the expressions of u3 , and v3 , then Eqs. (12) and (13) of Arora and Kumar [28], Eq. (52) of Yao [33], and Eqs. (20) and (21) of Zhang [38] can be retrieved.
2. Taking C4=0 in u3,v3 , u10,v10 , u5,v5 , u6,v6 , u8,v8 , u4,v4 , so the results represented by Eqs. (35)-(48) of Arnous et al. [29] can be obtained respectively.
3. Eqs. (3.21) and (3.23) of Khan et al. [26], Eq. (52) of Yao [33], and Eq. (37) of Shi et al. [42] can be obtained by replacing C4=C12=0 in u8 , and v8 of this study.
4. Putting C4=C12=0 in u8 and v8 of this study, one can get Eqs. (3.21) and (3.23) of Khan et al. [26], Eq. (52) of Yao [33], and Eq. (37) of Shi et al. [42].
5. For C4=C11=0 in u7 and v7 , then u2(ξ) and v2(ξ) in Eq. (36) of Zhang and Zhao [49] can be obtained.
6. Taking C4=C8=0 in u4 , and v4 , then Eq. (36) u3(ξ) , v3(ξ) of Zhang and Zhao [49], and u12 , v23,24 of Akbar and Ali [30] can be deduced.
7. Eqs. (28) and (29) of Zhang [38] can be proved replacing C4=C9=0 in u5 and v5 respectively.

5. Physical analysis of solutions and discussions

The graphical illustration of the mathematical expressions (21), (22), (29)-(40) and (42) of invariant solutions might make them more significant. These solutions are different from those available in the existing literature [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50]. Using the MATLAB code in the space range 20x20 , and selecting appropriate values of arbitrary variables and parameters such as a=0.9134 , k=0.1419 , and C1=0.4218 . Variations in u , and v can be observed via dominant behaviour frames of animations created with respect to changes in space and time.

6. Conclusions

Lie symmetry analysis is exploited successfully to produce one-dimensional optimum sub algebras of CDSWEs (1) in this research. The killing form has been used to derive invariants. The analytical solutions of CDSWEs (1) are represented by Eqs. (21), (22), (29)-(40) and (42). Some of the existing results [26], [28], [29], [30], [33], [38], [42], [49] can be derived from our solutions, demonstrating the originality of the results. The remaining solutions in this study are completely different from those established in the published works in Singh et al. [25], Khan et al. [26], Lu et al. [27], Arora and Kumar [28], Arnous et al. [29], Akbar and Ali [30], Chen et al. [31], Tariq and Seadawy [32], Yao [33], Lu et al. [34], Wen et al. [35], Javeed et al. [36], Shehata et al. [37], Zhang [38], Bibi and Tauseef [39], Bhatter et al. [40], Matjila et al. [41], Shi et al. [42], Khan et al. [43], Akbar et al. [44], Zhao and Zhi [45], Ullah et al. [46], Inc [47], Morris and Kara [48], Zhang and Zhao [49], Zhao et al. [50]. Asymptotic, bell-shaped, bright and dark soliton, bright soliton, parabolic, bright and kink, kink, and periodic nature of the solutions are shown in Fig. 1, Fig. 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7, Fig. 8, Fig. 9, Fig. 10, Fig. 11, Fig. 12, Fig. 13. Over a horizontal seafloor, CDSWEs (1) depict the dispersion of nonlinear surface gravity waves whose amplitude varies with time and space. Coastal engineers may get benefit from the analytical solutions of CDSWEs (1) with physical nature found in this study.
Fig.1 Asymptotic nature of u2 and v2 with c=0.6324 and d=0.0975 .
Fig.2 Bell-shaped profiles of the expressions u3 and v3 with constant C7=0.9157 .
Fig.3 Bright and dark solitons for u4 and v4 at arbitrary constant C8=0.9595 .
Fig.4 Parabolic profiles of the u5 and v5 taking C9=0.1712 .
Fig.5 Bright and kink profiles of the solutions represented as u6 and v6 with choice of C10=0.0971 .
Fig.6 Bright soliton for u7 and v7 with C11=0.8235 .
Fig.7 Kink profile of the u8 and v8 with C12=0.8147 .
Fig.8 Bright and kink profiles for u9 and v9 with C13=0.1270 .
Fig.9 Bright and kink soliton profile for u10 and v10 at C14=0.2785 .
Fig.10 Periodic nature of u11 and v11 choosing C15=0.9134 .
Fig.11 Periodic of u12 and v12 with C16=0.7431 .
Fig.12 Asymptotic nature profiles of u13 and v13 taking C17=0.7577 .
Fig.13 Asymptotic of u14 and v14 treating C18=0.6787 .

Availability of data and materials

Data sharing does not apply to this article.

Funding

Not applicable.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

[1]
M. Boiti, J.J.-P. Leon, F. Pempinelli. Inverse Probl., 3 (1987), pp. 371-387, DOI: 10.1088/0266-5611/3/3/007.

[2]
T.H.C. Herbers, S. Elgar, N.A. Sarap, R.T. Guza. J. Phys. Oceanogr., 32 (2002), pp. 1181-1193

[3]
R. Hirota, B. Grammaticos, A. Ramani. J. Math. Phys., 27 (1986), pp. 1498-1505, DOI: 10.1063/1.527110.

[4]
V.G. Drinfeld, V.V. Sokolov. Sov. Math. Dokl., 23 (1981), pp. 457-462

[5]
R. Hirota, J. Satsuma. Phys. Lett. A., 85 (8-9) (1981), pp. 407-408

[6]
J. Satsuma, R. Hirota. J. Phys. Soc. Jpn., 51 (1982), pp. 3390-3397

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

[8]
M. Niwas, S. Kumar, H. Kharbanda. J. Ocean Eng. Sci., 7 (2022), pp. 188-201

[9]
M.R. Ali, W.-X. Ma, R. Sadat. J. Ocean Eng. Sci. (2021), DOI: 10.1016/j.joes.2021.08.006.

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

[11]
A. Akbulut, M. Kaplan, M.K.A. Kaabar. J. Ocean Eng. Sci.(2021), DOI: 10.1016/j.joes.2021.09.010.

[12]
A.M. Wazwaz. Ocean Eng., 60 (2013), pp. 95-98, DOI: 10.1016/j.oceaneng.2012.12.049.

[13]
S. Kumar, M. Niwas, S.K. Dhiman. J. Ocean Eng. Sci.(2021), DOI: 10.1016/j.joes.2021.10.009.

[14]
A.C. Varsoliwala, T.R. Singh. J. Ocean Eng. Sci.(2021), DOI: 10.1016/j.joes.2021.03.003.

[15]
A.C. Varsoliwala, T.R. Singh. J. Ocean Eng. Sci., 6 (2021), pp. 367-375

[16]
S. Kumar, A. Kumar, B. Mohan. J. Ocean Eng. Sci. (2021), DOI: 10.1016/j.joes.2021.11.002.

[17]
S. Kumar, A. Kumar. J. Ocean Eng. Sci., 7 (2022), pp. 178-187

[18]
L.J.F. Broer. Appl. Sci. Res., 31 (1975), pp. 377-395, DOI: 10.1007/BF00418048.

[19]
K.U. Tariq, M. Inc, A. Pashrashid, M. Zubair, L. Akinyemi. J. Ocean Eng. Sci. (2022), DOI: 10.1016/j.joes.2022.01.011.

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

[21]
K. Debin, H. Rezazadeh, N. Ullah, J. Vahidi, K.U. Tariq, L. Akinyem. J. Ocean Eng. Sci. (2022), DOI: 10.1016/j.joes.2022.03.007.

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

[23]
N. Gedika, E. Irtem, S. Kabdasli. Ocean Eng., 32 (2005), pp. 513-528

[24]
M.M. Mousa. J. Ocean Eng. Sci., 3 (2018), pp. 303-309

[25]
P.K. Singh, K. Vishal, T. Som. Appl. Math., 10 (2015), pp. 460-472 DOI: 10.1504/IJICBM.2015.069647

[26]

[27]
D. Lu, D. Kang, B. Hong. J. Inf. Comput. Sci., 10 (18) (2013), pp. 5955-5962, DOI: 10.12733/jics20102632.

[28]
R. Arora, A. Kumar. Adv. Sci. Eng. Med., 5 (2013), pp. 1-7, DOI: 10.1166/asem.2013.1399.

[29]
A.H. Arnous, M. Mirzazadeh, M. Eslami. Pramana, 86 (2016), pp. 1153-1160, DOI: 10.1007/s12043-015-1179-1.

[30]
M.A. Akbar, N.M. Ali. Appl. Interdiscip. Math., 4 (2017), pp. 1282577-1282596, DOI: 10.1080/23311835.2017.1282577.

[31]
X. Chen, W. Shen, L. Wang, F. Wang. IAENG Int. J. Appl. Math.,47 (2) (2017), pp. 06-13 DOI: 10.37549/ar2385

[32]
K.U. Tariq, A.R. Seadawy. J. King Saud Univ. Sci., 4 (2018), pp. 1-21, DOI: 10.1016/j.jksus.2018.03.019.

[33]

[34]
D. Lu, A.R. Seadawy, M. Iqbal. Results Phys., 4 (2018), pp. 1-24, DOI: 10.1016/j.rinp.2018.11.014.

[35]

[36]
S. Javeed, S. Saif, D. Baleanu. Adv. Differ. Equ.,2018 (2018), pp. 459-464, DOI: 10.1186/s13662-018-1913-3.

[37]
A.R. Shehata, E.M. Kamal, H.A. Kareem. Int. J. Pure Appl. Math., 101 (4) (2015), pp. 477-487, DOI: 10.12732/ijpam.v101i4.3.

[38]
W.M. Zhang. Appl. Math. Sci., 5 (38) (2011), pp. 1887-1894

[39]

[40]
S. Bhatter, A. Mathur, D. Kumar, J. Singh. Phys. A., 537 (2020), pp. 122578-122591, DOI: 10.1016/j.physa.2019.122578.

[41]
C. Matjila, B. Muatjetjeja, C.M. Khalique. Abstr. Appl. Anal., 2014 (2014), pp. 6-13, DOI: 10.1155/2014/271960.

[42]
Y. Shi, X. Li, B.G. Zhang. Adv. Math. Phys., 2018 (2018), pp. 8-17, DOI: 10.1155/2018/8583418.

[43]
K. Khan, M.A. Akbar, H. Koppelaar. R. Soc. Open Sci., 2 (2015), pp. 140406-140419, DOI: 10.1098/rsos.140406.

[44]
M.A. Akbar, N.H.M. Ali, S.T. Mohyud-Din. Springerplus, 2 (2013), pp. 327-343

[45]
X.Q. Zhao, H.Y. Zhi. Commun. Theor. Phys., 50 (2008), pp. 309-314

[46]
H. Ullah, S. Islam, M. Idrees, R. Nawaz. Math. Probl. Eng., 6 (2013), pp. 8-17, DOI: 10.1155/2013/362816.

[47]

[48]
R. Morris, A.H. Kara. Appl. Math. Comput., 219 (2013), pp. 6473-6483

[49]
Y. Zhang, Z. Zhao. Bound. Value Probl., 154 (2017), pp. 154-171, DOI: 10.1186/s13661-017-0885-7.

[50]
Z. Zhao, Y. Zhang, Z. Han. Eur. Phys. J. Plus, 129 (2014), pp. 143-150

[51]
M. Kumar, R. Kumar. J. Phys. Conf. Ser., 574 (2015), pp. 012038-012043, DOI: 10.1088/1742-6596/574/1/012038.

[52]
S. Sahoo, G. Garai, S.S. Ray. Nonlinear Dyn., 87 (3) (2016), pp. 1995-2000, DOI: 10.1007/s11071-016-3169-3.

[53]
S. Sahoo, S.S. Ray. Comput. Math. Appl., 70 (2) (2015), pp. 158-166

[54]
M. Kumar, A.K. Tiwari, R. Kumar. Proc. Natl. Acad. Sci., India, Sect. A, 89 (4) (2018), pp. 747-755, DOI: 10.1007/s40010-018-0527-1.

[55]
R. Kumar. IOSR-JM, 12 (4) (2016), pp. 144-147, DOI: 10.9790/5728-120402144147.

[56]
R. Kumar, A. Kumar. Appl. Math. Comput., 422 (2022), p. 26976, ( 1-18)

[57]
R. Kumar, M. Kumar, A.K. Tiwari. Int. J. Comput. Methods Eng., 22 (3) (2021), pp. 225-234, DOI: 10.1080/15502287.2021.1916693.

[58]
R. Kumar, A. Kumar. Int. J. Appl. Comput. Math., 8 (11) (2022), pp. 1-22, DOI: 10.1007/s40819-021-01211-2.

[59]

[60]
R. Kumar. IJSET, 4 (2015), pp. 1430-1437

[61]
L. Akinyemi, K.S. Nisar, C. Ahamed Saleel, H. Rezazadeh, P. Veeresha, M.M.A. Khater, M. Inc. Results Phys., 31 (2021), p. 104958

[62]
M.M.A. Khater, M.A. Akbar, L. Akinyemi, A. Jhangeer, H. Rezazadeh, M. Inc. Opt. Quantum Electron., 53 (2021), p. 507 ( 1-16)

[63]
L. Akinyemi, H. Rezazadeh, Q.-H. Shi, M. Inc, M.M.A. Khater, H. Ahmad, A. Jhangeer, M.A. Akbar.Results Phys., 29 (2021), p. 104656

[64]
M.S. Mohamed, L. Akinyemi, S.A. Najati, S.K. Elagan. Opt.Quantum Electron., 54 (2022), p. 141

[65]
M. Mirzazadeh, A. Akbulut, F. TaÅŸcan, L. Akinyemi. Optik, 252 (2022), p. 168529

[66]
L. Akinyemi, M. Mirzazadeh, S.A. Badri, K. Hosseini. J. Mod. Opt., 69 (3) (2022), pp. 172-182

[67]
R. Kumar, R.S. Verma. Nonlinear Dyn. (2022), DOI: 10.1007/s11071-022-07389-4.

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

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

Outlines

/