Probing the diversity of soliton phenomena within conformable Estevez-Mansfield-Clarkson equation in shallow water

Mohammad Alqudah; Safyan Mukhtar; Haifa A. Alyousef; Sherif M. E. Ismaeel; S. A. El-Tantawy; Fazal Ghani; Mohammad Alqudah; Safyan Mukhtar; Haifa A. Alyousef; Sherif M. E. Ismaeel; S. A. El-Tantawy; Fazal Ghani

doi:10.3934/math.20241030

AIMS Mathematics

2024, Volume 9, Issue 8: 21212-21238. doi: 10.3934/math.20241030

Previous Article Next Article

Research article Special Issues

Probing the diversity of soliton phenomena within conformable Estevez-Mansfield-Clarkson equation in shallow water

1.
Department of Basic Sciences, School of Electrical Engineering & Information Technology, German Jordanian University, Amman 11180, Jordan
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 Physics, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
5.
Department of Physics, College of Science and Humanities in Al-Kharj, Prince Sattam bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
6.
Department of Physics, Faculty of Science, Ain Shams University, Cairo, Egypt
7.
Department of Physics, Faculty of Science, Port Said University, Port Said 42521, Egypt
8.
Department of Physics, Faculty of Science, Al-Baha University, Al-Baha P.O. Box 1988, Saudi Arabia
9.
Department of mathematics, Abdul Wali Khan University Mardan, Pakistan

Received: 31 January 2024 Revised: 26 May 2024 Accepted: 03 June 2024 Published: 01 July 2024
MSC : 34G20, 35A20, 35A22, 35R11

This study aims to employ the extended direct algebraic method (EDAM) to generate and evaluate soliton solutions to the nonlinear, space-time conformable Estevez Mansfield-Clarkson equation (CEMCE), which is utilized to simulate shallow water waves. The proposed method entails transforming nonlinear fractional partial differential equations (NFPDEs) into nonlinear ordinary differential equations (NODEs) under the assumption of a finite series solution by utilizing Riccati ordinary differential equations. Various mathematical structures/solutions for the current model are derived in the form of rational, exponential, trigonometric, and hyperbolic functions. The wide range of obtained solutions allows for a thorough analysis of their actual wave characteristics. The 3D and 2D graphs are used to illustrate that these behaviors consistently manifest as periodic, dark, and bright kink solitons. Notably, the produced soliton solutions offer new and critical insights into the intricate behaviors of the CEMCE by illuminating the basic mechanics of the wave's interaction and propagation. By analyzing these solutions, academics can better understand the model's behavior in various settings. These solutions shed light on complicated issues such as configuration dispersion in liquid drops and wave behavior in shallow water.

Keywords:

soliton solutions,
conformable Estevez-Mansfield-Clarkson equation,
conformable derivative,
shallow water,
extended direct algebraic method

Citation: Mohammad Alqudah, Safyan Mukhtar, Haifa A. Alyousef, Sherif M. E. Ismaeel, S. A. El-Tantawy, Fazal Ghani. Probing the diversity of soliton phenomena within conformable Estevez-Mansfield-Clarkson equation in shallow water[J]. AIMS Mathematics, 2024, 9(8): 21212-21238. doi: 10.3934/math.20241030

Related Papers:

[1]	M. Mossa Al-Sawalha, Saima Noor, Saleh Alshammari, Abdul Hamid Ganie, Ahmad Shafee . Analytical insights into solitary wave solutions of the fractional Estevez-Mansfield-Clarkson equation. AIMS Mathematics, 2024, 9(6): 13589-13606. doi: 10.3934/math.2024663
[2]	Naher Mohammed A. Alsafri . Solitonic behaviors in the coupled Drinfeld-Sokolov-Wilson system with fractional dynamics. AIMS Mathematics, 2025, 10(3): 4747-4774. doi: 10.3934/math.2025218
[3]	Mohammad Alqudah, Safyan Mukhtar, Albandari W. Alrowaily, Sherif. M. E. Ismaeel, S. A. El-Tantawy, Fazal Ghani . Breather patterns and other soliton dynamics in (2+1)-dimensional conformable Broer-Kaup-Kupershmit system. AIMS Mathematics, 2024, 9(6): 13712-13749. doi: 10.3934/math.2024669
[4]	Mahmoud Soliman, Hamdy M. Ahmed, Niveen Badra, Taher A. Nofal, Islam Samir . Highly dispersive gap solitons for conformable fractional model in optical fibers with dispersive reflectivity solutions using the modified extended direct algebraic method. AIMS Mathematics, 2024, 9(9): 25205-25222. doi: 10.3934/math.20241229
[5]	Chun Huang, Zhao Li . New soliton solutions of the conformal time derivative generalized $q$ -deformed sinh-Gordon equation. AIMS Mathematics, 2024, 9(2): 4194-4204. doi: 10.3934/math.2024206
[6]	Yanlin Li, Mohd Danish Siddiqi, Meraj Ali Khan, Ibrahim Al-Dayel, Maged Zakaria Youssef . Solitonic effect on relativistic string cloud spacetime attached with strange quark matter. AIMS Mathematics, 2024, 9(6): 14487-14503. doi: 10.3934/math.2024704
[7]	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
[8]	M. Mossa Al-Sawalha, Safyan Mukhtar, Azzh Saad Alshehry, Mohammad Alqudah, Musaad S. Aldhabani . Kink soliton phenomena of fractional conformable Kairat equations. AIMS Mathematics, 2025, 10(2): 2808-2828. doi: 10.3934/math.2025131
[9]	Yanlin Li, Dipen Ganguly, Santu Dey, Arindam Bhattacharyya . Conformal $\eta$ -Ricci solitons within the framework of indefinite Kenmotsu manifolds. AIMS Mathematics, 2022, 7(4): 5408-5430. doi: 10.3934/math.2022300
[10]	Chun Huang, Zhao Li . Soliton solutions of conformable time-fractional perturbed Radhakrishnan-Kundu-Lakshmanan equation. AIMS Mathematics, 2022, 7(8): 14460-14473. doi: 10.3934/math.2022797

Abstract

1. Introduction

Nonlinear fractional partial differential equations (NFPDEs) play a crucial role in the fields of applied science, engineering, and mathematics, as well as in modeling various physical nonlinear phenomena. The fields of fluid dynamics, fluid mechanics, neurons, optical fibers, electric circuits, water waves, plasma waves, capillary-gravity waves, chemical physics, and plasma physics have garnered significant attention due to the various fractional models developed ^[1,2]. Understanding the behaviors of distinct elements within various scientific fields is of the utmost significance, underscoring the necessity of exploring approaches to address non-functional partial differential equations (PDEs) ^[3,4,5].

The examination of soliton solutions to NFPDEs continues to hold academic significance, owing to the enhanced level of intricacy and comprehensiveness they provide compared to conventional solutions. A soliton can be defined as a solitary wave packet that exhibits self-reinforcing properties and traverses a medium without undergoing any changes in its shape or velocity. Fermion stability and robustness make them useful in many branches of physics and engineering. They facilitate the effective transmission of information and preserve coherence across long distances in nonlinear systems. Mathematicians have employed robust and various methodologies to pursue novel soliton solution outcomes. The notable methods include the Khater Ⅱ method ^[6], the modified Khater method ^[7], the ( $G^{\prime}/G$ )-expansion method ^[8], the $(G^{\prime2})$ -expansion method ^[9], the Kudryashov method ^[10,11], the Poincaré-Lighthill-Kuo method ^[12], the exp-function method ^[13], the fractional sub-equation method ^[14], the Sardar sub-equation method ^[15], the extended direct algebraic method (EDAM) ^{[16,17,18,19]}, and others.

In 1997, Mansfield and Clarkson introduced the Estevez-Mansfield-Clarkson Equation (EMCE), which is a nonlinear evolution equation of the fourth order ^[20]. This equation has proven advantageous in examining wave dynamics in shallow water, as it was initially employed in their investigation of pattern dispersion in liquid droplets. The EMCE is expressed in its formal form as follows:

$\begin{equation} \psi_{yyyt}+\varpi\psi_{yt}\psi_{y}+\varpi\psi_{yy}\psi_{t}+\psi _{tt} = 0. \end{equation}$

(1.1)

Here, $\psi\equiv\psi(y, t)$ indicates the wave function, and $\varpi$ is a non-zero constant. This study aims to create soliton solutions for the conformable EMCE (CEMCE) using the EDAM. This EMCE extension model substitutes integer-order derivatives with conformable fractional derivatives (CFDs). Numerous potential uses are introduced through the integration of CFDs in this framework. For example, CFDs, enable a more profound representation of the fractional-order kinetics present in some systems of matter that can enhance the precision of sims and forecasts. CFDs can help us better grasp the memory processes and non-local impacts prevalent in low-depth water interactions, which will help us simulate patterning generation, the dissemination of waves, and wave couplings in the EMCE. Implications for this improved system include tsunami modeling, coastal technology, research in the field of complicated hydrodynamics, and weather forecasting. Scholars can more effectively address real-world issues in related domains and get fresh knowledge into the behavior of shallow water swells by employing CFDs. The model is represented by the following equation ^[21]:

$\begin{equation} D_{t}^{\delta}(D_{y}^{\beta}(D_{y}^{\beta}(D_{y}^{\beta}\psi)))+\varpi D_{y}^{\beta}\psi D_{y}^{\beta}(D_{t}^{\delta}\psi)+\varpi D_{t}^{\delta}\psi D_{y}^{\beta}(D_{y}^{\beta}\psi)+D_{t}^{\delta}(D_{t}^{\delta}\psi ) = 0, \end{equation}$

(1.2)

where $D_{y}^{\beta}(\cdot)$ and $D_{t}^{\delta}(\cdot)$ are CFDs and $0 < \beta, \delta\leq1$ . CFDs provide an effective operator to explain the physical behavior of phenomena have a non-integer behavior, including occurrences that show effects of memory alongside their non-local and local behavior. In conventional math, localized behavior wherein the degree of shift at a location solely relies on its current surroundings is captured by conventional derivatives. Nonetheless, many systems in the real world exhibit non-local behavior, in which the system's past or remote locations impact the current condition of the structure at a given time. Fractional derivatives, including CFDs, extend calculus to explain such a behavior. For example, fractional interactions are a property of viscoelastic substances in the field of materials, wherein the tension or strain at a place depends on the material's distortion behavior and its current environment ^[22,23]. CFDs, which consider the entire system's behavior throughout an uninterrupted variety of past events, precisely represent this memory impact. Furthermore, the propagation of molecules in the diffusion process frequently exhibits non-local patterns. This results in fractional diffusing problems, where the spread of molecules across a non-local region determines the degree to which intensity changes occur at a given place. To accurately simulate these kinds of occurrences, CFDs offer a mathematical tool that bridges the gap between local and non-local behaviors and considers memory-related factors that are part of the system's behavior ^[24,25,26].

Prior to the present study, numerous scholars have dealt with the EMCE and the fractional EMCE (FEMCE). For example, the Kudryashov approach ^[21] yielded solutions to the nonlinear fractional space-time FEMCE, which formed exponential functions. Kinks and periodic waves characterized the behaviors found in the solution. The space-time fractional Ablowitz-Kaup-Newell-Segur (AKNS) water equation and the exact traveling wave solutions FEMCE were investigated by Phoosree & Chinviriyasit using the ( $G^{\prime}/G$ )-expansion approach in terms of hyperbolic functions, trigonometric functions, and rational functions ^[27]. Ultimately, Phoosree and Thadee ^[28] were able to obtain traveling wave solutions for the FEMCE by employing Jumarie's Riemann-Liouville derivative, thereby establishing the solution in the finite series, along with the simple equation technique and the Bernoulli equation. However in the present study, we will construct soliton solutions for the FEMCE using a novel EDAM approach. The recommended method transforms NFPDEs into nonlinear ordinary differential equations (NODEs) using Riccati Ordinary Differential Equations (ODE) and a finite series solution. This investigation has led to the discovery of soliton solutions that exhibit various mathematical structures, including trigonometric, exponential, rational, and hyperbolic functions. These structures allow for a detailed examination of their actual wave behaviors. In particular, we show that these behaviors are periodic, shock, and kink solitons. This is supported by showing how they look in contour, 3D, and 2D graphs. Interestingly, our results produce a wider variety of results than previous solutions, thus providing new and significant insights into the underlying phenomena of the model under study.

The rest of this investigation is organized in the following sections: Section 2 describes the materials and methodology for the EDAM; Section 3 shows the soliton solutions for the CEMCE; Section 4 includes a discussion and the corresponding graphics; and in Section 5, the most important findings are summarized in brief points.

2. Methodology & materials

2.1. The basic definition for the conformable fractional derivatives (CFDs)

By using the distinct advantages CFDs offer over other fractional derivative operators, it is possible to explicitly solve NFPDEs. Remarkably, different formulations of the fractional derivative for Eq (1.2) do not produce traveling wave and soliton solutions because they violate the chain rule ^[29,30]. Consequently, Eq (1.2) was modified to include the CFDs operator of order $\beta$ as follows ^[31]:

$\begin{equation} D_{\vartheta}^{\beta}w(\vartheta) = \lim\limits_{\gamma\rightarrow0}\frac {w(\gamma\vartheta^{1-\beta}+\vartheta)-w(\vartheta)}{\gamma}, \quad\beta \in(0, 1]. \end{equation}$

(2.1)

Here, we employ the following attributes of this derivative:

$\begin{equation} D_{\vartheta}^{\beta}\vartheta^{a} = a\vartheta^{a-\beta}, \end{equation}$

(2.2)

$\begin{equation} D_{\vartheta}^{\beta}(a_{1}\rho(\vartheta)\pm a_{2}\eta(\vartheta )) = a_{1}D_{\vartheta}^{\beta}(\rho(\vartheta))\pm a_{2}D_{\vartheta}^{\beta }(\eta(\vartheta)), \end{equation}$

(2.3)

$\begin{equation} D_{\vartheta}^{\beta}\chi\lbrack\zeta(\vartheta)] = \chi_{\zeta}^{\prime} (\zeta(\vartheta))D_{\vartheta}^{\beta}\zeta(\vartheta), \end{equation}$

(2.4)

where $\rho(\vartheta)$ , $\eta(\vartheta)$ , $\chi(\vartheta)$ , and $\zeta(\vartheta)$ are arbitrary differentiable functions, whereas $\left(a, a_{1}, a_{2}\right)$ are constants.

2.2. The working procedure of the EDAM

Here, we outline the operating strategy of the EDAM. Consider the following conditions and factors when working with the NFPDE:

$\begin{equation} R(\psi, D_{t}^{\alpha}\psi, D_{y_{1}}^{\beta}\psi, D_{y_{2}}^{\gamma}\psi, \psi D_{y_{1}}^{\beta}\psi, \ldots) = 0, \ \ 0 < \alpha, \beta, \gamma\leq1, \end{equation}$

(2.5)

where $\psi\equiv\psi(t, y_{1}, y_{2}, y_{3}, \ldots, y_{r})$ .

The following procedures are applied to solve Eq (2.5):

(1) We begin with a variable transformation of the type $\psi(t, y_{1}, y_{2}, y_{3}, \dots, y_{r})\equiv\Psi(\vartheta)$ , where $\vartheta$ can be written in several ways. This transformation converts the NFPDE (2.5) into the following NODE:

$\begin{equation} Q(\Psi, \Psi^{\prime}\Psi, \Psi^{\prime}, \dots) = 0, \end{equation}$

(2.6)

where $\Psi^{\prime} = \frac{d\Psi}{d\vartheta}$ .

Equation (2.6) may occasionally be integrated to make the NODE suitable for the homogenous balance principle.

(2) After that, by using the Riccati ODE, we consider the subsequent series-based solution to the NODE in Eq (2.6):

$\begin{equation} \Psi(\vartheta) = \sum\limits_{\sigma = -q}^{q}p_{\sigma}\zeta(\vartheta)^{\sigma}. \end{equation}$

(2.7)

In this context, $p_{\sigma}(\sigma = -q, ..., q)$ denotes the unknown constants, and $\zeta(\vartheta)$ is the general solution of the resulting ODE:

$\begin{equation} \zeta^{\prime}(\vartheta) = \lambda+\mu\zeta(\vartheta)+\nu\zeta(\vartheta)^{2}, \end{equation}$

(2.8)

where $\left(\lambda, \mu, \nu\right)$ are invariable.

(3) By creating a homogeneous balance for Eq (2.6) between the largest nonlinearity and the highest-order derivative, the positive integer $q$ introduced in Eq (2.7) is finally obtained.

(4) After that, we insert Eq (2.7) into either Eq (2.6) or the equation that arises from the integration of Eq (2.6) and combine the terms of $\zeta(\vartheta)$ into equivalent orders. This method results in an equation that includes $\zeta(\vartheta)$ . Then, the coefficients in this expression are set to zero, producing an algebraic system of equations for the variables $p_{\sigma}(\sigma = -q, ..., q)$ and other associated parameters.

(5) We apply MAPLE to solve this system of nonlinear algebraic equations.

(6) Then, the analytical soliton solutions for Eq (2.5) are then derived by computing and substituting the unknown values into Eq (2.7), together with the $\zeta(\vartheta)$ (the solution of Eq (2.8)). We can create the resulting families of soliton solutions using the general solution of Eq (2.8).

Family 1: For $M\left( = \mu^{2}-4\nu\lambda\right) < 0$ & $\nu\neq0$ , we obtain the following:

$\begin{align*} \zeta_{1}(\vartheta) & = -{\frac{\mu}{2\nu}}+{\frac{\sqrt{-M}}{2\nu}} \tan\left( \frac{1}{2}\, \sqrt{-M}\vartheta\right), \\ \zeta_{2}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{-M}}{2\nu}} \cot\left( \frac{1}{2}\, \sqrt{-M}\vartheta\right) , \\ \zeta_{3}(\vartheta) & = -{\frac{\mu}{2\nu}}+{\frac{\sqrt{-M}}{2\nu}}\left[ \tan\left( \sqrt{-M}\vartheta\right) +\sec\left( \sqrt{-M}\vartheta\right) \right] , \\ \zeta_{4}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{-M}}{2\nu}}\left[ \cot\left( \sqrt{-M}\vartheta\right) +\csc\left( \sqrt{-M}\vartheta\right) \right] , \\ \zeta_{5}(\vartheta) & = -{\frac{\mu}{2\nu}}+{\frac{\sqrt{-M}}{4\nu}}\left[ \tan\left( \frac{1}{4}\, \sqrt{-M}\vartheta\right) -\cot\left( \frac{1} {4}\, \sqrt{-M}\vartheta\right) \right] . \end{align*}$

Family 2: For $M > 0$ & $\nu\neq0$ , we obtain the following:

$\begin{align*} \zeta_{6}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{M}}{2\nu}} \tanh\left( \frac{1}{2}\, \sqrt{M}\vartheta\right) , \\ \zeta_{7}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{M}}{2\nu}} \coth\left( \frac{1}{2}\, \sqrt{M}\vartheta\right) , \\ \zeta_{8}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{M}}{2\nu}}\left[ \tanh\left( \sqrt{M}\vartheta\right) +\mathit{sech}\left( \sqrt{M} \vartheta\right) \right] , \\ \zeta_{9}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{M}}{2\nu}}\left[ \coth\left( \sqrt{M}\vartheta\right) +\mathit{csch}\left( \sqrt{M} \vartheta\right) \right] , \\ \zeta_{10}(\vartheta) & = -{\frac{\mu}{2\nu}}-{\frac{\sqrt{M}}{4\nu}}\left[ \tanh\left( \frac{1}{4}\, \sqrt{M}\vartheta\right) -\coth\left( \frac{1} {4}\, \sqrt{M}\vartheta\right) \right] . \end{align*}$

Family 3: For $\nu\lambda > 0$ & $\mu = 0$ , we obtain the following:

$\begin{align*} \zeta_{11}(\vartheta) & = \sqrt{{\frac{\lambda}{\nu}}}\tan\left( \sqrt{\lambda\nu}\vartheta\right) , \\ \zeta_{12}(\vartheta) & = -\sqrt{{\frac{\lambda}{\nu}}}\cot\left( \sqrt{\lambda\nu}\vartheta\right) , \\ \zeta_{13}(\vartheta) & = \sqrt{{\frac{\lambda}{\nu}}}\left[ \tan\left( 2\, \sqrt{\lambda\nu}\vartheta\right) +\sec\left( 2\, \sqrt{\lambda\nu }\vartheta\right) \right] , \\ \zeta_{14}(\vartheta) & = -\sqrt{{\frac{\lambda}{\nu}}}\left[ \cot\left( 2\, \sqrt{\lambda\nu}\vartheta\right) +\csc\left( 2\, \sqrt{\lambda\nu }\vartheta\right) \right] , \\ \zeta_{15}(\vartheta) & = \frac{1}{2}\sqrt{{\frac{\lambda}{\nu}}}\left[ \tan\left( \frac{1}{2}\, \sqrt{\lambda\nu}\vartheta\right) -\cot\left( \frac{1}{2}\, \sqrt{\lambda\nu}\vartheta\right) \right] . \end{align*}$

Family 4: For $\lambda\nu < 0$ and $\mu = 0$ , we obtain the following:

$\begin{align*} \zeta_{16}(\vartheta) & = -\sqrt{-{\frac{\lambda}{\nu}}}\tanh\left( \sqrt{-\lambda\nu}\vartheta\right) , \\ \zeta_{17}(\vartheta) & = -\sqrt{-{\frac{\lambda}{\nu}}}\coth\left( \sqrt{-\lambda\nu}\vartheta\right) , \\ \zeta_{18}(\vartheta) & = -\sqrt{-{\frac{\lambda}{\nu}}}\left[ \tanh\left( 2\, \sqrt{-\lambda\nu}\vartheta\right) +i\mathit{sech}\left( 2\, \sqrt {-\lambda\nu}\vartheta\right) \right] , \\ \zeta_{19}(\vartheta) & = -\sqrt{-{\frac{\lambda}{\nu}}}\left[ \coth\left( 2\, \sqrt{-\lambda\nu}\vartheta\right) +\mathit{csch}\left( 2\, \sqrt {-\lambda\nu}\vartheta\right) \right] , \\ \zeta_{20}(\vartheta) & = -\frac{1}{2}\sqrt{-{\frac{\lambda}{\nu}}}\left[ \tanh\left( \frac{1}{2}\, \sqrt{-\lambda\nu}\vartheta\right) +\coth\left( \frac{1}{2}\, \sqrt{-\lambda\nu}\vartheta\right) \right] . \end{align*}$

Family 5: For $\lambda = \nu$ & $\mu = 0$ , we obtain the following:

$\zeta_{21}(\vartheta) = \tan\left( \lambda\vartheta\right) ,$

$\zeta_{22}(\vartheta) = -\cot\left( \lambda\vartheta\right) ,$

$\zeta_{23}(\vartheta) = \tan\left( 2\, \lambda\vartheta\right) +\sec\left( 2\, \lambda\vartheta\right) ,$

$\zeta_{24}(\vartheta) = -\cot\left( 2\, \lambda\vartheta\right) +\csc\left( 2\, \lambda\vartheta\right) ,$

and

$\zeta_{25}(\vartheta) = \frac{1}{2}\tan\left( \frac{1}{2}\, \lambda \vartheta\right) -\frac{1}{2}\, \cot\left( \frac{1}{2}\, \lambda \vartheta\right) .$

Family 6: For $\nu = -\lambda$ & $\mu = 0$ , we obtain the following:

$\begin{align*} \zeta_{26}(\vartheta) & = -\tanh\left( \lambda\vartheta\right) , \\ \zeta_{27}(\vartheta) & = -\coth\left( \lambda\vartheta\right) , \\ \zeta_{28}(\vartheta) & = -\tanh\left( 2\, \lambda\vartheta\right) +i\mathit{sech}\left( 2\, \lambda\vartheta\right) , \\ \zeta_{29}(\vartheta) & = -\coth\left( 2\, \lambda\vartheta\right) +\mathit{csch}\left( 2\, \lambda\vartheta\right) , \\ \zeta_{30}(\vartheta) & = -\frac{1}{2}\tanh\left( \frac{1}{2}\, \lambda \vartheta\right) -\frac{1}{2}\, \coth\left( \frac{1}{2}\, \lambda \vartheta\right) . \end{align*}$

Family 7: For $M = 0$ , we obtain the following:

$\zeta_{31}(\vartheta) = -2\, {\frac{\lambda\left( \mu\vartheta\, +2\right) }{{\mu}^{2}\vartheta\, }}.$

Family 8: For $\mu = \tau$ , $\lambda = h\tau(h\neq0)$ & $\nu = 0$ , we obtain the following:

$\zeta_{32}(\vartheta) = {e}^{\tau\, \vartheta}-h.$

Family 9: For $\mu = \nu = 0$ , we obtain the following:

$\zeta_{33}(\vartheta) = \vartheta\, \lambda.$

Family 10: For $\mu = \lambda = 0$ , we obtain the following:

$\zeta_{34}(\vartheta) = -{\frac{1}{\vartheta\, \nu}}.$

Family 11: For $\lambda = 0$ , $\mu\neq0$ & $\nu\neq0$ , we obtain the following:

$\zeta_{35}(\vartheta) = -{\frac{\mu}{\nu\left( 1+\cosh\left( \mu \vartheta\right) -\sinh\left( \mu\vartheta\right) \right) }},$

and

$\zeta_{36}(\vartheta) = -{\frac{\mu\left[ \cosh\left( \mu\vartheta\right) +\sinh\left( \mu\vartheta\right) \right] }{\nu\left[ \cosh\left( \mu\vartheta\right) +\sinh\left( \mu\vartheta\right) +1\right] }}.$

Family 12: For $\mu = \tau$ , $\nu = h\tau(h\neq0)$ & $\lambda = 0$ , we obtain the following:

$\zeta_{37}(\vartheta) = {\frac{{e}^{\upsilon\, \vartheta}}{1-h{e}^{\upsilon \vartheta}}}.$

3. Execution of EDAM to CEMCE

In this study phase, the EDAM is employed to create soliton solutions for the CEMCE, as depicted in Eq (1.2). To initiate the technique, one must perform the following variable transformation:

$\begin{equation} \psi(x, y, t)\equiv\Psi(\vartheta)\;\text{with }\vartheta = \frac{\epsilon x^{\alpha}}{\alpha}+\frac{\kappa y^{\beta}}{\beta}-\frac{\omega t^{\delta} }{\delta}, \end{equation}$

(3.1)

which, with the help of Eqs (2.2) and (2.4), we have the following:

$\begin{split} & D_{t}^{\delta}\psi = D_{t}^{\delta}\Psi(\vartheta) = \Psi_{t}^{\prime\delta }(\vartheta) = -\omega\Psi^{\prime}(\vartheta), \\ & D_{t}^{\delta}(D_{t}^{\delta}\psi) = D_{t}^{\delta}(-\omega\Psi^{\prime }(\vartheta)) = -\omega\Psi_{t}^{\prime\prime\delta}\vartheta = \omega^{2} \Psi^{\prime\prime}(\vartheta), \\ & D_{y}^{\beta}\psi = D_{y}^{\beta}\Psi(\vartheta) = \Psi_{y}^{\prime\beta }(\vartheta) = \kappa\Psi^{\prime}(\vartheta), \\ & D_{y}^{\beta}(D_{y}^{\beta}\psi) = D_{y}^{\beta}(\kappa\Psi^{\prime }(\vartheta)) = \kappa\Psi_{y}^{\prime\prime\beta}\vartheta = \kappa^{2} \Psi^{\prime\prime}(\vartheta), \\ & D_{y}^{\beta}(D_{y}^{\beta}(D_{y}^{\beta}\psi)) = D_{y}^{\beta}(\kappa ^{2}\Psi^{\prime\prime2}\Psi_{y}^{\prime\prime\prime\beta}\vartheta = \kappa ^{3}\Psi^{\prime\prime\prime}(\vartheta), \\ & D_{y}^{\beta}(D_{y}^{\beta}(D_{y}^{\beta}(D_{y}^{\beta}\psi))) = D_{y} ^{\beta}(\kappa^{3}\Psi^{\prime\prime\prime3}\Psi^{(iv)}(\vartheta )D_{y}^{\beta}\vartheta = \kappa^{4}\Psi^{(iv)}(\vartheta), \\ & D_{t}^{\delta}(D_{y}^{\beta}\psi) = D_{t}^{\delta}(\kappa\Psi^{\prime }(\vartheta)) = \kappa\Psi_{t}^{\prime\prime\delta}\vartheta = -\omega\kappa \Psi^{\prime\prime}(\vartheta), \\ & D_{t}^{\delta}(D_{y}^{\beta}(D_{y}^{\beta}(D_{y}^{\beta}\psi))) = D_{t} ^{\delta}(\kappa^{3}\Psi^{\prime\prime\prime3}\Psi^{(iv)}(\vartheta )D_{t}^{\delta}\vartheta = -\omega\kappa^{3}\Psi^{(iv)}(\vartheta). \end{split}$

Now, by collecting all the above terms and adding them in Eq (1.2), the following NODE is obtained:

$\begin{equation} -\kappa^{3}\Psi^{(iv)}-2\kappa^{2}\varpi\Psi^{\prime}\Psi^{\prime\prime }+\omega\Psi^{\prime\prime} = 0. \end{equation}$

(3.2)

It is essential to observe that the variable transformation (3.1) is expressed in a two-dimensional format, although our problem is one-dimensional. Therefore, we can disregard the second dimension when performing our computations, (i.e., $\vartheta = \frac{\kappa y^{\beta}}{\beta }-\frac{\omega t^{\delta}}{\delta}$ ).

The result obtained from a single integration of Eq (3.2) with a constant of integration of zero is as follows:

$\begin{equation} -\kappa^{3}\Psi^{\prime\prime\prime}-\kappa^{2}\varpi\left( \Psi^{\prime }\right) ^{2}+\omega\Psi^{\prime} = 0. \end{equation}$

(3.3)

Upon constructing a homogeneous balancing condition between $\Psi ^{\prime\prime\prime}\longmapsto$ ( $3q$ ) and $(\left(\Psi^{\prime}\right) ^{2}$ $\longmapsto$ ( $q+2$ ), we conclude that $q$ $= 1$ . We obtain the following series-based solution for Eq (3.3) by substituting $q = 1$ into Eq (2.7):

$\begin{equation} \Psi(\vartheta) = \sum\limits_{\sigma = -1}^{1}p_{\sigma}\zeta(\vartheta)^{\sigma} = p_{-1}\zeta(\vartheta)^{-1}+p_{0}+p_{1}\zeta(\vartheta). \end{equation}$

(3.4)

We obtain an expression in $\zeta(\vartheta)$ by inserting Eq (3.4) into Eq (3.3) and accumulating the terms with similar powers of $\zeta (\vartheta)$ . The procedure produces a set of algebraic nonlinear equations when the coefficients are set to zero. When MAPLE is employed to solve this system, the two sets of solutions provided are as follows:

Case 1.

$\begin{equation} p_{{0}} = p_{{0}}, p_{{1}} = 0, p_{{-1}} = 3\, {\frac{\kappa\, \lambda}{\varpi}} , \omega = {\kappa}^{3}M, \kappa = \kappa. \end{equation}$

(3.5)

Case 2.

$\begin{equation} p_{{0}} = p_{{0}}, p_{{1}} = p_{{1}}, p_{{-1}} = 0, \omega = {\frac{-{\varpi}^{3}{p_{{1} }}^{3}M}{27{\nu}^{3}}}, \kappa = {\frac{-\varpi\, p_{{1}}}{3\nu}}. \end{equation}$

(3.6)

Now, we find all solutions for the two above cases.

Given case 1, we create the families of soliton solutions for the CEMCE (1.2) by employing Eqs (3.1) and (3.4) and the associated solution of Eq (2.8):

Family 1.1: For $M < 0\quad\nu\neq0$ , we obtain the following:

$\begin{equation} \psi_{1, 1} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}+\frac{1}{2}\, {\frac{\sqrt{-M}\tan\left( \frac{1}{2}\, \sqrt {-M}\vartheta\right) }{\nu}}\right] ^{-1}+p_{{0}}, \end{equation}$

(3.7)

$\begin{equation} \psi_{1, 2} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{2}\, {\frac{\sqrt{-M}\cot\left( \frac{1}{2}\, \sqrt {-M}\vartheta\right) }{\nu}}\right] ^{-1}+p_{{0}}, \end{equation}$

(3.8)

$\begin{equation} \psi_{1, 3} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}+\frac{1}{2}\, {\frac{\sqrt{-M}\left( \tan\left( \sqrt{-M} \vartheta\right) +\sec\left( \sqrt{-M}\vartheta\right) \right) }{\nu} }\right] ^{-1}+p_{{0}}, \end{equation}$

(3.9)

$\begin{equation} \psi_{1, 4} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{2}\, {\frac{\sqrt{-M}\left( \cot\left( \sqrt{-M} \vartheta\right) +\csc\left( \sqrt{-M}\vartheta\right) \right) }{\nu} }\right] ^{-1}+p_{{0}}, \end{equation}$

(3.10)

and

$\begin{equation} \psi_{1, 5} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}+\frac{1}{4}\, {\frac{\sqrt{-M}\left( \tan\left( \frac{1}{4} \, \sqrt{-M}\vartheta\right) -\cot\left( \frac{1}{4}\, \sqrt{-M} \vartheta\right) \right) }{\nu}}\right] ^{-1}+p_{{0}}. \end{equation}$

(3.11)

Family 1.2: For $M > 0\quad\nu\neq0$ , we obtain the following:

$\begin{equation} \psi_{1, 6} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{2}\, {\frac{\sqrt{M}\tanh\left( \frac{1}{2}\, \sqrt {M}\vartheta\right) }{\nu}}\right] ^{-1}+p_{{0}}, \end{equation}$

(3.12)

$\begin{equation} \psi_{1, 7} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{2}\, {\frac{\sqrt{M}\coth\left( \frac{1}{2}\, \sqrt {M}\vartheta\right) }{\nu}}\right] ^{-1}+p_{{0}}, \end{equation}$

(3.13)

$\begin{equation} \psi_{1, 8} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{2}\, {\frac{\sqrt{M}\left( \tanh\left( \sqrt{M} \vartheta\right) +i\mathit{sech}\left( \sqrt{M}\vartheta\right) \right) }{\nu}}\right] ^{-1}+p_{{0}}, \end{equation}$

(3.14)

$\begin{equation} \psi_{1, 9} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{2}\, {\frac{\sqrt{M}\left( \coth\left( \sqrt{M} \vartheta\right) +\mathit{csch}\left( \sqrt{M}\vartheta\right) \right) }{\nu}}\right] ^{-1}+p_{{0}}, \end{equation}$

(3.15)

and

$\begin{equation} \psi_{1, 10} = \frac{3\, \kappa\, \lambda}{\varpi}\left[ -\frac{1}{2}\, {\frac{\mu }{\nu}}-\frac{1}{4}\, {\frac{\sqrt{M}\left[ \tanh\left( \frac{1}{4}\, \sqrt {M}\vartheta\right) -\coth\left( \frac{1}{4}\, \sqrt{M}\vartheta\right) \right] }{\nu}}\right] ^{-1}+p_{{0}}. \end{equation}$

(3.16)

Family 1.3: For $\lambda\nu > 0$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{1, 11} = \frac{3\, \kappa\, }{\varpi}\sqrt{\lambda\nu}\left( \tan\left( \sqrt{\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.17)

$\begin{equation} \psi_{1, 12} = -\frac{3\, \kappa\, }{\varpi}\sqrt{\lambda\nu}\left( \cot\left( \sqrt{\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.18)

$\begin{equation} \psi_{1, 13} = \frac{3\, \kappa\, }{\varpi}\sqrt{\lambda\nu}\left( \tan\left( 2\, \sqrt{\lambda\, \nu}\vartheta\right) +\sec\left( 2\, \sqrt{\lambda\, \nu }\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.19)

$\begin{equation} \psi_{1, 14} = -\frac{3\, \kappa\, }{\varpi}\sqrt{\lambda\nu}\left( \cot\left( 2\, \sqrt{\lambda\, \nu}\vartheta\right) +\csc\left( 2\, \sqrt{\lambda\, \nu }\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.20)

and

$\begin{equation} \psi_{1, 15} = \frac{6\, \kappa\, }{\varpi}\sqrt{\lambda\nu}\left( \tan\left( \frac{1}{2}\, \sqrt{\lambda\, \nu}\vartheta\right) -\cot\left( \frac{1} {2}\, \sqrt{\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}. \end{equation}$

(3.21)

Family 1.4: For $\lambda\nu < 0$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{1, 16} = -\frac{3\, \kappa\, }{\varpi}\sqrt{-\lambda\nu}\left( \tanh\left( \sqrt{-\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.22)

$\begin{equation} \psi_{1, 17} = -\frac{3\, \kappa\, }{\varpi}\sqrt{-\lambda\nu}\left( \coth\left( \sqrt{-\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.23)

$\begin{equation} \psi_{1, 18} = -\frac{3\, \kappa\, }{\varpi}\sqrt{-\lambda\nu}\left( \tanh\left( 2\, \sqrt{-\lambda\, \nu}\vartheta\right) +i\mathit{sech}\left( 2\, \sqrt {-\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.24)

$\begin{equation} \psi_{1, 19} = -\frac{3\, \kappa\, }{\varpi}\sqrt{-\lambda\nu}\left( \coth\left( 2\, \sqrt{-\lambda\, \nu}\vartheta\right) +\mathit{csch}\left( 2\, \sqrt {-\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}, \end{equation}$

(3.25)

and

$\begin{equation} \psi_{1, 20} = \frac{-6\, \kappa\, }{\varpi}\sqrt{-\lambda\nu}\left( \tanh\left( \frac{1}{2}\, \sqrt{-\lambda\, \nu}\vartheta\right) +\coth\left( \frac{1} {2}\, \sqrt{-\lambda\, \nu}\vartheta\right) \right) ^{-1}+p_{{0}}. \end{equation}$

(3.26)

Family 1.5: For $\nu = \lambda$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{1, 21} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \tan\left( \lambda \, \vartheta\right) }}+p_{{0}}, \end{equation}$

(3.27)

$\begin{equation} \psi_{1, 22} = -3\, {\frac{\kappa\, \lambda}{\varpi\, \cot\left( \lambda \, \vartheta\right) }}+p_{{0}}, \end{equation}$

(3.28)

$\begin{equation} \psi_{1, 23} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \left( \tan\left( 2\, \lambda\, \vartheta\right) +\sec\left( 2\, \lambda\, \vartheta\right) \right) }}+p_{{0}}, \end{equation}$

(3.29)

$\begin{equation} \psi_{1, 24} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \left( -\cot\left( 2\, \lambda\, \vartheta\right) -\csc\left( 2\, \lambda\, \vartheta\right) \right) }}+p_{{0}}, \end{equation}$

(3.30)

and

$\begin{equation} \psi_{1, 25} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \left( \frac{1}{2} \, \tan\left( \frac{1}{2}\, \lambda\, \vartheta\right) -\frac{1}{2} \, \cot\left( \frac{1}{2}\, \lambda\, \vartheta\right) \right) }}+p_{{0}}. \end{equation}$

(3.31)

Family 1.6: For $\nu = -\lambda$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{1, 26} = -3\, {\frac{\kappa\, \lambda}{\varpi\, \tanh\left( \lambda \, \vartheta\right) }}+p_{{0}}, \end{equation}$

(3.32)

$\begin{equation} \psi_{1, 27} = -3\, {\frac{\kappa\, \lambda}{\varpi\, \coth\left( \lambda \, \vartheta\right) }}+p_{{0}}, \end{equation}$

(3.33)

$\begin{equation} \psi_{1, 28} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \left( -\tanh\left( 2\, \lambda\, \vartheta\right) -i\mathit{sech}\left( 2\, \lambda\, \vartheta \right) \right) }}+p_{{0}}, \end{equation}$

(3.34)

$\begin{equation} \psi_{1, 29} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \left( -\coth\left( 2\, \lambda\, \vartheta\right) -\mathit{csch}\left( 2\, \lambda\, \vartheta \right) \right) }}+p_{{0}}, \end{equation}$

(3.35)

and

$\begin{equation} \psi_{1, 30} = 3\, {\frac{\kappa\, \lambda}{\varpi\, \left( -\frac{1}{2} \, \tanh\left( \frac{1}{2}\, \lambda\, \vartheta\right) -\frac{1}{2} \, \coth\left( \frac{1}{2}\, \lambda\, \vartheta\right) \right) }}+p_{{0}}. \end{equation}$

(3.36)

Family 1.7: For $M = 0$ , we get

$\begin{equation} \psi_{1, 31} = -3/2\, {\frac{\kappa\, \lambda\, {\mu}^{2}\vartheta}{\varpi \, \lambda\left( \mu\, \vartheta+2\right) }}+p_{{0}}. \end{equation}$

(3.37)

Family 1.8: For $\mu = \tau$ , $\lambda = h\tau(h\neq0)$ and $\nu = 0$ ,

$\begin{equation} \psi_{1, 32} = {\frac{3\, \kappa\, h\tau+p_{{0}}\varpi\, {\mathrm{e}^{\tau \, \vartheta}}-p_{{0}}\varpi\, h}{\varpi\, \left( {\mathrm{e}^{\tau\, \vartheta} }-h\right) }}. \end{equation}$

(3.38)

Family 1.9: For $\mu = \nu = 0$ , we obtain the following:

$\begin{equation} \psi_{1, 33} = 3\, {\frac{\kappa}{\varpi\, \vartheta}}+p_{{0}}. \end{equation}$

(3.39)

According to case 1, the value of $\vartheta$ reads

$\vartheta = \frac{\kappa y^{\beta}}{\beta}-\frac{({\kappa}^{3}M)t^{\delta} }{\delta}.$

For case 2, and utilizing Eqs (3.1) and (3.4) with the corresponding solution of Eq (2.8), we construct the following families of traveling wave solutions for Eq (1.2):

Family 2.1: For $M < 0\quad\nu\neq0$ , we obtain the following:

$\begin{equation} \psi_{2, 1} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}+\frac{1} {2}\, {\frac{\sqrt{-M}\tan\left( \frac{1}{2}\, \sqrt{-M}\vartheta\right) } {\nu}}\right) , \end{equation}$

(3.40)

$\begin{equation} \psi_{2, 2} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {2}\, {\frac{\sqrt{-M}\cot\left( \frac{1}{2}\, \sqrt{-M}\vartheta\right) } {\nu}}\right) , \end{equation}$

(3.41)

$\begin{equation} \psi_{2, 3} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}+\frac{1} {2}\, {\frac{\sqrt{-M}\left( \tan\left( \sqrt{-M}\vartheta\right) +\sec\left( \sqrt{-M}\vartheta\right) \right) }{\nu}}\right) , \end{equation}$

(3.42)

$\begin{equation} \psi_{2, 4} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {2}\, {\frac{\sqrt{-M}\left( \cot\left( \sqrt{-M}\vartheta\right) +\csc\left( \sqrt{-M}\vartheta\right) \right) }{\nu}}\right) , \end{equation}$

(3.43)

and

$\begin{equation} \psi_{2, 5} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}+\frac{1} {4}\, {\frac{\sqrt{-M}\left( \tan\left( \frac{1}{4}\, \sqrt{-M}\vartheta \right) -\cot\left( 1/4\, \sqrt{-M}\vartheta\right) \right) }{\nu}}\right) . \end{equation}$

(3.44)

Family 2.2: For $M > 0\quad\nu\neq0$ , we obtain the following:

$\begin{equation} \psi_{2, 6} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {2}\, {\frac{\sqrt{M}\tanh\left( \frac{1}{2}\, \sqrt{M}\vartheta\right) }{\nu }}\right) , \end{equation}$

(3.45)

$\begin{equation} \psi_{2, 7} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {2}\, {\frac{\sqrt{M}\coth\left( \frac{1}{2}\, \sqrt{M}\vartheta\right) }{\nu }}\right) , \end{equation}$

(3.46)

$\begin{equation} \psi_{2, 8} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {2}\, {\frac{\sqrt{M}\left( \tanh\left( \sqrt{M}\vartheta\right) +i\mathit{sech}\left( \sqrt{M}\vartheta\right) \right) }{\nu}}\right) , \end{equation}$

(3.47)

$\begin{equation} \psi_{2, 9} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {2}\, {\frac{\sqrt{M}\left( \coth\left( \sqrt{M}\vartheta\right) +\mathit{csch}\left( \sqrt{M}\vartheta\right) \right) }{\nu}}\right) , \end{equation}$

(3.48)

and

$\begin{equation} \psi_{2, 10} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, {\frac{\mu}{\nu}}-\frac{1} {4}\, {\frac{\sqrt{M}\left( \tanh\left( \frac{1}{4}\, \sqrt{M}\vartheta \right) -\coth\left( 1/4\, \sqrt{M}\vartheta\right) \right) }{\nu}}\right) . \end{equation}$

(3.49)

Family 2.3: For $\lambda\nu > 0$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 11} = p_{{0}}+p_{{1}}\sqrt{{\frac{\lambda}{\nu}}}\tan\left( \sqrt{\lambda\, \nu}\vartheta\right) , \end{equation}$

(3.50)

$\begin{equation} \psi_{2, 12} = p_{{0}}-p_{{1}}\sqrt{{\frac{\lambda}{\nu}}}\cot\left( \sqrt{\lambda\, \nu}\vartheta\right) , \end{equation}$

(3.51)

$\begin{equation} \psi_{2, 13} = p_{{0}}+p_{{1}}\sqrt{{\frac{\lambda}{\nu}}}\left( \tan\left( 2\, \sqrt{\lambda\, \nu}\vartheta\right) +\sec\left( 2\, \sqrt{\lambda\, \nu }\vartheta\right) \right) , \end{equation}$

(3.52)

$\begin{equation} \psi_{2, 14} = p_{{0}}-p_{{1}}\sqrt{{\frac{\lambda}{\nu}}}\left( \cot\left( 2\, \sqrt{\lambda\, \nu}\vartheta\right) +\csc\left( 2\, \sqrt{\lambda\, \nu }\vartheta\right) \right) , \end{equation}$

(3.53)

and

$\begin{equation} \psi_{2, 15} = p_{{0}}+\frac{1}{2}\, p_{{1}}\sqrt{{\frac{\lambda}{\nu}}}\left( \tan\left( \frac{1}{2}\, \sqrt{\lambda\, \nu}\vartheta\right) -\cot\left( \frac{1}{2}\, \sqrt{\lambda\, \nu}\vartheta\right) \right) . \end{equation}$

(3.54)

Family 2.4: For $\lambda\nu < 0$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 16} = p_{{0}}-p_{{1}}\sqrt{-{\frac{\lambda}{\nu}}}\tanh\left( \sqrt{-\lambda\, \nu}\vartheta\right) , \end{equation}$

(3.55)

$\begin{equation} \psi_{2, 17} = p_{{0}}-p_{{1}}\sqrt{-{\frac{\lambda}{\nu}}}\coth\left( \sqrt{-\lambda\, \nu}\vartheta\right) , \end{equation}$

(3.56)

$\begin{equation} \psi_{2, 18} = p_{{0}}-p_{{1}}\sqrt{-{\frac{\lambda}{\nu}}}\left( \tanh\left( 2\, \sqrt{-\lambda\, \nu}\vartheta\right) +i\mathit{sech}\left( 2\, \sqrt {-\lambda\, \nu}\vartheta\right) \right) , \end{equation}$

(3.57)

$\begin{equation} \psi_{2, 19} = p_{{0}}-p_{{1}}\sqrt{-{\frac{\lambda}{\nu}}}\left( \coth\left( 2\, \sqrt{-\lambda\, \nu}\vartheta\right) +\mathit{csch}\left( 2\, \sqrt {-\lambda\, \nu}\vartheta\right) \right) , \end{equation}$

(3.58)

and

$\begin{equation} \psi_{2, 20} = p_{{0}}-\frac{1}{2}\, p_{{1}}\sqrt{-{\frac{\lambda}{\nu}}}\left( \tanh\left( \frac{1}{2}\, \sqrt{-\lambda\, \nu}\vartheta\right) +\coth\left( \frac{1}{2}\, \sqrt{-\lambda\, \nu}\vartheta\right) \right) . \end{equation}$

(3.59)

Family 2.5: For $\nu = \lambda$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 21} = p_{{0}}+p_{{1}}\tan\left( \lambda\, \vartheta\right) , \end{equation}$

(3.60)

$\begin{equation} \psi_{2, 22} = p_{{0}}-p_{{1}}\cot\left( \lambda\, \vartheta\right) , \end{equation}$

(3.61)

$\begin{equation} \psi_{2, 23} = p_{{0}}+p_{{1}}\left( \tan\left( 2\, \lambda\, \vartheta\right) +\sec\left( 2\, \lambda\, \vartheta\right) \right) , \end{equation}$

(3.62)

$\begin{equation} \psi_{2, 24} = p_{{0}}+p_{{1}}\left( -\cot\left( 2\, \lambda\, \vartheta\right) -\csc\left( 2\, \lambda\, \vartheta\right) \right) , \end{equation}$

(3.63)

and

$\begin{equation} \psi_{2, 25} = p_{{0}}+p_{{1}}\left( \frac{1}{2}\, \tan\left( \frac{1} {2}\, \lambda\, \vartheta\right) -\frac{1}{2}\, \cot\left( \frac{1}{2} \, \lambda\, \vartheta\right) \right) . \end{equation}$

(3.64)

Family 2.6: For $\nu = -\lambda$ and $\mu = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 26} = p_{{0}}-p_{{1}}\tanh\left( \lambda\, \vartheta\right) , \end{equation}$

(3.65)

$\begin{equation} \psi_{2, 27} = p_{{0}}-p_{{1}}\coth\left( \lambda\, \vartheta\right) , \end{equation}$

(3.66)

$\begin{equation} \psi_{2, 28} = p_{{0}}+p_{{1}}\left( -\tanh\left( 2\, \lambda\, \vartheta\right) -i\mathit{sech}\left( 2\, \lambda\, \vartheta\right) \right) , \end{equation}$

(3.67)

$\begin{equation} \psi_{2, 29} = p_{{0}}+p_{{1}}\left( -\coth\left( 2\, \lambda\, \vartheta\right) -\mathit{csch}\left( 2\, \lambda\, \vartheta\right) \right) , \end{equation}$

(3.68)

and

$\begin{equation} \psi_{2, 30} = p_{{0}}+p_{{1}}\left( -\frac{1}{2}\, \tanh\left( \frac{1} {2}\, \lambda\, \vartheta\right) -\frac{1}{2}\, \coth\left( \frac{1} {2}\, \lambda\, \vartheta\right) \right) . \end{equation}$

(3.69)

Family 2.7: For $M = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 31} = p_{{0}}-2\, {\frac{p_{{1}}\lambda\left( \mu\, \vartheta+2\right) }{{\mu}^{2}\vartheta}}. \end{equation}$

(3.70)

Family 2.8: For $\lambda = \mu = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 32} = p_{{0}}-{\frac{p_{{1}}}{\nu\, \vartheta}}. \end{equation}$

(3.71)

Family 2.9: For $\lambda = 0$ , $\mu\neq0$ and $\nu\neq0$ , we obtain the following:

$\begin{equation} \psi_{2, 33} = p_{{0}}-{\frac{p_{{1}}\mu}{\nu\, \left( \cosh\left( \mu\, \vartheta\right) -\sinh\left( \mu\, \vartheta\right) +1\right) }}, \end{equation}$

(3.72)

and

$\begin{equation} \psi_{2, 34} = p_{{0}}-{\frac{p_{{1}}\mu\, \left( \cosh\left( \mu\, \vartheta \right) +\sinh\left( \mu\, \vartheta\right) \right) }{\nu\, \left( \cosh\left( \mu\, \vartheta\right) +\sinh\left( \mu\, \vartheta\right) +1\right) }}. \end{equation}$

(3.73)

Family 2.10: For $\mu = \tau$ , $\nu = h\tau(h\neq0)$ and $\lambda = 0$ , we obtain the following:

$\begin{equation} \psi_{2, 35} = p_{{0}}+{\frac{p_{{1}}\mathrm{e^{\tau\, \vartheta}}} {1-h\mathrm{e^{\tau\, \vartheta}}}}. \end{equation}$

(3.74)

According to case 2, the value of $\vartheta$ reads obtain the following:

$\vartheta = \frac{({\frac{-\varpi\, p_{{1}}}{3\nu}})y^{\beta}}{\beta} -\frac{({\frac{-{\varpi}^{3}{p_{{1}}}^{3}M}{27{\nu}^{3}}})t^{\delta}}{\delta}.$

4. Discussion and graphs

This section provides graphic representations of various wave patterns for the model under consideration. The obtained solutions for different wave structures have been numerically examined in 3D and 2D forms using suitable values of the related parameters. We show that the obtained soliton solutions appear as periodic and kink wave solitons. Carefully choosing the parameters results in unique and educational visuals. Furthermore, the study's results are novel because, to the best of our knowledge, there is no prior literature that utilized this technique's results and applied them to the CEMCE. Moreover, it is essential to note that our analytical findings can be related to those derived from other approaches, such as the Riccati ODE expanding techniques, the F-expansion method, the tanh-method, and other related approaches, within specific limitations. Our work might be interpreted as a generalization of existing techniques, providing a more comprehensive framework from which these particular techniques can be drawn or comprehended.

The results obtained from the numerical investigations indicate that the soliton solutions developed can be classified into two main categories: Kink waves are distinguished by a restricted disruption or deflection in the wave pattern, which frequently retains its form throughout propagation; a periodic soliton is an autonomous wave propagating through a medium while preserving its shape and magnitude at regular intervals. Periodic, kink, and shock waves can display characteristics similar to wave breaking in shallow water waves. Wave breaking is a phenomenon that typically occurs near shorelines or other water barriers, and results in the sudden conversion of a wave into turbulent water. As a result, the generated soliton solutions provide valuable insights into the complex dynamics of the CEMCE and illustrate the fundamental principles that govern the interaction and propagation of the wave. By examining these solutions, academic researchers can enhance their comprehension of the model's performance in many contexts. These provide insights into phenomena such as the dispersion of configurations in liquid droplets and the behavior of waves in shallow water.

In this instance, we have conducted a graphical analysis of several derived aforementioned solutions to gain insights into the wave propagation mechanics and to enhance the reader's comprehension of the concept. For case (1), we numerically analyzed some of its solutions, as illustrated in –. In these figures, we studied and examined the effect of the fractional parameters $(\beta, \delta)$ on the wave profile described by these solutions. For instance, solution (3.12) was numerically analyzed, as shown in . The periodic wave solution (3.16) was analyzed, and the influence of the time-fractional parameter $\delta$ on the wave's behavior was investigated, as depicted in . Additionally, the impact of $\delta$ on the profile of the shock wave was investigated. Moreover, an investigation was conducted to examine the impact of $\delta$ on the characteristics of shock waves and solitons as specified by the solution (3.24). It is essential to acknowledge that the solution (3.24) exhibits complexity since it is comprised of both real and imaginary components. illustrates the shock wave, represented by the real part, and illustrates the solitons, represented by the imaginary part. Moreover, the absolute value of the solution (3.24) against $\delta$ was investigated, as shown in . The results obtained from the numerical analysis demonstrate that the time-fractional parameter exerts a substantial influence on the characteristics of both the shock waves and the soliton profiles. Additionally, another solution (3.26) for case (1) was examined, as depicted in . In this figure, we also studied the effect of $\delta$ on the profile of the shock waves described by solution (3.26).

Figure 1. The shock wave solution (3.12)

$\psi_{1, 6}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.7\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.85\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the solution (3.12) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 0.01,$

$\mu = 4,$

$\nu = 0.1,$

$\varpi = 1,$

$p_{{0}} = 10$ , and

$k = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. The periodic solution (3.16)

$\psi_{1, 10}$ is numerically examined: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.2\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the solution (3.16) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = \nu = 0.1,$

$\mu = 0,$

$\varpi = 1,$

$p_{{0}} = 10$ , and

$k = 10$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. The real part of solution (3.24) Re

$\left[\psi_{1, 18}\right]$ is numerically examined: (a) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.2\right)$ , (b) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (c) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the real part of shock wave solution (3.24) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = -2, \nu = 1,$

$\mu = 0,$

$\varpi = 3,$

$p_{{0}} = 0.1$ , and

$k = 0.25$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. The imaginary part of solution (3.24) Im

$\left[\psi_{1, 18}\right]$ is numerically examined: (a) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.2\right)$ , (b) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (c) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the imaginary part of shock wave solution (3.24) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = -2, \nu = 1,$

$\mu = 0,$

$\varpi = 3,$

$p_{{0}} = 0.1$ , and

$k = 0.25$ .

DownLoad: Full-Size Img PowerPoint

Figure 5. The absolute value of solution (3.24)

$\left\vert \psi_{1, 18}\right\vert$ is numerically examined: (a) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.2\right)$ , (b) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (c) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the absolute value of shock wave solution (3.24) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = -2, \nu = 1,$

$\mu = 0,$

$\varpi = 3,$

$p_{{0}} = 0.1$ , and

$k = 0.25$ .

DownLoad: Full-Size Img PowerPoint

Figure 6. Solution (3.26)

$\psi_{1, 20}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.1\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.2\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the solution (3.26) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = -2,$

$\mu = 0,$

$\nu = 1,$

$\varpi = 3,$

$p_{{0} } = 0.10$ , and

$k = 0.25$ .

DownLoad: Full-Size Img PowerPoint

To investigate the impact of the space-fractional parameter $\beta$ on the behavior of various solutions, we analyzed solutions (3.12), (3.16), (3.24), and (3.26), as depicted in . This diagram illustrates the susceptibility of these solutions to $\beta$ , which significantly influences the alteration of the waves behavior represented by these solutions. Furthermore, in this investigation, an examination was conducted on several solutions of case (2), as depicted in Figures 8–. For instance, we analyzed solution (3.45) and examined the impact of the time-fractional parameter on the properties of the shock waves described by this solution, as depicted in . Additionally, we analyzed the complex solution (3.47) and studied the effect of $\delta$ on the profile of the real and imaginary parts of this solution, as well as its absolute value, as illustrated in Figures 9–11, respectively. Moreover, we analyzed solutions (3.55), (3.65), (3.72), and (3.73) and examined the impact of the time-fractional parameter on the characteristics of shock waves, as depicted in Figures 12–15, respectively.

Figure 7. The impact of space-fractional parameter

$\beta$ on the profile of solutions (3.12), (3.16), (3.24), and (3.26) is numerically investigated: (a) solution (3.12) for

$\lambda = 0.01,$

$\mu = 4,$

$\nu = 0.1,$

$\varpi = 1,$

$p_{{0}} = 10$ ,

$k = 1$ ,

$y = 0.1$ , and

$\delta = 1$ , (b) solution (3.16) for

$\lambda = \nu = 0.1,$

$\mu = 0,$

$\varpi = 1,$

$p_{{0}} = 10$ ,

$k = 10$ ,

$y = 0.1$ , and

$\delta = 1$ , (c) solution (3.24) for

$\lambda = -2, \nu = 1,$

$\mu = 0,$

$\varpi = 3,$

$p_{{0}} = 0.1$ ,

$k = 0.25$ ,

$y = 0.1$ , and

$\delta = 1$ , (d) solution (3.26) for

$\lambda = -2,$

$\mu = 0,$

$\nu = 1,$

$\varpi = 3,$

$p_{{0}} = 0.10$ ,

$k = 0.25$ ,

$y = 0.1$ , and

$\delta = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 8. The shock wave solution (3.45)

$\psi_{2, 6}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.7\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.85\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the solution (3.45) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 1,$

$\mu = 4,$

$\nu = 1,$

$\varpi = 2,$ and

$p_{{0} } = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 9. The real part of solution (3.47) Re

$\left[\psi_{2, 8}\right]$ is numerically examined: (a) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the real part of shock wave solution (3.47) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 1, \nu = 1,$

$\mu = 4,$

$\varpi = 2,$

$p_{{0}} = 0.1$ , and

$p_{1} = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 10. The imaginary part of solution (3.47) Im

$\left[\psi_{2, 8}\right]$ is numerically examined: (a) 3D graph for solitons in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the imaginary part of solution (3.47) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 1, \nu = 1,$

$\mu = 4,$

$\varpi = 2,$

$p_{{0}} = 0.1$ , and

$p_{1} = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 11. The absolute value of solution (3.47)

$\left\vert \psi_{2, 8}\right\vert$ is numerically examined: (a) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph for shock waves in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the imaginary part of shock wave solution (3.47) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 1, \nu = 1,$

$\mu = 4,$

$\varpi = 2,$

$p_{{0}} = 0.1$ , and

$p_{1} = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 12. The shock wave solution (3.55)

$\psi_{2, 16}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the shock wave solution (3.55) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 1,$

$\mu = 0,$

$\nu = -1,$

$\varpi = 2,$

$p_{{0}} = 3$ , and

$p_{1} = 2$ .

DownLoad: Full-Size Img PowerPoint

Figure 13. The shock wave solution (3.65)

$\psi_{2, 26}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the shock wave solution (3.65) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 1,$

$\mu = 0,$

$\nu = -1,$

$\varpi = 2,$

$p_{{0}} = 3$ , and

$p_{1} = 2$ .

DownLoad: Full-Size Img PowerPoint

Figure 14. The shock wave solution (3.72)

$\psi_{2, 33}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the shock wave solution (3.72) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 0,$

$\mu = 1,$

$\nu = 1,$

$\varpi = 2,$

$p_{{0}} = 3$ , and

$p_{1} = 2$ .

DownLoad: Full-Size Img PowerPoint

Figure 15. The shock wave solution (3.73)

$\psi_{2, 34}$ is numerically investigated: (a) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.4\right)$ , (b) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 0.6\right)$ , (c) 3D graph in

$\left(y, t\right) -$ plane for

$\left(\beta, \delta\right) = \left(1, 1\right)$ , and (d) the shock wave solution (3.73) is plotted against

$y$ for different values of time-fractional parameter

$\delta$ and for

$\left(t, \beta\right) = \left(0.1, 1\right)$ . Here,

$\lambda = 0,$

$\mu = 1,$

$\nu = 1,$

$\varpi = 2,$

$p_{{0}} = 3$ , and

$p_{1} = 2$ .

DownLoad: Full-Size Img PowerPoint

Figure 16. The impact of space-fractional parameter

$\beta$ on the profile of solutions (3.55), (3.65), (3.72), and (3.73) is numerically investigated: (a) solution (3.55) for

$\lambda = 1,$

$\mu = 0,$

$\nu = -1,$

$\varpi = 2,$

$p_{{0}} = 3$ ,

$p_{1} = 2$ ,

$y = 0.1$ , and

$\delta = 1$ , (b) solution (3.65) for

$\lambda = 1,$

$\mu = 0,$

$\nu = -1,$

$\varpi = 2,$

$p_{{0} } = 3$ ,

$p_{1} = 2$ ,

$y = 0.1$ , and

$\delta = 1$ , (c) solution (3.72) for

$\lambda = 0,$

$\mu = 1,$

$\nu = 1,$

$\varpi = 2,$

$p_{{0}} = 3$ ,

$p_{1} = 2$ ,

$y = 0.1$ , and

$\delta = 1$ , (d) solution (3.73) for

$\lambda = 0,$

$\mu = 1,$

$\nu = 1,$

$\varpi = 2,$

$p_{{0}} = 3$ , and

$p_{1} = 2$ ,

$y = 0.1$ , and

$\delta = 1$ .

DownLoad: Full-Size Img PowerPoint

To see how $\beta$ changed the behavior of different solutions, we numerically examined the solutions (3.55), (3.65), (3.72), and (3.73) and saw how this parameter changed their behavior, as shown in . This figure shows the sensitivity of these solutions to $\beta$ , which has a strong effect on changing the behavior of the waves described by these solutions.

5. Conclusions and future work

The EDAM has been utilized to study and solve the nonlinear, space-time CEMCE. This approach used a suitable methodology to transform CEMCE into NODEs under the assumption of a finite series solution by utilizing Riccati ODEs. By solving the obtained equation, we received a series of different traveling wave solutions derived in the from of rational, exponential, trigonometric, and hyperbolic functions. We effectively and numerically analyzed some derived solutions using suitable values for the related parameters. Our analysis revealed that these solutions manifest as shock waves, solitons, periodic waves (cnoidal waves), and many different types of traveling wave solutions. The solutions obtained from the soliton simulations demonstrated the underlying principles that governed the interaction and propagation of the wave. These solutions offer valuable insights into the intricate dynamics of the CEMCE. By examining these solutions, scholars can enhance their understanding of the model's behavior in various settings, thus shedding light on issues such as wave dynamics in shallow water and the dispersion of configurations in liquid droplets. Therefore, the results of our study are of the utmost importance in addressing the practical challenges of studying different nonlinear structures that appear in shallow water. They contribute to our comprehension of the dynamics within these settings and have the potential to be utilized in predicting and manipulating essential fluid mechanics phenomena.

Furthermore, the suggested EDAM solves NFPDEs without requiring linearization or iterative techniques, and it is both user-friendly and efficient. While the EDAM has substantially increased our understanding of soliton dynamics and their impact on the models of interest, realizing the technique's limits is vital, mainly when the nonlinear term and the largest derivative are not homogeneously balanced. Despite this limitation, the work highlights the gaps in our understanding of soliton dynamics and nonlinear behaviors and presents novel opportunities for further investigations in future studies.

The obtained results demonstrated that the methodology used exhibited a notable abundance of solutions, thus enabling several researchers to employ this approach to model diverse physical and engineering issues. This method is anticipated to effectively analyze and model various evolution equations derived from fluid equations for diverse plasma systems. For instance, it can be utilized to analyze the family of fractional KdV-type equations ^[32,33], thus enabling an investigation into the fractional parameters' influence on the derived solutions' profile. Moreover, this method can examine the higher-order nonlinearity and derivative evolution equations that emerge in plasma, fluid mechanics, and nonlinear optics, such as the family of fractional Kawahara-type equations ^[34,35]. Additionally, we expect this method to successfully analyze fractional Schrödinger-type equations ^[36,37] and derive new modulated envelope solutions.

Author contributions

All authors of this article have been contributed equally. All authors have read and approved the final version of the manuscript for publication.

Use of AI tools declaration

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

Funding

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R17), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (KFU241071).

Acknowledgments

Conflict of interest

The authors declare that they have no conflicts of interest.

References

[1]	S. M. Zheng, Nonlinear evolution equations, New York: Chapman and Hall/CRC, 2004. https://doi.org/10.1201/9780203492222
[2]	R. Racke, Lectures on nonlinear evolution equations: Initial value problems, Cham: Birkhäuser, 2015. https://doi.org/10.1007/978-3-319-21873-1
[3]	C. Y. Zhu, M. Al-Dossari, S. Rezapour, S. Shateyi, On the exact soliton solutions and different wave structures to the modified Schrodinger's equation, Results Phys., 54 (2023), 107037. https://doi.org/10.1016/j.rinp.2023.107037 doi: 10.1016/j.rinp.2023.107037
[4]	C. Y. Zhu, M. Al-Dossari, N. S. A. El-Gawaad, S. A. M. Alsallami, S. Shateyi, Uncovering diverse soliton solutions in the modified Schrodinger's equation via innovative approaches, Results Phys., 54 (2023), 107100. https://doi.org/10.1016/j.rinp.2023.107100 doi: 10.1016/j.rinp.2023.107100
[5]	C. Y. Zhu, S. A. O. Abdallah, S. Rezapour, S. Shateyi, On new diverse variety analytical optical soliton solutions to the perturbed nonlinear Schrodinger equation, Results Phys., 54 (2023), 107046. https://doi.org/10.1016/j.rinp.2023.107046 doi: 10.1016/j.rinp.2023.107046
[6]	M. M. A. Khater, Physics of crystal lattices and plasma; analytical and numerical simulations of the Gilson-Pickering equation, Results Phys., 44 (2023), 106193. https://doi.org/10.1016/j.rinp.2022.106193 doi: 10.1016/j.rinp.2022.106193
[7]	M. M. Al-Sawalha, R. Shah, A. Khan, O. Y. Ababneh, T. Botmart, Fractional view analysis of Kersten-Krasil'shchik coupled KdV-mKdV systems with non-singular kernel derivatives, AIMS Mathematics, 7 (2022), 18334–18359. https://doi.org/10.3934/math.20221010 doi: 10.3934/math.20221010
[8]	H. Khan, S. Barak, P. Kumam, M. Arif, Analytical solutions of fractional Klein-Gordon and gas dynamics equations, via the (G'/G)-expansion method, Symmetry, 11 (2019), 566. https://doi.org/10.3390/sym11040566 doi: 10.3390/sym11040566
[9]	S. Behera, N. H. Aljahdaly, Nonlinear evolution equations and their traveling wave solutions in fluid media by modified analytical method, Pramana, 97 (2023), 130. https://doi.org/10.1007/s12043-023-02602-4 doi: 10.1007/s12043-023-02602-4
[10]	A. R. Adem, B. Muatjetjeja, T. S. Moretlo, An extended (2+1)-dimensional coupled burgers system in fluid mechanics: Symmetry reductions; Kudryashov method; conservation laws, Int. J. Theor. Phys., 62 (2023), 38. https://doi.org/10.1007/s10773-023-05298-9 doi: 10.1007/s10773-023-05298-9
[11]	K. J. Wang, Multi-wave complexiton, multi-wave, interaction-wave and the travelling wave solutions to the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli equation for the incompressible fluid, Pramana, 98 (2024), 47. https://doi.org/10.1007/s12043-024-02725-2 doi: 10.1007/s12043-024-02725-2
[12]	M. M. Bhatti, D. Q. Lu, An application of Nwogus Boussinesq model to analyze the head-on collision process between hydroelastic solitary waves, Open Phys., 17 (2019), 177–191. https://doi.org/10.1515/phys-2019-0018 doi: 10.1515/phys-2019-0018
[13]	J. H. He, X. H. Wu, Exp-function method for nonlinear wave equations, Chaos Soliton Fract., 30 (2006), 700–708. https://doi.org/10.1016/j.chaos.2006.03.020 doi: 10.1016/j.chaos.2006.03.020
[14]	J. F. Alzaidy, Fractional sub-equation method and its applications to the space-time fractional differential equations in mathematical physics, Brit. J. Math. Comput. Sci., 3 (2013), 153–163.
[15]	M. Cinar, A. Secer, M. Ozisik, M. Bayram, Derivation of optical solitons of dimensionless Fokas-Lenells equation with perturbation term using Sardar sub-equation method, Opt. Quant. Electron., 54 (2022), 402. https://doi.org/10.1007/s11082-022-03819-0 doi: 10.1007/s11082-022-03819-0
[16]	H. Yasmin, N. H. Aljahdaly, A. M. Saeed, R. Shah, Probing families of optical soliton solutions in fractional perturbed Radhakrishnan-Kundu-Lakshmanan model with improved versions of extended direct algebraic method, Fractal Fract., 7 (2023), 512. https://doi.org/10.3390/fractalfract7070512 doi: 10.3390/fractalfract7070512
[17]	M. Alqhtani, K. M. Saad, R. Shah, W. M. Hamanah, Discovering novel soliton solutions for (3+1)-modified fractional Zakharov-Kuznetsov equation in electrical engineering through an analytical approach, Opt. Quant. Electron., 55 (2023), 1149. https://doi.org/10.1007/s11082-023-05407-2 doi: 10.1007/s11082-023-05407-2
[18]	H. Yasmin, N. H. Aljahdaly, A. M. Saeed, R. Shah, Investigating families of soliton solutions for the complex structured coupled fractional Biswas-Arshed model in birefringent fibers using a novel analytical technique, Fractal Fract., 7 (2023), 491. https://doi.org/10.3390/fractalfract7070491 doi: 10.3390/fractalfract7070491
[19]	M. M. Al-Sawalha, H. Yasmin, R. Shah, A. H. Ganie, K. Moaddy, Unraveling the dynamics of singular stochastic solitons in stochastic fractional Kuramoto-Sivashinsky equation, Fractal Fract., 7 (2023), 753. https://doi.org/10.3390/fractalfract7100753 doi: 10.3390/fractalfract7100753
[20]	E. L. Mansfield, P. A. Clarkson, Symmetries and exact solutions for a 2+1-dimensional shallow water wave equation, Math. Comput. Simul., 43 (1997), 39–55. https://doi.org/10.1016/S0378-4754(96)00054-7 doi: 10.1016/S0378-4754(96)00054-7
[21]	W. Thadee, A. Chankaew, S. Phoosree, Effects of wave solutions on shallow-water equation, optical-fibre equation and electric-circuit equation, Maejo Int. J. Sci. Tech., 16 (2022), 262–274.
[22]	C. Y. Zhu, M. Al-Dossari, S. Rezapour, B. Gunay, On the exact soliton solutions and different wave structures to the (2+1) dimensional Chaffee-Infante equation, Results Phys., 57 (2024), 107431. https://doi.org/10.1016/j.rinp.2024.107431 doi: 10.1016/j.rinp.2024.107431
[23]	S. Lin, J. Zhang, C. Qiu, Asymptotic analysis for one-stage stochastic linear complementarity problems and applications, Mathematics, 11 (2023), 482. https://doi.org/10.3390/math11020482 doi: 10.3390/math11020482
[24]	Y. Kai, J. Ji, Z. Yin, Study of the generalization of regularized long-wave equation, Nonlinear Dyn., 107 (2022), 2745–2752. https://doi.org/10.1007/s11071-021-07115-6 doi: 10.1007/s11071-021-07115-6
[25]	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
[26]	W. Liu, X. Bai, H. Yang, R. Bao, J. Liu, Tendon driven bistable origami flexible gripper for high-speed adaptive grasping, IEEE Rob. Autom. Lett., 9 (2024), 5417–5424. https://doi.org/10.1109/LRA.2024.3389413 doi: 10.1109/LRA.2024.3389413
[27]	S. Phoosree, S. Chinviriyasit, New analytic solutions of some fourth-order nonlinear space-time fractional partial differential equations by G'/G-expansion method, Songklanakarin J. Sci. Technol., 43 (2021), 795–801.
[28]	S. Phoosree, W. Thadee, Wave effects of the fractional shallow water equation and the fractional optical fiber equation, Front. Appl. Math. Stat., 8 (2022), 900369. https://doi.org/10.3389/fams.2022.900369 doi: 10.3389/fams.2022.900369
[29]	V. E. Tarasov, On chain rule for fractional derivatives, Commun. Nonlinear Sci. Numer. Simul., 30 (2016), 1–4. https://doi.org/10.1016/j.cnsns.2015.06.007 doi: 10.1016/j.cnsns.2015.06.007
[30]	J. H. He, S. K. Elagan, Z. B. Li, Geometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus, Phys. Lett. A, 376 (2012), 257–259. https://doi.org/10.1016/j.physleta.2011.11.030 doi: 10.1016/j.physleta.2011.11.030
[31]	M. Z. Sarikaya, H. Budak, F. Usta, On generalized the conformable fractional calculus, TWMS J. Appl. Eng. Math., 9 (2019), 792–799.
[32]	S. A. Almutlak, S. Parveen, S. Mahmood, A. Qamar, B. M. Alotaibi, S. A. El-Tantawy, On the propagation of cnoidal wave and overtaking collision of slow shear Alfvén solitons in low $\beta$ -magnetized plasmas, Phys. Fluids, 35 (2023), 075130. https://doi.org/10.1063/5.0158292 doi: 10.1063/5.0158292
[33]	T. Hashmi, R. Jahangir, W. Masood, B. M. Alotaibi, S. M. E. Ismaeel, S. A. El-Tantawy, Head-on collision of ion-acoustic (modified) Korteweg-de Vries solitons in Saturn's magnetosphere plasmas with two temperature superthermal electrons, Phys. Fluids, 35 (2023), 103104. https://doi.org/10.1063/5.0171220 doi: 10.1063/5.0171220
[34]	B. S. Kashkari, S. A. El-Tantawy, A. H. Salas, L. S. El-Sherif, Homotopy perturbation method for studying dissipative nonplanar solitons in an electronegative complex plasma, Chaos Soliton Fract., 130 (2020), 109457. https://doi.org/10.1016/j.chaos.2019.109457 doi: 10.1016/j.chaos.2019.109457
[35]	S. A. El-Tantawy, A. M. Wazwaz, Anatomy of modified Korteweg-de Vries equation for studying the modulated envelope structures in non-Maxwellian dusty plasmas: Freak waves and dark soliton collisions, Phys. Plasmas, 25 (2018), 092105. https://doi.org/10.1063/1.5045247 doi: 10.1063/1.5045247
[36]	R. A. Alharbey, W. R. Alrefae, H. Malaikah, E. Tag-Eldin, S. A. El-Tantawy, Novel approximate analytical solutions to the nonplanar modified Kawahara equation and modeling nonlinear structures in electronegative plasmas, Symmetry, 15 (2023), 97. https://doi.org/10.3390/sym15010097 doi: 10.3390/sym15010097
[37]	S. A. El-Tantawy, A. H. Salas, H. A. Alyouse, M. R. Alharthi, Novel exact and approximate solutions to the family of the forced damped Kawahara equation and modeling strong nonlinear waves in a plasma, Chinese J. Phys., 77 (2022), 2454–2471. https://doi.org/10.1016/j.cjph.2022.04.009 doi: 10.1016/j.cjph.2022.04.009

This article has been cited by:

1.	Azzh Saad Alshehry, Safyan Mukhtar, Ali M. Mahnashi, Optical fractals and Hump soliton structures in integrable Kuralay-Ⅱ system, 2024, 9, 2473-6988, 28058, 10.3934/math.20241361
2.	Nauman Ahmed, Jorge E. Macías-Díaz, Shazia Umer, Muhammad Z. Baber, Muhammad Jawaz, Héctor Vargas-Rodríguez, On some solitary wave solutions of the Estevez--Mansfield--Clarkson equation with conformable fractional derivatives in time, 2024, 22, 2391-5471, 10.1515/phys-2024-0109

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.4

Metrics

Article views(1090) PDF downloads(61) Cited by(2)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(16)

AIMS Mathematics

Probing the diversity of soliton phenomena within conformable Estevez-Mansfield-Clarkson equation in shallow water

Related Papers:

Abstract

1. Introduction

2. Methodology & materials

2.1. The basic definition for the conformable fractional derivatives (CFDs)

2.2. The working procedure of the EDAM

3. Execution of EDAM to CEMCE

4. Discussion and graphs

5. Conclusions and future work

Author contributions

Use of AI tools declaration

Funding

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

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

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Mathematics

Probing the diversity of soliton phenomena within conformable Estevez-Mansfield-Clarkson equation in shallow water

Related Papers:

Abstract

1. Introduction

2. Methodology & materials

2.1. The basic definition for the conformable fractional derivatives (CFDs)

2.2. The working procedure of the EDAM

3. Execution of EDAM to CEMCE

4. Discussion and graphs

5. Conclusions and future work

Author contributions

Use of AI tools declaration

Funding

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

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

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog