Investigation of time fractional nonlinear KdV-Burgers equation under fractional operators with nonsingular kernels

Asif Khan; Tayyaba Akram; Arshad Khan; Shabir Ahmad; Kamsing Nonlaopon; Asif Khan; Tayyaba Akram; Arshad Khan; Shabir Ahmad; Kamsing Nonlaopon

doi:10.3934/math.2023063

AIMS Mathematics

2023, Volume 8, Issue 1: 1251-1268. doi: 10.3934/math.2023063

Previous Article Next Article

Research article

Investigation of time fractional nonlinear KdV-Burgers equation under fractional operators with nonsingular kernels

1.
Department of Mathematics, University of Malakand, Chakdara, Dir Lower, Khyber Pakhtunkhwa, Pakistan
2.
Department of Mathematics, COMSATS Institute of Information Technology, Lahore, 54000, Pakistan
3.
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand

Received: 14 July 2022 Revised: 30 September 2022 Accepted: 10 October 2022 Published: 19 October 2022
MSC : 35R11

In this manuscript, the Korteweg-de Vries-Burgers (KdV-Burgers) partial differential equation (PDE) is investigated under nonlocal operators with the Mittag-Leffler kernel and the exponential decay kernel. For both fractional operators, the existence of the solution of the KdV-Burgers PDE is demonstrated through fixed point theorems of $\alpha$ -type $\digamma$ contraction. The modified double Laplace transform is utilized to compute a series solution that leads to the exact values when fractional order equals unity. The effectiveness and reliability of the suggested approach are verified and confirmed by comparing the series outcomes to the exact values. Moreover, the series solution is demonstrated through graphs for a few fractional orders. Lastly, a comparison between the results of the two fractional operators is studied through numerical data and diagrams. The results show how consistently accurate the method is and how broadly applicable it is to fractional nonlinear evolution equations.

Keywords:

Citation: Asif Khan, Tayyaba Akram, Arshad Khan, Shabir Ahmad, Kamsing Nonlaopon. Investigation of time fractional nonlinear KdV-Burgers equation under fractional operators with nonsingular kernels[J]. AIMS Mathematics, 2023, 8(1): 1251-1268. doi: 10.3934/math.2023063

Related Papers:

[1]	Ihsan Ullah, Aman Ullah, Shabir Ahmad, Hijaz Ahmad, Taher A. Nofal . A survey of KdV-CDG equations via nonsingular fractional operators. AIMS Mathematics, 2023, 8(8): 18964-18981. doi: 10.3934/math.2023966
[2]	Rong Rong, Hui Liu . The Burgers-KdV limit in one-dimensional plasma with viscous dissipation: A study of dispersion and dissipation effects. AIMS Mathematics, 2024, 9(1): 1248-1272. doi: 10.3934/math.2024062
[3]	Musawa Yahya Almusawa, Hassan Almusawa . Numerical analysis of the fractional nonlinear waves of fifth-order KdV and Kawahara equations under Caputo operator. AIMS Mathematics, 2024, 9(11): 31898-31925. doi: 10.3934/math.20241533
[4]	Yanxia Hu, Qian Liu . On traveling wave solutions of a class of KdV-Burgers-Kuramoto type equations. AIMS Mathematics, 2019, 4(5): 1450-1465. doi: 10.3934/math.2019.5.1450
[5]	M. Mossa Al-Sawalha, Khalil Hadi Hakami, Mohammad Alqudah, Qasem M. Tawhari, Hussain Gissy . Novel Laplace-integrated least square methods for solving the fractional nonlinear damped Burgers' equation. AIMS Mathematics, 2025, 10(3): 7099-7126. doi: 10.3934/math.2025324
[6]	Shanhao Yuan, Yanqin Liu, Yibin Xu, Qiuping Li, Chao Guo, Yanfeng Shen . Gradient-enhanced fractional physics-informed neural networks for solving forward and inverse problems of the multiterm time-fractional Burger-type equation. AIMS Mathematics, 2024, 9(10): 27418-27437. doi: 10.3934/math.20241332
[7]	Hayman Thabet, Subhash Kendre, James Peters . Travelling wave solutions for fractional Korteweg-de Vries equations via an approximate-analytical method. AIMS Mathematics, 2019, 4(4): 1203-1222. doi: 10.3934/math.2019.4.1203
[8]	Jian-Gen Liu, Jian Zhang . A new approximate method to the time fractional damped Burger equation. AIMS Mathematics, 2023, 8(6): 13317-13324. doi: 10.3934/math.2023674
[9]	Yousef Jawarneh, Humaira Yasmin, Abdul Hamid Ganie, M. Mossa Al-Sawalha, Amjid Ali . Unification of Adomian decomposition method and ZZ transformation for exploring the dynamics of fractional Kersten-Krasil'shchik coupled KdV-mKdV systems. AIMS Mathematics, 2024, 9(1): 371-390. doi: 10.3934/math.2024021
[10]	Xiaoli Wang, Lizhen Wang . Traveling wave solutions of conformable time fractional Burgers type equations. AIMS Mathematics, 2021, 6(7): 7266-7284. doi: 10.3934/math.2021426

Abstract

1. Introduction

The Korteweg–de Vries–Burgers (KdV-Burgers) model, which comes up in several practical situations, such as the turbulence of undular bores in shallow water ^[1], the transport of liquids carrying gas bubbles ^[2], the waves that go through an elastic pipe that is filled with a viscous liquid ^[3], and weakly nonlinear plasma waves that have specific dissipative properties ^[4], has received a lot of interest over the last few decades. It may also be utilized as a nonlinear model in ferroelectricity theory, turbulence, circuit theory, and other fields ^[5,6]. The typical version of the KdV-Burger's equation is

$\begin{eqnarray} \frac{\partial \mathcal{U}}{\partial t}+\sigma\frac{\partial ^{3} \mathcal{U}}{\partial x^{3}}+\zeta \frac{\partial ^{2} \mathcal{U}}{\partial x^{2}}+\eta\mathcal{U}\frac{\partial\mathcal{U}}{\partial x} = 0,\; \; \; \; t > 0. \end{eqnarray}$

(1.1)

The KdV ^[7] and Burgers models ^[8] are commonly believed to be combined in Eq (1.1). Johnson ^[9] discovered that a specific limit of the matter led to the proposed model, where $\mathcal{U}(x, t)$ is proportional to the radial disturbance of the pipe wall, and $x$ and $t$ are the characterizing and temporal variables, in a field of wave propagation in fluid-filled elastic pipes. The model (1.1) was correct in the far-field of a near-field solution that was originally linear (small amplitude). Nonlinearity ( $\mathcal{U}\frac{\partial\mathcal{U}}{\partial x}$ ), dispersion ( $\frac{\partial ^{3} \mathcal{U}}{\partial x^{3}}$ ), and dissipation ( $\frac{\partial ^{2} \mathcal{U}}{\partial x^{2}}$ ) all exist in this equation, which is the basic version of a wave model. Fractional differential equations (FDEs) are extensions of differential equations (DEs) having integer order. FDEs have ample applications in different domains of sciences ^{[10,11,12,13,14]}. Due to the wide applications of FDEs, several operators have been defined in the literature ^[15,16]. The recent operators that are frequently used for studying DEs are the Caputo-Fabrizio (CF) ^[17] and Atangana-Baleanu (AB) ^[18] operators. These operators are dependent on the exponential and Mittag-Leffler kernels, respectively. The literature has several applications for the CF and the AB operators. For instance, HIV-1 infection has been investigated via the CF operator in ^[19]. Ahmad et al. studied the fractional-order Ambartsumian equation through the CF operator ^[20]. The $\Phi$ 4-model has been investigated using the CF and AB operators by Rahman et al. ^[21]. More applications can be found in the literature ^[22,23,24]. The fractional non-linear KdV-Burger's equation is taken into consideration in the form:

$\begin{eqnarray} \frac{\partial^{\alpha} \mathcal{U}}{\partial t^{\alpha}}+\sigma\frac{\partial ^{3} \mathcal{U}}{\partial x^{3}}+\zeta \frac{\partial ^{2} \mathcal{U}}{\partial x^{2}}+\eta\mathcal{U}\frac{\partial\mathcal{U}}{\partial x} = 0,\; \; \; 0 \; < \alpha \leq 1 ,\; \; t > 0, \end{eqnarray}$

(1.2)

with

$\mathcal{U}(x,0) = R(x),$

where $\mathcal{U}$ is a function of $x$ and $t$ , $x$ represents the space variable, $t$ represents the time variable and $\eta$ is a positive constant. We analyze Eq (1.2) in two ways: writing it first in the Atangana-Baleanu-Caputo (ABC) sense and then in the CF sense.

One of the most significant areas of research for FDEs is the quest for accurate and numerical solutions to FDEs. To date, many strategies for obtaining numerical and precise solutions of FDEs have been developed. A number of FDEs have been examined using these approaches. For example, Ahmad et al. ^[25] used the Laplace transform to find series of third order dispersive fractional PDEs. A generalized differential transform approach has been developed to solve fractional order PDEs by Odibat and Momani ^[26]. The homotopy perturbation technique has been proposed to solve the KdV-Burger's fractional PDE by Wang ^[27]. The Laplace transform was also observed to have a number of advantages, including its convergence to an exact solution of a problem after a certain iteration and that it does not allow any perturbation or discretization. Here, we utilize the double Laplace (DL) transform to compute a series solution of the considered equation.

The rest of article is organized as follows: Section 2 contains some basic definitions and a remark. The existence and uniqueness of the IVPs are presented in Section 3. In Section 4, the proposed techniques are presented. Section 5 consists of the application and comparison between results and diagrams of the proposed method. Finally, the conclusion is presented in Section 6.

2. Preliminaries

In this part, we provide some definitions, remarks, and lemmas about fractional calculus. Additionally, we provide a definition of the DL transformation and decomposition technique.

Definition 1. ^[18] Let $\mathcal{U}\in H^1(c, d)$ , $c > d$ and $\beta\in(0, 1]$ . Then the ABC operator is expressed as

$\begin{eqnarray*} ^{ABC}D^\beta \mathcal{U}(t) = \frac{B(\beta)}{(1-\beta)}\int_{c}^{t}\mathcal{U}'(s)E_{\beta}\left(\frac{-\beta(t-s)^\beta}{(1-\beta)}\right)\,ds, \label{eq2.1} \end{eqnarray*}$

where $B(\beta)$ is a normalizing factor with the conditions $B(0) = B(1) = 1$ , and $E_{\beta}$ is the Mittag-Leffler function

$\begin{equation*} E_{\beta}(x) = \sum\limits^{\infty}_{k = 0}\frac{x^{k}}{\Gamma(\beta k+1)},\; \; \; \; 0 < \beta < 1. \end{equation*}$

Definition 2. ^[17] Let $\mathcal{U}\in H^1(c, d)$ , $d > c$ and $\beta\in(0, 1]$ . Then the CF operator is written as

$\begin{eqnarray*} ^{CF}D^\beta \mathcal{U}(t) = \frac{E(\beta)}{1-\beta}\int_{c}^{t}\mathcal{U}'(t)\exp\left(\frac{-\beta(t-s)}{1-\beta}\right)\,ds, \label{eq2.2} \end{eqnarray*}$

where $E(\beta)$ is the normalizing factor such that $E(0) = E(1) = 1$ . When $\mathcal{U}(t)$ $\notin$ $H^1(c, d)$ then the above equation can be written for $u\ \epsilon\, L^{-1}(-\infty, d)$ and any $\beta\ \epsilon\, (0, 1]$ as

$^{CF}D^\beta \mathcal{U}(t) = \frac{\beta E(\beta)}{1-\beta}\int_{-\infty}^{t}(\mathcal{U}'(t)-\mathcal{U}(s))-\exp\left(\frac{-\beta(t-s)}{1-\beta}\right)\,ds. \label{eq2.3}$

Remark 1. For the above definitions, $n = [\beta]+1, \ [\beta]$ is the greatest integer not greater than $\beta$ , and " $\Gamma$ " is the well-known gamma function that can be calculated as

$\Gamma(\beta) = \int_0^\infty e^{-s} s^{\beta-1} ds. \label{eq2.4}$

Definition 3. Suppose that $\mathcal{U}$ is a function for $x, t > 0$ . The DL transformation of $\mathcal{U}$ is expressed as ^[28]

$\begin{equation*} {\mathrm{L}}_{x}{\mathrm{L}}_{t}\left[\mathcal{U}\right] = \int_0^{\infty}e^{-px}\int_0^{\infty} e^{-s{t}}\,\mathcal{U}\,dt\,dx, \label{eq2.5} \end{equation*}$

where $p$ and $s$ are complex numbers.

Definition 4. Application of the DL transform on the ABC operator is as follows:

$\begin{equation*} {\mathrm{L}}_x{\mathrm{L}}_t\left\{^{ABC}D_x^\theta\,\mathcal{U}\right\} = \frac{B(\theta)}{(1-\theta)(p^\theta+\frac{\theta}{(1-\theta)})}\left[p^\theta\,\bar{\mathcal{U}}(p,s)-\sum\limits_{k = 0}^{n-1}p^{\theta-1-k} {\mathrm{L}}_t\left\{\frac{\partial^k{\mathcal{U}(0,t)}}{\partial{x^k}} \right\}\right], \label{eq2.6} \end{equation*}$

and

$\begin{equation*} {\mathrm{L}}_x{\mathrm{L}}_t\left\{^{ABC}D_t^\vartheta\,\mathcal{U}\right\} = \frac{B(\vartheta)}{(1-\vartheta)(s^\vartheta+\frac{\vartheta}{(1-\vartheta)})}\left[s^\vartheta\,\bar{\mathcal{U}}(p,s)-\sum\limits_{k = 0}^{m-1}s^{\vartheta-1-k} {\mathrm{L}}_x\left\{\frac{\partial^k{\mathcal{U}(x,0)}}{\partial{t^k}} \right\}\right], \label{eq2.7} \end{equation*}$

where, $n = [\theta]+1, m = [\vartheta]+1$ .

Definition 5. Application of the DL transformation on the CF operator is as follows:

$\begin{equation*} {\mathrm{L}}_x{\mathrm{L}}_t\left\{^{CF}D_x^{\theta+n}\,\mathcal{U}\right\} = \frac{E(\theta)}{p+(1-\theta)p}\left[p^{n+1}\,\bar{\mathcal{U}}(p,s)-\sum\limits_{i = 0}^{n}p^{n-i} {\mathrm{L}}_t\left\{\frac{\partial^i{\mathcal{U}(0,t)}}{\partial{x^i}} \right\}\right], \label{eq2.8} \end{equation*}$

and

$\begin{equation*} {\mathrm{L}}_x{\mathrm{L}}_t\left\{^{CF}D_t^{\vartheta+m}\,\mathcal{U}\right\} = \frac{E(\vartheta)}{s+(1-\vartheta)s}\left[s^{m+1}\,\bar{\mathcal{U}}(p,s)-\sum\limits_{i = 0}^{m}s^{m-i} {\mathrm{L}}_x\left\{\frac{\partial^i{\mathcal{U}(x,0)}}{\partial{t^i}} \right\}\right], \label{eq2.9} \end{equation*}$

where, $n = [\theta]+1, m = [\vartheta]+1$ .

From the interpretation provided above, it is clear that

$\begin{equation*} {\mathrm{L}}_x{\mathrm{L}}_t{\mathcal{U}(x)\,v(t)} = \bar{\mathcal{U}}(p)\bar{v}(s) = {\mathrm{L}}_x{\mathcal{U}(x)}{\mathrm{L}}_t{v(t)}. \label{eq2.10} \end{equation*}$

A complex double-integral formulation is used to represent the inverse DL transform ${\mathrm{L}}_x^{-1}{\mathrm{L}}_t^{-1}\{\bar{\mathcal{U}}\} = \mathcal{U}$ :

$\begin{equation*} {\mathrm{L}}_x^{-1}{\mathrm{L}}_t^{-1}\Big\{\bar{\mathcal{U}}\Big\} = \frac{1}{2\pi\,i}\int\limits_{a-i\infty}^{a+i\infty}e^{s t}\int\limits_{b-i\infty}^{b+i\infty}e^{px} \bar{\mathcal{U}}(p,s)\,dp\,ds, \label{eq2.11} \end{equation*}$

where, $\bar{ \mathcal{U}}\left(p, s\right)$ , is an analytic function for all $p$ and $s$ that are described in the region $Re$ $(p)\geq a$ and $Re$ $(s)\geq b$ , where $a, b \in \mathbb{R}$ to be chosen appropriately.

3. Existence of the initial value problems (IVPs)

The existence and uniqueness of the IVPs are studied in this part employing $\alpha$ -type $\digamma$ -contraction. For this purpose, assume that $(\mathcal{Z}, d)$ is a complete metric space, and $\daleth$ is the collection of strictly increasing functions $\digamma:\mathcal{R_{+}}\to\mathcal{R}$ having the following required characteristics:

● $\lim\limits_{n\to\infty}\digamma(c_{n}) = -\infty$ if and only if, for each $\{c_{n}\}$ , $\lim\limits_{n\to\infty}(c_{n}) = 0;$

● there exists $\upsilon \in (0, 1)$ such that $\lim\limits_{c\to0^{+}}c^{\upsilon}\digamma(c) = 0.$

Definition 6. ^[29] Let $\mathcal{Q}:\mathcal{Z}\to\mathcal{Z}$ be self mapping with $\alpha:\mathcal{Z}\times\mathcal{Z}\to[0, \infty)$ . If

$\begin{eqnarray*} \alpha( \mathcal{X}, \mathcal{W}) \geq 1&\Rightarrow & \alpha(\mathcal{Q} \mathcal{X},\mathcal{Q} \mathcal{W})\geq1, \end{eqnarray*}$

for all $\mathcal{X}, \mathcal{W}\in\mathcal{Z}$ , then $\mathcal{Q}$ is referred to as $\alpha$ -admissible.

Definition 7. ^[30] Suppose that $\mathcal{Q}:\mathcal{Z}\to\mathcal{Z}$ , $\alpha:\mathcal{Z}\times\mathcal{Z}\to\{-\infty\}\cup[0, \infty)$ , and there exists $\omega > 0$ such that

$\begin{eqnarray*} \omega+\alpha(\mathcal{X},\mathcal{W})\digamma(d(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W})) &\leq& \digamma(d(\mathcal{X},\mathcal{W})) \end{eqnarray*}$

for each $\mathcal{X}, \mathcal{W}\in\mathcal{Z}$ with $d(\mathcal{Q}\mathcal{X}, \mathcal{Q}\mathcal{W}) > 0$ . Then, $\mathcal{Q}$ is called an $\alpha$ -type $\digamma$ -contraction.

Theorem 1. ^[30] Let $(\mathcal{Z}, d)$ be a complete metric space and $\mathcal{Q}:\mathcal{Z}\to\mathcal{Z}$ be an $\alpha$ -type $\digamma$ -contraction such that

1) there exists $\mathcal{X}_{\circ}\in \mathcal{Z}$ such that $\alpha(\mathcal{X}_{\circ}, \mathcal{Q}\mathcal{X}_{\circ})\geq1;$

2) if there exist $\{\mathcal{X}_{n}\}\subseteq \mathcal{Z}$ with $\alpha(\mathcal{X}_{n}, \mathcal{X}_{n+1})\geq1$ and $\mathcal{X}_{n}\to \mathcal{X}$ , then $\alpha(\mathcal{X}_{n}, \mathcal{X})\geq1$ for all $n\in \mathcal{N};$

3) $\digamma$ is continuous.

Then, $\mathcal{Q}$ has a fixed point $\mathcal{X}^{*}\in \mathcal{Z}$ . Also for $\mathcal{X}_{\circ}\in \mathcal{Z}$ , the sequence $\{\mathcal{Q}^{n}\mathcal{X}_{\circ}\}_{n\in \mathcal{N}}$ is convergent to $\mathcal{X}_{\circ}$

Let $\mathcal{Z} = \mathcal{C}([0, 1]^{2}, \mathcal{R})$ , where $\mathcal{C}$ is the space of all continuous functions $\mathcal{X}: [0, 1]\times [0, 1]\to \mathcal{R}$ , and $d(\mathcal{X}(x, t), \mathcal{W}(x, t)) = \sup_{x, t\in [0, 1]}\{|\mathcal{X}(x, t)-\mathcal{W}(x, t)|\}$ . Then we can write the IVP (1.1) in the CF fractional derivative sense as

$\begin{eqnarray} ^{CF}D_t^\alpha\,\mathcal{X}(x,t)& = & \mathcal{F}(x,t,\mathcal{X}(x,t)) ,\, \,0 < \alpha\leq 1 , \end{eqnarray}$

(3.1)

with initial condition

$\begin{eqnarray*} \mathcal{X}(x,0) = g(x),\, \, \end{eqnarray*}$

where $\mathcal{F}(x, t, \mathcal{X}(x, t)) = -\sigma\mathcal{U}_{xxx}-\zeta\mathcal{U}_{xx}-\eta\mathcal{U}\mathcal{U}_{x}.$

The following theorem demonstrates the existence of a solution of the problem (3.1).

Theorem 2. There exists $\mathcal{G}:\mathcal{R}^{2}\to\mathcal{R}$ such that

1) $|\mathcal{F}(x, t, \mathcal{X})-\mathcal{F}(x, t, \mathcal{W})|\leq\frac{2-\gamma M(\gamma)}{2} e^{b}|\mathcal{X}(x, t)-\mathcal{W}(x, t)|$ for $(x, t)\in[0, 1]^{2}$ and $\mathcal{X}, \mathcal{W}\in \mathcal{R};$

2) there exists $\mathcal{X}_{1}\in \mathcal{Z}$ such that $\mathcal{G}(\mathcal{X_{1}}, \mathcal{Q}\mathcal{X}_{1})\geq0$ , where $\mathcal{Q}:\mathcal{Z}\to\mathcal{Z}$ is defined by

$\mathcal{Q}\mathcal{X} = \mathcal{X}_{\circ}+{}^{CF}_{0}I^{\gamma}\mathcal{F}(x,t,\mathcal{X}(x,t));$

3) for $\mathcal{X}, \mathcal{W}\in \mathcal{Z}, \; \mathcal{G}(\mathcal{X}, \mathcal{W})\geq0$ implies that $\mathcal{G}(\mathcal{Q}\mathcal{X}, \mathcal{Q}\mathcal{W})\geq0$ ;

4) $\{\mathcal{X}_{n}\}\subseteq\mathcal{Z}, \lim\limits_{n\to\infty}\mathcal{X}_{n} = \mathcal{X}$ , where $\mathcal{X}\in \mathcal{Z}$ and $\mathcal{G}(\mathcal{X}_{n}, \mathcal{X}_{n+1})\geq0$ implies that $\mathcal{G}(\mathcal{X}_{n}, \mathcal{X})\geq0,$ for all $n\in \mathcal{N}$ .

Then, there exists at least one fixed point of $\mathcal{Q}$ that is the solution of the given model (3.1).

Proof. To prove that $\mathcal{Q}$ has a fixed point, we consider

$\begin{eqnarray*} |\mathcal{Q}\mathcal{X}-\mathcal{Q}\mathcal{W}||\mathcal{Q}\mathcal{X}-\mathcal{Q}\mathcal{W}+1|& = &|{}^{CF}I[\mathcal{F}(x,t\mathcal{X})-\mathcal{F}(x,t\mathcal{W})]||{}^{CF}I[\mathcal{F}(x,t\mathcal{X})-\mathcal{F}(x,t\mathcal{W})]+1| \\ &\leq& \big(\frac{2(1-\gamma)}{(2-\gamma)\mathcal{M}(\gamma)}|\mathcal{F}(x,t\mathcal{X})-\mathcal{F}(x,t\mathcal{W})|\\&&+\frac{2\gamma}{(2-\gamma)\mathcal{M}(\gamma)}\int_{0}^{\tau}|\mathcal{F}(x,t\mathcal{X})-\mathcal{F}(x,t\mathcal{W})|d\tau\big)\\&&\times\big(\frac{2(1-\gamma)}{(2-\gamma)\mathcal{M}(\gamma)}|\mathcal{F}(x,t\mathcal{X})-\mathcal{F}(x,t\mathcal{W})|\\&&+\frac{2\gamma}{(2-\gamma)\mathcal{M}(\gamma)}\int_{0}^{\tau}|\mathcal{F}(x,t\mathcal{X})-\mathcal{F}(x,t\mathcal{W})|d\tau+1\big) \big) \\ &\leq&\big(\frac{2(1-\gamma)}{(2-\gamma)\mathcal{M}(\gamma)}.\frac{(2-\gamma\mathcal{M}(\gamma))}{2}e^{-b}|\mathcal{X}-\mathcal{W}|\\&&+\frac{2\gamma}{(2-\gamma)\mathcal{M}(\gamma)}.\frac{2-\gamma\mathcal{M}(\gamma)}{2}\int_{0}^{\tau}e^{-b}|\mathcal{X}-\mathcal{W}|d\tau\big)\\&&\times\big(\frac{2(1-\gamma)}{(2-\gamma)\mathcal{M}(\gamma)}.\frac{(2-\gamma\mathcal{M}(\gamma))}{2}e^{-b}|\mathcal{X}-\mathcal{W}|\\&&+\frac{2\gamma}{(2-\gamma)\mathcal{M}(\gamma)}.\frac{2-\gamma\mathcal{M}(\gamma)}{2}\int_{0}^{\tau}e^{-b}|\mathcal{X}-\mathcal{W}|d\tau+1\big) \\ &\leq& \big(e^{-b}\sup\limits_{x,t\in[0,1]}|\mathcal{X}(x,t)-\mathcal{W}(x,t)|\big)\big(e^{-b}\sup\limits_{x,t\in[0,1]}|\mathcal{X}(x,t)-\mathcal{W}(x,t)|+1\big) \\ & = &\big( e^{-b}d(\mathcal{X},\mathcal{W})\big)\big(e^{-b}d(\mathcal{X},\mathcal{W})+1\big)\\ & = &e^{-b}[e^{-b}(d(\mathcal{X},\mathcal{W}))^{2}+d(\mathcal{X},\mathcal{W})] \\ &\leq& e^{-b}[(d(\mathcal{X},\mathcal{W}))^{2}+d(\mathcal{X},\mathcal{W})]. \end{eqnarray*}$

Thus, for $\mathcal{X}, \mathcal{W}\in \mathcal{Z}$ with $\mathcal{G}(\mathcal{X}, \mathcal{W})\geq0,$ we obtain

$(d(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W}))^{2}+d(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W})$

$\leq e^{-b}[(d(\mathcal{X},\mathcal{W}))^{2}+d(\mathcal{X},\mathcal{W})].$

Taking the $\ln$ on both sides, we have

$\begin{eqnarray*} b+\ln[d(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W}))^{2}+d(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W}] &\leq& \ln[d(\mathcal{X},\mathcal{W}))^{2}+d(\mathcal{X},\mathcal{W}]. \end{eqnarray*}$

If $\mathcal{F}:[0, \infty)\to\mathcal{R}$ is defined by $\mathcal{F}(\mathfrak{u}) = \ln[\mathfrak{u}^{2}+\mathfrak{u}], \; \; \mathfrak{u} > 0$ , then $\mathcal{F}\in \delta$ .

Now, define $\alpha:\mathcal{Z}\times\mathcal{Z}\to \{-\infty\}\cup[0, \infty)$ as

$\begin{eqnarray*} \alpha(\mathcal{X},\mathcal{W}) & = & \begin{cases} 1, & {{\rm{if }}}\; \mathcal{G}(\mathcal{X}(x,t),\mathcal{W}(x,t))\geq0 \; \; for\; \; all\; \; x,t\in [0,1], \\ -\infty, & {{\rm{otherwise}}}. \end{cases} \end{eqnarray*}$

Then,

$\begin{eqnarray*} b+\alpha(\mathcal{X},\mathcal{W})\mathcal{F}\big(d(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W})\big) &\leq& \mathcal{F}(d(\mathcal{X},\mathcal{W})), \end{eqnarray*}$

for $\mathcal{X}, \mathcal{W}\in\mathcal{Z}$ with $d(\mathcal{Q}\mathcal{X}, \mathcal{Q}\mathcal{W}) > 0$ . Now, by $\mathcal{G}3,$

$\begin{eqnarray*} \alpha(\mathcal{X},\mathcal{W}) \geq 1\; \; &\Rightarrow &\; \; \mathcal{G}(\mathcal{X},\mathcal{W})\geq0\; \; \Rightarrow\; \; \mathcal{G}(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W})\geq0\\ &\Rightarrow & \alpha(\mathcal{Q}\mathcal{X},\mathcal{Q}\mathcal{W})\geq1, \end{eqnarray*}$

for all $\mathcal{X}, \mathcal{W}\in \mathcal{Z}$ . From $\mathcal{G}2$ , there exists $\mathcal{X}_{\circ}\in \mathcal{Z}$ such that $\alpha(\mathcal{X}_{\circ}, \mathcal{Q}\mathcal{X}_{\circ})\geq1.$ Therefore by $\mathcal{G}4$ and Theorem 1, there exists $\mathcal{X}^{*}\in \mathcal{Z}$ such that $\mathcal{X}^{*} = \mathcal{Q}\mathcal{X}^{*}$ . Hence, $\mathcal{X}^{*}$ is the solution of the problem (3.1)

Similarly, we can write the IVP (1.1) in the ABC sense as

$\begin{eqnarray} ^{ABC}D_t^\alpha\,\mathcal{X}(x,t)& = & \mathcal{F}(x,t,\mathcal{X}(x,t)) ,\, \,0 < \alpha\leq 1 , \end{eqnarray}$

(3.2)

with initial condition

$\begin{eqnarray*} \mathcal{X}(x,0) = g(x),\, \, \end{eqnarray*}$

where $\mathcal{F}(x, t, \mathcal{X}(x, t)) = -\sigma\mathcal{U}_{xxx}-\zeta\mathcal{U}_{xx}-\eta\mathcal{U}\mathcal{U}_{x}.$

The following theorem shows the existence of a solution of the problem (3.2).

Theorem 3. There exists $\mathcal{G}:\mathcal{R}^{2}\to\mathcal{R}$ such that

1) $|\mathcal{F}(x, t, \mathcal{X})-\mathcal{F}(x, t, \mathcal{W})|\leq\frac{\Gamma\gamma M(\gamma)}{(1-\gamma)\Gamma\gamma+1} e^{\frac{-b}{2}}|\mathcal{X}(x, t)-\mathcal{W}(x, t)|$ for $(x, t)\in[0, 1]^{2}$ and $\mathcal{X}, \mathcal{W}\in \mathcal{R};$

2) there exists $\mathcal{X}_{1}\in \mathcal{Z}$ such that $\mathcal{G}(\mathcal{Y_{1}}, \mathcal{Q}\mathcal{X}_{1})\geq0$ , where $\mathcal{Q}:\mathcal{Z}\to\mathcal{Z}$ is defined by

$\mathcal{Q}\mathcal{X} = \mathcal{X}_{\circ}+{}^{ABC}_{0}I^{\gamma}\mathcal{F}(x,t,\mathcal{X}(x,t));$

3) for $\mathcal{X}, \mathcal{W}\in \mathcal{Z}, \; \mathcal{G}(\mathcal{X}, \mathcal{W})\geq0$ implies that $\mathcal{G}(\mathcal{Q}\mathcal{X}, \mathcal{Q}\mathcal{W})\geq0$ ;

Then, there exists at least one fixed point of $\mathcal{Q}$ which is the solution of the problem (3.2).

Proof.

$\begin{array}{l} |\mathcal{Q}\mathcal{X}-\mathcal{Q}\mathcal{W}|^{2} = |{}^{AB}_{0}I^{\gamma}[\mathcal{F}(x,t,\mathcal{X}(x,t))-\mathcal{F}(x,t,\mathcal{W}(x,t))]|^{2} \\ \leq |\{\frac{1-\gamma}{\mathcal{M}(\gamma)}[\mathcal{F}(x,t,\mathcal{X})-\mathcal{F}(x,t,\mathcal{W})]+\\ \frac{\gamma}{\mathcal{M}(\gamma)}{}_{0}I^{\gamma}[ \mathcal{F}(x,t, \mathcal{X}(x,t))- \mathcal{F}(x,t, \mathcal{W}(x,t))]\}|^{2}\\ \leq\{\frac{1-\gamma}{\mathcal{M}(\gamma)}|\mathcal{F}(x,t,\mathcal{X})-\mathcal{F}(x,t,\mathcal{W})| \\+\frac{\gamma}{\mathcal{M}(\gamma)}{}_{0}I^{\gamma}| \mathcal{F}(x,t, \mathcal{X}(x,t))- \mathcal{F}(x,t, \mathcal{W}(x,t))|\}^{2}\\ \leq \{\frac{1-\gamma}{\mathcal{M}(\gamma)}.\frac{\mathcal{M}(\gamma)\Gamma\gamma}{(1-\gamma)\Gamma\gamma+1}e^{\frac{-b}{2}}\sqrt{| \mathcal{X}- \mathcal{W}|^{2}}\\ +\frac{\gamma}{\mathcal{M}(\gamma)}\frac{\mathcal{M}(\gamma)\Gamma\gamma}{(1-\gamma)\Gamma\gamma+1}{}_{0}I^{\gamma}1.e^{\frac{-b}{2}}\sqrt{| \mathcal{X}- \mathcal{W}|^{2}}\}^{2}\\ = \{\frac{\mathcal{M}(\gamma)\Gamma\gamma}{(1-\gamma)\Gamma\gamma+1}e^{\frac{-b}{2}}\sqrt{| \mathcal{X}- \mathcal{W}|^{2}}\}^{2}\\ \{\frac{1-\gamma}{\mathcal{M}(\gamma)}+\frac{\gamma}{\mathcal{M}(\gamma)\gamma\Gamma\gamma}\}^{2}\\ \leq\{\frac{\mathcal{M}(\gamma)\Gamma\gamma}{(1-\gamma)\Gamma\gamma+1}e^{\frac{-b}{2}}\sqrt{\sup\limits_{x,t\in[0,1]}| \mathcal{X}(x,t)- \mathcal{W}(x,t)|^{2}}\}^{2}\\ \{\frac{1-\gamma}{\mathcal{M}(\gamma)}+\frac{\gamma}{\mathcal{M}(\gamma)\gamma\Gamma\gamma}\}^{2}\\ = \{\frac{\mathcal{M}(\gamma)\Gamma\gamma}{(1-\gamma)\Gamma\gamma+1}e^{\frac{-b}{2}}\sqrt{d( \mathcal{X}, \mathcal{W})}\}^{2}\\ \{\frac{1-\gamma}{\mathcal{M}(\gamma)}+\frac{\gamma}{\mathcal{M}(\gamma)\Gamma\gamma}\}^{2}\\ = e^{-b}d( \mathcal{X}, \mathcal{W}). \end{array}$

Consequently,

$\begin{eqnarray*} d(\mathcal{Q} \mathcal{X},\mathcal{Q} \mathcal{W}) &\leq&e^{-b}d( \mathcal{X}, \mathcal{W}). \end{eqnarray*}$

Applying " $\ln$ '' on both sides, we have

$\ln(d(\mathcal{Q} \mathcal{X},\mathcal{Q} \mathcal{W}))\leq\ln (e^{-b}d( \mathcal{X}, \mathcal{W})),$

and

$b+\ln(d(\mathcal{Q} \mathcal{X},\mathcal{Q} \mathcal{W})\leq\ln (d( \mathcal{X}, \mathcal{W})).$

Let $\digamma:[0, \infty)\to \mathcal{R}$ be defined by $\digamma(\lambda) = \ln\lambda,$ where $\lambda > 0$ . Then it is easy to show that $\digamma\in \daleth$ .

Now, define $\alpha:\mathcal{Z}\times\mathcal{Z}\to\{-\infty\}\cup[0, \infty)$ by

$\alpha( \mathcal{X}, \mathcal{W}) = \begin{cases} 1, & {{\rm{if }}} \;\mathcal{G}( \mathcal{X}(x,t), \mathcal{W}(x,t))\geq0\; \; for\; all\; x,t \in [0,1], \\ -\infty, & {{\rm{otherwise}}}. \end{cases}$

Thus, $b+\alpha(\mathcal{X}, \mathcal{W})\digamma(d(\mathcal{Q} \mathcal{X}, \mathcal{Q} \mathcal{W}))\leq\digamma(d(\mathcal{X}, \mathcal{W}))$ for $\mathcal{X}, \mathcal{W}\in \mathcal{Z}$ with $d(\mathcal{Q} \mathcal{X}, \mathcal{Q} \mathcal{W})\geq0$ . Therefore, $\mathcal{Q}$ is an $\alpha$ -type $\digamma$ -contraction. From $(\mathcal{G}3)$ , we have

$\begin{eqnarray*} \alpha( \mathcal{X}, \mathcal{W})\geq1 &\Rightarrow & \mathcal{G}( \mathcal{X}, \mathcal{W})\geq0\; \; \Rightarrow\; \; \mathcal{G}(\mathcal{Q} \mathcal{X},\mathcal{Q} \mathcal{W})\\ &\Rightarrow & \alpha(\mathcal{Q} \mathcal{X},\mathcal{Q} \mathcal{W})\geq1, \end{eqnarray*}$

for all $x, t\in[0, 1]$ . Thus, $\mathcal{Q}$ is $\alpha$ -admissible. From $(\mathcal{G}2)$ , there exists $\mathcal{X}_{\circ}\in \mathcal{Z}$ with $\alpha(\mathcal{X}_{\circ}, \mathcal{Q} \mathcal{X}_{\circ})\geq1.$ From $(\mathcal{G}4)$ and Theorem^[29], there exists $\mathcal{X}^{*}\in\mathcal{Z}$ such that $\mathcal{Q} \mathcal{X}^{*}$ . Hence, $\mathcal{X}^{*}$ is the solution of the IVP (3.2).

4. Modified double Laplace transform decomposition method (MDLDM)

Here, we briefly introduce the suggested technique MDLDM, which combines the DL transformation and Adomian decomposition method (ADM) for obtaining a series solution of non-linear ordinary DEs and PDEs. It is a very effective approach for obtaining the approximate values of dynamical problems such as KdV-Burgers, Sine-Gordon and KdV type equations. Here, we first introduce the technique and then the application of the given method on Eq (1.1).

Take the following form:

$\begin{equation} L\,\mathcal{U}+R\,\mathcal{U}+N\,\mathcal{U} = f(x,t),\; \; \; \forall\; \; t\in \mathbb{R}, \end{equation}$

(4.1)

where $L$ is linear, $R$ and $N$ are operators containing linear and non-linear terms, respectively, and $f(x, t)$ is some external function. $A_n$ , the well-known Adomian polynomials ^[31] of the functions $\mathcal{U}_0, \mathcal{U}_1, \mathcal{U}_2, \dots$ , can be described as

$\begin{eqnarray} A_n = \frac{1}{n!}\frac{d^n}{d\lambda^n}\left[\sum\limits_{k = 0}^n \lambda^k\mathcal{U}_k\mathcal{U}_{kx}\right]_{\lambda} = 0. \end{eqnarray}$

(4.2)

4.1. The suggested model in ABC sense

Here, first take Eq (1.1) in the form of ABC sense as:

$\begin{eqnarray} \ {^{ABC}D_{t}^{\alpha}\mathcal{U}}+\sigma\mathcal{U}_{xxx}+\zeta\mathcal{U}_{xx}+\eta\mathcal{U}\mathcal{U}_{x} = 0\; , \end{eqnarray}$

(4.3)

with initial conditions

$\mathcal{U}(x,0) = R(x). \label{eq4.2}$

By applying the DL transform method on both sides of Eq (4.3), we obtain

$\begin{eqnarray*} {\mathrm{L}}_{x}{\mathrm{L}}_{t}{^{ABC}D_{t}^{\alpha}\mathcal{U}} = -1\bigg[\sigma{\mathrm{L}}_{x}{\mathrm{L}}_{t}\mathcal{U}_{xxx}+\zeta{\mathrm{L}}_{x}{\mathrm{L}}_{t}\mathcal{U}_{xx} +\eta {\mathrm{L}}_{x}{\mathrm{L}}_{t}\mathcal{U}\mathcal{U}_{x}\bigg]. \label{eq4.3} \end{eqnarray*}$

Using initial condition and after some calculation, we have

$\begin{eqnarray*} {\mathrm{L}}_{x}{\mathrm{L}}_{t}\mathcal{U} = \frac{1}{s} {\mathrm{L}}_{x}R(x)-\bigg(1-\alpha+ \frac{\alpha}{s^{\alpha}}\bigg){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg[\sigma\mathcal{U}_{xxx}+\zeta\mathcal{U}_{xx}+\eta\mathcal{U}\mathcal{U}_{x}\bigg]. \label{eq4.4} \end{eqnarray*}$

Consider the series form

$\begin{eqnarray} \mathcal{U} = \sum\limits_{n = 0}^{\infty}\mathcal{U}_{n}. \end{eqnarray}$

(4.4)

The non-linear term $\mathcal{U}\mathcal{U}_{x}$ can be calculated using "ADM"; we obtain

$\mathcal{U}\mathcal{U}_{x} = \sum\limits_{n = 0}^{\infty}A_{n},$

$\begin{equation} \begin{aligned} {\mathrm{L}}_{x}{\mathrm{L}}_{t}\sum\limits_{n = 0}^{\infty}\mathcal{U}_{n} = \frac{1}{s} {\mathrm{L}}_{x}R(x)-\bigg(1-\alpha+ \frac{\alpha}{s^{\alpha}}\bigg)\bigg[{\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma\sum\limits_{n = 0}^{\infty} \mathcal{U}_{nxxx}+\zeta\sum\limits_{n = 0}^{\infty} \mathcal{U}_{nxx}+\eta\sum\limits_{n = 0}^{\infty}A_{n} \bigg)\bigg]. \end{aligned} \end{equation}$

(4.5)

Applying the double Laplace inverse transform and Equating terms of both sides,

$\begin{eqnarray*} \mathcal{U}_0& = &R(x),\\ \mathcal{U}_1& = &- {\mathrm{L}}_{x}^{-1}{\mathrm{L}}_{t}^{-1}\left(1-\alpha+\frac{\alpha}{s^{\alpha}}\right){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma\mathcal{U}_{0xxx}+\zeta\mathcal{U}_{0xx}+\eta A_{0}\bigg),\\ \mathcal{U}_2& = &- {\mathrm{L}}_{x}^{-1}{\mathrm{L}}_{t}^{-1}\left(1-\alpha+\frac{\alpha}{s^{\alpha}}\right){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma\mathcal{U}_{1xxx}+\zeta \mathcal{U}_{1xx}+\eta A_{1}\bigg),\\ \mathcal{U}_3& = &- {\mathrm{L}}_{x}^{-1}{\mathrm{L}}_{t}^{-1}\left(1-\alpha+\frac{\alpha}{s^{\alpha}}\right){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma \mathcal{U}_{2xxx}+\zeta \mathcal{U}_{2xx}+\eta A_{2}\bigg),\\ \\&\vdots. \label{eq4.6} \end{eqnarray*}$

The final solution can be calculated as follows:

$\begin{equation} \mathcal{U} = \sum\limits_{n = 0}^{\infty}\mathcal{U}_{n}. \end{equation}$

(4.6)

4.2. The proposed model in CF Form

Now, we take Eq (1.1) in CF fractional derivative sense as:

$\begin{equation} \begin{aligned} \ {^{CF}D_{t}^{\alpha}\mathcal{U}}+\sigma\mathcal{U}_{xxx}+\zeta\mathcal{U}_{xx}+\eta\mathcal{U}\mathcal{U}_{x} = 0\; , \end{aligned} \end{equation}$

(4.7)

with initial conditions

$\label{eq5.2} \mathcal{U}(x,0) = R(x).$

Using (3.1), we obtain the following series solutions:

$\begin{eqnarray*} \mathcal{U}_0& = & R(x).\\ \mathcal{U}_1& = &- {\mathrm{L}}_{x}^{-1}{\mathrm{L}}_{t}^{-1}\left(\frac{(1-s)\alpha+s}{s}\right){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma\mathcal{U}_{0xxx}+\zeta\mathcal{U}_{0xx}+\eta A_{0}\bigg),\\ \mathcal{U}_2& = &- {\mathrm{L}}_{x}^{-1}{\mathrm{L}}_{t}^{-1}\left(\frac{(1-s)\alpha+s}{s}\right){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma\mathcal{U}_{1xxx}+\zeta \mathcal{U}_{1xx}+\eta A_{1}\bigg),\\ \mathcal{U}_3& = &- {\mathrm{L}}_{x}^{-1}{\mathrm{L}}_{t}^{-1}\left(\frac{(1-s)\alpha+s}{s}\right){\mathrm{L}}_{x}{\mathrm{L}}_{t}\bigg(\sigma \mathcal{U}_{2xxx}+\zeta \mathcal{U}_{2xx}+\eta A_{2}\bigg),\\ &\vdots. \label{eq5.3} \end{eqnarray*}$

The final form can be expressed as

$\begin{equation} \mathcal{U} = \sum\limits_{n = 0}^\infty \mathcal{U}_n. \end{equation}$

(4.8)

5. Applications

Here, we consider a numerical example in the form of time fractional derivatives in ABC and CF operators.

Example 1. Consider Eq (1.1) in the following form:

(5.1)

with initial condition

$\begin{eqnarray*} \mathcal{U}_0 = {\rm{a}}_0-\frac{3 \zeta ^2 \tanh ^2\left(\frac{\zeta x}{10 \sigma }\right)}{25 \eta \sigma }+\frac{6 \zeta ^2 \tanh \left(\frac{\zeta x}{10 \sigma }\right)}{25 \eta \sigma }. \end{eqnarray*}$

The exact solution of Eq (5.1) is:

$\begin{array}{l} \mathcal{U} = {\rm{a}}_0+\frac{6 \zeta ^2}{25 \eta \sigma } \tanh \left(\frac{\zeta}{10 \sigma } \left(x+\frac{ \left(3 \zeta ^2-25 {\rm{a}}_0 \eta \sigma \right)t}{25 \sigma }\right)\right)\\-\frac{3 \zeta ^2}{25 \eta \sigma } \tanh^{2} \left(\frac{\zeta}{10 \sigma } \left(x+\frac{ \left(3 \zeta ^2-25 {\rm{a}}_0 \eta \sigma \right)t}{25 \sigma }\right)\right). \label{eq11} \end{array}$

Here, we will discuss two cases.

Case Ⅰ:

Consider Eq (5.1) in ABC form as follows:

$\begin{eqnarray} \ {^{ABC}D_{t}^{\alpha}\mathcal{U}}+\sigma\mathcal{U}_{xxx}+\zeta\mathcal{U}_{xx}+\eta\mathcal{U}\mathcal{U}_{x} = 0. \end{eqnarray}$

(5.2)

The approximate solutions of Eq (5.2) by using the techniques discussed in Section 4 are obtained in the series up to $\mathcal{O}(3)$ and given by

$\begin{eqnarray*} \mathcal{U}_0& = &{\rm{a}}_0-\frac{3 \zeta ^2 \tanh ^2\left(\frac{\zeta x}{10 \sigma }\right)}{25 \eta \sigma }+\frac{6 \zeta ^2 \tanh \left(\frac{\zeta x}{10 \sigma }\right)}{25 \eta \sigma },\\ \mathcal{U}_1& = &-\frac{3 \zeta ^3 \left(t^{\alpha}+\Gamma (\alpha )-\alpha\Gamma (\alpha )\right) \left(3 \zeta ^2-25 {\rm{a}}_0 \eta \sigma \right) \left(\tanh \left(\frac{\zeta x}{10 \sigma }\right)-1\right) {\rm{sech}}^2\left(\frac{\zeta x}{10 \sigma }\right)}{3125 \eta \sigma ^3 \Gamma (\alpha )},\\ \mathcal{U}_2& = &\frac{3 \zeta^{4}}{781250 \eta \sigma ^5 \Gamma(\alpha) \Gamma(1+2\alpha)} \left((t^{2\alpha}\alpha+2(\alpha-1)^2\Gamma(2\alpha))\Gamma(1+\alpha) -2t^{\alpha}(\alpha-1)\Gamma(1+2\alpha)\right)\\ &\; &\left(3 \zeta ^2-25 {\rm{a}}_0 \eta \sigma \right)^2 \left(\tanh \left(\frac{\zeta x}{10 \sigma }\right)-1\right) \left(3 \tanh \left(\frac{\zeta x}{10 \sigma }\right)+1\right) {\rm{sech}}^2\left(\frac{\zeta x}{10 \sigma }\right),\\ \end{eqnarray*}$

$\begin{array}{l} \mathcal{U}_3 = \frac{-3\zeta^{5}(3\zeta^{2}-25 {\rm{a}}_0 \eta \sigma)^{2}{\rm{sech}}^6\left(\frac{x\zeta}{10\sigma}\right)\left( \tanh \left(\frac{\zeta x}{10\sigma }\right)-1\right)}{781250000\sigma^{7}\eta(\Gamma (\alpha))^2\\\Gamma(1+\alpha)\Gamma(1+2\alpha)\Gamma(1+3\alpha)} [24\left(-1+\alpha\right)^{2}(\zeta)^{2}(\Gamma(\alpha))^{2}(t^{\alpha}\alpha-(-1+\alpha)\\\Gamma(1+\alpha))\Gamma(1+2\alpha)\Gamma(1+3\alpha) \left(\cosh\left(\frac{x\zeta}{10\sigma}\right)\\-\sinh\left(\frac{x\zeta}{10\sigma}\right)\right) \left(\cosh\left(\frac{x\zeta}{10\sigma}\right)+3\sinh\left(\frac{x\zeta}{10\sigma}\right)\right)\\ -24t^{2\alpha}\zeta^{2}\Gamma(1+\alpha)\Gamma(1+2\alpha) (t^{\alpha}\alpha\Gamma(1+2\alpha)\\-(-1+\alpha)\Gamma(1+3\alpha)) \left(-2+\cosh\left(\frac{x\zeta}{5\sigma}\right)-\sinh\left(\frac{x\zeta}{5\sigma}\right)\right)\\ +\Gamma(\alpha)\Gamma(1+\alpha)(2t^{\alpha}(-1+\alpha)^{2}\Gamma(1+3\alpha) (\alpha\Gamma(2\alpha)(-123\zeta^{2}\\+225 {\rm{a}}_0 \eta\sigma+4(9\zeta^{2}+25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{x\zeta}{5\sigma}\right)\\ +5(3\zeta^{2}-25 {\rm{a}}_0 \eta \sigma)\cosh\left(\frac{2x\zeta}{5\sigma}\right)-30(\zeta^{2}+5 {\rm{a}}_0 \eta\sigma)\sinh\left(\frac{x\zeta}{5\sigma}\right)\\+9\zeta^{2}\sinh\left(\frac{2x\zeta}{5\sigma}\right)\\- 75 {\rm{a}}_0 \eta\sigma\sinh\left(\frac{2x\zeta}{5\sigma}\right))\\ +\Gamma(1+2\alpha)(-75\zeta^{2}+225 {\rm{a}}_0 \eta\sigma+4(3\zeta^{2}+25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{x\zeta}{5\sigma}\right) \\+5(3\zeta^{2}-25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{2x\zeta}{5\sigma}\right)\\ -6(\zeta^{2}+25{\rm{a}}_0\eta\sigma)\sinh\left(\frac{x\zeta}{5\sigma}\right) +9\zeta^{2}\sinh\left(\frac{2x\zeta}{5\sigma}\right)\\- 75{\rm{a}}_0\eta\sigma\sinh\left(\frac{2x\zeta}{5\sigma}\right))) +\Gamma(1+\alpha)(t^{3\alpha}\alpha^{2}\Gamma(1+2\alpha)(-123\zeta^{2}\\ +225 {\rm{a}}_0\eta\sigma+4(9\zeta^{2}+25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{x\zeta}{5\sigma}\right)\\+5(3\zeta^{2}-25 {\rm{a}}_0\eta\sigma)\cosh\left(\frac{2x\zeta}{5\sigma}\right) -30(\zeta^{2}+5 {\rm{a}}_0\eta\sigma)\sinh\left(\frac{x\zeta}{5\sigma}\right)\\ +9\zeta^{2}\sinh\left(\frac{2x\zeta}{5\sigma}\right)- 75{\rm{a}}_0\eta\sigma\sinh\left(\frac{2x\zeta}{5\sigma}\right)) \\-(-1+\alpha)\Gamma(1+3\alpha)(3t^{2\alpha}\alpha(-91\zeta^{2} +225 {\rm{a}}_0\eta\sigma\\ +20(\zeta^{2}+5{\rm{a}}_0\eta\sigma)\cosh\left(\frac{x\zeta}{5\sigma}\right) \\+5(3\zeta^{2}-25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{2x\zeta}{5\sigma}\right) -14\zeta^{2}\sinh\left(\frac{x\zeta}{5\sigma}\right) -150{\rm{a}}_0\eta\sigma\sinh\left(\frac{x\zeta}{5\sigma}\right)\\ +9\zeta^{2}\sinh\left(\frac{2x\zeta}{5\sigma}\right)- 75{\rm{a}}_0\eta\sigma\sinh\left(\frac{2x\zeta}{5\sigma}\right)) \\+2(-1+\alpha)^{2}\Gamma(2\alpha)(-123\zeta^{2} +225 {\rm{a}}_0\eta\sigma\\ +4(9\zeta^{2}+25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{x\zeta}{5\sigma}\right) \\+5(3\zeta^{2}-25{\rm{a}}_0\eta\sigma)\cosh\left(\frac{2x\zeta}{5\sigma}\right) -30(\zeta^{2}+5 {\rm{a}}_0\eta\sigma)\sinh\left(\frac{x\zeta}{5\sigma}\right)\\ +9\zeta^{2}\sinh\left(\frac{2x\zeta}{5\sigma}\right)- 75{\rm{a}}_0\eta\sigma\sinh\left(\frac{2x\zeta}{5\sigma}\right)))))]\\ \vdots. \end{array}$

The simplified form of $\mathcal{U}_3$ is:

$\begin{array}{l} \mathcal{U}_3 = \frac{3\zeta^{5}(3\zeta^{2}-25 {\rm{a}}_0 \eta \sigma)^{2}{\rm{sech}}^4\left(\frac{x\zeta}{10\sigma}\right)\left( \tanh \left(\frac{\zeta x}{10\sigma }\right)-1\right)}{1562500000\sigma^{7}\eta}\\\Bigg[\frac{4}{3(1+\exp^{\frac{x\zeta}{5\sigma}})^2}\bigg(\exp^{\frac{-x\zeta}{5\sigma}}\bigg( 6(\alpha-1)^{3}(3(1+3\exp^{\frac{x\zeta}{5\sigma}}-25\exp^{\frac{2x\zeta}{5\sigma}}\\+\exp^{\frac{3x\zeta}{5\sigma}}+4\exp^{\frac{4x\zeta}{5\sigma}} )\zeta^{2}-25 {\rm{a}}_0 (1+\exp^{\frac{x\zeta}{5\sigma}})^{2}(1-7\exp^{\frac{x\zeta}{5\sigma}}+4\exp^{\frac{2x\zeta}{5\sigma}})\eta \sigma)+\\t^{2\alpha}\alpha\bigg(\frac{9(\alpha-1)}{\Gamma(2\alpha)}((3+17\exp^{\frac{x\zeta}{5\sigma}}-91\exp^{\frac{2x\zeta}{5\sigma}}\\ +3\exp^{\frac{3x\zeta}{5\sigma}}+12\exp^{\frac{4x\zeta}{5\sigma}})\zeta^{2}-25{\rm{a}}_0 (1+\exp^{\frac{x\zeta}{5\sigma}})^{2}(1-7\exp^{\frac{x\zeta}{5\sigma}}+4\exp^{\frac{2x\zeta}{5\sigma}})\eta \sigma)-\\\frac{2t^{\alpha}\alpha}{\Gamma(3\alpha)}(3(1+11\exp^{\frac{x\zeta}{5\sigma}}-41\exp^{\frac{2x\zeta}{5\sigma}} +\exp^{\frac{3x\zeta}{5\sigma}}\\ +4\exp^{\frac{4x\zeta}{5\sigma}})\zeta^{2} -25{\rm{a}}_0(1+\exp^{\frac{x\zeta}{5\sigma}})^{2}(1 -7\exp^{\frac{x\zeta}{5\sigma}}+4\exp^{\frac{2x\zeta}{5\sigma}})\eta \sigma)\bigg)\bigg)\\-\frac{48(2\exp^{\frac{x\zeta}{5\sigma}}-1)t^{2\alpha}\zeta^{2}(-3(\alpha-1)\Gamma(3\alpha)+t^{\alpha}\Gamma(1+2\alpha))} {(\Gamma(\alpha))^{2}\Gamma(3\alpha)}\bigg)\\ -\frac{24t^{\alpha}(\alpha-1)}{\Gamma(\alpha-1)}\bigg(-9\zeta^{2}+175{\rm{a}}_0\eta\sigma +5(3\zeta^{2}-25{\rm{a}}_0\eta\sigma)\cosh\bigg[\frac{x\zeta}{5\sigma}\bigg]\\-18\zeta^{2}sech^{2}\bigg[\frac{x\zeta}{5\sigma}\bigg] +3(3\zeta^{2}-25{\rm{a}}_0\eta\sigma)\sinh\bigg[\frac{x\zeta}{5\sigma}\bigg]\\ -12\zeta^{2}\tanh\bigg[\frac{x\zeta}{5\sigma}\bigg]\bigg)\Bigg]. \end{array}$

Therefore,

$\begin{eqnarray} \mathcal{U} = \mathcal{U}_0+\mathcal{U}_1+\mathcal{U}_2+\mathcal{U}_3..... \end{eqnarray}$

(5.3)

Case Ⅱ:

Consider Eq (1.2) in the CF sense as

$\begin{eqnarray} \ {^{CF}D_{t}^{\alpha}\mathcal{U}}+\sigma\mathcal{U}_{xxx}+\zeta\mathcal{U}_{xx}+\eta\mathcal{U}\mathcal{U}_{x} = 0\; . \end{eqnarray}$

(5.4)

The approximate solutions of Eq (5.4) obtained by using the approaches discussed in Section 4 up to $\mathcal{O}(3)$ are given by

$\begin{array}{l} \mathcal{U}_0 = = {\rm{a}}_0-\frac{\left(3 \zeta ^2\right) \tanh ^2\left(\frac{\zeta x}{10 \sigma }\right)}{25 \eta \sigma }+\frac{\left(6 \zeta ^2\right) \tanh \left(\frac{\zeta x}{10 \sigma }\right)}{25 \eta \sigma },\\ \mathcal{U}_1 = -\frac{3 \zeta ^3 (1+(t-1)\alpha) \left(3 \zeta ^2-25 {\rm{a}}_0 \eta \sigma \right) \left(\tanh \left(\frac{\zeta x}{10 \sigma }\right)-1\right) {\rm{sech}}^2\left(\frac{\zeta x}{10 \sigma }\right)}{3125 \eta \sigma ^3},\\ \mathcal{U}_2 =\\ \frac{3 \zeta ^4 \left(\alpha ^2 (t^2-4t+2)+4 \alpha (t-1)+2\right) \left(3 \zeta ^2-25 {\rm{a}}_0 \eta \sigma \right)^2 \left(\tanh \left(\frac{\zeta x}{10 \sigma }\right)-1\right) \left(3 \tanh \left(\frac{\zeta x}{10 \sigma }\right)+1\right) {\rm{sech}}^2\left(\frac{\zeta x}{10 \sigma }\right)}{1562500 \eta \sigma ^5},\\ \end{array}$

$\begin{array}{l} \mathcal{U}_3 = \frac{\zeta^{5}(3\zeta^{2}-25 {\rm{a}}_0 \eta \sigma)^{2}{\rm{sech}}^6(\frac{x\zeta}{10\sigma})}{1562500000 \eta \sigma^{7}}[4(3(-6-18(-1+t)\alpha\\-3(6-12t+t^{2})\alpha^{2}+ (6-18t+3t^{2}+t^{3})\alpha^{3})\zeta^{2}\\ -25 {\rm{a}}_0(6+18(-1+t)\alpha+9(2-4t+t^{2})\alpha^{2} +(-6+18t-9t^{2}+t^{3})\alpha^{3})\eta \sigma)\\\cosh(\frac{x\zeta}{5\sigma})-5(6+18(-1+t)\alpha\\ +9(2-4t+t^{2})\alpha^{2}\\+(-6+18t-9t^{2}+t^{3})\alpha^{3}) (3\zeta^{2}-25{\rm{a}}_0\eta \sigma)\\\cosh(\frac{2x\zeta}{5\sigma})-3(-3(50+150(-1+t)\alpha\\ +(150-300t+59t^{2})\alpha^{2}+(-50+150t-59t^{2}+3t^{3})\alpha^{3})\zeta^{2} \\+75 {\rm{a}}_0(6+18(-1+t)\alpha+9(2-4t+t^{2})\alpha^{2}\\ +(-6+18t-9t^{2}+t^{3})\alpha^{3})\eta \sigma +2(3(-2-6(-1+t)\alpha\\+(-6+12t+t^{2})\alpha^{2} +(2-6t-t^{2}+t^{2})\alpha^{3})\zeta^{2}\\ -25 {\rm{a}}_0(6+18(-1+t)\alpha+9(2-4t+t^{2})\alpha^{2} \\+(-6+18t-9t^{2}+t^{3})\alpha^{3})\eta \sigma)\sinh(\frac{x\zeta}{5\sigma})+(6+18(-1+t)\alpha\\ +9(2-4t+t^{2})\alpha^{2}+(-6+18t-9t^{2}+t^{3})\alpha^{3}) \\(3\zeta^{2}-25{\rm{a}}_0\eta\sigma)\sinh(\frac{2x\zeta}{5\sigma}))] (-1+\tanh(\frac{x\zeta}{10\sigma}))\\ \vdots. \end{array}$

The simplified form is

$\begin{array}{l} \mathcal{U}_3 = \frac{\zeta^{5}(3\zeta^{2}-25 {\rm{a}}_0 \eta \sigma)^{2}{\rm{sech}}^6(\frac{x\zeta}{10\sigma})(\tanh(\frac{x\zeta}{10\sigma})-1)}{1562500000 \eta \sigma^{7}} \\\Bigg[9\bigg(50+\alpha\bigg(150t(\alpha-1)^{2}-59t^{2}(\alpha-1)\alpha+3t^{3}\alpha^{2}\\ -50(3+(\alpha-3)\alpha)\bigg)\bigg)\zeta^{2} -225{\rm{a}}_0\bigg(6+18(t-1)\alpha\\+9(2+(t-4)t)\alpha^{2}+(-6+(t-6)(t-3)t)\alpha^{3}\bigg)\eta\sigma\\ +12\bigg(-6-18(t-1)\alpha-3(6+(t-12)t)\alpha^{2} \\+(6+(t-3)t(t+6))\alpha^{3}\bigg)\zeta^{2}\cosh\bigg[\frac{x\zeta}{5\sigma}\bigg]+\bigg(6+18(t-1)\alpha\\ +9(2+(t-4)t)\alpha^{2}+(-6+(t-6)(t-3)t)\alpha^{3}\bigg)\\(3\zeta^{2} -25{\rm{a}}_0\eta\sigma)\bigg(-5\cosh\bigg[\frac{2x\zeta}{5\sigma}\bigg]-3\sinh\bigg[\frac{2x\zeta}{5\sigma}\bigg]\bigg)\\ -18\bigg(-2+\alpha(6-6t+(-6+t(t+12))\alpha+(2+(t-3)t(t+2))\alpha^{2})\bigg)\\\zeta^{2}\sinh\bigg[\frac{x\zeta}{5\sigma}\bigg] +50{\rm{a}}_0\bigg(6+18(-1+t)\alpha\\ +9(2-4t+t^{2})\alpha^{2} +(-6+18t-9t^{2}+t^{3})\alpha^{3}\bigg)\sigma\eta\\\bigg(-2\cosh\bigg[\frac{x\zeta}{5\sigma}\bigg]+3\sinh\bigg[\frac{x\zeta}{5\sigma}\bigg]\bigg)\Bigg] \end{array}$

The final result in the series form up to $\mathcal{O}(3)$ is given by

$\begin{eqnarray} \mathcal{U} = \mathcal{U}_0+\mathcal{U}_1+\mathcal{U}_2+\mathcal{U}_3..... \end{eqnarray}$

(5.5)

The errors between the exact and approximate solutions in the ABC and CF senses are shown in the Tables 1 and 2, respectively.

Table 1. The comparison of approximated ABC and exact solutions for

$\alpha = 1$ at various points for Example 1.

(x, t)	$\mathcal{U}$	Exact	$\mid$ Exact $-$ $\mathcal{U}$ $\mid$	(x, t)	$\mathcal{U}$	Exact	$\mid$ Exact $-$ $\mathcal{U}$ $\mid$
( $-$ 4, 0.01)	$-$ 3.5967	$-$ 3.5967	4.437 $\times 10^{-11}$	( $-$ 2, 0.01)	$-$ 3.42484	$-$ 3.42484	1.488 $\times 10^{-9}$
( $-$ 1, 0.01)	$-$ 2.50239	$-$ 2.50239	7.202 $\times 10^{-9}$	(0, 0.01)	0.02862	0.02862	1.666 $\times 10^{-8}$
(1, 0.01)	1.13462	1.13462	4.376 $\times 10^{-9}$	(2, 0.01)	1.19852	1.19852	2.930 $\times 10^{-10}$
(4, 0.01)	1.2000	1.2000	1.184 $\times 10^{-13}$	( $-$ 4, 0.1)	$-$ 3.5959	$-$ 3.5959	4.635 $\times 10^{-7}$
( $-$ 2, 0.1)	$-$ 3.3840	$-$ 3.3840	1.515 $\times 10^{-5}$	( $-$ 1, 0.1)	$-$ 2.2941	$-$ 2.2942	7.761 $\times 10^{-5}$
(0, 0.1)	0.2693	0.2695	1.725 $\times 10^{-4}$	(1, 0.1)	1.1556	1.1556	4.302 $\times 10^{-5}$
(2, 0.1)	1.1990	1.1990	2.715 $\times 10^{-6}$	(4, 0.1)	1.2000	1.2000	1.085 $\times 10^{-9}$
( $-$ 4, 0.05)	$-$ 3.5963	$-$ 3.5963	2.827 $\times 10^{-8}$	( $-$ 2, 0.05)	$-$ 3.4077	$-$ 3.4077	9.378 $\times 10^{-7}$
( $-$ 1, 0.05)	$-$ 2.4131	$-$ 2.4131	4.656 $\times 10^{-6}$	(0, 0.05)	0.1395	0.1395	1.06 $\times 10^{-5}$
(1, 0.05)	1.1448	1.1448	2.716 $\times 10^{-6}$	(2, 0.05)	1.1987	1.1987	1.769 $\times 10^{-7}$
(4, 0.05)	1.2000	1.2000	7.099 $\times 10^{-11}$	( $-$ 4, 0.2)	$-$ 3.5948	$-$ 3.5948	7.795 $\times 10^{-6}$
( $-$ 2, 0.2)	$-$ 3.3282	$-$ 3.3279	2.4641 $\times 10^{-4}$	( $-$ 1, 0.2)	$-$ 2.0304	$-$ 2.0317	1.3403 $\times 10^{-3}$
(0, 0.2)	0.4958	0.4986	2.8191 $\times 10^{-3}$	(1, 0.2)	1.1720	1.1713	6.7111 $\times 10^{-4}$
(2, 0.2)	1.1994	1.1994	4.007 $\times 10^{-5}$	(4, 0.2)	1.2000	1.2000	1.5911 $\times 10^{-8}$

| Show Table

DownLoad: CSV

Table 2. The comparison of approximated the CF and exact solution for

$\alpha = 1$ at various points for Example 1.

(x, t)	$\mathcal{U}$	Exact	$\mid$ Exact $-$ $\mathcal{U}$ $\mid$	(x, t)	$\mathcal{U}$	Exact	$\mid$ Exact $-$ $\mathcal{U}$ $\mid$
( $-$ 4, 0.01)	$-$ 3.5967	$-$ 3.5967	4.437 $\times 10^{-11}$	( $-$ 2, 0.01)	$-$ 3.4248	$-$ 3.4248	1.488 $\times 10^{-9}$
( $-$ 1, 0.01)	$-$ 2.5023	$-$ 2.5023	7.202 $\times 10^{-9}$	(0, 0.01)	0.02862	0.02862	1.666 $\times 10^{-8}$
(1, 0.01)	1.13462	1.13462	4.376 $\times 10^{-9}$	(2, 0.01)	1.19852	1.19852	2.930 $\times 10^{-10}$
(4, 0.01)	1.2000	1.2000	1.184 $\times 10^{-13}$	( $-$ 4, 0.1)	$-$ 3.5959	$-$ 3.5959	4.635 $\times 10^{-7}$
( $-$ 2, 0.1)	$-$ 3.3840	$-$ 3.3840	1.515 $\times 10^{-5}$	( $-$ 1, 0.1)	$-$ 2.2941	$-$ 2.2942	7.761 $\times 10^{-5}$
(0, 0.1)	0.2693	0.2695	1.725 $\times 10^{-4}$	(1, 0.1)	1.1556	1.1556	4.302 $\times 10^{-5}$
(2, 0.1)	1.1990	1.1990	2.715 $\times 10^{-6}$	(4, 0.1)	1.2000	1.2000	1.085 $\times 10^{-9}$
( $-$ 4, 0.05)	$-$ 3.5963	$-$ 3.5963	2.827 $\times 10^{-8}$	( $-$ 2, 0.05)	$-$ 3.4077	$-$ 3.4077	9.378 $\times 10^{-7}$
( $-$ 1, 0.05)	$-$ 2.4131	$-$ 2.4131	4.656 $\times 10^{-6}$	(0, 0.05)	0.1395	0.1395	1.06 $\times 10^{-5}$
(1, 0.05)	1.1448	1.1448	2.716 $\times 10^{-6}$	(2, 0.05)	1.1987	1.1987	1.769 $\times 10^{-7}$
(4, 0.05)	1.2000	1.2000	7.099 $\times 10^{-11}$	( $-$ 4, 0.2)	$-$ 3.5948	$-$ 3.5948	7.795 $\times 10^{-6}$
( $-$ 2, 0.2)	$-$ 3.3282	$-$ 3.3279	2.4641 $\times 10^{-4}$	( $-$ 1, 0.2)	$-$ 2.0304	$-$ 2.0317	1.3403 $\times 10^{-3}$
(0, 0.2)	0.4958	0.4986	2.8191 $\times 10^{-3}$	(1, 0.2)	1.1720	1.1713	6.7111 $\times 10^{-4}$
(2, 0.2)	1.1994	1.1994	4.007 $\times 10^{-5}$	(4, 0.2)	1.2000	1.2000	1.5911 $\times 10^{-8}$

| Show Table

DownLoad: CSV

5.1. Discussion

Take the parameteric values as $\eta = 1,$ ${\rm{a}}_0 = 0$ and $\zeta$ and $\sigma = 0.1$ for numerical computations. Tables 1 and show the errors between the exact and approximated values in the ABC and CF senses, while the exact and computational values are also shown in these tables. depicts the comparison of approximated solutions in ABC and CF forms for different values of $\alpha$ at $t = 1$ . The surfaces in show the exact and approximated solutions of the ABC and the CF forms, respectively, for Eq (1.2) at $\alpha = 1$ . and present the approximate solutions when $\alpha = 0.95$ and $\alpha = 0.90$ , respectively. These graphical representations show the behaviour of ABC and CF operators respectively for the proposed problem. and depict the behaviours of surfaces for various values of $\alpha$ at different time levels. The first representation depicts the ABC form, while the second shows the CF form. It is straightforward to deduce that, as the fractional parameter $\alpha$ decreases, the wave response bifurcates into a wave but only for small values of $x$ , and we observe that the amplitude of the wave grows over time $t$ .

Figure 1. Comparison plots of approximated solutions of ABC and CF, respectively, for various values

$\alpha$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. Comparison plots of exact solutions, approximated solutions of ABC and CF for

$\alpha = 1$ respectively.

DownLoad: Full-Size Img PowerPoint

Figure 3. Comparison plots of approximated solutions of ABC and CF, respectively, for

$\alpha = 0.95$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. Comparison plots of approximated solutions of ABC and CF, respectively, for

$\alpha = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 5. Comparison plots of approximated solutions of ABC and CF, respectively, for

$\alpha = 0.5$ .

DownLoad: Full-Size Img PowerPoint

Figure 6. Comparison plots of approximated solutions of ABC and CF, respectively, for

$\alpha = 0.25$ .

DownLoad: Full-Size Img PowerPoint

6. Conclusions

In this manuscript, the time fractional KdV-Burger's model with initial conditions under two non-local operators with exponential kernel and Mittag-Leffler kernel has been investigated. The existence of the solution for both operators has been demonstrated through fixed point results of $\alpha$ -type $\digamma$ contraction. The MDLDM was utilized to compute a series solutions that tends to the exact values for the special case $\alpha = 1$ . As a result, we found highly accurate computed solutions to the fractional KdV-Burger's equation. The solutions show how consistently accurate the technique is and how broadly applicable it is to fractional nonlinear evolution problems. Quick convergence is seen when solutions are simulated numerically. Furthermore, neither linearization nor perturbation are needed. Therefore, it offers more authentic series solutions that typically converge quickly in actual physical problems.

Conflict of interest

The authors declare that they have no conflict of interest.

References

[1]	R. S. Johnson, Shallow water waves on a viscous fluid–the undular bore, Phys. Fluids, 15 (1972), 1693–1699. https://doi.org/10.1063/1.1693764 doi: 10.1063/1.1693764
[2]	L. van Wijngaarden, On the motion of gas bubbles in a perfect fluid, Annu. Rev. Fluid Mech., 4 (1972), 369–373. https://doi.org/10.1146/annurev.fl.04.010172.002101 doi: 10.1146/annurev.fl.04.010172.002101
[3]	R. S. Johnson, A nonlinear equation incorporating damping and dispersion, J. Fluid Mech., 42 (1970), 49–60. https://doi.org/10.1017/S0022112070001064 doi: 10.1017/S0022112070001064
[4]	P. N. Hu, Collisional theory of shock and nonlinear waves in a plasma, Phys. Fluids, 15 (1972), 854–864. https://doi.org/10.1063/1.1693994 doi: 10.1063/1.1693994
[5]	G. Gao, A theory of interaction between dissipation and dispersion of turbulence, Sci. China Ser. A, 28 (1985), 616–627.
[6]	S. D. Liu, S. K. Liu, KdV–Burgers equation modelling of turbulence, Sci. China Ser. A, 35 (1992), 576–586.
[7]	D. J. Korteweg, G. de Vries, On the change of form of long waves advancing in a rectangular channel, and on a new type of long stationary wave, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 39 (1895), 422–443. https://doi.org/10.1080/14786449508620739 doi: 10.1080/14786449508620739
[8]	J. M. Burgers, Mathematical examples illustrating relations occurring in the theory of turbulent fluid motion, Kon. Ned. Akad. Wet., 17 (1939), 1–53.
[9]	R. S. Johnson, A modern introduction to the mathematical theory of water waves, Cambridge: Cambridge University Press, 1997. https://doi.org/10.1017/CBO9780511624056
[10]	T. Akram, M. Abbas, M. B. Riaz, A. I. Ismail, N. M. Ali, Development and analysis of new approximation of extended cubic B-spline to the non-linear time fractional Klein-Gordon equation, Fractals, 28 (2020), 2040039. https://doi.org/10.1142/S0218348X20400393 doi: 10.1142/S0218348X20400393
[11]	T. Akram, M. Abbas, A. Ali, A numerical study on time fractional Fisher equation using an extended cubic B-spline approximation, J. Math. Comput. Sci., 22 (2020), 85–96. http://doi.org/10.22436/jmcs.022.01.08 doi: 10.22436/jmcs.022.01.08
[12]	R. T. Alqahtani, S. Ahmad, A. Akgül, Mathematical analysis of biodegradation model under nonlocal operator in Caputo sense, Mathematics, 9 (2021), 2787. https://doi.org/10.3390/math9212787 doi: 10.3390/math9212787
[13]	S. Ahmad, A. Ullah, A. Akgül, D. Baleanu, On solution of fractional model of human liver under hybrid fractional derivative, Discontinuity, Nonlinearity, and Complexity, 11 (2022), 183–190. https://doi.org/10.5890/DNC.2022.03.015 doi: 10.5890/DNC.2022.03.015
[14]	S. Saifullah, A. Ali, M. Irfan, K. Shah, Time-fractional Klein-Gordon equation with solitary/shock waves solutions, Math. Probl. Eng., 2021 (2021), 6858592. https://doi.org/10.1155/2021/6858592 doi: 10.1155/2021/6858592
[15]	T. Akram, M. Abbas, A. I. Izani, An extended cubic B-spline collocation scheme for time fractional sub-diffusion equation, AIP Conference Proceedings, 2184 (2019), 060017. https://doi.org/10.1063/1.5136449 doi: 10.1063/1.5136449
[16]	T. Akram, M. Abbas, A. I. Izani, Numerical solution of fractional cable equation via extended cubic B-spline, AIP Conference Proceedings, 2138 (2019), 030004. https://doi.org/10.1063/1.5121041 doi: 10.1063/1.5121041
[17]	M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Progress in Fractional Differentiation and Applications, 1 (2015), 73–85.
[18]	A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769. https://doi.org/10.2298/TSCI160111018A doi: 10.2298/TSCI160111018A
[19]	S. Ahmad, A. Ullah, M. Partohaghighi, S. Saifullah, A. Akgül, F. Jarad, Oscillatory and complex behaviour of Caputo-Fabrizio fractional order HIV-1 infection model, AIMS Mathematics, 7 (2022), 4778–4792. https://doi.org/10.3934/math.2022265 doi: 10.3934/math.2022265
[20]	S. Ahmad, A. Ullah, A. Akgül, M. De la Sen, A study of fractional order Ambartsumian equation involving exponential decay kernel, AIMS Mathematics, 6 (2021), 9981–9997. https://doi.org/10.3934/math.2021580 doi: 10.3934/math.2021580
[21]	F. Rahman, A. Ali, S. Saifullah, Analysis of time-fractional $\Phi$ 4-equation with singular and non-singular kernels, Int. J. Appl. Comput. Math., 7 (2021), 192. https://doi.org/10.1007/s40819-021-01128-w doi: 10.1007/s40819-021-01128-w
[22]	D. Kumar, J. Singh, D. Baleanu, On the analysis of vibration equation involving a fractional derivative with Mittag‐Leffler law, Math. Method. Appl. Sci., 43 (2020), 443–457. https://doi.org/10.1002/mma.5903 doi: 10.1002/mma.5903
[23]	M. Rahman, M. Arfan, W. Debani, P. Kumam, Z. Shah, Analysis of time-fractional Kawahara equation under Mittag-Leffler Power Law, Fractals, 30 (2022), 2240021. https://doi.org/10.1142/S0218348X22400217 doi: 10.1142/S0218348X22400217
[24]	K. A. Abro, I. Khan, K. S. Nisar, Use of Atangana–Baleanu fractional derivative in helical flow of a circular pipe, Fractals, 28 (2020), 2040049. https://doi.org/10.1142/S0218348X20400496 doi: 10.1142/S0218348X20400496
[25]	S. Ahmad, A. Ullah, K. Shah, A. Akgül, Computational analysis of the third order dispersive fractional PDE under exponential-decay and Mittag-Leffler type kernels, Numer. Meth. Part. Differ. Equ., in press. https://doi.org/10.1002/num.22627
[26]	Z. Odibat, S. Momani, A generalized differential transform method for linear partial differential equations of fractional order, Appl. Math. Lett., 21 (2008), 194–199. https://doi.org/10.1016/j.aml.2007.02.022 doi: 10.1016/j.aml.2007.02.022
[27]	Q. Wang, Homotopy perturbation method for fractional KdV-Burgers equation, Chaos Soliton. Fract., 35 (2008), 843–850. https://doi.org/10.1016/j.chaos.2006.05.074 doi: 10.1016/j.chaos.2006.05.074
[28]	I. N. Sneddon, The use of integral transforms, McGraw Hill, 1974.
[29]	B. Samet, C. Vetro, P. Vetro, Fixed point theorems for $\alpha$ - $\psi$ -contractive type mappings, Nonlinear Anal. Theor., 75 (2012), 2154–2165. https://doi.org/10.1016/j.na.2011.10.014 doi: 10.1016/j.na.2011.10.014
[30]	D. Gopal, M. Abbas, D. P. Kumar, C. Vetro, Fixed points of $\alpha$ -type $F$ -contractive mappings with an application to nonlinear fractional differential equation, Acta Math. Sci., 36 (2016), 957–970. https://doi.org/10.1016/S0252-9602(16)30052-2 doi: 10.1016/S0252-9602(16)30052-2
[31]	J. Duan, New recurrence algorithms for the nonclassic Adomian polynomials, Comput. Math. Appl., 62 (2011), 2961–2977. https://doi.org/10.1016/j.camwa.2011.07.074 doi: 10.1016/j.camwa.2011.07.074

This article has been cited by:

1.	Sajjad Ali, Aman Ullah, Shabir Ahmad, Kamsing Nonlaopon, Ali Akgül, Analysis of Kink Behaviour of KdV-mKdV Equation under Caputo Fractional Operator with Non-Singular Kernel, 2022, 14, 2073-8994, 2316, 10.3390/sym14112316
2.	Sadia Arshad, Iram Saleem, Ali Akgül, Jianfei Huang, Yifa Tang, Sayed M Eldin, A novel numerical method for solving the Caputo-Fabrizio fractional differential equation, 2023, 8, 2473-6988, 9535, 10.3934/math.2023481
3.	Vahisht K. Tamboli, Priti V. Tandel, Solution of the non-linear time-fractional Kudryashov–Sinelshchikov equation using fractional reduced differential transform method, 2024, 30, 1405-213X, 10.1007/s40590-024-00602-x
4.	Khan Shehzada, Aman Ullah, Sayed Saifullah, Ali Akgül, Fractional generalized perturbed KdV equation with a power Law kernel: A computational study, 2023, 12, 26667207, 100298, 10.1016/j.rico.2023.100298
5.	Amjad E. Hamza, Osman Osman, Muhammad Umair Sarwar, Khaled Aldwoah, Hicham Saber, Manel Hleili, Exploring Solitons Solutions of a (3+1)-Dimensional Fractional mKdV-ZK Equation, 2024, 8, 2504-3110, 498, 10.3390/fractalfract8090498
6.	Asif Khan, Abid Ullah Khan, Shabir Ahmad, Investigation of fractal fractional nonlinear Korteweg-de-Vries-Schrödinger system with power law kernel, 2023, 98, 0031-8949, 085202, 10.1088/1402-4896/ace08b
7.	Laila A. Al-Essa, Mati ur Rahman, A survey on fractal fractional nonlinear Kawahara equation theoretical and computational analysis, 2024, 14, 2045-2322, 10.1038/s41598-024-57389-0
8.	Sadia Yasin, Meraj Ali Khan, Shabir Ahmad, Saud Fahad Aldosary, Abundant new optical solitary waves of paraxial wave dynamical model with kerr media via new extended direct algebraic method, 2024, 56, 1572-817X, 10.1007/s11082-024-06845-2
9.	Wafa Alfwzan, Shao-Wen Yao, F.M. Allehiany, Shabir Ahmad, Sayed Saifullah, Mustafa Inc, Analysis of fractional non-linear tsunami shallow-water mathematical model with singular and non singular kernels, 2023, 52, 22113797, 106707, 10.1016/j.rinp.2023.106707
10.	Khalid K. Ali, M. Maneea, New approximation solution for time-fractional Kudryashov-Sinelshchikov equation using novel technique, 2023, 72, 11100168, 559, 10.1016/j.aej.2023.04.027
11.	Abdulrahman B. M. Alzahrani, Ghadah Alhawael, Novel Computations of the Time-Fractional Coupled Korteweg–de Vries Equations via Non-Singular Kernel Operators in Terms of the Natural Transform, 2023, 15, 2073-8994, 2010, 10.3390/sym15112010
12.	Sayed Saifullah, Nahid Fatima, Shaimaa A. M. Abdelmohsen, Meznah M. Alanazi, Shabir Ahmad, Dumitru Baleanu, Analysis of a conformable generalized geophysical KdV equation with Coriolis effect, 2023, 73, 11100168, 651, 10.1016/j.aej.2023.04.058
13.	Miguel Vivas-Cortez, Majeed A. Yousif, Bewar A. Mahmood, Pshtiwan Othman Mohammed, Nejmeddine Chorfi, Alina Alb Lupas, High-Accuracy Solutions to the Time-Fractional KdV–Burgers Equation Using Rational Non-Polynomial Splines, 2024, 17, 2073-8994, 16, 10.3390/sym17010016

Reader Comments

Your name:*

Email:*
© 2023 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(1954) PDF downloads(117) Cited by(13)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(6) / Tables(2)

AIMS Mathematics

Investigation of time fractional nonlinear KdV-Burgers equation under fractional operators with nonsingular kernels

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Existence of the initial value problems (IVPs)

4. Modified double Laplace transform decomposition method (MDLDM)

4.1. The suggested model in ABC sense

4.2. The proposed model in CF Form

5. Applications

5.1. Discussion

6. Conclusions

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

Investigation of time fractional nonlinear KdV-Burgers equation under fractional operators with nonsingular kernels

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Existence of the initial value problems (IVPs)

4. Modified double Laplace transform decomposition method (MDLDM)

4.1. The suggested model in ABC sense

4.2. The proposed model in CF Form

5. Applications

5.1. Discussion

6. Conclusions

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