Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators

M. Mossa Al-Sawalha; Osama Y. Ababneh; Rasool Shah; Amjad khan; Kamsing Nonlaopon; M. Mossa Al-Sawalha; Osama Y. Ababneh; Rasool Shah; Amjad khan; Kamsing Nonlaopon

doi:10.3934/math.2023120

AIMS Mathematics

2023, Volume 8, Issue 1: 2308-2336. doi: 10.3934/math.2023120

Previous Article Next Article

Research article

Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators

1.
Department of Mathematics, College of Science, University of Ha'il, Ha'il 2440, Saudi Arabia
2.
Department of Mathematics, Faculty of Science, Zarqa University, Zarqa 13110, Jordan
3.
Department of Mathematics, Abdul Wali khan university Mardan 23200, Pakistan
4.
Department of Mathematics and statistics, Bacha khan University, Palosa, Charsadda, Khyber Pakhtunkhwa, Pakistan
5.
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand

Received: 15 September 2022 Revised: 17 October 2022 Accepted: 19 October 2022 Published: 01 November 2022
MSC : 34A08, 35A20, 35R11

This paper solves a fractional system of non-linear Whitham-Broer-Kaup equations using a natural decomposition technique with two fractional derivatives. Caputo-Fabrizio and Atangana-Baleanu fractional derivatives were applied in a Caputo-manner. In addition, the results of the suggested method are compared to those of well-known analytical techniques such as the Adomian decomposition technique, the Variation iteration method, and the optimal homotopy asymptotic method. Two non-linear problems are utilized to demonstrate the validity and accuracy of the proposed methods. The analytical solution is then utilized to test the accuracy and precision of the proposed methodologies. The acquired findings suggest that the method used is very precise, easy to implement, and effective for analyzing the nature of complex non-linear applied sciences.

Keywords:

Citation: M. Mossa Al-Sawalha, Osama Y. Ababneh, Rasool Shah, Amjad khan, Kamsing Nonlaopon. Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators[J]. AIMS Mathematics, 2023, 8(1): 2308-2336. doi: 10.3934/math.2023120

Related Papers:

[1]	Maysaa Al Qurashi, Saima Rashid, Sobia Sultana, Fahd Jarad, Abdullah M. Alsharif . Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory. AIMS Mathematics, 2022, 7(7): 12587-12619. doi: 10.3934/math.2022697
[2]	Hussain Gissy, Abdullah Ali H. Ahmadini, Ali H. Hakami . The travelling wave phenomena of the space-time fractional Whitham-Broer-Kaup equation. AIMS Mathematics, 2025, 10(2): 2492-2508. doi: 10.3934/math.2025116
[3]	Faeza Lafta Hasan, Mohamed A. Abdoon, Rania Saadeh, Ahmad Qazza, Dalal Khalid Almutairi . Exploring analytical results for (2+1) dimensional breaking soliton equation and stochastic fractional Broer-Kaup system. AIMS Mathematics, 2024, 9(5): 11622-11643. doi: 10.3934/math.2024570
[4]	M. Mossa Al-Sawalha, Rasool Shah, Adnan Khan, Osama Y. Ababneh, Thongchai Botmart . Fractional view analysis of Kersten-Krasil'shchik coupled KdV-mKdV systems with non-singular kernel derivatives. AIMS Mathematics, 2022, 7(10): 18334-18359. doi: 10.3934/math.20221010
[5]	Ibrahim-Elkhalil Ahmed, Ahmed E. Abouelregal, Doaa Atta, Meshari Alesemi . A fractional dual-phase-lag thermoelastic model for a solid half-space with changing thermophysical properties involving two-temperature and non-singular kernels. AIMS Mathematics, 2024, 9(3): 6964-6992. doi: 10.3934/math.2024340
[6]	Ahu Ercan . Comparative analysis for fractional nonlinear Sturm-Liouville equations with singular and non-singular kernels. AIMS Mathematics, 2022, 7(7): 13325-13343. doi: 10.3934/math.2022736
[7]	Saima Rashid, Rehana Ashraf, Fahd Jarad . Strong interaction of Jafari decomposition method with nonlinear fractional-order partial differential equations arising in plasma via the singular and nonsingular kernels. AIMS Mathematics, 2022, 7(5): 7936-7963. doi: 10.3934/math.2022444
[8]	Thongchai Botmart, Ravi P. Agarwal, Muhammed Naeem, Adnan Khan, Rasool Shah . On the solution of fractional modified Boussinesq and approximate long wave equations with non-singular kernel operators. AIMS Mathematics, 2022, 7(7): 12483-12513. doi: 10.3934/math.2022693
[9]	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
[10]	Changdev P. Jadhav, Tanisha B. Dale, Vaijanath L. Chinchane, Asha B. Nale, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On solutions of fractional differential equations for the mechanical oscillations by using the Laplace transform. AIMS Mathematics, 2024, 9(11): 32629-32645. doi: 10.3934/math.20241562

Abstract

1. Introduction

In recent years, fractional calculus (FC) has been successfully applied to describe the mathematical issues in real materials ^[1,2]. For instance, ^[3] provided the representation of the constitutive relationship (RCR) of the fractional mechanical element with fractional derivative (FD) of the Riemann-Liouville type, and ^[4] presented its generalised form. In ^[5], the RCR with FC of the Liouville-Caputo type was taken into consideration. It has been demonstrated that the local fractional derivative is a wonderful mathematical tool for addressing fractal issues and some challenging natural phenomena ^[6,7]. The Caputo-Fabrizio type's RCR through fractional derivative was introduced in ^[8]. The study of fractional differential equations has expanded owing to the Caputo-Fabrizio fractional derivative. The nonsingular kernel of the new derivative ^[9], is what makes it so beautiful. The Caputo-Fabrizio derivative has the same additional motivating characteristics of heterogeneity and configuration ^[10,11] with different scales as it does in the Caputo and Riemann-Liouville fractional derivatives despite being created through the convolution of an ordinary derivative and an exponential function. In the past two years, numerous studies related to the the new Caputo-Fabrizio fractional derivative have been published. It has been shown that modelling using the fractional derivative of the Atangana-Baleanu has a brief random walk. Additionally, it has been found that the Mittag-Leffler function is a more significant and practical filter tool than the power and exponential law functions, making the Atangana-Baleanu fractional derivative in the sense of Caputo an effective mathematical tool for simulating more difficult real-world problems ^[12,13].

Using partial differential equations, it is possible to express a specific relationship between its partial derivatives (PDEs) and an unknown function. PDEs may be found in almost every field of engineering and research. In recent years, the application of PDEs in fields such as finance, biology, image graphics and processing, and social sciences has risen. As a result, when certain independent variables interact with one another in each of the above-mentioned fields, appropriate functions in these variables may be established, allowing for the modelling of a variety of processes via the use of equations for the related functions ^[14,15,16]. The study of PDEs has various elements. The traditional method, which dominated the nineteenth century, was to develop procedures for identifying explicit answers ^[17,18,19]. Theoretical PDE analysis offers a wide range of applications. It's worth noting that there are certain really difficult equations that even supercomputers can't solve. All one can do in these situations is attempt to get qualitative data. Moreover, the formulation of the equation, as well as its associated side conditions, is significant. A model of a physical or technical issue is frequently used to generate the equation. It is not immediately clear that the model is continuous in the view that it results to a solved PDE ^[20,21,22]. Furthermore, in most circumstances, it is preferable that the answer be one of a kind and stable under minor data disruptions. Theoretical comprehension of the equation aids in determining whether or not these requirements are satisfied. Several methods for solving classical PDEs have already been suggested, and several solutions have been discovered ^{[23,24,25,26,27,28,29]}.

The propagation of shallow water is represented by numerous well-known integral models, such as the Boussinesq equation, KdV equation, Whitham-Broer-Kaup equation, and others. Whitham, Broer, and Kaup ^{[30,31,32,33,34]} created non-linear Whitham-Broer-Kaup equations employing the Boussinesq approximation

$\begin{equation} \left\{\begin{array}{l} \psi_{\tau}+\psi \psi_{\upsilon}+\varphi_{\upsilon}+q \psi_{\upsilon\upsilon } = 0 \\ \varphi_{\tau}+\varphi \psi_{\upsilon}+\psi \varphi_{\upsilon}-q \varphi_{\upsilon\upsilon}+p \psi_{\upsilon\upsilon\upsilon} = 0, \end{array}\right. \end{equation}$

(1.1)

where $\psi = \psi(\upsilon, \tau), \varphi = \varphi(\upsilon, \tau)$ represents the velocity of horizontal and fluid height, which fluctuates substantially from equilibriums, and $q, p$ are constants made up of many diffusion power. Wang and Zheng ^[32] implemented riccati sub-equation method to obtain the result of fractional order Whitham-Broer-Kaup equations. The analytical and numerical methods have been applied to find the solutions of Whitham-Broer-Kaup equations, such as reduced differential transform method, residual power series method (1.1), finite element method ^[35], power-series method^[36], the finite difference approach ^[37], the exponential-function technique ^[38], the variation iteration technique, the homotopy perturbation technique, the homotopy analysis technique and others ^[39,40,41].

Adomian introduced the Adomian decomposition approach in 1980, which is an effective approach for locating numeric and explicit solutions to a wide variety of differential equations that describe physical conditions. This method is applicable to initial value problems, boundary value problems, partial and ordinary differential equations, including linear and nonlinear equations, and stochastic systems. The Natural transform decomposition method is constructed by combining the Adomian decomposition method with the Natural transform method (NTDM). NTDM has also been used to analyze fractional-order non-linear partial differential equations numerically in a number of articles ^[42,43,44].

In this research, we use NTDM in combination with two different derivatives to investigate the general as well as numerical solution of the coupled system of Whitham-Broer-Kaup equations of fractional order, as inspired by the above papers. NTDM is a simple and effective technique that does not require any perturbation. We compare our suggested method's results to those of other well-known approaches like the variational iteration method, Adomain decomposition method, and Optimal homoptoy analysis method. We can observe that the proposed method is superior to the previously discussed methods for finding nonlinear fractional order partial differential equation solutions. We use Maple to perform the calculations. The suggested techniques' convergence is also ensured by extending the concept mentioned in ^[45,46].

2. Basic definitions

Fractional derivatives and integrals have a great number of properties and definitions. We propose changes to some basic fractional calculus definitions and preliminaries used in this study.

Definition 2.1. The Riemann-Liouville integral of order fraction for a function $j\in C_{\wp}, \wp\geq-1$ , is given as ^[47,48,49]

$\begin{equation} \begin{split} I^{{{\varrho}}}j(\eta) = &\frac{1}{\Gamma({{\varrho}})}\int_{0}^{\eta}(\eta-\wp)^{{{\varrho}}-1}j(\wp)d\wp, \ \ {{\varrho}} > 0, \ \ \eta > 0.\\ and \ \ I^{0}j(\eta) = &j(\eta). \end{split} \end{equation}$

(2.1)

Definition 2.2. The derivative with order fraction for $j(\eta)$ in Caputo manner is stated as ^[47,48,49]

$\begin{equation} \begin{split} D_{\eta}^{{{\varrho}}}j(\eta) = I^{{m}-{{\varrho}}}D^{{m}}j(\eta) = \frac{1}{{m}-{{\varrho}}}\int_{0}^{\eta}(\eta-\wp)^{{m}-{{\varrho}}-1}j^{({m})}(\wp)d\wp, \end{split} \end{equation}$

(2.2)

for ${m}-1 < {{\varrho}}\leq {m}, \ \ {m}\in N, \ \ \eta > 0, j\in C^{{m}}_{\wp}, \wp\geq-1$ .

Definition 2.3. The derivative with order fraction for $j(\eta)$ in CF sense is stated as ^[47]

$\begin{equation} \begin{split} D_{\eta}^{{{\varrho}}}j(\eta) = \frac{F({{\varrho}})}{1-{{\varrho}}}\int^{\eta}_{0}\exp\left(\frac{-{{\varrho}}(\eta-\wp)}{1-{{\varrho}}}\right)D(j(\wp))d\wp, \end{split} \end{equation}$

(2.3)

having $0 < {{\varrho}} < 1$ and $F({{\varrho}})$ is a normalization function having $F(0) = F(1) = 1$ .

Definition 2.4. The derivative with order fraction for $j(\eta)$ in ABC sense is stated as ^[47]

$\begin{equation} \begin{split} D_{\eta}^{{{\varrho}}}j(\eta) = \frac{B({{\varrho}})}{1-{{\varrho}}}\int^{\eta}_{0}E_{{{\varrho}}}\left(\frac{-{{\varrho}}(\eta-\wp)}{1-{{\varrho}}}\right)D(j(\wp))d\wp, \end{split} \end{equation}$

(2.4)

having $0 < {{\varrho}} < 1$ , $B(\varrho)$ is normalization function and $E_{{{\varrho}}}(z) = \sum_{m = 0}^{\infty}\frac{z^{m}}{\Gamma({m}{{\varrho}}+1)}$ is the Mittag-Leffler function.

Definition 2.5. For a function $\psi(\tau)$ , the natural transform (NT) is stated as

$\begin{equation} \begin{split} \mathbf{N}(\psi(\tau)) = \mathcal{U}({\kappa}, {{\ell}}) = \int_{-\infty}^{\infty}e^{-\kappa\tau}\psi(\ell\tau)d\tau, \ \ \kappa , \ell\in (-\infty, \infty), \end{split} \end{equation}$

(2.5)

and for $\tau\in(0, \infty)$ , the NT of $\psi(\tau)$ is stated as

$\begin{equation} \begin{split} \mathbf{N}(\psi(\tau)H(\tau)) = \mathbf{N}^{+} = \mathcal{U}^{+}({\kappa}, {{\ell}}) = \int_{0}^{\infty}e^{-\kappa\tau}\psi(\ell\tau)d\tau, \ \ \kappa , \ell\in (0, \infty), \end{split} \end{equation}$

(2.6)

where $H(\tau)$ is the Heaviside function.

Definition 2.6. For the function $\psi({\kappa}, {{\ell}})$ , the inverse NT is stated as

$\begin{equation} \begin{split} \mathbf{N}^{-1}[\mathcal{U}({\kappa}, {{\ell}})] = \psi(\tau), \ \ \forall \tau\geq0. \end{split} \end{equation}$

(2.7)

Lemma 2.7. If the NT of $\psi_{1}(\tau)$ and $\psi_{2}(\tau)$ are $\mathcal{U}_{1}({\kappa}, {{\ell}})$ and $\mathcal{U}_{2}({\kappa}, {{\ell}})$ respectively, so

$\begin{equation} \begin{split} \mathbf{N}[c_{1}\psi_{1}(\tau)+c_{2}\psi_{2}(\tau)] = c_{1}\mathbf{N}[\psi_{1}(\tau)]+c_{2}\mathbf{N}[\psi_{2}(\tau)] = c_{1}\mathcal{U}_{1}({\kappa}, {{\ell}})+c_{2}\mathcal{U}_{2}({\kappa}, {{\ell}}), \end{split} \end{equation}$

(2.8)

having $c_{1}$ and $c_{2}$ constants.

Lemma 2.8. If $\psi_{1}(\tau)$ and $\psi_{2}(\tau)$ are the inverse NT of $\psi_{1}({\kappa}, {{\ell}})$ and $\psi_{2}({\kappa}, {{\ell}})$ respectively, so

$\begin{equation} \mathcal\{N\}^{-1}[c_{1}\mathcal{U}_{1}({\kappa}, {{\ell}})+c_{2}\mathcal{U}_{2}({\kappa}, {{\ell}})] = c_{1}\mathbf{N}^{-1}[\mathcal{U}_{1}({\kappa}, {{\ell}})]+c_{2}\mathbf{N}^{-1}[\mathcal{U}_{2}({\kappa}, {{\ell}})] = c_{1}\psi_{1}(\tau)+c_{2}\psi_{2}(\tau), \end{equation}$

(2.9)

having $c_{1}$ and $c_{2}$ constants.

Definition 2.9. The NT of $D_{\tau}^{{{\varrho}}}\psi(\tau)$ in Caputo sense is stated as ^[47]

$\begin{equation} \mathbf{N}[D_{\tau}^{{{\varrho}}}\psi(\tau)] = \left(\frac{\kappa}{{\ell}}\right)^{{{\varrho}}}\left(\mathbf{N}[\psi(\tau)]-\left(\frac{1}{\kappa}\right)\psi(0)\right). \end{equation}$

(2.10)

Definition 2.10. The NT of $D_{\tau}^{{{\varrho}}}\psi(\tau)$ in CF sense is stated as ^[47]

$\begin{equation} \mathbf{N}[D_{\tau}^{{{\varrho}}}\psi(\tau)] = \frac{1}{1-{{\varrho}}+{{\varrho}}(\frac{{\ell}}{\kappa})}\left(\mathbf{N}[\psi(\tau)]-\left(\frac{1}{\kappa}\right)\psi(0)\right). \end{equation}$

(2.11)

Definition 2.11. The NT of $D_{\tau}^{{{\varrho}}}\psi(\tau)$ in ABC sense is stated as ^[47]

$\begin{equation} \begin{split} \mathbf{N}[D_{\tau}^{{{\varrho}}}\psi(\tau)] = \frac{M[{{\varrho}}]}{1-{{\varrho}}+{{\varrho}}(\frac{{\ell}}{\kappa})^{{{\varrho}}}}\left(\mathbf{N}[\psi(\tau)]-\left(\frac{1}{\kappa}\right)\psi(0)\right), \end{split} \end{equation}$

(2.12)

with $M[{{\varrho}}]$ denoting a normalization function.

3. General implementation of method

Consider the fractional partial differential equation

$\begin{equation} \begin{split} D_{\tau}^{{{\varrho}}}\psi(\upsilon, \tau) = \mathcal{L}(\psi(\upsilon, \tau))+{N}(\psi(\upsilon, \tau))+h(\upsilon, \tau) = M(\upsilon, \tau), \end{split} \end{equation}$

(3.1)

with the initial condition

$\begin{equation} \begin{split} \psi(\upsilon, 0) = \phi(\upsilon), \end{split} \end{equation}$

(3.2)

which $\mathcal{L}$ is linear, ${N}$ non-linear functions and sources function $h(\upsilon, \tau)$ .

3.1. Case I $(NTDM_{CF})$

The fractional CF derivative and natural transformation, Eq (3.1) can be stated as,

$\begin{equation} \begin{split} \frac{1}{p({{\varrho}}, {{\ell}}, {\kappa})}\left(\mathbf{N}[\psi(\upsilon, \tau)]-\frac{\phi(\upsilon)}{\kappa}\right) = \mathbf{N}[M(\upsilon, \tau)], \end{split} \end{equation}$

(3.3)

with

$\begin{equation} \begin{split} p({{\varrho}}, {{\ell}}, {\kappa}) = 1-{{\varrho}}+{{\varrho}}(\frac{{\ell}}{\kappa}). \end{split} \end{equation}$

(3.4)

On employing natural inverse transformation, we get

$\begin{equation} \begin{split} \psi(\upsilon, \tau) = \mathbf{N}^{-1}\left(\frac{\phi(\upsilon)}{\kappa}+p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[M(\upsilon, \tau)]\right). \end{split} \end{equation}$

(3.5)

Thus, for $\psi(\upsilon, \tau)$ , the series type solution is given as

$\begin{equation} \begin{split} \psi(\upsilon, \tau) = \sum\limits_{i = 0}^{\infty}\psi_{i}(\upsilon, \tau), \end{split} \end{equation}$

(3.6)

and ${N}(\psi(\upsilon, \tau))$ can be decomposed as

$\begin{equation} \begin{split} {N}(\psi(\upsilon, \tau)) = \sum\limits_{i = 0}^{\infty}A_{i}(\psi_{0}, ..., \psi_{i}), \end{split} \end{equation}$

(3.7)

the Adomian polynomials $A_{i}$ is given as

$A_{n} = \frac{1}{n!} \frac{d^n}{d\varepsilon^n}{N} (t, \Sigma_{k = 0}^n \varepsilon^k \psi_k)|_{\varepsilon = 0.}$

Putting Eqs (3.7) and (3.6) into (3.5), we obtain

$\begin{equation} \begin{split} \sum\limits_{i = 0}^{\infty}\psi_{i}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(\frac{\phi(\upsilon)}{\kappa}+p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[h(\upsilon, \tau)]\right)\\ &+\mathbf{N}^{-1}\left(p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\sum\limits_{i = 0}^{\infty}\mathcal{L}(\psi_{i}(\upsilon, \tau))+A_{\tau}\right]\right). \end{split} \end{equation}$

(3.8)

From (3.8), we get,

$\begin{equation} \begin{split} \psi_{0}^{CF}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(\frac{\phi(\upsilon)}{\kappa}+p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[h(\upsilon, \tau)]\right), \\ \psi_{1}^{CF}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\mathcal{L}(\psi_{0}(\upsilon, \tau))+A_{0}\right]\right), \\ &\vdots\\ \psi_{l+1}^{CF}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\mathcal{L}(\psi_{l}(\upsilon, \tau))+A_{l}\right]\right), \ \ l = 1, 2, 3, \cdots.\\ \end{split} \end{equation}$

(3.9)

Thus, we obtain the result of (3.1) by putting (3.9) into (3.6) applying $NTDM_{CF},$

$\begin{equation} \begin{split} \psi^{CF}(\upsilon, \tau) = \psi_{0}^{CF}(\upsilon, \tau)+\psi_{1}^{CF}(\upsilon, \tau)+\psi_{2}^{CF}(\upsilon, \tau)+\cdots. \end{split} \end{equation}$

(3.10)

3.2. Case II $(NTDM_{ABC})$

The fractional ABC derivative and natural transformation, Eq (3.1) is given as,

$\begin{equation} \begin{split} \frac{1}{q({{\varrho}}, {{\ell}}, {\kappa})}\left(\mathbf{N}[\psi(\upsilon, \tau)]-\frac{\phi(\upsilon)}{\kappa}\right) = \mathbf{N}[M(\upsilon, \tau)], \end{split} \end{equation}$

(3.11)

with

$\begin{equation} \begin{split} q({{\varrho}}, {{\ell}}, {\kappa}) = \frac{1-{{\varrho}}+{{\varrho}}(\frac{{\ell}}{\kappa})^{{{\varrho}}}}{B({{\varrho}})}. \end{split} \end{equation}$

(3.12)

On employing natural inverse transformation, we get

$\begin{equation} \begin{split} \psi(\upsilon, \tau) = \mathbf{N}^{-1}\left(\frac{\phi(\upsilon)}{\kappa}+q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[M(\upsilon, \tau)]\right). \end{split} \end{equation}$

(3.13)

The Adomain decomposition, we obtain as

$\begin{equation} \begin{split} \sum\limits_{i = 0}^{\infty}\psi_{i}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(\frac{\phi(\upsilon)}{\kappa}+q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[h(\upsilon, \tau)]\right)\\ &+\mathbf{N}^{-1}\left(q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\sum\limits_{i = 0}^{\infty}\mathcal{L}(\psi_{i}(\upsilon, \tau))+A_{\tau}\right]\right). \end{split} \end{equation}$

(3.14)

From (3.8), we get,

$\begin{equation} \begin{split} \psi_{0}^{ABC}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(\frac{\phi(\upsilon)}{\kappa}+q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[h(\upsilon, \tau)]\right), \\ \psi_{1}^{ABC}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\mathcal{L}(\psi_{0}(\upsilon, \tau))+A_{0}\right]\right), \\ &\vdots\\ \psi_{l+1}^{ABC}(\upsilon, \tau) = &\mathbf{N}^{-1}\left(q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\mathcal{L}(\psi_{l}(\upsilon, \tau))+A_{l}\right]\right), \ \ l = 1, 2, 3, \cdots.\\ \end{split} \end{equation}$

(3.15)

Thus, we obtain the result of (3.1), by applying $NTDM_{ABC}$

$\begin{equation} \begin{split} \psi^{ABC}(\upsilon, \tau) = \psi_{0}^{ABC}(\upsilon, \tau)+\psi_{1}^{ABC}(\upsilon, \tau)+\psi_{2}^{ABC}(\upsilon, \tau)+\cdots. \end{split} \end{equation}$

(3.16)

4. Convergence analysis

In this part, we give the uniqueness and convergence of the $NTDM_{ABC}$ and $NTDM_{CF}$ .

Theorem 4.1. Let $|\mathcal{L}(\psi)-\mathcal{L}(\psi^{*})| < \gamma_{1}|\psi-\psi^{*}|$ and $|{N}(\psi)-{N}(\psi^{*})| < \gamma_{2}|\psi-\psi^{*}|$ , where $\psi: = \psi(\wp, \tau)$ and $\psi^{*}: = \psi^{*}(\wp, \tau)$ are values of two different functions and $\gamma_{1}$ , $\gamma_{2}$ are Lipschitz constants.

$\mathcal{L}$ and ${N}$ are the operators stated in (3.1). Then for $NTDM_{CF}$ the solution of (3.1) is unique when $0 < (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}} \tau) < 1$ for all $\tau.$

Proof. Let $H = (C[J], ||.||)$ having norm $||\phi(\tau)|| = {max}_{\tau\in J}|\phi(\tau)|$ is the Banach space, $\forall$ continuous function over the interval $J = [0, T].$ Let $I:H\rightarrow H$ a non-linear mapping, with

$\psi^{C}_{l+1} = \psi^{C}_{0}+\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[\mathcal{L}(\psi_{l}(\wp, \tau))+{N}(\psi_{l}(\wp, \tau))]], \ \ l\geq0.$

$\begin{equation} \begin{split} ||I(\psi)-I(\psi^{*})||&\leq max_{\tau\in J}|\mathbf{N}^{-1}\Big[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[\mathcal{L}(\psi)-\mathcal{L}(\psi^{*})]\\ &+p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[{N}(\psi)-{N}(\psi^{*})]|\Big]\\ &\leq max_{\tau\in J}\Big[\gamma_{1}\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[|\psi-\psi^{*}|]]\\ &+\gamma_{2}\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[|\psi-\psi^{*}|]]\Big]\\ &\leq max_{\tau\in J}(\gamma_{1}+\gamma_{2})\left[\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}|\psi-\psi^{*}|]\right]\\ &\leq(\gamma_{1}+\gamma_{2})\left[\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}||\psi-\psi^{*}||]\right]\\ & = (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}} \tau)||\psi-\psi^{*}||. \end{split} \end{equation}$

(4.1)

So, I is contraction as $0 < (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}} \tau) < 1$ . The solution of (3.1) is unique as according to Banach fixed point theorem.

Theorem 4.2. As according to the above theorem, the solution of (3.1) is unique for $NTDM_{ABC}$ when $0 < (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}}\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}) < 1$ for all $\tau.$

Proof. As from above theorem let $H = (C[J], ||.||)$ be the Banach space, $\forall$ continuous function over the interval $J.$ Let $I:H\rightarrow H$ is non-linear mapping, with

$\psi^{C}_{l+1} = \psi^{C}_{0}+\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[\mathcal{L}(\psi_{l}(\wp, \tau))+{N}(\psi_{l}(\wp, \tau))]], \ \ l\geq0.$

$\begin{equation} \begin{split} ||I(\psi)-I(\psi^{*})||&\leq max_{\tau\in J}|\mathbf{N}^{-1}\Big[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[\mathcal{L}(\psi)-\mathcal{L}(\psi^{*})]\\ &+q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[{N}(\psi)-{N}(\psi^{*})]|\Big]\\ &\leq max_{\tau\in J}\Big[\gamma_{1}\mathbf{N}^{-1}[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[|\psi-\psi^{*}|]]\\ &+\gamma_{2}\mathbf{N}^{-1}[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[|\psi-\psi^{*}|]]\Big]\\ &\leq max_{\tau\in J}(\gamma_{1}+\gamma_{2})\left[\mathbf{N}^{-1}[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}|\psi-\psi^{*}|]\right]\\ &\leq(\gamma_{1}+\gamma_{2})\left[\mathbf{N}^{-1}[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}||\psi-\psi^{*}||]\right]\\ & = (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}}\frac{\tau^{{{\varrho}}}}{\Gamma{{({\varrho}}+1)}})||\psi-\psi^{*}||. \end{split} \end{equation}$

(4.2)

So, I is contraction as $0 < (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}}\frac{\tau^{{{\varrho}}}}{\Gamma{{({\varrho}}+1)}}) < 1$ . The solution of (3.1) is unique as according to Banach fixed point theorem.

Theorem 4.3. Let $\mathcal{L}$ and ${N}$ are a Lipschitz functions as stated in the above theorems, then the $NTDM_{CF}$ solution of (3.1) is convergent.

Proof. Suppose $H$ be the Banach space determined above and let $\psi_{m} = \sum_{r = 0}^{m}\psi_{r}(\wp, \tau)$ . To prove that $\psi_{m}$ is a Cauchy sequence in H. Let,

$\begin{equation} \begin{split} ||\psi_{m}-\psi_{n}||& = max_{\tau\in J}|\sum\limits_{r = n+1}^{m}\psi_{r}|, \ \ n = 1, 2, 3, \cdots\\ &\leq max_{\tau\in J} \left|\mathbf{N}^{-1}\left[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\sum\limits_{r = n+1}^{m}(\mathcal{L}(\psi_{r-1})+{N}(\psi_{r-1}))\right]\right]\right|\\ & = max_{\tau\in J} \left|\mathbf{N}^{-1}\left[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\sum\limits_{r = n+1}^{m-1}(\mathcal{L}(\psi_{r})+{N}(\psi_{r}))\right]\right]\right|\\ &\leq max_{\tau\in J}|\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[(\mathcal{L}(\psi_{m-1})-\mathcal{L}(\psi_{n-1})+{N}(\psi_{m-1})-{N}(\psi_{n-1}))]]|\\ &\leq\gamma_{1}max_{\tau\in J}|\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[(\mathcal{L}(\psi_{m-1})-\mathcal{L}(\psi_{n-1}))]]|\\ &+\gamma_{2}max_{\tau\in J}|\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[({N}(\psi_{m-1})-{N}(\psi_{n-1}))]]|\\ & = (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}} \tau)||\psi_{m-1}-\psi_{n-1}||. \end{split} \end{equation}$

(4.3)

Let $m = n+1$ , then

$\begin{equation} \begin{split} ||\psi_{n+1}-\psi_{n}||\leq \gamma||\psi_{n}-\psi_{n-1}||\leq \gamma^{2}||\psi_{n-1}\psi_{n-2}||\leq \cdots \leq \gamma^{n}||\psi_{1}-\psi_{0}||, \end{split} \end{equation}$

(4.4)

where $\gamma = (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}} \tau)$ . Thus, we get

$\begin{equation} \begin{split} ||\psi_{m}-\psi_{n}||&\leq||\psi_{n+1}-\psi_{n}||+||\psi_{n+2}\psi_{n+1}||+\cdots+||\psi_{m}-\psi_{m-1}||, \\ &(\gamma^{n}+\gamma^{n+1}+\cdots+\gamma^{m-1})||\psi_{1}-\psi_{0}||\\ &\leq \gamma^{n}\left(\frac{1-\gamma^{m-n}}{1-\gamma}\right)||\psi_{1}||. \end{split} \end{equation}$

(4.5)

As $0 < \gamma < 1$ , we have $1-\gamma^{m-n} < 1$ . Thus,

$\begin{equation} \begin{split} ||\psi_{m}-\psi_{n}||\leq \frac{\gamma^{n}}{1-\gamma}max_{\tau\in J}||\psi_{1}||. \end{split} \end{equation}$

(4.6)

Since $||\psi_{1}|| < \infty, \ \ ||\psi_{m}-\psi_{n}||\rightarrow 0$ when $n \rightarrow \infty$ . As a result, $\psi_m$ is a Cauchy sequence in H, stated that the series $\psi_m$ is convergent.

Theorem 4.4. Let $\mathcal{L}$ and ${N}$ are a Lipschitz functions as stated in the above theorems, then the $NTDM_{ABC}$ solution of (3.1) is convergent.

Proof. Suppose $\psi_{m} = \sum_{r = 0}^{m}\psi_{r}(\wp, \tau)$ . To prove that $\psi_{m}$ is a Cauchy sequence in H. Let,

$\begin{align} \begin{split} ||\psi_{m}-\psi_{n}||& = max_{\tau\in J}|\sum\limits_{r = n+1}^{m}\psi_{r}|, \ \ n = 1, 2, 3, \cdots\\ &\leq max_{\tau\in J} \left|\mathbf{N}^{-1}\left[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\sum\limits_{r = n+1}^{m}(\mathcal{L}(\psi_{r-1})+{N}(\psi_{r-1}))\right]\right]\right|\\ & = max_{\tau\in J} \left|\mathbf{N}^{-1}\left[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}\left[\sum\limits_{r = n+1}^{m-1}(\mathcal{L}(\psi_{r})+{N}(u_{r}))\right]\right]\right|\\ &\leq max_{\tau\in J}|\mathbf{N}^{-1}[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[(\mathcal{L}(\psi_{m-1})-\mathcal{L}(\psi_{n-1})+{N}(\psi_{m-1})-{N}(\psi_{n-1}))]]|\\ &\leq\gamma_{1}max_{\tau\in J}|\mathbf{N}^{-1}[q({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[(\mathcal{L}(\psi_{m-1})-\mathcal{L}(\psi_{n-1}))]]|\\ &+\gamma_{2}max_{\tau\in J}|\mathbf{N}^{-1}[p({{\varrho}}, {{\ell}}, {\kappa})\mathbf{N}[({N}(\psi_{m-1})-{N}(\psi_{n-1}))]]|\\\nonumber & = (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}}\frac{\tau^{{{\varrho}}}}{\Gamma({{\varrho}}+1)})||\psi_{m-1}-\psi_{n-1}||. \end{split} \end{align}$

(4.7)

Suppose $m = n+1$ , thus

$\begin{equation} \begin{split} ||\psi_{n+1}-\psi_{n}||\leq \gamma||\psi_{n}-\psi_{n-1}||\leq \gamma^{2}||\psi_{n-1}-\psi_{n-2}||\leq \cdots \leq \gamma^{n}||\psi_{1}-\psi_{0}||, \end{split} \end{equation}$

(4.8)

with $\gamma = (\gamma_{1}+\gamma_{2})(1-{{\varrho}}+{{\varrho}}\frac{\tau^{{{\varrho}}}}{\Gamma({{\varrho}}+1)})$ . Thus, we get

$\begin{equation} \begin{split} ||\psi_{m}-\psi_{n}||&\leq||\psi_{n+1}-\psi_{n}||+||\psi_{n+2}-\psi_{n+1}||+\cdots+||\psi_{m}-\psi_{m-1}||, \\ &\leq(\gamma^{n}+\gamma^{n+1}+\cdots+\gamma^{m-1})||\psi_{1}-\psi_{0}||\\ &\leq \gamma^{n}\left(\frac{1-\gamma^{m-n}}{1-\gamma}\right)||\psi_{1}||. \end{split} \end{equation}$

(4.9)

As $0 < \gamma < 1$ , we have $1-\gamma^{m-n} < 1$ . Thus,

$\begin{equation} \begin{split} ||\psi_{m}-\psi_{n}||\leq \frac{\gamma^{n}}{1-\gamma}max_{\tau\in J}||\psi_{1}||. \end{split} \end{equation}$

(4.10)

Since $||\psi_{1}|| < \infty, \ \ ||\psi_{m}-\psi_{n}||\rightarrow 0$ when $n \rightarrow \infty$ . As a result, $\psi_m$ is a Cauchy sequence in H, stated that the series $\psi_m$ is convergent.

5. Applications

Example 1. Consider the fractional-order system of Whitham-Broer-Kaup equations is given as

$\begin{array}{l} D_{{\tau}}^{\varrho}\psi(\upsilon, \tau)+{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}} = 0, \ \ 0 < \varrho \leq 1, \quad -1 < {\tau} \leq 1, \quad -10 \leq {\upsilon} \leq 10, \\ D_{{\tau}}^{\varrho}\varphi(\upsilon, \tau)+{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}} = 0, \end{array}$

(5.1)

the initial condition

$\begin{equation} \begin{split} &\psi({\upsilon}, 0) = \frac{1}{2}-8\tanh(-2{\upsilon}), \\ &\varphi({\upsilon}, 0) = 16-16\tanh^2(-2{\upsilon}). \end{split} \end{equation}$

(5.2)

On employing natural transform, we have

$\begin{equation} \begin{split} &{\mathbf{N}}[D_{\tau}^{\varrho}\psi({\upsilon}, {\tau})] = -{\mathbf{N}}\left[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\right], \\ &{\mathbf{N}}[D_{\tau}^{\varrho}\varphi({\upsilon}, {\tau})] = -{\mathbf{N}}\left[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\right], \end{split} \end{equation}$

(5.3)

Thus, we have

$\begin{equation} \begin{split} \frac{1}{\kappa^{\varrho}}{\mathbf{N}}[\psi({\upsilon}, {\tau})]-{\kappa}^{2-\varrho}{{\psi}({\upsilon}, 0)}& = -{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg], \\ \frac{1}{\kappa^{\varrho}}{\mathbf{N}}[\varphi({\upsilon}, {\tau})]-{\kappa}^{2-\varrho}{{\psi}({\upsilon}, 0)}& = -{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\\ &+3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg].\\ \end{split} \end{equation}$

(5.4)

On simplification we have

$\begin{array}{l} {\mathbf{N}}[\psi({\upsilon}, {\tau})] = {{\kappa}^{2}}\Bigg[\frac{1}{2}-8\tanh(-2{\upsilon})\Bigg]-\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\\ +\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg], \\ {\mathbf{N}}[\varphi({\upsilon}, {\tau})] = {{\kappa}^{2}}\Bigg[16-16\tanh^2(-2{\upsilon})\Bigg]-\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\\ +3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg].\end{array}$

(5.5)

Using the inverse natural transform, we get

$\begin{equation} \begin{split} \psi({\upsilon}, {\tau})& = \Bigg[\frac{1}{2}-8\tanh(-2{\upsilon})\Bigg]\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg\}\Bigg], \\ \varphi({\upsilon}, {\tau})& = \Bigg[16-16\tanh^2(-2{\upsilon})\Bigg]\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\\ &+3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg\}\Bigg].\\ \end{split} \end{equation}$

(5.6)

5.1. $NTDM_{CF}$ solution

The series form solution for the unknown term $\psi(\upsilon, \tau)$ and $\varphi(\upsilon, \tau)$ is given as

$\begin{equation} \psi({\upsilon}, {\tau}) = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}) \quad and \quad \varphi({\upsilon}, {\tau}) = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}) . \end{equation}$

(5.7)

The non-linear functions define by Adomian polynomials are given as $\psi\psi_{\upsilon} = \sum_{m = 0}^{\infty}\mathcal{A}_m$ , $\psi\varphi_{\upsilon} = \sum_{m = 0}^{\infty}\mathcal{B}_m$ and $\varphi\psi_{\upsilon} = \sum_{m = 0}^{\infty}\mathcal{C}_m$ , thus by means of these functions Eq (5.6) can be calculated as

$\begin{equation} \begin{split} \sum\limits_{{l} = 0}^{\infty}\psi_{{l}+1}({\upsilon}, {\tau})& = \frac{1}{2}-8\tanh(-2{\upsilon})\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal A_{{l}}\\ &+\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg\}\Bigg], \\ \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}+1}({\upsilon}, {\tau})& = 16-16\tanh^2(-2{\upsilon})\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal B_{{l}}+\sum\limits_{{l} = 0}^{\infty}\mathcal C_{{l}}\\ &+3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg\}\Bigg].\\ \end{split} \end{equation}$

(5.8)

Both side comparison of Eq (5.8), we achieve

$\begin{equation*} \begin{split} &\psi_{0}({\upsilon}, {\tau}) = \frac{1}{2}-8\tanh(-2{\upsilon}), \\ &\varphi_{0}({\upsilon}, {\tau}) = 16-16\tanh^2(-2{\upsilon}), \\ \end{split} \end{equation*}$

$\begin{equation} \begin{split} &\psi_{1}({\upsilon}, {\tau}) = -8\sec h^2(-2{\upsilon})\left(\varrho(\tau-1)+1\right), \\ &\varphi_{1}({\upsilon}, {\tau}) = -32{\rm{sech}}^2(-2{\upsilon})\tanh(-2{\upsilon})\left(\varrho(\tau-1)+1\right), \\ \end{split} \end{equation}$

(5.9)

$\begin{equation} \begin{split} \psi_{2}({\upsilon}, {\tau})& = -16{\rm{sech}}^2(-2{\upsilon})\Big(4{\rm{sech}}^2(-2{\upsilon})\\ &-8\tanh^2(-2{\upsilon})+3\tanh(-2{\upsilon})\Big)\Big((1-\varrho)^{2}\\ &+2\varrho(1-\varrho)\tau+\frac{\varrho^{2}\tau^2}{2}\Big), \\ \varphi_{2}({\upsilon}, {\tau})& = -32\sec h^2(-2{\upsilon})\{40\sec h^2(-2{\upsilon})\tanh(-2{\upsilon})\\ &+96\tanh(-2{\upsilon})-2\tanh^2(-2{\upsilon}) -32\tanh^3(-2{\upsilon})\\ &-25\sec h^2(-2{\upsilon})\}\Big((1-\varrho)^{2}+2\varrho(1-\varrho)\tau+\frac{\varrho^{2}\tau^2}{2}\Big).\\ \end{split} \end{equation}$

(5.10)

Thus, for $\psi_{{l}}$ and $\varphi_{{l}}$ with $({l}\geq3)$ the remaining terms are simply calculate. So, the result in series type is as

(5.11)

5.2. $NDM_{ABC}$ solution

The series form solution for the unknown term $\psi(\upsilon, \tau)$ and $\varphi(\upsilon, \tau)$ is given as

$\begin{equation} \begin{split} &\psi({\upsilon}, {\tau}) = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}), \\ &\varphi({\upsilon}, {\tau}) = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}).\\ \end{split} \end{equation}$

(5.12)

The non-linear functions by mean of Adomian polynomials is given as $\psi\psi_{{\upsilon}} = \sum_{{l} = 0}^{\infty}\mathcal A_{{l}}$ and $\psi^{2}\psi_{{\upsilon}} = \sum_{{l} = 0}^{\infty}\mathcal B_{{l}}$ , thus by mean of these functions Eq (5.6) can be calculated as

$\begin{equation} \begin{split} \sum\limits_{{l} = 0}^{\infty}\psi_{{l}+1}({\upsilon}, {\tau})& = \frac{1}{2}-8\tanh(-2{\upsilon})\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{\ell^\varrho(\kappa^\varrho+\varrho(\ell^\varrho-{\kappa}^\varrho))}{{\kappa}^{2\varrho}}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal A_{{l}}+\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg\}\Bigg], \\ \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}+1}({\upsilon}, {\tau})& = 16-16\tanh^2(-2{\upsilon})\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{\ell^\varrho(\kappa^\varrho+\varrho(\ell^\varrho-{\kappa}^\varrho))}{{\kappa}^{2\varrho}}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal B_{{l}}+\sum\limits_{{l} = 0}^{\infty}\mathcal C_{{l}}+3\frac{\partial^3{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^3}}-\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg\}\Bigg].\\ \end{split} \end{equation}$

(5.13)

Both side comparison of Eq (5.13), we get

$\begin{equation*} \begin{split} &\psi_{0}({\upsilon}, {\tau}) = \frac{1}{2}-8\tanh(-2{\upsilon}), \\ &\varphi_{0}({\upsilon}, {\tau}) = 16-16\tanh^2(-2{\upsilon}), \\ \end{split} \end{equation*}$

$\begin{equation} \begin{split} &\psi_{1}({\upsilon}, {\tau}) = -8{\rm{sech}}^2(-2{\upsilon})\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right), \\ &\varphi_{1}({\upsilon}, {\tau}) = -32{\rm{sech}}^2(-2{\upsilon})\tanh(-2{\upsilon})\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right), \\ \end{split} \end{equation}$

(5.14)

$\begin{equation} \begin{split} \psi_{2}({\upsilon}, {\tau})& = -16\sec h^2(-2{\upsilon})\left(4\sec h^2(-2{\upsilon})-8\tanh^2(-2{\upsilon})+3\tanh(-2{\upsilon}) \right) \Bigg[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}\\ &+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\Bigg], \\ \varphi_{2}({\upsilon}, {\tau})& = -32\sec h^2(-2{\upsilon})\{40\sec h^2(-2{\upsilon})\tanh(-2{\upsilon})+96\tanh(-2{\upsilon})-2\tanh^2(-2{\upsilon})\\ &-32\tanh^3(-2{\upsilon})-25\sec h^2(-2{\upsilon})\}\left[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\right].\\ \end{split} \end{equation}$

(5.15)

Thus, for $\psi_{{l}}$ with $({l}\geq3)$ the remaining component are simply calculated. So, the series form solution is given as

$\begin{equation} \begin{split} \psi({\upsilon}, {\tau})& = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}) = \psi_{0}({\upsilon}, {\tau})+\psi_{1}({\upsilon}, {\tau})+\psi_{2}({\upsilon}, {\tau})+\cdots, \\ \psi({\upsilon}, {\tau})& = \frac{1}{2}-8\tanh(-2{\upsilon})-8{\rm{sech}}^2(-2{\upsilon})\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right)\\ &-16\sec h^2(-2{\upsilon})\Bigg(4\sec h^2(-2{\upsilon})-8\tanh^2(-2{\upsilon})\\ &+3\tanh(-2{\upsilon}) \Bigg)\left[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\right]+\cdots.\\ \varphi({\upsilon}, {\tau})& = \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}}({\upsilon}, {\tau}) = \varphi_{0}({\upsilon}, {\tau})+\varphi_{1}({\upsilon}, {\tau})+\varphi_{2}({\upsilon}, {\tau})+\cdots, \\ \varphi({\upsilon}, {\tau})& = 16-16\tanh^2(-2{\upsilon})-32{\rm{sech}}^2(-2{\upsilon})\tanh(-2{\upsilon})\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right)\\ &-32\sec h^2(-2{\upsilon})\{40\sec h^2(-2{\upsilon})\tanh(-2{\upsilon})+96\tanh(-2{\upsilon})-2\tanh^2(-2{\upsilon})\\ &-32\tanh^3(-2{\upsilon})-25\sec h^2(-2{\upsilon})\}\Bigg[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\Bigg]+\cdots.\\ \end{split} \end{equation}$

(5.16)

The exact result is

$\begin{equation} \begin{split} &\psi({\upsilon}, {\tau}) = \frac{1}{2}-8\tanh\left\{-2\left({\upsilon}-\frac{{\tau}}{2}\right)\right\}, \\ &\varphi({\upsilon}, {\tau}) = 16-16\tanh^2\left\{-2\left({\upsilon}-\frac{{\tau}}{2}\right)\right\}.\\ \end{split} \end{equation}$

(5.17)

Example 2. Consider the fractional-order system of Whitham-Broer-Kaup equations is given as

$\begin{equation} \begin{split} &D_{{\tau}}^{\varrho}\psi(\upsilon, \tau)+{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}} = 0, \\ &D_{{\tau}}^{\varrho}\varphi(\upsilon, \tau)+{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}} = 0, \\ &0 < \varrho \leq 1, \quad -1 < {\tau} \leq 1, \quad -10 \leq {\upsilon} \leq 10, \end{split} \end{equation}$

(5.18)

with the initial condition

$\begin{equation} \begin{split} &\psi({\upsilon}, 0) = \lambda-\kappa\coth[\kappa({\upsilon}+\theta)], \\ &\varphi({\upsilon}, 0) = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]. \end{split} \end{equation}$

(5.19)

On employing natural transform, we have

$\begin{equation} \begin{split} &{\mathbf{N}}[D_{\tau}^{\varrho}\psi({\upsilon}, {\tau})] = -{\mathbf{N}}\left[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\right], \\ &{\mathbf{N}}[D_{\tau}^{\varrho}\varphi({\upsilon}, {\tau})] = -{\mathbf{N}}\left[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\right]. \end{split} \end{equation}$

(5.20)

Thus, we have

$\begin{equation} \begin{split} &\frac{1}{\kappa^{\varrho}}{\mathbf{N}}[\psi({\upsilon}, {\tau})]-{\kappa}^{2-\varrho}{{\psi}({\upsilon}, 0)} = -{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg], \\ &\frac{1}{\kappa^{\varrho}}{\mathbf{N}}[\varphi({\upsilon}, {\tau})]-{\kappa}^{2-\varrho}{{\psi}({\upsilon}, 0)} = -{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg].\\ \end{split} \end{equation}$

(5.21)

On simplification we have

$\begin{equation} \begin{split} {\mathbf{N}}[\psi({\upsilon}, {\tau})]& = {{\kappa}^{2}}\Bigg[\lambda-\kappa\coth[\kappa({\upsilon}+\theta)]\Bigg]-\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\\ &+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg], \\ {\mathbf{N}}[\varphi({\upsilon}, {\tau})]& = {{\kappa}^{2}}\Bigg[-\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\Bigg]-\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}{\mathbf{N}}\Bigg[{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\\ &+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg].\\ \end{split} \end{equation}$

(5.22)

Using the inverse natural transform, we have

$\begin{equation} \begin{split} \psi({\upsilon}, {\tau})& = \Bigg[\lambda-\kappa\coth[\kappa({\upsilon}+\theta)]\Bigg]\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg\}\Bigg], \\ \varphi({\upsilon}, {\tau})& = \Bigg[-\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\Bigg]\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{{{\psi({\upsilon}, {\tau})}}\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+{{\varphi({\upsilon}, {\tau})}}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg\}\Bigg].\\ \end{split} \end{equation}$

(5.23)

5.3. $NDM_{CF}$ solution

The series form solution for the unknown term $\psi(\upsilon, \tau)$ and $\varphi(\upsilon, \tau)$ is given as

(5.24)

The non-linear functions by Adomian polynomials is given as $\psi\psi_{\upsilon} = \sum_{m = 0}^{\infty}\mathcal{A}_m$ , $\psi\varphi_{\upsilon} = \sum_{m = 0}^{\infty}\mathcal{B}_m$ and $\varphi\psi_{\upsilon} = \sum_{m = 0}^{\infty}\mathcal{C}_m$ , thus by mean of these functions Eq (5.23) can be calculated as

$\begin{equation} \begin{split} \sum\limits_{{l} = 0}^{\infty}\psi_{{l}+1}({\upsilon}, {\tau})& = \lambda-\kappa\coth[\kappa({\upsilon}+\theta)]\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal A_{{l}}+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg\}\Bigg], \\ \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}+1}({\upsilon}, {\tau})& = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\\ &-{\mathbf{N}}^{-1}\Bigg[\frac{{{\varrho}}(\kappa-\varrho({\kappa}-{{\varrho}}))}{{\kappa}^2}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal B_{{l}}+\sum\limits_{{l} = 0}^{\infty}\mathcal C_{{l}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg\}\Bigg].\\ \end{split} \end{equation}$

(5.25)

Both sides comparisons of Eq (5.25), we achieve

$\begin{equation*} \begin{split} &\psi_{0}({\upsilon}, {\tau}) = \lambda-\kappa\coth[\kappa({\upsilon}+\theta)], \\ &\varphi_{0}({\upsilon}, {\tau}) = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)], \\ \end{split} \end{equation*}$

$\begin{equation} \begin{split} &\psi_{1}({\upsilon}, {\tau}) = -\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\left(\varrho(\tau-1)+1\right), \\ &\varphi_{1}({\upsilon}, {\tau}) = -\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\coth[\kappa({\upsilon}+\theta)]\left(\varrho(\tau-1)+1\right), \end{split} \end{equation}$

(5.26)

$\begin{equation} \begin{split} \psi_{2}({\upsilon}, {\tau})& = \lambda\kappa^4 {csch}^2[\kappa({\upsilon}+\theta)]\Big\{{2\lambda\kappa}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+(1-{\varrho})^{3}+\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}\\ &-{(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Big((1-\varrho)^{2}+2\varrho(1-\varrho)\tau+\frac{\varrho^{2}\tau^2}{2}\Big)\Big\}, \\ \varphi_{2}({\upsilon}, {\tau})& = [2\lambda\kappa^5 {csch}^2[\kappa({\upsilon}+\theta)]]\left[{\lambda\kappa {csch}^2(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+\right.\\ &\left.+(1-{\varrho})^{3} \frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}+\frac{2\lambda\kappa {csch}^2\coth^2([\kappa({\upsilon}+\theta)]){\tau}^{3\varrho}}{\Gamma(\varrho+1)\Gamma(3\varrho+1)}\right.\\ &\left.-{2\lambda\coth(3 {csch}^2([\kappa({\upsilon}+\theta)]-1))}\Big((1-\varrho)^{2}+2\varrho(1-\varrho)\tau+\frac{\varrho^{2}\tau^2}{2}\Big)\right]. \end{split} \end{equation}$

(5.27)

Thus, for $\psi_{{l}}$ and $\varphi_{{l}}$ with $({l}\geq3)$ the remaining terms are simply calculate. So, the series form solution is as

$\begin{equation} \begin{split} \psi({\upsilon}, {\tau})& = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}) = \psi_{0}({\upsilon}, {\tau})+\psi_{1}({\upsilon}, {\tau})+\psi_{2}({\upsilon}, {\tau})+\cdots, \\ \psi({\upsilon}, {\tau})& = \lambda-\kappa\coth[\kappa({\upsilon}+\theta)]-\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\left(\varrho(\tau-1)+1\right)\\ &+\lambda\kappa^4 {csch}^2[\kappa({\upsilon}+\theta)]\Big\{{2\lambda\kappa}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+(1-{\varrho})^{3}+\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}\\ &-{(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Big((1-\varrho)^{2}+2\varrho(1-\varrho)\tau+\frac{\varrho^{2}\tau^2}{2}\Big)\Big\}+\cdots.\\ \varphi({\upsilon}, {\tau})& = \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}}({\upsilon}, {\tau}) = \varphi_{0}({\upsilon}, {\tau})+\varphi_{1}({\upsilon}, {\tau})+\varphi_{2}({\upsilon}, {\tau})+\cdots, \\ \varphi({\upsilon}, {\tau})& = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]-\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\coth[\kappa({\upsilon}+\theta)]\left(\varrho(\tau-1)+1\right)\\ &+[2\lambda\kappa^5 {csch}^2[\kappa({\upsilon}+\theta)]]\Bigg[{\lambda\kappa {csch}^2(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}\\ &+(1-{\varrho})^{3}+\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}+\frac{2\lambda\kappa {csch}^2\coth^2([\kappa({\upsilon}+\theta)]){\tau}^{3\varrho}}{\Gamma(\varrho+1)\Gamma(3\varrho+1)}\\ &-{2\lambda \coth(3 {csch}^2([\kappa({\upsilon}+\theta)]-1))}\Big((1-\varrho)^{2}+2\varrho(1-\varrho)\tau+\frac{\varrho^{2}\tau^2}{2}\Big)\Bigg]+\cdots. \end{split} \end{equation}$

(5.28)

5.4. $NDM_{ABC}$ solution

The series form solution for the unknown term $\psi(\upsilon, \tau)$ and $\varphi(\upsilon, \tau)$ is given as

(5.29)

The non-linear function by mean of Adomian polynomials are define as $\psi\psi_{{\upsilon}} = \sum_{{l} = 0}^{\infty}\mathcal A_{{l}}$ and $\psi^{2}\psi_{{\upsilon}} = \sum_{{l} = 0}^{\infty}\mathcal B_{{l}}$ , thus by mean of these function Eq (5.23) can be calculated as

$\begin{equation} \begin{split} \sum\limits_{{l} = 0}^{\infty}\psi_{{l}+1}({\upsilon}, {\tau})& = \lambda-\kappa\coth[\kappa({\upsilon}+\theta)]\\ &+{\mathbf{N}}^{-1}\Bigg[\frac{\ell^\varrho(\kappa^\varrho+\varrho(\ell^\varrho-{\kappa}^\varrho))}{{\kappa}^{2\varrho}}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal A_{{l}}+\frac{1}{2}\frac{\partial{{\psi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}+\frac{\partial{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}}}\Bigg\}\Bigg], \\ \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}+1}({\upsilon}, {\tau})& = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\\ &+{\mathbf{N}}^{-1}\Bigg[\frac{\ell^\varrho(\kappa^\varrho+\varrho(\ell^\varrho-{\kappa}^\varrho))}{{\kappa}^{2\varrho}}\mathbf{N}\Bigg\{\sum\limits_{{l} = 0}^{\infty}\mathcal B_{{l}}+\sum\limits_{{l} = 0}^{\infty}\mathcal C_{{l}}-\frac{1}{2}\frac{\partial^2{{\varphi({\upsilon}, {\tau})}}}{\partial{{\upsilon}^2}}\Bigg\}\Bigg].\\ \end{split} \end{equation}$

(5.30)

Both side comparison of Eq (5.30), we achieve

$\begin{equation} \begin{split} &\psi_{1}({\upsilon}, {\tau}) = -\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right), \\ &\varphi_{1}({\upsilon}, {\tau}) = -\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\coth[\kappa({\upsilon}+\theta)]\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right), \\ \end{split} \end{equation}$

(5.31)

$\begin{equation} \begin{split} \psi_{2}({\upsilon}, {\tau})& = \lambda\kappa^4 {csch}^2[\kappa({\upsilon}+\theta)]\Big\{{2\lambda\kappa}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+(1-{\varrho})^{3}+\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}\\ &-{(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\left[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\right]\Big\}, \\ \varphi_{2}({\upsilon}, {\tau})& = [2\lambda\kappa^5 {csch}^2[\kappa({\upsilon}+\theta)]]\Bigg[{\lambda\kappa {csch}^2(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+(1-{\varrho})^{3}\\ &+\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}+\frac{2\lambda\kappa {csch}^2\coth^2([\kappa({\upsilon}+\theta)]){\tau}^{3\varrho}}{\Gamma(\varrho+1)\Gamma(3\varrho+1)}\\ &-{2\lambda \coth(3 {csch}^2([\kappa({\upsilon}+\theta)]-1))}\Bigg[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\Bigg]\Bigg].\\ \end{split} \end{equation}$

(5.32)

Thus, for $\psi_{{l}}$ and $\varphi_{{l}}$ with $({l}\geq3)$ , the remaining terms are simply calculated. So, the series form solution is as

$\begin{array}{l}\psi({\upsilon}, {\tau}) = \sum\limits_{{l} = 0}^{\infty}\psi_{{l}}({\upsilon}, {\tau}) = \psi_{0}({\upsilon}, {\tau})+\psi_{1}({\upsilon}, {\tau})+\psi_{2}({\upsilon}, {\tau})+\cdots, \\ \psi({\upsilon}, {\tau}) = \lambda-\kappa\coth[\kappa({\upsilon}+\theta)]-\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right)\\ +\lambda\kappa^4 {csch}^2[\kappa({\upsilon}+\theta)]\Big\{{2\lambda\kappa}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+(1-{\varrho})^{3}+\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}\\ -{(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Bigg[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\Bigg]\Big\}+\cdots.\\ \varphi({\upsilon}, {\tau}) = \sum\limits_{{l} = 0}^{\infty}\varphi_{{l}}({\upsilon}, {\tau}) = \varphi_{0}({\upsilon}, {\tau})+\varphi_{1}({\upsilon}, {\tau})+\varphi_{2}({\upsilon}, {\tau})+\cdots, \\ \varphi({\upsilon}, {\tau}) = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]-\lambda\kappa^2 {csch}^2[\kappa({\upsilon}+\theta)]\coth[\kappa({\upsilon}+\theta)]\left(1-\varrho+\frac{\varrho\tau^\varrho}{\Gamma(\varrho+1)}\right)\\ +[2\lambda\kappa^5 {csch}^2[\kappa({\upsilon}+\theta)]]\Bigg[{\lambda\kappa {csch}^2(3\coth^2([\kappa({\upsilon}+\theta)]-1))}\Big\{(1-{\varrho})^{2}3{\varrho}{{\tau}}+(1-{\varrho})^{3}\\ +\frac{3{\varrho}^{2}(1-{\varrho}){{\tau}}^{2}}{2}+\frac{{\varrho}^{3}{{\tau}}^{3}}{3!}\Big\}+\frac{2\lambda\kappa {csch}^2\coth^2([\kappa({\upsilon}+\theta)]){\tau}^{3\varrho}}{\Gamma(\varrho+1)\Gamma(3\varrho+1)}\\ -{2\lambda \coth(3 {csch}^2([\kappa({\upsilon}+\theta)]-1))}\Bigg[\frac{\varrho^{2}\tau^{2\varrho}}{\Gamma(2\varrho+1)}+2\varrho(1-\varrho)\frac{\tau^{\varrho}}{\Gamma(\varrho+1)}+(1-\varrho)^{2}\Bigg]\Bigg]+\cdots. \end{array}$

(5.33)

We achieve the series type result at integer-order $\varrho = 1, \kappa = 0.1$ , $\lambda = 0.005$ , $\theta = 10$ , as

$\begin{aligned} \psi(\upsilon, \tau) = & 0.005-0.1 \operatorname{coth}(0.1 \upsilon+10)-0.0005 {csch}^{2}(0.1 \upsilon+10) \tau\\ &+5 \times 10^{-7} {csch}^{2}(0.1 \upsilon+10) 0.003 \tau^{3} \\ &\left.-0.5\left(3 \operatorname{coth}^{2}(0.1 \upsilon+10)-1 .\right)\right) \tau^{2}, \\ \varphi(\upsilon, \tau) = &-0.01 {csch}^{2}(0.1 \upsilon+10)-0.000010 {csch}^{2}(0.1 \upsilon+10) \\ &\times \operatorname{coth}(0.1 \upsilon+10) \tau+1.0 \times 10^{-7} {csch}^{2}(0.1 \upsilon+10)\\ &\times \left[8.3 \times 10^{-5} \tau^{3} {csch}^{2}(0.1 \upsilon+10)(3 \operatorname{coth}(0.1 \upsilon+10)-1) \right.\\ &\left.-\tau^{2} \operatorname{coth}(0.1 \upsilon+10)\left(3 {csch}^{2}(0.1 \upsilon+10)-1\right) \right.\\ &\left.+1.6 \times 10^{-4} \tau^{3} {cosec} h^{2}(0.1 \upsilon+10) \operatorname{coth}(0.1 \upsilon+10)\right] . \end{aligned}$

The exact result of Eq (5.18) at $\varrho = 1$ , and taking $\kappa = 0.1$ , $\lambda = 0.005$ , $\theta = 10$ ,

$\begin{equation} \begin{split} &\psi({\upsilon}, {\tau}) = =\lambda-\kappa\coth[\kappa({\upsilon}+\theta-\lambda {\tau})], \\ &{\varphi({\upsilon}, {\tau})} = -\kappa^2 {csch}^2[\kappa({\upsilon}+\theta-\lambda {\tau})]. \end{split} \end{equation}$

(5.34)

6. Numerical results and discussion

In this study, we investigated the numerical solution of coupled fractional Whitham-Broer-Kaup equation systems using two novel approaches. For any order and different values of the space and time variables, numerical data for the coupled fractional Whitham-Broer-Kaup equations can be found using Maple. We create numerical simulations for system 1 at various $\upsilon$ and $\tau$ values in Tables 1 and 2. The Adomian decomposition approach, variational iteration method, optimal homotopy asymptotic method, and natural decomposition method are numerically compared in Tables 3 and 4 in terms of absolute error. and display the results of calculations performed for the coupled system taken into account in Problem 2. We can conclude that the Natural decomposition approach yields more accurate results based on the data in the tables above. For $\psi(\upsilon, \tau)$ of Problem 1, shows the behavior of the exact and Natural decomposition approach result, while shows the behavior of the analytical result at various fractional-orders of $\varrho$ . shows the solution to the equation $\psi(\upsilon, \tau)$ in various fractional orders, whereas shows the absolute error. The behavior of the exact and analytical solutions for $\varphi(\upsilon, \tau)$ is shown in , while the behavior of the analytical result at various fractional orders of $\varrho$ is shown in and . The behavior of the exact and analytical solution for $\psi(\upsilon, \tau)$ is depicted in , while the absolute error for Problem 2 is shown in . Similarly, the behavior of the exact and analytical solution for $\varphi(\upsilon, \tau)$ is depicted in Figure 10, while the absolute error is shown in Figure 11 for Problem 2.

Table 1. The suggested techniques result for

$\psi({\upsilon}, {\tau})$ at various fractional-orders of Example 1.

$({\upsilon}, {\tau})$	$\psi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.6	$\psi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.8	$(NTDM_{ABC})$ at ${ \varrho }=$ 1	$(NTDM_{CF})$ at ${ \varrho }=$ 1	Exact solution
(0.2, 0.01)	3.538802	3.538856	3.538907	3.538907	3.538907
(0.4, 0.02)	5.811722	5.811755	5.811846	5.811846	5.811846
(0.6, 0.03)	7.168821	7.168876	7.168992	7.168992	7.168992
(0.2, 0.01)	3.525702	3.525764	3.525891	3.525891	3.525891
(0.4, 0.02)	5.803176	5.803245	5.803337	5.803337	5.803337
(0.6, 0.03)	7.164207	7.164278	7.164348	7.164348	7.164348
(0.2, 0.01)	3.537413	3.537498	3.537537	3.537537	3.537537
(0.4, 0.02)	5.810812	5.810856	5.810952	5.810952	5.810952
(0.6, 0.03)	7.168423	7.168479	7.168504	7.168504	7.168504
(0.2, 0.01)	3.536701	3.536746	3.536853	3.536853	3.536853
(0.4, 0.02)	5.810411	5.810466	5.810504	5.810504	5.810504
(0.6, 0.03)	7.168112	7.168160	7.168260	7.168260	7.168260
(0.2, 0.01)	3.536076	3.536134	3.536168	3.536168	3.536168
(0.4, 0.02)	5.810002	5.810010	5.810057	5.810057	5.810057
(0.6, 0.03)	7.168000	7.168002	7.168016	7.168016	7.168016

| Show Table

DownLoad: CSV

Table 2. The suggested techniques solution for

$\varphi({\upsilon}, {\tau})$ at various fractional-orders of Example 1.

$({\upsilon}, {\tau})$	$\varphi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.6	$\varphi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.8	$(NTDM_{ABC})$ at ${ \varrho }=$ 1	$(NTDM_{CF})$ at ${ \varrho }=$ 1	Exact solution
(0.2, 0.01)	13.691122	13.691178	13.691260	13.691260	13.691260
(0.4, 0.02)	8.946002	8.946023	8.946070	8.946070	8.946070
(0.6, 0.03)	4.881079	4.881103	4.881133	4.881133	4.881133
(0.2, 0.01)	13.710837	13.710898	13.710995	13.710995	13.710995
(0.4, 0.02)	8.968516	8.968579	8.968653	8.968653	8.968653
(0.6, 0.03)	4.896526	4.896588	4.896615	4.896615	4.896615
(0.2, 0.01)	13.693268	13.693302	13.693340	13.693340	13.693340
(0.4, 0.02)	8.948378	8.948412	8.948446	8.948446	8.948446
(0.6, 0.03)	4.882679	4.882705	4.882761	4.882761	4.882761
(0.2, 0.01)	13.694302	13.694323	13.694380	13.694380	13.694380
(0.4, 0.02)	8.949587	8.949612	8.949634	8.949634	8.949634
(0.6, 0.03)	4.883501	4.883525	4.883575	4.883575	4.883575
(0.2, 0.01)	13.695389	13.695402	13.695420	13.695420	13.695420
(0.4, 0.02)	8.950736	8.950778	8.950823	8.950823	8.950823
(0.6, 0.03)	4.884288	4.884326	4.884389	4.884389	4.884389

| Show Table

DownLoad: CSV

Table 3. Absolute Error (AE) comparison for

$\psi({\upsilon}, {\tau})$ at

${\varrho = 1}$ obtained by different techniques.

$({\upsilon}, {\tau})$	AE Of ADM ^[50]	AE Of VIM ^[51]	AE Of OHAM ^[52]	AE of $(NTDM_{ABC})$	AE of $(NTDM_{CF})$
(0.1, 0.2)	1.05983 ${ \times }10^{-5}$	1.34144 ${ \times }10^{-6}$	1.18169 ${ \times }10^{-7}$	1.45621 ${ \times }10^{-10}$	1.45621 ${ \times }10^{-10}$
(0.1, 0.4)	9.75585 ${ \times }10^{-6}$	3.78688 ${ \times }10^{-6}$	3.15656 ${ \times }10^{-7}$	2.95673 ${ \times }10^{-09}$	2.95673 ${ \times }10^{-09}$
(0.1, 0.6)	8.77423 ${ \times }10^{-6}$	6.27984 ${ \times }10^{-6}$	4.92412 ${ \times }10^{-7}$	1.10432 ${ \times }10^{-8}$	1.10432 ${ \times }10^{-8}$
(0.2, 0.2)	4.37319 ${ \times }10^{-5}$	1.38978 ${ \times }10^{-6}$	1.12395 ${ \times }10^{-7}$	1.44301 ${ \times }10^{-10}$	1.44301 ${ \times }10^{-10}$
(0.2, 0.4)	3.82189 ${ \times }10^{-5}$	3.51189 ${ \times }10^{-6}$	2.86457 ${ \times }10^{-6}$	2.35213 ${ \times }10^{-09}$	2.35213 ${ \times }10^{-09}$
(0.2, 0.6)	3.51272 ${ \times }10^{-5}$	6.10117 ${ \times }10^{-5}$	4.51245 ${ \times }10^{-6}$	1.02145 ${ \times }10^{-08}$	1.02145 ${ \times }10^{-08}$
(0.3, 0.2)	9.62833 ${ \times }10^{-5}$	1.25698 ${ \times }10^{-5}$	1.13664 ${ \times }10^{-6}$	2.27511 ${ \times }10^{-10}$	2.27511 ${ \times }10^{-10}$
(0.3, 0.4)	8.84418 ${ \times }10^{-5}$	3.61977 ${ \times }10^{-5}$	2.62353 ${ \times }10^{-6}$	2.02314 ${ \times }10^{-09}$	2.02314 ${ \times }10^{-09}$
(0.3, 0.6)	8.33563 ${ \times }10^{-5}$	5.96721 ${ \times }10^{-5}$	4.46642 ${ \times }10^{-6}$	1.03541 ${ \times }10^{-08}$	1.03541 ${ \times }10^{-08}$
(0.4, 0.2)	1.86687 ${ \times }10^{-4}$	1.24938 ${ \times }10^{-5}$	9.24537 ${ \times }10^{-5}$	1.32601 ${ \times }10^{-10}$	1.32601 ${ \times }10^{-10}$
(0.4, 0.4)	1.72542 ${ \times }10^{-4}$	3.52859 ${ \times }10^{-5}$	2.63564 ${ \times }10^{-5}$	1.11427 ${ \times }10^{-09}$	1.11427 ${ \times }10^{-09}$
(0.4, 0.6)	1.58687 ${ \times }10^{-4}$	5.81821 ${ \times }10^{-5}$	4.65446 ${ \times }10^{-5}$	1.87965 ${ \times }10^{-08}$	1.87965 ${ \times }10^{-08}$
(0.5, 0.2)	2.88628 ${ \times }10^{-4}$	1.21847 ${ \times }10^{-5}$	9.72736 ${ \times }10^{-5}$	1.78145 ${ \times }10^{-10}$	1.78145 ${ \times }10^{-10}$
(0.5, 0.4)	2.47825 ${ \times }10^{-4}$	3.44373 ${ \times }10^{-5}$	2.33457 ${ \times }10^{-5}$	1.43901 ${ \times }10^{-09}$	1.43901 ${ \times }10^{-09}$
(0.5, 0.6)	2.47295 ${ \times }10^{-4}$	5.47346 ${ \times }10^{-5}$	4.38895 ${ \times }10^{-5}$	1.43421 ${ \times }10^{-08}$	1.43421 ${ \times }10^{-08}$

| Show Table

DownLoad: CSV

Table 4. Absolute Error (AE) comparison for

$\varphi({\upsilon}, {\tau})$ at

${\varrho = 1}$ obtained by different techniques.

$({\upsilon, \tau})$	AE Of ADM ^[50]	AE Of VIM ^[51]	AE Of OHAM ^[52]	AE of $(NTDM_{ABC})$	AE of $(NTDM_{CF})$
(0.1, 0.2)	6.52318 ${ \times }10^{-4}$	1.23581 ${ \times }10^{-5}$	5.72451 ${ \times }10^{-6}$	1.73281 ${ \times }10^{-10}$	1.73281 ${ \times }10^{-10}$
(0.1, 0.4)	5.87694 ${ \times }10^{-4}$	3.53456 ${ \times }10^{-5}$	3.24632 ${ \times }10^{-6}$	2.00931 ${ \times }10^{-09}$	2.00931 ${ \times }10^{-09}$
(0.1, 0.6)	5.72618 ${ \times }10^{-4}$	5.63261 ${ \times }10^{-5}$	3.38923 ${ \times }10^{-6}$	1.08631 ${ \times }10^{-08}$	1.08631 ${ \times }10^{-08}$
(0.2, 0.2)	1.44292 ${ \times }10^{-3}$	1.18127 ${ \times }10^{-5}$	5.45771 ${ \times }10^{-6}$	1.11621 ${ \times }10^{-10}$	1.11621 ${ \times }10^{-10}$
(0.2, 0.4)	1.33452 ${ \times }10^{-3}$	3.34512 ${ \times }10^{-5}$	2.86341 ${ \times }10^{-5}$	2.07941 ${ \times }10^{-09}$	2.07941 ${ \times }10^{-09}$
(0.2, 0.6)	1.25527 ${ \times }10^{-3}$	5.47838 ${ \times }10^{-5}$	2.82545 ${ \times }10^{-5}$	1.03361 ${ \times }10^{-08}$	1.03361 ${ \times }10^{-08}$
(0.3, 0.2)	2.14563 ${ \times }10^{-3}$	1.14848 ${ \times }10^{-5}$	5.36746 ${ \times }10^{-5}$	1.16031 ${ \times }10^{-10}$	1.16031 ${ \times }10^{-10}$
(0.3, 0.4)	1.84963 ${ \times }10^{-3}$	3.22828 ${ \times }10^{-5}$	2.74231 ${ \times }10^{-5}$	2.006841 ${ \times }10^{-09}$	2.006841 ${ \times }10^{-09}$
(0.3, 0.6)	1.72318 ${ \times }10^{-3}$	5.32558 ${ \times }10^{-4}$	2.66463 ${ \times }10^{-5}$	2.86931 ${ \times }10^{-08}$	2.86931 ${ \times }10^{-08}$
(0.4, 0.2)	2.98211 ${ \times }10^{-3}$	1.11468 ${ \times }10^{-4}$	5.23838 ${ \times }10^{-5}$	4.09741 ${ \times }10^{-10}$	4.09741 ${ \times }10^{-10}$
(0.4, 0.4)	2.59845 ${ \times }10^{-3}$	3.13456 ${ \times }10^{-4}$	2.72338 ${ \times }10^{-3}$	3.16439 ${ \times }10^{-09}$	3.16439 ${ \times }10^{-09}$
(0.4, 0.6)	2.61896 ${ \times }10^{-3}$	5.15382 ${ \times }10^{-3}$	2.54328 ${ \times }10^{-3}$	1.00693 ${ \times }10^{-09}$	1.00693 ${ \times }10^{-09}$
(0.5, 0.2)	3.84384 ${ \times }10^{-3}$	9.86396 ${ \times }10^{-3}$	4.83832 ${ \times }10^{-3}$	1.00631 ${ \times }10^{-10}$	1.00631 ${ \times }10^{-10}$
(0.5, 0.4)	3.58728 ${ \times }10^{-3}$	2.84228 ${ \times }10^{-3}$	2.84563 ${ \times }10^{-3}$	2.12692 ${ \times }10^{-09}$	2.12692 ${ \times }10^{-09}$
(0.5, 0.6)	3.35348 ${ \times }10^{-3}$	4.72446 ${ \times }10^{-3}$	2.52741 ${ \times }10^{-3}$	1.00471 ${ \times }10^{-09}$	1.00471 ${ \times }10^{-09}$

| Show Table

DownLoad: CSV

Table 5. The suggested techniques result for

$\psi({\upsilon}, {\tau})$ at various fractional-orders of Example 2.

$({\upsilon}, {\tau})$	$\psi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.5	$\psi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.75	$(NTDM_{ABC})$ at ${ \varrho }=$ 1	$(NTDM_{CF})$ at ${ \varrho }=$ 1	Exact result
(0.2, 0.01)	-0.124902	-0.124896	-0.124892	-0.124892	-0.124892
(0.4, 0.01)	-0.123585	-0.123568	-0.123553	-0.123553	-0.123553
(0.6, 0.01)	-0.122297	-0.122292	-0.122280	-0.122280	-0.122280
(0.2, 0.02)	-0.124908	-0.124898	-0.124892	-0.124892	-0.124892
(0.4, 0.02)	-0.123597	-0.123578	-0.123553	-0.123553	-0.123553
(0.6, 0.02)	-0.122299	-0.122293	-0.122280	-0.122280	-0.122280
(0.2, 0.03)	-0.124909	-0.124899	-0.124892	-0.124892	-0.124892
(0.4, 0.03)	-0.123586	-0.123568	-0.123553	-0.123553	-0.123553
(0.6, 0.03)	-0.122299	-0.122289	-0.122280	-0.122280	-0.122280
(0.2, 0.04)	-0.124908	-0.124896	-0.124892	-0.124892	-0.124892
(0.4, 0.04)	-0.123589	-0.123576	-0.123553	-0.123553	-0.123553
(0.6, 0.04)	-0.122296	-0.122288	-0.122280	-0.122280	-0.122280
(0.2, 0.05)	-0.124905	-0.124897	-0.124892	-0.124892	-0.124892
(0.4, 0.05)	-0.123578	-0.123564	-0.123553	-0.123553	-0.123553
(0.6, 0.05)	-0.122298	-0.122290	-0.122280	-0.122280	-0.122280

| Show Table

DownLoad: CSV

Table 6. The suggested techniques result for

$\varphi({\upsilon}, {\tau})$ at various fractional-orders of Example 2.

$({\upsilon}, {\tau})$	$\varphi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.5	$\varphi({\upsilon}, {\tau})$ at ${ \varrho }=$ 0.75	$(NTDM_{ABC})$ at ${ \varrho }=$ 1	$(NTDM_{CF})$ at ${ \varrho }=$ 1	Exact result
(0.2, 0.01)	-0.006894	-0.006886	-0.006872	-0.006872	-0.006872
(0.4, 0.01)	-0.006546	-0.006537	-0.006525	-0.006525	-0.006525
(0.6, 0.01)	-0.006219	-0.006213	-0.006200	-0.006200	-0.006200
(0.2, 0.02)	-0.006894	-0.006885	-0.006872	-0.006872	-0.006872
(0.4, 0.02)	-0.006542	-0.006534	-0.006525	-0.006525	-0.006525
(0.6, 0.02)	-0.006219	-0.006213	-0.006200	-0.006200	-0.006200
(0.2, 0.03)	-0.006889	-0.006882	-0.006872	-0.006872	-0.006872
(0.4, 0.03)	-0.006543	-0.006536	-0.006525	-0.006525	-0.006525
(0.6, 0.03)	-0.006216	-0.006208	-0.006200	-0.006200	-0.006200
(0.2, 0.04)	-0.006889	-0.006881	-0.006872	-0.006872	-0.006872
(0.4, 0.04)	-0.006541	-0.006532	-0.006525	-0.006525	-0.006525
(0.6, 0.04)	-0.006218	-0.006207	-0.006200	-0.006200	-0.006200
(0.2, 0.05)	-0.006889	-0.006883	-0.006872	-0.006872	-0.006872
(0.4, 0.05)	-0.006546	-0.006538	-0.006525	-0.006525	-0.006525
(0.6, 0.05)	-0.006223	-0.006213	-0.006200	-0.006200	-0.006200

| Show Table

DownLoad: CSV

Figure 1. The actual and analytic solutions for

$\psi({\upsilon}, {\tau})$ at

$\varrho = 1$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. The analytic result for

$\psi({\upsilon}, {\tau})$ at

$\varrho = 0.8, 0.6$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. The analytic solution graph at different value of

$\varrho$ for

$\psi({\upsilon}, {\tau})$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. The absolute error for

$\psi({\upsilon}, {\tau})$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 5. The actual and analytic result for

$\varphi({\upsilon}, {\tau})$ at

$\varrho = 1$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 6. The analytic result graph for

$\varphi({\upsilon}, {\tau})$ at

$\varrho = 0.8, 0.6$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 7. The analytic result graph at different value of

$\varrho$ for

$\varphi({\upsilon}, {\tau})$ of Example

$1$ .

DownLoad: Full-Size Img PowerPoint

Figure 8. The exact and analytic solutions graph for

$\psi({\upsilon}, {\tau})$ at

$\varrho = 1$ of Example

$2$ .

DownLoad: Full-Size Img PowerPoint

Figure 9. The absolute error for

$\psi({\upsilon}, {\tau})$ of Example

$2$ .

DownLoad: Full-Size Img PowerPoint

Figure 10. The exact and analytic solutions for

$\varphi({\upsilon}, {\tau})$ at

$\varrho = 1$ of Example

$2$ .

DownLoad: Full-Size Img PowerPoint

Figure 11. The absolute error for

$\varphi({\upsilon}, {\tau})$ of Example

$2$ .

DownLoad: Full-Size Img PowerPoint

7. Conclusions

The Atangana-Baleanu and Caputo-Fabrizio operators are used in this work to attempt a semi-analytic solution of fractional Whitham-Broer-Kaup equations. To show and confirm the efficiency of the recommended technique, two examples are solved. The numerical results show that the proposed method for solving time-fractional Whitham-Broer-Kaup equations is quite effective and accurate. According to numerical data, the method is very effective and reliable for getting close solutions for nonlinear fractional partial differential equations. Compared to other analytical methods, the proposed method is a quick and easy way to look into the numerical solution of nonlinear coupled systems of fractional partial differential equations. The proposed method gives solutions in the form of a series that is more accurate and take less time to figure out. The calculated results have been displayed graphically and in tables. For both considerably coupled systems, computations were done to determine the absolute error. Several computational solutions are contrasted with well-known analytical methods and the precise results at $\varrho = 1$ . Fewer calculations and more precision are two advantages of the present approaches. Lastly, we can say that the proposed approaches are very effective and useful and that they can be used to study any nonlinear problems that come up in complex phenomena.

Conflict of interest

The authors declare that they have no competing interests.

References

[1]	F. Mainardi, Fractional calculus and waves in linear viscoelasticity, Singapore: World Scientific, 2010. https://doi.org/10.1142/p614
[2]	A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Preface, North-Holland Math. Stud., 204 (2006), 7–10. https://doi.org/10.1016/s0304-0208(06)80001-0
[3]	G. W. Scott Blair, The role of psychophysics in rheology, J. Colloid Sci., 2 (1947), 21–32. https://doi.org/10.1016/0095-8522(47)90007-x doi: 10.1016/0095-8522(47)90007-x
[4]	D. Alpay, H. Dym, On a new class of reproducing kernel spaces and a new generalization of the iohvidov laws, Linear Algebra Appl., 178 (1993), 109–183. https://doi.org/10.1016/0024-3795(93)90339-p doi: 10.1016/0024-3795(93)90339-p
[5]	M. Caputo, F. Mainardi, Linear models of dissipation in Anelastic solids, La Riv. Nuovo Cim., 1 (1971), 161–198. https://doi.org/10.1007/bf02820620 doi: 10.1007/bf02820620
[6]	K. L. Wang, Exact travelling wave solution for the local fractional Camassa-Holm-Kadomtsev-Petviashvili equation, Alex. Eng. J., 63 (2023), 371–376. https://doi.org/10.1016/j.aej.2022.08.011 doi: 10.1016/j.aej.2022.08.011
[7]	K. J. Wang, J. Si, On the non-differentiable exact solutions of the (2+1)-dimensional local fractional breaking soliton equation on Cantor sets, Math. Method. Appl. Sci., 2022. https://doi.org/10.1002/mma.8588
[8]	F. Gao, X. J. Yang, Fractional Maxwell fluid with fractional derivative without singular kernel, Therm. Sci., 20 (2016), 871–877. https://doi.org/10.2298/tsci16s3871g doi: 10.2298/tsci16s3871g
[9]	M. Al-Refai, K. Pal, New aspects of caputo-fabrizio fractional derivative, Prog. Fract. Differ. Appl., 5 (2019), 157–166. https://doi.org/10.18576/pfda/050206 doi: 10.18576/pfda/050206
[10]	A. Atangana, B. S. Alkahtani, Extension of the resistance, inductance, capacitance electrical circuit to fractional derivative without singular kernel, Adv. Mech. Eng., 7 (2015). https://doi.org/10.1177/1687814015591937
[11]	J. F. Gómez-Aguilar, M. G. López-Lópezb, V. M. Alvarado-Martínezb, J. Reyes-Reyesb, M. Adam-Medina, Modeling diffusive transport with a fractional derivative without singular kernel, Phy. A, 447 (2016), 467–481. https://doi.org/10.1016/j.physa.2015.12.066 doi: 10.1016/j.physa.2015.12.066
[12]	K. M. Owolabi, A. Atangana, Robustness of fractional difference schemes via the Caputo subdiffusion-reaction equations, Chaos Solitons Fract., 111 (2018), 119–127. https://doi.org/10.1016/j.chaos.2018.04.019 doi: 10.1016/j.chaos.2018.04.019
[13]	K. W. Owolabi, Robust and adaptive techniques for numerical simulation of nonlinear partial differential equations of fractional order, Commun. Nonlinear Sci., 44 (2017), 304–317. https://doi.org/10.1016/j.cnsns.2016.08.021 doi: 10.1016/j.cnsns.2016.08.021
[14]	H. X. Chen, M. M. Liu, Y. T. Chen, S. Y. Li, Y. Z. Miao, Nonlinear lamb wave for structural incipient defect detection with sequential probabilistic ratio test, Secur. Commun. Netw., 12 (2022), 9851533. https://doi.org/10.1155/2022/9851533 doi: 10.1155/2022/9851533
[15]	H. Chen, S. Li, Collinear nonlinear mixed-frequency ultrasound with FEM and experimental method for structural health prognosis, Processes, 10 (2022), 656. https://doi.org/10.3390/pr10040656 doi: 10.3390/pr10040656
[16]	X. Qin, L. Zhang, L. Yang, S. Cao, Heuristics to sift extraneous factors in Dixon resultants, J. Symb. Comput., 112 (2021), 105–121. https://doi.org/10.1016/j.jsc.2022.01.003 doi: 10.1016/j.jsc.2022.01.003
[17]	H. Beyer, S. Kempfle, Definition of physically consistent damping laws with fractional derivatives, ZAMM J. Appl. Math. Mech., 75 (1995), 623–635. https://doi.org/10.1002/zamm.19950750820 doi: 10.1002/zamm.19950750820
[18]	W. Okrasinski, S. Vila, Approximations of solutions to some second order nonlinear differential equations, Nonlinear Anal.-Theor., 35 (1999), 1061–1072. https://doi.org/10.1016/s0362-546x(99)80001-7 doi: 10.1016/s0362-546x(99)80001-7
[19]	M. Caputo, Linear models of dissipation whose Q is almost frequency independent–II, Geophys. J. Int., 13 (1967), 529–539. https://doi.org/10.1111/j.1365-246x.1967.tb02303.x doi: 10.1111/j.1365-246x.1967.tb02303.x
[20]	V. N. Kovalnogov, R. V. Fedorov, T. V. Karpukhina, T. E. Simos, C. Tsitouras, Runge-kutta pairs of orders 5(4) trained to best address keplerian type orbits, Mathematics, 9 (2021), 2400. https://doi.org/ 10.3390/math9192400 doi: 10.3390/math9192400
[21]	V. N. Kovalnogov, R. V. Fedorov, T. V. Karpukhina, T. E. Simos, C. Tsitouras, Sixth order numerov- type methods with coefficients trained to perform best on problems with oscillating solutions, Mathematics, 9 (2021), 2756. https://doi.org/10.3390/math9212756 doi: 10.3390/math9212756
[22]	R. Ye, P. Liu, K. Shi, B. Yan, State damping control: a novel simple method of rotor UAV with high performance, IEEE Access, 8 (2020), 214346–214357. https://doi.org/10.1109/ACCESS.2020.3040779 doi: 10.1109/ACCESS.2020.3040779
[23]	K. Nonlaopon, A. M. Alsharif, A. M. Zidan, A. Khan, Y. S. Hamed, R. Shah, Numerical investigation of fractional-order Swift-Hohenberg equations via a Novel transform, Symmetry, 13 (2021), 1263. https://doi.org/10.3390/sym13071263 doi: 10.3390/sym13071263
[24]	M. K. Alaoui, R. Fayyaz, A. Khan, R. Shah, M. S. Abdo, Analytical investigation of Noyes-Field model for time-fractional Belousov-Zhabotinsky reaction, Complexity, 2021 (2021), 1–21. https://doi.org/10.1155/2021/3248376 doi: 10.1155/2021/3248376
[25]	M. Areshi, A. Khan, R. Shah, K. Nonlaopon, Analytical investigation of fractional-order Newell-Whitehead-Segel Equations via a novel transform, AIMS Math., 7 (2022), 6936–6958. https://doi.org/10.3934/math.2022385 doi: 10.3934/math.2022385
[26]	A. M. Zidan, A. Khan, R. Shah, M. K. Alaoui, W. Weera, Evaluation of time-fractional Fisher's equations with the help of Analytical Methods, AIMS Math., 7 (2022), 18746–18766. https://doi.org/10.3934/math.20221031 doi: 10.3934/math.20221031
[27]	M. M. Al-Sawalha, A. S. Alshehry, K. Nonlaopon, R. Shah, O. Y. Ababneh, Fractional view analysis of delay differential equations via numerical method, AIMS Math., 7 (2022), 20510–20523. https://doi.org/10.3934/math.20221123 doi: 10.3934/math.20221123
[28]	M. M. Al-Sawalha, A. S. Alshehry, K. Nonlaopon, R. Shah, O. Y. Ababneh, Approximate analytical solution of time-fractional vibration equation via reliable numerical algorithm, AIMS Math., 7 (2022), 19739–19757. https://doi.org/10.3934/math.20221082 doi: 10.3934/math.20221082
[29]	M. M. Al-Sawalha, R. P. Agarwal, R. Shah, O. Y. Ababneh, W. Weera, A reliable way to deal with fractional-order equations that describe the unsteady flow of a polytropic gas, Mathematics, 10 (2022), 2293. https://doi.org/10.3390/math10132293 doi: 10.3390/math10132293
[30]	G. B. Whitham, Variational methods and applications to water waves, Proc. Roy. Soc. London. Ser. A Math. Phy. Sci., 299 (1967), 6–25. https://doi.org/10.1098/rspa.1967.0119 doi: 10.1098/rspa.1967.0119
[31]	D. J. Kaup, A higher-order water-wave equation and the method for solving it, Prog. Theor. Phys., 54 (1975), 396–408. https://doi.org/10.1143/ptp.54.396 doi: 10.1143/ptp.54.396
[32]	B. Zheng, C. Wen, Exact solutions for fractional partial differential equations by a new fractional sub-equation method, Adv. Differ. Equ., 2013 (2013). https://doi.org/10.1186/1687-1847-2013-199
[33]	S. M. El-Sayed, D. Kaya, Exact and numerical traveling wave solutions of Whitham-Broer-Kaup equations, Appl. Math. Comput., 167 (2005), 1339–1349. https://doi.org/10.1016/j.amc.2004.08.012 doi: 10.1016/j.amc.2004.08.012
[34]	S. Saha Ray, A novel method for travelling wave solutions of fractional Whitham-Broer-Kaup, fractional modified Boussinesq and fractional approximate long wave equations in shallow water, Math. Method. Appl. Sci., 38 (2014), 1352–1368. https://doi.org/10.1002/mma.3151 doi: 10.1002/mma.3151
[35]	M. Cui, Compact finite difference method for the fractional diffusion equation, J. Comput. Phy., 228 (2009), 7792–7804. https://doi.org/10.1016/j.jcp.2009.07.021 doi: 10.1016/j.jcp.2009.07.021
[36]	B. Ren, J. Lin, Symmetry reduction related by nonlocal symmetry and explicit solutions for the Whitham-Broer-Kaup system, J. Korean Phy. Soc., 73 (2018), 538–546. https://doi.org/10.3938/jkps.73.538 doi: 10.3938/jkps.73.538
[37]	Q. Huang, G. Huang, H. Zhan, A finite element solution for the fractional advection-dispersion equation, Adv. Water Resour., 31 (2008), 1578–1589. https://doi.org/10.1016/j.advwatres.2008.07.002 doi: 10.1016/j.advwatres.2008.07.002
[38]	B. Zheng, EXP-function method for solving fractional partial differential equations, Sci. World J., 2013 (2013), 1–8. https://doi.org/10.1155/2013/465723 doi: 10.1155/2013/465723
[39]	C. Chiu, F. C. Hoppensteadt, Mathematical models and simulations of bacterial growth and chemotaxis in a diffusion gradient chamber, J. Math. Biol., 42 (2001), 120–144. https://doi.org/10.1007/s002850000069 doi: 10.1007/s002850000069
[40]	H. N. Hassan, M. A. El-Tawil, A new technique of using homotopy analysis method for solving high-order nonlinear differential equations, Math. Method. Appl. Sci., 34 (2010), 728–742. https://doi.org/10.1002/mma.1400 doi: 10.1002/mma.1400
[41]	Z. Zhang, X. Yong, Y. Chen, Symmetry analysis for Whitham-Broer-Kaup equations, J. Nonlinear Math. Phy., 15 (2008), 383. https://doi.org/10.2991/jnmp.2008.15.4.3 doi: 10.2991/jnmp.2008.15.4.3
[42]	N. Shah, Y. Hamed, K. Abualnaja, J. Chung, R. Shah, A. Khan, A comparative analysis of fractional-order Kaup-Kupershmidt equation within different operators, Symmetry, 14 (2022), 986. https://doi.org/10.3390/sym14050986 doi: 10.3390/sym14050986
[43]	M. M. Al-Sawalha1, 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 Math., 7 (2022), 18334–18359. https://doi.org/10.3934/math.20221010 doi: 10.3934/math.20221010
[44]	A. S. Alshehry, M. Imran, A. Khan, R. Shah, W. Weera, Fractional view analysis of Kuramoto-Sivashinsky equations with non-singular kernel operators, Symmetry, 14 (2022), 1463. https://doi.org/10.3390/sym14071463 doi: 10.3390/sym14071463
[45]	M. K. Alaoui, K. Nonlaopon, A. M. Zidan, A. Khan, R. Shah, Analytical investigation of fractional-order Cahn-Hilliard and Gardner equations using two novel techniques, Mathematics, 10 (2022), 1643. https://doi.org/10.3390/math10101643 doi: 10.3390/math10101643
[46]	T. Botmart, R. P. Agarwal, M. Naeem, A. Khan, R. Shah, On the solution of fractional modified Boussinesq and approximate long wave equations with non-singular kernel operators, AIMS Math., 7 (2022), 12483–12513. https://doi.org/10.3934/math.2022693 doi: 10.3934/math.2022693
[47]	M. X. Zhou, A. S. V. Ravi Kanth, K. Aruna, K. Raghavendar, H. Rezazadeh, M. Inc, et al., Numerical solutions of time fractional Zakharov-Kuznetsov equation via natural transform decomposition method with nonsingular kernel derivatives, J. Funct. Space., 2021 (2021), 1–17. https://doi.org/10.1155/2021/9884027 doi: 10.1155/2021/9884027
[48]	A. A. Alderremy, S. Aly, R. Fayyaz, A.Khan, R. Shah, N. Wyal, The analysis of fractional-order nonlinear systems of third order kdv and burgers equations via a novel transform, Complexity, 2022 (2022), 1–24. https://doi.org/10.1155/2022/4935809 doi: 10.1155/2022/4935809
[49]	P. Sunthrayuth, H. A. Alyousef, S. A. El-Tantawy, A. Khan, N. Wyal, Solving fractional-order diffusion equations in a plasma and fluids via a novel transform, J. Funct. Space., 2022 (2022), 1–19. https://doi.org/10.1155/2022/1899130 doi: 10.1155/2022/1899130
[50]	S. M. El-Sayed, D. Kaya, Exact and numerical traveling wave solutions of Whitham-Broer-Kaup equations, Appl. Math. Comput., 167 (2005), 1339–1349. https://doi.org/10.1016/j.amc.2004.08.012 doi: 10.1016/j.amc.2004.08.012
[51]	M. Rafei, H. Daniali, Application of the variational iteration method to the whitham-broer-kaup equations, Comput. Math. Appl., 54 (2007), 1079–1085. https://doi.org/10.1016/j.camwa.2006.12.054 doi: 10.1016/j.camwa.2006.12.054
[52]	S. Haq, M. Ishaq, Solution of coupled whitham-broer-KAUP equations using optimal homotopy asymptotic method, Ocean Eng., 84 (2014), 81–88. https://doi.org/10.1016/j.oceaneng.2014.03.031 doi: 10.1016/j.oceaneng.2014.03.031

This article has been cited by:

1.	Benoumran Telli, Mohammed Said Souid, Ivanka Stamova, Boundary-Value Problem for Nonlinear Fractional Differential Equations of Variable Order with Finite Delay via Kuratowski Measure of Noncompactness, 2023, 12, 2075-1680, 80, 10.3390/axioms12010080
2.	M. Mossa Al-Sawalha, Osama Y. Ababneh, Rasool Shah, Nehad Ali Shah, Kamsing Nonlaopon, Combination of Laplace transform and residual power series techniques of special fractional-order non-linear partial differential equations, 2022, 8, 2473-6988, 5266, 10.3934/math.2023264
3.	M. Mossa Al-Sawalha, Rasool Shah, Kamsing Nonlaopon, Imran Khan, Osama Y. Ababneh, Fractional evaluation of Kaup-Kupershmidt equation with the exponential-decay kernel, 2023, 8, 2473-6988, 3730, 10.3934/math.2023186
4.	M. Mossa Al-Sawalha, Rasool Shah, Kamsing Nonlaopon, Osama Y. Ababneh, Numerical investigation of fractional-order wave-like equation, 2022, 8, 2473-6988, 5281, 10.3934/math.2023265
5.	Aisha Abdullah Alderremy, Rasool Shah, Nehad Ali Shah, Shaban Aly, Kamsing Nonlaopon, Comparison of two modified analytical approaches for the systems of time fractional partial differential equations, 2023, 8, 2473-6988, 7142, 10.3934/math.2023360
6.	Jagdev Singh, Arpita Gupta, Dumitru Baleanu, Computational analysis for fractional model of coupled Whitham-Broer-Kaup equation, 2025, 110, 11100168, 613, 10.1016/j.aej.2024.09.061
7.	Mohsen Ebrahimi, Esmat Sadat Alaviyan Shahri, Alireza Alfi, A graphical method-based Kharitonov theorem for robust stability analysis of incommensurate fractional-order uncertain systems, 2024, 43, 2238-3603, 10.1007/s40314-024-02610-z

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(2400) PDF downloads(176) Cited by(7)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(11) / Tables(6)

AIMS Mathematics

Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators

Related Papers:

Abstract

1. Introduction

2. Basic definitions

3. General implementation of method

3.1. Case I $(NTDM_{CF})$

3.2. Case II $(NTDM_{ABC})$

4. Convergence analysis

5. Applications

5.1. $NTDM_{CF}$ solution

5.2. $NDM_{ABC}$ solution

5.3. $NDM_{CF}$ solution

5.4. $NDM_{ABC}$ solution

6. Numerical results and discussion

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

Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators

Related Papers:

Abstract

1. Introduction

2. Basic definitions

3. General implementation of method

3.1. Case I (NTDMCF) (NTDM_{CF})

3.2. Case II (NTDMABC) (NTDM_{ABC})

4. Convergence analysis

5. Applications

5.1. NTDMCF NTDM_{CF} solution

5.2. NDMABC NDM_{ABC} solution

5.3. NDMCF NDM_{CF} solution

5.4. NDMABC NDM_{ABC} solution

6. Numerical results and discussion

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

3.1. Case I $(NTDM_{CF})$

3.2. Case II $(NTDM_{ABC})$

5.1. $NTDM_{CF}$ solution

5.2. $NDM_{ABC}$ solution

5.3. $NDM_{CF}$ solution

5.4. $NDM_{ABC}$ solution