Approximate Solution of Nonlinear Time-Fractional Klein-Gordon Equations Using Yang Transform

Abstract

The algebras of the symmetry operators for the Klein–Gordon equation are important for a charged test particle, moving in an external electromagnetic field in a space time manifold on the isotropic hydrosulphate. In this paper, we develop an analytical and numerical approach for providing the solution to a class of linear and nonlinear fractional Klein–Gordon equations arising in classical relativistic and quantum mechanics. We study the Yang homotopy perturbation transform method ($\mathbb{Y}$HPTM), which is associated with the Yang transform ($\mathbb{Y}$T) and the homotopy perturbation method (HPM), where the fractional derivative is taken in a Caputo–Fabrizio (CF) sense. This technique provides the solution very accurately and efficiently in the form of a series with easily computable coefficients. The behavior of the approximate series solution for different fractional-order ℘ values has been shown graphically. Our numerical investigations indicate that $\mathbb{Y}$HPTM is a simple and powerful mathematical tool to deal with the complexity of such problems.

1. Introduction

Recently, fractional calculus has grown in popularity due to its significant prospective applications in physics and engineering such as biology, mathematics, chemistry, fluid mechanics, physics, and nonlinear optics [1,2]. Fractional partial differential Equations (FPDEs) are a contemporary tool in calculus that can be used to simulate a wide range of classifications in applied sciences and engineering [3,4,5].

The Klein–Gordon (KG) equation performs a significant role in mathematical physics and many other scientific studies such as quantum field theory, nonlinear optics, and solidstate physic [6,7,8,9,10]. On the other hand, the fractional-order KG equation is derived from the classical KG equation by substituting the time order derivative with the fractional derivative of order ℘. The fractional-order KG equation can be illustrated as below

$$\begin{array}{c}\hfill D_q^{\wp }\vartheta (ϵ,q)-D_{ϵ}^2\vartheta (ϵ,q)+a_1\vartheta (ϵ,q)+a_{2}G\left(\vartheta (ϵ,q)\right)=f(ϵ,q),\end{array}$$

(1)

with initial conditions

$$\begin{array}{c}\hfill \vartheta (ϵ,0)=f_1\left(ϵ\right),{\vartheta }_q(ϵ,0)=f_2\left(ϵ\right),\end{array}$$

(2)

where $D_q^{\wp }$ represents the Caputo fractional time derivative, $a_1$ and $a_2$ are real constants, $f(ϵ,q)$, $f_1\left(ϵ\right)$ and $f_2\left(ϵ\right)$ are known as analytical functions, whereas $G\left(\vartheta \right(ϵ,q\left)\right)$ is a nonlinear, and $\vartheta$ is an unknown function of $ϵ$ and q.

Various authors [11,12,13,14,15] have investigated different analytical and numerical strategies to examine the solution to the KG equation but with some restrictions and lacks. Tamsir and Srivastava [16] used fractional reduced differential transform to obtain the analytical solution of linear and nonlinear KG equation with time-fractional order. Bansu and Kumar [17] used a radial basis approach, and Kurulay [18] applied the homotopy analysis method to evaluate the numerical solution of the space-time fractional KG equation. Later, Khader and Adel [19] applied a hybridization scheme to achieve the solution of the fractional KG equation. Zhmud and Dimitrov [20] developed the fractional integration method, which is based on extrapolation using a series of integrating and differentiating links with a time constant that changes symmetrically from one step to another. In order to obtain the solution of FPDEs, several valuable strategies have been considered, such as the generalized differential transform method [21], the adomian decomposition method [22], the homotopy analysis method [23], the variational iteration method [24], the homotopy perturbation method [25], the Elzaki transform decomposition method [26], the fractional wavelet method [27,28] and the residual power series method [29,30].

In this paper, we present the Yang homotopy perturbation transform method ($\mathbb{Y}$HPTM), which is a composition of $\mathbb{Y}$T and HPM. The primary objective of this approach is to investigate the approximate solution of fractional KG equations and minimize the computational work that overcomes nonlinear problems easily. Next, this scheme can promptly deal with the nonlinear KG equation. Finally, this method can reduce the range of the computations and generate an approximate solution with elegantly computed expressions, which is its most impressive advantage. The design of this paper is framed as follows. In Section 2, we start with some primary definitions of Caputo–Fabrizio. In Section 3, we formulate the idea of the Yang homotopy perturbation transform method. In Section 4, we perform this scheme on some illustrative examples to show its capability and efficiency. Concluding remarks are given in Section 5.

2. Preliminaries and Concepts

Definition 1.

The CF derivative is described as [31]

$$\begin{array}{c}\hfill {}^{CF}{D}_q^{\wp }\vartheta (ϵ,q)=\frac{S(\wp )}{1-\wp }{\int }_0^q\left[Q^{\prime }\left(\varrho \right)K(q,\varrho )\right]d\varrho ,n-1<\wp \le n\end{array}$$

(3)

$S(\wp )$ is the normalization function with $S\left(0\right)=S\left(1\right)=1$, and then, Equation (3) becomes as

$$\begin{array}{c}\hfill {}^{CF}{D}_q^{\wp }\vartheta (ϵ,q)=\frac{S(\wp )}{1-\wp }{\int }_0^q[Q\left(q\right)-Q\left(\varrho \right)]K(q,\varrho )d\varrho ,n-1<\wp \le n\end{array}$$

(4)

Definition 2.

The fractional CF integral is stated as [32]

$$\begin{array}{c}\hfill {}^{CF}{I}_q^{\wp }\vartheta (ϵ,q)=\frac{1-\wp }{S(\wp )}Q\left(q\right)+\frac{\wp }{S(\wp )}{\int }_0^{q}Q\left(\varrho \right)d\varrho ,q\ge 0,\wp ϵ(0,1].\end{array}$$

(5)

Definition 3.

For $S(\wp )=1$, the Laplace transform of the CF derivative is [33]

$$\begin{array}{c}\hfill L\left[{}^{CF}{D}_q^{\wp }Q\left[\left(q\right)\right]\right]=\frac{vL\left[Q\right(q)-Q(0\left)\right]}{v+\wp (1-v)}.\end{array}$$

(6)

Definition 4.

The $\mathbb{Y}$T of $Q\left(q\right)$ is framed as [34]

$$\begin{array}{c}\hfill \mathbb{Y}\left[Q\left(q\right)\right]=\chi \left(v\right)={\int }_0^{\infty }Q\left(q\right)e^{-\frac{q}{v}}dq.q>0\end{array}$$

(7)

Remarks

The $\mathbb{Y}$T of some helpful expressions are as follows:

$$\begin{array}{c}\hfill \mathbb{Y}\left[1\right]=v;\\ \hfill \mathbb{Y}\left[q\right]=v^2;\\ \hfill \mathbb{Y}\left[q^i\right]=\mathsf{\Gamma }(i+1)v^{i+1}.\end{array}$$

Lemma 1.

Let the Laplace transform of $Q\left(q\right)$ be $F\left(v\right)$, and then $\chi \left(v\right)=F(1/v)$ [35].

Proof. 

From Equation (7), we can obtain the Yang transform by putting $q/v=\zeta$ as

$$\begin{array}{c}\hfill L\left[Q\left(q\right)\right]={\int }_0^{\infty }Q\left(v\zeta \right)e^{\zeta }d\zeta ,\zeta >0\end{array}$$

(8)

since $L\left[Q\right(q\left)\right]=F\left(v\right)$, which implies that

$$\begin{array}{c}\hfill F\left(v\right)=L\left[Q\left(q\right)\right]={\int }_0^{\infty }Q\left(q\right)e^{-vq}dq.\end{array}$$

(9)

Putting q = $\zeta /v$ in Equation (9), we obtain

$$\begin{array}{c}\hfill F\left(v\right)=\frac{1}{v}{\int }_0^{\infty }Q\left(\frac{\zeta }{v}\right)e^{\zeta }d\zeta .\end{array}$$

(10)

Thus, from Equation (8), we obtain:

$$\begin{array}{c}\hfill F\left(v\right)=\chi \left(\frac{1}{v}\right).\end{array}$$

(11)

Furthermore, from Equations (7) and (9), we obtain

$$\begin{array}{c}\hfill F\left(\frac{1}{v}\right)=\chi \left(v\right).\end{array}$$

(12)

The links between Equations (11) and (12) represent the duality connection among the Laplace and Yang transforms.  □

Lemma 2.

Let $Q\left(q\right)$ be a function, then $\mathbb{Y}$T of CF derivatives of $Q\left(q\right)$ is [35]

$$\begin{array}{c}\hfill \mathbb{Y}\left[Q\left(q\right)\right]=\frac{\mathbb{Y}\left[Q\left(q\right)-vQ\left(0\right)\right]}{v+\wp (v-1)}.\end{array}$$

(13)

Proof. 

The fractional Laplace transform of CF is defined as in Equation (13)

$$\begin{array}{c}\hfill L\left[Q\left(q\right)\right]=\frac{L\left[vQ\left(q\right)-Q\left(0\right)\right]}{v+\wp (1-v)}.\end{array}$$

(14)

However, we have a correlation among the $\mathbb{Y}$T and Laplace properties, namely $\chi \left(v\right)=F(1/v)$, so put $1/v$ for v in Equation (14), and we obtain

$$\begin{array}{cc}\hfill \mathbb{Y}\left[Q\left(q\right)\right]& =\frac{\mathbb{Y}\left[\frac{1}{v}Q\left(q\right)-Q\left(0\right)\right]}{\frac{1}{v}+\wp (1-\frac{1}{v})},\hfill \\ \hfill \mathbb{Y}\left[Q\left(q\right)\right]& =\frac{\mathbb{Y}\left[Q\left(q\right)-vQ\left(0\right)\right]}{1+\wp (v-1)}.\hfill \end{array}$$

(15)

Thus, the proof is satisfied.  □

3. Idea of Yang Homotopy Perturbation Transform Method (${Y}$HPTM)

In this part, we will demonstrate the concept of $\mathbb{Y}$HPTM. Let us assume a nonlinear fractional-order PDE, such as

$$\begin{array}{c}\hfill {}^{CF}{D}_q^{\wp }\vartheta (ϵ,q)+R\vartheta (ϵ,q)+N\vartheta (ϵ,q)=g(ϵ,q),\end{array}$$

(16)

$$\begin{array}{c}\hfill \vartheta (ϵ,0)=h\left(ϵ\right),\end{array}$$

(17)

where $g(ϵ,q)$ is called the source function. Applying the $\mathbb{Y}$T to Equation (16),

$$\begin{array}{c}\hfill \frac{1}{v^{\wp }}\mathbb{Y}\left[\vartheta (ϵ,q)-v\vartheta (ϵ,0)\right]=-\mathbb{Y}[R\left(\vartheta (ϵ,q)\right)+N\left(\vartheta (ϵ,q)\right)+\mathbb{Y}\left[g(ϵ,q)\right]],\\ \hfill \mathbb{Y}\left[\vartheta (ϵ,q)\right]=vh\left(ϵ\right)-v^{\wp }\left[\mathbb{Y}\left[R\right(\vartheta (ϵ,q))+N(\vartheta (ϵ,q)\left)\right]\right]+\mathbb{Y}\left[g(ϵ,q)\right].\end{array}$$

By using inverse $\mathbb{Y}$T,

$$\begin{array}{c}\hfill \vartheta (ϵ,q)=\vartheta (ϵ,0)-{\mathbb{Y}}^{-1}\left[v^{\wp }\left[\mathbb{Y}\left[R\right(\vartheta (ϵ,q))+N(\vartheta (ϵ,q)\left)\right]\right]+\mathbb{Y}\left[g(ϵ,q)\right]\right].\end{array}$$

(18)

However, HPM is stated as

$$\begin{array}{c}\hfill \vartheta (ϵ,q)=\sum _{i=0}^{\infty }{p}^i{\vartheta }_i(ϵ,q),\end{array}$$

(19)

and

$$\begin{array}{c}\hfill N\vartheta (ϵ,q)=\sum _{i=0}^{\infty }{p}^i{H}_i\vartheta (ϵ,q).\end{array}$$

(20)

The following strategy can be operated to acquire the He’s polynomials,

$$\begin{array}{c}\hfill H_i({\vartheta }_0+{\vartheta }_1+\cdots +{\vartheta }_i)=\frac{1}{n!}\frac{{\partial }^i}{\partial p^i}{\left(N(\sum _{i=0}^{\infty }{p}^i{\vartheta }_i)\right)}_{p=0}.n=0,1,2,\cdots \end{array}$$

With the help of Equations (19) and (20), we can obtain Equation (18), such as

$$\begin{array}{c}\hfill \sum _{i=0}^{\infty }{p}^i{\vartheta }_i(ϵ,q)=\vartheta (ϵ,0)-p{\mathbb{Y}}^{-1}\left[v^{\wp }\mathbb{Y}\left\{R\left(\sum _{i=0}^{\infty }{p}^i{\vartheta }_i(ϵ,q)\right)+\sum _{i=0}^{\infty }{p}^i{H}_n{\vartheta }_i(ϵ,q)\right\}\right].\end{array}$$

We can obtain the following terms by evaluating the p components:

$$\begin{array}{cc}\hfill p^0& ={\vartheta }_0(ϵ,q)=\vartheta (ϵ,0),\hfill \\ \hfill p^1& ={\vartheta }_1(ϵ,q)=-{\mathbb{Y}}^{-1}\left[v^{\wp }\mathbb{Y}\left\{R{\vartheta }_0(ϵ,q)+H_0\left(\vartheta \right)\right\}\right],\hfill \\ \hfill p^2& ={\vartheta }_2(ϵ,q)=-{\mathbb{Y}}^{-1}\left[v^{\wp }\mathbb{Y}\left\{R{\vartheta }_1(ϵ,q)+H_1\left(\vartheta \right)\right\}\right],\hfill \\ \hfill p^3& ={\vartheta }_3(ϵ,q)=-{\mathbb{Y}}^{-1}\left[v^{\wp }\mathbb{Y}\left\{R{\vartheta }_2(ϵ,q)+H_2\left(\vartheta \right)\right\}\right],\hfill \\ \hfill & ⋮\hfill \\ \hfill p^i& ={\vartheta }_i(ϵ,q)=-{\mathbb{Y}}^{-1}\left[v^{\wp }\mathbb{Y}\left\{R{\vartheta }_i(ϵ,q)+H_i\left(\vartheta \right)\right\}\right],\hfill \end{array}$$

(21)

Thus, we can summarize the set of Equations (21) in the series form, such as

$$\begin{array}{cc}\hfill \vartheta (ϵ,q)& ={\vartheta }_0(ϵ,q)+{\vartheta }_1(ϵ,q)+{\vartheta }_2(ϵ,q)+\cdots \hfill \\ \hfill \vartheta (ϵ,q)& =\underset{N\to \infty }{lim}\sum _{n=0}^N{\vartheta }_n(ϵ,q)\hfill \end{array}$$

(22)

4. Numerical Applications

4.1. Example 1

Consider a linear time-fractional KG problem

$$\begin{array}{c}\hfill D_q^{\wp }\vartheta (ϵ,q)-D_{ϵ}^2\vartheta (ϵ,q)-\vartheta (ϵ,q)=0,\end{array}$$

(23)

with the initial condition

$$\begin{array}{c}\hfill \vartheta (ϵ,0)=1+sin\left(ϵ\right).\end{array}$$

(24)

Taking $\mathbb{Y}$T of Equation (23), we obtain

$$\begin{array}{c}\hfill \mathbb{Y}\left[\frac{{\partial }^{\wp }\vartheta }{\partial q^{\wp }}\right]=\mathbb{Y}\left[\frac{{\partial }^2\vartheta }{\partial {ϵ}^2}+\vartheta \right].\end{array}$$

Executing the differential property of $\mathbb{Y}$T, we obtain

$$\begin{array}{c}\hfill \frac{1}{v^{\wp }}\mathbb{Y}\left[\vartheta (ϵ,q)-v\vartheta (ϵ,0)\right]=\mathbb{Y}\left[\frac{{\partial }^2\vartheta }{\partial q^2}+\vartheta \right],\end{array}$$

$$\begin{array}{c}\hfill \mathbb{Y}\left[\vartheta (ϵ,q)\right]=v\vartheta (ϵ,0)+v^{\wp }\mathbb{Y}\left[\frac{{\partial }^2\vartheta }{\partial q^2}+\vartheta \right].\end{array}$$

The inverse $\mathbb{Y}$T indicates

$$\begin{array}{c}\hfill \vartheta (ϵ,q)=\vartheta (ϵ,0)+{\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}\left(\frac{{\partial }^2\vartheta }{\partial q^2}+\vartheta \right)\right\}\right].\end{array}$$

Employing HPM such as

$$\begin{array}{c}\hfill \vartheta (ϵ,q)={\vartheta }_0+p{\vartheta }_1+p^2{\vartheta }_2+\cdots ,\end{array}$$

$$\begin{array}{c}\hfill \sum _{i=0}^{\infty }{p}^i{\vartheta }_i(ϵ,q)=1+sin\left(q\right)+p\left({\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}\left(\sum _{i=0}^{\infty }{p}^i\frac{{\partial }^2{\vartheta }_i}{\partial q^2}+\sum _{i=0}^{\infty }{p}^i{\vartheta }_i\right)\right\}\right]\right),\end{array}$$

on comparing the identical of p, we obtain

$$\begin{array}{cc}\hfill p^0& ={\vartheta }_0(ϵ,q)=\vartheta (ϵ,0),\hfill \\ \hfill p^1& ={\vartheta }_1(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_0}{\partial q^2}+{\vartheta }_0)\right\}\right],\hfill \\ \hfill p^2& ={\vartheta }_2(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_1}{\partial q^2}+{\vartheta }_1)\right\}\right],\hfill \\ \hfill p^3& ={\vartheta }_3(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_2}{\partial q^2}+{\vartheta }_2)\right\}\right],\hfill \\ \hfill & ⋮\hfill \end{array}$$

With help of Equation (24), we gain the iterations successively ${\vartheta }_i\left(ϵ\right)$, $i=1,2,3,\cdots$, as follows:

$$\begin{array}{cc}\hfill {\vartheta }_0(ϵ,q)=& 1+sin\left(ϵ\right),\hfill \\ \hfill {\vartheta }_1(ϵ,q)=& \frac{1}{\mathsf{\Gamma }(1+\wp )}{q}^{\wp },\hfill \end{array}$$

$$\begin{array}{cc}\hfill {\vartheta }_2(ϵ,q)=& \frac{1}{\mathsf{\Gamma }(1+2\wp )}{q}^{2\wp },\hfill \\ \hfill {\vartheta }_3(ϵ,q)=& \frac{1}{\mathsf{\Gamma }(1+3\wp )}{q}^{3\wp },\hfill \\ \hfill & ⋮\hfill \\ \hfill {\vartheta }_i(ϵ,q)=& \frac{1}{\mathsf{\Gamma }(1+i\wp )}{q}^{i\wp },\hfill \end{array}$$

Thus, the approximate solution can be obtained by:

$$\begin{array}{cc}\hfill \vartheta (ϵ,q)=& 1+sin\left(ϵ\right)+\frac{1}{\mathsf{\Gamma }(1+\wp )}{q}^{\wp }+\frac{1}{\mathsf{\Gamma }(1+2\wp )}{q}^{2\wp }+\frac{1}{\mathsf{\Gamma }(1+3\wp )}{q}^{3\wp }+\cdots \hfill \\ \hfill =& 1+sin\left(ϵ\right)+\sum _{i=0}^{\infty }{p}^i\frac{q^{i\wp }}{\mathsf{\Gamma }(1+i\wp )},\hfill \end{array}$$

(25)

which implies the exact solution of Equation (23), In particular, at $\wp =1$, we obtain

$$\begin{array}{cc}\hfill \vartheta (ϵ,q)=& 1+sin\left(ϵ\right),\hfill \end{array}$$

(26)

which is in full agreement.

Figure 1a–d indicate the physical behavior of the obtained solution at $ϵ\in [0,4]$ and $q\in [0,0.8]$. From these figures, it can be observed that the solution graphs of the problem show the friendly touch with each other. Figure 1a–d demonstrate that the solution achieved by $\mathbb{Y}$HPTM approaches the precise solution very rapidly with more iterations. In Figure 2, we have plotted the graph of $\vartheta (ϵ,q)$ at different fractional order of $\wp =0.25,0.50,0.75,1$ and $ϵ\in [0,2\pi ]$ with different values of q.

4.2. Example 2

Assume a nonlinear time-fractional KG problem

$$\begin{array}{c}\hfill D_q^{\wp }\vartheta (ϵ,q)-D_{ϵ}^2\vartheta (ϵ,q)+{\vartheta }^2(ϵ,q)=0,\end{array}$$

(27)

with the initial condition

$$\begin{array}{c}\hfill \vartheta (ϵ,0)=1+sin\left(ϵ\right).\end{array}$$

(28)

Taking the Yang transform of Equation (27), we obtain

$$\begin{array}{c}\hfill \mathbb{Y}\left[\frac{{\partial }^{\wp }\vartheta }{\partial q^{\wp }}\right]=\mathbb{Y}\left[\frac{{\partial }^2\vartheta }{\partial {ϵ}^2}-{\vartheta }^2\right].\end{array}$$

Executing the differential property of $\mathbb{Y}$T, we obtain

$$\begin{array}{c}\hfill \frac{1}{v^{\wp }}\mathbb{Y}\left[\vartheta (ϵ,q)-v\vartheta (ϵ,0)\right]=\mathbb{Y}\left[\frac{{\partial }^2\vartheta }{\partial q^2}-{\vartheta }^2\right],\end{array}$$

$$\begin{array}{c}\hfill \mathbb{Y}\left[\vartheta (ϵ,q)\right]=v\vartheta (ϵ,0)+v^{\wp }\mathbb{Y}\left[\frac{{\partial }^2\vartheta }{\partial q^2}-{\vartheta }^2\right].\end{array}$$

The inverse $\mathbb{Y}$T indicates

$$\begin{array}{c}\hfill \vartheta (ϵ,q)=\vartheta (ϵ,0)+{\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}\left(\frac{{\partial }^2\vartheta }{\partial q^2}-{\vartheta }^2\right)\right\}\right].\end{array}$$

Employing HPM such as

$$\begin{array}{c}\hfill \sum _{i=0}^{\infty }{p}^i{\vartheta }_i(ϵ,q)=1+sin\left(q\right)+p\left({\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}\left(\sum _{i=0}^{\infty }{p}^i\frac{{\partial }^2{\vartheta }_i}{\partial q^2}-\sum _{i=0}^{\infty }{p}^i{\vartheta }_i^2\right)\right\}\right]\right),\end{array}$$

on comparing the identical of p, we obtain

$$\begin{array}{cc}\hfill p^0& ={\vartheta }_0(ϵ,q)=\vartheta (ϵ,0),\hfill \\ \hfill p^1& ={\vartheta }_1(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_0}{\partial q^2}-{\vartheta }_0^2)\right\}\right],\hfill \\ \hfill p^2& ={\vartheta }_2(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_1}{\partial q^2}-2{\vartheta }_0{\vartheta }_1)\right\}\right],\hfill \\ \hfill p^3& ={\vartheta }_3(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_2}{\partial q^2}-{\vartheta }_1^2-2{\vartheta }_0{\vartheta }_2)\right\}\right],\hfill \\ \hfill & ⋮\hfill \end{array}$$

With help of Equation (28), we gain the iterations successively ${\vartheta }_i\left(ϵ\right)$, $i=1,2,3,\cdots$, as follows:

$$\begin{array}{cc}\hfill {\vartheta }_0(ϵ,q)=& 1+sin\left(ϵ\right),\hfill \\ \hfill {\vartheta }_1(ϵ,q)=& \frac{-q^{\wp }}{\mathsf{\Gamma }(1+\wp )}\left(1+3sin(ϵ)+{sin}^2(ϵ)\right),\hfill \\ \hfill {\vartheta }_2(ϵ,q)=& \frac{q^{2\wp }}{\mathsf{\Gamma }(1+2\wp )}\left(11sin(ϵ)+12{sin}^2(ϵ)+2{sin}^3(ϵ)\right),\hfill \\ \hfill {\vartheta }_3(ϵ,q)=& \frac{q^{3\wp }}{\mathsf{\Gamma }(1+3\wp )}\left(18-57sin(ϵ)-160{sin}^2(ϵ)-82{sin}^3(ϵ)-10sin(4ϵ)\right),\hfill \\ \hfill & ⋮\hfill \end{array}$$

Thus, the approximate solution can be obtained by:

$$\begin{array}{cc}\hfill \vartheta (ϵ,q)=& 1+sin\left(ϵ\right)-\frac{q^{\wp }}{\mathsf{\Gamma }(1+\wp )}\left(1+3sin\left(ϵ\right)+{sin}^2\left(ϵ\right)\right)+\frac{q^{2\wp }}{\mathsf{\Gamma }(1+2\wp )}\left(11sin\left(ϵ\right)+12{sin}^2\left(ϵ\right)+2{sin}^3\left(ϵ\right)\right)\hfill \\ & +\frac{q^{3\wp }}{\mathsf{\Gamma }(1+3\wp )}\left(18-57sin\left(ϵ\right)-160{sin}^2\left(ϵ\right)-82{sin}^3\left(ϵ\right)-10sin\left(4ϵ\right)\right)+\cdots \hfill \end{array}$$

Figure 3a–d indicate the physical behavior of the obtained solution at $ϵ\in [0,1]$ and $q\in [0,1]$. From these figures, it can be observed that with the increase in the value of ℘, the approximate solution become close to the exact solution at $\wp =1$. In Figure 4, we have plotted the graph of $\vartheta (ϵ,q)$ with different fractional orders of $\wp =0.25,0.50,0.75,1$ at $ϵ\in [0,2\pi ]$ with different values of q. It is obvious that this approximation can only be employed numerically, even though a closed form solution is not accessible. It can be seen that our approximate solution using $\mathbb{Y}$HPTM in

Table 1. Comparison between the value $\vartheta (ϵ,q)$ for the solution of the KG equation.

Sr. No.	$\mathit{q}=0.1$			$\mathit{q}=0.2$			$\mathit{q}=0.3$
$ϵ$	[36]	[37]	$\mathbb{Y}$HPTM	[36]	[37]	$\mathbb{Y}$HPTM	[36]	[37]	$\mathbb{Y}$HPTM
0.0	0.9949999861	0.9950000249	0.903	0.9799991162	0.9800015775	0.824	0.9549900052	0.9550176534	0.781
0.1	1.093291132	1.093291179	0.976100	1.073723730	1.073726319	0.871321	1.073723730	1.073726319	0.792208
0.2	1.190502988	1.190503087	1.04725	1.166134875	1.166138050	0.915126	1.125945576	1.125974851	0.794835
0.3	1.285668610	1.285668848	1.11584	1.256326130	1.256331032	0.955409	1.208114007	1.208147932	0.789972
0.4	1.377844211	1.377844710	1.18132	1.343423788	1.343432104	0.992136	1.287043874	1.287088824	0.778571
0.5	1.466118315	1.466119219	1.24317	1.426594492	1.426608263	1.0252	1.362025218	1.362089477	0.761295
0.6	1.549620480	1.549621939	1.3009	1.505052082	1.505073495	1.05442	1.432404521	1.432497282	0.738476
0.7	1.627529538	1.627531694	1.35406	1.578063673	1.578094808	1.07951	1.497587424	1.497717706	0.710192
0.8	1.699081273	1.699084244	1.40223	1.644954933	1.644997540	1.0023	1.557040327	1.557215916	0.676451
0.9	1.763575490	1.763579356	1.44504	1.705114628	1.705169916	1.11635	1.610291023	1.610517519	0.63744
1.0	1.820382425	1.820387216	1.48219	1.757998450	1.758066925	1.12781	1.656928567	1.657208637	0.593784

is more significant than that obtained in [36,37].

4.3. Example 3

Consider another nonlinear time-fractional KG problem

$$\begin{array}{c}\hfill D_q^{\wp }\vartheta (ϵ,q)-D_{ϵ}^2\vartheta (ϵ,q)+\vartheta (ϵ,q)-{\vartheta }^3(ϵ,q)=0,\end{array}$$

(29)

with the initial condition

$$\begin{array}{c}\hfill \vartheta (ϵ,0)=-sech\left(ϵ\right).\end{array}$$

(30)

According to the idea of $\mathbb{Y}$HPTM, we can obtain the following relation

$$\begin{array}{c}\hfill \sum _{i=0}^{\infty }{p}^i{\vartheta }_i(ϵ,q)=sech\left(ϵ\right)+p\left({\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}\left(\sum _{i=0}^{\infty }{p}^i\frac{{\partial }^2{\vartheta }_i}{\partial q^2}-\sum _{i=0}^{\infty }{p}^i{\vartheta }_i+\sum _{i=0}^{\infty }{p}^i{\vartheta }_i^2\right)\right\}\right]\right),\end{array}$$

when the coefficients of like powers of p are compared, we obtain

$$\begin{array}{cc}\hfill p^0& ={\vartheta }_0(ϵ,q)=\vartheta (ϵ,0),\hfill \\ \hfill p^1& ={\vartheta }_1(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_0}{\partial q^2}-{\vartheta }_0+{\vartheta }_0^3)\right\}\right],\hfill \\ \hfill p^2& ={\vartheta }_2(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_1}{\partial q^2}-{\vartheta }_1+3{\vartheta }_0^2{\vartheta }_1)\right\}\right],\hfill \\ \hfill p^3& ={\vartheta }_3(ϵ,q)={\mathbb{Y}}^{-1}\left[v^{\wp }\left\{\mathbb{Y}(\frac{{\partial }^2{\vartheta }_2}{\partial q^2}-{\vartheta }_2+3{\vartheta }_0{\vartheta }_1^2+3{\vartheta }_0^2{\vartheta }_2)\right\}\right],\hfill \\ \hfill & ⋮\hfill \end{array}$$

with help of Equation (30), we gain the iterations successively ${\vartheta }_i\left(ϵ\right)$, $i=1,2,3,\cdots$, as follows:

$$\begin{array}{cc}\hfill {\vartheta }_0(ϵ,q)=& -sech\left(ϵ\right),\hfill \\ \hfill {\vartheta }_1(ϵ,q)=& -\frac{q^{\wp }}{\mathsf{\Gamma }(1+\wp )}\left(2sech\left(ϵ\right)-3{sech}^3\left(ϵ\right)\right),\hfill \\ \hfill {\vartheta }_2(ϵ,q)=& -\frac{q^{2\wp }}{\mathsf{\Gamma }(1+2\wp )}\left(3sech\left(ϵ\right)-34{sech}^3\left(ϵ\right)-18{sech}^5\left(ϵ\right)\right),\hfill \\ \hfill {\vartheta }_3(ϵ,q)=& -\frac{q^{3\wp }}{\mathsf{\Gamma }(1+3\wp )}\left(64{sech}^3\left(x\right)-288{sech}^5\left(ϵ\right)+240{sech}^7\left(ϵ\right)\right),\hfill \\ \hfill & ⋮\hfill \end{array}$$

Thus, the approximate solution can be obtained by:

$$\begin{array}{cc}\hfill \vartheta (ϵ,q)=& -sech\left(ϵ\right)-\frac{q^{\wp }}{\mathsf{\Gamma }(1+\wp )}\left(2sech\left(ϵ\right)-3{sech}^3\left(ϵ\right)\right)-\frac{q^{2\wp }}{\mathsf{\Gamma }(1+2\wp )}\left(3sech\left(ϵ\right)-34{sech}^3\left(ϵ\right)-18{sech}^5\left(ϵ\right)\right)\hfill \\ & -\frac{q^{3\wp }}{\mathsf{\Gamma }(1+3\wp )}\left(64{sech}^3\left(ϵ\right)-288{sech}^5\left(ϵ\right)+240{sech}^7\left(ϵ\right)\right)+\cdots \hfill \end{array}$$

Figure 5a–d indicates the physical behavior of the obtained solution at $ϵ\in [-2,2]$ and $q\in [0,0.1]$. From these figures, it can be observed that with increase in the value of ℘, the approximate solution graph comes close to the exact exact solution at $\wp =1$. In Figure 6, we have plotted the graph of $\vartheta (ϵ,q)$ at different fractional orders of $\wp =0.25,0.50,0.75,1$ and $ϵ\in [0,2\pi ]$ with different values of q. We compared our graphical results obtained by $\mathbb{Y}$HPTM, which converges to the exact solution very rapidly with a small amount of q compared to [38] at $\wp =1$.

5. Conclusions

In this study, $\mathbb{Y}$HPTM has been utilized to achieve an approximate solution of nonlinear time-fractional KG equations. We demonstrate some illustrations to show the validity of the method, which reveals that the obtained findings are satisfactory. It is important to note that in order to improve the accuracy, a greater number of iterations with excessive order of p are required. The best advantage of $\mathbb{Y}$HPTM is that it generates the approximate solution in a quickly convergent power series form. As a result, this strategy can be adopted to elucidate a wide classification of nonlinear challenges in science and engineering with no need for linearization, discretization or perturbation. The proposed strategy has the privilege of being able to tackle linear and nonlinear time-fractional KG problems simultaneously. Mathematica package 11.0.1 has been operated for the graphical analysis as well as for the computations in this paper. We recommend the readers to consider this problem for the Atangana–Baleanu fractional derivative and many others in the place of the Caputo–Fabrizio operator.