# Novel Investigation of Fractional-Order Cauchy-Reaction Diffusion Equation Involving Caputo-Fabrizio Operator

Authors: Meshari Alesemi; Naveed Iqbal; Mohammed S. Abdo
Journal: Journal of Function Spaces
 (2022)

---

## Abstract

In this article, the new iterative transform technique and homotopy perturbation transform method are applied to calculate the fractional-order Cauchy-reaction diffusion equation solution. Yang transformation is mixed with the new iteration method and homotopy perturbation method in these methods. The fractional derivative is considered in the sense of Caputo-Fabrizio operator. The convection-diffusion models arise in physical phenomena in which energy, particles, or other physical properties are transferred within a physical process via two processes: diffusion and convection. Four problems are evaluated to demonstrate, show, and verify the present methods’ efficiency. The analytically obtained results by the present method suggest that the method is accurate and simple to implement.


---

## Body

## 1. Introduction

The convection-diffusion equation is a mixture of convection and diffusion equations and identifies physical processes where energy, particles, or other physical properties are transmitted inside a physical process due to two process steps: diffusion and convection. In standard form, the convection-diffusion model is written as follows:(1)∂U∂I=∇·D·∇U−∇·v⟶U+R,where D is the diffusivity, U is the variable term, such as thermal diffusivity for heat flow or mass diffusion coefficient for particle, and v⟶ is the average velocity that the volume is travelling. For instance, in convection, u might be the density of river in salt and then the flow velocity of water v⟶. For example, in a calm lake, v⟶ would be the average velocity of bubbles rising to the surface due to buoyancy, and U would be the concentration of small bubbles. R defines “sinks” or “sources” of the quantity U. For a chemical species, R>0 indicates that a chemical reaction is increasing the number of the species, while R>0 indicates that a chemical reaction is decreasing the number of the species. If thermal energy is generated by friction, R>0 may occur in heat transport. ∇ denotes gradient, while ∇· denotes divergence. Previously, different techniques have been applied to investigate these models such as Adomian’s decomposition technique [1], variational iteration technique [2], Bessel collocation technique [3], and homotopy perturbation technique [4].In recent decades, fractional derivatives have been used to interpret many physical problems mathematically, and these representations have produced excellent results in modelling real-world issues. Many basic definitions of fractional operators were given by Riesz, Riemann-Liouville, Hadamard, Weyl, Grunwald-Letnikov, Liouville-Caputo, Caputo-Fabrizio, and Atangana-Baleanu, among others [5–8]. Over the last few years, many nonlinear equations have been developed and widely used in nonlinear physical sciences like chemistry, biology, mathematics, and different branches of physics like plasma physics, condensed matter physics, fluid mechanics, field theory, and nonlinear optics. The exact outcome of nonlinear equations is crucial in determining the characteristics and behaviour of physical processes. Still, it is impossible to find exact results when dealing with linear equations. Many useful methods have been applied to investigate nonlinear fractional partial differential equations, for example, analytical solutions with the help of natural decomposition method of fractional-order heat and wave equations [9], fractional-order partial differential equations with proportional delay [10], fractional-order hyperbolic telegraph equation [11] and fractional-order diffusion equations [12], the variational iterative transform method [13], the homotopy perturbation transform method [14, 15], the homotopy analysis transform method [16, 17], reduced differential transform method [18, 19], q-homotopy analysis transform method [20–24], the finite element technique [25], the finite difference technique [26], and so on [27–30].Daftardar-Gejji and Jafari developed a new iterative method of analysis for solving nonlinear equations in 2006 [31, 32]. It is the first application of Laplace transformation in iterative technique by Jafari et al. Iterative Laplace transformation method [33] was introduced as a simple method for estimating approximate effects of the fractional partial differential equation system. Iterative Laplace transformation method (NITM) is used to solve linear and nonlinear partial differential equations such as fractional-order Fornberg Whitham equations [34], time-fractional Zakharov Kuznetsov equation [35], and fractional-order Fokker Planck equations [36].In 1999, He developed the homotopy perturbation method (HPM) [37], which combines the homotopy technique, and the standard perturbation method has been broadly utilized to both linear and nonlinear models [38–40]. The homotopy perturbation method is important because it eliminates the need for a small parameter in the model, eliminating the disadvantages of traditional perturbation techniques. The main goal of this paper is to use HPM to solve nonlinear fractional-order Cauchy-reaction diffusion equation using a newly introduced integral transformation known as the “Yang transform” [41]. The suggested technique is applied to analyse two well-known nonlinear partial differential equations. In the context of a quickly convergent series, we obtain a power series solution, and only a few iterations are required to obtain very efficient solutions. There is no need for a discretization technique or linearization for the nonlinear equations, and just a few few can yield a result that can be quickly estimated to utilize these methods.


## 2. Basic Definitions

We provide the fundamental definitions that will be used throughout the article. For the purpose of simplification, we write the exponential decay kernel as,KI,ϱ=e−℘I−ϱ/1−℘.Definition 1.
The Caputo-Fabrizio derivative is given as follows [42]:
(2)CFDI℘ℙI=N℘1−℘∫0Iℙ′ϱKI,ϱdϱ,n−1<℘≤n.
N℘ is the normalization function with N0=N1=1.
(3)CFDI℘ℙI=N℘1−℘∫0IℙI−ℙϱKI,ϱdϱ.Definition 2.
The fractional integral Caputo-Fabrizio is given as [42]
(4)CFII℘ℙI=1−℘N℘ℙI+℘N℘∫0Iℙϱdϱ,I≥0,℘∈0,1.Definition 3.
ForN℘=1, the following result shows the Caputo-Fabrizio derivative of Laplace transformation [42]:
(5)LCFDI℘ℙI=vLℙI−ℙ0v+℘1−v.Definition 4.
The Yang transformation ofℙI is expressed as [42]
(6)YℙI=χv=∫0∞ℙIe−IvdI.I>0,Remark 5.
Yang transformation of few useful functions is defined as below.(7)Y1=v,YI=v2,YIi=Γi+1vi+1.Lemma 6 (Laplace-Yang duality).
Let the Laplace transformation ofℙI be Fv, and then, χv=F1/v [43].Proof.
From Equation (5), we can achieve another type of the Yang transformation by putting I/v=ζ as
(8)LℙI=χv=v∫0∞ℙvζeζdζ.ζ>0.
SinceLℙI=Fv, this implies that
(9)Fv=LℙI=∫0∞ℙIe−vIdI.
PutI=ζ/v in (8), and we have
(10)Fv=1v∫0∞ℙζveζdζ.
Thus, from Equation (7), we achieve
(11)Fv=χ1v.
Also from Equations (5) and (8), we achieve
(12)F1v=χv.
The connections (10) and (11) represent the duality link between the Laplace and Yang transformation.Lemma 7.
LetℙI be a continuous function; then, the Caputo-Fabrizio derivative Yang transformation of ℙI is define by [43]
(13)YℙI=YℙI−vℙ01+℘v−1.Proof.
The Caputo-Fabrizio fractional Laplace transformation is given by(14)LℙI=LvℙI−ℙ0v+℘1−v.
Also, we have that the connection among Laplace and Yang property, i.e.,χv=F1/v. To achieve the necessary result, we substitute v by 1/v in Equation (13), and we get
(15)YℙI=1/vYℙI−ℙ01/v+℘1−1/v,YℙI=YℙI−vℙ01+℘v−1.
The proof is completed.


## 3. Algorithm of the HPTM

The procedure of general nonlinear Caputo-Fabrizio fractional partial differential equations is through HPTM. Let us take a general nonlinear Caputo-Fabrizio partial differential equations with nonlinear functionNUφ,I and linear fractional LUφ,I as [43]
(16)CFDIρVφ,I+LVφ,I+NVφ,I=gφ,I,Vφ,0=hφ,where the term gφ,I shows the source function. Using Yang transformation to Equation (16), one can obtain
(17)YVφ,I−vVφ,01+℘v−1=−YLVφ,I+NVφ,I+Ygφ,I,YVφ,I=vhφ−1+℘v−1YLVφ,I+NVφ,I+Ygφ,I.Implementing inverse Yang transformation, we obtain(18)Vφ,I=Vφ,0−Y−11+℘v−1YLVφ,I+NVφ,I+Ygφ,I,where the term Vφ,I shows the source function and with the initial condition. Now, we apply homoptopy perturbation method.
(19)Vφ,I=∑i=0∞ρiViφ,I.We decompose the nonlinear termNVφ,I as
(20)NVφ,I=∑i=0∞ρiHiV,where HiV represents the He’s polynomial and is calculated through the following formula:
(21)HiV1,V2,V3,⋯,Vi=1Γi+1∂i∂ρiN∑i=0∞ρiViρ=0,i=1,2,3.Substituting Equations (19) and (20) in Equation (18), we obtain
(22)∑i=0∞ρiViφ,I=Vφ,I−ρY−11+℘v−1YL∑i=0∞ρiiVφ,I+N∑i=0∞ρiHiV.We obtain the following terms by coefficients comparing ofρ in (22):
(23)ρ0:V0φ,I=Vφ,I,ρ1:V1φ,I=Y−11+℘v−1YLV0φ,I+H0V,ρ2:V2φ,I=Y−11+℘v−1YLV1φ,I+H1V,ρ3:V3φ,I=Y−11+℘v−1YLV2φ,I+H2V,⋮ρi:Viφ,I=Y−11+℘v−1YLViφ,I+HiV.As a result, the obtained solution of Equation (16) can be written as follows:
(24)Vφ,I=V0φ,I+V1φ,I+⋯.


## 4. Error Analysis and Convergence

The following theorems are fundamental on the techniques address the original models [16] error analysis and convergence.Theorem 8.
LetVφ,I be the actual result of (16), and let Viφ,I∈H and σ∈0,1, where H denotes the Hilbert space. Then, the achieved result ∑i=0∞Viφ,I will converge Vφ,I if Viφ,I≤Vi−1φ,I∀i>A, i.e., for any ω>0∃A>0, such that Vi+nφ,I≤β,∀i,n∈N [43].Proof.
We make a sequence of∑i=0∞Viφ,I.(25)C0φ,I=V0φ,I,C1φ,I=V0φ,I+V1φ,I,C2φ,I=V0φ,I+V1φ,I+V2φ,I,C3φ,I=V0φ,I+V1φ,I+V2φ,I+V3φ,I,⋮Ciφ,I=V0φ,I+V1φ,I+V2φ,I+⋯+Viφ,I.
To provide the correct outcome, we have to demonstrate thatCiφ,I forms a “Cauchy sequence.” Take, for example,
(26)Ci+1φ,I−Ciφ,I=Vi+1φ,I≤σViφ,I≤σ2Vi−1φ,I≤σ3Vi−2φ,I⋯≤σi+1V0φ,I.
Fori,n∈N, we acquire
(27)Ciφ,I−Cnφ,I=Vi+nφ,I=Ciφ,I−Ci−1φ,I+Ci−1φ,I−Ci−2φ,I+Ci−2φ,I−Ci−3φ,I+⋯+Cn+1φ,I−Cnφ,I≤Ciφ,I−Ci−1φ,I+Ci−1φ,I−Ci−2φ,I+Ci−2φ,I−Ci−3φ,I+⋯+Cn+1φ,I−Cnφ,I≤σiV0φ,I+σi−1V0φ,I+⋯+σi+1V0φ,I=V0φ,Iσi+σi−1+σi+1=V0φ,I1−σi−n1−σi+1σn+1.
Since0<σ<1, and V0φ,I is bounded, let us take β=1−σ/1−σi−nσn+1V0φ,I. Thus, Viφ,Ii=0∞ forms a “Cauchy sequence” in H. It follows that the sequence Viφ,Ii=0∞ is a convergent sequence with the limit limi⟶∞Viφ,I=Vφ,I for ∃Vφ,I∈H. Hence, this ends the proof.Theorem 9.
Let∑h=0kVhφ,I is finite and Vφ,I represents the obtained series solution. Let σ>0 such that Vh+1φ,I≤Vhφ,I; then, the following relation gives the maximum absolute error [43].
(28)Vφ,I−∑h=0kVhφ,I<σk+11−σV0φ,I.Proof.
Since∑h=0kVhφ,I is finite, this implies that ∑h=0kVhφ,I<∞.
Consider(29)Vφ,I−∑h=0kVhφ,I=∑h=k+1∞Vhφ,I≤∑h=k+1∞Vhφ,I≤∑h=k+1∞σhV0φ,I≤σk+11+σ+σ2+⋯V0φ,I≤σk+11−σV0φ,I.
This ends the theorem’s proof.


## 5. The General Procedure of NITM

The general solution of fractional-order partial differential equation is as follows:(30)CFDI℘Vφ,I+NVφ,I+MVφ,I=hφ,I,i∈N,i−1<℘≤i,where N is nonlinear and M linear functions.With the initial condition(31)Vkφ,0=gkφ,k=0,1,2,⋯,i−1,implementing the Yang transformation of Equation (30), we get
(32)YDI℘Vφ,I+YNVφ,I+MVφ,I=Yhφ,I.Applying the Yang differentiation is given to(33)YVφ,I=vVφ,0+1+℘v−1Yhφ,I−1+℘v−1YNVφ,I+MVφ,I.Using inverse Yang transformation Equation (32), we get
(34)Vφ,I=Y−1vVφ,0+1+℘v−1Yhφ,I−Y−11+℘v−1YNVφ,I+MVφ,I.By iterative method, we get(35)Vφ,I=∑i=0∞Viφ,I,(36)N∑i=0∞Viφ,I=∑i=0∞NViφ,I.The nonlinear termN is identified as
(37)N∑i=0∞Viφ,I=V0φ,I+N∑i=0∞Viφ,I−M∑i=0∞Viφ,I.Putting Equations (35)–(37) in Equation (34), we have obtain the following solution:
(38)∑i=0∞Viφ,I=Y−11+℘v−1∑i=0∞s2−φ+iuiφ,0+Yhφ,I−Y−11+℘v−1YN∑i=0∞Viφ,I−M∑i=0∞Viφ,I,V0φ,I=Y−1vVφ,0+1+℘v−1Ygφ,I,V1φ,I=−Y−11+℘v−1YNV0φ,I+MV0φ,I,Vm+1φ,I=−Y−11+℘v−1Y−N∑i=0iViφ,I−M∑i=0iViφ,I,m≥1.Lastly, Equations (30) and (31) provide the i-term solution in series form which is expressed as
(39)Vφ,I≅V0φ,I+V1φ,I+V2φ,I+⋯..,+Viφ,I,i=1,2,⋯.Example 10.
Consider fractional-order Cauchy-reaction diffusion equation as [44]
(40)CFDI℘Vφ,I=DI2Vφ,I−Vφ,I,0<℘≤1,φ,I∈Ω⊂R2,with initial and boundary conditions
(41)Vφ,0=e−φ+φ=gφ,V0,I=1=f0I,∂V0,I∂I=e−I−1=f1I,φ,I∈R.The methodology consists of applying Yang transformation first on both side in (40) and utilizing the differentiation property of Yang transformation, and we have
(42)YVφ,I=ve−φ+φ+1+℘v−1YDI2V−V.Using Yang inverse transform, we get(43)Vφ,I=e−φ+φ+Y−11+℘v−1YDI2V−V.Now, we apply the new iterative transform method(44)V0φ,I=e−φ+φ,V1φ,I=Y−11+℘v−1YDI2V0−V0=−φ1+℘I−℘,V2φ,I=Y−11+℘v−1YDI2V1−V1=φ1−℘2℘I+1−℘2+℘2I22,V3φ,I=Y−11+℘v−1YDI2V2−V2=−φ1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!⋮The series type solution is given as(45)Vφ,I=V0φ,I+V1φ,I+V2φ,I+V3φ,I+⋯Viφ,I.The approximate solution is achieved as(46)Vφ,I=e−φ+φ1−1+℘I−℘+1−℘2℘I+1−℘2+℘2I22−1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯.Now applying the HPTM, we get(47)∑i=0∞piViφ,I=e−φ+φ+pY−11+℘v−1Y∑i=0∞piHiV,where the polynomials represent the nonlinear functions are HiV. For instance, the terms of He’s polynomials are achieved through the recursive relationship HiV=DI2Vi−Vi, ∀n∈N. Now, as the correspond power coefficients of p is comparison on both sides, the following solution is obtained as follows:
(48)p0:V0φ,I=e−φ+φ,p1:V1φ,I=Y−11+℘v−1YH0V=−φ1+℘I−℘,p2:V2φ,I=Y−11+℘v−1YH1V=φ1−℘2℘I+1−℘2+℘2I22,p3:V3φ,I=Y−11+℘v−1YH2V=−φ1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!,⋮Then, the homotopy perturbation method series form solution is defined as(49)Vφ,I=∑i=0∞piViφ,I.The analytical result of the above equation is defined as(50)Vφ,I=e−φ+φ1−1+℘I−℘+1−℘2℘I+1−℘2+℘2I22−1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯=φ∑i=0∞I℘iΓi℘+1,Vφ,I=e−φ+φE℘I℘.The exact result of the above equation is(51)Vφ,I=e−φ+φe−I.Figure1 shows the analytical solution of two methods at different fractional-order ℘=1 and 0.8, and Figure 2 shows separate fractional-order at ℘=0.6 and 0.4 with close contact with each other. In Figure 3, the graph shows the different fractional-order ℘ of Example 10.Example 11.
Consider fractional-order Cauchy-reaction diffusion equation as [44]
(52)CFDI℘Vφ,I=DI2Vφ,I−1+4φ2Vφ,I,0<℘≤1,φ,t∈Ω⊂R2,with initial condition
(53)Vφ,0=eφ2.and the exact result is given as
(54)Vφ,I=eφ2+1.Figure 1
(a)℘=1 and (b) the fractional-order ℘=0.8 of Example 10.Figure 2
Different fractional-order of℘=0.6 and 0.4 of Example 10.Figure 3
The different fractional-order℘ of Example 10.Now, we apply the new iterative transform method(55)V0φ,I=eφ2,V1φ,I=Y−11+℘v−1YDI2V0φ,I−1+4φ2V0φ,I=eφ21+℘I−℘,V2φ,I=Y−11+℘v−1YDI2V1φ,I−1+4φ2V1φ,I=eφ21−℘2℘I+1−℘2+℘2I22,V3φ,I=Y−11+℘v−1YDI2V2φ,I−1+4φ2V2φ,I=eφ21−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!,⋮The series type solution is given as(56)Vφ,I=V0φ,I+V1φ,I+V2φ,I+V3φ,I+⋯Viφ,I.The approximate solution of the above equation is defined as(57)Vφ,I=eφ21+1+℘I−℘+1−℘2℘I+1−℘2+℘2I22+1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯,Vφ,I=eφ2E℘I℘.Now by applying homotopy perturbation transform method, we get(58)∑i=0∞piViφ,I=eφ2+pY−11+℘v−1Y∑i=0∞piHiw.Comparing the coefficients of powerp, we get
(59)p0:V0φ,I=eφ2,p1:V1φ,I=Y−11+℘v−1YH0w=eφ21+℘I−℘,p2:V2φ,I=Y−11+℘v−1YH1w=eφ21−℘2℘I+1−℘2+℘2I22,p3:V3φ,I=Y−11+℘v−1YH2w=eφ21−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!,⋮The HPTM series solution is given as(60)Vφ,I=∑i=0∞piViφ,I,Vφ,I=eφ21+1+℘I−℘+1−℘2℘I+1−℘2+℘2I22+1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯,Vφ,I=eφ2E℘I℘.Now℘=1; then, the actual result of Equation (52) is Vφ,I=eφ2+I.Figure4 shows the analytical solution of two methods at different fractional-order ℘=1 and 0.8, and Figure 5 shows the separate fractional-order at ℘=0.6 and 0.4 with close contact with each other. In Figure 6, the graph shows the different fractional-order ℘ of Example 11.Example 12.
Consider fractional-order Cauchy-reaction diffusion equation [44]
(61)CFDI℘Vφ,I=DI2Vφ,I+2IVφ,I,0<℘≤1,φ,I∈Ω⊂R2,with initial condition
(62)Vφ,0=eφ.
The exact result is(63)Vφ,I=eφ+I+I2.
By using the Yang transformation, we get(64)Vφ,I=eφ+Y−11+℘v−1YDI2Vφ,I+2IVφ,I.Figure 4
(a)℘=1 and (b) the fractional-order ℘=0.8 of Example 10.Figure 5
The different fractional-order of℘=0.6 and 0.4 of Example 10.Figure 6
The different fractional-order℘ of Example 10.Now, we apply the new iterative transform method(65)V0φ,I=eφ,V1φ,I=Y−11+℘v−1YDI2V0φ,I+2IV0φ,I=eφ1+℘I−℘+1−℘2℘I+1−℘2+℘2I22,V2φ,I=Y−11+℘v−1YDI2V1φ,I+2IV1φ,I=eφ1−℘2℘I+1−℘2+℘2I22+1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!,⋮The series type solution is given as(66)Vφ,I=V0φ,I+V1φ,I+V2φ,I+V3φ,I+⋯Viφ,I.The approximate solution of the above equation is defined as(67)Vφ,I=eφ+eφ1+℘I−℘+1−℘2℘I+1−℘2+℘2I22+eφ1−℘2℘I+1−℘2+℘2I22+1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+.Now, using HPM, we get(68)∑i=0∞piViφ,I=eφ+pY−11+℘v−1Y∑i=0∞piHiw.Comparing the coefficients of powerp, we get
(69)p0:V0φ,I=eφ,p1:V1φ,I=Y−11+℘v−1YH0w=eφ1+℘I−℘+1−℘2℘I+1−℘2+℘2I22,p2:V2φ,I=Y−11+℘v−1YH1w=eφ1−℘2℘I+1−℘2+℘2I22+1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!.Proceeding in this path, the rest of theVnφ,I for n≥3 component can be completely recovered and the series solution can therefore be absolutely determined. Eventually, we approximate the numerical solution Vφ,I to the truncated series.
(70)Vφ,I=limN⟶∞∑n=1NViφ,I,Vφ,I=eφ+eφ1+℘I−℘+1−℘2℘I+1−℘2+℘2I22+eφ1−℘2℘I+1−℘2+℘2I22+1−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯.Now for℘=1, the closed form of the above series is
(71)Vφ,I=eφ+I+I2.Figure7 shows the analytical solution of two methods at different fractional-order ℘=1 and 0.8, and Figure 8 shows the separate fractional-order at ℘=0.6 and 0.4 with close contact with each other. In Figure 9, the graph shows the different fractional-order ℘ of Example 12.Example 13.
Consider fractional-order Cauchy-reaction diffusion equation as [44]
(72)CFDI℘Vφ,I=DI2Vφ,I−4φ2−2I+2Vφ,I,0<℘≤1,φ,I∈Ω⊂R2,with initial condition
(73)Vφ,0=eφ2.The exact result is
(74)Vφ,I=eφ2+I2.Figure 7
(a)℘=1 and (b) the fractional-order ℘=0.8 of Example 10.Figure 8
The different fractional-order of℘=0.6 and 0.4 of Example 10.Figure 9
The different fractional-order℘ of Example 10.Now, we apply the new iterative transform method(75)V0φ,I=eφ2,V1φ,I=Y−11+℘v−1YDI2V0φ,I−4φ2−2I+2V0φ,I=eφ21+℘I−℘,V2φ,I=Y−11+℘v−1YDI2V1φ,I−4φ2−2I+2V1φ,I=eφ21−℘2℘I+1−℘2+℘2I22,V3φ,I=Y−11+℘v−1YDI2V2φ,I−4φ2−2I+2V2φ,I=eφ21−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!,⋮The series type solution is given as(76)Vφ,I=V0φ,I+V1φ,I+V2φ,I+V3φ,I+⋯Viφ,I.The approximate solution of the above example is(77)Vφ,I=eφ2+eφ21+℘I−℘+eφ21−℘2℘I+1−℘2+℘2I22+eφ21−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯.Now, using the HPM, we get(78)∑i=0∞piViφ,I=eφ2+pY−11+℘v−1Y∑i=0∞piHiw.Comparing the coefficients of powerp, we have
(79)p0:V0φ,I=eφ2,p1:V1φ,I=Y−11+℘v−1YH0w=eφ21+℘I−℘,p2:V2φ,I=Y−11+℘v−1YH1w=eφ21−℘2℘I+1−℘2+℘2I22,p3:V3φ,I=Y−11+℘v−1YH2w=eφ21−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!.Similarly, the remainder of theViφ,I components for n≥4 can be completely achieved, thereby fully evaluating the series solutions. Finally, we estimate the approximate result Vφ,I by truncated sequence
(80)Vφ,I=limN⟶∞∑n=1NViφ,I,Vφ,I=eφ2+eφ21+℘I−℘+eφ21−℘2℘I+1−℘2+℘2I22+eφ21−℘23℘I+1−℘3+3℘21−℘I22+℘3I33!+⋯.Figure10 shows the analytical solution of two methods at different fractional-order ℘=1, 0.8, 0.6, and 0.4 of Example 13. The special case for ℘=1, and the above problem close form is given as
(81)Vφ,I=eφ2+I2.Figure 10
The different fractional-order℘ of Example 13.


## 6. Conclusion

The homotopy perturbation transform technique and the Iterative transform method are used in this article to obtain numerical solutions for the fractional-order Cauchy-reaction diffusion equation, which is broadly used in applied sciences as a problem for spatial effects. In physical models, the techniques produce a series of form results that converge quickly. The obtained results in this article are expected to be useful for further analysis of complicated nonlinear physical problems. The calculations for these techniques are very simple and straightforward. As a result, we can conclude that these techniques can be used to solve a variety of nonlinear fractional-order partial differential equation schemes.



---

*Source: 4284060-2022-01-27.xml*