 Research Note
 Open access
 Published:
A hybrid yang transform adomian decomposition method for solving timefractional nonlinear partial differential equation
BMC Research Notes volumeÂ 17, ArticleÂ number:Â 226 (2024)
Abstract
Nonlinear timefractional partial differential equations (NTFPDEs) play a great role in the mathematical modeling of realworld phenomena like traffic models, the design of earthquakes, fractional stochastic systems, diffusion processes, and control processing. Solving such problems is reasonably challenging, and the nonlinear part and fractional operator make them more problematic. Thus, developing suitable numerical methods is an active area of research. In this paper, we develop a new numerical method called Yang transform Adomian decomposition method (YTADM) by mixing the Yang transform and the Adomian decomposition method for solving NTFPDEs. The derivative of the problem is considered in sense of Caputo fractional order. The stability and convergence of the developed method are discussed in the Banach space sense. The effectiveness, validity, and practicability of the method are demonstrated by solving four examples of NTFPEs. The findings suggest that the proposed method gives a better solution than other compared numerical methods. Additionally, the proposed scheme achieves an accurate solution with a few numbers of iteration, and thus the method is suitable for handling a wide class of NTFPDEs arising in the application of nonlinear phenomena.
Introduction
Different natural phenomena that appear in various areas of science and engineering can be successfully demonstrated by using the notion of fractional calculus [1]. Fractional calculus plays a great role in the mathematical modeling of realworld phenomena, for example, reactionâ€“diffusion processes [2, 3], decentralized wireless networks [4, 5], water wave movement [6, 7], signal processing [8, 9], population growth [10, 11], design of earthquake [12, 13], traffic models with fractional derivatives [14], diffusion processes [15], fractional stochastic systems [16,17,18,19], control processing [20, 21], medical sciences [22] and many other physical processes [23,24,25,26].
In most cases, fractional differential equations (FDEs) are considered a generalization of differential equations (DEs) because they describe functional values at every continuous point without losing any memory or hereditary behaviors of natural phenomena. For this reason, many authors have studied FDEs for modeling and deep understanding of realworld natural phenomena such as business models with fractional derivations, evaluation of viscoelastic surface characteristics, human behavior in the field of mathematical psychology [27, 28], models of love between humans [29], and models of happiness and in other fields of applications [23,24,25, 30, 31].
Nonlinear timefractional partial differential equations (NTFPDEs) have played an alternative role in diverse fields such as applied science, mathematics, physics, medicine, biology, geology, chemistry, and engineering due to the nonlocality nature of the powerlaw structure with arbitrary positive order. In the modern age, it is impossible to model many realworld problems without using nonlinear FPDEs. Fractional calculus can be called this centuryâ€™s calculus because of the diversity of applications in every discipline of applied science, technology, and engineering [32,33,34].
In applied mathematics and many other areas of science, it is important to obtain accurate solutions to FPDEs. In general, FPDEs do not have closedform exact solutions for most problems, and it is very difficult to obtain exact solutions [35]. For this reason, various numerical methods have been developed for solving linear and nonlinear FPDEs. For instance, Mohamed and Torky [33] proposed the Legendre multiwavelet Galerkin method to solve the conformable fractional modified KdVâ€“ZK equation. Singh et al. [36] announced the homotopy perturbation Sumudu transform method for nonlinear FPDEs arising in the spatial diffusion of biological populations. Wang and Liu [37] familiarized themselves with the new Sumudu transform iterative method for the timeâ€‘fractional Cauchy reactionâ€“diffusion equation. Ziane et al. [38] introduced the Elzaki transform combined with the variational iteration method for PDEs of fractional order. Dehestani et al. [39] proposed Legendreâ€“Laguerre functions based on Legendre and Laguerre polynomials for solving a class of linear and nonlinear timeâ€“space fractional PDEs with variable coefficients. Dehestani et al. [40] presents the discrete scheme based on Genocchi polynomials and fractional Laguerre functions to solve multiterm variableorder timefractional PDEs in the large interval. Dehestani et al. [41] familiarized themselves with the fractionalorder Genocchiâ€“Petrovâ€“Galerkin method to investigate the approximate solution of timeâ€“space fractional Fokkerâ€“Planck equations. Wang et al. [42] proposed two analytical methods, the residual power series method and the homotopy analysis transform method, to solve NTFPDEs with proportional delay.
Recently, researchers have used several numerical methods to solve different types of NTFPDEs. For example, Bekela et al. [32] applied the hybrid numerical method of Laplacelike transform and variational theory for solving NTFPDEs with proportional delay. Malyk et al. [43] employed the Yangâ€“Abdelâ€“Cattani derivative operator for solving the nonlinear fractional diffusion equation. The numerical method is presented by using a novel spline technique for solving fourthorder timefractional evolution problems [44]. AlDeiakeh et al. [45] introduced an approximation for nonlinear FPDE using the combination of the multiLaplace transform and ADM. Mallick et al. [46] proposed an iterative method for solving timefractional PDEs with proportional delays. A wide application of NTFPDEs in many real life phenomena and the growing of to search for suitable numerical methods motivate us to propose a hybrid numerical method called the Yang transform Adomian decomposition method (YTADM) for solving NTFPDEs.
George Adomian introduced a modern mathematical method to solve nonlinear DEs in the 1980s, described as ADM [47,48,49]. Similarly, another powerful method found by Yang to solve linear PDEs was described as YT, which transforms the initial DEs into an algebraic equation [50,51,52]. However, it does not handle nonlinear terms in NTFPDEs because of the difficulties that are caused by nonlinear terms. ADM is a powerful numerical method for solving nonlinear DEs, but it does not properly decompose the fractional order parts. The main objective of this study is to combine two methods for solving the NTFPDEs: YT is used to decompose FDEs into algebraic equations, and ADM is used to decompose the nonlinear terms in NTFPDEs into a series solution. The stability and convergence conditions of the proposed method were analyzed. In addition, some test examples of NTFPDEs are solved to illustrate the capability, accuracy, and simplicity of the proposed method.
Furthermore, to demonstrate the versatility and robustness of the YTADM in different areas of science, we take four examples that have a wide range of applications in various fields of science and technology. Accordingly, first we considered the nonlinear timefractional NewellWhiteheadSegel equation (NWSE), which predicts the existence of traveling wave patterns. These kinds of systems can be observed in a variety of natural systems, including physical, biological, and chemical systems. For example, in medicine: the mechanisms of spreading epidemic diseases, healing wounded tissue, chemical reactions of multiple substances, and neutron diffusion theory all belong to reactionâ€“diffusion systems [49].
Then, we considered the nonlinear fractional Cauchy reactionâ€“diffusion equation (CRDE), which is used to describe the evolution of a system over time with reaction and diffusion processes in various fields. For instance, in biology and medicine, these equations are used to model the spreading of diseases, tumor growth, etc. In chemistry, CRDEs are employed to understand chemical reactions that exhibit complex kinetics and diffusion processes. In population dynamics, these equations are employed to model the behavior of ecological and natural populations. In material science, CRDEs are expended for modeling transport marvels in porous media whereas in finance and economics, they are used for modeling the financial markets with memory goods and longrange dependence [47, 48], and references therein.
Hereafter, we solve the nonlinear timefractional FornbergWhitham equation (FWE). The FWE is a type of traveling wave equation and have variety of applications in physics and engineering. Such as in the propagation of electrical signals and optimization of guided communication systems. In mathematical physics to investigate how nonlinear dispersive water waves break [50]. Final, we considered one problem from nonlinear timefractional Fokkerâ€“Planck equation (FPE). The nonlinear timefractional FPE arises in various fields like in chemistry, natural science, including astrophysics problems, biological applications, chemical physics, and other fields [51].
Here, the numerical results obtained by YTADM are summarized in tables and figures to illustrate easily understanding for various fractional order \(\alpha\). For the fractional order \(\alpha =1\), the obtained numerical solution have excellent agreement with the exact solution. It is evident from the theory of fractional calculus that the approximate solution of the problem with derivative \(1\) continuously tends to the exact solution when the fractional derivative \(\alpha (m1<\alpha \le m)\) tends to positive integer \(m=1\)[52].
The remaining paper is organized in the following way: Sect.Â "Preliminaries" of this work gives detailed concepts related to fractional integrals, derivatives, and YT. Sect.Â "A hybrid numerical method for solving NTFPDEs" familiarizes reader with the mathematical formulation of YTADM using YT and ADM. Sects.Â "Stability analysis of YTADM to solve NTFPDEs" and "Convergence analysis of YTADM for solving NTFPDEs" offer stability and convergence of the proposed method in the Banach space sense. Sect.Â "Numerical results and discussions" presents a numerical simulation of the method on some NTFPDEs. Finally, the conclusion of the method is presented in Sect.Â "Conclusion".
Preliminaries
This part is devoted to some basic concepts and definitions of YT, fractional integrals, and fractional derivatives, which are essential for accepting the remainder of the monograph.
Definition 2.1
[6, 7]. A real function \(f\left(t\right), t>0\) is in the space \({C}_{\tau }, \tau \in R\) if there exists a real number \(\left(p>\tau \right)\), such that \(f\left(t\right)={t}^{p}{f}_{1}\left(t\right)\), where \({f}_{1}\left(t\right)\in C\left[0,\infty )\right.\), and it is said to be in the space \({C}_{\tau }^{n}\) if \({f}^{\left(n\right)}\in {C}_{\tau }, n\in N\), clearly \({C}_{\tau }\subset {C}_{\rho }\) if \(\rho \le \tau\).
Definition 2.2
[24, 25]. The left sides Riemannâ€“Liouville fractional integral of order \(\alpha \ge 0\), of a function \(f\left(t\right)\in {C}_{\tau },\tau \ge 1\) is defined as.
where \(\Gamma \left(.\right)\) is the gamma function.
Definition 2.3
[7, 55]. The fractional derivative of \(f(t)\) in the Caputo sense is defined as.
where \({{}_{a}{}^{C}D}_{t}^{\alpha }\) is the Caputo fractional derivative operator of order \(\alpha\), \(n1<\alpha \le n, n\in N\).
Definition 2.4
[56, 57]. The YT is defined over the set of function.
where \(C\) is constant and \({m}_{1},{m}_{2}\) are either finite or infinite by the following integral formula
It converges if the limit of the integral exists, and diverges if not. Accordingly, the following results follows directly from the above integral formula of YT [56, 59]:

1.
\(Y\left[1\right]=v\),

2.
\(Y\left[t\right]={v}^{2}\),

3.
\(Y\left[{t}^{n}\right]=n!{v}^{n+1}=\Gamma \left(n+1\right){v}^{n+1}, n\in N,\)

4.
\(Y\left[{t}^{\alpha }\right]=\Gamma \left(\alpha +1\right){v}^{\alpha +1}, \alpha >0,\)

5.
\(Y\left[af\left(t\right)+bg\left(t\right)\right]=aY\left[f\left(t\right)\right]+bY\left[g\left(t\right)\right],\) linearity property.
Definition 2.5
[56, 59]. The YT of the \({n}^{th}\) derivative of the function \(f\left(t\right)\) is defined as
Definition 2.6
[58]. The YT of the function \(f\left(t\right)\) in the Caputo fractional derivative is defined as.
for \(T\left(v\right)=Y\left[f\left(t\right)\right]\) and \(n1<\alpha \le n, n\in N\).
Definition 2.7
([32]) Let \(P:X\to X\) be a mapping of a set \(X\) into itself. The fixed point of \(P\) is \(x\in X\) which is mapped onto itself, that is,
the image \(Px\) coincides with \(x\).
Definition 2.8
(Contraction [32]) Let \(X=(X, d)\) be a metric space. A mapping \(P:X\to X\) is called a contraction on \(X\) if there is \(0\le \varepsilon <1\) such that, for all \(x,y\in X\),
Theorem 2.1
(Banachâ€™s fixed point theorem [32]) Consider a metric space \(X=(X, d)\), where \(x\ne \varnothing\). Suppose that \(X\) is complete, and let \(P:X\to X\) be a contraction on \(X\). Then \(P\) has a unique fixed point.
Theorem 2.2
[60]. Let \(\left(X,\Vert .\Vert \right)\) be a Banach space and \(P\) be a selfmap of \(X\) (\(P:X\to X\)), satisfying.
then \(P\) is Picard, \(\text{P}\)stable.
A hybrid numerical method for solving NTFPDEs
In this section, we derive a hybrid numerical method that uses YT and ADM. We named this developed method YTADM. To illustrate the basic idea of this method, we consider the following general NTFPDEs with a source term [38, 58]
with the initial condition
where \({}_{0}{}^{C}{D}_{t}^{\alpha }=\frac{{\partial }^{\alpha }}{\partial {t}^{\alpha }}\) is the Caputo fractional derivative of order \(\alpha\) with respect to \(t\), \(R\) is the linear differential operator, \(N\) represents the general nonlinear differential operator and \(f\left(x,t\right)\) is the source term. To develop YTADM for solving (1) and (2), we follow the following steps.
Step 1: Apply YT on both sides of (1) to obtain
By using Definition 2.6 we transform the fractional derivative on the left side of (3) in to
Note that in this step, YT changes the fractional derivative into an algebraic equation, and using the given initial condition, we can rewrite (4) as
where \(T\left(x,v\right)=Y\left[u\left(x,t\right)\right]\).
By the linearity property of the YT, the above equation becomes
Step 2: Apply the inverse YT on both sides of (5), to obtain
where \(u\left(x,t\right)=\,{Y}^{1}\left[T\left(x,v\right)\right].\)
Step 3: In order to decompose the nonlinear term in NTFPDE and generate a series solution of the proposed method, the ADM is applied. The ADM defines the solution in an infinite series of linear terms as
and the nonlinear term in the problem expressed in terms of the Adomian polynomial as follows:
where
Step 4: Substituting the obtained values in (7) and (8) into (6), we obtain
Step 5: To obtain the component solutions of the proposed method, we equating the terms on both sides of (10), and get the following relation
Step 6: Following component approximate solution, we obtain the general solution in a series form as
Note that it is impossible to find infinite component of solution. Hence, after \({\left(m+1\right)}^{th}\) truncation, the iteration formula of YTADM is given by
which is the \({\left(m+1\right)}^{th}\) approximate solution of the NTFPDE in (1) for \(0<\alpha \le 1\).
Stability analysis of YTADM to solve NTFPDEs
In this section, we discuss the stability analysis of the method presented in Sect.Â "A hybrid numerical method for solving NTFPDEs". For this, we state and prove an essential condition concerning the stability of YTADM when it is applied to solve NTFPDEs. To reveal Picard stability, it enough to show that the mapping associated with YTADM in (13) satisfies the conditions of Theorem 2.2.
Theorem 3.1
Let \(\left(\text{X},\Vert .\Vert \right)\) be a Banach space and \(P:X\to X\) be a selfmap of \(X\). Then, the iteration procedure of YTADM defined by.
is Pstable if.

i)
\(\Vert R{u}_{m}\left(x,t\right)R{u}_{n}\left(x,t\right)\Vert \le {\varepsilon }_{0}\Vert {u}_{m}\left(x,t\right){u}_{n}\left(x,t\right)\Vert\) for some \({\varepsilon }_{0}\in {R}^{+}\),

ii)
\(\Vert {A}_{m}\left(u\right){A}_{n}\left(u\right)\Vert \le {\varepsilon }_{1}\Vert {u}_{m}\left(x,t\right){u}_{n}\left(x,t\right)\Vert\) for some \({\varepsilon }_{1}\in {R}^{+}\),

iii)
\(\varepsilon =\left({\varepsilon }_{0}+{\varepsilon }_{1}\right)\Vert \frac{{t}^{\alpha }}{\Gamma \left(\alpha +1\right)}\Vert <1\).
Proof.
First we show that \(P\) has a fixed point. To do this, for \(n,m\in N\), we have.
By subtracting (15) from (14), we obtain
Taking the norm on both sides of (16), without loss of generality, we have
Using the linearity property of the YT and its inverse, we obtain.
Using the properties of the norm, we proceed.
Now, assuming that.
\(\Vert R{u}_{m}\left(x,t\right)R{u}_{n}\left(x,t\right)\Vert \le {\varepsilon }_{0}\Vert {u}_{m}\left(x,t\right){u}_{n}\left(x,t\right)\Vert\) and.
\(\Vert {A}_{m}\left(u\right){A}_{n}\left(u\right)\Vert \le {\varepsilon }_{1}\Vert {u}_{m}\left(x,t\right){u}_{n}\left(x,t\right)\Vert\) for some \({\varepsilon }_{0},{\varepsilon }_{1}\in {R}^{+}\).
Then, (17) becomes
However, from the properties of YT, we have.
Therefore, from (18), we obtain
where \(\varepsilon = \left( {\varepsilon _{0} + \varepsilon _{1} } \right)\,\left\ {\frac{{t^{\alpha } }}{{\Gamma \left( {\alpha + 1} \right)}}} \right\\).
Hence, the selfmapping \(P\) has a fixed point. Now, we show \(P\) satisfies the condition in TheoremÂ \(3.1\). For this, we have
for \(C = 0,\varepsilon = \left( {\varepsilon _{0} + \varepsilon _{1} } \right)\,\left\ {\frac{{t^{\alpha } }}{{\Gamma \left( {\alpha + 1} \right)}}} \right\{\text{ < }}\,1.\) This shows that the conditions of Theorem 3.1 hold for selfmapping \(P\). Hence, by Theorem 2.2, YTADM is Picard Tstable if \(\varepsilon <1\).
Convergence analysis of YTADM for solving NTFPDEs
In this section, we prove the convergence of YTADM when applied to NTFPDEs.
Theorem 3.2.
Let \(\left(X,\Vert .\Vert \right)\) be a Banach space and \(P: X\to X\) be a mapping associated with YTADM be defined by (13). Then, \(P\) has a unique fixed point and the sequence solution \({\left\{{u}_{n}\left(x,t\right)\right\}}_{n=0}^{\infty }\) converges to the fixed point of \(X\) with an initial value \({u}_{0}\in X\), if \(0<\varepsilon <1\) such that \(\Vert {u}_{n+1}\left(x,t\right)\Vert \le \varepsilon \Vert {u}_{n}\left(x,t\right)\Vert\).
Proof.
First, we show the existence of a fixed point of \(P\). For this define that \(\left\{{s}_{n}\right\}\) is the sequence of partial sums of the series Eq.Â (13) as,
Now, we show that \({\left\{{s}_{n}\right\}}_{n=0}^{\infty }\) is a Cauchy sequence in Banach space \(\left(X,\Vert .\Vert \right)\). For this purpose, we consider.
Now, for every \(m,n\in N\) with \(n>m\), then by using (20) and triangle inequality, we obtain
Since \(0<\varepsilon <1\), the sum \(1+\varepsilon +{\varepsilon }^{2}+\dots +{\varepsilon }^{nm}\) represents a finite geometric sequence whose total sum is \(\frac{1{\varepsilon }^{nm}}{1\varepsilon }\). Now, we have \(1{\varepsilon }^{nm}<1\); then, (21) becomes.
Since \({u}_{0}\left(x,t\right)\) is bounded,
This implies that \({s}_{n}={s}_{m}\). Therefore, \({\left\{{s}_{n}\right\}}_{n=0}^{\infty }\) is a Cauchy sequence and is convergent. Hence, \(P\) has a fixed point.
Next, to complete this proof, let \(\left\{{u}_{n}\left(x,t\right)\right\}\) converge to \(u\in X\) and we need to show that \(u\) is the unique fixed point of \(P\). To achieve this, let \(l\) be another fixed point of \(P\). Then, by (19), we have
Since,\(\left(1\varepsilon \right)<0\) for \(0<\varepsilon <1\), the above inequality can only hold if \(\Vert ul\Vert =0\) implies \(u=l\). Hence, \(u\) is the unique fixed point of \(P\), which completes the proof.
Numerical results and discussions
In this section, four examples of NTFPDEs are solved by using the YTADM. The validity of the proposed numerical method is observed in terms of its absolute errors. To show the behaviors of the corresponding examples, some figures are also plotted for fractional orders \(\alpha\). All the results are calculated using MATLAB.
Example 4.1
[61]. Consider the following nonlinear timefractional NewellWhiteheadSegel equation.
with initial conditions
The exact solution of this problem is \(u\left(x,t\right)=\frac{\frac{2}{3}\lambda {e}^{2t}}{\frac{2}{3}+\lambda \lambda {e}^{2t}}.\)for the special case \(\alpha =1\). For simplicity, we drop \(\left(x,t\right)\) from \(u\left(x,t\right)\) and \(\frac{{\partial }^{2}}{\partial {x}^{2}}u\left(x,t\right).\)
Step 1: Applying YT on both sides of (23) and using the differentiation property, we have
In our case \(0<\alpha \le 1\), (24) is reduced to
Using the given initial condition, we obtain
Step 2: Applying the inverse YT on both sides of (25), we obtain.
This implies that.
Step 3: ADM defines the solution by an infinite series of linear terms as
and the nonlinear term as
where \({A}_{n}\left(u\right)\) are the Adomian polynomials to be determined.
Step 4: substituting (28) and (27) into (26) by using the linearity property of the YT, to obtain
Step 5: Equating the terms on both sides of (29), we have the following recurrence relation.
By using (9), the few components of the Adomian polynomials for the nonlinear term \({u}^{2}\) have been derived in the form
Consequently, solving the above equations, the first few components of the YTADM solution for (30) are derived as follows:
Using the same technique,
and so on.
Step 6: The numerical solution of YTADM is given as
To verify the solution approximated by the YTADM graphically, we depict the results in Fig.Â 1 for different values of the fractional order (\(\alpha =0.4, 0.6, \text{0.8,1}\)) and in Fig.Â 2 for \(\alpha =1\).
As shown in Figs.Â 1a and b, the approximate solutions of the pattern behavior gradually decrease as the fractional order values get closer and closer to. The same as in the first example, for comparison purposes, we consider only the special case because the exact solution exists. For this reason, in Fig.Â 2, we compare the approximate solution with the exact solution.
As seen from Fig.Â 2, the approximate solution obtained is similar to the exact solution. This finding implies that the results obtained by YTADM are in good agreement with the exact results. The numerical solution to this problem is also compared with the results of the Laplace transform decomposition method (LTDM) [64], the Aboodh transform homotopy perturbation method (ATHPM) [65], and the qhomotopy Shehu analysis transform method (qHSATM) [61] in TableÂ 1 at the same number of iterations and domain discretization.
Table 1 shows the numerical approximate and comparison solutions of YTADM and other numerical methods. The absolute errors of YTADM at various points in the corresponding domain are presented and compared with LTDM [64], qHSATM [61], ATHPM [65], and the exact solution at \(\alpha =1\), it yields far more robust results than all the other numerical methods. Hence, we conclude that the absolute error determined in the table shows that the method is too accurate for treating NTFPDEs.
Example 4.2
[63] Consider the nonlinear fractional Cauchy reactionâ€“diffusion equation.
with initial condition
The exact solution of this problem is \(u(x,t)={e}^{x+t}\) for the special case \(\alpha =1\). For simplicity we drop \(\left(x,t\right)\) from each term of (32). To solve this problem by YTADM, we follow the procedure stated in Sect.Â "A hybrid numerical method for solving NTFPDEs".
Step 1: Apply YT on both sides of (32) to obtain
In our case \(0<\alpha \le 1\), (33) is reduced to
Step 2: Taking the inverse YT on both sides of (34),
Step 3: The ADM defines the solution by an infinite series of linear terms as
and the nonlinear term as
where \({A}_{n}\) and \({B}_{n}\) are the Adomian polynomials to be determined.
Step 4: Substituting (36), (37) and (38) into (35) by using the linearity property of the YT, to obtain
Step 5: To obtain the component solutions, we equating the terms on both sides of (39) to get
By using (9), the few components of the Adomian polynomials for the nonlinear terms \(u{u}_{xx}\) and \({u}^{2}\) have been derived in the form
\({A}_{3}\left(u\right)={u}_{0}{u}_{3x}+{u}_{1}{u}_{2x}+{u}_{2}{u}_{1x}+{u}_{3}{u}_{0x},\) and so on.
\({B}_{0}\left(u\right)={u}_{0}^{2},\) \({B}_{1}\left(u\right)=2{u}_{0}{u}_{1},\) \({B}_{2}\left(u\right)=2{u}_{0}{u}_{2}+{u}_{1}^{2},\) \({B}_{3}\left(u\right)=2{u}_{0}{u}_{3}+2{u}_{1}{u}_{2},\) and so on.
Since \({u}_{0}\left(x,t\right)={e}^{x}\), the first iteration solution is given by
Using the same technique, we obtain.
\({u}_{2}\left(x,t\right)=\frac{{e}^{x}{t}^{2\alpha }}{\Gamma \left(2\alpha +1\right)}\), \({u}_{3}\left(x,t\right)=\frac{{e}^{x}{t}^{3\alpha }}{\Gamma \left(3\alpha +1\right)}\), \({u}_{4}\left(x,t\right)=\frac{{e}^{x}{t}^{4\alpha }}{\Gamma \left(4\alpha +1\right)}\) and so on.
Step 6: Following component approximate solution, we obtain the general solution as
where \({\text{E}}_{\alpha } \left( . \right)\) is wellknown as the Mittagâ€“Leffler function defined in [7].
To show the solution behavior of the approximate solution obtained by YTADM, we consider different values of the fractional order \(\alpha\) \(\left(\alpha =0.4, 0.6, 0.8, 1\right)\) and depict the results in Figs.Â 3 and 4.
The numerical solutions obtained by YTADM for distinct fractional orders are depicted in Figs.Â 3a and b. We have observed from this figure that the solution behaviors of the timefractional derivatives are clearly decreasing as the values of fractional order increase. In the case of ordinary derivatives, the solution behavior of such a DE is not clearly shown. Therefore, arbitrariness in fractionalorder derivatives introduces more degrees of freedom in the design and study of realtime events. It is observed from Figs.Â 4b and d that the approximate solution obtained by YTADM is almost identical to the exact solution at. From Fig.Â 4, it is found that exact and approximate solutions are in complete agreement. The absolute errors illustrated in Figs.Â 4c and e indicate that the designed numerical method is a suitable one for solving nonlinear fractional Cauchy reactionâ€“diffusion equations that arise in various fields of science. We have also compared YTADM with the Aboodh variational iteration method (AVIM) [63] using their absolute errors. For this reason, we consider the time discretization points as in AVIM and record the obtained results in TableÂ 2.
Table 2 shows the numerical solution of Example 4.2 for different points in the domain with different values of fractional order. For the special case of fractional order \(\alpha =1\), we compared the numerical solution obtained by YTADM with the solution in the [63] and the exact solution. The obtained result indicates that the solution obtained by YTADM has a better agreement with the exact solution than the approximate solution obtained by AVIM. One of the drawbacks of an AVIM in [63] is that the author only considered the first four iterations. It is well known that to get accurate solution considering more number of iterations is advisable. For this example, we considered the first six terms of the solution components and we obtained a more accurate solution than the one found in [63] as shown in TableÂ 2.
Example 4.3
[50] Consider the following nonlinear timefractional FornbergWhitham equation:
with initial conditions
The exact solution of this problem is \(u(x,t)={e}^{\frac{x}{2}\frac{2t}{3}}\) for the special case \(\alpha =1\).
Step 1: Applying YT on both sides of (41) and using the differentiation property, we have.
In our case \(0<\alpha \le 1,\) (42) is reduced to
Step 2: Applying the inverse YT on both sides of (43), we obtain.
Step 3: ADM defines the solution by an infinite series of linear terms as
and the nonlinear term as
where \({A}_{n}, {B}_{n}\) and \({C}_{n}\) are the Adomian polynomials to be determined.
Step 4: Using the linearity property of the YT and substituting (45)(48) into (44), we obtain
Step 5: Equating the terms on both sides of (49), we have the following relation
By using (9), the few components of the Adomian polynomials for the nonlinear terms \(u{u}_{xxx}, u{u}_{x}\) and \({u}_{x}{u}_{xx}\) have been derived in the form
\({A}_{3}\left(u\right)={u}_{0}{u}_{3xxx}+{u}_{1}{u}_{2xxx}+{u}_{2}{u}_{1xxx}+{u}_{3}{u}_{0xxx},\) and so on.
\({B}_{3}\left(u\right)={u}_{0}{u}_{3x}+{u}_{1}{u}_{2x}+{u}_{2}{u}_{1x}+{u}_{3}{u}_{0x},\) and so on.
\({C}_{3}\left(u\right)={u}_{0x}{u}_{3xx}+{u}_{1x}{u}_{2xx}+{u}_{2x}{u}_{1xx}+{u}_{3x}{u}_{0xx},\) and so on.
Then the iteration solution of YTADM is obtained as follows
In a similar manner, we obtain.
\({u}_{4}\left(x,t\right)=\frac{1}{16}{e}^\frac{x}{2}\frac{{t}^{4\alpha }}{\Gamma \left(4\alpha +1\right)}\frac{3}{32}{e}^\frac{x}{2}\frac{{t}^{4\alpha 1}}{\Gamma \left(4\alpha \right)}+\frac{3}{64}{e}^\frac{x}{2}\frac{{t}^{4\alpha 2}}{\Gamma \left(4\alpha 1\right)}\frac{1}{128}{e}^\frac{x}{2}\frac{{t}^{4\alpha 3}}{\Gamma \left(4\alpha 2\right)},\) and so on.
Step 6: The numerical solution of YTADM is given as
Here we also investigated the behavior of the approximated solution by up to three iterations of YTADM given in (50) by varying the values of the timefractional order (\(\alpha =0.75, 0.85, 0.95, 1\)), and the results are presented in Fig.Â 5, and the comparison in Fig.Â 6 at \(\alpha =1\)
FigureÂ 6 shows the comparison of the numerical solution obtained by YTADM and the exact solution in the 3D and 2D plots of Example 4.3. The 3D and 2D plots have confirmed the closed contact between the YTADM and the exact solutions of Example 4.3.
In TableÂ 3 we compared the obtained results of YTADM with the results of the residual power series method (RPSM) found in [50]. The results show that, with error extremely close to zero, our technique provides a superior answer than the numerical method in [50].
Example 4.4
[62] Consider the following nonlinear timefractional Fokkerâ€“Planck equation:
with initial conditions
The exact solution of this problem is \(u(x,t)={x}^{2}{e}^{t}\) for the special case \(\alpha =1\).
Step 1: Applying YT on both sides of (51), we obtain.
In our case,\(0<\alpha \le 1\), then, (52) is reduced to.
Step 2: Applying the inverse YT on both sides of (53), to get.
Step 3: The ADM defines the solution by an infinite series of linear terms as
and the nonlinear terms as
where \({A}_{n}\) is the Adomian polynomial to be determined.
Step 4: Substituting (55) and (56) into (54), we obtain.
Step 5: Equating the terms on both sides of (57), we have the following relation.
By using (9), the few components of the Adomian polynomials for the nonlinear term \({u}^{2}\) have been derived as
Then, the iteration solution of YTADM is obtained as follows
In a similar manner, we obtain.
\({u}_{3}\left(x,t\right)={x}^{2}\frac{{t}^{3\alpha }}{\Gamma \left(3\alpha +1\right)},\) And so on.
Step 6: The numerical solution of YTADM is given as
Here, we investigate the behavior graphical representations of the approximated solution of YTADM in Fig.Â 7, 8. The approximate solution presented in Fig.Â 7a and b for various values of fractional order \(\alpha \left(\alpha =0.75, 0.85, 0.95, 1\right)\) is decreases as the fractional order \(\alpha\) increases. From the comparison in Fig.Â 8, we conclude that the obtained results of YTADM have a nice relationship with the exact solution. We have also compared the obtained results of YTADM with the results taken in [66] and the exact solution in TableÂ 4.
Table 4 shows the exact solution, the approximate solution of YTADM for fractional order \(\alpha =0.9\) and \(1\), and a comparison of the exact solution and the approximate solution obtained by HPTM, ADM, and VIM found in the reference [66] for \(\alpha =1\). The approximate solution found by YTADM is generally in considerably better agreement with the exact result than the numerical approach, in [66] as can be seen from the TableÂ 4. Additionally, it is evident from the table that we get more accurate result as \(\alpha \to 1\), which indicates that the other approximate solutions are likely to be valid.
Conclusion
In this manuscript, we investigate the numerical solution of NTFPDEs using YTADM. The procedure is understandable to the readers because it consists of the direct implementation of the YT on the portion containing the fractional derivative of the given problem to change it into the algebraic form. Finally, the ADM is applied to decompose the nonlinear portion and provides a series solution to the given problem. We offered the recommended method's stability and convergence criteria along with its proof. The theoretical explanation of the suggested strategy was supported by the presentation of four illustrative instances. The findings presented in terms of figures confirm that the results obtained by the present method are in good agreement with the exact solutions for especial case \(\alpha =1\) and the proposed approach gives better solution with error much near to zero. It is noted that the behavior of the approximate solution values at various \(\alpha\)'s is identical to that of the values obtained with the exact solution, where \(\alpha =1\). This demonstrates that the approximate solution is efficient. Tables 1, 2, 3 and 4 illustrate the findings, which indicate that the suggested YTADM approach outperforms the numerical methods that have been examined and published in the literature. Consequently, we draw the conclusion that YTADM is highly effective and potent in locating numerical solutions for a broad range of NTFPDEs. Only time fractional nonlinear PDEs in onedimensional spaces were used in the current work. This may be extended to the spaceâ€“time fractional PDEs with multidimensional spaces as well as the space fractional PDEs. Furthermore, for the purpose of addressing nonlinear realworld problems, the fractional order derivative is taken in the Caputo meanings. This may also be extended to other fractional derivatives such as CaputoFabrizio, conformable, and AtanganaBaleanu (ABC) fractional derivative orders.
Availability of data and materials
No datasets were generated or analysed during the current study.
References
Mainardi F. Fractional calculus: some basic problems in continuum and statistical mechanics. Fractals Fract Calc Contin Mech. 1997;378:291â€“348.
DelCastilloNegrete D, Carreras BA, Lynch VE. Front dynamics in reactiondiffusion systems with Levy flights: a fractional diffusion approach. Phys Rev Lett. 2003; 91(1).
Singh J, Rashidi MM, Kumar D, Swroop R. A fractional model of a dynamical Brusselator reactiondiffusion system arising in triple collision and enzymatic reactions. Nonlinear Eng. 2016;5(4):277â€“85.
Jindal N, Weber S, Andrews JG. Fractional power control for decentralized wireless networks. IEEE Trans Wirel Commun. 2008;7:5482â€“92.
Zappone A, Jorswieck E. Energy efficiency in wireless networks via fractional programming theory. Found Trends Commun Inf Theory. 2014;11:185â€“396.
Liu J, Zhang Y. Analytical study of exact solutions of the nonlinear Kortewegde Vries equation with spacetime fractional derivatives. Mod Phys Lett. 2018: 1850012.
Mamun AA, Lu C, NaharAnanna S, Uddin MM. Dynamical behavior of water wave phenomena for the 3D fractional WBBM equations using rational sineGordon expansion method. Sci Rep. 2024;14:6455.
Ferdi Y. Some applications of fractional order calculus to design digital filters for biomedical signal processing. J Mech Med Biol. 2012; 12(02).
Meena SC, Lata P, Meena M. A study of signal processing with fractional calculus. Eur Chem Bull. 2023;12(9):2858â€“69.
Jafari H, Ganji RM, Nkomo NS, Lv YP. A numerical study of fractional order population dynamics model. Results Phys. 2012;27: 104456.
Awadalla M, Noupoue YY, Abuasbeh K. Population growth modeling via RayleighCaputo fractional derivative. J Stat Appl Pro. 2021;10(1):11â€“6.
Pelap FB, Tanekou GB, Fogang CF, Kengne R. Fractionalorder stability analysis of earthquake dynamics. J Geophys Eng. 2018;15:1673â€“87.
Bisong MS, Mofor A, Fomekong F, Tsasse LC, Mando AK, Manfouo F, Tanekou GB, Motchongom MT, Afungchui D. Delayed and fractionalorder dynamics of a model for earthquake induced by stress perturbations. Phys Scr. 2024;99: 035205.
Zhukov D, Otradnov K, Kalinin V. Fractionaldifferential models of the time series evolution of sociodynamic processes with possible selforganization and memory. Mathematics. 2024;12(3):484.
Momani S, Yildirim A. Analytical approximate solutions of the fractional convectiondiffusion equation with nonlinear source term by Heâ€™s homotopy perturbation method. Int J Comput Math. 2010;87:1057â€“65.
Sakthivel R, Suganya S, Anthoni SM. Approximate controllability of fractional stochastic evolution equations. Comput Math Appl. 2012;2012(63):660â€“8.
Sadeghian H, Salarieh H, Alasty A, Meghdari A. Controllability of linear fractional stochastic systems. Sci Iran. 2015;22(1):264â€“71.
Lima LDS. Fractional stochastic differential equation approach for spreading of diseases. Entropy. 2022;24(5):719.
Luoc MTD. Existence and finitetime stability results for impulsive Caputotype fractional stochastic differential equations with time delays. Math Slovaca. 2023;73(2):387â€“406.
Dinh TN, Kamal S, Pandey RK. Fractionalorder system: control theory and applications. Fractal Fract. 2023;7(48):1â€“4.
Oprzedkiewicz K, Dziedzic K. Fractional discretecontinuous model of heat transfer process. Arch Control Sci. 2021;31(2):287â€“306.
Kumar D, Singh J, Kumar D, Singh J. Fractional calculus in medical and health science. 1st Edition, 2020.
Hilfer R. Applications of Fractional Calculus in Physics. River Edge, NJ, USA: World Scientific Publishing; 2000.
Podlubny I. Fractional Differential Equations: An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. Elsevier: Amsterdam, Netherlands. 1998; 198.
Miller KS, Ross B. An introduction to the fractional calculus and fractional differential equations, 1st ed.; Wiley: London, UK. 1993, p. 384.
Mohamed AM, Torky SM. Approximate solution of fractional nonlinear partial differential equations by the Legendre Multiwavelet Galerkin method. J. Appl. Math. 2014; vol. 2014, Article ID 192519, 12 pages.
Ogunmiloro OM. Mathematical analysis and approximate solution of a fractional order Caputo fascioliasis disease model. Chaos Solitons Fractals. 2021;2021(146): 110851.
Ogunmiloro OM. A fractional order mathematical model of teenage pregnancy problems and rehabilitation in Nigeria. Math model control. 2022;2(4):139â€“52.
Talabi AT, Odunaike RK, Odetunde OS, Omoteso KA. A modified fractional differential love model. Moroccan J Pure Appl Anal. 2019;5(1):37â€“45.
Song L, Xu S, Yang J. Dynamical models of happiness with fractional order. Commun Nonlinear Sci Numer Simul. 2010;15(3):616â€“28.
Ahmad WM, ElKhazali R. Fractionalorder dynamical models of love. Chaos Solitons Fractals. 2007;33(4):1367â€“75.
Bekela AS, Belachew MT, Wole GA. A numerical method using Laplacelike transform and variational theory for solving timefractional nonlinear partial differential equations with proportional delay. Adv Differ Equ. 2020;2020:586.
Cenesiz Y, Tasbozan O, Kurt A. Functional variable method for conformable fractional modified KdVâ€“ZK equation and Maccari system. Tbil Math J. 2017;10(1):118â€“26.
Jaradat I, Alquran M, Katatbeh Q, Yousef F, Momani S, Baleanu D. An avantgarde handling of temporalspatial fractional physical models. 2020; 21(2): 183â€“194.
Javeed S, Baleanu D, Waheed A, Khan MS, Affan H. Analysis of homotopy perturbation method for solving fractional order differential equations. Mathematics. 2019; 7(40).
Singh J, Kumar D, KiliÃ§man A. Numerical solutions of nonlinear fractional partial differential equations arising in spatial diffusion of biological populations. Abstr. Appl. Anal. 2014; vol. 2014, Article ID 535793, 12 pages.
Wang K, Liu S. A new Sumudu transform iterative method for timefractional Cauchy reactionâ€“diffusion equation. Springer Plus. 2016:2â€“20.
Ziane D, Elzaki TM, Cherif MH. Elzaki transform combined with variational iteration method for partial differential equations of fractional order. Fundam J Math Appl. 2018;1(1):102â€“8.
Dehestani H, Ordokhani Y, Razzaghi M. Fractionalorder LegendreLaguerre functions and their applications in fractional partial differential equations. Appl Math Comput Journal. 2018;336:433â€“53.
Dehestani H, Ordokhani Y, Razzaghi M. Application of the modified operational matrices in multiterm variableorder timefractional partial differential equations. Math Meth Appl Sci. 2019;42:7296â€“313.
Dehestani H, Ordokhani Y, Razzaghi M. Fractionalorder Genocchiâ€“Petrovâ€“Galerkin method for solving timeâ€“space fractional FokkerPlanck equations arising from the physical phenomenon. Int J Appl Comput Math. 2020;6:100.
Wang L, Wu Y, Ren Y, Chen X. Two analytical methods for fractional partial differential equations with proportional delay. IAENG Int J Appl Math. 2019;49(1):1â€“6.
Malyk I, Mykola G, Chaudhary A, Sharma S, Dubey RS. Numerical solution of nonlinear fractional diffusion equation in the framework of the Yangâ€“Abdelâ€“Cattani derivative operator. Fractal Fract. 2021;5:64.
Akram G, Abbas M, Tariq H, Sadaf M, Abdeljawad T, Alqudah MA. Numerical approximations for the solutions of fourth order time fractional evolution problems using a novel spline technique. Fractal Fract. 2022;6:170.
AlDeiakeh R, Ali M, Alquran MA, Sulaiman T, Momani S, AlSmadi M. On finding closedform solutions to some nonlinear fractional systems via the combination of multiLaplace transform and the Adomian decomposition method. Rom Rep Phys. 2022; 74(111).
Mallick B, Sahu PK, Routaray M. An iterative method for solving timefractional partial differential equations with proportional delays. TWMS J App Eng Math. 2023;13(1):74â€“85.
Kumar M. An iterative approach for solving fractional order Cauchy reactiondiffusion equations. J Appl Math Comput Mech. 2023;22(3):19â€“32.
Brikaa MG. An analytic algorithm for the spacetime fractional reactiondiffusion equation. J Interpolat Approx Sci Comput. 2015;2015(2):112â€“27.
Nadeem M, He J, He C, Sedighi HM, Shirazi AH. A numerical solution of nonlinear fractional NewellWhiteheadSegel equation using natural transform. J Pure Appl Math. 2022;13(2):168â€“82.
Khan H, Kumam P, Nawaz A, Khan Q, Khan S. The fractional investigation of FornbergWhitham equation using an efficient technique. Comp Mod Eng Sci. 2023;135(1):260â€“73.
Firoozjaee MA, Youse SA, Jafari H. A numerical approach to FokkerPlanck equation with space and timefractional and nonfractional derivatives. MATCH Commun Math Comput Chem. 2015;2015(74):449â€“64.
Alkhaled K. Numerical solution of timefractional partial differential equations using Sumudu decomposition method. Rom.6Journ.oPhys. 2015; 60(1â€“2): 99â€“110.
Adomian G. A review of the decomposition method in applied mathematics. J Math Anal Appl. 1988;135(2):501â€“44.
Wazwaz AM, Partial differential equations and solitary waveâ€™s theory. Springer Science & Business Media. 2010.
Akinola E, Oladejo J, Akinpelu F, Owolabi J. On the application of Sumudu transform series decomposition method and oscillation equations. 2017.
Yang XJ. A new integral transform method for solving steady heattransfer problem. Therm Sci. 2016;20(3):639â€“42.
Bhuvaneswari R, Bhuvaneswari K. Application of yang transform in cryptography. Int J Eng Sci Math. 2020;9(3):2320â€“0294.
Aruldoss R, Jasmine G. Numerical solutions of time fractional nonlinear partial differential equations using Yang transform combined with variational iteration method. Glob J Pure Appl Math. 2020;16(2):249â€“60.
Dattu KU. New integral transform: fundamental properties, investigations and applications. IAETSD J adv res appl sci. 2018;5(4):534â€“9.
Khan H, Khan A, Chen W, Shah K. Stability analysis and a numerical scheme for fractional kleingordon equations. Math Methods Appl Sci. 2019;42(2):723â€“32.
Bektas U, Anac H. A hybrid method to solve a fractionalorder Newellâ€“Whiteheadâ€“Segel equation. Bound Value Probl. 2024;2024(38):1â€“13.
Prakash A, Kumar, M. Numerical solution of timefractional order FokkerPlanck equation. 2019; 9(3): 446454.
Khalouta A. Exact solution of nonlinear timefractional reactiondiffusionconvection equation via a new coupling method. Int J Nonlinear Anal Appl. 2022;13(2):333â€“44.
Ayata M, Ozkan O. A new application of conformable Laplace decomposition method for fractional NewellWhiteheadSegel equation. AIMS Math. 2020;5(6):7402â€“12.
Jani HP, Singh TR. Aboodh transform homotopy perturbation method for solving fractionalorder NewellWhiteheadSegel equation. Math Methods Appl Sci. 2022. https://doi.org/10.1002/mma.8886.
Saravanan A, Magesh N. An efficient computational technique for solving the FokkerPlanck equation with space and time fractional derivatives. J King Saud Univ Sci. 2016;28:160â€“6.
Acknowledgements
Not applicable.
Funding
There is no any fund support in this research work.
Author information
Authors and Affiliations
Contributions
ATD formulated the problem, designed, and drafted the manuscript. ASB carried out scheme development, numerical experimentation and produced tables and figures, and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bekela, A.S., Deresse, A.T. A hybrid yang transform adomian decomposition method for solving timefractional nonlinear partial differential equation. BMC Res Notes 17, 226 (2024). https://doi.org/10.1186/s13104024068777
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104024068777