New Numerical Results for the Time-Fractional Phi-Four Equation Using a Novel Analytical Approach

Abstract

This manuscript investigates the fractional Phi-four equation by using $q$-homotopy analysis transform method ($q$-HATM) numerically. The Phi-four equation is obtained from one of the special cases of the Klein-Gordon model. Moreover, it is used to model the kink and anti-kink solitary wave interactions arising in nuclear particle physics and biological structures for the last several decades. The proposed technique is composed of Laplace transform and $q$-homotopy analysis techniques, and fractional derivative defined in the sense of Caputo. For the governing fractional-order model, the Banach’s fixed point hypothesis is studied to establish the existence and uniqueness of the achieved solution. To illustrate and validate the effectiveness of the projected algorithm, we analyze the considered model in terms of arbitrary order with two distinct cases and also introduce corresponding numerical simulation. Moreover, the physical behaviors of the obtained solutions with respect to fractional-order are presented via various simulations.

1. Introduction

Recently, many new fractional-order operator types have been introduced to the literature by some scholars, and established its fundamental properties [1,2,3,4,5,6]. Fractional theory is widely used to explain many properties of phenomena such as nanotechnology [7], optics [8], human diseases [9], chaos theory [10], and others [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]. As an example, the tumor-immune surveillance model has been investigated in [33]. Some important properties of dengue fever have been effectively investigated in [34,35]. The analytical as well as numerical solutions for the equations illustrating the above cited phenomena, have an important role in describing the behavior of nonlinear models ascends [36,37,38,39,40]. For the differential equation, symmetry is a transformation that keeps its family of solutions invariant, and further, its analysis can be applied to examine and illustrate various classes of differential equations. The study of fractional calculus associated to symmetry recently attracted many researchers from different disciplines in order to present their viewpoints while analyzing real-word problems.

The Klein-Gordon (KG) equation is derived by physicists Klein and Gordon to investigate the behavior of relativistic electrons. As a specific form of the KG equation, the Phi-four equation is defined to illustrate wave interactions [41,42] which is considered as in fractional order, namely, fractional Phi-four (FPF) equation

$$D^{\mu }v\left(x,t\right)=v_{xx}\left(x,t\right)+{\lambda }_{1}v\left(x,t\right)+{\lambda }_2{v}^3\left(x,t\right), 1<\mu \le 2,$$

(1)

with

$$v\left(x,0\right)=f_1\left(x\right) \mathrm{and}{v}_t\left(x,0\right)=f_2\left(x\right).$$

Here, ${\lambda }_1$ and ${\lambda }_2$ are real-valued constants. Here,$v\left(x, t\right)$ is the normalized propagation of distance and retracted time $t$.

Many physicists and mathematicians recently proposed very accurate and more effective algorithms to find and examine the solution for equations describing nonlinear mechanisms arising in science and engineering. As an approximate method, the homotopy analysis method (HAM) has been introduced by L. Shijun [43,44], which is formed to describe deformation from zero to 1. Recently, it has been efficiently and beneficially employed to examine the nature of nonlinear problems without linearization and perturbation. However, HAM necessitates huge time and computer memory for computational work. Therefore, there is an essence of the consolidation of this scheme with familiar transform algorithms. The projected solution procedure is the mixture of $q$-HAM with LT [45]. Since the future scheme is a modified method of HAM, it does not kernel discretization, perturbation or linearization [46,47,48,49,50,51,52,53,54,55,56].

The considered nonlinear problem recently fascinated authors from various areas of science. Due to the numerous applications of the proposed model and the important role it plays in describing various nonlinear phenomena, several researchers found and investigated the solution numerically and also analytically; for instance, researchers in [57] find the new exact solution for the fractional case of considered nonlinear problem, and the spectral collocation method is considered in [58] to find the solution for PF equation. The considered nonlinear PF equation has attracted the attention of researchers, and, in order to present their viewpoint, many techniques are considered [59,60,61,62,63,64,65], presenting some interesting results. Recently many researchers from all over the world have investigated Fractional Calculus (FC) as an effective tool by comparing with other operators some important physical problems [66,67,68,69,70,71,72,73,74,75,76,77,78,79]. For instance, authors in [80] presented the reflection symmetry in fractional calculus, and proposed some interesting properties and corresponding consequences.

The main aim of this paper is to investigate the numerical solutions for the considered fractional Phi-Four equation using the novel technique. The proposed algorithm provides more liberty to consider the diverse class complex as well as nonlinear problems and the initial guess.

In this paper, we try to capture in the physical behaviors obtained numerical solutions with distinct arbitrary order and the parameters accessible by the projected algorithm. Moreover, in order to present accuracy and effectiveness, we present the numerical simulations. By the help of these results, we try to explain the diverse model exemplifying numerous phenomena arising in daily life. The rest of the paper is organized as follows: the basic properties of FC are presented in Section 2, the basic solution algorithm of the projected method for the considered problem is defined in Section 3, the convergence analysis of the considered scheme is illustrated in Section 4, the solution for FPF equation and its corresponding consequences in terms plots and also numerical simulation are respectively cited in Section 5 and Section 6. Finally, some important findings of the paper are presented as a conclusion section.

2. Preliminaries

In this subsection of the paper, it is presented the basic definitions of fractional calculus.

Definition 1.

Fractional integral of$v\left(x,t\right)\in C_{\delta }\left(\delta \ge -1\right)$with the order$\mu$is defined in RL as

$$\begin{gathered} J^{\mu }v\left(x,t\right)=\frac{1}{\Gamma \left(\mu \right)}{{\displaystyle \int }}_0^t{\left(t-\vartheta \right)}^{\mu -1}v\left(x,\vartheta \right)d\vartheta , \\ {J}^{0}v\left(x,t\right)=v\left(x,t\right). \end{gathered}$$

(2)

Definition 2.

Fractional derivative of$f\in C_{-1}^n$is given in the sense of Caputo as the following:

$$D_t^{\mu }v\left(x,t\right)=\{\begin{array}{ll}\frac{d^{n}v\left(x,t\right)}{d{t}^n}, & \mu =n\in \mathbb{N},\\ \frac{1}{\Gamma \left(n-\mu \right)}{{\displaystyle \int }}_0^t{\left(x, t-\vartheta \right)}^{n-\mu -1}{v}^{\left(n\right)}\left(x,\vartheta \right)d\vartheta ,& n-1<\mu <n ,n\in \mathbb{N}.\end{array}$$

(3)

Definition 3.

The Laplace transform (LT) of$D_t^{\mu }v\left(x,t\right)$for a Caputo fractional derivative is introduced as below:

$$L\left[D_t^{\mu }v\left(x,t\right)\right]=s^{\mu }V\left(x,s\right)-{{\displaystyle \sum }}_{r=0}^{n-1}{s}^{\mu -r-1}{v}^{\left(r\right)}\left(x,0^{+}\right), \left(n-1<\mu \le n\right),$$

(4)

where$V\left(x, s\right)$is LT of$v\left(x, t\right)$.

3. Solution Procedure for Fractional Phi-Four Equation

Here, we consider the nonlinear FPF equation to illustrate the basic solution algorithm of the projected method with associated initial conditions as follows:

$$D_t^{\mu }v\left(x,t\right)=\frac{{\partial }^{2}v\left(x,t\right)}{\partial x^2}+{\lambda }_{1}v\left(x,t\right)+{\lambda }_2{v}^3\left(x,t\right), 1<\mu \le 2, t>0,$$

(5)

and

$$v\left(x,0\right)=g_1\left(x\right) \mathrm{and} v_t\left(x,0\right)=g_2\left(x\right).$$

(6)

Here, $D_t^{\mu }v\left(x,t\right)$ signifies the fractional Caputo-derivative of $v\left(x,t\right).$ Here, $v\left(x,t\right)$ is a bounded function. On using the LT on Equation (5), and then by the help of Equation (6), it is given as following

$$L\left[v\left(x,t\right)\right]-\frac{1}{s}\left[g_1\left(x\right)\right]-\frac{1}{s^2}\left[g_2\left(x\right)\right]-\frac{1}{s^{\mu }}L\left\{\frac{{\partial }^{2}v}{\partial x^2}+{\lambda }_1 v+{\lambda }_2{v}^3\right\}=0.$$

(7)

By the assist of Equation (7), the nonlinear operator is given as follows:

$$\begin{array}{ll}\mathcal{N}\left[\phi \left(x,t;q\right)\right]& =L\left[\phi \left(x,t;q\right)\right]-\frac{1}{s}\left[g_1\left(x\right)\right]-\frac{1}{s^2}\left[g_2\left(x\right)\right]\\ & -\frac{1}{s^{\mu }}L\left\{\frac{{\partial }^2\phi \left(x,t;q\right)}{\partial x^2}+{\lambda }_1\phi \left(x,t;q\right)+{\lambda }_2{\phi }^3\left(x,t;q\right)\right\},\end{array}$$

(8)

Now, the homotopy structure is defined by

$$\left(1-nq\right)L\left[\phi \left(x, t;q\right)-v_0\left(x,t\right)\right]=\hslash q\mathcal{N}\left[\phi \left(x,t;q\right)\right],$$

(9)

where $q\in \left[0,\frac{1}{n}\right]\left(n\ge 1\right)$ and $\hslash \ne 0$ are respectively the embedding and auxiliary parameter. Further, these parameters help us to control and adjust the convergence region of the obtained solution. For the proper choice of $n$ and $\hslash$, the obtained solution quickly converges to the exact solution in an admissible domain. The following results are respectively true for $q=0$ and $q=\frac{1}{n}$

$$\phi \left(x,t;0\right)=v_0\left(x,t\right), \phi \left(x,t;\frac{1}{n}\right)=v\left(x,t\right).$$

(10)

Near to $q$, we define $\phi \left(x,t;q\right)$ in series form by help of the Taylor theorem and then one can obtain

$$\phi \left(x,t;q\right)=v_0\left(x,t\right)+{{\displaystyle \sum }}_{m=1}^{\infty }{v}_m\left(x,t\right)q^m,$$

(11)

where

$$v_m\left(x,t\right)=\frac{1}{m!}\frac{{\partial }^m\phi \left(x, t;q\right)}{\partial q^m}{|}_{q=0}.$$

(12)

For the proper chaise of $v_0\left(x,t\right),n$ and $\hslash$; the series (11) converges at $q=\frac{1}{n}$. Then

$$v\left(x,t\right)=v_0\left(x,t\right)+{{\displaystyle \sum }}_{m=1}^{\infty }{v}_m\left(x,t\right){\left(\frac{1}{n}\right)}^m.$$

(13)

Dividing by $m!$ after $m$-times differentiating Equation (9) with $q$ and then putting $q=0$, it simplifies to

$$L\left[v_m\left(x,t\right)-k_m{v}_{m-1}\left(x,t\right)\right]=\hslash {\Re }_m\left({\overrightarrow{v}}_{m-1}\right),$$

(14)

where $k_m$ and $v_m$ are defined by

$$k_m=\{\begin{array}{c}0, m\le 1,\\ n, m>1.\end{array}$$

(15)

and

$${\overrightarrow{v}}_m=\left\{v_0\left(x,t\right),v_1\left(x,t\right),\dots ,v_m\left(x,t\right)\right\}.$$

(16)

On applying inverse LT on Equation (14), produces

$$v_m\left(x,t\right)=k_m{v}_{m-1}\left(x,t\right)+\hslash L^{-1}\left[{\Re }_m\left({\overrightarrow{v}}_{m-1}\right)\right],$$

(17)

where

$$\begin{gathered} {\Re }_m\left[{\overrightarrow{v}}_{m-1}\right]=L\left[v_{m-1}\left(x,t\right)\right]-\left(1-\frac{k_m}{n}\right)\left(\frac{1}{s}\left[g_1\left(x\right)\right]+\frac{1}{s^2}\left[g_2\left(x\right)\right]\right) \\ -\frac{1}{s^{\mu }}L\left\{\frac{{\partial }^2{v}_{m-1}}{\partial x^2}+{\lambda }_1{v}_{m-1}+{\lambda }_2{\displaystyle \sum }_{i=0}^{m-1}\left({\displaystyle \sum }_{j=0}^i{v}_j{v}_{i-j}\right)v_{m-1-i}\right\}. \end{gathered}$$

(18)

With the help of Equations (17) and (18), we have

$$\begin{array}{ll}{v}_m\left(x,t\right)& =\left(k_m+\hslash \right)v_{m-1}\left(x,t\right)-\left(1-\frac{k_m}{n}\right)L^{-1}\left\{\frac{1}{s}\left[g_1\left(x\right)\right]+\frac{1}{s^2}\left[g_2\left(x\right)\right]\right\}\\ & -\hslash L^{-1}\left[\frac{1}{s^{\mu }}L\left\{\frac{{\partial }^2{v}_{m-1}}{\partial x^2}+{\lambda }_1{v}_{m-1}+{\lambda }_2{{\displaystyle \sum }}_{i=0}^{m-1}\left({{\displaystyle \sum }}_{j=0}^i{v}_j{v}_{i-j}\right)v_{m-1-i}\right\}\right].\end{array}$$

(19)

Using Equation (19), we achieve the series of $v_m\left(x,t\right)$. Lastly, the series $q$-HATM solution is defined as

$$v\left(x,t\right)={\displaystyle \sum }_{m=0}^{\infty }{v}_m\left(x,t\right).$$

(20)

4. Convergence Analysis of the Proposed Method

In this part of the paper, for the considered problem, we introduce the convergence analysis as follows.

Theorem 1 (Convergence theorem).

Let$F:E\to E$is mapping (nonlinear) with Banach space$E$. Presume that

$$\Vert F\left(u\right)-F\left(v\right)\Vert \le \alpha \Vert u-v\Vert , \forall u, v\in E,$$

(21)

then there is a fixed point for$F$with the help of Banach’s fixed point theory [81,82,83]. Then, for $u_0, v_0\in E$ the sequence obtained by $q$ -HATM solution converges to a specified point of $F$, and

$$\Vert v_m-v_p\Vert \le \frac{{\alpha }^p}{1-\alpha }\Vert v_1-v_0\Vert .$$

(22)

Proof. 

Let us consider a Banach space $\left(C\left[J\right],\Vert .\Vert \right)$ with the norm defined as $\Vert g\left(t\right)\Vert =\underset{t\in J}{\mathit{max}}\left|g\left(t\right)\right|$ for all continuous function on $J$. Now, we verify $\left\{v_p\right\}$ is Cauchy sequence in $\left(C\left[J\right],\Vert .\Vert \right)$. For that, let

$$\begin{array}{l}\Vert v_m-v_p\Vert =\underset{t\in J}{\mathit{max}}\left|v_m-v_p\right|\\ =\underset{t\in J}{\mathit{max}}|\left(k_m+\hslash \right)\left(v_{m-1}-v_{p-1}\right)-\hslash L^{-1}\left(\frac{1}{s^{\mu }}L\left[\left(\frac{{\partial }^2{v}_{m-1}}{\partial x^2}-\frac{{\partial }^2{v}_{p-1}}{\partial x^2}\right)+{\lambda }_1\left(v_{m-1}-v_{p-1}\right)+{\lambda }_2\left(v_{m-1}^3-v_{p-1}^3\right)\right]\right)|\\ \le \underset{t\in J}{\mathit{max}}\{\left(k_m+\hslash \right)\left|\left(v_{m-1}-v_{p-1}\right)\right|-\hslash L^{-1}(\frac{1}{s^{\mu }}L[\left|\frac{{\partial }^2{v}_{m-1}}{\partial x^2}-\frac{{\partial }^2{v}_{p-1}}{\partial x^2}\right|\\ +{\lambda }_1\left|v_{m-1}-v_{p-1}\right|+{\lambda }_2\left|v_{m-1}^3-v_{p-1}^3\right|])\}.\end{array}$$

For LT, with the help of the convolution theorem, we get

$$\begin{array}{l}\Vert v_m-v_p\Vert \le \underset{t\in J}{max}[\left(k_m+\hslash \right)\left|\left(v_{m-1}-v_{p-1}\right)\right|\\ -\hslash {{\displaystyle \int }}_0^t\left(\left|\frac{{\partial }^2{v}_{m-1}}{\partial x^2}-\frac{{\partial }^2{v}_{p-1}}{\partial x^2}\right|+{\lambda }_1\left|v_{m-1}-v_{p-1}\right|+{\lambda }_2\left|v_{m-1}^3-v_{p-1}^3\right|\right)]\frac{{\left(t-\xi \right)}^{\mu }}{\Gamma \left(\mu +1\right)}d\xi \\ \le \underset{t\in J}{max}[\left(k_m+\hslash \right)\left|v_{m-1}-v_{p-1}\right|-\hslash {{\displaystyle \int }}_0^t({\delta }^2\left|v_{m-1}-v_{p-1}\right|+{\lambda }_1\left|\left(v_{m-1}-v_{p-1}\right)\right|\\ +{\lambda }_2\left|\left(v_{m-1}-v_{p-1}\right)\left(P^2+Q^2+PQ\right)\right|)]\frac{{\left(t-\xi \right)}^{\mu }}{\Gamma \left(\mu +1\right)}d\xi .\end{array}$$

Then, the above inequality reduces to

$$\begin{array}{l}\Vert v_m-v_p\Vert \le \underset{t\in J}{max}[\left(k_m+\hslash \right)\left|v_{m-1}-v_{p-1}\right|+\hslash ({\delta }^2\left|v_{m-1}-v_{p-1}\right|\\ +{\lambda }_1\left|\left(v_{m-1}-v_{p-1}\right)\right|+{\lambda }_2\left|\left(v_{m-1}-v_{p-1}\right)\left(P^2+Q^2+PQ\right)\right|)\mathcal{T}],\\ \Vert v_m-v_p\Vert \le \alpha \Vert v_{m-1}-v_{p-1}\Vert .\end{array}$$

Setting $m=p+1$, it yields

$$\Vert v_{p+1}-v_p\Vert \le \alpha \Vert v_p-v_{p-1}\Vert \le {\alpha }^2\Vert v_{p-1}-v_{p-2}\Vert \le \dots \le {\alpha }^p\Vert v_1-v_0\Vert .$$

On using triangular inequality, we have

$$\begin{array}{l}\Vert v_m-v_p\Vert \le \Vert v_{p+1}-v_p\Vert +\Vert v_{p+2}-v_{p+1}\Vert +\cdots +\Vert v_m-v_{m-1}\Vert \\ \le \left[{\alpha }^p+{\alpha }^{p+1}+\dots +{\alpha }^{m-1}\right]\Vert v_1-v_0\Vert \\ \le {\alpha }^p\left[1+\alpha +{\alpha }^2+\dots +{\alpha }^{m-p-1}\right]\Vert v_1-v_0\Vert \\ \le {\alpha }^p\left[\frac{1-{\alpha }^{m-p-1}}{1-\alpha }\right]\Vert v_1-v_0\Vert .\end{array}$$

As $0<\alpha <1$, so $1-{\alpha }^{m-p-1}<1$, then we have

$$\Vert v_m-v_p\Vert \le \frac{{\alpha }^p}{1-\alpha }\Vert v_1-v_0\Vert .$$

But $\Vert v_1-v_0\Vert <\infty$, consequently as $m\to \infty$ than $\Vert v_m-v_p\Vert \to 0$, which proves $\left\{v_p\right\}$ is Cauchy sequence in $C\left[J\right]$. It completes our required proof. □

Theorem 2 (Uniqueness theorem).

For Equation (5), the$q$-HATM solution is unique wherever$0<\alpha <1$, where$\alpha =\left(k_m+\hslash \right)+\hslash \left({\delta }^2+{\lambda }_1+{\lambda }_2\left(P^2+Q^2+PQ\right)\right)\mathcal{T}$.

Proof. 

For fractional PF equation cited in Equation (5), the solution is illustrated as following

$$v\left(x,t\right)={\sum }_{m=0}^{\infty }{v}_m\left(x,t\right)$$

where

$$\begin{gathered} v_m\left(x,t\right)=\left(k_m+\hslash \right)v_{m-1}\left(x,t\right)-\left(1-\frac{k_m}{n}\right)L^{-1}\left\{\frac{1}{s}\left[f_1\left(x\right)\right]+\frac{1}{s^2}\left[f_2\left(x\right)\right]\right\} \\ -\hslash L^{-1}\left[\frac{1}{s^{\mu }}L\left\{\frac{{\partial }^2{v}_{m-1}}{\partial x^2}+{\lambda }_1{v}_{m-1}+{\lambda }_2{{\displaystyle \sum }}_{i=0}^{m-1}\left({{\displaystyle \sum }}_{j=0}^i{v}_j{v}_{i-j}\right)v_{m-1-i}\right\}\right]. \end{gathered}$$

Suppose $v$ and $v^{\star }$ be the two different solutions for Equation (5) such that $\left|v\right|\le P$ and $\left|v^{\star }\right|\le Q$, then using the above relation, it is obtained as following

$$\left|v-v^{\star }\right|=\left|\left(k_m+\hslash \right)\left(v-v^{\star }\right)-\hslash L^{-1}\left\{\frac{1}{s^{\mu }}L\left(\left(\frac{{\partial }^{2}v}{\partial x^2}-\frac{{\partial }^2{v}^{\star }}{\partial x^2}\right)+{\lambda }_1\left(v-v^{\star }\right)+{\lambda }_2\left(v^3-v^{\star }{}^3\right)\right)\right\}\right|.$$

(23)

For Laplace transform, we have, by the help of the convolution theorem,

$$\begin{gathered} \left|v-v^{\star }\right|=\left(k_m+\hslash \right)\left|v-v^{\star }\right|+\hslash {{\displaystyle \int }}_0^t\left(\left|\frac{{\partial }^{2}v}{\partial x^2}-\frac{{\partial }^2{v}^{\star }}{\partial x^2}\right|+\left|{\lambda }_1\left(v-v^{\star }\right)\right|+\left|{\lambda }_2\left(v^3-v^{\star }{}^3\right)\right|\right)\frac{{\left(t-\xi \right)}^{\mu }}{\Gamma \left(\mu +1\right)}d\xi , \\ \le \left(k_m+\hslash \right)\left|v-v^{\star }\right|+\hslash {{\displaystyle \int }}_0^t(\frac{{\partial }^2}{\partial x^2}\left|v-v^{\star }\right|+{\lambda }_1\left|\left(v-v^{\star }\right)\right| \\ +{\lambda }_2\left|\left(v-v^{\star }\right)\left(v^2+v^{\star }{}^2+v{v}^{\star }\right)\right|\frac{{\left(t-\xi \right)}^{\mu }}{\Gamma \left(\mu +1\right)}d\xi , \\ \le \left(k_m+\hslash \right)\left|v-v^{\star }\right|+\hslash {{\displaystyle \int }}_0^t({\delta }^2\left|v-v^{\star }\right|+{\lambda }_1\left|\left(v-v^{\star }\right)\right| \\ +{\lambda }_2\left|\left(v-v^{\star }\right)\left(P^2+Q^2+PQ\right)\right|\frac{{\left(t-\xi \right)}^{\mu }}{\Gamma \left(\mu +1\right)}d\xi , \end{gathered}$$

where ${\delta }^2=\frac{{\partial }^2}{\partial x^2}$. Using integral mean value, we have

$$\begin{gathered} \left|v-v^{\star }\right|\le \left(k_m+\hslash \right)\left|v-v^{\star }\right|+\hslash ({\delta }^2\left|v-v^{\star }\right|+{\lambda }_1\left|\left(v-v^{\star }\right)\right| \\ +{\lambda }_2\left|\left(v-v^{\star }\right)\left(P^2+Q^2+PQ\right)\right|)\mathcal{T} \\ \left|v-v^{\star }\right|\le \alpha \left|v-v^{\star }\right| \\ \left(1-\alpha \right)\left|v-v^{\star }\right|\le 0. \end{gathered}$$

Since $0<\alpha <1$, therefore $\left|v-v^{\star }\right|=0$, which gives $v=v^{\star }$. This completes the required proof. □

5. Solution for Fractional PF Equation

In order to validate the applicability and efficiency of the future technique, here we consider the two distinct examples.

Example 1.

Consider the time-fractional Phi-four equation [84,85,86]

$$D_t^{\mu }v\left(x,t\right)=\frac{{\partial }^{2}v\left(x,t\right)}{\partial x^2}-{\lambda }_{1}v\left(x,t\right)-{\lambda }_2{v}^3\left(x,t\right),$$

(24)

subjected to the initial conditions

$$v\left(x, 0\right)=\sqrt{\frac{-{\lambda }_1{}^2}{{\lambda }_2}}\tan h\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}x\right)\mathrm{and}{v}_t\left(x, 0\right)=-{\lambda }_1\eta \sqrt{\frac{-{\lambda }_1{}^2}{2{\lambda }_2\left({\eta }^2-1\right)}}{\sec h}^2\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}x\right).$$

(25)

By using the proposed algorithm and Equation (25), the Equation (24) becomes

$$\begin{gathered} v_m\left(x,t\right)=\left(k_m+\hslash \right)v_{m-1}\left(x,t\right)-\left(1-\frac{k_m}{n}\right)L^{-1}\{\frac{1}{s}\left(\sqrt{\frac{-{\lambda }_1{}^2}{{\lambda }_2}}\tan h\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}x\right)\right) \\ -\frac{1}{s^2}\left({\lambda }_1\eta \sqrt{\frac{-{\lambda }_1{}^2}{2{\lambda }_2\left({\eta }^2-1\right)}}{\sec h}^2\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}x\right)\right)\}-\hslash L^{-1}[\frac{1}{s^{\mu }}L\{\frac{{\partial }^2{v}_{m-1}}{\partial x^2}-{\lambda }_{1}v \\ -{\lambda }_2{{\displaystyle \sum }}_{i=0}^{m-1}\left({{\displaystyle \sum }}_{j=0}^i{v}_j{v}_{i-j}\right)v_{m-1-i}\}]. \end{gathered}$$

(26)

We get the iterative terms of $v_m\left(x,t\right)$ by solving Equation (26), and as follows

$$\begin{array}{c}{v}_0\left(x,t\right)=\sqrt{\frac{-{\lambda }_1{}^2}{{\lambda }_2}}\tan h\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}x\right)-{\lambda }_1\eta t\sqrt{\frac{-{\lambda }_1{}^2}{2{\lambda }_2\left({\eta }^2-1\right)}}{\sec h}^2\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}x\right)\\ v_1\left(x,t\right)=\frac{{\lambda }_1{}^2{\eta }^2\hslash t^{\mu }}{8\left(-1+{\eta }^2\right)\Gamma \left[\mu +1\right]}{\sec h}^6\left(\frac{{\lambda }_{1}x\sqrt{\frac{1}{-1+{\eta }^2}}}{\sqrt{2}}\right)(-3\sqrt{2}{\lambda }_{1}t\eta \sqrt{\frac{{\lambda }_1{}^2}{{\lambda }_2-{\eta }^2{\lambda }_2}}\\ +2\sqrt{2}{\lambda }_1{}^3{t}^3\eta \sqrt{\frac{{\lambda }_1{}^2}{{\lambda }_2-{\eta }^2{\lambda }_2}}\\ -2\sqrt{2}{\lambda }_{1}t\eta \sqrt{\frac{{\lambda }_1{}^2}{{\lambda }_2-{\eta }^2{\lambda }_2}}\cos h\left(\sqrt{2}{\lambda }_{1}x\sqrt{\frac{1}{-1+{\eta }^2}}\right)\\ +\sqrt{2}{\lambda }_{1}t\eta \sqrt{\frac{{\lambda }_1{}^2}{{\lambda }_2-{\eta }^2{\lambda }_2}}\cos h\left(2\sqrt{2}{\lambda }_{1}x\sqrt{\frac{1}{-1+{\eta }^2}}\right)\\ +2\sqrt{-\frac{{\lambda }_1{}^2}{{\lambda }_2}}\sin h\left(\sqrt{2}{\lambda }_{1}x\sqrt{\frac{1}{-1+{\eta }^2}}\right)+6t^2{\left(-\frac{{\lambda }_1{}^2}{{\lambda }_2}\right)}^{3/2}{\lambda }_2\sinh \left(\sqrt{2}{\lambda }_{1}x\sqrt{\frac{1}{-1+{\eta }^2}}\right)\\ +\sqrt{-\frac{{\lambda }_1{}^2}{{\lambda }_2}}\sinh \left(2\sqrt{2}{\lambda }_{1}x\sqrt{\frac{1}{-1+{\eta }^2}}\right)),\end{array}$$

Similarly, the remaining part of the series can be calculated (See Figure 1, Figure 2 and Figure 3 and

Table 1. Numerical simulation conducted for FPF equation defined in Example1 with different $x$ and $t$ at ${\lambda }_1=1, {\lambda }_2=-1, \eta =3, n=1, \mu =2$ and $\hslash =-1$.

$\left|v_{Exact}-v_{q-HATM}^{\left(3\right)}\right|$
$x/t$	$0.01$	$0.02$	$0.03$	$0.04$	$0.05$
$-5$	$9.18570\times {10}^{-8}$	$7.38629\times {10}^{-7}$	$2.50566\times {10}^{-6}$	$5.96979\times {10}^{-6}$	$1.17195\times {10}^{-5}$
$-3$	$3.68683\times {10}^{-8}$	$3.07718\times {10}^{-7}$	$1.08179 \times {10}^{-6}$	$2.66706\times {10}^{-6}$	$5.41063\times {10}^{-6}$
$-1$	$2.15263\times {10}^{-7}$	$1.70951 \times {10}^{-6}$	$5.72612 \times {10}^{-6}$	$1.34678\times {10}^{-5}$	$2.60945\times {10}^{-5}$
$1$	$2.18306\times {10}^{-7}$	$1.75819 \times {10}^{-6}$	$5.97254 \times {10}^{-6}$	$1.42464\times {10}^{-5}$	$2.79946\times {10}^{-5}$
$3$	$3.36893\times {10}^{-8}$	$2.56849\times {10}^{-7}$	$8.24228\times {10}^{-7}$	$1.85290\times {10}^{-6}$	$3.42245\times {10}^{-6}$
$5$	$9.09190\times {10}^{-8}$	$7.23617\times {10}^{-7}$	$2.42963\times {10}^{-6}$		$1.11322\times {10}^{-5}$

). Now, the series form of the $q$-HATM solution is as follows:

$$v\left(x,t\right)=v_0\left(x, t\right)+{\displaystyle \sum }_{m=1}^{\infty }{v}_m\left(x,t\right).$$

For $n=1, \mu =2$ and $\hslash =-1$, the $q$-HATM solution converges to analytical solution

$$v\left(x,t\right)=\sqrt{\frac{-{\lambda }_1{}^2}{{\lambda }_2}}\tan h\left({\lambda }_1\sqrt{\frac{1}{2\left({\eta }^2-1\right)}}\left(x-\eta t\right)\right).$$

Example 2.

Consider the FPF equation at${\lambda }_1=-1$ and ${\lambda }_2=1$ [57]

$$D_t^{\mu }v\left(x,t\right)=\frac{{\partial }^{2}v\left(x,t\right)}{\partial x^2}+v\left(x,t\right)-v^3\left(x,t\right),$$

(27)

subjected to the initial conditions

$$v\left(x, 0\right)=\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)\mathrm{and}{v}_t\left(x, 0\right)=\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right).$$

(28)

By using the proposed algorithm and Equation (28), the Equation (27) becomes

$$v_m\left(x,t\right)=\left(k_m+\hslash \right)v_{m-1}\left(x,t\right)-\left(1-\frac{k_m}{n}\right)L^{-1}\left\{\left(\frac{1}{s}+\frac{1}{s^2}\right)\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)\right\}-\hslash L^{-1}\left[\frac{1}{s^{\mu }}L\left\{\frac{{\partial }^2{v}_{m-1}}{\partial x^2}+v-{\displaystyle \sum }_{i=0}^{m-1}\left({\displaystyle \sum }_{j=0}^i{v}_j{v}_{i-j}\right)v_{m-1-i}\right\}\right].$$

(29)

We evaluate the iterative terms of $v_m\left(x,t\right)$ with the assistance of Equation (29), and as follows

$$\begin{gathered} v_0\left(x,t\right)=\left(1+t\right)\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right), \\ {v}_1\left(x,t\right)=\frac{\hslash t^{\mu }}{\Gamma \left[\mu +1\right]}\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)\left(-1-t-\frac{\left(1+t\right){\sec h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)}{k^2-1}+{\left(1+t\right)}^3{\tan h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)\right), \\ {v}_2\left(x,t\right)=\frac{\hslash \left(n+\hslash \right)t^{\mu }}{\Gamma \left[\mu +1\right]}\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)\left(-1-t-\frac{\left(1+t\right){\sec h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)}{k^2-1}+{\left(1+t\right)}^3{\tan h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)\right) \\ +\frac{{\hslash }^2{t}^{2\mu }}{\Gamma \left[2\mu +1\right]}\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)(1+t+\frac{\left(1+t\right)}{-1+k^2}{\sec h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)+\frac{2\left(1+t\right)\left(k^2{\left(1+t\right)}^2-t\left(2+t\right)\right)}{{\left(-1+k^2\right)}^2} \\ \times {\sec h}^4\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)-\frac{\left(1+t\right)\left(k^2{\left(1+t\right)}^2-t\left(2+t\right)\right)\left(-2+\cos h\left(\sqrt{2}\sqrt{\frac{1}{1-k^2}}x\right)\right)}{{\left(-1+k^2\right)}^2}{\sec h}^4\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right) \\ -{\left(1+t\right)}^3{\tan h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)+\frac{2{\sec h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)}{2-2k^2}(-1-t-\frac{\left(1+t\right)\mathrm{Sec}h{\left[\sqrt{\frac{1}{2-2k^2}}x\right]}^2}{-1+k^2} \\ +{\left(1+t\right)}^3{\tan h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right))+3{\left(1+t\right)}^2{\tan h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)(-1-t \\ -\frac{\left(1+t\right){\sec h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right)}{-1+k^2}+{\left(1+t\right)}^3{\tan h}^2\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}x\right))), \end{gathered}$$

Now, the series form of $q$-HATM solution is as follows:

$$v\left(x,t\right)=v_0\left(x, t\right)+{\displaystyle \sum }_{m=1}^{\infty }{v}_m\left(x,t\right).$$

For $n=1, \mu =2$ and $\hslash =-1$, the $q$-HATM solution converges to analytical solution $v\left(x,t\right)=\tan h\left(\sqrt{\frac{1}{2\left(1-k^2\right)}}\left(x-kt\right)\right)$ (See Figure 4, Figure 5 and Figure 6 and

Table 2. Numerical simulation conducted for fractional Phi-four (FPF) equation defined in Example2 with different $x$ and $t$ at $k=0.1, \mu =2, n=1$ and $\hslash =-1$.

$\left|v_{Exact}-v_{q-HATM}^{\left(3\right)}\right|$
$x/t$	$0.01$	$0.02$	$0.03$	$0.04$	$0.05$
$-5$	$1.96030\times {10}^{-5}$	$9.96170\times {10}^{-3}$	$1.99430\times {10}^{-2}$	$2.99243\times {10}^{-2}$	$3.99056\times {10}^{-2}$
$-3$	$3.17168\times {10}^{-4}$	$9.36662\times {10}^{-3}$	$1.90505\times {10}^{-2}$	$2.87343\times {10}^{-2}$	$3.84183\times {10}^{-2}$
$-1$	$4.33502\times {10}^{-3}$	$1.33124\times {10}^{-3}$	$6.99789\times {10}^{-3}$	$1.26649\times {10}^{-2}$	$1.83324\times {10}^{-2}$
$1$	$1.65573\times {10}^{-2}$	$2.31142\times {10}^{-2}$	$2.96714\times {10}^{-2}$	$3.62290\times {10}^{-2}$	$4.27869\times {10}^{-2}$
$3$	$1.97624\times {10}^{-2}$	$2.95239\times {10}^{-2}$	$3.92856\times {10}^{-2}$	$4.90472\times {10}^{-2}$	$5.88089\times {10}^{-2}$
$5$	$1.99868\times {10}^{-2}$	$2.99728\times {10}^{-2}$	$3.99588\times {10}^{-2}$	$4.99447\times {10}^{-2}$	$5.99307\times {10}^{-2}$

).

6. Numerical Results and Discussion

Here, we have demonstrated the numerical simulations of the Phi-four equation having fractional-order by using $q$-HATM. The fourth order solution (i.e., up to $v_4\left(x,t\right)$) we consider to illustrate the behavior of the obtained solution in terms of plots and tables. In Table 1 and Table 2, the numerical studies have been conducted to ensure the exactness of the projected algorithm, and further from the tables, it is confirmed that the considered model conspicuously depends on the time. The surfaces of the $q$-HATM solution, exact solution and absolute error are captured in Figure 1 for the FPF equation considered in Example 1. To show the projected problem remarkably depends on arbitrary order and the response of $q$-HATM solution for different $\mu$, the plots have been drowned and are cited in Figure 2 and Figure 5 for Examples 1 and 2., respectively. Similarly, nature obtained and exact solutions in association with absolute error for Example 2, are presented in Figure 4. The $\hslash$-curves have been drowned to analyze the behavior of the achieved solution related homotopy parameter $\hslash$ with different $\mu$ for both examples and have been respectively captured in Figure 3 and Figure 6. These curves can help us to adjust and regulate the region of convergence for the $q$-HATM solution. For unsuitable $\hslash$, the $q$-HATM solution swiftly converges to exact solution. The demonstrated plots help us to better understand the nature of the fractional Phi-four equation when temporal-spatial variables vary in comparison with arbitrary order.

7. Conclusions

In the present framework, $q$-HATM has been employed to find the numerical solutions for the fractional Phi-four equation. With the help of Banach’s fixed point theory, the convergence analysis projecting the nonlinear problem has been presented. To present the efficiency as well as the applicability of the projected algorithm, we have considered two distinct cases. The present study confirms that the projected nonlinear problem is remarkable with the time instant as well as the time history and these can be effectively exemplified by employing the concept of fractional calculus. The governing model plays a vibrant role while analyzing many physical phenomena, and thus, for future work, this can be examined by using recently proposed and nurtured effective and accurate methods [87,88]. With the help of the obtained results, we can capture more interesting consequences. Finally, it has been observed that the considered scheme is highly effective, more accurate and extremely methodical, and this can be employed to exemplify the various classes of nonlinear models that exist in science and technology. The results obtained demonstrate that the considered method is more effective and easy to employ to scrutinize the behaviors of fractional differential equations with multi-dimensions arising in associated areas of science and technology.