The Numerical Investigation of a Fractional-Order Multi-Dimensional Model of Navier–Stokes Equation via Novel Techniques

Abstract

In this study, numerical results of a fractional-order multi-dimensional model of the Navier–Stokes equations will be achieved via adoption of two analytical methods, i.e., the Adomian decomposition transform method and the q-Homotopy analysis transform method. The Caputo–Fabrizio operator will be used to define the fractional derivative. The proposed methods will be implemented to provide the series form results of the given models. The series form results of proposed techniques will be validated with the exact results available in the literature. The proposed techniques will be investigated to be efficient, straightforward, and reliable for application to many other scientific and engineering problems.

1. Introduction

A renowned mathematician, Leibnitz, introduced the concept of fractional derivative in 1695. Fractional calculus is associated with non-integral differential and integral operators. The fractional-order differential operator is a non-local operator, implying that its present and prior states determine a system’s subsequent state. The popularity of fractional calculus increases day by day due to its various implementations in a broad area of non-linear complex systems occurring in viscoelasticity, fluid mechanics, life sciences, mathematical biology, physics, and electrochemistry [1,2,3]. The improvements in fractional differential equations have also received a great deal of attention in recent years [4,5]. There is beauty in symmetry analysis, particularly in the study of partial differential equations, and more specifically those equations coming from the Mathematics of Finance. The secret of nature is symmetry, but most observations in nature do not exhibit symmetry. A profound way to hide symmetry is the phenomenon of spontaneous symmetry-breaking. There are two types of symmetries: finite and infinitesimal. Finite symmetries can be discrete or continuous. Parity and time reversal are discrete symmetries of nature, while space, on the other hand, is a continuous transformation. Mathematicians have forever been fascinated by patterns. Classifications of planar patterns and spatial patterns began seriously in the nineteenth century. Unfortunately, finding an accurate solution to non-linear fractional differential equations has proven quite difficult. Fractional derivatives, for example, can be used to describe non-linear earthquake oscillation, and fractional derivatives can help a fluid dynamic traffic model overcome the inadequacies imposed by the continuum traffic-flow assumption [6,7,8,9,10]. Fractional differential equations have piqued the interest of researchers owing to their precise representation of non-linear processes, particularly in nano-hydrodynamics, where the continuum assumption fails, and the fractional model is the best choice. These discoveries sparked a surge of interest in fractional calculus research across many scientific and technical disciplines [11,12,13].

In 1822, the Navier–Stokes (NS) equation, which governs the motion of the viscous fluid flow, was developed [14,15,16]. The equation is a mixture of the energy equation, the continuity equation, and the momentum equation, and may be thought of as Newton’s second law of motion for fluid substances. This equation may represent a variety of physical phenomena, including liquid flow in pipes, ocean currents, airflow, and blood flow over aircraft wings [17,18,19,20]. Salem and El-Shahed [21] were the first to perform fractional modeling of NS equations in 2005. The authors [21] used the Laplace transform, finite Fourier Sine transform, and finite Hankel transforms to generalize the classical NS equations. Kumar et al. [22] solved a non-linear fractional model of the NS equation analytically by linking the Laplace transform and the homotopy perturbation method. Ganji et al. [23] and Ragab et al. [24] used HAM to solve the non-linear fractional NS equation. For numerical calculation of the fractional NS equation, Maitama [25] and Momani and Birajdar [26] used the Adomian decomposition method. Kumar et al. [27] achieved an approximate result of the time-fractional NS problem by combining the Laplace transform and the Adomian decomposition method, whereas Kumar and Chaurasia [28] solved the same equation by coupling the Laplace transform and finite Hankel transform.

Because most non-linear FDEs do not have precise solutions, numerical approaches are needed to estimate their numerical solution. Modeling the dimensions of equations is vital, but so is the reliability of solution methods. It is self-evident that coupling a technique with a transform [29,30] eliminates time-consuming issues and reduces the amount of CPU time required to investigate numerical solutions to non-linear problems. The q-homotopy analysis transform method (q-HATM) [31,32] is a beautiful combination of the Laplace transform and the q-homotopy analysis method. It has the benefit of including powerful computational approaches for investigating FDEs. By properly selecting, it provides a more effortless technique to regulate the convergence area of the series solution in a broad permitted domain. The series solution’s exact sequence and grid point provide more acceptable results. The efficacy of a key in the convergent zone is demonstrated by h and n-curves. The q-HATM has the advantages of not requiring perturbations, linearization, discretization, or any restrictive assumptions, promising a large convergence region, significantly reducing mathematical calculations, not requiring the computation of complicated polynomials providing a non-local effect, and physical parameters [33,34].

One of the most effective analytical strategies for solving linear and non-linear equations is the Laplace decomposition method [35,36]. The Laplace decomposition method offers benefits over other approximation techniques such as perturbation, since it is free of tiny or big parameters. The Laplace decomposition method does not require any linearization and discretization, unlike other analytical approaches. As a result, the LDM outputs are more realistic and efficient. Approximate solutions to a class of non-linear partial and ordinary differential equations have been obtained using this technique [37,38]. The Klein–Gordon equation [39] and the diffusion-wave equations [40] are two examples. This paper applied the q-homotopy analysis transform method and the Adomian decomposition method combined with a Yang transform Caputo–Fabrizio operator for the first time.

The rest of this article is organized as follows. In Section 2, we present some basic definitions and properties. In Section 3, we give the description of the Adomian decomposition transform method for solving fractional partial differential equations and in Section 4, the existence and uniqueness solution for the Adomian decomposition transform technique. Then, in Section 5, we apply this method to establish a two-dimensional NS equation. In Section 6, we discuss the q-homotopy analysis transform method and graphical discussion. The conclusions are presented at the end of the article.

2. Preliminary Concepts

In this part, we address several key ideas, conceptions, and terminologies related to fractional derivative operators involving index and exponential decay as a kernel, as well as the Yang transform’s specific repercussions.

Definition 1.

The Caputo fractional derivative $\left(CFD\right)$ is described as follows [41,42]:

${}_0{}^c\mathbf{D}_{\tau }^{\varpi }\mathbb{W}\left(\tau \right)=\left\{\begin{array}{cc}\frac{1}{\mathsf{\Gamma }(r-\varpi )}{\int }_0^{\tau }\frac{{\mathbb{W}}^{\left(r\right)}\left({\mathbf{y}}_1\right)}{{\left(\tau -{\mathbf{y}}_1\right)}^{\varpi +1-r}}d{\mathbf{y}}_1,\hfill & r-1<\varpi <r\\ \frac{d^r}{d{\tau }^r}\mathbb{W}\left(\tau \right),\hfill & \varpi =r\end{array}\right.$

Definition 2.

The Caputo–Fabrizio fractional derivative operator is described as follows [41,42]:

${}_{}{}^{CF}\mathbf{D}_{\tau }^{\varpi }\left(\mathbb{W}\left(\tau \right)\right)=\frac{(2-\varpi )\mathbb{B}\left(\varpi \right)}{2(1-\varpi )}{\int }_0^{\tau }exp\left(-\frac{\varpi \left(\tau -{\mathbf{y}}_1\right)}{1-\varpi }\right){\mathbb{W}}^{\prime }\left(\tau \right)d\tau$

where $\mathbb{W}\in {\mathbf{H}}^1(\mathbf{a},\mathbf{b})$ (Sobolev space), $\mathbf{a}<\mathbf{b},\varpi \in [0,1]$ and $\mathbb{B}\left(\varpi \right)$ signifies a normalization function as $\mathbb{B}\left(\varpi \right)=\mathbb{B}\left(0\right)=\mathbb{B}\left(1\right)=1.$

Definition 3.

The fractional integral of the Caputo–Fabrizio operator is defined as [41,42]:

${}^{CF}{\mathcal{I}}_{\tau }^{\varpi }\left(\mathbb{W}\left(\tau \right)\right)=\frac{2(1-\varpi )}{(2-\varpi )\mathbb{B}\left(\varpi \right)}\mathbb{W}\left(\tau \right)+\frac{2\varpi }{(2-\varpi )\mathbb{B}\left(\varpi \right)}{\int }_0^{\tau }\mathbb{W}\left({\mathbf{y}}_1\right)d{\mathbf{y}}_1$

Definition 4.

The Yang transform is described as follows [41,42]:

$\mathbb{Y}\left[\mathbb{W}\left(\phi \right)\right]=\mathbb{Y}\left(\omega \right)={\int }_0^{\infty }\mathbb{W}\left(\phi \right)exp\left(-\frac{\phi }{\omega }\right)d\phi ,\phi >0.$

The Yang transform of a range of vital expressions is as follows:

$\begin{array}{cc}& \mathbb{Y}\left[1\right]=\omega \hfill \\ & \mathbb{Y}\left[\phi \right]={\omega }^2\hfill \\ & ⋮\hfill \\ & \mathbb{Y}\left[\frac{{\phi }^{\varpi }}{\mathsf{\Gamma }(\varpi +1)}\right]={\omega }^{\varpi +1}\hfill \end{array}$

Definition 5.

The inverse Yang transform ${\mathbb{Y}}^{-1}$ is defined by

${\mathbb{Y}}^{-1}\left[\mathbb{Y}\left(\omega \right)\right]=h\left(\tau \right)=\frac{1}{2\pi \iota }{\int }_{\alpha -\iota \infty }^{\alpha +\iota \infty }h\left(\frac{1}{\omega }\right)e^{\omega \tau }\omega d\omega =\mathsf{\Sigma }residuesofh\left(\frac{1}{\omega }\right)e^{\omega \Im }\omega .$

Definition 6.

The Yang transform of the CFD operator is mentioned as [41,42]:

$\mathbb{Y}\left\{{}_0^c{\mathbf{D}}_{\tau }^{\varpi }\left(\mathbb{W}\left(\tau \right)\right),\mathfrak{s}\right\}={\phi }^{-\varpi }\mathbf{Q}\left(\mathfrak{s}\right)-\sum _{\kappa =0}^{\varpi -1}{\phi }^{1-\varpi -\kappa }\left(\mathfrak{s}\right){\mathbb{W}}^{\left(\kappa \right)}\left(0\right),r-1<\varpi <r,\phi >0.$

Definition 7.

The Yang transform of the Caputo–Fabrizio fractional derivative operator is stated as [41,42]:

$\mathbb{Y}\left\{{}_0^{CF}{\mathbf{D}}_e^{\varpi }\left(\mathbb{W}\left(\phi \right)\right),\omega \right\}=\frac{\mathbb{Y}\left[\mathbb{W}\right(\phi )-\omega \mathbb{W}(0\left)\right]}{1+\varpi (\omega -1)}$

3. The Procedure of Adomian Decomposition Transform Method (ADTM)

In this section, we presented the procedure of ADTM for fractional PDEs [40,43].

(1)$\begin{array}{cc}\hfill & {}^{CF}{D}_{\tau }^{\varpi }\mu (\chi ,\tau )+{\mathcal{R}}_1(\mu ,\nu )+{\mathcal{N}}_1(\mu ,\nu )-{\mathcal{P}}_1(\chi ,\tau )=0,\hfill \\ \hfill & {}^{CF}{D}_{\tau }^{\varpi }\nu (\chi ,\tau )+{\mathcal{R}}_2(\mu ,\nu )+{\mathcal{N}}_2(\mu ,\nu )-{\mathcal{P}}_2(\chi ,\tau )=0,0<\varpi \le 1,\hfill \end{array}$

with the initial condition

(2)$\mu (\chi ,0)=g_1\left(\chi \right),\nu (\chi ,0)=g_2\left(\chi \right).$

where ${}^{CF}{D}_{\tau }^{\varpi }=\frac{{\partial }^{\varpi }}{\partial {\tau }^{\varpi }}$ is the Caputo–Fabrizio derivative of fractional-order $\varpi$, ${\mathcal{R}}_1$, ${\mathcal{R}}_2$ and ${\mathcal{N}}_1$, ${\mathcal{N}}_2$ are linear and non-linear terms, respectively, and ${\mathcal{P}}_1$, ${\mathcal{P}}_2$ are source functions.

Using the Yang transformation to Equation (1), we obtain

(3)$\begin{array}{cc}\hfill & \mathbb{Y}{[}^{CF}{D}_{\tau }^{\varpi }\mu (\chi ,\tau )]+\mathbb{Y}[{\mathcal{R}}_1(\mu ,\nu )+{\mathcal{N}}_1(\mu ,\nu )-{\mathcal{P}}_1(\chi ,\tau )]=0,\hfill \\ \hfill & \mathbb{Y}{[}^{CF}{D}_{\tau }^{\varpi }\nu (\chi ,\tau )]+\mathbb{Y}[{\mathcal{R}}_2(\mu ,\nu )+{\mathcal{N}}_2(\mu ,\nu )-{\mathcal{P}}_2(\chi ,\tau )]=0.\hfill \end{array}$

(4)$\begin{array}{cc}\hfill & \frac{1}{(1+\varpi (v-1\left)\right)}\mathbb{Y}\left[\mu (\chi ,\tau )\right]-v\mu (\chi ,0)=\mathbb{Y}\left[{\mathcal{P}}_1(\chi ,\tau )\right]-\mathbb{Y}\{{\mathcal{R}}_1(\mu ,\nu )+{\mathcal{N}}_1(\mu ,\nu )\}],\hfill \\ \hfill & \frac{1}{(1+\varpi (v-1\left)\right)}\mathbb{Y}\left[\nu (\chi ,\tau )\right]-v\nu (\chi ,0)=\mathbb{Y}\left[{\mathcal{P}}_2(\chi ,\tau )\right]-\mathbb{Y}\{{\mathcal{R}}_2(\mu ,\nu )+{\mathcal{N}}_2(\mu ,\nu )\}],\hfill \end{array}$

(5)$\begin{array}{cc}\hfill & \mathbb{Y}\left[\mu (\chi ,\tau )\right]=v\mu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left[{\mathcal{P}}_1(\chi ,\tau )\right]-(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_1(\mu ,\nu )+{\mathcal{N}}_1(\mu ,\nu )\}],\hfill \\ \hfill & \mathbb{Y}\left[\nu (\chi ,\tau )\right]=v\nu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left[{\mathcal{P}}_2(\chi ,\tau )\right]-(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_2(\mu ,\nu )+{\mathcal{N}}_2(\mu ,\nu )\}],\hfill \end{array}$

ADTM defines the result of infinite series $\mu (\chi ,\tau )$ and $\nu (\chi ,\tau )$,

(6)$\begin{array}{c}\hfill \mu (\chi ,\tau )=\sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\tau ),\nu (\chi ,\tau )=\sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\tau ),\end{array}$

The non-linear functions defined by Adomian polynomials ${\mathcal{N}}_1$ and ${\mathcal{N}}_2$ are expressed as

(7)${\mathcal{N}}_1(\mu ,\nu )=\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell },{\mathcal{N}}_2(\mu ,\nu )=\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell },$

The Adomian polynomials can be expressed as

(8)$\begin{array}{c}\hfill {\mathcal{A}}_{\ell }=\frac{1}{\ell !}{\left[\frac{{\partial }^{\ell }}{\partial {\lambda }^{\ell }}\left\{{\mathcal{N}}_1\left(\sum _{k=0}^{\infty }{\lambda }^k{\mu }_k,\sum _{k=0}^{\infty }{\lambda }^k{\nu }_k\right)\right\}\right]}_{\lambda =0},\\ \hfill {\mathcal{B}}_{\ell }=\frac{1}{\ell !}{\left[\frac{{\partial }^{\ell }}{\partial {\lambda }^{\ell }}\left\{{\mathcal{N}}_2\left(\sum _{k=0}^{\infty }{\lambda }^k{\mu }_k,\sum _{k=0}^{\infty }{\lambda }^k{\nu }_k\right)\right\}\right]}_{\lambda =0},\end{array}$

Putting Equations (6) and (8) into Equation (5), gives

(9)$\begin{array}{cc}\hfill & \mathbb{Y}\left[\sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\tau )\right]=v\mu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left\{{\mathcal{P}}_1(\chi ,\tau )\right\}-(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_1(\sum _{\ell =0}^{\infty }{\mu }_{\ell },\sum _{\ell =0}^{\infty }{\nu }_{\ell })+\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell }\},\hfill \\ \hfill & \mathbb{Y}\left[\sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\tau )\right]=v\nu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left\{{\mathcal{P}}_2(\chi ,\tau )\right\}-(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_2(\sum _{\ell =0}^{\infty }{\mu }_{\ell },\sum _{\ell =0}^{\infty }{\nu }_{\ell })+\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell }\},\hfill \end{array}$

Using the inverse Yang transform on Equation (9),

(10)$\begin{array}{cc}\hfill & \sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\tau )={\mathbb{Y}}^{-1}[v\mu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left\{{\mathcal{P}}_1(\chi ,\tau )\right\}]-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_1(\sum _{\ell =0}^{\infty }{\mu }_{\ell },\sum _{\ell =0}^{\infty }{\nu }_{\ell })+\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell }\}\right],\hfill \\ \hfill & \sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\tau )={\mathbb{Y}}^{-1}[v\nu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left\{{\mathcal{P}}_2(\chi ,\tau )\right\}]-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_2(\sum _{\ell =0}^{\infty }{\mu }_{\ell },\sum _{\ell =0}^{\infty }{\nu }_{\ell })+\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell }\}\right],\hfill \end{array}$

we expressed the following terms,

(11)$\begin{array}{cc}\hfill & {\mu }_0(\chi ,\tau )={\mathbb{Y}}^{-1}[v\mu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left\{{\mathcal{P}}_1(\chi ,\tau )\right\}],\hfill \\ \hfill & {\nu }_0(\chi ,\tau )={\mathbb{Y}}^{-1}[v\mu (\chi ,0)+(1+\varpi (v-1))\mathbb{Y}\left\{{\mathcal{P}}_2(\chi ,\tau )\right\}],\hfill \end{array}$

$\begin{array}{cc}\hfill & {\mu }_1(\chi ,\tau )=-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_1({\mu }_0,{\nu }_0)+{\mathcal{A}}_0\}\right],\hfill \\ \hfill & {\nu }_1(\chi ,\tau )=-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_2({\mu }_0,{\nu }_0)+{\mathcal{B}}_0\}\right],\hfill \end{array}$

the general for $\ell \ge 1$, is given by

$\begin{array}{cc}\hfill & {\mu }_{\ell +1}(\chi ,\tau )=-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_1({\mu }_{\ell },{\nu }_{\ell })+{\mathcal{A}}_{\ell }\}\right],\hfill \\ \hfill & {\nu }_{\ell +1}(\chi ,\tau )=-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\{{\mathcal{R}}_2({\mu }_{\ell },{\nu }_{\ell })+{\mathcal{B}}_{\ell }\}\right],\hfill \end{array}$

4. Solution of ADTM

In this section, we apply ADTM coupled with Yang transformation for the Caputo–Fabrizio fractional derivative to solve fractional-order two-dimensional NS equations.

Example 1.

Consider the fractional-order two-dimensional NS equation [31]

(12) $\begin{array}{cc}\hfill & {}^{CF}{D}_{\tau }^{\varpi }\left(\mu \right)+\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }=\rho \left[\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right]+q,\hfill \\ \hfill & {}^{CF}{D}_{\tau }^{\varpi }\left(\nu \right)+\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }=\rho \left[\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right]-q,\hfill \end{array}$

with initial conditions

(13) $\left\{\begin{array}{cc}\hfill & \mu (\chi ,\xi ,0)=-sin\left(\chi +\xi \right),\hfill \\ \hfill & \nu (\chi ,\xi ,0)=sin\left(\chi +\xi \right).\hfill \end{array}\right.$

Using the Yang transform on Equation (12), we get

$\begin{array}{cc}\hfill & \mathbb{Y}\left\{\frac{{}^{CF}{\partial }^{\varpi }\mu (\chi ,\xi ,\tau )}{\partial {\tau }^{\varpi }}\right\}=-\mathbb{Y}\left[\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }-\rho \left[\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right]+q\right],\hfill \\ \hfill & \mathbb{Y}\left\{\frac{{}^{CF}{\partial }^{\varpi }\nu (\chi ,\xi ,\tau )}{\partial {\tau }^{\varpi }}\right\}=-\mathbb{Y}\left[\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }-\rho \left[\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right]-q\right],\hfill \end{array}$

$\begin{array}{cc}\hfill & \frac{1}{(1+\varpi (v-1\left)\right)}\mathbb{Y}\left\{\mu (\chi ,\xi ,\tau )-v\mu (\chi ,\xi ,0)\right\}=-\mathbb{Y}\left[\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }-\rho \left[\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right]+q\right],\hfill \\ \hfill & \frac{1}{(1+\varpi (v-1\left)\right)}\mathbb{Y}\left\{\nu (\chi ,\xi ,\tau )-v\nu (\chi ,\xi ,0)\right\}=-\mathbb{Y}\left[\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }-\rho \left[\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right]-q\right].\hfill \end{array}$

The above equations can be written as

(14) $\begin{array}{cc}\hfill & \mathbb{Y}\left\{\mu (\chi ,\xi ,\tau )\right\}=v\left\{\mu (\chi ,\xi ,0)\right\}-(1+\varpi (v-1))\mathbb{Y}\left[\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }-\rho \left[\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right]+q\right],\hfill \\ \hfill & \mathbb{Y}\left\{\nu (\chi ,\xi ,\tau )\right\}=v\left\{\nu (\chi ,\xi ,0)\right\}-(1+\varpi (v-1))\mathbb{Y}\left[\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }-\rho \left[\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right]-q\right].\hfill \end{array}$

Using the inverse Yang transform, we have

(15) $\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=\mu (\chi ,\xi ,0)-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[q\right]\right]-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }-\rho \left(\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right)\right\}\right],\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=\nu (\chi ,\xi ,0)-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[q\right]\right]-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }-\rho \left(\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right)\right\}\right].\hfill \end{array}$

Suppose that the unknown functions $\mu (\chi ,\xi ,\tau )$ and $\nu (\chi ,\xi ,\tau )$ infinite series solution is as follows:

$\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau ),and\hfill \\ & \nu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi ,\tau ).\hfill \end{array}$

Note that $\mu {\mu }_{\chi }={\sum }_{\ell =0}^{\infty }{\mathcal{A}}_{\ell }$, $\mu {\nu }_{\xi }={\sum }_{\ell =0}^{\infty }{\mathcal{B}}_{\ell }$, $\nu {\mu }_{\chi }={\sum }_{\ell =0}^{\infty }\mathcal{C}m$ and $\nu {\nu }_{\xi }={\sum }_{\ell =0}^{\infty }{\mathcal{D}}_{\ell }$ are the Adomian polynomials and the non-linear terms were described. Using some term, Equation (15) can be rewritten in the form

$\begin{array}{cc}\hfill & \sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau )=\mu (\chi ,\xi ,0)+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{q\right\}\right]\hfill \\ \hfill & +{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell }\right)+\rho \left\{\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right\}\right]\right],\hfill \\ \hfill & \sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi ,\tau )=\nu (\chi ,\xi ,0)-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{q\right\}\right]\hfill \\ \hfill & +{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{C}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{D}}_{\ell }\right)+\rho \left\{\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right\}\right]\right].\hfill \end{array}$

(16) $\begin{array}{cc}\hfill \sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau )=& -sin\left(\chi +\xi \right)+q\left\{1+\varpi \tau -\varpi \right\}+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell }\right)\right]\right]\hfill \\ \hfill & +{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[\rho \left\{\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\mu }_{\ell }}{\partial {\chi }^2}+\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\mu }_{\ell }}{\partial {\xi }^2}\right\}\right]\right],\hfill \\ \hfill \sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi ,\tau )=& sin\left(\chi +\xi \right)-q\left\{1+\varpi \tau -\varpi \right\}+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{C}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{D}}_{\ell }\right)\right]\right]+\hfill \\ \hfill & {\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[\rho \left\{\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\nu }_{\ell }}{\partial {\chi }^2}+\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\nu }_{\ell }}{\partial {\xi }^2}\right\}\right]\right].\hfill \end{array}$

According to Equation (8), the Adomian polynomials can be expressed as

$\begin{array}{cc}\hfill & {\mathcal{A}}_0={\mu }_0\frac{\partial {\mu }_0}{\partial \chi },{\mathcal{A}}_1={\mu }_0\frac{\partial {\mu }_1}{\partial \chi }+{\mu }_1\frac{\partial {\mu }_0}{\partial \chi },{\mathcal{B}}_0={\nu }_0\frac{\partial {\mu }_0}{\partial \chi },{\mathcal{B}}_1={\nu }_0\frac{\partial {\mu }_1}{\partial \xi }+{\nu }_1\frac{\partial {\mu }_0}{\partial \xi },\hfill \\ \hfill & {\mathcal{C}}_0={\mu }_0\frac{\partial {\nu }_0}{\partial \chi },{\mathcal{C}}_1={\mu }_0\frac{\partial {\nu }_1}{\partial \chi }+{\mu }_1\frac{\partial {\nu }_0}{\partial \chi },{\mathcal{D}}_0={\nu }_0\frac{\partial {\nu }_0}{\partial \chi },{\mathcal{D}}_1={\nu }_0\frac{\partial {\nu }_1}{\partial \chi }+{\nu }_1\frac{\partial {\nu }_0}{\partial \chi }.\hfill \end{array}$

Thus, we can easily achieve the recursive relationship Equation (16)

$\begin{array}{cc}\hfill & {\mu }_0(\chi ,\xi ,\tau )=-sin\left(\chi +\xi \right)+q\left\{1+\varpi \tau -\varpi \right\},{\nu }_0(\chi ,\xi ,\tau )=sin\left(\chi +\xi \right)-q\left\{1+\varpi \tau -\varpi \right\}.\hfill \end{array}$

For $\ell =0$

$\begin{array}{cc}\hfill & {\mu }_1(\chi ,\xi ,\tau )=sin\left(\chi +\xi \right)2\rho \left\{1+\varpi \tau -\varpi \right\},{\nu }_1(\chi ,\xi ,\tau )=-sin\left(\chi +\xi \right)2\rho \left\{1+\varpi \tau -\varpi \right\}.\hfill \end{array}$

For $\ell =1$

$\begin{array}{cc}\hfill & {\mu }_2(\chi ,\xi ,\tau )=-sin\left(\chi +\xi \right){\left(2\rho \right)}^2\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\},\hfill \\ & {\nu }_2(\chi ,\xi ,\tau )=sin\left(\chi +\xi \right){\left(2\rho \right)}^2\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}.\hfill \end{array}$

For $\ell =2$

$\begin{array}{cc}\hfill & {\mu }_3(\chi ,\xi ,\tau )=sin\left(\chi +\xi \right){\left(2\rho \right)}^3\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\},\hfill \\ & {\nu }_3(\chi ,\xi ,\tau )=-sin\left(\chi +\xi \right){\left(2\rho \right)}^3\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\}.\hfill \\ \hfill & ⋮\hfill \end{array}$

In same method, the remaining ${\mu }_{\ell }$ and ${\nu }_{\ell }$ $(\ell \ge 3)$ components of the YDM solution can be obtained seamlessly. Consequently, we describe the series of alternative solutions as

$\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi )={\mu }_0(\chi ,\xi )+{\mu }_1(\chi ,\xi )+{\mu }_2(\chi ,\xi )+{\mu }_3(\chi ,\xi )+\cdots .\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi )={\nu }_0(\chi ,\xi )+{\nu }_1(\chi ,\xi )+{\nu }_2(\chi ,\xi )+{\nu }_3(\chi ,\xi )+\cdots .\hfill \end{array}$

$\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=-sin\left(\chi +\xi \right)+q\left\{1+\varpi \tau -\varpi \right\}+sin\left(\chi +\xi \right)2\rho \left\{1+\varpi \tau -\varpi \right\}\hfill \\ \hfill & -sin\left(\chi +\xi \right){\left(2\rho \right)}^2\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}\hfill \\ \hfill & +sin\left(\chi +\xi \right){\left(2\rho \right)}^3\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\}-\cdots .\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=sin\left(\chi +\xi \right)-q\left\{1+\varpi \tau -\varpi \right\}-sin\left(\chi +\xi \right)2\rho \left\{1+\varpi \tau -\varpi \right\}\hfill \\ \hfill & +sin\left(\chi +\xi \right){\left(2\rho \right)}^2\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}\hfill \\ \hfill & -sin\left(\chi +\xi \right){\left(2\rho \right)}^3\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\}+\cdots .\hfill \end{array}$

The exact result of Equation (12) at $\varpi =1$ and $q=0$,

(17) $\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=-e^{-2\rho \tau }sin\left(\chi +\xi \right),\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=e^{-2\rho \tau }sin\left(\chi +\xi \right).\hfill \end{array}$

Example 2.

Consider the fractional-order two-dimensional NS equation [31]

(18) $\begin{array}{cc}\hfill & {}^{CF}{D}_{\tau }^{\varpi }\left(\mu \right)+\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }=\rho \left[\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right]+q,\hfill \\ \hfill & {}^{CF}{D}_{\tau }^{\varpi }\left(\nu \right)+\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }=\rho \left[\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right]-q,\hfill \end{array}$

with the initial conditions

(19) $\left\{\begin{array}{cc}\hfill & \mu (\chi ,\xi ,0)=-e^{\chi +\xi },\hfill \\ \hfill & \nu (\chi ,\xi ,0)=e^{\chi +\xi }.\hfill \end{array}\right.$

Using the Yang transform on Equation (18), we get

$\begin{array}{cc}\hfill & \mathbb{Y}\left\{\frac{{}^{CF}{\partial }^{\varpi }\mu }{\partial {\tau }^{\varpi }}\right\}=\mathbb{Y}\left[-\left(\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right\}+q\right],\hfill \\ \hfill & \mathbb{Y}\left\{\frac{{}^{CF}{\partial }^{\varpi }\nu }{\partial {\tau }^{\varpi }}\right\}=\mathbb{Y}\left[-\left(\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right\}-q\right],\hfill \end{array}$

$\begin{array}{cc}\hfill & \frac{1}{(1+\varpi (v-1\left)\right)}\mathbb{Y}\left\{\mu (\chi ,\xi ,\tau )-v\mu (\chi ,\xi ,0)\right\}=\mathbb{Y}\left[-\left(\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right\}+q\right],\hfill \\ \hfill & \frac{1}{(1+\varpi (v-1\left)\right)}\mathbb{Y}\left\{\nu (\chi ,\xi ,\tau )-v\nu (\chi ,\xi ,0)\right\}=\mathbb{Y}\left[-\left(\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right\}-q\right].\hfill \end{array}$

The above equations can be written as

(20) $\begin{array}{cc}\hfill & \mathbb{Y}\left\{\mu (\chi ,\xi ,\tau )\right\}=v\left\{\mu (\chi ,\xi ,0)\right\}+(1+\varpi (v-1))\mathbb{Y}\left[-\left(\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right\}+q\right],\hfill \\ \hfill & \mathbb{Y}\left\{\nu (\chi ,\xi ,\tau )\right\}=v\left\{\nu (\chi ,\xi ,0)\right\}=(1+\varpi (v-1))\mathbb{Y}\left[-\left(\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right\}-q\right].\hfill \end{array}$

Applying the inverse Yang transformation, we obtain

(21) $\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=\mu (\chi ,\xi ,0)+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{q\right\}\right]+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right\}\right]\right],\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=\nu (\chi ,\xi ,0)-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{q\right\}\right]+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }\right)+\rho \left\{\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right\}\right]\right].\hfill \end{array}$

Suppose that the unknown functions $\mu (\chi ,\xi ,\tau )$ and $\nu (\chi ,\xi ,\tau )$ infinite series result as follows:

$\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau ),and\hfill \\ & \nu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi ,\tau ).\hfill \end{array}$

Note that $\mu {\mu }_{\chi }={\sum }_{\ell =0}^{\infty }{\mathcal{A}}_{\ell }$, $\nu {\mu }_{\xi }={\sum }_{\ell =0}^{\infty }{\mathcal{B}}_{\ell }$, $\mu {\nu }_{\chi }={\sum }_{\ell =0}^{\infty }{\mathcal{C}}_{\ell }$ and $\nu {\nu }_{\xi }={\sum }_{\ell =0}^{\infty }{\mathcal{D}}_{\ell }$ are the Adomian polynomials and the non-linear terms were described. Using some term, Equation (21) can be rewritten in the form

$\begin{array}{cc}\hfill & \sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau )=\mu (\chi ,\xi ,0)+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{q\right\}\right]\hfill \\ \hfill & +{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell }\right)+\rho \left\{\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right\}\right]\right],\hfill \\ \hfill & \sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi ,\tau )=\nu (\chi ,\xi ,0)-{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left\{q\right\}\right]\hfill \\ \hfill & +{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{C}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{D}}_{\ell }\right)+\rho \left\{\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right\}\right]\right].\hfill \end{array}$

(22) $\begin{array}{cc}\hfill \sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau )=& -sin\left(\chi +\xi \right)+q\left\{1+\varpi \tau -\varpi \right\}+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{A}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{B}}_{\ell }\right)\right]\right]\hfill \\ \hfill & +{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[\rho \left\{\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\mu }_{\ell }}{\partial {\chi }^2}+\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\mu }_{\ell }}{\partial {\xi }^2}\right\}\right]\right],\hfill \\ \hfill \sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi ,\tau )=& sin\left(\chi +\xi \right)-q\left\{1+\varpi \tau -\varpi \right\}+{\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[-\left(\sum _{\ell =0}^{\infty }{\mathcal{C}}_{\ell }+\sum _{\ell =0}^{\infty }{\mathcal{D}}_{\ell }\right)\right]\right]+\hfill \\ \hfill & {\mathbb{Y}}^{-1}\left[(1+\varpi (v-1))\mathbb{Y}\left[\rho \left\{\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\nu }_{\ell }}{\partial {\chi }^2}+\sum _{\ell =0}^{\infty }\frac{{\partial }^2{\nu }_{\ell }}{\partial {\xi }^2}\right\}\right]\right].\hfill \end{array}$

According to Equation (8), the Adomian polynomials can be expressed as

$\begin{array}{cc}\hfill & {\mathcal{A}}_0={\mu }_0\frac{\partial {\mu }_0}{\partial \chi },{\mathcal{A}}_1={\mu }_0\frac{\partial {\mu }_1}{\partial \chi }+{\mu }_1\frac{\partial {\mu }_0}{\partial \chi },{\mathcal{B}}_0={\nu }_0\frac{\partial {\mu }_0}{\partial \beta },{\mathcal{B}}_1={\nu }_0\frac{\partial {\mu }_1}{\partial \beta }+{\nu }_1\frac{\partial {\mu }_0}{\partial \beta },\hfill \\ \hfill & {\mathcal{C}}_0={\mu }_0\frac{\partial {\nu }_0}{\partial \chi },{\mathcal{C}}_1={\mu }_0\frac{\partial {\nu }_1}{\partial \chi }+{\mu }_1\frac{\partial {\nu }_0}{\partial \chi },{\mathcal{D}}_0={\nu }_0\frac{\partial {\nu }_0}{\partial \chi },{\mathcal{D}}_1={\nu }_0\frac{\partial {\nu }_1}{\partial \chi }+{\nu }_1\frac{\partial {\nu }_0}{\partial \chi }.\hfill \end{array}$

Thus, we can easily achieve the recursive relationship Equation (22)

$\begin{array}{cc}\hfill & {\mu }_0(\chi ,\xi ,\tau )=-e^{\chi +\xi }+q\left\{1+\varpi \tau -\varpi \right\},{\nu }_0(\chi ,\xi ,\tau )=e^{\chi +\xi }-q\left\{1+\varpi \tau -\varpi \right\}.\hfill \end{array}$

For $\ell =0$

$\begin{array}{cc}\hfill & {\mu }_1(\chi ,\xi ,\tau )=e^{\chi +\xi }2\rho \left\{1+\varpi \tau -\varpi \right\},{\nu }_1(\chi ,\xi ,\tau )=-e^{\chi +\xi }2\rho \left\{1+\varpi \tau -\varpi \right\}.\hfill \end{array}$

For $\ell =1$

$\begin{array}{cc}\hfill & {\mu }_2(\chi ,\xi ,\tau )=-e^{\chi +\xi }{\left(2\rho \right)}^2\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\},\hfill \\ \hfill & {\nu }_2(\chi ,\xi ,\tau )=e^{\chi +\xi }{\left(2\rho \right)}^2\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}.\hfill \end{array}$

For $\ell =2$

$\begin{array}{cc}\hfill & {\mu }_3(\chi ,\xi ,\tau )=e^{\chi +\xi }\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\},\hfill \\ \hfill & {\nu }_3(\chi ,\xi ,\tau )=-e^{\chi +\xi }\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\}.\hfill \\ \hfill & ⋮\hfill \end{array}$

Using the same method, the remaining ${\mu }_{\ell }$ and ${\nu }_{\ell }$ $(\ell \ge 3)$ components of the YDM solution can be obtained seamlessly. Consequently, we describe the series of alternative solutions as

$\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\mu }_{\ell }(\chi ,\xi )={\mu }_0(\chi ,\xi )+{\mu }_1(\chi ,\xi )+{\mu }_2(\chi ,\xi )+{\mu }_3(\chi ,\xi )+\cdots .\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=\sum _{\ell =0}^{\infty }{\nu }_{\ell }(\chi ,\xi )={\nu }_0(\chi ,\xi )+{\nu }_1(\chi ,\xi )+{\nu }_2(\chi ,\xi )+{\nu }_3(\chi ,\xi )+\cdots .\hfill \end{array}$

$\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=-e^{\chi +\xi }+q\left\{1+\varpi \tau -\varpi \right\}+e^{\chi +\xi }2\rho \left\{1+\varpi \tau -\varpi \right\}-e^{\chi +\xi }{\left(2\rho \right)}^2\hfill \\ \hfill & \left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}+e^{\chi +\xi }{\left(2\rho \right)}^3\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\}-\cdots .\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=e^{\chi +\xi }-q\left\{1+\varpi \tau -\varpi \right\}-e^{\chi +\xi }2\rho \left\{1+\varpi \tau -\varpi \right\}+e^{\chi +\xi }{\left(2\rho \right)}^2\hfill \\ \hfill & \left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}-e^{\chi +\xi }{\left(2\rho \right)}^3\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\}+\cdots .\hfill \end{array}$

The exact result of Equation (18) at $\varpi =1$ and $q=0$,

(23) $\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )=-e^{\chi +\xi +2\rho \tau },\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=e^{\chi +\xi +2\rho \tau }.\hfill \end{array}$

5. The Methodology of q-Homotopy Analysis Transform Method

Consider a non-linear, non-homogeneous fractional partial differential equation [31];

(24)${}^{CF}{D}_{\tau }^{\varpi }\mu (\chi ,\xi ,\tau )+R\mu (\chi ,\xi ,\tau )+N\mu (\chi ,\xi ,\tau )=f(\chi ,\xi ,\tau ),n-1<\varpi \le n.$

Here, ${}^{CF}{D}_{\tau }^{\varpi }\mu$ is the Caputo–Fabrizio derivative and R represents the linear and N non-linear terms, respectively. $f(\chi ,\xi ,\tau )$ is the source function.

Applying the Yang transform on Equation (24), we get

(25)$\mathbb{Y}\left[\mu \right(\chi ,\xi ,\tau \left)\right]-v\mu (\chi ,\xi ,0)+(1+\varpi (v-1\left)\right)\left\{\mathbb{Y}\right[R\mu (\chi ,\xi ,\tau )]+\mathbb{Y}[NL(\chi ,\xi ,\tau )]-\mathbb{Y}[f(\chi ,\xi ,\tau )\left]\right\}=0.$

The non-linear function is

(26)$\begin{array}{cc}\hfill N\left[\varphi (\chi ,\xi ,\tau ;\overline{q})\right]=& \mathbb{Y}\left[\varphi (\chi ,\xi ,\tau ;\overline{q})\right]-v\varphi (\chi ,\xi ,\tau ;\overline{q})\left(0^{+}\right)+(1+\varpi (v-1))\{\mathbb{Y}\left[R\varphi (\chi ,\xi ,\tau ;\overline{q})\right]\hfill \\ \hfill & +\mathbb{Y}\left[N\varphi (\chi ,\xi ,\tau ;\overline{q})\right]-\mathbb{Y}\left[f(\chi ,\xi ,\tau )\right]\}.\hfill \end{array}$

Here, $\varphi (\chi ,\xi ,\tau ;\overline{q})$ is an unknowns function and $\overline{q}\in \left[0,\frac{1}{4}\right]$ are the embedded parameters, $n\ge 1$. Construct a homotopy as

(27)$\begin{array}{cc}\hfill & (1-n\overline{q})\mathbb{Y}\left[\varphi (\chi ,\xi ,\tau ;\overline{q})-{\mu }_0(\chi ,\xi ,\tau )\right]=ћqH(\chi ,\xi ,\tau )N\left[\varphi (\chi ,\xi ,\tau ;\overline{q})\right]\hfill \end{array}$

where ${\mu }_0$ is an initial condition and $ћ\ne 0$ is an auxiliary parameter.

(28)$\begin{array}{cc}\hfill & \varphi (\chi ,\xi ,\tau ;0)={\mu }_0(\chi ,\xi ,\tau ),\hfill \\ \hfill & \varphi \left(\chi ,\xi ,\tau ;\frac{1}{n}\right)=\mu (\chi ,\xi ,\tau ).\hfill \end{array}$

By calculating $q,\varphi$ convergence type $U_0$ to U and intensifying $\varphi$ about q by Taylor’s theorem, we get

(29)$\begin{array}{c}\hfill \varphi (\chi ,\xi ,\tau ;\overline{q})=U_0+\sum _{\ell =1}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau ){\overline{q}}^m,\end{array}$

where

(30)$\begin{array}{c}\hfill {\mu }_{\ell }={\left.\frac{1}{m!}\frac{{\partial }^m\varphi (\chi ,\xi ,\tau ;\overline{q})}{\partial {\overline{q}}^m}\right|}_{\overline{q}=0}.\end{array}$

With an appropriated selection of auxiliary linear terms, $U_0,n,ћ$ and H, series (29) convergence at $\overline{q}=1/n$, thereby provides a solution

(31)$\begin{array}{c}\hfill \mu (\chi ,\xi ,\tau )={\mu }_0+\sum _{\ell =1}^{\infty }{\mu }_{\ell }(\chi ,\xi ,\tau ){\left(\frac{1}{n}\right)}^m.\end{array}$

Now, differentiate Equation (27) m times, divide by $m!$ and take $\overline{q}=0$

$\mathbb{Y}\left[{\mu }_{\ell }(\chi ,\xi ,\tau )-k_{\ell }{\mu }_{m-1}(\chi ,\xi ,\tau )\right],$

where the vector is described as

(32)$\begin{array}{c}\hfill {\overrightarrow{\mu }}_{\ell }=\left\{{\mu }_0(\chi ,\xi ,\tau ),{\mu }_1(\chi ,\xi ,\tau ),\cdots ,{\mu }_{\ell }(\chi ,\xi ,\tau )\right\}.\end{array}$

Applying the inverse transform on Equation (32), we obtain

(33)$\begin{array}{cc}\hfill {\mu }_{\ell }(\chi ,\xi ,\tau )=& k_{\ell }{\mu }_{\ell -1}(\chi ,\xi ,\tau )+ћ{\mathbb{Y}}^{-1}\left[H(\chi ,\xi ,\tau ){\Re }_{\ell }\left({\overrightarrow{\mu }}_{\ell -1}\right)\right].\hfill \end{array}$

Here,

(34)$\begin{array}{cc}\hfill & {\Re }_{\ell }\left({\overrightarrow{\mu }}_{\ell -1}\right)={\left.\frac{1}{(\ell -1)!}\frac{{\partial }^{\ell -1}N\left[\varphi (\chi ,\xi ,\tau ;\overline{q})\right]}{\partial {\overline{q}}^{\ell -1}}\right|}_{\overline{q}=0},\hfill \end{array}$

and

(35)$\begin{array}{c}\hfill k_r=\left\{\begin{array}{cc}0,\hfill & r\le 1\hfill \\ n,\hfill & r>1\hfill \end{array}\right.\end{array}$

Lastly, by solving Equation (33), the q-homotopy analysis transform method results elements are readily available.

Example 3.

Consider the fractional-order two-dimensional NS equation [31]

(36) $\begin{array}{c}\hfill \left\{\begin{array}{c}{}^{CF}{D}_{\tau }^{\varpi }\mu +\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }={\rho }_0\left(\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right)+g,\hfill \\ {}^{CF}{D}_{\tau }^{\varpi }\nu +\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }={\rho }_0\left(\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right)-g,\hfill \end{array}0<\varpi \le 1,\right.\end{array}$

with initial conditions

(37) $\begin{array}{cc}\hfill & \nu (\chi ,\xi ,0)=sin(\chi +\xi ),\hfill \\ \hfill & \mu (\chi ,\xi ,0)=-sin(\chi +\xi ).\hfill \end{array}$

Applying the Yang transform on Equation (36) and using Equation (37), we get

(38) $\begin{array}{c}\hfill \mathbb{Y}\left[\mu (\chi ,\xi ,\tau )\right]+vsin(\chi +\xi )+(1+\varpi (v-1))\mathbb{Y}\left\{\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right)-g\right\}=0,\\ \hfill \mathbb{Y}\left[\nu (\chi ,\xi ,\tau )\right]-vsin(\chi +\xi )+(1+\varpi (v-1))\mathbb{Y}\left\{\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right)+g\right\}=0.\end{array}$

Define the non-linear operators:

(39) $\begin{array}{cc}\hfill & N^{\mathrm{l}}\left[{\varphi }_1(\chi ,\xi ,\tau ;\overline{q}),{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})\right]=\mathbb{Y}\left[{\varphi }_1(\chi ,\xi ,\tau ;\overline{q})\right]+vsin(\chi +\xi )+(1+\varpi (v-1))\mathbb{Y}\left\{{\varphi }_1(\chi ,\xi ,\tau ;\overline{q})\frac{\partial {\varphi }_1(\chi ,\xi ,\tau ;\overline{q})}{\partial \chi }\right.\hfill \\ \hfill & \left.+{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})\frac{\partial {\varphi }_1(\chi ,\xi ,\tau ;\overline{q})}{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2{\varphi }_1(\chi ,\xi ,\tau ;\overline{q})}{\partial {\chi }^2}\frac{{\partial }^2{\varphi }_1(\chi ,\xi ,\tau ;\overline{q})}{\partial {\xi }^2}\right)-g\right\},\hfill \\ \hfill & N^2\left[{\varphi }_1(\chi ,\xi ,\tau ;\overline{q}),{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})\right]=\mathbb{Y}\left[{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})\right]-vsin(\chi +\xi )+(1+\varpi (v-1))\mathbb{Y}\left\{{\varphi }_1(\chi ,\xi ,\tau ;\overline{q})\frac{\partial {\varphi }_2(\chi ,\xi ,\tau ;\overline{q})}{\partial \chi }\right.\hfill \\ \hfill & \left.+{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})\frac{\partial {\varphi }_2(\chi ,\xi ,\tau ;\overline{q})}{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})}{\partial {\chi }^2}+\frac{{\partial }^2{\varphi }_2(\chi ,\xi ,\tau ;\overline{q})}{\partial {\xi }^2}\right)+g\right\},\hfill \end{array}$

and the Yang operators as

(40) $\begin{array}{cc}\hfill & \mathbb{Y}\left[{\mu }_{\ell }(\chi ,\xi ,\tau )-k_{\ell }{\mu }_{\ell -1}(\chi ,\xi ,\tau )\right]=ћR_{1,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right],\hfill \\ \hfill & \mathbb{Y}\left[{\nu }_{\ell }(\chi ,\xi ,\tau )-k_{\ell }{\nu }_{\ell -1}(\chi ,\xi ,\tau )\right]=ћR_{2,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right],\hfill \end{array}$

(41) $\begin{array}{cc}\hfill & R_{1,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]=\mathbb{Y}\left[{\mu }_{\ell -1}(\chi ,\xi ,\tau )\right]+\left(1-\frac{k_{\ell }}{n}\right)vsin(\chi +\xi )+(1+\varpi (v-1))\mathbb{Y}\left\{\sum _{i=0}^{\ell -1}{\mu }_i\frac{\partial {\mu }_{\ell }-1-i}{\partial \chi }\right.\hfill \\ \hfill & \left.+\sum _{i=0}^{\ell -1}\nu \frac{\partial {\mu }_{\ell -1}-i}{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2{\mu }_{\ell -1}}{\partial {\chi }^2}+\frac{{\partial }^2{\mu }_{\ell -1}}{\partial {\xi }^2}\right)-g\right\}.\hfill \end{array}$

(42) $\begin{array}{cc}\hfill & R_{2,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]=\mathbb{Y}\left[{\nu }_{\ell -1}(\chi ,\xi ,\tau )\right]-\left(1-\frac{k_{\ell }}{n}\right)vsin(\chi +\xi )+(1+\varpi (v-1))\mathbb{Y}\left\{\sum _{i=0}^{\ell -1}{\mu }_i\frac{\partial {\nu }_{\ell -1}-i}{\partial \chi }\right.\hfill \\ \hfill & \left.+\sum _{i=0}^{\ell -1}{\nu }_i\frac{\partial {\nu }_{\ell -1}-i}{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2{\nu }_{\ell -1}}{\partial {\chi }^2}+\frac{{\partial }^2{\nu }_{\ell -1}}{\partial {\xi }^2}\right)+g\right\}.\hfill \end{array}$

Using the inverse Yang transformation on Equation (40), we have

(43) $\begin{array}{cc}\hfill {\mu }_{\ell }(\chi ,\xi ,\tau )=& k_{\ell }{\mu }_{\ell -1}+ћ{\mathbb{Y}}^{-1}\left\{R_{1,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]\right\},\hfill \\ \hfill {\nu }_{\ell }(\chi ,\xi ,\tau )=& k_{\ell }{\nu }_{\ell -1}+ћ{\mathbb{Y}}^{-1}\left\{R_{2,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overline{\nu }}_{\ell -1}\right]\right\}.\hfill \end{array}$

Using ${\mu }_0$ and ${\nu }_0$ in Equation (43), we get

$\begin{array}{cc}\hfill {\mu }_1=& -2{\rho }_0ћsin(\chi +\xi )\left\{1+\varpi \tau -\varpi \right\},\hfill \\ \hfill {\nu }_1=& 2{\rho }_0ћsin(\chi +\xi )\left\{1+\varpi \tau -\varpi \right\},\hfill \\ \hfill {\mu }_2=& -2(n+n){\rho }_0ћsin(\chi +\xi )\left\{1+\varpi \tau -\varpi \right\}-4{\rho }_0^2{ћ}^{2}sin(\chi +\xi )\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\},\hfill \\ \hfill {\nu }_2=& 2(n+ћ){\rho }_0ћsin(\chi +\xi )\left\{1+\varpi \tau -\varpi \right\}+4{\rho }_0^2{ћ}^{2}sin(\chi +\xi )\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\},\hfill \end{array}$

$\begin{array}{cc}\hfill {\mu }_3=& -2{(n+ћ)}^2{\rho }_0ћsin(\chi +\xi )\left\{1+\varpi \tau -\varpi \right\}-8(n+ћ){\rho }_0^2{ћ}^{2}sin(\chi +\xi )\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}\hfill \\ & -8{\rho }_0^3{ћ}^{3}sin(\chi +\xi )\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\},\hfill \\ \hfill {\nu }_3=& 2{(n+ћ)}^2{\rho }_0ћsin(\chi +\xi )\left\{1+\varpi \tau -\varpi \right\}+8(n+ћ){\rho }_0^2{ћ}^{2}sin(\chi +\xi )\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}\hfill \\ & +8{\rho }_0^3{ћ}^{3}sin(\chi +\xi )\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\},\hfill \end{array}$

as well as others. The remaining components are found in the same manner. The following is the q-HATM result of Equation (36):

(44) $\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )={\mu }_0+\sum _{\ell =1}^{\infty }{\mu }_{\ell }{\left(\frac{1}{n}\right)}^{\ell },\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )=2_0+\sum _{w=1}^{\infty }{\nu }_{\ell }{\left(\frac{1}{n}\right)}^{\ell }.\hfill \end{array}$

For $\varpi =1,n=-1,n=1$ and $g=0$, results ${\sum }_{\ell =1}^N{\mu }_{\ell }(\chi ,\xi ,\tau ){(1/n)}^{\ell }$ and ${\sum }_{\ell =1}^N{\nu }_{\ell }(\chi ,\xi ,\tau ){(1/n)}^{\ell }$ convergent to exact solutions as $N\to \infty$:

$\begin{array}{cc}\hfill \mu (\chi ,\xi ,\tau )=& -sin(\chi +\xi )\left[1-\frac{2{\rho }_0\tau }{1!}+\frac{{\left(2{\rho }_0\tau \right)}^2}{2!}-\frac{{\left(2{\rho }_0\tau \right)}^3}{3!}+\cdots \right]\hfill \\ \hfill =& -{\mathrm{e}}^{-2\rho 0\tau }sin(\chi +\xi ),\hfill \\ \hfill \nu (\chi ,\xi ,\tau )=& sin(\chi +\xi )\left[1-\frac{2{\rho }_0\tau }{1!}+\frac{{\left(2{\rho }_0\tau \right)}^2}{2!}-\frac{{\left(2{\rho }_0\tau \right)}^3}{3!}+\cdots \right]\hfill \\ \hfill =& {\mathrm{e}}^{-2{\rho }_0\tau }sin(\chi +\xi ).\hfill \end{array}$

Example 4.

Consider Equation (36) and take the initial conditions

(45) $\begin{array}{c}\hfill \nu (\chi ,\xi ,0)={\mathrm{e}}^{\chi +\xi },\mu (\chi ,\xi ,0)=-{\mathrm{e}}^{\chi +\xi }.\end{array}$

Using the Yang transformation on Equation (36) and applying Equation (45), we have

(46) $\begin{array}{c}\hfill \mathbb{Y}\left[\mu \right]+v{\mathrm{e}}^{\chi +\xi }+(1+\varpi (v-1))\mathbb{Y}\left\{\mu \frac{\partial \mu }{\partial \chi }+\nu \frac{\partial \mu }{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2\mu }{\partial {\chi }^2}+\frac{{\partial }^2\mu }{\partial {\xi }^2}\right)-g\right\}=0,\\ \hfill \mathbb{Y}\left[\nu \right]-v{\mathrm{e}}^{\chi +\xi }+(1+\varpi (v-1))\mathbb{Y}\left\{\mu \frac{\partial \nu }{\partial \chi }+\nu \frac{\partial \nu }{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2\nu }{\partial {\chi }^2}+\frac{{\partial }^2\nu }{\partial {\xi }^2}\right)+g\right\}=0.\end{array}$

Define non-linear operators as

(47) $\begin{array}{cc}\hfill N^1\left[{\varphi }_1,{\varphi }_2\right]=& \mathbb{Y}\left[{\varphi }_1\right]+v{\mathrm{e}}^{\chi +\xi }+(1+\varpi (v-1))\mathbb{Y}\left\{{\varphi }_1\frac{\partial {\varphi }_1}{\partial \chi }+{\varphi }_2\frac{\partial {\varphi }_1}{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2{\varphi }_1}{\partial {\chi }^2}+\frac{{\partial }^2{\varphi }_1}{\partial {\xi }^2}\right)-g\right\},\hfill \\ \hfill N^2\left[{\varphi }_1,{\varphi }_2\right]=& \mathbb{Y}\left[{\varphi }_2\right]-v{\mathrm{e}}^{\chi +\xi }+(1+\varpi (v-1))\mathbb{Y}\left\{{\varphi }_1\frac{\partial {\varphi }_2}{\partial \chi }+{\varphi }_2\frac{\partial {\varphi }_2}{\partial \xi }-{\rho }_0\left(\frac{{\partial }^2{\varphi }_2}{\partial {\chi }^2}+\frac{{\partial }^2{\varphi }_2}{\partial {\xi }^2}\right)+g\right\},\hfill \end{array}$

and Yang operators as

(48) $\begin{array}{cc}\hfill & \mathbb{Y}\left[{\mu }_{\ell }(\chi ,\xi ,\tau )-k_{\ell }{\mu }_{\ell -1}(\chi ,\xi ,\tau )\right]=ћR_{1,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right],\hfill \\ \hfill & \mathbb{Y}\left[{\nu }_{\ell }(\chi ,\xi ,\tau )-k_{\ell }{\nu }_{\ell -1}(\chi ,\xi ,\tau )\right]=ћR_{2,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right],\hfill \end{array}$

where

(49) $\begin{array}{cc}\hfill & R_{1,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]=\mathbb{Y}\left[{\mu }_{\ell -1}\right]+\left(1-\frac{k_{\ell }}{n}\right)\frac{e^{\chi +\xi }}{s}+(1+\varpi (v-1))\mathbb{Y}\left\{\sum _{i=0}^{\ell -1}{\mu }_i\frac{\partial {\mu }_{\ell }-1-i}{\partial \chi }+\sum _{i=0}^{\ell -1}{\nu }_i\frac{\partial {\mu }_{\ell -1}-i}{\partial \xi }\right.\hfill \\ \hfill & \left.-{\rho }_0\left(\frac{{\partial }^2{\mu }_{\ell -1}}{\partial {\chi }^2}+\frac{{\partial }^2{\mu }_{\ell -1}}{\partial {\xi }^2}\right)-g\right\},\hfill \\ \hfill & R_{2,\mu }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]=\mathbb{Y}\left[{\nu }_{\ell -1}\right]-\left(1-\frac{k_{\ell }}{n}\right)v{\mathrm{e}}^{\chi +\xi }+(1+\varpi (v-1))\mathbb{Y}\left\{\sum _{i=0}^{\ell -1}{\mu }_i\frac{\partial {\nu }_{\ell -1-i}}{\partial \chi }+\sum _{i=0}^{\ell -1}{\nu }_i\frac{\partial {\nu }_{\ell -1-i}}{\partial \xi }\right.\hfill \\ \hfill & \left.-{\rho }_0\left(\frac{{\partial }^2{\nu }_{\ell }-1}{\partial {\chi }^2}+\frac{{\partial }^2{\nu }_{\ell }-1}{\partial {\xi }^2}\right)+g\right\}.\hfill \end{array}$

Using the inverse Yang transformation on Equation (48), we get

(50) $\begin{array}{cc}\hfill {\mu }_{\ell }(\chi ,\xi ,\tau )=& k_{\ell }{\mu }_{\ell -1}+ћ{\mathbb{Y}}^{-1}\left\{R_{1,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]\right\},\hfill \\ \hfill {\nu }_{\ell }(\chi ,\xi ,\tau )=& k_{\ell }{\nu }_{\ell -1}+ћ{\mathbb{Y}}^{-1}\left\{R_{2,\ell }\left[{\overrightarrow{\mu }}_{\ell -1},{\overrightarrow{\nu }}_{\ell -1}\right]\right\}.\hfill \end{array}$

Using ${\mu }_0$ and ${\nu }_0$, we get from Equation (50),

${\mu }_1=2{\rho }_0ћ{\mathrm{e}}^{\chi +\xi }{\tau }^{\varpi }\left\{1+\varpi \tau -\varpi \right\},{\nu }_1=-2{\rho }_0ћ{\mathrm{e}}^{\chi +\xi }{\tau }^{\varpi }\left\{1+\varpi \tau -\varpi \right\},$

${\mu }_2=2(n+ћ){\rho }_0ћe^{\chi +\xi }\left\{1+\varpi \tau -\varpi \right\}-4{\rho }_0^2{ћ}^2{e}^{\chi +\xi }\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\},$

${\nu }_2=-2(n+ћ){\rho }_0ћ{\mathrm{e}}^{\chi +\xi }\left\{1+\varpi \tau -\varpi \right\}+4{\rho }_0^2{ћ}^2{\mathrm{e}}^{\chi +\xi }\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\},$

$\begin{array}{cc}\hfill {\mu }_3=& 2{(n+ћ)}^2{\rho }_0ћ{\mathrm{e}}^{\chi +\xi }\left\{1+\varpi \tau -\varpi \right\}-8(n+ћ){\rho }_0^2{ћ}^2{e}^{\chi +\xi }\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}\hfill \\ \hfill & +8{\rho }_0^3{ћ}^3{e}^{\chi +\xi }\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\},\hfill \end{array}$

$\begin{array}{cc}\hfill {\nu }_3=& -2{(n+ћ)}^2{\rho }_0ћe^{\chi +\xi }\left\{1+\varpi \tau -\varpi \right\}+8(n+ћ){\rho }_0^2{ћ}^2{e}^{\chi +\xi }\left\{(1-\varpi )2\varpi \tau +{(1-\varpi )}^2+\frac{{\varpi }^2{\tau }^2}{2}\right\}\hfill \\ & -8{\rho }_0^3{ћ}^3{e}^{\chi +\xi }\left\{{(1-\varpi )}^{2}3\varpi \tau +{(1-\varpi )}^3+\frac{3{\varpi }^2(1-\varpi ){\tau }^2}{2}+\frac{{\varpi }^3{\tau }^3}{3!}\right\},\hfill \end{array}$

as well as others. The remaining components are found in the same manner. The following is the q-HATM result of Equation (36):

(51) $\begin{array}{cc}\hfill & \mu (\chi ,\xi ,\tau )={\mu }_0+\sum _{\ell =1}^{\infty }{\mu }_{\ell }{\left(\frac{1}{n}\right)}^{\ell },\hfill \\ \hfill & \nu (\chi ,\xi ,\tau )={\nu }_0+\sum _{i=1}^{\infty }{\nu }_{\ell }{\left(\frac{1}{n}\right)}^{\ell }.\hfill \end{array}$

For $\varpi =1=n,ћ=-1$ and $g=0$, results ${\sum }_{\ell =1}^N{\mu }_{\ell }{(1/n)}^{\ell }$ and ${\sum }_{\ell =1}^N{\nu }_{\ell }{(1/n)}^{\ell }$ convergence to the actual solutions as $N\to \infty$

$\begin{array}{cc}& \mu (\chi ,\xi ,\tau )=-{\mathrm{e}}^{\chi +\xi }\left[1+\frac{2{\rho }_0\tau }{1!}+\frac{{\left(2{\rho }_0\tau \right)}^2}{2!}+\frac{{\left(2{\rho }_0\tau \right)}^3}{3!}+\cdots \right]\hfill \\ & =-{\mathrm{e}}^{\chi +\xi +2{\rho }_0\tau },\hfill \end{array}$

$\begin{array}{cc}& \nu (\chi ,\xi ,\tau )={\mathrm{e}}^{\chi +\xi }\left[1+\frac{2{\rho }_0\tau }{1!}+\frac{{\left(2{\rho }_0\tau \right)}^2}{2!}+\frac{{\left(2{\rho }_0\tau \right)}^3}{3!}+\cdots \right]\hfill \\ & ={\mathrm{e}}^{\chi +\xi +2{\rho }_0\tau }.\hfill \end{array}$

6. Results and Discussion

The aim of the this paper is to investigate an analytical solution of fractional-order Navier–Stokes equations, applying efficient analytical methods. The Adomian decomposition transform method and q-homotopy analysis transform method are applied to solve the given examples. The Caputo–Fabrizio definition of fractional derivative is applied to define the fractional derivative. To check the accuracy of the proposed methods, the solution of some illustrative problems are presented. Solutions figures are plotted for both integer and fractional-order problems. In Figure 1, the exact and analytical solution $\mu (\chi ,\xi ,\tau )$ of Example 1 is shown. It is seen that the approximate analytical solution obtained by proposed methods decreases very rapidly with the increases in $\tau$ at the value of $\chi =1$. Similarly, in Figure 2, the figures of the exact and analytical solution $\nu (\chi ,\xi ,\tau )$ are represented at $\varpi =1$. It is observed that the exact, ADTM, and q-HATM solutiosn are in close contact with the exact results of the Examples. In Figure 3, Figure 4, Figure 5 and Figure 6 the ADTM and q-HATM solutions of Example 1 are also calculated at different fractional-order $\varpi =0.8$, $0.6$, $0.4$. It is confirmed that ADTM and q-HATM solutions are in strong agreement with each other. Similarly, in Figure 7, the exact and analytical solution $\mu (\chi ,\xi ,\tau )$ of Example 2 is shown. In Figure 8, the graph of exact and analytical solution $\nu (\chi ,\xi ,\tau )$ is represented at $\varpi =1$. It is observed that the exact, ADTM, and q-HATM solutions are in close contact with the exact results of the Examples. Additionally, in Figure 9, Figure 10, Figure 11 and Figure 12 the ADTM and q-HATM solutions of Example 1 are calculated at different fractional-order $\varpi =0.8$, $0.6$, $0.4$. In Table 1 and Table 2 confirmed that ADTM and q-HATM solutions are in strong agreement with each other. The same convergence phenomena of the fractional-order solutions towards integer-order solutions are observed.

7. Conclusions

This paper evaluates a result of the fractional scheme of Navier–Stokes equations calculating numerical results utilizing the proposed q-homotopy analysis transform method and the Yang decomposition method. In a rapid convergence series, the solution is attained. The efficacy and effectiveness of the method are demonstrated by the test samples presented. The suggested approach incorporates a parameter ℏ that controls the convergence zone of the series solution. Because the q-homotopy analysis transform method and the Yang decomposition approach do not require minor perturbations, linearization, or discretization, computation times are greatly reduced. The q-homotopy analysis transform method and Yang decomposition method are competent tools for obtaining mathematical results of system non-linear fractional partial differential equations when compared to other methodologies.

Author Contributions

Formal analysis, S.M.; Methodology, R.S.; Project administration, S.N.; Software, R.S.; Supervision, S.M.; Validation, S.N; Writing—original draft, R.S.; Writing—review & editing, S.M. All authors have equal contribution. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The numerical data used to support the findings of this study are included within the article.

Acknowledgments

This work was supported through the annual funding track by the Deanship of Scientific Research, Vice Presidency of Graduate Studies and Scientific Research King Faisal University Saudi Arabia (project No. AN000377).

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. The figure of actual and q-HATM/ADTM solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 1 and 3.

Figure 2. The figure of actual and q-HATM/ADTM solutions of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 1 and 3.

Figure 3. The various fractional-order solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 1 and 3.

Figure 4. The various fractional-order solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 1 and 3.

Figure 5. The figure of analytic solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] of example 1 and 3.

Figure 6. The figure of analytic solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] of example 1 and 3.

Figure 7. The figure of actual and q-HATM/ADTM solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 2 and 4.

Figure 8. The figure of actual and q-HATM/ADTM solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 2 and 4.

Figure 9. The various fractional-order solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 2 and 4.

Figure 10. The various fractional-order solution of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] of example 2 and 4.

Figure 11. The figure of analytical result of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] of example 2 and 4.

Figure 12. The figure of analytical result of [Forumla omitted. See PDF.] at [Forumla omitted. See PDF.] of example 2 and 4.

Table 1

Comparative study between LDM [27], $q-HATM$ and $ADTM$ for the numerical result $\mu (\chi ,\xi ,\tau )$ of Example 1.

$\mathit{\tau }$	$(\mathit{\chi },\mathit{\xi })$	$|{\mathit{\nu }}_{Exact}-{\mathit{\nu }}_{LDM}|$	$|{\mathit{\nu }}_{Exact}-{\mathit{\nu }}_{ADTM}|$	$|{\mathit{\nu }}_{Exact}-{\mathit{\nu }}_{\mathit{q}-HATM}|$
	0.1	4.3210$\times {10}^{-9}$	4.2464$\times {10}^{-10}$	4.2464$\times {10}^{-10}$
	0.2	5.890$\times {10}^{-9}$	4.5486$\times {10}^{-10}$	4.5486$\times {10}^{-10}$
0.1	0.3	2.7301$\times {10}^{-8}$	2.3578$\times {10}^{-9}$	2.3578$\times {10}^{-9}$
	0.4	5.9700$\times {10}^{-8}$	6.3267$\times {10}^{-9}$	6.3267$\times {10}^{-9}$
	0.5	2.9771$\times {10}^{-7}$	3.2436$\times {10}^{-9}$	3.2436$\times {10}^{-9}$
	0.1	2.4400$\times {10}^{-7}$	4.7421$\times {10}^{-9}$	4.7421$\times {10}^{-9}$
	0.2	8.3100$\times {10}^{-7}$	3.1235$\times {10}^{-8}$	3.1235$\times {10}^{-8}$
0.2	0.3	2.8500$\times {10}^{-6}$	4.5682$\times {10}^{-8}$	4.5682$\times {10}^{-8}$
	0.4	9.7940$\times {10}^{-6}$	3.5223$\times {10}^{-7}$	3.5223$\times {10}^{-7}$
	0.5	3.2012$\times {10}^{-7}$	2.9315$\times {10}^{-7}$	2.9315$\times {10}^{-7}$
	0.1	2.2981$\times {10}^{-6}$	3.2245$\times {10}^{-7}$	3.2245$\times {10}^{-7}$
	0.2	5.4602$\times {10}^{-6}$	4.2659$\times {10}^{-7}$	4.2659$\times {10}^{-7}$
0.3	0.3	2.5432$\times {10}^{-5}$	1.5348$\times {10}^{-6}$	1.5348$\times {10}^{-6}$
	0.4	6.4229$\times {10}^{-4}$	8.2374$\times {10}^{-6}$	8.2374$\times {10}^{-6}$
	0.5	2.8364$\times {10}^{-4}$	4.1975$\times {10}^{-5}$	4.1975$\times {10}^{-5}$
	0.1	5.5428$\times {10}^{-5}$	2.1351$\times {10}^{-6}$	2.1351$\times {10}^{-6}$
	0.2	2.4133$\times {10}^{-5}$	2.6276$\times {10}^{-6}$	2.6276$\times {10}^{-6}$
0.4	0.3	6.3743$\times {10}^{-5}$	2.2334$\times {10}^{-5}$	2.2334$\times {10}^{-5}$
	0.4	2.9070$\times {10}^{-4}$	1.2035$\times {10}^{-5}$	1.2035$\times {10}^{-5}$
	0.5	6.9763$\times {10}^{-4}$	2.2145$\times {10}^{-4}$	2.2145$\times {10}^{-4}$
	0.1	2.2529$\times {10}^{-5}$	2.3223$\times {10}^{-6}$	2.3223$\times {10}^{-6}$
	0.2	4.9868$\times {10}^{-5}$	3.2721$\times {10}^{-5}$	3.2721$\times {10}^{-5}$
0.5	0.3	4.1932$\times {10}^{-4}$	3.0767$\times {10}^{-5}$	3.0767$\times {10}^{-5}$
	0.4	5.5568$\times {10}^{-4}$	2.3742$\times {10}^{-4}$	2.3742$\times {10}^{-4}$
	0.5	2.4350$\times {10}^{-3}$	1.3223$\times {10}^{-3}$	1.3223$\times {10}^{-3}$

Table 2

Comparative study between LDM [27], $q-HATM$ and $ADTM$ for the numerical result $\nu (\chi ,\xi ,\tau )$ of Example 1.

$\mathit{\tau }$	$(\mathit{\chi },\mathit{\xi })$	$|{\mathit{\nu }}_{Exact}-{\mathit{\nu }}_{LDM}|$	$|{\mathit{\nu }}_{Exact}-{\mathit{\nu }}_{ADTM}|$	$|{\mathit{\nu }}_{Exact}-{\mathit{\nu }}_{\mathit{q}-HATM}|$
	0.1	9.0202$\times {10}^{-9}$	1.3770$\times {10}^{-11}$	8.6253$\times {10}^{-11}$
	0.2	5.4060$\times {10}^{-9}$	4.8036$\times {10}^{-10}$	3.1054$\times {10}^{-10}$
0.1	0.3	2.7960$\times {10}^{-8}$	1.6734$\times {10}^{-9}$	1.1992$\times {10}^{-9}$
	0.4	6.5902$\times {10}^{-8}$	5.8013$\times {10}^{-9}$	5.5827$\times {10}^{-9}$
	0.5	3.9762$\times {10}^{-7}$	1.9773$\times {10}^{-8}$	3.6150$\times {10}^{-8}$
	0.1	3.3610$\times {10}^{-8}$	2.3548$\times {10}^{-9}$	2.8894$\times {10}^{-9}$
	0.2	9.1810$\times {10}^{-8}$	8.2143$\times {10}^{-8}$	1.0171$\times {10}^{-8}$
0.2	0.3	1.9482$\times {10}^{-7}$	2.8611$\times {10}^{-8}$	3.6538$\times {10}^{-8}$
	0.4	9.8750$\times {10}^{-7}$	9.2053$\times {10}^{-7}$	1.4010$\times {10}^{-7}$
	0.5	4.2127$\times {10}^{-6}$	3.3727$\times {10}^{-7}$	6.3855$\times {10}^{-7}$
	0.1	2.3872$\times {10}^{-7}$	1.2779$\times {10}^{-8}$	2.3189$\times {10}^{-8}$
	0.2	5.3523$\times {10}^{-7}$	4.4573$\times {10}^{-8}$	8.1227$\times {10}^{-8}$
0.3	0.3	2.5565$\times {10}^{-6}$	1.5522$\times {10}^{-7}$	2.8704$\times {10}^{-7}$
	0.4	6.3787$\times {10}^{-6}$	5.3749$\times {10}^{-6}$	1.0445$\times {10}^{-6}$
	0.5	2.8365$\times {10}^{-5}$	1.8244$\times {10}^{-6}$	4.1512$\times {10}^{-6}$
	0.1	5.2272$\times {10}^{-6}$	4.3423$\times {10}^{-7}$	1.0363$\times {10}^{-7}$
	0.2	2.6215$\times {10}^{-6}$	1.5145$\times {10}^{-7}$	3.6229$\times {10}^{-7}$
0.4	0.3	6.3642$\times {10}^{-6}$	5.2723$\times {10}^{-6}$	1.2717$\times {10}^{-6}$
	0.4	2.9203$\times {10}^{-5}$	1.8244$\times {10}^{-6}$	4.5250$\times {10}^{-6}$
	0.5	6.9958$\times {10}^{-5}$	6.1751$\times {10}^{-5}$	1.6814$\times {10}^{-5}$
	0.1	2.2532$\times {10}^{-6}$	1.1434$\times {10}^{-7}$	3.3630$\times {10}^{-7}$
	0.2	4.8935$\times {10}^{-6}$	3.9879$\times {10}^{-6}$	1.1745$\times {10}^{-6}$
0.5	0.3	2.4921$\times {10}^{-5}$	1.3880$\times {10}^{-6}$	4.1074$\times {10}^{-6}$
	0.4	5.8486$\times {10}^{-5}$	4.7977$\times {10}^{-5}$	1.4434$\times {10}^{-5}$
	0.5	1.6542$\times {10}^{-4}$	1.6182$\times {10}^{-5}$	5.1527$\times {10}^{-5}$

References

1. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl.; 2015; 1, pp. 1-13.

2. Sheikh, N.A.; Ali, F.; Khan, I.; Saqib, M. A modern approach of Caputo-Fabrizio time-fractional derivative to MHD free convection flow of generalized second-grade fluid in a porous medium. Neural Comput. Appl.; 2018; 30, pp. 1865-1875. [DOI: https://dx.doi.org/10.1007/s00521-016-2815-5]

3. Atangana, A.; Alqahtani, R.T. Numerical approximation of the space-time Caputo-Fabrizio fractional derivative and application to groundwater pollution equation. Adv. Differ. Equ.; 2016; 2016, 156. [DOI: https://dx.doi.org/10.1186/s13662-016-0871-x]

4. Aljahdaly, N.; Akgül, A.; Shah, R.; Mahariq, I.; Kafle, J. A Comparative Analysis of the Fractional-Order Coupled Korteweg–De Vries Equations with the Mittag–Leffler Law. J. Math.; 2022; 2022, pp. 1-30. [DOI: https://dx.doi.org/10.1155/2022/8876149]

5. Shamshuddin, M.D.; Mishra, S.R.; Beg, O.A.; Kadir, A. Numerical study of heat transfer and viscous flow in a dual rotating extendable disk system with a non-Fourier heat flux model. Heat Transf. Asian Res.; 2019; 48, pp. 435-459. [DOI: https://dx.doi.org/10.1002/htj.21392]

6. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Bilinear forms through the binary Bell polynomials, N solitons and Backlund transformations of the Boussinesq-Burgers system for the shallow water waves in a lake or near an ocean beach. Commun. Theor. Phys.; 2020; 72, 095002. [DOI: https://dx.doi.org/10.1088/1572-9494/aba23d]

7. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Optical waves/modes in a multicomponent inhomogeneous optical fiber via a three-coupled variable-coefficient nonlinear Schrodinger system. Appl. Math. Lett.; 2021; 120, 107161. [DOI: https://dx.doi.org/10.1016/j.aml.2021.107161]

8. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Symbolic computation on a (2+1)-dimensional generalized variable-coefficient Boiti-Leon-Pempinelli system for the water waves. Chaos Solitons Fractals; 2021; 150, 111066. [DOI: https://dx.doi.org/10.1016/j.chaos.2021.111066]

9. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Looking at an open sea via a generalized (2+1)-dimensional dispersive long-wave system for the shallow water: Scaling transformations, hetero-Backlund transformations, bilinear forms and N solitons. Eur. Phys. J. Plus; 2021; 136, pp. 1-9. [DOI: https://dx.doi.org/10.1140/epjp/s13360-021-01773-6]

10. Gao, X.Y.; Guo, Y.J.; Shan, W.R.; Yin, H.M.; Du, X.X.; Yang, D.Y. Certain electromagnetic waves in a ferromagnetic film. Commun. Nonlinear Sci. Numer. Simul.; 2022; 105, 106066. [DOI: https://dx.doi.org/10.1016/j.cnsns.2021.106066]

11. Saad, K.M.; Atangana, A.; Baleanu, D. New fractional derivatives with non-singular kernel applied to the Burgers equation. Chaos Interdiscip. J. Nonlinear Sci.; 2018; 28, 063109. [DOI: https://dx.doi.org/10.1063/1.5026284]

12. Nonlaopon, K.; Naeem, M.; Zidan, A.; Shah, R.; Alsanad, A.; Gumaei, A. Numerical Investigation of the Time-Fractional Whitham–Broer–Kaup Equation Involving without Singular Kernel Operators. Complexity; 2021; 2021, pp. 1-21. [DOI: https://dx.doi.org/10.1155/2021/7979365]

13. Baleanu, D.; Machado, J.A.T.; Luo, A.C. Fractional Dynamics and Control; Springer: New York, NY, USA, 2012.

14. Poincare, H. Memoires et observations. Sur l’equilibre d’une masse fluide animee d’un mouvement de rotation. Bull. Astron. Ser. I; 1885; 2, pp. 109-118. [DOI: https://dx.doi.org/10.3406/bastr.1885.2592]

15. Adomian, G. Analytical solution of Navier-Stokes flow of a viscous compressible fluid. Found. Phys. Lett.; 1995; 8, pp. 389-400. [DOI: https://dx.doi.org/10.1007/BF02187819]

16. Krasnoschok, M.; Pata, V.; Siryk, S.V.; Vasylyeva, N. A subdiffusive Navier-Stokes-Voigt system. Phys. D Nonlinear Phenom.; 2020; 409, 132503. [DOI: https://dx.doi.org/10.1016/j.physd.2020.132503]

17. Wang, Y.; Zhao, Z.; Li, C.; Chen, Y.Q. Adomian’s method applied to Navier-Stokes equation with a fractional order. Proceedings of the ASME 2009 IDETC/CIE; San Diego, CA, USA, 30 August–2 September 2009; pp. 1047-1054.

18. Yu, Q.; Song, J.; Liu, F.; Anh, V.; Turner, I. An approximate solution for the Rayleigh-Stokes problem for a heated generalized second grade fluid with fractional derivative model using the Adomian decomposition method. J. Algorithms Comput. Technol.; 2009; 3, pp. 553-571. [DOI: https://dx.doi.org/10.1260/174830109789621310]

19. Krasnoschok, M.; Pata, V.; Siryk, S.V.; Vasylyeva, N. Equivalent definitions of Caputo derivatives and applications to subdiffusion equations. Dyn. PDE; 2020; 17, pp. 383-402. [DOI: https://dx.doi.org/10.4310/DPDE.2020.v17.n4.a4]

20. Bazhlekova, E.; Jin, B.; Lazarov, R.; Zhou, Z. An analysis of the Rayleigh-Stokes problem for a generalized second-grade fluid. Numer. Math.; 2015; 131, pp. 1-31. [DOI: https://dx.doi.org/10.1007/s00211-014-0685-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28615736]

21. El-Shahed, M.; Salem, A. On the generalized Navier-Stokes equations. Appl. Math. Comput.; 2004; 156, pp. 287-293. [DOI: https://dx.doi.org/10.1016/j.amc.2003.07.022]

22. Kumar, D.; Singh, J.; Kumar, S. A fractional model of Navier-Stokes equation arising in unsteady flow of a viscous fluid. J. Assoc. Arab. Univ. Basic Appl. Sci.; 2015; 17, pp. 14-19. [DOI: https://dx.doi.org/10.1016/j.jaubas.2014.01.001]

23. Ganji, Z.Z.; Ganji, D.D.; Ganji, A.D.; Rostamian, M. Analytical solution of time-fractional Navier-Stokes equation in polar coordinate by homotopy perturbation method. Numer. Methods Partial. Differ. Equ. Int. J.; 2010; 26, pp. 117-124. [DOI: https://dx.doi.org/10.1002/num.20420]

24. Ragab, A.A.; Hemida, K.M.; Mohamed, M.S.; Abd El Salam, M.A. Solution of time-fractional Navier-Stokes equation by using homotopy analysis method. Gen. Math. Notes; 2012; 13, pp. 13-21.

25. Maitama, S. Analytical solution of time-fractional Navier-Stokes equation by natural homotopy perturbation method. Prog. Fract. Differ. Appl.; 2018; 4, pp. 123-131. [DOI: https://dx.doi.org/10.18576/pfda/040206]

26. Birajdar, G.A. Numerical solution of time fractional Navier-Stokes equation by discrete Adomian decomposition method. Nonlinear Eng.; 2014; 3, pp. 21-26. [DOI: https://dx.doi.org/10.1515/nleng-2012-0004]

27. Kumar, S.; Kumar, D.; Abbasbandy, S.; Rashidi, M.M. Analytical solution of fractional Navier-Stokes equation by using modified Laplace decomposition method. Ain Shams Eng. J.; 2014; 5, pp. 569-574. [DOI: https://dx.doi.org/10.1016/j.asej.2013.11.004]

28. Chaurasia, V.B.L.; Kumar, D. Solution of the time-fractional Navier-Stokes equation. Gen. Math. Notes; 2011; 4, pp. 49-59.

29. Prakash, A.; GAlaoui, M.; Fayyaz, R.; Khan, A.; Shah, R.; Abdo, M. Analytical Investigation of Noyes–Field Model for Time-Fractional Belousov–Zhabotinsky Reaction. Complexity; 2021; 2021, pp. 1-21.

30. El-Tawil, M.A.; Huseen, S.N. On convergence of the q-homotopy analysis method. Int. J. Contemp. Math. Sci.; 2013; 8, pp. 481-497. [DOI: https://dx.doi.org/10.12988/ijcms.2013.13048]

31. Prakash, A.; Veeresha, P.; Prakasha, D.G.; Goyal, M. A new efficient technique for solving fractional coupled Navier-Stokes equations using q-homotopy analysis transform method. Pramana; 2019; 93, pp. 1-10. [DOI: https://dx.doi.org/10.1007/s12043-019-1763-x]

32. Arafa, A.A.; Hagag, A.M.S. Q-homotopy analysis transform method applied to fractional Kundu-Eckhaus equation and fractional massive Thirring model arising in quantum field theory. Asian-Eur. J. Math.; 2019; 12, 1950045. [DOI: https://dx.doi.org/10.1142/S1793557119500451]

33. Jena, R.M.; Chakraverty, S. Q-Homotopy Analysis Aboodh Transform Method based solution of proportional delay time-fractional partial differential equations. J. Interdiscip. Math.; 2019; 22, pp. 931-950. [DOI: https://dx.doi.org/10.1080/09720502.2019.1694742]

34. Iqbal, N.; Akgül, A.; Shah, R.; Bariq, A.; Mossa Al-Sawalha, M.; Ali, A. On Solutions of Fractional-Order Gas Dynamics Equation by Effective Techniques. J. Funct. Spaces; 2022; 2022, pp. 1-14. [DOI: https://dx.doi.org/10.1155/2022/3341754]

35. Khuri, S.A. A Laplace decomposition algorithm applied to a class of nonlinear differential equations. J. Appl. Math.; 2001; 1, pp. 141-155. [DOI: https://dx.doi.org/10.1155/S1110757X01000183]

36. Khuri, S.A. A new approach to Bratus problem. Appl. Math. Comput.; 2004; 147, pp. 131-136.

37. Mohamed, M.A.; Torky, M.S. Numerical solution of nonlinear system of partial differential equations by the Laplace decomposition method and the Pade approximation. Am. J. Comput. Math.; 2013; 3, 175. [DOI: https://dx.doi.org/10.4236/ajcm.2013.33026]

38. Ghazi, F.F.; Tawfiq, L.N.M. Coupled Laplace-decomposition method for solving Klein-Gordon equation. Int. J. Mod. Math. Sci.; 2020; 18, pp. 31-41.

39. Hosseinzadeh, H.; Jafari, H.; Roohani, M. Application of Laplace decomposition method for solving Klein-Gordon equation. World Appl. Sci. J.; 2010; 8, pp. 809-813.

40. Khan, H.; Farooq, U.; Shah, R.; Baleanu, D.; Kumam, P.; Arif, M. Analytical Solutions of (2+Time Fractional Order) Dimensional Physical Models, Using Modified Decomposition Method. Appl. Sci.; 2019; 10, 122. [DOI: https://dx.doi.org/10.3390/app10010122]

41. Yang, X.J. A new integral transform method for solving steady heat-transfer problem. Therm. Sci.; 2016; 20, (Suppl. 3), pp. 639-642. [DOI: https://dx.doi.org/10.2298/TSCI16S3639Y]

42. Caputo, M.; Fabrizio, M. On the singular kernels for fractional derivatives. some applications to partial differential equations. Progr. Fract. Differ. Appl.; 2021; 7, pp. 1-4.

43. Shah, R.; Khan, H.; Baleanu, D. Fractional Whitham–Broer–Kaup Equations within Modified Analytical Approaches. Axioms; 2019; 8, 125. [DOI: https://dx.doi.org/10.3390/axioms8040125]