Original article

An extensive analytical and numerical study of the generalized q -deformed Sinh-Gordon equation

  • Khalid K. Ali , a, * ,
  • Abdel-Haleem Abdel-Aty , b, c
Expand
  • a Mathematics Department, Faculty of Science, Al-Azhar University, Nasr-City, Cairo, Egypt
  • b Department of Physics, College of Sciences, University of Bisha, PO Box 344, Bisha 61922, Saudi Arabia
  • c Physics Department, Faculty of Science, Al-Azhar University, Assiut 71524, Egypt
*E-mail addresses: (K.K. Ali);
(A.-H. Abdel-Aty).

Received date: 2022-04-29

  Revised date: 2022-05-22

  Accepted date: 2022-05-29

  Online published: 2022-06-01

Abstract

In this paper, the q-deformed Sinh-Gordon equation is solved analytically using a new general form based on the extended tanh approach. The numerical solutions of the equation is obtained using a b-spline finite element method. Also, we present numerous figures to demonstrate the various solitons propagation patterns. This type of equation has not been previously dealt with in such ways, whether analytical or numerical. This study is very useful in studying several physical systems that have lost their symmetry.

Highlights

● The q-deformed Sinh-Gordon equation was investigated analytically using the new general form of the extended tanh approach and numerically using a b-spline finite element method in this paper.

● We also present numerous figures to demonstrate the various solitons propagation patterns.

● This type of equation has not been previously dealt with in such ways, whether analytical or numerical.

● We believe that this study will be very useful in studying physical systems that have lost their symmetry.

Cite this article

Khalid K. Ali , Abdel-Haleem Abdel-Aty . An extensive analytical and numerical study of the generalized q -deformed Sinh-Gordon equation[J]. Journal of Ocean Engineering and Science, 2024 , 9(3) : 232 -239 . DOI: 10.1016/j.joes.2022.05.034

1. Introduction

Nonlinear systems attracted much interest due to it's importance in explaining the dynamics of many systems in natures and science such as, gravitational waves, acoustics plasma waves, shallow water dynamics, no-linear optics, fluid dynamics and surface ocean waves [1], [2], [3], [4], [5], [6]. The microscopic systems that are governed by the principles of quantum physics is considered as a best example of q-deformed Sinh-Gordon equation [7], [8]. Many mathematical models describing some physical systems have been extensively studied by a group of researchers [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24].
Quantum calculus, often known as calculus without bounds, is the same as regular infinitesimal calculus but without the concept of bounds. It defines ”q -calculus” and ”h -calculus”, where h ostensibly stands for Planck's constant while q stands for quantum. The two parameters are related by the formula [25]
q=eih=e2iπ,
where =h2π . The study of the equations containing q -calculus has not been covered except in limited papers such as [26], [27].
q -calculus, this area of research has several applications there are several developments and applications of the q -calculus in mathematical physics, the theory of relativity and special functions [28], [29], [30], [31], [32], [33].
The study of the generalized q -deformed Sinh-Gordon equation (Eleuch equation) is the subject of this article [26], [27]. Eleuch equation is written as:
2Ωx22Ωt2=[sinhq(βΩα)]mϑ.
where α,m are constantes and q(0,1) . For q=1 gives the standard Sinh-Gordon equation.
In this work, we use the extended tanh approach to examine the soliton solutions of (1.2) [34], [35]. Besides, we solve this equation numerically using the b-spline finite element method [36], [37], [38], [39]. We study two cases for this equation.
This paper has been organized as follows: The second section presents some definitions and properties of q -calculus. In the third section, we present the methodology of the analytical method. In the fourth section, we present the solutions obtained. In the fifth section, we present the numerical solutions to the equation. In the sixth section, we present some different forms of solutions. In the seventh section, we present a conclusion of what we have done in the work.

2. q-Calculus

Now, we introduce some basic definitions of the q -calculus as follows: [33]
● For any real number a , let us define [a]q by [a]q=1qa1q,aR.
● In the q -calculus, differential of function is defined as:
dq(f(t))=f(qt)f(t).
● The q -gamma function is defined by
Γq(t)=G(qt)(1q)t1G(q),tR{0,1,2,},
where G(qt)=1(qt;q). Or, equivalently, Γq(t)=1q(t1)1qt1, and satisfies Γq(t)+1)=[t]qΓq(t).
● The q -derivative of a function f is given by (Dqf)(t)=dq(f(t))dq(t)=f(t)f(tq)tqt,(Dqf)(0)=limt(Dqf)(t).
● The sinhq is defined by: sinhq(t)=etqet2,
● The coshq is defined by: coshq(t)=et+qet2,
● We list below some simple and useful relations for q -deformed functions
(coshq(t))2(sinhq(t))2=q,
sinhqt)=qsinh1q(t), coshqt)=qcosh1q(t),
ddtcoshq(t)=sinhq(t), ddtsinhq(t)=coshq(t),

3. Mathematical analysis of the model

The following transformation is used to determine the traveling wave solution of (1.2),
Ω(x,t)=u(δ),
where
δ=axct+θ.
where c denotes the speed of traveling wave. Using (3.1) and (3.2), then (1.2) becomes,
(a2c2)u(δ)+ϑ[sinhqu(δ)α)]m=0.
Now, we take two cases for (3.3).
Case one: Suppose m=α=1,ϑ=0.
Then, (3.3) can be written as:
(a2c2)u(δ)sinhq(u(δ))=0.
We can multiply both side of (3.4) by u(δ) and after the integration we get
12(a2c2)(u(δ))2coshq(u(δ))c1=0,
where c1 is a constant of the integration.
Let
u(δ)=ln(v(δ)).
Thus, (3.5) becomes,
(c2a2)v2(δ)+2C1v2(δ)+qv(δ)+v3(δ)=0.
Thus, we can solve (3.7) using our method and from (3.6) and (3.1) we can get the solution of (1.2) in the first case.
Case two: Suppose α=1,m=2,ϑ=q2.
Then, (3.3) can be written as:
(a2c2)u(δ)(sinhq(u(δ)))2+q2=0.
After simplify (3.8) we get
(a2c2)u(δ)12coshq2(2u(δ))=0.
Let
u(δ)=12ln(v(δ)).
Thus, (3.10) becomes,
2(a2c2)v2(δ)2(a2c2)v(δ)v(δ)+q2v(δ)+v3(δ)=0.
Thus, we can solve (3.11) using our method and from (3.10) and (3.1) we can get the solution of (1.2) in the second case.

4. The strategy of the extended tanh function method

Postulate the governing equation in the form
F(Ω,Ωxx,Ωtt,)=0,
where: F is a polynomial Ω=Ω(x,t) and its partial derivatives. Applying a traveling wave transformations to convert (3.1) to an ordinary differential equation:
H(u,u,)=0.
The essential steps of the new generalization for the Extended tanh function method are:
Step 1: The modified extended tanh method present the solution of (4.2) by the following finite series:
u(δ)=A0+i=1N(Aiϕ(δ)i+Biϕ(δ)i),
where ϕ=ϕ(δ) , satisfy the Riccati equation
dϕdδ=ϱ+ϕ2.
The solution of the Riccati Eq. (4.4) has the following cases of solutions: If ϱ<0 then
ϕ(δ)=ϱtanh(ϱδ),ϕ(δ)=ϱcoth(ϱδ).
If ϱ=0 then
ϕ(δ)=1δ.
If ϱ>0 then
ϕ(δ)=ϱtan(ϱδ),ϕ(δ)=ϱcot(ϱδ).
Step 2: N can be determined by balancing the highest order derivative term with the highest power nonlinear term in (4.2).
Step 3: Substituting (4.3) and (4.4) into (4.2) then gather all coefficients of the same powers of ϕ(δ)i and put them equal to zero, we get a system of algebraic equations for ϱ,A0,,AN,B1,,BN , solving this equations we get all constants.

5. The analytical solution for the model

In this section, the generalization of kudryashov method applies to find the analytic solution to the two cases that were imposed for the Eq. (1.2)
The analytical solution of case one at m=α=1,ϑ=0 :
Applying the balance principle in (3.7) between v2 and v3 we get 2N+2=3NN=2 . From (4.3), the solution of (3.7) can be presented as:
v(δ)=A0+i=12(Aiϕ(δ)i+Biϕ(δ)i).
Substituting (5.1) into (3.7), setting the coefficient of like power of ϕ(δ) equal to zero, we acquire the following system: 4a2A1B1ϱ+16a2A2B2ϱa2A12ϱ2a2B124A1B1c2ϱ16A2B2c2ϱ +4A1B1c1+4A2B2c1+6A1A0B1+6A2A0B2+3A2B12+3A12B2+A12c2ϱ2 +2A02c1+A0q+A03+B12c2=0,
8a2A2B1ϱ+4a2A1B24a2A1A2ϱ28A2B1c2ϱ4A1B2c2+4A2B1c1+3A12B1+6A0A2B1+6A1A2B2+4A1A2c2ϱ2+4A0A1c1 +A1+q+3A02A1=0, 2a2A1B1+8a2A2B24a2A22ϱ22a2A12ϱ2A1B1c28A2B2c2 +6A1A2B1+3A22B2+4A22c2ϱ2+2A12c2ϱ+2A12c1 +4A0A2c1+A2q+3A0A12+3A02A2=0,
4a2A2B18a2A1A2ϱ4A2B1c2+3A22B1+8A1A2c2ϱ +4A1A2c1+A13+6A0A1A2=0,
4a2A1A2+4A1A2c2+3A1A22=0, 8a2A22ϱa2A12+8A22c2ϱ+A12c2+2A22c1++3A2A12+3A0A22=0, 4a2A22+4A22c2+A23=0,
4a2A2B1ϱ2+8a2A1B2ϱ4a2B1B24A2B1c2ϱ28A1B2c2ϱ +4A0B1c1+4A1B2c1+3A1B12+3A02B1+6A0A1B2+6A2B1B2 +4B1B2c2+B1+q=0, +2a2A1B1ϱ2+8a2A2B2ϱ22a2B12ϱ4a2B222A1B1c2ϱ2 8A2B2c2ϱ2+4A0B2c1+3A0B12+3A2B22+3A02B2 +6A1B1B2+2B12c2ϱ+4B22c2+2B12c1+B2q=0, a2B12ϱ28a2B22ϱ+3A0B22+B12c2ϱ2+8B22c2ϱ+2B22c1+3B2B12=0, 4a2A1B2ϱ28a2B1B2ϱ4A1B2c2ϱ2 +3A1B22+6A0B1B2+8B1B2c2ϱ+4B1B2c1+B13=0, 4a2B1B2ϱ2+4B1B2c2ϱ2+3B1B22=0, 4a2B22ϱ2+4B22c2ϱ2+B23=0. Solving the previous set of equations with the aid of Mathematica program, we obtain the following set of solutions.
Class 1:
A0=q2,A1=0,A2=q4ϱ,B1=0,B2=qϱ4,a=c2q16ϱ,c1=q.
Substituting from (5.2) into (5.1) with (4.5), (4.6), (4.7), (3.6) and (3.1), we get the solutions of (1.2) at m=α=1,ϑ=0 :
If ϱ<0 , then
Ω1,2(x,t)=ln(A2(ϱcoth((axct+θ)ϱ))2+A0+B2(ϱcoth((axct+θ)ϱ))2).
Ω3,4(x,t)=ln(A2(ϱtanh((axct+θ)ϱ))2+A0+B2(ϱtanh((axct+θ)ϱ))2).
If ϱ>0 , then
Ω5,6(x,t)=ln(A2(ϱcot((axct+θ)ϱ))2+A0+B2(ϱcot((axct+θ)ϱ))2).
Ω7,8(x,t)=ln(A2(ϱtan((axct+θ)ϱ))2+A0+B2(ϱtan((axct+θ)ϱ))2).
Class 2:
A0=0,A1=0,A2=qϱ,B1=0,B2=0,a=c2q4ϱ,c1=q.
Substituting from (5.7) into (5.1) with (4.5), (4.6), (4.7), (3.6) and (3.1), we get the solutions of (1.2) at m=α=1,ϑ=0 :
If ϱ<0 , then
Ω9,10(x,t)=ln(A2(ϱcoth((axct+θ)ϱ))2).
Ω11,12(x,t)=ln(A2(ϱtanh((axct+θ)ϱ))2).
If ϱ>0 , then
Ω13,14(x,t)=ln(A2(ϱcot((axct+θ)ϱ))2).
Ω15,16(x,t)=ln(A2(ϱtan((axct+θ)ϱ))2).
Class 3:
A0=0,A1=0,A2=0,B1=0,B2=qϱ,a=c2+q4ϱ,c1=q.
Substituting from (5.12) into (5.1) with (4.5), (4.6), (4.7), (3.6) and (3.1), we get the solutions of (1.2) at m=α=1,ϑ=0 :
If ϱ<0 , then
Ω9,10(x,t)=ln(B2(ϱcoth((axct+θ)ϱ))2).
Ω11,12(x,t)=ln(B2(ϱtanh((axct+θ)ϱ))2).
If ϱ>0 , then
Ω13,14(x,t)=ln(B2(ϱcot((axct+θ)ϱ))2).
Ω15,16(x,t)=ln(B2(ϱtan((axct+θ)ϱ))2).
The analytical solution of case two at m=2,α=1,ϑ=q2 :
Applying the balance principle in (3.11) between vv and v3 we get 2N+2=3NN=2 . From (4.3), the solution of (3.7) can be presented as:
v(δ)=A0+i=12(Aiϕ(δ)i+Biϕ(δ)i).
Substituting (5.17) into (3.11), setting the coefficient of like power of ϕ(δ) equal to zero, we acquire the following system:
16a2A1B1ϱ64a2A2B2ϱ4a2A0B24a2A2A0ϱ2+2a2A12ϱ2 +2a2B12+16A1B1c2ϱ+64A2B2c2ϱ+4A0B2c2+6A1A0B1 +6A2A0B2+3A2B12+3A12B2+4A2A0c2ϱ22A12c2ϱ2 +A0q2+A032B12c2=0,
36a2A2B1ϱ16a2A1B2+4a2A1A2ϱ24a2A0A1ϱ +36A2B1c2ϱ+16A1B2c2+3A12B1+6A0A2B1+6A1A2B2 4A1A2c2ϱ2+4A0A1c2ϱ+A1+q2+3A02A1=0, 8a2A1B132a2A2B2+4a2A22ϱ216a2A0A2ϱ+8A1B1c2 +32A2B2c2+6A1A2B1+3A22B24A22c2ϱ2 +16A0A2c2ϱ+A2q2++3A0A12+3A02A2=0,
20a2A2B14a2A2A1ϱ4a2A0A1+20A2B1c2+3A22B1 +4A2A1c2ϱ+4A0A1c2+A13+6A0A2A1=0,
2a2A1212a2A0A2+2A12c2 +12A0A2c2+3A0A22+3A12A2=0,
8a2A1A2+8A1A2c2+3A1A22=0, 4a2A22+4A22c2+A23=0,
16a2A2B1ϱ24a2A0B1ϱ36a2A1B2ϱ+4a2B1B2 +16A2B1c2ϱ2+4A0B1c2ϱ+36A1B2c2ϱ+3A1B12
+3A02B1+6A0A1B2+6A2B1B24B1B2c2+B1q2=0, 8a2A1B1ϱ232a2A2B2ϱ216a2A0B2ϱ+4a2B22 +8A1B1c2ϱ2+32A2B2c2ϱ2+16A0B2c2ϱ+3A0B12
+3A2B22+3A02B2+6A1B1B24B22c2+B2q2=0,
4a2A0B1ϱ220a2A1B2ϱ24a2B2B1ϱ+4A0B1c2ϱ2
+20A1B2c2ϱ2+6A0B2B1+3A1B22+4B2B1c2ϱ+B13=0,
12a2A0B2ϱ22a2B12ϱ2+12A0B2c2ϱ2 +3A0B22+2B12c2ϱ2+3B12B2=0, 8a2B1B2ϱ2+8B1B2c2ϱ2+3B1B22=0,
4a2B22ϱ2+4B22c2ϱ2+B23=0.
Solving the previous set of equations with the aid of Mathematica program, we obtain the following set of solutions.
Class 1:
A0=12(iq),A1=0,A2=iq4ϱ,B1=0,B2=iqϱ4,a=16c2ϱ+iq4ϱ.
Substituting from (5.18) into (5.17) with (4.5), (4.6), (4.7), (3.6) and (3.1), we get the solutions of (1.2) at m=α=1,ϑ=0 :
If ϱ<0 , then
Ω1,2(x,t)=ln(A2(ϱcoth((axct+θ)ϱ))2+A0+B2(ϱcoth((axct+θ)ϱ))2).
Ω3,4(x,t)=ln(A2(ϱtanh((axct+θ)ϱ))2+A0+B2(ϱtanh((axct+θ)ϱ))2).
If ϱ>0 , then
Ω5,6(x,t)=ln(A2(ϱcot((axct+θ)ϱ))2+A0+B2(ϱcot((axct+θ)ϱ))2).
Ω7,8(x,t)=ln(A2(ϱtan((axct+θ)ϱ))2+A0+B2(ϱtan((axct+θ)ϱ))2).
Class 2:
A0=0,A1=0,A2=iqϱ,B1=0,B2=0,a=4c2ϱ+iq2ϱ.
Substituting from (5.23) into (5.17) with (4.5), (4.6), (4.7), (3.6) and (3.1), we get the solutions of (1.2) at m=α=1,ϑ=0 : If ϱ<0 , then
Ω9,10(x,t)=ln(A2(ϱcoth((axct+θ)ϱ))2).
Ω11,12(x,t)=ln(A2(ϱtanh((axct+θ)ϱ))2).
If ϱ>0 , then
Ω13,14(x,t)=ln(A2(ϱcot((axct+θ)ϱ))2).
Ω15,16(x,t)=ln(A2(ϱtan((axct+θ)ϱ))2).
Class 3:
A0=0,A1=0,A2=0,B1=0,B2=iqϱ,a=4c2ϱ+iq2ϱ.
Substituting from (5.28) into (5.17) with (4.5), (4.6), (4.7), (3.6) and (3.1), we get the solutions of (1.2) at m=α=1,ϑ=0 : If ϱ<0 , then
Ω9,10(x,t)=ln(B2(ϱcoth((axct+θ)ϱ))2).
Ω11,12(x,t)=ln(B2(ϱtanh((axct+θ)ϱ))2).
If ϱ>0 , then
Ω13,14(x,t)=ln(B2(ϱcot((axct+θ)ϱ))2).
Ω15,16(x,t)=ln(B2(ϱtan((axct+θ)ϱ))2).

6. The numerical solution for the model

In this section, we shall verify the results obtained in the last section using the cubic B-spline method. First, we approximate the variables of the space an time which are x and t with their derivatives as in [36], [37], [38], [39]. Next, assuming that the value of the function Ω(x,t) which is the exact solution of the model at the points on the grid (xi,tj) and Ωi,j to be the same as the approximate solution at these points. The required values of Ωi,j and its first and the second derivatives, Ωi,j and Ωi,j , at nodal points xi are identified in terms of ci as
Ωi,j=ci+1,j+4ci,j+ci1,j,Ωx=Ωi,j=3h(ci+1,jci1,j),Ωxx=Ωi,j=6h2(ci1,j+ci+1,j2ci,j),
and if the time derivative is discretized using finite differences, we have where
Ωtt=ci,j1+ci,j+12ci,jk2.
Substituting (6.1) and (6.2) into (1.2) a we obtain the system of difference equations that can be solved to obtain the numerical results.

6.1. The numerical results

Now, we will introduce some numerical results for the generalized q -deformed Sinh-Gordon equation by using the cubic b-spline method. As we studied the analytic solution for two special cases of the generalized q -deformed Sinh-Gordon equation, we also present the numerical solution for these two cases they are as follows:
The numerical results for case one: m=α=1,ϑ=0
In Table 1 and Fig. 5 we introduce comparison between the numerical results with the analytical solution (5.5) for (1.2) at h=0.2,k=0.1,θ=1,q=0.01,c=0.2,ϱ=0.1 .
Table 1 Comparison between the numerical results with the analytical solution.
x Numerical solutionAnalytical solutionAbsolute error
-10.02.215272.215277.62745 E-9
-8.02.231272.231279.27982 E-9
-6.02.314822.314822.74825 E-8
-4.03.173303.173301.03583 E-8
4.02.535822.535826.64020 E-8
6.02.314822.314822.74823 E-8
8.02.231272.231279.28340 E-9
10.02.215272.215277.62753 E-9
Looking at the numerical results in the Table 1, which we obtained through the application of the cubic b-spline method, and comparing these results with the analytical solutions that we also obtained through the application of the extended tanh approach. We can see the extent to which the analytical results are in agreement with the numerical results. Which gives a clear indication of the two methods used in this work.
The numerical results for case two: m=2,α=1,ϑ=q2
In Table 2 and Fig. 6 we introduce comparison between the numerical results with the analytical solution (5.21) for (1.2) at h=0.2,k=0.1,θ=1,q=0.4,c=0.8,ϱ=0.01 .
Table 2 Comparison between the numerical results with the analytical solution.
x Numerical solutionAnalytical solutionAbsolute error
-10.00.9195510.9195511.11712 E-7
-8.00.9910310.9910312.30162 E-7
-6.01.0598901.0598906.37558 E-7
-4.00.9352640.9352651.78687 E-6
4.00.6022470.6022471.06871 E-7
6.00.9280880.9280886.04087 E-7
8.00.9829970.9829972.46975 E-7
10.00.9442640.9442641.10615 E-7
From taking a look at the numerical findings in the Table 2, which were produced using the cubic b-spline method, and compare them to the analytical results found using the extended tanh approach. The degree to which the analytical and numerical results are in agreement can be shown.
Now, after what we presented of analytical solutions and numerical solutions to two special cases of the proposed equation, we present in the next section some figures through which we explain our good results.

7. Graphical illustrations

Herein, we present some figures in the two-dimensional, three-dimensional to clarify the solutions that we presented. Some of the analytical and numerical solutions are presented in Fig. 1, Fig. 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6. In Fig. 1, we introduce the graph of (5.5) using our method at θ=1,q=0.1,c=0.2,ϱ=0.1 . Also, the graph of (5.6) θ=1,q=0.001,c=0.8,ϱ=0.001 , is presented in Fig. 2. Moreover, we present the graph of (5.20) at θ=1,q=0.5,c=0.5,ϱ=0.04 in Fig. 3. The graph of (5.21) θ=1,q=0.4,c=0.8,ϱ=0.01 , is presented in Fig. 4. Furthermore, in Fig. 5, we present comparison between the numerical results of (1.2) with the analytical solution (5.5) at h=0.2,k=0.1,θ=1,q=0.01,c=0.2,ϱ=0.1 . Finally, the comparison between the numerical results of (1.2) with the analytical solution (5.21) at h=0.2,k=0.1,θ=1,q=0.4,c=0.8,ϱ=0.01 , is presented in Fig. 5.
Fig.1 Graph of (5.5) at θ=1,q=0.1,c=0.2,ϱ=0.1 .
Fig.2 Graph of (5.6) at θ=1,q=0.001,c=0.8,ϱ=0.001 .
Fig.3 Graph of (5.20)using The new Kudryashov's method at θ=1,q=0.5,c=0.5,ϱ=0.04 .
Fig.4 Graph of (5.21)using The new Kudryashov's method at θ=1,q=0.4,c=0.8,ϱ=0.01 .
Fig.5 Comparison between the numerical results of (1.2) with the analytical solution (5.5) at h=0.2,k=0.1,θ=1,q=0.01,c=0.2,ϱ=0.1 .
Fig.6 Comparison between the numerical results of (1.2) with the analytical solution (5.21) at h=0.2,k=0.1,θ=1,q=0.4,c=0.8,ϱ=0.01 .

8. Conclusion

Finally, we presented an extended tanh approach method to study the generalized q -deformed Sinh-Gordon equation. We have studied the model in this way and presented figures showing the correctness of what we hoped to reach from the proposed method. In addition, we presented an extensive numerical study of this model using the finite element method. In addition, the analytical solutions and the numerical solutions were compared. Through what we have reached, we can say that our results are a clear contribution to this field. At the end of our article, we present a future plan for our study, as we can solve more than one model in this way, and we can also develop other ways to solve different models. The proposed equation has opened up new possibilities for modeling physical systems in which symmetry has been compromised.

Declaration of Competing Interest

There is no conflict between the authors.

[1]
G. Vilasi,J. Phys., 87 (2007), p. 012017 DOI: 10.1088/1742-6596/87/1/012017

[2]
S.D. Odintsov, V.K. Oikonomou,R. Myrzakulov. Symmetry, 14 (2022), p. 729 DOI: 10.3390/sym14040729

[3]
A. Goswami, J. Singh, D. Kumar, S. Gupta, Sushila. J. Ocean Eng. Sci., 4 (2) (2019), pp. 85-99

[4]
R.M. El-Shiekh, M. Gaballah. J. Ocean Eng. Sci., 5 (2) (2020), pp. 180-185

[5]
T. Islam, A. Akbar, H. Rezazadeh, A. Bekir. J. Ocean Eng. Sci. (2022), ISSN 2468-0133.

[6]
M. Cai, D. Li. Nonlinear Anal., 74 (13) (2011), pp. 4581-4588

[7]
A. Yusuf, T.A. Sulaiman, M. Mirzazadeh, K. Hosseini. Opt Quantum Electron., 53 (10) (2021), pp. 1-17

[8]
A.A. Kilbas, H.M. Srivastava, J.J. Trujillo, Theory and applications of fractional differential equations. Elsevier, Amsterdam (2006)

[9]
H. Durur, E. Ilhan,H. Bulut. Fractal Fract., 4 (3) (2020), p. 41 DOI: 10.3390/fractalfract4030041

[10]
H. Almusawa, K.A. Khalid, A.-M. Wazwaz, M. S. Mehanna, D. Baleanu, M. S. Osman,Results Phys., 31 (2021), p. 104933

[11]
K.A. Khalid, M. S. Mehanna,Results Phys., 28 (2021), p. 104596

[12]
K.A. Khalid, M. S. Mehanna. Opt,Quantum Electron., 53 (2021), p. 723

[13]
K.K. Ali, M. S. Mehanna. Alex, Eng. J., 60 (2021), pp. 5275-5285

[14]
A. Saha, K.A. Khalid, H. Rezazadeh,Y. Ghatani. Opt. Quantum Electron., 53 (2021), p. 150

[15]
A. Saha, K.A. Khalid, H. Rezazadeh,Y. Ghatani. Opt. Quantum Electron., 53 (2021), p. 150

[16]
A. Zafar, K.A. Khalid, M. Raheel, N. Jafar, K.S. Nisar. Eur.Phys, J.Plus, 135 (2020), p. 726

[17]
A. M. S. Mahdy. J. Ocean Eng. Sci. (2022), DOI:

[18]
A.M.S. Mahdy, K. Lotfy, W. Hassan, A.A. El-Bary, Waves Random Comp. Media, 31 (6) (2021), pp. 2040-2057 DOI: 10.1080/17455030.2020.1717673

[19]
M. Khader, N.H. Swetlam, A.M.S. Mahdy, Wseas Trans. Math., 13 (2014), pp. 31-38

[20]
A.M.S. Mahdy, K. Lotfy, A. El-Bary, H.H. Sarhan, Eur. Phy. J. Plus, 136 (5) (2021), pp. 1-17

[21]
S.B. GaziKarakoc, K.A. Khalid, Appl. Math. Comput., 371 (2020), p. 124933

[22]
A numerical. J. Ocean Eng. Sci., 3 (1) (2018), pp. 67-75

[23]
Ak Turgut S. Karakoc B. Gazi B, Anjan. J. Comput. Theor. Nanos., 13 (10) (2016), pp. 1-9

[24]
S.B. Karakoc, K.K. Ali. Comput, Methods Differ. Equ., 9 (3) (2021), pp. 670-691

[25]
International Bureau of Weights and Measures. SI Brochure: The International System of Units (SI) ( PDF) (9th ed.) (2020), p. 131 ISBN 978-92-822-2272-0.

[26]
H. Eleuch, Hindawi Adv. Math. Phys., 7 (2018), p. 5242757, DOI:

[27]
N. Raza, S. Arshed, H.I. Alrebdi, A.-H. Abdel-Aty,H. Eleuch. Results Phys., 35 (2022), p. 105358

[28]
T. Abdeljawad, D. Baleanu. Commun. Nonlinear Sci. Numer. Simul., 16 (12) (2011), pp. 4682-4688

[29]
R.P. Agarwal, Mathematical Proceedings of the Cambridge Philosophical Society, vol. 66, Cambridge University Press (1969), pp. 365-370

[30]
R.I. Butt, T. Abdeljawad, M.A. Alqudah, M.u. Rehman. J, Inequal. Appl., 2019 (1) (2019), pp. 1-13

[31]
M. Jleli, M. Mursaleen, B. Samet. Electron. J. Differ. Equ., 2016 (17) (2016), pp. 1-14

[32]
V.G. Kac, P. Cheung. Quantum calculus. Springer (2002), p. 113

[33]
K. Victor, C. Pokman, Q. calculus. Universitext. Springer-Verlag (2002) ISBN 0-387-95341-8.

[34]
K.A. Khalid, M. S. Mehanna. Opt,Quantum Electron., 53 (2021), p. 723

[35]
K.K. Ali, M. S. Mehanna. Alex, Eng. J., 60 (2021), pp. 5275-5285

[36]
K.R. Raslan,K.K. Ali. Nonlinear Eng., 9 (2020), pp. 432-448 DOI: 10.1515/nleng-2020-0027

[37]
K.R. Raslan, K.K. Ali, H.M.Y. Al-Bayatti. Appl. Math. Inf. Sci., 15 (5) (2021), pp. 599-611

[38]
K.R. Raslan, K.K. Ali, M.S. Mohamed, A.R. Hadhoud. Adv. Differ. Equ., 442 (1) (2021), p. 2021

[39]
K.R. Raslan, K.K. Ali, M.A. Shaalan, H.K. Al-Jeaid. Int. J, Appl. Comput. Math, 8 (2022), p. 91

Outlines

/