Fractional view evaluation system of Schrödinger-KdV equation by a comparative analysis

Rasool Shah; Abd-Allah Hyder; Naveed Iqbal; Thongchai Botmart; Rasool Shah; Abd-Allah Hyder; Naveed Iqbal; Thongchai Botmart

doi:10.3934/math.20221087

AIMS Mathematics

2022, Volume 7, Issue 11: 19846-19864. doi: 10.3934/math.20221087

Previous Article Next Article

Research article Special Issues

Fractional view evaluation system of Schrödinger-KdV equation by a comparative analysis

1.
Department of Mathematics, Abdul Wali khan university Mardan 23200, Pakistan
2.
Department of Mathematics, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
3.
Department of Engineering Mathematics and Physics, Faculty of Engineering, Al-Azhar University, Cairo, Egypt
4.
Department of Mathematics, Faculty of Science, University of Ha'il, Ha'il 81451, Saudi Arabia
5.
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand

Received: 15 June 2022 Revised: 23 August 2022 Accepted: 30 August 2022 Published: 08 September 2022
MSC : 34A34, 35A20, 35A22, 44A10, 33B15

The time-fractional coupled Schrödinger-KdV equation is an interesting mathematical model because of its wide and significant application in mathematics and applied sciences. A fractional coupled Schrödinger-KdV equation in the sense of Caputo derivative is investigated in this article. Namely, we provide a comparative study of the considered model using the Adomian decomposition method and the homotopy perturbation method with Shehu transform. Approximate solutions obtained using the Adomian decomposition and homotopy perturbation methods were numerically evaluated and presented in graphs and tables. Then, these solutions were compared to the exact solutions, demonstrating the simplicity, effectiveness, and good accuracy of the applied method. To demonstrate the accuracy and efficiency of the suggested techniques, numerical problem are provided.

Keywords:

Citation: Rasool Shah, Abd-Allah Hyder, Naveed Iqbal, Thongchai Botmart. Fractional view evaluation system of Schrödinger-KdV equation by a comparative analysis[J]. AIMS Mathematics, 2022, 7(11): 19846-19864. doi: 10.3934/math.20221087

Related Papers:

[1]	M. Hafiz Uddin, M. Ali Akbar, Md. Ashrafuzzaman Khan, Md. Abdul Haque . New exact solitary wave solutions to the space-time fractional differential equations with conformable derivative. AIMS Mathematics, 2019, 4(2): 199-214. doi: 10.3934/math.2019.2.199
[2]	Gulnur Yel, Haci Mehmet Baskonus, Wei Gao . New dark-bright soliton in the shallow water wave model. AIMS Mathematics, 2020, 5(4): 4027-4044. doi: 10.3934/math.2020259
[3]	Chander Bhan, Ravi Karwasra, Sandeep Malik, Sachin Kumar, Ahmed H. Arnous, Nehad Ali Shah, Jae Dong Chung . Bifurcation, chaotic behavior and soliton solutions to the KP-BBM equation through new Kudryashov and generalized Arnous methods. AIMS Mathematics, 2024, 9(4): 8749-8767. doi: 10.3934/math.2024424
[4]	Khudhayr A. Rashedi, Saima Noor, Tariq S. Alshammari, Imran Khan . Lump and kink soliton phenomena of Vakhnenko equation. AIMS Mathematics, 2024, 9(8): 21079-21093. doi: 10.3934/math.20241024
[5]	F. A. Mohammed . Soliton solutions for some nonlinear models in mathematical physics via conservation laws. AIMS Mathematics, 2022, 7(8): 15075-15093. doi: 10.3934/math.2022826
[6]	Ye Zhao, Chunfeng Xing . Orbital stability of periodic traveling waves to some coupled BBM equations. AIMS Mathematics, 2023, 8(9): 22225-22236. doi: 10.3934/math.20231133
[7]	Miao Yang, Lizhen Wang . Lie symmetry group, exact solutions and conservation laws for multi-term time fractional differential equations. AIMS Mathematics, 2023, 8(12): 30038-30058. doi: 10.3934/math.20231536
[8]	M. Mossa Al-Sawalha, Humaira Yasmin, Ali M. Mahnashi . Kink phenomena of the time-space fractional Oskolkov equation. AIMS Mathematics, 2024, 9(11): 31163-31179. doi: 10.3934/math.20241502
[9]	Nauman Raza, Maria Luz Gandarias, Ghada Ali Basendwah . Symmetry reductions and conservation laws of a modified-mixed KdV equation: exploring new interaction solutions. AIMS Mathematics, 2024, 9(4): 10289-10303. doi: 10.3934/math.2024503
[10]	Abdelkader Moumen, Khaled A. Aldwoah, Muntasir Suhail, Alwaleed Kamel, Hicham Saber, Manel Hleili, Sayed Saifullah . Investigation of more solitary waves solutions of the stochastics Benjamin-Bona-Mahony equation under beta operator. AIMS Mathematics, 2024, 9(10): 27403-27417. doi: 10.3934/math.20241331

Abstract

1. Introduction

Nonlinear partial differential equations (NPDEs) are crucial due to their ability to decode a wide-ranging of phenomena, as well as wave bending, fluid mechanics, hydrodynamics, organic molecular dispersion, magnetism, thermal conductivity, and many more. These phenomena are found across multiple scientific areas, such as mathematics, physics, engineering, biology, and finance. The broad range of NLPDEs is shown by a variety of equations, such as the Higgs system, advection equation, Boussinesq equation, Fisher's equation, and Burger's equation, among others ^[1,2,3,4,5].

Several asymptotic methods have been proposed in other studies to explore the internal dynamic behaviors of nonlinear partial differential equations (NPDEs) and fractional NPDEs with investigations on propagating solitons and other travelling wave solutions ^[6,7,8,9,10]. A soliton is an auto-oscillating wave-packet that preserves both its shape and speed during its propagation due to the fine-tuning between dispersion and nonlinearity. Kink waves, shock waves, lump waves, damped waves, periodic waves, etc are a few examples of solitons are present, which are described by soliton theory ^[11,12,13]. Despite the abundant availability of numerical solutions, analysts often resort to analytical methods because of the ability to explain the flow of physics and the accuracy of estimate of the behaviors of the system ^[14,15,16]. As a result, work on developing analytical treatments to investigate the solitonic behavior of NPDEs and fractional NPDEs is still ongoing, and several analytical methods have been created to investigate soliton solutions. $F$ -expansion technique ^[17], sech-tanh technique ^[18], Sardar sub-equation technique ^[19], $(G'/G)$ -expansion technique ^{[20,21,22,23]}, exp-function technique ^[24], sub-equation technique ^[25], tanh technique ^[26], $(G'/G^2)$ -expansion technique ^[27], Hirota bilinear technique ^[28], Kudryashov technique ^[29], Poincaré-Lighthill Kuo approach ^[30], unified technique ^[31], Riccati-Bernoulli Sub-ODE technique ^[32], extended direct algebraic technique ^{[33,34,35,36,37]}, auxiliary equation method ^[38], simple equation method ^[39], and RMESEM ^[40] are a few of these methods.

To discover cases of propagating solitons and other traveling waves, results for NPDEs and NFPDEs, this work establishes a new modification in a novel analytical technique called RMESEM with extended Riccati equation. This approach uses a variable-form wave transformation to transform the NPDE or NFPDE into an integer-order NODE. The resulting NODE is considered to have a series-based result (incorporating the solution of the Riccati equation). The substitution of the supposed solution in the resultant NODE transforms it into a set of algebraic equations. The soliton solutions for the relevant NPDEs and NFPDEs are acquired by more thoroughly finding the solutions of algebraic equations. The families of soliton solutions produced by RMESEM with the extended Riccati equation can help us comprehend the fundamental physical mechanisms and behavioral patterns of the nonlinear system.

To showcase the effectiveness of the improved RMESEM, the method is utilized to acquire soliton solutions for the (2+1)-dimensional GZK-BBME. In 1972, Benjamin et al. ^[41] investigated the issue of long waves with small and finite amplitudes and put out the Benjamin-Bona-Mahony equation (BBME), which took the following form:

$\begin{equation} v_t+\alpha v_x+\beta v v_x-\gamma v_{txx} = 0. \end{equation}$

(1.1)

To study weakly nonlinear ion-acoustic oscillations in low-pressure magnetized plasma, Zakharov and Kuznetsov ^[42] extended the Korteweg-de Vries (KdV) equation, resulting in the formulation of the Zakharov-Kuznetsov equation (ZKE) as:

$\begin{equation} v_t+\alpha v v_x+(v_{xx}+v_{yy}+v_{zz})_x = 0, \end{equation}$

(1.2)

where $v$ denotes ion velocity with the magnetic fields and is non-dimensional. Wazwaz ^[43] merged the ZKE and BBME equations in 2005 to create the (2+1) dimensional GZK-BBM equation ^[40]:

$\begin{equation} v_t+v_x+\alpha (v^\delta)_x+\beta(v_{xt}+v_{yy})_x = 0, \quad \delta > 1, \end{equation}$

(1.3)

where $v = v(t, x, y)$ , the coefficients $\alpha$ and $\beta$ are the relative dispersion and nonlinear parameters respectively, and $x$ and $y$ are the propagating & transverse coordinates. The nonlinear influence is introduced into the equation by the term $(v^\delta)_x$ . This formula is applied to the ZK-model analysis of long waves with finite amplitude. This equation's odd-order derivatives, $v_{xxt}$ and $v_{yyx}$ , take the dispersion effect into account. Examining the way the nonlinear and dispersion effects interact with the problem is one goal of the generalized GZK-BBME investigation. For instance, the tanh and sine-cosine approach was used to solve (3) ^[44]. By using a modified simple equation approach, Khan et al. ^[45] were able to get solitary wave solutions to the GZK-BBME for $n = 3$ . In 2015, Guner et al. ^[46] revealed the dark and bright soliton solution of the GZK-BBME. Finally, Patel and Kumar acquired numerical and semi-analytical solutions for GZK-BBME for $n = 2$ and $n = 3$ by the Adomian decomposition method and the variational iteration method, respectively, under some initial conditions ^[40]. The remaining task of the present analysis is to constructs and examine the propagation of the solitons in the context of GZK-BBME for $n = 3$ articulated in (2.1).

The rest of the study is structured as follows: In Section 2, we detail the working process of RMESEM. Following this approach, we develop soliton solutions for the GZK-BBME in Section 3. The depictions and graphical discussion are presented in Section 4. Lastly, Section 5 provides the conclusion in our research proposal.

2. The operational methodology of RMESEM

This section outlines the operational mechanism of RMESEM for constructing soliton solutions for NPDEs. Suppose the following general NPDE:

$\begin{equation} P( v , v_t, v_{x}, v_{y}, v v_{t}, \ldots) = 0, \end{equation}$

(2.1)

where $v = v(t, x, y)$ is an unknown function, $P$ is the polynomial of $v(t, x, y)$ , while the subscripts signify partial derivatives.

The main steps of the proposed RMESEM are as follows:

Step 1. Take into consideration the ensuing wave transformation

$\begin{equation} v(t, x, y) = V(\zeta), \quad \text{where} \quad \zeta = x+y- \omega t, \end{equation}$

(2.2)

where $\omega$ denotes wave speed. Equation (2.1) is converted into the following NODE using the above wave transformation:

$\begin{equation} Q(V, V'V, V', \dots) = 0, \end{equation}$

(2.3)

where the primes indicate the ordinary derivatives of $V$ with respect to $\zeta$ , and $Q$ is a polynomial of $V$ and its derivatives.

Step 2. Equation (2.3) is sometimes integrated term by term to be made conformable for the homogeneous balancing rule.

Step 3. Following that, we suppose that a closed-form wave solution for Eq (2.3) can be expressed in the ensuing form:

$\begin{equation} V(\zeta) = \sum\limits_{j = 0}^{\gamma}k_j \left(\frac{\Phi'(\zeta)}{\Phi(\zeta)}\right)^j+ \sum\limits_{\epsilon = 0}^{\gamma-1}s_\epsilon \left(\frac{\Phi'(\zeta)}{\Phi(\zeta)}\right)^\epsilon \cdot \left(\frac{1}{\Phi(\zeta)} \right), \end{equation}$

(2.4)

where $k_j (j = 0, ..., \gamma)$ and $s_\epsilon (\epsilon = 0, ..., \gamma-1)$ represent the unknown constants that need to be determined later and $\Phi(\zeta)$ satisfies the subsequent 1st order Riccati equation:

$\begin{equation} \Phi'(\zeta) = p+q \Phi(\zeta)+r(\Phi(\zeta))^2, \end{equation}$

(2.5)

where $p, q$ and $r$ are constants.

Step 4. To calculate the integer $\gamma$ presented in Eq (2.4), we take the homogeneous balance between the highest nonlinear term and the highest-order derivative term in Eq (2.3).

Step 5. Substituting the value of $\gamma$ got in step 4 into Eq (2.4) and substituting the result together with $\Phi(\zeta)$ raised to the exponent equal to the index of the integral on the left-hand side of Eq (2.3) or substituting the result with $\Phi(\zeta)$ when we have integrated Eq (2.3) gives an expression in terms of $\Phi(\zeta)$ . New additional construction by comparison of coefficient expression gives a system of algebraic equations of $k_j (j = 0, ..., \gamma)$ and $s_\epsilon (\epsilon = 0, ..., \gamma-1)$ with other associated parameters that are introduced.

Step 6. Solving the algebraic system obtained in Step 5 with the help of algebraic software Maple yields values of $k_j (j = 0, ..., \gamma)$ and $s_\epsilon (\epsilon = 0, ..., \gamma-1)$ with additional associated parameters.

Step 7. Finally, soliton solutions to Eq (2.1) are derived by determining and substituting the calculated values of parameters in Eq (2.4) with the solutions of Eq (2.5) that are given in Table 1.

Table 1. The solutions

$\Phi(\zeta)$ that meet the particular Riccati equation in (2.4) and the structure of

$\bigg(\frac{\Phi'(\zeta)}{\Phi(\zeta)}\bigg)$ , where

$\kappa = q^2-4r p$ and

$\wp = \cosh \left(\frac{1}{4}\, \sqrt {\kappa}\zeta \right) \sinh \left(\frac{1}{4}\, \sqrt {\kappa}\zeta \right)$ .

S. No.	Family	Condition(s)	$\Phi(\zeta)$	$\bigg(\frac{\Phi'(\zeta)}{\Phi(\zeta)}\bigg)$
1	Trigonometric Solutions	$\kappa < 0, \quad r\neq0$
	Trigonometric Solutions		$-{\frac {q}{2r}}+{\frac {\sqrt {-\kappa}\tan \left(\frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{2r}}$ ,	$-\frac{1}{2}\, {\frac {\kappa\, \left(1+ \left(\tan \left(\frac{1}{2}\, \sqrt {-\kappa}\zeta \right) \right) ^{2} \right) }{-q+\sqrt {-\kappa}\tan \left(\frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }}$ ,
			$-{\frac {q}{2r}}-{\frac {\sqrt {-\kappa}\cot \left(\frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{2r}}$ ,	$\frac{1}{2}\, {\frac { \left(1+ \left(\cot \left(\frac{1}{2}\, \sqrt {-\kappa}\zeta \right) \right) ^{2} \right) \kappa}{q+\sqrt {-\kappa}\cot \left(\frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }}$ ,
			$-{\frac {q}{2r}}+{\frac {\sqrt {-\kappa} \left(\tan \left(\sqrt{-\kappa}\zeta \right) + \left(\sec \left(\sqrt {-\kappa}\zeta \right) \right) \right) }{2r}}$ ,	$-{\frac {\kappa\, \left(1+\sin \left(\sqrt {-\kappa}\zeta \right) \right) \sec \left(\sqrt {-\kappa}\zeta \right)}{ -q\cos \left(\sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa}\sin \left(\sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa} }}$ ,
			$-{\frac {q}{2r}}+{\frac {\sqrt {-\kappa} \left(\tan \left(\sqrt{-\kappa}\zeta \right) - \left(\sec \left(\sqrt {-\kappa}\zeta \right) \right) \right) }{2r}}$ .	${\frac {\kappa\, \left(\sin \left(\sqrt {-\kappa}\zeta \right) -1 \right) \sec \left(\sqrt {-\kappa}\zeta \right)}{ -q\cos \left(\sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa}\sin \left(\sqrt {-\kappa}\zeta \right) -\sqrt {-\kappa} }}$ .
2	Rational Solutions
	Rational Solutions	$\kappa=0$	$-2\, {\frac {p \left(q\zeta\, +2 \right) }{{q}^{2}\zeta\, }}$ ,	$-2\, {\frac {1}{\zeta\, \left(q\zeta+2 \right) }}$ ,
		$\kappa=0$ , & $q=r=0$	$\zeta\, p$ ,	$\frac {1}{\zeta\, }$ ,
		$\kappa=0$ , & $q=p=0$	$-{\frac {1}{\zeta\, r}}$ .	$-\frac {1}{\zeta\, }$ .
3	Hyperbolic Solutions	$\kappa > 0, \quad r\neq0$
	Hyperbolic Solutions		$-{\frac {q}{2r}}-{\frac {\sqrt {\kappa}\tanh \left(\frac{1}{2}\, \sqrt {\kappa}\zeta \right) }{2r}}$ ,	$-\frac{1}{2}\, {\frac { \left(-1+ \left(\tanh \left(\frac{1}{2}\, \sqrt {\kappa}\zeta \right) \right) ^{2} \right) \kappa}{q+\sqrt {\kappa}\tanh \left(\frac{1}{2}\, \sqrt {\kappa}\zeta \right) }}$ ,
			$-{\frac {q}{2r}}-{\frac {\sqrt {\kappa} \left(\tanh \left(\sqrt{\kappa}\zeta \right) + i \left({\rm sech} \left(\sqrt {\kappa}\zeta \right) \right) \right) }{2r}}$ ,	$-{\frac {\kappa\, \left(-1+i\sinh \left(\sqrt {\kappa}\zeta \right) \right) }{\cosh \left(\sqrt {\kappa}\zeta \right) \left(q\cosh \left(\sqrt {\kappa}\zeta \right) +\sqrt {\kappa}\sinh \left(\sqrt { \kappa}\zeta \right) +i\sqrt {\kappa} \right) }}$ ,
			$-{\frac {q}{2r}}-{\frac {\sqrt {\kappa} \left(\tanh \left(\sqrt{\kappa}\zeta \right) - i \left({\rm sech} \left(\sqrt {\kappa}\zeta \right) \right) \right) }{2r}}$ ,	$-{\frac {\kappa\, \left(1+i\sinh \left(\sqrt {\kappa}\zeta \right) \right) }{\cosh \left(\sqrt {\kappa}\zeta \right) \left(-q\cosh \left(\sqrt {\kappa}\zeta \right) -\sqrt {\kappa}\sinh \left(\sqrt { \kappa}\zeta \right) +i\sqrt {\kappa} \right) }},$
			$-{\frac {q}{2r}}-{\frac {\sqrt {\kappa} \left(\coth \left(\sqrt{\kappa}\zeta \right) + \left({\rm csch} \left(\sqrt {\kappa}\zeta \right) \right) \right) }{2r}}$ .	$-\frac{1}{4}\, {\frac {\kappa\, \left(2\, \left(\cosh \left(\frac{1}{4}\, \sqrt { \kappa}\zeta \right) \right) ^{2}-1 \right) }{\wp \left(-2\, q\wp +\sqrt {\kappa} \right) }}.$
4	Rational-Hyperbolic Solutions
		$p=0$ , & $q\neq0$ , $r\neq 0$	$-{\frac {\lambda\, b}{r \left(\cosh \left(b \zeta \right) -\sinh \left(b \zeta \right) +\lambda \right) }}$ ,	${\frac {q \left(\sinh \left(q\zeta \right) -\cosh \left(q\zeta \right) \right) }{-\cosh \left(q\zeta \right) +\sinh \left(q\zeta \right) -\lambda}}$ ,
			$-{\frac {b \left(\cosh \left(b \zeta \right) +\sinh \left(b \zeta \right) \right) }{r \left(\cosh \left(b \zeta \right) +\sinh \left(b \zeta \right) +\mu \right) }}$ .	${\frac {q\mu}{\cosh \left(q\zeta \right) +\sinh \left(q\zeta \right) + \mu}}$ .
5	Exponential Solutions
	Exponential Solutions	$r=0$ , & $q=\Upsilon$ , $p=h\Upsilon$	${e}^{\Upsilon\, \zeta}-h$ ,	${\frac {\Upsilon\, {{\rm e}^{\Upsilon\, \zeta}}}{{{\rm e}^{\Upsilon\, \zeta}}-h}}$ ,
		$p=0$ , & $q=\Upsilon$ , $r=h\Upsilon$	${\frac {{e}^{\Upsilon\, \zeta}}{1-h{e}^{\Upsilon \zeta}}}$ .	$-{\frac {\Upsilon}{-1+h{{\rm e}^{\Upsilon\, \zeta}}}}$ .

| Show Table

DownLoad: CSV

3. Execution of RMESEM

This section employs the proposed RMESEM for the establishment of new plethora of soliton solutions for GZK-BBME with $\delta = 3$ of the form:

$\begin{equation} v_t+v_x+\alpha (v^3)_x+\beta(v_{xt}+v_{yy})_x = 0. \end{equation}$

(3.1)

We proceed by performing the wave transformation given in Eq (2.2), which transforms Eq (3.1) into the ensuing NODE:

$\begin{equation} -\omega V'+V'+\alpha (V^3)'+\beta(-\omega V''+V'')' = 0. \end{equation}$

(3.2)

Upon integrating Eq (3.2) with zero constant of integration, we obtain the following NODE:

$\begin{equation} (1-\omega) V+\alpha V^3+ (1-\omega) \beta V''. \end{equation}$

(3.3)

Establishing the principle of homogenous balance between terms $V''$ and $V^3$ in Eq (3.3) suggests that $2 +\gamma = 3 \gamma$ which applies $\gamma = 1$ . Substituting $\gamma = 1$ in Eq (2.4) presents the series form closed solutions for Eq (3.3):

$\begin{equation} V(\zeta) = \sum\limits_{j = 0}^{1}k_j \left(\frac{\Phi'(\zeta)}{\Phi(\zeta)}\right)^j+ \left(\frac{s_0 }{\Phi(\zeta)} \right). \end{equation}$

(3.4)

An expression in $\Phi(\zeta)$ is generated by entering Eq (3.4) into Eq (3.3) and gathering all terms with the equal exponents of $\Phi(\zeta)$ . The achieved expressions can be reduced to the ensuing scheme of seven nonlinear algebraic equations by putting the coefficients to zero:

$\begin{equation*} \begin{split} 2\, \beta\, k_{{1}}{r}^{3}-2\, \beta\, \omega\, k_{{1}}{r}^{3}+\alpha\, {k_{{1} }}^{3}{r}^{3} = 0, \end{split} \end{equation*}$

$\begin{equation*} \begin{split} 3\, \alpha\, {k_{{1}}}^{3}q{r}^{2}-3\, \beta\, \omega\, k_{{1}}q{r}^{2}+3\, \alpha\, k_{{0}}{k_{{1}}}^{2}{r}^{2}+3\, \beta\, k_{{1}}q{r}^{2} = 0, \end{split} \end{equation*}$

$\begin{equation*} \begin{split} & -2\, \beta\, \omega\, k_{{1}}{r}^{2}p+2\, \beta\, k_{{1}}{r}^{ 2}p+3\, \alpha\, {k_{{0}}}^{2}k_{{1}}r-\beta\, \omega\, k_{{1}}{q}^{2}r+3\, \alpha\, {k_{{1}}}^{3}p{r}^{2}\\ &-\omega\, k_{{1}}r+k_{{1}}r+3\, \alpha\, {k_{{1}}}^{3}{q}^{ 2}r+6\, \alpha\, k_{{0}}{k_{{1}}}^{2}qr+\beta\, k_{{1}}{q}^{2}r+3\, \alpha \, {k_{{1}}}^{2}{r}^{2}s_{{0}} = 0, \end{split} \end{equation*}$

$\begin{equation*} \begin{split} & 6\, \alpha\, {k_{{1}}}^{3}pqr+\alpha\, {k_{{0}}}^{3}+2\, \beta\, k_{{1}}qpr +6\, \alpha\, k_{{0}}k_{{1}}rs_{{0}}+k_{{1}}q-\omega\, k_{{1}}q+6\, \alpha\, k_{{0}}{k_{{1}}}^{2}pr+3\, \alpha\, {k_{{0}}}^{2}k_{{1}}q\\ &+6\, \alpha\, {k_{{1}}}^{2}qrs_{{0}}+k_{{0}}+\beta\, s_{{0}}qr-\beta\, \omega\, s_ {{0}}qr+\alpha\, {k_{{1}}}^{3}{q}^{3}+3\, \alpha\, k_{{0}}{k_{{1}}}^{2} {q}^{2}-\omega\, k_{{0}}-2\, \beta\, \omega\, k_{{1}}qpr = 0, \end{split} \end{equation*}$

$\begin{equation*} \begin{split} & -\omega\, s_{{0}}+\beta\, s_{{0}}{q}^{2}+6\, \alpha\, {k_{{1}}}^{2}prs_{{0} }+6\, \alpha\, k_{{0}}{k_{{1}}}^{2}pq+3\, \alpha\, {k_{{1}}}^{3}{p}^{2}r -\omega\, k_{{1}}p+k_{{1}}p+3\, \alpha\, {k_{{0}}}^{2}s_{{0}}\\ &-2\, \beta\, \omega\, k_{{1}}r{p}^{2}+s_{{0}}-2\, \beta\, \omega\, s_{{0}}pr+6\, \alpha\, k_{{0}}k_{{1}}qs_{{0}}+3\, \alpha\, {k_{{0}}}^{2}k_{{1}}p+2\, \beta\, s_{{0 }}pr+\beta\, k_{{1}}{q}^{2}p\\ &+3\, \alpha\, k_{{1}}r{s_{{0}}}^{2}+3\, \alpha \, {k_{{1}}}^{3}p{q}^{2}-\beta\, \omega\, s_{{0}}{q}^{2}+3\, \alpha\, {k_{{1 }}}^{2}{q}^{2}s_{{0}}+2\, \beta\, k_{{1}}r{p}^{2}-\beta\, \omega\, k_{{1}}{q}^ {2}p = 0, \end{split} \end{equation*}$

$\begin{equation*} \begin{split} & 3\, \alpha\, {k_{{1}}}^{3}{p}^{2}q+3\, \beta\, s_{{0}}pq+3\, \alpha\, k_{{1} }q{s_{{0}}}^{2}-3\, \beta\, \omega\, s_{{0}}pq+6\, \alpha\, k_{{0}}k_{{1}}ps _{{0}}\\ &-3\, \beta\, \omega\, k_{{1}}q{p}^{2}+6\, \alpha\, {k_{{1}}}^{2}pqs_{{0 }}+3\, \alpha\, k_{{0}}{k_{{1}}}^{2}{p}^{2}+3\, \beta\, k_{{1}}q{p}^{2}+3\, \alpha\, k_{{0}}{s_{{0}}}^{2} = 0, \end{split} \end{equation*}$

and

$\begin{equation*} \begin{split} -2\, \beta\, \omega\, s_{{0}}{p}^{2}+\alpha\, {k_{{1}}}^{3}{p}^{3}+3\, \alpha\, {k_{{1}}}^{2}{p}^{2}s_{{0}}+\alpha\, {s_{{0}}}^{3}-2\, \beta\, \omega\, k_{{1}}{p}^{3}+2\, \beta\, s_{{0}}{p}^{2}+3\, \alpha\, k_{{1}}p{s_{ {0}}}^{2}+2\, \beta\, k_{{1}}{p}^{3} = 0. \end{split} \end{equation*}$

When tackling the result scheme with Maple, the followings three sorts of results become available:

Case 1.

$\begin{equation} \begin{split} k_{{1}} = 0, \; \; s_{{0}} = s_{{0}}, \; \; k_{{0}} = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}, \; \; \beta = \frac{2\, }{ \kappa}, \; \; \alpha = 4\, {\frac {{p}^{2} \left( \omega-1 \right) }{ \kappa\, {s_{{0}}}^{2}}}, \; \; \omega = \omega. \end{split} \end{equation}$

(3.5)

Case 2.

$\begin{equation} \begin{split} k_{{1}} = k_{{1}}, \; \; s_{{0}} = s_{{0}}, \; \; k_{{0}} = k_{{0}}, \; \; \beta = \beta, \; \; \alpha = 0, \; \; \omega = 1. \end{split} \end{equation}$

(3.6)

Case 3.

$\begin{equation} \begin{split} k_{{1}} = k_{{1}}, \; \; s_{{0}} = -k_{{1}}p, \; \; k_{{0}} = -\frac{1}{2}\, k_{{1}}q, \; \; \beta = \frac{2\, }{ \kappa}, \; \; \alpha = 4\, {\frac {\omega-1}{\kappa\, {k_{{1}}}^{2}}}, \; \; \omega = \omega. \end{split} \end{equation}$

(3.7)

When Case 1 is assumed and Eqs (2.2) and (3.4), together with the consistent general result of Eq (2.5) given in Table 1. The following cases of soliton results for GZK-BBME expressed in Eq (3.1) result:

Set. 1.1. With $\kappa < 0, \; r\neq0$ ,

$\begin{equation} \begin{split} v_{1, 1}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, { \frac {\sqrt {-\kappa}\tan \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{r }} \right) ^{-1} , \end{split} \end{equation}$

(3.8)

$\begin{equation} \begin{split} v_{1, 2}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, { \frac {\sqrt {-\kappa}\cot \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{r }} \right) ^{-1} , \end{split} \end{equation}$

(3.9)

$\begin{equation} \begin{split} v_{1, 3}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, { \frac {\sqrt {-\kappa} \left( \tan \left( \sqrt {-\kappa}\zeta \right) +\sec \left( \sqrt {-\kappa}\zeta \right) \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.10)

and

$\begin{equation} \begin{split} v_{1, 4}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, { \frac {\sqrt {-\kappa} \left( \tan \left( \sqrt {-\kappa}\zeta \right) -\sec \left( \sqrt {-\kappa}\zeta \right) \right) }{r}} \right) ^{-1} . \end{split} \end{equation}$

(3.11)

Set. 1.2. With $\kappa > 0, \; r\neq0$ ,

$\begin{equation} \begin{split} v_{1, 5}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, { \frac {\sqrt {\kappa}\tanh \left( \frac{1}{2}\, \sqrt {\kappa}\zeta \right) }{r} } \right) ^{-1} , \end{split} \end{equation}$

(3.12)

$\begin{equation} \begin{split} v_{1, 6}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, { \frac {\sqrt {\kappa} \left( \tanh \left( \sqrt {\kappa}\zeta \right) + i{\rm sech} \left( \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^ {-1} , \end{split} \end{equation}$

(3.13)

$\begin{equation} \begin{split} v_{1, 7}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, { \frac {\sqrt {\kappa} \left( \tanh \left( \sqrt {\kappa}\zeta \right) - i{\rm sech} \left( \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^ {-1} , \end{split} \end{equation}$

(3.14)

and

$\begin{equation} \begin{split} v_{1, 8}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}q}{p}}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{4}\, { \frac {\sqrt {\kappa} \left( \tanh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) -\coth \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^{-1} . \end{split} \end{equation}$

(3.15)

Set. 1.3. With $q = \Upsilon$ , $p = h\Upsilon(h\neq0)$ and $r = 0$ ,

$\begin{equation} \begin{split} v_{1, 9}(x, y, t) = \frac{1}{2}\, {\frac {s_{{0}}}{h}}+{\frac {s_{{0}}}{{{\rm e}^{\Upsilon\, \zeta}}-h}} . \end{split} \end{equation}$

(3.16)

In above solutions $\zeta = x+y-\omega t$ .

When Case 2 is assumed and Eqs (2.2) and (3.4) together with the consistent general result of Eq (2.5) given in Table 1. The following cases of soliton results for GZK-BBME expressed in Eq (3.1) result:

Set. 2.1. With $\kappa < 0, \; r\neq0$ ,

$\begin{equation} \begin{split} v_{2, 1}(x, y, t) = k_{{0}}-\frac{1}{2}\, {\frac {k_{{1}}\kappa\, \left( 1+ \left( \tan \left( \frac{1}{2} \, \sqrt {-\kappa}\zeta \right) \right) ^{2} \right) }{-q+\sqrt {- \kappa}\tan \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, {\frac {\sqrt {-\kappa}\tan \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.17)

$\begin{equation} \begin{split} v_{2, 2}(x, y, t) = k_{{0}}+\frac{1}{2}\, {\frac {k_{{1}}\kappa\, \left( 1+ \left( \cot \left( \frac{1}{2} \, \sqrt {-\kappa}\zeta \right) \right) ^{2} \right) }{q+\sqrt {-\kappa }\cot \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }}+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt {-\kappa}\cot \left( \frac{1}{2}\, \sqrt {- \kappa}\zeta \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.18)

$\begin{equation} \begin{split} v_{2, 3}(x, y, t) = & k_{{0}}-{\frac {k_{{1}}\kappa\, \left( 1+\sin \left( \sqrt {-\kappa} \zeta \right) \right) }{\cos \left( \sqrt {-\kappa}\zeta \right) \left( -q\cos \left( \sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa}\sin \left( \sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa} \right) }}\\ &+s_{{0} } \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, {\frac {\sqrt {-\kappa} \left( \tan \left( \sqrt {-\kappa}\zeta \right) +\sec \left( \sqrt {-\kappa}\zeta \right) \right) }{r}} \right) ^{-1}+ k_{{0}} , \end{split} \end{equation}$

(3.19)

and

$\begin{equation} \begin{split} v_{2, 4}(x, y, t) = & {\frac {k_{{1}}\kappa\, \left( \sin \left( \sqrt {-\kappa}\zeta \right) -1 \right) }{\cos \left( \sqrt {-\kappa}\zeta \right) \left( -q\cos \left( \sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa}\sin \left( \sqrt {-\kappa}\zeta \right) -\sqrt {-\kappa} \right) }}\\ &+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, {\frac {\sqrt {-\kappa} \left( \tan \left( \sqrt {-\kappa}\zeta \right) -\sec \left( \sqrt {-\kappa}\zeta \right) \right) }{r}} \right) ^{-1} + k_{{0}} . \end{split} \end{equation}$

(3.20)

Set. 2.2. With $\kappa > 0, \; r\neq0$ ,

$\begin{equation} \begin{split} v_{2, 5}(x, y, t) = k_{{0}}-\frac{1}{2}\, {\frac {k_{{1}}\kappa\, \left( -1+ \left( \tanh \left( 1/ 2\, \sqrt {\kappa}\zeta \right) \right) ^{2} \right) }{q+\sqrt {\kappa} \tanh \left( \frac{1}{2}\, \sqrt {\kappa}\zeta \right) }}+s_{{0}} \left( -\frac{1}{2}\, { \frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt {\kappa}\tanh \left( \frac{1}{2}\, \sqrt { \kappa}\zeta \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.21)

$\begin{equation} \begin{split} v_{2, 6}(x, y, t) = & -{\frac {k_{{1}}\kappa\, \left( -1+i\sinh \left( \sqrt {\kappa} \zeta \right) \right) }{\cosh \left( \sqrt {\kappa}\zeta \right) \left( q\cosh \left( \sqrt {\kappa}\zeta \right) +\sqrt {\kappa}\sinh \left( \sqrt {\kappa}\zeta \right) +i\sqrt {\kappa} \right) }}\\ &+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt {\kappa} \left( \tanh \left( \sqrt {\kappa}\zeta \right) +i{\rm sech} \left( \sqrt {\kappa} \zeta \right) \right) }{r}} \right) ^{-1} + k_{{0}} , \end{split} \end{equation}$

(3.22)

$\begin{equation} \begin{split} v_{2, 7}(x, y, t) = & -{\frac {k_{{1}}\kappa\, \left( 1+i\sinh \left( \sqrt {\kappa} \zeta \right) \right) }{\cosh \left( \sqrt {\kappa}\zeta \right) \left( -q\cosh \left( \sqrt {\kappa}\zeta \right) -\sqrt {\kappa}\sinh \left( \sqrt {\kappa}\zeta \right) +i\sqrt {\kappa} \right) }}\\ &+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt {\kappa} \left( \tanh \left( \sqrt {\kappa}\zeta \right) -i{\rm sech} \left( \sqrt {\kappa} \zeta \right) \right) }{r}} \right) ^{-1} + k_{{0}} , \end{split} \end{equation}$

(3.23)

and

$\begin{equation} \begin{split} v_{2, 8}(x, y, t) = & -\frac{1}{4}\, {\frac {k_{{1}}\kappa\, \left( 2\, \left( \cosh \left( \frac{1}{4}\sqrt {\kappa}\zeta \right) \right) ^{2}-1 \right) }{\cosh \left( \frac{1}{4}\sqrt {\kappa}\zeta \right) \sinh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \left( -2\, q\cosh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \sinh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) +\sqrt {\kappa} \right) } }\\ &+s_{{0}} \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{4}\, {\frac {\sqrt {\kappa} \left( \tanh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) -\coth \left( \frac{1}{4} \, \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^{-1}+ k_{{0}} . \end{split} \end{equation}$

(3.24)

Set. 2.3. With $\kappa = 0, \quad q\neq0$ ,

$\begin{equation} \begin{split} v_{2, 9}(x, y, t) = k_{{0}}-2\, {\frac {k_{{1}}}{\zeta\, \left( q\zeta+2 \right) }}-\frac{1}{2}\, { \frac {s_{{0}}{q}^{2}\zeta}{p \left( q\zeta+2 \right) }} . \end{split} \end{equation}$

(3.25)

Set. 2.4. With $\kappa = 0$ , in case when $q = r = 0$ ,

$\begin{equation} \begin{split} v_{2, 10}(x, y, t) = k_{{0}}+{\frac {k_{{1}}}{\zeta}}+{\frac {s_{{0}}}{p\zeta}} . \end{split} \end{equation}$

(3.26)

Set. 2.5. With $\kappa = 0$ , in case when $q = p = 0$ ,

$\begin{equation} \begin{split} v_{2, 11}(x, y, t) = k_{{0}}-{\frac {k_{{1}}}{\zeta}}-s_{{0}}r\zeta . \end{split} \end{equation}$

(3.27)

Set. 2.6. With $q = \Upsilon$ , $p = h\Upsilon(h\neq0)$ and $r = 0$ ,

$\begin{equation} \begin{split} v_{2, 12}(x, y, t) = k_{{0}}+{\frac {k_{{1}}\Upsilon\, {{\rm e}^{\Upsilon\, \zeta}}}{{{\rm e}^{\Upsilon\, \zeta}}-h}}+{\frac {s_{{0}}}{{{\rm e}^{\Upsilon\, \zeta}}-h}} . \end{split} \end{equation}$

(3.28)

Set. 2.7. With $q = \Upsilon$ , $r = h\Upsilon(h\neq0)$ and $p = 0$ ,

$\begin{equation} \begin{split} v_{2, 13}(x, y, t) = k_{{0}}-{\frac {k_{{1}}\Upsilon}{-1+h{{\rm e}^{\Upsilon\, \zeta}}}}+{\frac {s_{{0 }} \left( 1-h{{\rm e}^{\Upsilon\, \zeta}} \right) }{{{\rm e}^{\Upsilon\, \zeta}}}} . \end{split} \end{equation}$

(3.29)

Set. 2.8. With $p = 0$ , $r\neq0$ and $q\neq0$ ,

$\begin{equation} \begin{split} v_{2, 14}(x, y, t) = k_{{0}}+{\frac {k_{{1}}q \left( \sinh \left( q\zeta \right) -\cosh \left( q\zeta \right) \right) }{-\cosh \left( q\zeta \right) +\sinh \left( q\zeta \right) -\lambda}}-{\frac {s_{{0}}r \left( \cosh \left( q\zeta \right) -\sinh \left( q\zeta \right) +\lambda \right) }{\lambda\, q}} , \end{split} \end{equation}$

(3.30)

and

$\begin{equation} \begin{split} v_{2, 15}(x, y, t) = k_{{0}}+{\frac {k_{{1}}q\mu}{\cosh \left( q\zeta \right) +\sinh \left( q\zeta \right) +\mu}}-{\frac {s_{{0}}r \left( \cosh \left( q\zeta \right) +\sinh \left( q\zeta \right) +\mu \right) }{q \left( \cosh \left( q\zeta \right) +\sinh \left( q\zeta \right) \right) }} . \end{split} \end{equation}$

(3.31)

In above solutions $\zeta = x+y-t$ .

When Case 3 is assumed and Eqs (2.2) and (3.4) together with the consistent general result of Eq (2.5) given in Table 1. The following cases of soliton results for GZK-BBME expressed in Eq (3.1) result:

Set. 3.1. With $\kappa < 0, \; r\neq0$ ,

$\begin{equation} \begin{split} v_{3, 1}(x, y, t) = -\frac{1}{2}\, k_{{1}}q-\frac{1}{2}\, {\frac {k_{{1}}\kappa\, \left( 1+ \left( \tan \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) \right) ^{2} \right) }{-q+ \sqrt {-\kappa}\tan \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }}-k_{{1}} p \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, {\frac {\sqrt {-\kappa}\tan \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.32)

$\begin{equation} \begin{split} v_{3, 2}(x, y, t) = -\frac{1}{2}\, k_{{1}}q+\frac{1}{2}\, {\frac {k_{{1}}\kappa\, \left( 1+ \left( \cot \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) \right) ^{2} \right) }{q+ \sqrt {-\kappa}\cot \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }}-k_{{1}} p \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt {-\kappa}\cot \left( \frac{1}{2}\, \sqrt {-\kappa}\zeta \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.33)

$\begin{equation} \begin{split} v_{3, 3}(x, y, t) = & -{\frac {k_{{1}}\kappa\, \left( 1+\sin \left( \sqrt {- \kappa}\zeta \right) \right) }{\cos \left( \sqrt {-\kappa}\zeta \right) \left( -q\cos \left( \sqrt {-\kappa}\zeta \right) +\sqrt {- \kappa}\sin \left( \sqrt {-\kappa}\zeta \right) +\sqrt {-\kappa} \right) }}\\ &-k_{{1}}p \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, {\frac {\sqrt {- \kappa} \left( \tan \left( \sqrt {-\kappa}\zeta \right) +\sec \left( \sqrt {-\kappa}\zeta \right) \right) }{r}} \right) ^{-1}-\frac{1}{2}\, k_{{1}}q , \end{split} \end{equation}$

(3.34)

and

$\begin{equation} \begin{split} v_{3, 4}(x, y, t) = & {\frac {k_{{1}}\kappa\, \left( \sin \left( \sqrt {- \kappa}\zeta \right) -1 \right) }{\cos \left( \sqrt {-\kappa}\zeta \right) \left( -q\cos \left( \sqrt {-\kappa}\zeta \right) +\sqrt {- \kappa}\sin \left( \sqrt {-\kappa}\zeta \right) -\sqrt {-\kappa} \right) }}\\ &-k_{{1}}p \left( -\frac{1}{2}\, {\frac {q}{r}}+\frac{1}{2}\, {\frac {\sqrt {- \kappa} \left( \tan \left( \sqrt {-\kappa}\zeta \right) -\sec \left( \sqrt {-\kappa}\zeta \right) \right) }{r}} \right) ^{-1}-\frac{1}{2}\, k_{{1}}q . \end{split} \end{equation}$

(3.35)

Set. 3.2. With $\kappa > 0, \; r\neq0$ ,

$\begin{equation} \begin{split} v_{3, 5}(x, y, t) = -\frac{1}{2}\, k_{{1}}q-\frac{1}{2}\, {\frac {k_{{1}}\kappa\, \left( -1+ \left( \tanh \left( \frac{1}{2}\, \sqrt {\kappa}\zeta \right) \right) ^{2} \right) }{q+ \sqrt {\kappa}\tanh \left( \frac{1}{2}\, \sqrt {\kappa}\zeta \right) }}-k_{{1}}p \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt {\kappa}\tanh \left( \frac{1}{2}\sqrt {\kappa}\zeta \right) }{r}} \right) ^{-1} , \end{split} \end{equation}$

(3.36)

$\begin{equation} \begin{split} v_{3, 6}(x, y, t) = &-{\frac {k_{{1}}\kappa\, \left( -1+i\sinh \left( \sqrt { \kappa}\zeta \right) \right) }{\cosh \left( \sqrt {\kappa}\zeta \right) \left( q\cosh \left( \sqrt {\kappa}\zeta \right) +\sqrt { \kappa}\sinh \left( \sqrt {\kappa}\zeta \right) +i\sqrt {\kappa} \right) }}\\ &-k_{{1}}p \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt { \kappa} \left( \tanh \left( \sqrt {\kappa}\zeta \right) +i{\rm sech} \left( \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^{-1}-\frac{1}{2}\, k_{{1}}q , \end{split} \end{equation}$

(3.37)

$\begin{equation} \begin{split} v_{3, 7}(x, y, t) = & -{\frac {k_{{1}}\kappa\, \left( 1+i\sinh \left( \sqrt { \kappa}\zeta \right) \right) }{\cosh \left( \sqrt {\kappa}\zeta \right) \left( -q\cosh \left( \sqrt {\kappa}\zeta \right) -\sqrt { \kappa}\sinh \left( \sqrt {\kappa}\zeta \right) +i\sqrt {\kappa} \right) }}\\ &-k_{{1}}p \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{2}\, {\frac {\sqrt { \kappa} \left( \tanh \left( \sqrt {\kappa}\zeta \right) -i{\rm sech} \left( \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^{-1}-\frac{1}{2}\, k_{{1}}q , \end{split} \end{equation}$

(3.38)

and

$\begin{equation} \begin{split} v_{3, 8}(x, y, t) = &-\frac{1}{4}\, {\frac {k_{{1}}\kappa\, \left( 2\, \left( \cosh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \right) ^{2}-1 \right) }{ \cosh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \sinh \left( \frac{1}{4}\, \sqrt { \kappa}\zeta \right) \left( -2\, q\cosh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \sinh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) +\sqrt {\kappa} \right) }}\\ &-k_{{1}}p \left( -\frac{1}{2}\, {\frac {q}{r}}-\frac{1}{4}\, {\frac {\sqrt { \kappa} \left( \tanh \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) -\coth \left( \frac{1}{4}\, \sqrt {\kappa}\zeta \right) \right) }{r}} \right) ^{-1} -\frac{1}{2}\, k_{{1}}q . \end{split} \end{equation}$

(3.39)

Set. 3.3. With $q = \Upsilon$ , $p = h\Upsilon(h\neq0)$ and $r = 0$ ,

$\begin{equation} \begin{split} v_{3, 9}(x, y, t) = -\frac{1}{2}\, k_{{1}}\Upsilon+{\frac {k_{{1}}\Upsilon\, {{\rm e}^{\Upsilon\, \zeta}}}{{ {\rm e}^{\Upsilon\, \zeta}}-h}}-{\frac {k_{{1}}h\Upsilon}{{{\rm e}^{\Upsilon\, \zeta}} -h}} . \end{split} \end{equation}$

(3.40)

Set. 3.4. With $q = \Upsilon$ , $r = h\Upsilon(h\neq0)$ and $p = 0$ ,

$\begin{equation} \begin{split} v_{3, 10}(x, y, t) = -\frac{1}{2}\, k_{{1}}\Upsilon-{\frac {k_{{1}}\Upsilon}{-1+h{{\rm e}^{\Upsilon\, \zeta}}}} . \end{split} \end{equation}$

(3.41)

Set. 3.5. With $p = 0$ , $r\neq0$ and $q\neq0$ ,

$\begin{equation} \begin{split} v_{3, 11}(x, y, t) = -\frac{1}{2}\, k_{{1}}q+{\frac {k_{{1}}q \left( \sinh \left( q\zeta \right) - \cosh \left( q\zeta \right) \right) }{-\cosh \left( q\zeta \right) + \sinh \left( q\zeta \right) -\lambda}} , \end{split} \end{equation}$

(3.42)

and

$\begin{equation} \begin{split} v_{3, 12}(x, y, t) = -\frac{1}{2}\, k_{{1}}q+{\frac {k_{{1}}q\mu}{\cosh \left( q\zeta \right) +\sinh \left( q\zeta \right) +\mu}} . \end{split} \end{equation}$

(3.43)

In above solutions $\zeta = x+y- \omega t$ .

4. Graphical discussion

We present depictions for the numerous wave forms found in the framework pursuant to assessment in this part of the paper. We compiled and graphically displayed waves such as dark solitary, bright, dark-bright, lump-like, dark, anti, and cuspon kinks in 2D, 3D, and contour modes through RMESEM. The results obtained are essential for interpreting the manner in which attributed physical phenomena operate. The objectives of the produced soliton solutions are to significantly expand our comprehension with regard to the theory of long waves with finite amplitude and the related field. Additionally, it has been graphically demonstrated that the solitons in the context of GZK-BBME take the shapes of kink solitons prominently.

Peculiar wave solutions in NPDEs with a supple, resilient, confined transition across the two asymptotic shifts are known as kink solitons. Such waves are seen in some NPDEs, including the GZK-BBME, which models wave propagation in a range of physical systems, including fluids and plasma.

Kink solitons are put into many groups according to the characteristics they exhibit, such as dark, lump-like, cuspon kink, dark-kink, grey kink, and dark-bright kink solitons. In contrast, bright kinks are concentrated, steady wave packets with an energy peak or accumulation that maintains its peak shape and dimension through transmission. Dark-bright kink solitons bring together the properties of both dark wave and bright kinks. Lump-like kinks display locally lump-shaped forms; during propagation, they can alter structure or orientation. Erratic, cusp-type field discontinuities separate cuspon kinks from smoother solitons. Lastly, A grey kink, which produces a waveform with a steeper amplification plunge compared to a black kink, is referred to as a confined, uniform transition with a non-zero deviance distinguishing two hyperbole phases. Being that they preserve their initial configuration as they propagate through the GZK-BBME, kink solitons are important for studying long-wave dispersal amplitude that is finite. Kink solitons persistence occurrence in a variety of media, such as water, provides insight into wave conduct, surf-to-wave interaction, and patterning integrity throughout time and space. Because kink solitons are stable, confined waves that transmit without altering form, they are useful for explaining fluid dynamics, nonlinear wave theory, and plasma physics scenarios. In particular, kink solitons may be used to describe ion-acoustic waveforms in polarized plasmas in plasma physics. Shallow waves of water influenced by both dispersive and nonlinear factors are described by them in fluid dynamics. Furthermore, they are perfect for researching the transmission of energy and communication in nonlinear optical systems and other dispersive media controlled by comparable dynamical equations due to their stability and durability. As demonstrated in this study, these applications highlight the significance of investigating their many forms and behaviors. This information is essential for describing & predicting wave motion in dynamics of fluids, plasma physics, and related purposes.

Remark 1: is plotted for $v_{ 1, 6}$ given in (3.13), which displays an anti-kink soliton profile.

Figure 1. The real part of anti-kink soliton solution

$v_{1, 6}$ , as described in (3.13), is represented in three dimension, with contours and in two dimension (

$y = -1$ ) for the following values of

$p : = 1; q : = 10; r : = 8; \omega : = -20; s_0 : = 10; t : = 2$ .

DownLoad: Full-Size Img PowerPoint

Remark 2: is plotted for $v_{1, 9 }$ given in (3.16), which displays an anti-kink soliton profile.

Figure 2. The anti-kink soliton solution

$v_{1, 9}$ , as described in (3.16), is represented in three dimension, with contours and in two dimension (

$y = -100$ ) for the following values of

$p : = 25; q : = 5; h : = 5; \Upsilon : = 5; r : = 0; \omega : = -10; s_0 : = 30; t : = 4$ .

DownLoad: Full-Size Img PowerPoint

Remark 3: is plotted for $v_{ 2, 5}$ given in (3.21), which displays an anti-kink soliton profile.

Figure 3. The anti-kink soliton solution

$v_{2, 5}$ , as described in (3.21), is represented in three dimension, with contours and in two dimension (

$y = 0$ ) for the following values of

$p : = 1; q : = 5; r : = 4; \omega : = 1; k_0 : = 1; k_1 : = 2; s_0 : = 5; t : = 0$ .

DownLoad: Full-Size Img PowerPoint

Remark 4: is plotted for $v_{ 2, 8}$ given in (3.24), the profile shows a bright kink soliton.

Figure 4. The bright kink soliton solution (also known as a hump kink)

$v_{2, 8}$ , as described in (3.24), is represented in three dimension, with contours and in two dimension (

$y = -1$ ) for the following values of

$p : = 4; q : = 10; r : = 4; \omega : = 1; k_0 : = 3; k_1 : = 6; s_0 : = 15; t : = 1$ .

DownLoad: Full-Size Img PowerPoint

Remark 5: is plotted for $v_{ 2, 9}$ given in (3.25), the profile shows a lump-like kink soliton.

Figure 5. The lump-type kink soliton solution

$v_{2, 9}$ , as described in (3.25), is represented in three dimension, with contours and in two dimension (

$y = -50$ ) for the following values of

$p : = 1; q : = 2; r : = 1; \omega : = 1; k_0 : = 4; k_1 : = 8; s_0 : = 20; t : = 10$ .

DownLoad: Full-Size Img PowerPoint

Remark 6: is plotted for $v_{ 2, 10}$ given in (3.26), the profile shows a lump-like kink soliton.

Figure 6. The lump-type kink soliton solution

$v_{2, 10}$ , as described in (3.26), is represented in three dimension, with contours and in two dimension (

$y = 1$ ) for the following values of

$p : = 15; q : = 0; r : = 0; \omega : = 1; k_0 : = 5; k_1 : = 10; s_0 : = 25; t : = 20$ .

DownLoad: Full-Size Img PowerPoint

Remark 7: is plotted for $v_{2, 12 }$ given in (3.28), the profile shows a cuspon anti-kink soliton.

Figure 7. The cuspon anti-kink soliton solution

$v_{2, 12}$ , as described in (3.28), is represented in three dimension, with contours and in two dimension (

$y = 50$ ) for the following values of

$p : = 6; q : = 3; h : = 2; \Upsilon : = 3; r : = 0; \omega : = 1; k_0 : = 10; k_1 : = 20; s_0 : = 50; t : = 100$ .

DownLoad: Full-Size Img PowerPoint

Remark 8: is plotted for $v_{ 3, 1}$ given in (3.32), the profile shows a bright-dark kink soliton.

Figure 8. The bright-dark kink soliton solution

$v_{3, 1}$ , as described in (3.32), is represented in three dimension, with contours and in two dimension (

$y = 0.3$ ) for the following values of

$p : = 6; q : = 1; r : = 5; \omega : = 10; k_0 : = 1; k_1 : = 5; s_0 : = 50; t : = 50$ .

DownLoad: Full-Size Img PowerPoint

Remark 9: is plotted for $v_{ 3, 5}$ given in (3.36), the profile shows a grey kink soliton.

Figure 9. The gray kink soliton solution

$v_{3, 5}$ , as described in (3.36), is represented in three dimension, with contours and in two dimension (

$y = 1000$ ) for the following values of

$p : = 4; q : = 10; r : = 4; \omega : = 20; k_1 : = 30; t : = 100$ .

DownLoad: Full-Size Img PowerPoint

Remark 10: is plotted for $v_{ 3, 10 }$ given in (3.41), the profile shows a solitary kink soliton.

Figure 10. The solitary kink soliton solution

$v_{3, 10}$ , as described in (3.41), is represented in three dimension, with contours and in two dimension (

$y = 45$ ) for the following values of

$p : = 0; q : = 10; r : = 20; \Upsilon : = 10; h : = 2; \omega : = 5; k_1 : = 10; t : = 50$ .

DownLoad: Full-Size Img PowerPoint

Remark 11: is plotted for $v_{ 3, 12}$ given in (3.43), the profile shows a solitary kink soliton.

Figure 11. The solitary kink soliton solution

$v_{3, 12}$ , as described in (3.43), is represented in three dimension, with contours and in two dimension (

$y = 50$ ) for the following values of

$p : = 0; q : = 5; r : = 5; \omega : = 1; k_0 : = 2; k_1 : = 3; s_0 : = 5; t : = 20; \mu : = 5$ .

DownLoad: Full-Size Img PowerPoint

5. Conclusions

The modernized RMESEM was established in this research to address a nonlinear model, namely GZK-BBME. With the help of the Riccati equation, the RMESEM was capable of arriving at a close form solution for the NODE that the model generated. The propagating soliton solutions that are significant to the problem's physical interpretation were subsequently obtained by shaping this solution into a system of nonlinear algebraic equations. It was shown that different travelling solitons, including dark-kink, lump-type, dark-bright, grey kink, and cuspon kink solitons, exist in kink soliton solutions by presenting multiple 3D, 2D, and contour graphs. The research highlights the implications for several practical applications in the linked fields of nonlinear GZK-BBME and demonstrates how the RMESEM may be utilized to build arrays of soliton solutions for difficult problems, particularly plasma physics and fluid dynamics. Thus, despite the fact that the analysis within the framework of the GZK-BBME provides insight into soliton dynamics that relate to the models of interest, it is constructive to also point out the drawbacks of using this method, especially when the largest derivative and nonlinear term are not equivalent. However, this limitation does not detract from the present study since this work clearly shows that the strategy used in this work is highly efficient, portable, and reliable for nonlinear problems of various natural science disciplines.

Appendix

Some of the above mentioned analytical methods rely on the Riccati equation. These methods are convenient to analyze soliton effects in nonlinear models since the equation of the Riccati type possesses solitary solutions ^[47]. Based on these applications of the Riccati hypothesis, the current study employed the Riccati equation comprising RMESEM ^[48] to generate and simulate soliton dynamics in GZK-BBME. This addition was useful as it generated five new families of kink soliton solutions for the targeted model: rational, hyperbolic, periodic, exponential and rational-hyperbolic. From the solutions obtained, significant progress was made towards the understanding of soliton behaviour and establishing a linkage between the events in the targeted model and the mentioned theories. Limiting our method's solutions results in some other strategies' solutions. The analogy is given in the part that follows:

Comparison with other analytical techniques

The outcomes of the other analytical techniques can be obtained using our procedure. As an instance:

Axiom 6.1.1. The following develops after $k_1 = 0$ is configured in (3.4):

$\begin{equation} {V}(\zeta) = {d}_{{0}} \left(\frac{1}{\Phi(\zeta)} \right). \end{equation}$

(A.1)

This shows that the closed-type result is associated with EDAM. Thus, attaining $k_1 = 0$ , our solutions can likewise lead to the results generated by EDAM.

Axiom 6.1.2. Similarly, the following develops after $s_0 = 0$ is configured in (3.4):

$\begin{equation} {V}(\zeta) = \sum\limits_{{j} = 0}^{{1}}{C}_{j} \left(\frac{\Phi'(\zeta)}{\Phi(\zeta)}\right)^{j}, \end{equation}$

(A.2)

This is the closed form solution obtained is applying the Riccati equation in the $(G'/G)$ -expansion approach.

As a result, the results of our study might potentially provide a wider range of results generated by the EDAM and $(G'/G)$ -expansion techniques.

Author contributions

N.M.A.A.; Conceptualization, H.Z.; formal analysis, H.Z.; investigation, N.M.A.A.; validation, H.Z.; visualization, N.M.A.A. and H.Z. funding; H.Z.; Data curation, H.Z.; resources, H.Z.; validation, N.M.A.A.; software, H.Z.; resources, N.M.A.A.; project administration, H.Z. writing-review & editing. All authors contributed equally. All authors have read and agreed to the published version of the manuscript.

Use of Generative-AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Conflict of interest

The authors declare that they have no conflicts of interest.

References

[1]	V. Prajapati, R. Meher, A robust analytical approach to the generalized Burgers-Fisher equation with fractional derivatives including singular and non-singular kernels, J. Ocean Eng. Sci., 2022. https://doi.org/10.1016/j.joes.2022.06.035 doi: 10.1016/j.joes.2022.06.035
[2]	L. Verma, R. Meher, Z. Avazzadeh, O. Nikan, Solution for generalized fuzzy fractional Kortewege-de Varies equation using a robust fuzzy double parametric approach, J. Ocean Eng. Sci., 2022. https://doi.org/10.1016/j.joes.2022.04.026 doi: 10.1016/j.joes.2022.04.026
[3]	M. Alqhtani, K. Saad, W. Weera, W. Hamanah, Analysis of the fractional-order local Poisson equation in fractal porous media, Symmetry, 14 (2022), 1323. https://doi.org/10.3390/sym14071323 doi: 10.3390/sym14071323
[4]	M. Al-Sawalha, R. Agarwal, R. Shah, O. Ababneh, W. Weera, A reliable way to deal with fractional-order equations that describe the unsteady flow of a polytropic gas, Mathematics, 10 (2022), 2293. https://doi.org/10.3390/math10132293 doi: 10.3390/math10132293
[5]	S. Mukhtar, S. Noor, The numerical investigation of a fractional-order multi-dimensional model of Navier-Stokes equation via novel techniques, Symmetry, 14 (2022), 1102. https://doi.org/10.3390/sym14061102 doi: 10.3390/sym14061102
[6]	N. Shah, Y. Hamed, K. Abualnaja, J. Chung, A. Khan, A comparative analysis of fractional-order Kaup-Kupershmidt equation within different operators, Symmetry, 14 (2022), 986. https://doi.org/10.3390/sym14050986 doi: 10.3390/sym14050986
[7]	N. Shah, H. Alyousef, S. El-Tantawy, J. Chung, Analytical investigation of fractional-order Korteweg-De-Vries-Type equations under Atangana-Baleanu-Caputo operator: Modeling nonlinear waves in a plasma and fluid, Symmetry, 14 (2022), 739. https://doi.org/10.3390/sym14040739 doi: 10.3390/sym14040739
[8]	N. Aljahdaly, A. Akgul, R. Shah, I. Mahariq, J. Kafle, A comparative analysis of the fractional-order coupled Korteweg-De Vries equations with the Mittag-Leffler law, J. Math., 2022 (2022), 1–30. https://doi.org/10.1155/2022/8876149 doi: 10.1155/2022/8876149
[9]	H. Ismael, H. Bulut, H. Baskonus, W-shaped surfaces to the nematic liquid crystals with three nonlinearity laws, Soft Comput., 25 (2020), 4513–4524. https://doi.org/10.1007/s00500-020-05459-6 doi: 10.1007/s00500-020-05459-6
[10]	M. Yavuz, T. Sulaiman, A. Yusuf, T. Abdeljawad, The Schrodinger-KdV equation of fractional order with Mittag-Leffler nonsingular kernel, Alex. Eng. J., 60 (2021), 2715–2724. https://doi.org/10.1016/j.aej.2021.01.009 doi: 10.1016/j.aej.2021.01.009
[11]	K. Safare, V. Betageri, D. Prakasha, P. Veeresha, S. Kumar, A mathematical analysis of ongoing outbreak COVID-19 in India through nonsingular derivative, Numer. Meth. Part. D. E., 37 (2020), 1282–1298. https://doi.org/10.1002/num.22579 doi: 10.1002/num.22579
[12]	H. Ismael, H. Baskonus, H. Bulut, Abundant novel solutions of the conformable Lakshmanan-Porsezian-Daniel model, Discrete Cont. Dyn.-S., 14 (2021), 2311. https://doi.org/10.3934/dcdss.2020398. doi: 10.3934/dcdss.2020398
[13]	H. Yasmin, N. Iqbal, A comparative study of the fractional coupled Burgers and Hirota-Satsuma KdV equations via analytical techniques, Symmetry, 14 (2022), 1364. https://doi.org/10.3390/sym14071364 doi: 10.3390/sym14071364
[14]	R. Alyusof, S. Alyusof, N. Iqbal, M. Arefin, Novel evaluation of the fractional acoustic wave model with the exponential-decay kernel, Complexity, 2022 (2022), 1–14. https://doi.org/10.1155/2022/9712388 doi: 10.1155/2022/9712388
[15]	M. Yavuz, T. Sulaiman, A. Yusuf, T. Abdeljawad, The Schrodinger-KdV equation of fractional order with Mittag-Leffler nonsingular kernel, Alex. Eng. J., 60 (2021), 2715–2724. https://doi.org/10.1016/j.aej.2021.01.009 doi: 10.1016/j.aej.2021.01.009
[16]	H. Triki, A. Biswas, Dark solitons for a generalized nonlinear Schrodinger equation with parabolic law and dual-power law nonlinearities, Math. Meth. Appl. Sci., 34 (2011), https://doi.org/10.1002/mma.1414 doi: 10.1002/mma.1414
[17]	L. Zhang, J. Si, New soliton and periodic solutions of (1+2)-dimensional nonlinear Schrodinger equation with dual-power law nonlinearity, Commun. Nonlinear Sci. Numer. Simul., 15 (2010), 2747–2754. https://doi.org/10.1016/j.cnsns.2009.10.028 doi: 10.1016/j.cnsns.2009.10.028
[18]	Q. Xu, S. Songhe, Numerical analysis of two local conservative methods for two-dimensional nonlinear Schrodinger equation, SCIENTIA SINICA Math., 48 (2017), 345. https://doi.org/10.1360/scm-2016-0308 doi: 10.1360/scm-2016-0308
[19]	J. C. Bronski, L. D. Carr, B. Deconinck, J. N. Kutz, Bose-Einstein condensates in standing waves: The cubic nonlinear Schrödinger equation with a periodic potential, Phys. Rev. Lett., 86 (2001), 1402–1405. https://doi.org/10.1103/PhysRevLett.86.1402 doi: 10.1103/PhysRevLett.86.1402
[20]	A. Trombettoni, A. Smerzi, Discrete solitons and breathers with dilute Bose-Einstein condensates, Phys. Rev. Lett., 86 (2001), 2353–2356. https://doi.org/10.1103/physrevlett.86.2353 doi: 10.1103/physrevlett.86.2353
[21]	A. Biswas, Quasi-stationary non-Kerr law optical solitons, Opt. Fiber Technol., 9 (2003), 224–259. https://doi.org/10.1016/s1068-5200(03)00044-0 doi: 10.1016/s1068-5200(03)00044-0
[22]	M. Eslami, M. Mirzazadeh, Topological 1-soliton solution of nonlinear Schrödinger equation with dual-power law nonlinearity in nonlinear optical fibers, Eur. Phys. J. Plus, 128 (2013). https://doi.org/10.1140/epjp/i2013-13140-y doi: 10.1140/epjp/i2013-13140-y
[23]	A. Hyder, The influence of the differential conformable operators through modern exact solutions of the double Schrödinger-Boussinesq system, Phys. Scripta, 2021. https://doi.org/10.1088/1402-4896/ac169f doi: 10.1088/1402-4896/ac169f
[24]	G. Adomian, R. Rach, Inversion of nonlinear stochastic operators, J. Math. Anal. Appl., 91 (1983), 39–46. https://doi.org/10.1016/0022-247x(83)90090-2 doi: 10.1016/0022-247x(83)90090-2
[25]	S. Alinhac, M. Baouendi, A counterexample to strong uniqueness for partial differential equations of Schrodinger's type, Commun. Part. Diff. Eq., 19 (1994), 1727–1733. https://doi.org/10.1080/03605309408821069 doi: 10.1080/03605309408821069
[26]	J. H. He, Homotopy perturbation technique, Comput. Meth. Appl. Mech. Eng., 178 (1999), 257–262. https://doi.org/10.1016/s0045-7825(99)00018-3 doi: 10.1016/s0045-7825(99)00018-3
[27]	J. H. He, A coupling method of homotopy technique and perturbation technique for nonlinear problems, Int. J. Non-Linear Mech., 35 (2000), 743. https://doi.org/10.1016/s0020-7462(98)00085-7 doi: 10.1016/s0020-7462(98)00085-7
[28]	D. D. Ganji, M. Rafei, Solitary wave solutions for a generalized Hirota Satsuma coupled KdV equation by homotopy perturbation method, Phys. Lett. A, 356 (2006), 131–137. https://doi.org/10.1016/j.physleta.2006.03.039 doi: 10.1016/j.physleta.2006.03.039
[29]	A. M. Siddiqui, R. Mahmood, Q. K. Ghori, Homotopy perturbation method for thin film flow of a fourth grade fluid down a vertical cylinder, Phys. Lett. A, 352 (2006), 404–410. https://doi.org/10.1016/j.physleta.2005.12.033 doi: 10.1016/j.physleta.2005.12.033
[30]	Y. Qin, A. Khan, I. Ali, M. Al Qurashi, H. Khan, R. Shah, D. Baleanu, An efficient analytical approach for the solution of certain fractional-order dynamical systems, Energies, 13 (2020), 2725. https://doi.org/10.3390/en13112725 doi: 10.3390/en13112725
[31]	M. Alaoui, R. Fayyaz, A. Khan, R. Shah, M. Abdo, Analytical investigation of Noyes-Field model for time-fractional Belousov-Zhabotinsky reaction, Complexity, 2021 (2021), 1–21. https://doi.org/10.1155/2021/3248376 doi: 10.1155/2021/3248376
[32]	J. H. He, Y. O. El-Dib, Homotopy perturbation method for Fangzhu oscillator, J. Math. Chem., 58 (2020), 2245–2253. https://doi.org/10.1007/s10910-020-01167-6 doi: 10.1007/s10910-020-01167-6
[33]	J. Honggang, Z. Yanmin, Conformable double Laplace-Sumudu transform decomposition method for fractional partial differential equations, Complexity, 2022 (2022), 1–8. https://doi.org/10.1155/2022/7602254 doi: 10.1155/2022/7602254
[34]	L. Akinyemi, O. S. Iyiola, Exact and approximate solutions of time-fractional models arising from physics via Shehu transform, Math. Meth. Appl. Sci., 43 (2020), 7442–7464. https://doi.org/10.1002/mma.6484 doi: 10.1002/mma.6484
[35]	L. Ali, R. Shah, W. Weera, Fractional view analysis of Cahn-Allen equations by new iterative transform method, Fractal Fract., 6 (2022), 293. https://doi.org/10.3390/fractalfract6060293 doi: 10.3390/fractalfract6060293
[36]	S. Maitama, W. Zhao, Homotopy perturbation Shehu transform method for solving fractional models arising in applied sciences, J. Appl. Math. Comput. Mech., 20 (2021), 71–82. https://doi.org/10.17512/jamcm.2021.1.07 doi: 10.17512/jamcm.2021.1.07
[37]	M. Liaqat, A. Khan, M. Alam, M. Pandit, S. Etemad, S. Rezapour, Approximate and Closed-Form solutions of Newell-Whitehead-Segel eEquations via modified conformable Shehu transform decomposition method, Math. Probl. Eng., 2022 (2022), 1–14. https://doi.org/10.1155/2022/6752455 doi: 10.1155/2022/6752455

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)