Original article

New-fashioned solitons of coupled nonlinear Maccari systems describing the motion of solitary waves in fluid flow

  • Tarikul Islam a ,
  • Ali Akbar b ,
  • Hadi Rezazadeh c ,
  • Ahmet Bekir , d
Expand
  • a Department of Mathematics, Hajee Mohammad Danesh Science and Technology University, Dinajpur, Bangladesh
  • b Department of Applied Mathematics, University of Rajshahi, Bangladesh
  • c Faculty of Modern Technologies Engineering, Amol University of Special Modern Technologies, Amol, Iran
  • d Neighbourhood of Akcaglan, Imarli Street, No. 28/4, Eskisehir 26030, Turkey

Received date: 2021-08-19

  Revised date: 2022-03-05

  Accepted date: 2022-03-06

  Online published: 2022-03-08

Abstract

The Schrodinger equation type nonlinear coupled Maccari system is a significant equation that flourished with the wide-ranging arena concerning fluid flow and the theory of deep-water waves, physics of plasma, nonlinear optics, etc. We exploit the enhanced tanh approach and the rational (G′/G) -expansion process to retrieve the soliton and dissimilar soliton solutions to the Maccari system in this study. The suggested systems of nonlinear equations turn into a differential equation of single variable through executing some operations of wave variable alteration. Thereupon, with the successful implementation of the advised techniques, a lot of exact soliton solutions are regained. The obtained solutions are depicted in 2D, 3D, and contour traces by assigning appropriate values of the allied unknown constants. These diverse graphical appearances assist the researchers to understand the underlying processes of intricate phenomena of the leading equations. The individual performances of the employed methods are praiseworthy which justify further application to unravel many other nonlinear evolution equations ascending in various branches of science and engineering.

Highlights

● We have obtained new solitons for the coupled nonlinear Maccari's system describing the motion of isolated waves in fluid-flow.

● This model describes propagation of solitons in the theory of deep-water waves.

● Abundant closed form wave solutions are successfully generated in terms of rational, trigonometric and hyperbolic functions.

● The acquired solutions alongside particular values of involved free parameters are figured out in 3D, 2D and contour profiles to depict diverse soliton patterns.

Cite this article

Tarikul Islam , Ali Akbar , Hadi Rezazadeh , Ahmet Bekir . New-fashioned solitons of coupled nonlinear Maccari systems describing the motion of solitary waves in fluid flow[J]. Journal of Ocean Engineering and Science, 2024 , 9(4) : 353 -363 . DOI: 10.1016/j.joes.2022.03.003

1. Introduction

Nonlinear fractional and classical differential equations are significant tools to model and depict the interior behavior of complicated physical phenomena arise throughout the natural world [1], [2], [3]. That is why the researchers have paid attention to explore the analytical and numerical solutions of these equations [4], [5], [6], [7]. The continued efforts of scholars to examine the results of nonlinear partial differential equations attained through diverse approaches are opening new doors of research [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21].
The nonlinear Maccari system [30], [31], [32] is the Schrodinger equation type coupled equations iqt+qxx+rq=0, ist+sxx+rs=0, int+nxx+rn=0,
rt+ry+(|q+s+n|2)x=0,
where q,s,n stand for complex scalar fields, r represents real scalar fields and all depends on space-variables x , y and time-variable t , has wide-spectral applications in the theory of deep-water waves, fluid flow, plasma physics, nonlinear optics, etc. This system has been investigated by some researchers to unravel for exact wave solutions, such as Ciancio et al [22]. utilized the modified exp(ϕ(η)) -expansion function procedure to find solutions relating to hyperbolic and trigonometric functions; Inc et al [23]. studied this system by the use of the generalized Riccati equation scheme and modified F-expansion method; Neirameh employed the improved (G/G) -expansion scheme to acquire certain closed-from wave solutions [24]; Baskonus et al [25]. the Maccari system (1.1) by means of sine-Gordon expansion technique and found only hyperbolic type solutions; the modified trial equation process has been adopted by Bulut et al [26]. which provide few traveling wave solutions; Rostamy and Zabihi [27] have applied the first integral method to construct exponential functions and hyperbolic functions wave solutions. Vahidi et al [28]. studied the Maccari system by means of rational sinh-cosh and extended rational sine-cosine techniques, which is closely related to Eq. (1.1); Akbar et al [29]. have also considered the same system to extract by a generalized auxiliary equation mapping algorithm. Now, we have a brief look at the literature to recall the original foundation of Maccari system. This system containing two nonlinear evolution equations (Maccari, [30]) stands out among the nonlinear models for understanding the movement of solitary waves in plasmas, nonlinear optics, quantum field theory. Several studies relating to this system is available in the literature for instance Jabbari et al [31]. used (G/G) -expansion and He's semi-inverse techniques which provide estimations in terms rational, trigonometric and hyperbolic form; Khater et al [32]. studied the system by using general exp-function approach and obtained standard soliton solutions; Cheemaa et al [33]. used a modified auxiliary equation approach and found different types closed-form solutions; Liu et al [34]. have studied the same system and obtained bright and dark N-soliton solutions, Hafez et al [35]. examined this system through the new (G/G) -expansion scheme.
Khater et al [36]. analyzed the nonlinear Kolmogorov-Petrovskii-Piskunov model of fractional order by the use of the Adomian decomposition procedure and numerical approach. Chu et al [37]. investigated the general (2 + 1)-dimensional wave model arise in shallow-water using the simplest calculation approach and the advanced Kudryashov scheme. Khater et al [38]. studied the semi-analytical and analytical solutions to the fractional HIV-1 infection of CD4+ T -cells, and also approximate wave solutions to the complex nonlinear Fokas-Lenells models [39] of short pulses in optical fibers, the semi-exact estimations of the nonlinear Phi-4 equation [40], the numerical simulation of the time-space fractional telegraphic nonlinear equation [41]. Khater and Ghanbari [42] investigated the soliton solutions to the Chaffee-Infante model using five techniques. Khater et al [43]. scrutinized numerical estimations to the modified BBM model and examined the stability property of the attained solutions, abundant novel soliton solutions to the Klein-Fock-Gordon equation [44], the analytical solutions to the perturbed Schrödinger nonlinear equation [45], the analytical and semi-exact estimations to the time-fractional Cahn-Allen model [46], the numerical estimations to the Klein-Gordon-Zakharov equation [47], the analytical accurate traveling wave solutions of the higher-dimensional Schrödinger nonlinear model [48], the semi-analytical and exact estimation to the fractional quadratic-cubic Schrödinger equation [49]. Khater [50] unraveled abundant breather solutions to the Vakhnenko-Parkes equation by means of a couple of analytical approaches. Khater et al [51]. investigated the numerical solutions of two bio-mathematical and statistical fractional order nonlinear equations. Khater and Lu [52] searched numerical and exact solutions to the space and time fractional nonlinear telegraphic model using the trigonometric-quantic-B-spline method. Khater and Ahmed [53] evaluated the solutions of the Klein-Gordon-Zakharov nonlinear model by exponential cubic B-spline schemes. Khater et al [54]. derived bright-dark soliton to the (1 + 1)-dimensional nonlinear Kaup-Kupershmidt model utilizing modified F-expansion approach. Yue et al [55]. perused analytical estimations to the fractional Hirota-Satsuma models via the modified Kudryashov and the enhanced simplest equation tools. Li et al [56]. probed numerical estimations to the (2 + 1)-dimensional KPBBM model by employing the modified direct algebraic method. Khater and Alabdali [57] studied numerical and analytical wave solutions of (2 + 1)-dimensional Fisher-KolmogorovPetrovskii-Piskunov model by imposing the modified Kudryashov, trigonometric-quantic B-spline, and modified direct algebraic schemes.
To the furthest of our sagaciousness, the Maccari system (1.1) has not been studied using the enhanced tanh approach and the method of rational (G/G) -expansion. In continuation of the above works fascinate us and thereupon we aim to extract some primal and standard soliton solutions to the Maccari system through the indicated methods from which a number of reputed solutions will be re-established as well as some effective and reliable solitary solutions will be extracted. The closed waveform solutions are received in hyperbolic, rational, and trigonometric types. The well-furnished solutions are interpreted in 2D, 3D, and contour traces to bring out their physical meanings. It is designated that this work is new in the literature, and the results might have a substantial impact on ocean engineering [22], [23], [24], [25], [26], [27]. The Journal of Ocean Engineering and Science has published some results of various investigations [58], [59], [60], [61], [62], [63], [64] that are compatible with our findings. The considered governing model is concerned with the phenomena of water waves, fluid flow, plasma physics etc. which are related to wave energy converters and hence oceanography.

2. Explanation of the methods

A nonlinear evolution equation is assumed as
Θ(u,ux,uy,ut,uxy,uxt,uyt,uxx,uyy,utt,)=0,
where Θ consists of u and it distinct derivatives. The new wave variable adaptation
u=u(ϕ),ϕ=ϕ(x,y,t,)
makes this equation to the following general differential equation relating to ϕ :
Λ(U,U,U,U,)=0.
As many times as feasible, we take the anti-derivative of Eq. (2.3). We can ignore the integral constant, since we are scrutinizing soliton solutions. The following are the key steps of the introduced techniques:

2.1. The improved tanh method

We estimate the solution of Eq. (2.3) in accordance with the improved tanh approach in the ensuing form
Q(ϕ)=ι0+i=1e(ιiΦi(ϕ)+τiΦi(ϕ))0+i=1e(iΦi(ϕ)+εiΦi(ϕ))
whose unknown parameters are evaluated hereafter [65]. The integral value of e is set in view of the balancing principle applying to Eq. (2.3). The Riccati equation satisfying by Φ=Φ(ϕ) is
Φ(ϕ)=δ+Φ2(ϕ),
where δ is a real arbitrary number Eq. (2.1.2). has the subsequent outcomes:
Φ(ϕ)=δtanh(δϕ)orΦ(ϕ)=δcoth(δϕ),δ<0
Φ(ϕ)=1/ϕorδ=0
Φ(ϕ)=δtan(δϕ)orΦ(ϕ)=δcot(δϕ),δ>0
When Eq. (2.3) is combined with Eqs. (2.1.1) and (2.1.2), a polynomial in Φ(ϕ) is obtained. We unravel the equations using computational software that equalizes the coefficients of this polynomial to zero for the values of the arbitrary numbers belonging to Eq. (2.1.1). We substitute the received values in Eq. (2.1.1), which, when combined with the solutions of Eq. (2.1.2), result the exact traveling wave solutions of Eq. (2.1).

2.2. The rational (G/G) -expansion method

In the following form,
Q(ϕ)=i=0eιi(G(ϕ)/G(ϕ))ii=0eτi(G(ϕ)/G(ϕ))i,
we estimate the solution of Eq. (2.3) using the enhanced tanh technique, where e is determined by equilibrium state the linear and nonlinear terms of the maximum order appearing in Eq. (2.3) and unknown constants ais and bis are calculated after some operations [66]. The function G(ϕ) satisfies
G(ϕ)+bG(ϕ)+cG(ϕ)=0
This ODE leaves,
G(ϕ)G(ϕ)=b+b24ctanh(0.5b24cϕ)/2,b24c>0
G(ϕ)G(ϕ)={b+b24ccoth(0.5b24cϕ)}/2,b24c>0
G(ϕ)G(ϕ)=1ϕb2,b24c=0
G(ϕ)G(ϕ)={b4cb2tan(0.54cb2ϕ)}/2,b24c<0
G(ϕ)G(ϕ)={b+4cb2cot(0.54cb2ϕ)}/2,b24c<0
With the help of Eqs. (2.2.1) and (2.2.2), from Eq. (2.3) it is resulted a polynomial in (G(ϕ)/G(ϕ)) with zero coefficients, yielding a set of algebraic equations. We use Maple to unravel the system and seek for the values of the arbitrary numbers occur in Eq. (2.2.1). The accurate traveling wave solutions of Eq. (2.1) are then obtained by combining these values into solution (2.2.1).

3. Formulation of solutions

In this part, the nonlinear coupled Maccari system is considered to be unraveled by means of two competent tools like the approach of rational (G/G) -expansion and improved tanh function method. We introduce the transformation, q(x,y,t)=Q(ϕ)eiψ,s(x,y,t)=S(ϕ)eiψ
n(x,y,t)=N(ϕ)eiψ,r(x,y,t)=R(ϕ)
where ϕ=x+y2κt , ψ=κx+αy+γt+ζ .
After specific operation on Eq. (1.1) by the support of Eq. (3.1) turns into the resulting equation:
(12κ)Q(12κ)(γ+κ2)Q(1+λ+μ)2Q3=0,
where S=λQ , N=μQ and R=(1+λ+μ)2Q22κ1 .
Balancing the terms Q and Q3 we determine e=1 . Thereupon, the advised techniques are adopted to seek for solitary wave solutions of Eq. (1.1) in appropriate form.

3.1. The improved tanh method

Q(ϕ)=ι0+ι1Φ(ϕ)+τ1Φ1(ϕ)0+1Φ(ϕ)+ε1Φ1(ϕ)
By the assistance of solution (3.1.1) and wave transformation (2.1.1), Eq. (3.2) becomes a polynomial in Φ(ϕ) . By setting zero to all coefficients of this polynomial and solving those with the computer package, we get the results
Set 1:
ι0=ι1=0,0=±2τ1(1+λ+μ)(γ+κ2)24κ,1=ε1=0,δ=γ+κ22
Set 2:
ι0=0,τ1=ι1(γ+κ2)4,0=±ι1(1+λ+μ)24κ,1=0,δ=γ+κ24
Set 3:
ι0=±ε124κ1+λ+μ,ι1=1=0,τ1=ε1(24κ)(2γ+2κ2)2(1+λ+μ),0=2ε12γ+2κ2,δ=γ+κ22
Set 4:
ι0=±ε124κ1+λ+μ,ι1=0,τ1=0(γ+κ2)(24κ)2(1+λ+μ),1=0,δ=γ+κ22
Set 5:
ι0=0,τ1=ι1(γ+κ2)8,0=±ι1(1+λ+μ)24κ,1=0,ε1=0,δ=γ+κ28
Set 6:
ι0=±ε124κ1+λ+μ,τ1=0,0=±ι1(1+λ+μ)24κ,1=0,δ=γ+κ22
Injecting the scores appeared in Eq. (3.1.2) into Eq. (3.1.1) and thereupon adopting Eqs. (2.1.3)-(2.1.5) provide the following solitary wave solutions:
For δ<0 , the solutions are received in hyperbolic function form as
q1(x,y,t)=±(γ+κ2)24κeiψcoth(δϕ)2τ1δ(1+λ+μ),
q2(x,y,t)=±(γ+κ2)24κeiψtanh(δϕ)2τ1δ(1+λ+μ),
s1(x,y,t)=±λ(γ+κ2)24κeiψcoth(δϕ)2τ1δ(1+λ+μ),
s2(x,y,t)=±λ(γ+κ2)24κeiψtanh(δϕ)2τ1δ(1+λ+μ),
n1(x,y,t)=±μ(γ+κ2)24κeiψcoth(δϕ)2τ1δ(1+λ+μ),
n2(x,y,t)=±μ(γ+κ2)24κeiψtanh(δϕ)2τ1δ(1+λ+μ),
r1(x,y,t)=(γ+κ2)22τ12δcoth2(δϕ),
r2(x,y,t)=(γ+κ2)22τ12δtanh2(δϕ),
where ϕ=x+y2κt and ψ=κx+αy+γt+ζ .
Under the assumption δ=0 the rational function solutions are
q3(x,y,t)=±(γ+κ2)24κϕeiψ2τ1(1+λ+μ),
s3(x,y,t)=±λ(γ+κ2)24κϕeiψ2τ1(1+λ+μ),
n3(x,y,t)=±μ(γ+κ2)24κϕeiψ2τ1(1+λ+μ),
r3(x,y,t)=(γ+κ2)2ϕ22τ12,
where ϕ=x+y2κt and ψ=κx+αy+γt+ζ .
For the postulate δ>0 the trigonometric solutions are appeared as
q4(x,y,t)=±(γ+κ2)24κeiψcot(δϕ)2τ1δ(1+λ+μ),
q5(x,y,t)=±(γ+κ2)24κeiψtan(δϕ)2τ1δ(1+λ+μ),
s4(x,y,t)=±λ(γ+κ2)24κeiψcot(δϕ)2τ1δ(1+λ+μ),
s5(x,y,t)=±λ(γ+κ2)24κeiψtan(δϕ)2τ1δ(1+λ+μ),
n4(x,y,t)=±μ(γ+κ2)24κeiψcot(δϕ)2τ1δ(1+λ+μ),
n5(x,y,t)=±μ(γ+κ2)24κeiψtan(δϕ)2τ1δ(1+λ+μ),
r4(x,y,t)=(γ+κ2)2cot2(δϕ)2τ12δ,
r5(x,y,t)=(γ+κ2)2tan2(δϕ)2τ12δ,
where ϕ=x+y2κt and ψ=κx+αy+γt+ζ .
Utilizing the values involved in Eq. (3.1.3) into Eq. (3.1.1) and then adopting Eqs. (2.1.3)-(2.1.5) deliver the closed form traveling wave solutions as bellow:
Considering δ<0 the subsequent hyperbolic form solutions are attained
q6(x,y,t)=±24κeiψ4δ(1+λ+μ)×{4δtanh(δϕ)+(γ+κ2)coth(δϕ)},
q7(x,y,t)=±24κeiψ4δ(1+λ+μ)×{4δcoth(δϕ)+(γ+κ2)tanh(δϕ)},
s6(x,y,t)=±λ24κeiψ4δ(1+λ+μ)×{4δtanh(δϕ)+(γ+κ2)coth(δϕ)},
s7(x,y,t)=±λ24κeiψ4δ(1+λ+μ)×{4δcoth(δϕ)+(γ+κ2)tanh(δϕ)},
n6(x,y,t)=±μ24κeiψ4δ(1+λ+μ)×{4δtanh(δϕ)+(γ+κ2)coth(δϕ)},
n7(x,y,t)=±μ24κeiψ4δ(1+λ+μ)×{4δcoth(δϕ)+(γ+κ2)tanh(δϕ)},
r6(x,y,t)={4δtanh(δϕ)+(γ+κ2)coth(δϕ)}28δ,
r7(x,y,t)={4δcoth(δϕ)+(γ+κ2)tanh(δϕ)}28δ,
where ϕ=x+y2κt and ψ=κx+αy+γt+ζ .
The rational type solutions with the assumption δ=0 are achieved as
q8(x,y,t)=±24κeiψ4(1+λ+μ)×{(γ+κ2)ϕ4ϕ1},
s8(x,y,t)=±λ24κeiψ4(1+λ+μ)×{(γ+κ2)ϕ4ϕ1},
n8(x,y,t)=±μ24κeiψ4(1+λ+μ)×{(γ+κ2)ϕ4ϕ1},
r8(x,y,t)={(γ+κ2)ϕ4ϕ1}28,
where ϕ=x+y2κt and ψ=κx+αy+γt+ζ .
Under the hypothesis δ>0 we retrieve the trigonometric function solutions as follows:
q9(x,y,t)=±24κeiψ4δ(1+λ+μ)×{4δtan(δϕ)(γ+κ2)cot(δϕ)},
q10(x,y,t)=±24κeiψ4δ(1+λ+μ)×{(γ+κ2)tan(δϕ)4δcot(δϕ)},
s9(x,y,t)=±λ24κeiψ4δ(1+λ+μ)×{4δtan(δϕ)(γ+κ2)cot(δϕ)},
s10(x,y,t)=±λ24κeiψ4δ(1+λ+μ)×{(γ+κ2)tan(δϕ)4δcot(δϕ)},
n9(x,y,t)=±μ24κeiψ4δ(1+λ+μ)×{4δtan(δϕ)(γ+κ2)cot(δϕ)},
n10(x,y,t)=±μ24κeiψ4δ(1+λ+μ)×{(γ+κ2)tan(δϕ)4δcot(δϕ)},
r9(x,y,t)={4δtan(δϕ)(γ+κ2)cot(δϕ)}28δ,
r10(x,y,t)={(γ+κ2)tan(δϕ)4δcot(δϕ)}28δ,
where ϕ=x+y2κt , ψ=κx+αy+γt+ζ .
The above procedure after assimilation Eqs. (3.1.1), (3.1.4)-(3.1.7) and (2.1.3)-(2.1.5) might provide abundant wave solutions in closed form as hyperbolic, rational, and trigonometric types which are not recorded in this text for simplicity.

3.2. The rational (G/G) -expansion method

Q(ϕ)=ι0+ι1(G(ϕ)/G(ϕ))τ0+τ1(G(ϕ)/G(ϕ)).
Eq. (3.2) is used to reduce a polynomial in (G(ϕ)/G(ϕ)) with the help of Eqs. (3.2.1) and (2.2.2). For an over determined system of algebraic equations, we collect each term and set it to zero. The following are the results of solving this system with Maple.
Set 1:
ι0=±(2cτ1bτ0)24κ2(1+λ+μ),ι1=±(bτ12τ0)24κ2(1+λ+μ),γ=4cb22κ22
Set 2:
ι0=±bτ024κ2(1+λ+μ),ι1=±τ024κ1+λ+μ,τ1=0,γ=4cb22κ22
Set 3:
ι0=±τ1(b24c)24κ4(1+λ+μ),ι1=0,τ0=τ1b2,γ=4cb22κ22
Merging Eqs. (3.2.1), (3.2.2) and (2.2.3)-(2.2.7), we formulate the succeeding exact solitary wave solutions in three types such as hyperbolic, rational and trigonometric functions:
For the assumption b24c>0 , the exact traveling wave solutions are obtained in hyperbolic form as
q1(x,y,t)=±24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+b24ctanh(0.5b24cϕ)}2τ0+τ1{b+b24ctanh(0.5b24cϕ)},
q2(x,y,t)=±24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+b24ccoth(0.5b24cϕ)}2τ0+τ1{b+b24ccoth(0.5b24cϕ)},
s1(x,y,t)=±λ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+b24ctanh(0.5b24cϕ)}2τ0+τ1{b+b24ctanh(0.5b24cϕ)},
s2(x,y,t)=±λ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+b24ccoth(0.5b24cϕ)}2τ0+τ1{b+b24ccoth(0.5b24cϕ)},
n1(x,y,t)=±μ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+b24ctanh(0.5b24cϕ)}2τ0+τ1{b+b24ctanh(0.5b24cϕ)},
n2(x,y,t)=±μ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+b24ccoth(0.5b24cϕ)}2τ0+τ1{b+b24ccoth(0.5b24cϕ)},
r1(x,y,t)=12(2(2cτ1bτ0)+(bτ12τ0){b+b24ctanh(0.5b24cϕ)}2τ0+τ1{b+b24ctanh(0.5b24cϕ)})2,
r2(x,y,t)=12(2(2cτ1bτ0)+(bτ12τ0){b+b24ccoth(0.5b24cϕ)}2τ0+τ1{b+b24ccoth(0.5b24cϕ)})2,
where ϕ=x+y2κt and ψ=κx+αy+4cb22κ22t+ζ .
When b24c=0 , the rational function solutions are gained as
q3(x,y,t)=±24κeiψ2(1+λ+μ)×(2cτ1bτ0)+(bτ12τ0)(1ϕb2)τ0+τ1(1ϕb2),
s3(x,y,t)=±λ24κeiψ2(1+λ+μ)×(2cτ1bτ0)+(bτ12τ0)(1ϕb2)τ0+τ1(1ϕb2),
n3(x,y,t)=±μ24κeiψ2(1+λ+μ)×(2cτ1bτ0)+(bτ12τ0)(1ϕb2)τ0+τ1(1ϕb2),
r3(x,y,t)=12((2cτ1bτ0)+(bτ12τ0)(1ϕb2)τ0+τ1(1ϕb2))2,
where ϕ=x+y2κt and ψ=κx+αy+4cb22κ22t+ζ .
Under the assumption b24c<0 the wave solutions in the form of trigonometric function are constructed as,
q4(x,y,t)=±24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b4cb2tan(0.54cb2ϕ)}2τ0+τ1{b4cb2tan(0.54cb2ϕ)},
q5(x,y,t)=±24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+4cb2cot(0.54cb2ϕ)}2τ0+τ1{b+4cb2cot(0.54cb2ϕ)},
s4(x,y,t)=±λ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b4cb2tan(0.54cb2ϕ)}2τ0+τ1{b4cb2tan(0.54cb2ϕ)},
s5(x,y,t)=±λ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+4cb2cot(0.54cb2ϕ)}2τ0+τ1{b+4cb2cot(0.54cb2ϕ)},
n4(x,y,t)=±μ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b4cb2tan(0.54cb2ϕ)}2τ0+τ1{b4cb2tan(0.54cb2ϕ)},
n5(x,y,t)=±μ24κeiψ2(1+λ+μ)×2(2cτ1bτ0)+(bτ12τ0){b+4cb2cot(0.54cb2ϕ)}2τ0+τ1{b+4cb2cot(0.54cb2ϕ)},
r4(x,y,t)=12(2(2cτ1bτ0)+(bτ12τ0){b4cb2tan(0.54cb2ϕ)}2τ0+τ1{b4cb2tan(0.54cb2ϕ)})2,
r5(x,y,t)=12(2(2cτ1bτ0)+(bτ12τ0){b+4cb2cot(0.54cb2ϕ)}2τ0+τ1{b+4cb2cot(0.54cb2ϕ)})2,
where ϕ=x+y2κt and ψ=κx+αy+4cb22κ22t+ζ .
We assimilate Eqs. (3.2.1), (3.2.3) and (2.2.3)-(2.2.7) to celebrate the following three types exact traveling wave solutions:
Imposing the condition b24c>0 the solutions as hyperbolic type are found as below,
q6(x,y,t)=±(24κ)(b24c)eiψ2(1+λ+μ)coth(0.5b24cϕ),
q7(x,y,t)=±(24κ)(b24c)eiψ2(1+λ+μ)tanh(0.5b24cϕ),
s6(x,y,t)=±λ(24κ)(b24c)eiψ2(1+λ+μ)coth(0.5b24cϕ),
s7(x,y,t)=±λ(24κ)(b24c)eiψ2(1+λ+μ)tanh(0.5b24cϕ),
n6(x,y,t)=±μ(24κ)(b24c)eiψ2(1+λ+μ)coth(0.5b24cϕ),
n7(x,y,t)=±μ(24κ)(b24c)eiψ2(1+λ+μ)tanh(0.5b24cϕ),
r6(x,y,t)=4cb22coth2(0.5b24cϕ),
r7(x,y,t)=4cb22tanh2(0.5b24cϕ),
where ϕ=x+y2κt and ψ=κx+αy+4cb22κ22t+ζ .
For b24c=0 , the following rational type solutions are gained:
q8(x,y,t)=±24κeiψ(1+λ+μ)ϕ,
s8(x,y,t)=±λ24κeiψ(1+λ+μ)ϕ,
n8(x,y,t)=±μ24κeiψ(1+λ+μ)ϕ,
r8(x,y,t)=2ϕ2,
where ϕ=x+y2κt and ψ=κx+αy+4cb22κ22t+ζ .
When b24c<0 , the ensuing solutions in trigonometric form are
q9(x,y,t)=±(24κ)(4cb2)eiψ2(1+λ+μ)cot(0.54cb2ϕ),
q10(x,y,t)=±(24κ)(4cb2)eiψ2(1+λ+μ)tan(0.54cb2ϕ),
s9(x,y,t)=±λ(24κ)(4cb2)eiψ2(1+λ+μ)cot(0.54cb2ϕ),
s10(x,y,t)=±λ(24κ)(4cb2)eiψ2(1+λ+μ)tan(0.54cb2ϕ),
n9(x,y,t)=±μ(24κ)(4cb2)eiψ2(1+λ+μ)cot(0.54cb2ϕ),
n10(x,y,t)=±μ(24κ)(4cb2)eiψ2(1+λ+μ)tan(0.54cb2ϕ),
r9(x,y,t)=b24c2cot2(0.54cb2ϕ),
r10(x,y,t)=b24c2tan2(0.54cb2ϕ),
where ϕ=x+y2κt and ψ=κx+αy+4cb22κ22t+ζ .
If we assemble Eqs. (3.2.1), (3.2.4) and (2.2.3)-(2.2.7) and also maintain the same procedure as above, then more exact wave solutions might be found which has been avoided to record in the text.
Remarks: The rational (G/G) -expansion approach and improved tanh scheme have successfully been employed to examine the coupled nonlinear Maccari system of Schrodinger type and gained numerous significant and interesting wave solutions. Ciancio et al. have applied the modified exp(ϕ(η)) -expansion function procedure and found trigonometric and hyperbolic function solutions [22]. Inc et al. have also studied this system by using the modified F-expansion process and the generalized Riccati equation method [23]. We claim that our obtained outcomes are comparatively different, novel and significant.

4. Results discussion and graphical representations

We have made a comparable study of the well-generated wave solutions with the earlier results discussed in the introduction section and claimed that this exploration provides novel and diverse exact analytic solutions in trigonometric, hyperbolic and rational forms [22], [23], [24], [25], [26], [27]. The acquired solutions are figured out in the profiles of 3D, 2D and contour to illustrate various soliton patterns. The portrayed outlines stand for bell shape, kink type, anti-bell shape, anti-kink type, singular bell shape, singular kink type, compacton, periodic, peakon, cuspon, singular periodic etc. The 2D profiles are considered to depict the wave propagation due to the change of the velocity of wave. It is worthy to observe that the nature of the solutions obtained depends on the wave velocity ( 2κ ) and affects the propagation of the wave. We bring out different 2D profiles with the several particular values of κ which clarify various patterns of wave due to the change of wave velocity. Few graphical representations are as follows:
Outline of solution (3.1.9) is appeared as kink type: Fig. 1(a), (b) stand for 3D, contour shapes for κ=0.25 , α=λ=0.5 , γ=0.8 , ζ=μ=1 , δ=1 , τ1=0.2 , y=0 in the range 10x,t10 while Fig. 1(c) represents a class of 2D profiles with the different values of κ (green (κ=0.1 ), red (κ=0.2 ), pink (κ=0.3 ), blue (κ=0.4 )) and t=2 . Plot of solution (3.1.15) takes bell shape: 3D, contour outlines are appeared in Fig. 2(a), (b) for κ=0.5 , γ=τ1=2 , y=0 , δ=1 in the interval 3x,t3
whereas various 2D shapes are shown in Fig. 2(c) by assigning the different values of κ (green (κ=0.25 ), red (κ=0.5 ), pink (κ=0.75 ), blue (κ=1 )) and t=1 . Scheme of solution (3.1.19) stands as compacton: the assumption of the parameters as κ=0.5 , γ=0.8 , τ1=2 , y=0 bring out the 3D, contour outlines shown in Fig. 3(a), (b) within the interval 10x,t10 and the distinct values of κ (green (κ=0.25 ), red (κ=0.5 ), pink (κ=0.75 ), blue (κ=1 )) and t=0.5 display diverse 2D profiles in Fig. 3(c). Sketch of solution (3.1.39) looks like anti-bell shape: Fig. 4(a), (b) represent 3D, contour schemes for κ=0.25 , γ=0.5 , y=0 within the range 15x,t15 while Fig. 4(c) stands for a class of 2D profiles under the different values of κ (green (κ=0.25 ), red (κ=0.5 ), pink (κ=0.75 ), blue (κ=1 )) and t=1 . Outline of solution (3.1.40) is appeared as periodic: 3D, contour shapes appeared in Fig. 5(a), (b) are drawn for the parameters values κ=0.25 , α=0.5 , γ=ζ=μ=1, δ=3, λ=1 , y=0 in the interval 15x,t15
whereas various 2D sketches under the several values of κ (green (κ=0.1 ), red (κ=0.2 ), pink (κ=0.3 ), blue (κ=0.4 )) and t=1 are displayed in Fig. 5(c). Plot of solution (3.2.5) stands for kink type: the supposition of arbitrary parameters as κ=0.25 , α=λ=τ0=0.5 , ζ=μ=1 , τ1=y=0 , b=4 , c=3 ensures 3D, contour shapes presented in Fig. 6(a), (b) within 10x,t10 and the different values of κ (green (κ=0.4 ), red (κ=0.3 ), pink (κ=0.2 ), blue (κ=0.1 )) and t=1 exhibit various 2D profiles in Fig. 6(c). Sketch of solution (3.2.24) is in periodic form: Fig. 7(a), (b) stand for 3D, contour are drawn with the particular values κ=0.5 , τ0=1 , τ1=2 , y=0 , b=4 , c=5 within the interval 13x,t13 while Fig. 7(c) is appeared for several 2D profiles along with the different values of κ (green (κ=0.25 ), red (κ=0.5 ), pink (κ=0.75 ), blue (κ=1 )) and t=1 . Outline in bell shape of solution (3.2.31): 3D, contour profiles in Fig. 8(a), (b) are obtained after assigning the parameters as κ=0.5 , y=0 , b=4 , c=3 in the range 3x,t3 and different 2D sketches in Fig. 8(c) are portrayed with the several values of κ (green (κ=0.5 ), red (κ=0.25 ), pink (κ=0.75 ), blue (κ=1 )) and t=1 . Singular kink type plot of solution (3.2.33): the assumption of the arbitrary constants as κ=0.25 , α=ζ=λ=0.5 , μ=c=1 , y=0 , b=2 provides 3D, contour shapes presented in Fig. 9(a), (b) within 15x,t15 while the different values of κ (green (κ=0.4 ), red (κ=0.3 ), pink (κ=0.2 ), blue (κ=0.1 )) and t=1 display various 2D outlines in Fig. 9(c). Anti-bell shape outline of solution (3.2.36): Fig. 10(a), (b) displaying 3D, contour profiles are brought out after assigning the unknown parameters as κ=0.25 , y=0 within the range 3x,t3 and different 2D sketches in Fig. 8(c) are depicted with the several values of κ (green (κ=0.25 ), red (κ=0.5 ), pink (κ=0.75 ), blue (κ=1 )) and t=1 . The diverse 3d, 2D and contour profiles are much more significant for describing the novel structures of solitons to the considered Maccari's system.
Fig.1 (a), (b) stand for 3D, contour shapes for κ=0.25, α=λ=0.5, γ=0.8, ζ=μ=1, δ=-1, τ1=0.2, y=0, in the range -10≤x,t≤10 while Fig 1(c) represents a class of 2D profiles with the different values of κ (green (κ=0.1), red (κ=0.2), pink (κ=0.3), blue (κ=0.4)) and t=2.
Fig.2 (a), 2(b) stand for 3D, contour shapes for κ=-0.5, γ=τ1=2, δ=-1, y=0, in the range -3≤x,t≤3 while Fig. 2(c) represents a class of 2D profiles with the different values of κ (green (κ=-0.25), red (κ=-0.5), pink (κ=-0.75), blue (κ=-1)) and t=1.
Fig.3 (a), (b) within the interval -10≤x,t≤10 and the distinct values of κ (green (κ=-0.25), red (κ=-0.5), pink (κ=-0.75), blue (κ=-1)) and t=0.5 display diverse 2D profiles in Fig. 3(c).
Fig.4 (a), (b) stand for 3D, contour shapes for κ=0.25, γ=0.5, y=0, in the range -15≤x,t≤15 while Fig. 4(c) represents a class of 2D profiles with the different values of κ (green (κ=0.25), red (κ=0.5), pink (κ=0.75), blue (κ=1)) and t=1.
Fig.5 (a), (b) stand for 3D, contour shapes for κ=0.25, α=0.5, λ=-1, γ=ζ=μ=1, y=0, in the range -15≤x,t≤15 while Fig. 5(c) represents a class of 2D profiles with the different values of κ (green (κ=0.1), red (κ=0.2), pink (κ=0.3), blue (κ=0.4)) and t=1.
Fig.6 (a), (b) within the interval -10≤x,t≤10 and the distinct values of κ (green (κ=0.4), red (κ=0.3), pink (κ=0.2), blue (κ=0.1)) and t=1 display diverse 2D profiles in Fig. 6(c).
Fig.7 (a), (b) stand for 3D, contour shapes for κ=-0.5, τ0=1, τ1=2, b=4, c=5, y=0, in the range -13≤x,t≤13 while Fig. 7(c) represents a class of 2D profiles with the different values of κ (green (κ=-0.25), red (κ=-0.5), pink (κ=-0.75), blue (κ=-1)) and t=1.
Fig.8 (a), (b) stand for 3D, contour shapes for κ=0.5, b=4, c=3, y=0, in the range -3≤x,t≤3 while Fig. 8(c) represents a class of 2D profiles with the different values of κ (green (κ=0.5), red (κ=0.25), pink (κ=0.75), blue (κ=1)) and t=1.
Fig.9 (a), (b) within the interval -15≤x,t≤15 and the distinct values of κ (green (κ=0.4), red (κ=0.3), pink (κ=0.2), blue (κ=0.1)) and t=1 display diverse 2D profiles in Fig. 9(c).
Fig.10 (a), (b) stand for 3D, contour shapes for κ=0.25, y=0, in the range -3≤x,t≤3 while Fig. 10(c) represents a class of 2D profiles with the different values of κ (green (κ=0.25), red (κ=0.5), pink (κ=0.75), blue (κ=1)) and t=1.

5. Conclusions

In this study, the coupled nonlinear Maccari system has been studied and extracted analytic solutions of by means of two competent techniques namely improved tanh method and rational (G/G) -expansion scheme. Abundant closed-form soliton solutions are successfully generated in terms of rational, hyperbolic and trigonometric functions. The acquired solutions alongside particular values of involved free parameters are illustrated in 2D, 3D and contour traces to depict diverse soliton patterns. The soliton patterns are appeared as kink type, bell shape, periodic, compacton, cuspon, peakon, etc. Different wave propagations due to the change of wave velocity are brought out through various 2D outlines. The relative study of the obtained outcomes and the results available in the knowledge has made us to claim the novelty and diversity of the achieved solution in this work. The adopted tools have performed as concise, efficient and productive to make available a heap of exact solutions which might convey a significant contribution to depict the intrinsic behaviors of nonlinear natural phenomena. The future research relevant to this study considering the mentioned techniques and any nonlinear evolution equations is just waiting for a few moments.
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]
A.M. Wazwaz. Partial Differential Equations: Method and Applications. Taylor and Francis ( 2002)

[2]
K.S. Miller, B. Ross. An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York (1993)

[3]
A.A. Kilbas, H.M. Srivastava, J.J. Trujillo. Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam (2006)

[4]
M. Bilal, A.R. Seadawy, M. Younis, S.T.R. Rizvi, K. El-Rashidy, S.F. Mahmoud. Res. Phys., 23 ( 2021), Article 103959

[5]
A. Yokus, H.M. Baskonus, T.A. Sulaiman, H. Bulut. Numer. Method Partial Differ. Eq., 34 ( 2018), pp. 211-227 DOI: 10.1002/num.22192

[6]
M.T. Islam,M.A. Akter. J. Res. Eng. Appl. Sci., 6 (2021), pp. 11-18 DOI: 10.1117/12.2596325

[7]
H. Ahmad, A.R. Seadawy, T.A. Khan,P. Thounthong. J. Taibah Univ. Sci., 14 (2020), pp. 346-358 DOI: 10.1080/16583655.2020.1741943

[8]
Y. Hu, Y. Luo, Z. Lu. J. Comput. Appl. Math., 215 (2008), pp. 220-229

[9]
S. Kumar, A. Kumar, A.M. Wazwaz. Eur. Phys. J. Plus., 135 (2020), pp. 1-17

[10]
A. Akbulut, M. Kaplan, A. Bekir. Int. J. Nonlinear Sci. Numer. Simul., 17 ( 2016), pp. 413-420 DOI: 10.1515/ijnsns-2016-0023

[11]
H. Bulut, H.M. Baskonus, Y. Pandir. Abst. Appl. Anal. (2013), Article 6368022013

[12]
H.Y. Martinez, J.F.G. Aguilar, A. Atangana. Math. Model. Nat. Phenom., 13 (2018), p. 14

[13]
M.A. Akbar, N.H.M. Ali, M.T. Islam.AIMS Math., 4 (2019), pp. 397-411 DOI: 10.32520/stmsi.v8i3.521

[14]
W.X. Ma, J.H. Lee. Chaos Solitons Fractals, 42 (2009), pp. 1356-1363

[15]
A. Yokus, H. Durur, S. Duran. Opt. Quantum Electron., 53 (402) (2021), pp. 1-15

[16]
M.T. Islam, M.A. Akbar, A.K. Azad. J. Ocean. Eng. Sci., 3 (2018), pp. 76-81

[17]
S.T.R. Rizvi, A.R. Seadawy, M. Younis, S. Iqbal, S. Althobaiti, A.M. El-Shehawi. Res. Phys., 23 ( 2021), Article 103998

[18]
M.T. Islam, M.E. Islam, M.A. Akbar. Arab J. Basic Appl. Sci.,28 (2021), pp. 64-72 DOI: 10.1080/25765299.2021.1886678

[19]
M. Younis, N. Cheemaa, S.A. Mehmood, S.T.R. Rizvi, A. Bekir. Com. Med., 30 (2018), pp. 490-499

[20]
A. Bhalraj, A. Azmi, M.H. Mohd. Int. J. Math. Com. Sci., 16 (2021), pp. 949-961

[21]
A.A. Al-Shawba, K.A. Gepreel, F.A. Abdullah, A. Azmi. Res. Phys., 9 (2018), pp. 337-343

[22]
A. Ciancio, H.M. Baskonus, T.A. Sulaiman, H. Bulut. Indian J. Phys., 92 ( 2018), pp. 1281-1290 DOI: 10.1007/s12648-018-1204-6

[23]
M. Inc, A.I. Aliyu, A. Yusuf, D. Baleanu, E. Nuray. Mod. Phys. Lett. B, 32 (2018), Article 1850014

[24]
A. Neirameh. Alex. Eng. J., 55 (2016), pp. 2839-2847

[25]
H.M. Baskonus, T.A. Sulaiman, H. Bulut. Opt. (Stuttg), 131 (2017), pp. 1036-1043

[26]
H. Bulut, G. Yel, H.M. Baskonus. Adv. Math. Mod. Appl., 2 (2017), pp. 14-19

[27]
D. Rostamy, F. Zabihi. Nonlinear Studies, 19 (2012), pp. 229-239

[28]
J. Vahidi, S.M. Zekavatmand, H. Rezazadeh, M. Inc, M.A. Akinlar,Y. Chu. Res. Phys., 21 (2021), Article 103801

[29]
Y. Akbar, H. Afsar, S. Abbas, M.W. Javed,N. Ullah. Chaos Solitons Fractals, 150 (2021), Article 111114

[30]
A. Maccari. J. Math. Phys., 12 (1996), pp. 6207-6212

[31]
A. Jabbari, H. Kheiri, A. Bekir. Comput. Math. Appl., 62 (2011), pp. 2177-2186

[32]
M.M.A. Khater, A.R. Seadawy, D. Lu. J. King Saud Univ. Sci., 30 (2018), pp. 417-423

[33]
N. Cheemaa, S. Chen, A.R. Seadawy. Res. Phys., 17 ( 2020), Article 102987

[34]
L. Liu, B. Tian, Y.Q. Yuan, Y. Sun. Eur. Phys. J. Plus., 133 (2018), pp. 1-14

[35]
M.G. Hafez, B. Zheng, M.A. Akbar. Egypt. J. Basic Appl. Sci., 2 (2015), pp. 206-220 https://www.sciencedirect.com/science/article/pii/S2314808X15000214/pdfft?md5=e82da4c639b347a349811b4ebead87fd&pid=1-s2.0-S2314808X15000214-main.pdf DOI: 10.1016/j.ejbas.2015.04.002

[36]
M.M.A. Khater, M.S. Mohamed, R.A.M. Attia.Chaos Solitons Fractals, 144 (2021), Article 110676

[37]
Y. Chu, M.M.A. Khater, Y.S. Hamed. AIP Adv., 11 (2021), Article 015223

[38]
M.M.A. Khater, A.E. Ahmed, M.A. El-Shorbagy.Results Phys., 22 (2021), Article 103890

[39]
M.M.A. Khater, A.E. Ahmed, S.H. Alfalqi, J.F. Alzaidi, S. Elbendary, A.M. Alabdali.Results Phys., 25 (2021), Article 104322

[40]
M.M.A. Khater, A.A. Mousa, M.A. El-Shorbagy, R.A.M. Attia.Results Phys., 22 (2021), Article 103954

[41]
M.M.A. Khater, K.S. Nisar, M.S. Mohamed. Math. Met. Appl. Sci., 44 (6) (2021), pp. 4598-4606 DOI: 10.1002/mma.7052

[42]
M.M.A. Khater, B. Ghanbari. Eur. Phys. J. Plus, 136 (447) (2021), pp. 1-28 DOI: 10. 3138/diaspora.21.1.2020-06-15

[43]
M.M.A. Khater, T.A. Nofal, H.A. Zinadah, M.S.M. Lotayife, D. Lu. Alex. Eng. J., 60 (2021), pp. 1797-1806

[44]
M.M.A. Khater, M.S. Mohamed, S.K. Elagan. Results Phys., 23 (2021), Article 104003

[45]
M.M.A. Khater, S. Anwar, K.U. Tariq, M.S. Mohamed. AIP Adv., 11 (2021), Article 025130

[46]
M.M.A. Khater, A. Bekir, D. Lu, R.A.M. Attia. Math. Met. Appl. Sci., 44 (3) (2021), pp. 2682-2691 DOI: 10.1002/mma.6951

[47]
M.M.A. Khater, A.A. Mousa, M.A. El-Shorbagy, R.A.M. Attia. Eur. Phys. J. Plus., 136 (604) ( 2021), pp. 1-11 DOI: 10. 3138/diaspora.21.1.2020-06-15

[48]
M.M.A. Khater, S.K. Elagan, A.A. Mousa, M.A. El-Shorbagy, S.H. Alfalqi, J.F. Alzaidi,D. Lu. Results Phys., 25 (2021), Article 104133

[49]
R.A.M. Attia, M.M.A. Khater, A.E. Ahmed, M.A. El-Shorbagy. AIP Adv., 11 (2021), Article 055105

[50]
M.M.A. Khater. Mod. Phys. Lett. B, 35 (22) ( 2021), Article 2150372 DOI: 10.1142/s0217984921503723

[51]
M.M.A. Khater, A.E. Ahmed, S.H. Alfalqi, J.F. Alzaidi. Phys. Scr., 96 (7) (2021), Article 075207

[52]
M.M.A. Khater, D. Lu. Mod. Phys. Lett. B, 35 (19) ( 2021), Article 2150324 DOI: 10.1142/s0217984921503243

[53]
M.M.A. Khater, A.E. Ahmed.AIMS Math., 6 (6) (2021), pp. 5896-5908 DOI: 10.3934/math.2021349

[54]
M.M.A. Khater, L. Akinyemi, S.K. Elagan, M.A. El-Shorbagy, S.H. Alfalqi, J.F. Alzaidi, N.A. Alshehri.Symmetry (Basel), 13 (2021), p. 963 DOI: 10.3390/sym13060963

[55]
C. Yue, D. Lu, M.M.A. Khater. Fluids, 6 (7) ( 2021), p. 235 DOI: 10.3390/fluids6070235

[56]
W. Li, L. Akinyemi, D. Lu, M.M.A. Khater.Symmetry (Basel), 13 (2021), p. 1085 DOI: 10.3390/sym13061085

[57]
M.M.A. Khater, A.M. Alabdali. Mathematics, 9 ( 2021), p. 1440 DOI: 10.3390/math9121440

[58]
P. Veeresha, D.G. Prakasha, N. Magesh, A.J. Christopher, D.U. Sarwe. J. Ocean Eng. Sci., 6 (3) (2021), pp. 265-275

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

[60]
A. Vithya, M.S.M. Rajan. J. Ocean Eng. Sci., 5 (3) (2020), pp. 205-213

[61]
M.M. Miah, A.R. Seadawy, H.M.S. Ali, M.A. Ali. J. Ocean Eng. Sci., 5 (3) (2020), pp. 269-278

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

[63]
M.P. Haghighi, J. Manafian. J. Ocean Eng. Sci., 5 (3) (2020), pp. 197-204

[64]
M.A. Shallal, K.K. Ali, K.R. Raslan, H. Rezazadeh, A. Bekir. J. Ocean Eng. Sci., 5 (3) (2020), pp. 223-229

[65]
M.T. Islam, M.A. Akter. Partial Differ. Eq. Appl. Math., 3 ( 2021), Article 100031

[66]
M.T. Islam, M.A. Akbar, M.A.K. Azad. Nonlinear Stud., 6 (2015), pp. 1-11

Options
Outlines

/