A new solitary wave solution of the fractional phenomena Bogoyavlenskii equation via Bäcklund transformation

Yousef Jawarneh; Humaira Yasmin; Ali M. Mahnashi; Yousef Jawarneh; Humaira Yasmin; Ali M. Mahnashi

doi:10.3934/math.20241678

AIMS Mathematics

2024, Volume 9, Issue 12: 35308-35325. doi: 10.3934/math.20241678

Previous Article Next Article

Research article Special Issues

A new solitary wave solution of the fractional phenomena Bogoyavlenskii equation via Bäcklund transformation

1.
Department of Mathematics, College of Science, University of Ha'il, Ha'il 2440, Saudi Arabia
2.
Department of Basic Sciences, General Administration of Preparatory Year, King Faisal University, P.O. Box 400, Al Ahsa 31982, Saudi Arabia
3.
Department of Mathematics and Statistics, College of Science, King Faisal University, P.O. Box 400, Al Ahsa 31982, Saudi Arabia
4.
Department of Mathematics, Faculty of Science, Jazan University, P.O. Box 2097, Jazan 45142, Kingdom of Saudi Arabia

Received: 20 September 2024 Revised: 04 December 2024 Accepted: 06 December 2024 Published: 18 December 2024
MSC : 34G20, 35A20, 35A22, 35R11

In this paper, we use the Riccati–Bernoulli sub-ODE method in conjunction with the Bäcklund transformation to find out the exact solutions of the nonlinear time–space fractional Bogoyavlenskii equation. The obtained solutions encompass multiple kink solitary wave solutions that are quite unique and important in addition to solutions presented in hyperbolic, trigonometric, and rational function forms. This equation describes central factors influencing its behavior including fluid dynamics in shallow water waves and plasma, which demonstrates our conclusions have broad applications for such systems. We also study the effect of the fractional order parameter ( $\alpha$ ) on solutions and plot their behavior using MATLAB in two dimensions. This work also contributes to the knowledge of the physical structures of the fractional Bogoyavlenskyi equation apart from showcasing the potential of the Riccati–Bernoulli sub-ODE method when applied to nonlinear fractional differential equations.

Keywords:

Citation: Yousef Jawarneh, Humaira Yasmin, Ali M. Mahnashi. A new solitary wave solution of the fractional phenomena Bogoyavlenskii equation via Bäcklund transformation[J]. AIMS Mathematics, 2024, 9(12): 35308-35325. doi: 10.3934/math.20241678

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]	M. Almheidat, H. Yasmin, M. Al Huwayz, R. Shah, S. A. El-Tantawy, A novel investigation into time-fractional multi-dimensional Navier–Stokes equations within Aboodh transform, Open Phys., 22 (2024), 20240081. https://doi.org/10.1515/phys-2024-0081 doi: 10.1515/phys-2024-0081
[2]	M. M. Al-Sawalha, H. Yasmin, S. Muhammad, Y. Khan, R. Shah, Optimal power management of a stand-alone hybrid energy management system: hydro-photovoltaic-fuel cell, Ain Shams Eng. J., 15 (2024), 103089. https://doi.org/10.1016/j.asej.2024.103089 doi: 10.1016/j.asej.2024.103089
[3]	A. H. Ganie, H. Yasmin, A. A. Alderremy, A. S. Alshehry, S. Aly, Fractional view analytical analysis of generalized regularized long wave equation, Open Phys., 22 (2024), 20240025. https://doi.org/10.1515/phys-2024-0025 doi: 10.1515/phys-2024-0025
[4]	A. S. Alshehry, H. Yasmin, R. Shah, A. Ali, I. Khan, Fractional-order view analysis of Fisher's and foam drainage equations within Aboodh transform, Eng. Computation., 41 (2024), 489–515. https://doi.org/10.1108/EC-08-2023-0475 doi: 10.1108/EC-08-2023-0475
[5]	H. Yasmin, A. S. Alshehry, A. H. Ganie, A. M. Mahnashi, R. Shah, Perturbed Gerdjikov–Ivanov equation: soliton solutions via Backlund transformation, Optik, 298 (2024), 171576. https://doi.org/10.1016/j.ijleo.2023.171576 doi: 10.1016/j.ijleo.2023.171576
[6]	H. Durur, A. Yokus, ( $1/G'$ )-Açılım metodunu kullanarak Sawada–Kotera denkleminin hiperbolik yürüyen dalga cözümleri, (Turkish), Afyon Kocatepe Üniversitesi Fen Ve Mühendislik Bilimleri Dergisi, 19 (2019), 615–619. https://doi.org/10.35414/akufemubid.559048 doi: 10.35414/akufemubid.559048
[7]	Sirendaoreji, S. Jiong, Auxiliary equation method for solving nonlinear partial differential equations, Phys. Lett. A, 309 (2003), 387–396. https://doi.org/10.1016/S0375-9601(03)00196-8 doi: 10.1016/S0375-9601(03)00196-8
[8]	S. Abbasbandy, F. S. Zakaria, Soliton solutions for the fifth-order KdV equation with the homotopy analysis method, Nonlinear Dyn., 51 (2008), 83–87. https://doi.org/10.1007/s11071-006-9193-y doi: 10.1007/s11071-006-9193-y
[9]	A. Korkmaz, O. E. Hepson, K. Hosseini, H. Rezazadeh, M. Eslami, Sine-Gordon expansion method for exact solutions to conformable time fractional equations in RLW-class, J. King Saud Univ. Sci., 32 (2020), 567–574. https://doi.org/10.1016/j.jksus.2018.08.013 doi: 10.1016/j.jksus.2018.08.013
[10]	Z.-Y. Zhang, Exact traveling wave solutions of the perturbed Klein–Gordon equation with quadratic nonlinearity in (1+1)-dimension, Part Ⅰ: without local inductance and dissipation effect, Turk. J. Phys., 37 (2013), 259–267. https://doi.org/10.3906/fiz-1205-13 doi: 10.3906/fiz-1205-13
[11]	Z.-Y. Zhang, Z.-H. Liu, X.-J. Miao, Y.-Z. Chen, New exact solutions to the perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity, Appl. Math. Comput., 216 (2010), 3064–3072. https://doi.org/10.1016/j.amc.2010.04.026 doi: 10.1016/j.amc.2010.04.026
[12]	Z.-Y. Zhang, Z.-H. Liu, X.-J. Miao, Y.-Z. Chen, Qualitative analysis and traveling wave solutions for the perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity, Phys. Lett. A, 375 (2011), 1275–1280. https://doi.org/10.1016/j.physleta.2010.11.070 doi: 10.1016/j.physleta.2010.11.070
[13]	M. N. Khan, Siraj-ul-Islam, I. Hussain, I. Ahmad, H. Ahmad, A local meshless method for the numerical solution of space-dependent inverse heat problems, Math. Method. Appl. Sci., 44 (2021), 3066–3079. https://doi.org/10.1002/mma.6439 doi: 10.1002/mma.6439
[14]	M. N. Khan, I. Ahmad, H. Ahmad, A radial basis function collocation method for space-dependent inverse heat problems, J. Appl. Comput. Mech., 6 (2020), 1187–1199. https://doi.org/10.22055/JACM.2020.32999.2123 doi: 10.22055/JACM.2020.32999.2123
[15]	K. R. Kamal, G. Rahmat, K. Shah, On the numerical approximation of three-dimensional time fractional convection-diffusion equations, Math. Probl. Eng., 2021 (2021), 4640467. https://doi.org/10.1155/2021/4640467 doi: 10.1155/2021/4640467
[16]	Y. Kai, Z. Yin, On the Gaussian traveling wave solution to a special kind of Schrödinger equation with logarithmic nonlinearity, Mod. Phys. Lett. B, 36 (2022), 2150543. https://doi.org/10.1142/S0217984921505436 doi: 10.1142/S0217984921505436
[17]	Y. Kai, Z. Yin, Linear structure and soliton molecules of Sharma-Tasso-Olver-Burgers equation, Phys. Lett. A, 452 (2022), 128430. https://doi.org/10.1016/j.physleta.2022.128430 doi: 10.1016/j.physleta.2022.128430
[18]	C. Zhu, S. A. Idris, M. E. M. Abdalla, S. Rezapour, S. Shateyi, B. Gunay, Analytical study of nonlinear models using a modified Schrödinger's equation and logarithmic transformation, Results Phys., 55 (2023), 107183. https://doi.org/10.1016/j.rinp.2023.107183 doi: 10.1016/j.rinp.2023.107183
[19]	T. A. A. Ali, Z. Xiao, H. Jiang, B. Li, A class of digital integrators based on trigonometric quadrature rules, IEEE Trans. Ind. Electron., 71 (2024), 6128–6138. https://doi.org/10.1109/TIE.2023.3290247 doi: 10.1109/TIE.2023.3290247
[20]	A. Zulfiqar, J. Ahmad, Soliton solutions of fractional modified unstable Schrödinger equation using Exp-function method, Results Phys., 19 (2020), 103476. https://doi.org/10.1016/j.rinp.2020.103476 doi: 10.1016/j.rinp.2020.103476
[21]	M. A. Akbar, N. H. M. Ali, M. T. Islam, Multiple closed form solutions to some fractional order nonlinear evolution equations in physics and plasma physics, AIMS Math., 4 (2019), 397–411. https://doi.org/10.3934/math.2019.3.397 doi: 10.3934/math.2019.3.397
[22]	T. Liu, Exact solutions to time-fractional fifth order KdV equation by trial equation method based on symmetry, Symmetry, 11 (2019), 742. https://doi.org/10.3390/sym11060742 doi: 10.3390/sym11060742
[23]	L. Bai, J. Qi, Y. Sun, Physical phenomena analysis of solution structures in a nonlinear electric transmission network with dissipative elements, Eur. Phys. J. Plus, 139 (2024), 9. https://doi.org/10.1140/epjp/s13360-023-04736-1 doi: 10.1140/epjp/s13360-023-04736-1
[24]	L. Bai, J. Qi, Y. Sun, Further physical study about solution structures for nonlinear q-deformed Sinh–Gordon equation along with bifurcation and chaotic behaviors, Nonlinear Dyn., 111 (2023), 20165–20199. https://doi.org/10.1007/s11071-023-08882-0 doi: 10.1007/s11071-023-08882-0
[25]	W. Chen, Time-space fabric underlying anomalous diffusion, Chaos Soliton. Fract., 28 (2006), 923–929. https://doi.org/10.1016/j.chaos.2005.08.199 doi: 10.1016/j.chaos.2005.08.199
[26]	W. Chen, F. Wang, B. Zheng, W. Cai, Non-Euclidean distance fundamental solution of Hausdorff derivative partial differential equations, Eng. Anal. Bound. Elem., 84 (2017), 213–219. https://doi.org/10.1016/j.enganabound.2017.09.003 doi: 10.1016/j.enganabound.2017.09.003
[27]	J.-H. He, A tutorial review on fractal spacetime and fractional calculus, Int. J. Theor. Phys., 53 (2014), 3698–3718. https://doi.org/10.1007/s10773-014-2123-8 doi: 10.1007/s10773-014-2123-8
[28]	J.-H. He, Seeing with a single scale is always unbelieving: from magic to two-scale fractal, Therm. Sci., 25 (2021), 1217–1219. https://doi.org/10.2298/TSCI2102217H doi: 10.2298/TSCI2102217H
[29]	M. Z. Sarikaya, H. Budak, F. Usta, On generalized the conformable fractional calculus, TWMS J. Appl. Eng. Math., 9 (2019), 792–799.
[30]	O. I. Bogoyavlenskii, Breaking solitons in 2+1-dimensional integrable equations, Russ. Math. Surv., 45 (1990), 1–86. https://doi.org/10.1070/RM1990v045n04ABEH002377 doi: 10.1070/RM1990v045n04ABEH002377
[31]	Y.-Z. Peng, M. Shen, On exact solutions of Bogoyavlenskii equation, Pramana J. Phys., 67 (2006), 449–456. https://doi.org/10.1007/s12043-006-0005-1 doi: 10.1007/s12043-006-0005-1
[32]	A. Malik, F. Chand, H. Kumar, S. C. Mishra, Exact solutions of the Bogoyavlenskii equation using the multiple $G'/G$ -expansion method, Comput. Math. Appl., 64 (2012), 2850–2859. https://doi.org/10.1016/j.camwa.2012.04.018 doi: 10.1016/j.camwa.2012.04.018
[33]	E. H. M. Zahran, M. M. A. Khater, Modified extended tanh-function method and its applications to the Bogoyavlenskii equation, Appl. Math. Model., 40 (2016), 1769–1775. https://doi.org/10.1016/j.apm.2015.08.018 doi: 10.1016/j.apm.2015.08.018
[34]	M. N. Alam, C. Tunc, An analytical method for solving exact solutions of the nonlinear Bogoyavlenskii equation and the nonlinear diffusive predator-prey system, Alex. Eng. J., 55 (2016), 1855–1865. https://doi.org/10.1016/j.aej.2016.04.024 doi: 10.1016/j.aej.2016.04.024
[35]	A. S. Alshehry, H. Yasmin, M. A. Shah, R. Shah, Analyzing fuzzy fractional Degasperis–Procesi and Camassa–Holm equations with the Atangana–Baleanu operator, Open Phys., 22 (2024), 20230191. https://doi.org/10.1515/phys-2023-0191 doi: 10.1515/phys-2023-0191
[36]	A. S. Alshehry, H. Yasmin, A. A. Khammash, R. Shah, Numerical analysis of dengue transmission model using Caputo–Fabrizio fractional derivative, Open Phys., 22 (2024), 20230169. https://doi.org/10.1515/phys-2023-0169 doi: 10.1515/phys-2023-0169
[37]	A. S. Alshehry, A. M. Mahnashi, Analyzing fractional PDE system with the Caputo operator and Mohand transform techniques, AIMS Math., 9 (2024), 32157–32181. https://doi.org/10.3934/math.20241544 doi: 10.3934/math.20241544
[38]	H. Yasmin, A. H. Almuqrin, Efficient solutions for time fractional Sawada-Kotera, Ito, and Kaup-Kupershmidt equations using an analytical technique, AIMS Math., 9 (2024), 20441–20466. https://doi.org/10.3934/math.2024994 doi: 10.3934/math.2024994
[39]	H. Yasmin, A. H. Almuqrin, Analytical study of time-fractional heat, diffusion, and Burger's equations using Aboodh residual power series and transform iterative methodologies, AIMS Math., 9 (2024), 16721–16752. https://doi.org/10.3934/math.2024811 doi: 10.3934/math.2024811
[40]	M. A. E. Abdelrahman, M. A. Sohaly, Solitary waves for the modified Korteweg-de Vries equation in deterministic case and random case, J. Phys. Math., 8 (2017), 214. https://doi.org/10.4172/2090-0902.1000214 doi: 10.4172/2090-0902.1000214
[41]	M. A. E. Abdelrahman, M. A. Sohaly, Solitary waves for the nonlinear Schrödinger problem with the probability distribution function in the stochastic input case, Eur. Phys. J. Plus, 132 (2017), 339. https://doi.org/10.1140/epjp/i2017-11607-5 doi: 10.1140/epjp/i2017-11607-5
[42]	X.-F. Yang, Z.-C. Deng, Y. Wei, A Riccati-Bernoulli sub-ODE method for nonlinear partial differential equations and its application, Adv. Differ. Equ., 2015 (2015), 117.
[43]	D. Lu, Q. Shi, New Jacobi elliptic functions solutions for the combined KdV-mKdV equation, Internatinal Journal of Nonlinear Science, 10 (2010), 320–325.
[44]	Y. Zhang, Solving STO and KD equations with modified Riemann–Liouville derivative using improved ( $G/G'$ )-expansion function method, International Journal of Applied Mathematics, 45 (2015), 16–22.
[45]	J. Lu, S. Shen, L. Chen, Variational approach for time-space fractal Bogoyavlenskii equation, Alex. Eng. J., 97 (2024), 294–301. https://doi.org/10.1016/j.aej.2024.04.031 doi: 10.1016/j.aej.2024.04.031

Reader Comments

Your name:*

Email:*
© 2024 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)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.1

Metrics

Article views(644) PDF downloads(34) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Mathematics

A new solitary wave solution of the fractional phenomena Bogoyavlenskii equation via Bäcklund transformation

Related Papers:

Abstract

1. Introduction

2. The operational methodology of RMESEM

3. Execution of RMESEM

4. Graphical discussion

5. Conclusions

Appendix

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Abstract

1. Introduction

2. The operational methodology of RMESEM

3. Execution of RMESEM

4. Graphical discussion

5. Conclusions

Appendix

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References

AIMS Mathematics

A new solitary wave solution of the fractional phenomena Bogoyavlenskii equation via Bäcklund transformation

Related Papers:

Abstract

1. Introduction

2. The operational methodology of RMESEM

3. Execution of RMESEM

4. Graphical discussion

5. Conclusions

Appendix

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. The operational methodology of RMESEM

3. Execution of RMESEM

4. Graphical discussion

5. Conclusions

Appendix

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References