# Analytic Solution for the Strongly Nonlinear Multi-Order Fractional Version of a BVP Occurring in Chemical Reactor Theory

Authors: Vedat Suat Erturk; A.K. Alomari; Pushpendra Kumar; Marina Murillo-Arcila
Journal: Discrete Dynamics in Nature and Society
 (2022)

---

## Abstract

This study is devoted to constructing an approximate analytic solution of the fractional form of a strongly nonlinear boundary value problem with multi-fractional derivatives that comes in chemical reactor theory. We construct the solution algorithm based on the generalized differential transform technique in four simple steps. The fractional derivative is defined in the sense of Caputo. We also mathematically prove the convergence of the algorithm. The applicability and effectiveness of the given scheme are justified by simulating the equation for given parameter values presented in the system and compared with existing published results in the case of standard derivatives. In addition, residual error computation is used to check the algorithm’s correctness. The results are presented in several tables and figures. The goal of this study is to justify the effects and importance of the proposed fractional derivative on the given nonlinear problem. The generalization of the adopted integer-order problem into a fractional-order sense which includes the memory in the system is the main novelty of this research.


---

## Body

## 1. Introduction

Chemical reactors are containers used in chemical engineering to contain chemical processes. Because of their numerous industrial uses, these reactors are crucial. Biological treatment, algae production, and gasoline production are some of the applications for tubular reactors. The mathematical model for an adiabatic tubular chemical reactor that performs an irreversible exothermic chemical reaction is examined in this work. The model may be simplified into the following nonlinear ordinary differential equation for steady-state solutions [1]:(1)d2udx2−λdudx+λμβ−uexpu=0,where λ,μ, and β are the Péclet number, Damköhler number, and adiabatic temperature rise, respectively. The relative boundary conditions are given by(2)u′0=λu0,u′1=0.In references [2, 3], the authors investigated the existence of a solution to equations (1) and (2). In [1], the researchers established the existence of numerous solutions. To solve the problem under specific evaluations, certain numerical approaches were used. Green’s function, for example, is used to turn the issue into a Hammerstein integral equation in [4]. After that, Adomian’s decomposition approach was used to solve the resultant equation. The problem was solved using the Chebyshev finite difference approach in [5]. The authors in [6] used a solution aligned on embedding Green’s function inside Krasnoselskii–Mann fixed point iteration method to solve the problem.In this article, we incorporate fractional order into (1). The following noninteger-order differential equation describes the new equation:(3)dqudxq−λdγudxγ+λμβ−uexpu=0,where dq/dxq and dγ/dxγ are the Caputo derivative operators along with fractional orders q∈1,2, γ∈0,1 subject to the boundary conditions. (1) can be considered as a particular form of (3) by fixing the orders q=2 and γ=1. This means the final solution of the fractional-order system must converge to the solution of the integer-order counterpart of the equation. Some recent results on boundary value problems in fractional-order sense can be seen in references [7, 8]. Mostly, the nature of complex dynamics cannot be better stated by integer-order differential equations. In the present case, the fractional-order model may strongly define the nature of the given system. In some circumstances, fractional models yield superior approximation results, according to Abbas et al. [9]. Iyiolaa et al. [10] have also demonstrated that for cancer tumours, the fractional model delivers a better approximation outcome than the integer-order one. Some recent studies related to the modeling in terms of fractional-order boundary value problems can be learned from references [11, 12]. Recently, a number of nonlinear fractional-order models have been proposed by the researchers to describe the dynamics of various real-life problems like AH1N1/09 influenza [13], childhood diseases [14], human liver dynamics [15], greenhouse gas effects on the population of aquatic animals [16], mosaic diseases in plants [17], maize streak virus [18], and so on.Fractional calculus has more than three century history and has progressed steadily to the present day. Rieman and Liouville were the first ones who defined the fractional-order differentiation notion in the nineteenth century. Fractional differential equations (FDEs) have been shown to be a valuable tool for representing a wide range of scientific and engineering phenomena. Many FDEs that describe any phenomenon have lack of analytic solutions. As a result of the absence of analytic solutions, a significant variety of techniques for solving FDEs have been devised [19]. FDEs have attracted much more attention as a part of fractional calculus. It is worth noting that a general solution strategy for fractional differential equations is yet to be developed. The majority of problem-solving strategies in this field have been created for certain categories of challenges. For this reason, a single standard technique for solving problems related to fractional calculus has not been found. As a result, identifying compelling and beneficial solution strategies in combination with quick application techniques is valuable and worthy of further investigation [20] (see [21–27] for further information).Under the best of our investigations, our study introduces the firstly produced numerical solution of (3). For this target, we are directed to find the approximate solution of (3) via generalized differential transform method (GDTM) [28–31].This paper is organized as follows. In Section2, a review of the GDTM [28] is given and some important preliminaries are given. In Section 3, the solution procedure is presented. Convergence theorem of the present solution is proved in Section 4. The solution approximations for equations (2) AND (3) are established in Section 5. Some conclusions are mentioned at the end.


## 2. The Generalized Differential Transform Method

For reader’s facility, this part covers a review of the generalized differential transform [28] as well as some fundamental fractional calculus ideas and terminology.The generalized differential transform of thekth derivative of the analytic function fx is given by(4)Fαk=1Γαk+1Dx0αkfxx=x0,where 0<α≤1,Dx0αk=Dx0α.Dx0α…Dx0α,k-times, and Dx0α denotes the Caputo fractional differential operator of order α given by(5)Daαfx=Jam−αDmfx.HereDm is the integer-order differential operator of order m and Jm is the Riemann–Liouville integral operator of order μ with μ>0, which is given by(6)Jaμfx=1Γμ∫axx−tμ−1ftdt,x>0.The generalized differential inverse transform ofFαk is defined as(7)fx=∑k=0∞Fαkx−x0αk,which practically can be approximated by the following finite series:(8)fx≅∑k=0MFαkx−x0αk.Because the initial conditions are represented as integer-order derivatives, the GDTM defines the transformation of the initial conditions as follows:(9)Fαk=1αk!dαkfxdxαk|x=x0,ifαk∈Z+,0,ifαk∉Z,fork=0,1,…,qα−1,where q is the order of considered FDEs.Putting (4) in (7) yields(10)fx=∑k=0∞Fαkx−x0αk.The idea of a generalized differential transform is obtained from generalized Taylor’s formula [28]. It is worth noting that the extended differential transform technique simplifies to the conventional differential transform method when α=1 [32]. Table 1 lists some of the essential features of GDTM derived from equations (4) and (5).Table 1
Basic properties of GDTM [28].
Original functionTransformed functionfx=gx±hxFαk=Gαk±Hαkfx=agxFαk=aGαkfx=gxhxFαk=∑l=0kGαlHαk−lfx=Dx0αgx, 0<α≤1Fαk=Γαk+1+1/Γαk+1Gαk+1fx=x−x0γFαk=δk−γ/α, δk=1,ifk=00,ifk≠0fx=Dx0βgx,m−1<β≤m,m∈Z+Fαk=Γαk+β+1/Γαk+1Gαk+β/αfx=expgxFαk=∑i=0k−1i+1/kGαi+1Fαk−i−1, where Fα0=expGα0


## 3. Solution Procedure

For solving (2) and (3), we proceed with the following algorithm steps.(1)
Choose proper value ofα that satisfies q/α,α/α∈Z+.(2)
Using the generalized differential transform on both sides of equation (3) and the characteristics mentioned in Table 1, we obtain the following recurrence relation:(11)Uαk+qα=βλΓαk+γ+1Γαk+1Uαk+γα−λμ∑m=0kWαmβδk−m−Uαk−m,wherek=0,1,2,…,β=Γαk+1/Γαk+q+1, and Wαk is the generalized differential transform of eux which is given by(12)Wαj=1j∑m=0j−1m+1Uαm+1Wαj−m−1,j=qα,qα+1,…(3)
The boundary conditions given in (2) for x=0 are transformed by employing (9) as follows:(13)Uα0=A,Uα1α=λA,Uαi=0,foriα∉ℤ+,i=1,…,qα−1,whereA=u0 is the initial condition. We can define(14)Wα0=eA,Wα1α=λAeA,Wαi=0,i=1,2…,1α−1,1α−1,…,qα−1.The conditions in (2) for x=1 are transformed by employing (7) as follows:(15)∑k=0NαkUαk=0.(4)
Equations (11) and (13) are utilized to find Uαk up to any N-terms. Then, by using (15), the value of A is evaluated. Also, the Nth order of approximation is(16)uNx=∑k=0NUαkxαk.For simplicity, we can generate the solution forq=1.9,γ=0.9 by assuming α=0.1. Then, applying (11) for k=0,1,2,…,15, we have(17)u34=∑i=034U0.1ixi/10=A+16Axλ6+xλ3+xλ+eAx19/10λ29A−β+10xA2+A−1+Aβλμ29Γ29/10.The condition in (15) gives(18)Aλ+Aλ2+Aλ32+19AeAλμ−eAβλμ10Γ29/10+292AeAλ2μ+A2eAλ2μ−eAβλ2μ−AeAβλ2μ10Γ39/10=0.By fixing the values ofλ,μ, and β, it is easy to solve the equation via the Newton–Raphson method.


## 4. Convergence Analysis

Lemma 1.
(see [33]). The standard power series ∑i=0∞Uiri, r∈ℝ, has a radius of convergence R if and only if the fractional one ∑i=0∞Uαirαi,r≥0, has a radius of R1/μ.Theorem 1.
The fractional power series:(19)∑i=0∞Uαirαi,where the coefficients are mentioned in (11), has a positive radius of convergence.Proof.
From equation (8), we have(20)Uαk≤θ1λμ∑m=0n−q/αWmUαn−m−qα−βδn−m−qα+θ2λΓnα+γ+1Γnα+1Uαγα+n−qα,where(21)θ1=Γn−q/αα+1Γq+n−q/αα+1,and(22)θ2=Γn−q/αα+γ+1Γq+n−q/αα+1Γnα+1Γnα+γ+1.
Now, consider the series(23)zr=∑n=0∞anrn,for which a0=Uα0, a1=Uα1/α, and(24)an=θ1λμ∑m=0n−q/αWmUαn−m−qα−βδn−m−qα+θ2λΓnα+γ+1Γnα+1Uαγα+n−qα,for n=q/α,q/α+1,… and an=0 for 1<n<q/α.
Then, define(25)Ψ=zr=a0+a1r+∑n=0∞an+2rnr2=a0+a1r+r2∑n=0∞θ1λμ∑m=0n−q/αWmYn−m−qα−βδn−m−qα+θ2λΓnα+γ+1Γnα+1Yγα+n−qαrn.
Now, we have the function of two variables:(26)Θr,Ψ=Ψ−a0−a1r−r2θ1λμeΨΨ−βrq+λθ2DγΨrq,which is analytic in the plane r,Ψ with the characteristics Θ0,a0=0 and ΘΨ0,a0=1≠0. Since zr is an analytic function in a neighborhood of the point 0,a0 of the r,Ψ-plane with a positive radius of convergence, then by implicit function theorem, the series in (10) is convergent by Lemma 1.


## 5. Numerical Experiments

This section derives the numerical experiments of the given procedure of Section3 for solving equations (2) and (3). Because the exact solution to the given problem is not known, we instead find the absolute residual error function, which justifies how accurately the numerical solution agrees to the solution of main problems (2) and (3). So, the absolute residual error function is(27)ERNx=dquNxdxq−λdγuNxdxγ+λμβ−uNxexpuNx,0≤x≤1.Firstly, we start to generate the results for standard fractional derivativeq=2,γ=1. By fixing α=1,N=21, Table 2 displays our numerical results and matches them with those of the mentioned outputs in [4] for the same iteration number. Numerical solution using default Mathematica package and the GDTM solution and its absolute error and residual error are given in Table 3. Moreover, Figure 1 shows the approximate solutions for μ=0.7,λ=5, and β=0.8. Figure 1 is in good agreement with Figure 2 given in [5].Table 2
Numerical solutions forμ=0.7,λ=5,β=0.8.
xGDTMMethod in [4]0.00.101646231514738255570.101646231064122758510.10.151608003475384056640.151608002444185139700.20.199690826938690141330.199690825987745986300.30.245647640713457210420.245647639198821209860.40.289195528463070561380.289195526526831751590.50.329973876788146658350.329973875007283831840.60.367466685033411122330.367466683425880139690.70.400861474839358078450.400861474627291976290.80.428795636271987135050.428795638521272575240.90.448902769569831106010.448902770897070312011.00.457005437961182814370.45700543763742257810Table 3
GDTM solution and numerical solution with its absolute and residual errors forλ=5,μ=0.7,β=0.8.
xGDTMNumerical solutionAbsolute errorResidual error0.00.10164623070.10164623158.085193209×10−101.665334537×10−160.10.15160800190.15160800351.624451518×10−96.106226635×10−160.20.1996908250.19969082691.943355393×10−94.996003611×10−160.30.24564763750.24564764073.187346082×10−95.828670879×10−160.40.28919552370.28919552854.774220019×10−91.276756478×10−150.50.32997387020.32997387686.62490951×10−98.398837181×10−140.60.36746667510.3674666859.92236443×10−91.042421705×10−110.70.40086146030.40086147481.456031457×10−86.394951235×10−100.80.42879561370.42879563632.26142593×10−82.315349923×10−80.90.44890272820.44890276964.136782861×10−85.61353519×10−71.00.45700537630.4570054386.16973152×10−89.929563344×10−6Figure 1
Graph of numerical outputs forλ=5,μ=0.7, and β=0.8.Figure 2
GDTM solution:q=1.9 and γ=1 (line); q=1.9 and γ=0.9 (dashed).Now, we explore the impact of the fractional derivative on the solution of the model. In Figure3, we show the convergence of obtained missing condition A with increasing N when q=1.9 and γ=0.9. It is clear that the value of A starts to be fixed when N>80.Figure 3
Variation ofA with N for the case λ=5,μ=0.05,β=0.53.Tables4 and 5 introduce the approximate solutions and the residual errors of problems (2) and (3) for different values of q,γ,λ,μ, and β. The residual error indicator demonstrates that the results are accurate for at least 6×10−5.Table 4
Numerical outputs forλ=5,μ=0.05,β=0.53.
xq=1.9,γ=1.0q=1.9,γ=0.9GDTMResidual errorGDTMResidual error0.00.005242585.55112×10−170.006128911.69136×10−170.10.007857492.77556×10−170.009051659.54098×10−180.20.01045811.5786×10−160.01193414.55365×10−170.30.0130392.43347×10−140.01489691.82536×10−150.40.01558944.46358×10−120.01799962.99132×10−130.50.01809012.4929×10−100.02128231.55918×10−110.60.02050546.57931×10−90.02477513.86786×10−100.70.02277071.03475×10−70.02850215.73491×10−90.80.0247681.11374×10−60.03248325.82065×10−80.90.02628228.96982×10−60.03673544.41167×10−71.00.02692095.7448×10−50.04127362.64861×10−6Table 5
Numerical solutions forλ=0.05,μ=0.5,β=0.6.
xq=1.9,γ=0.9q=1.9,γ=1.0GDTMResidual errorGDTMResidual error0.00.2334823.46945×10−180.2335280.0.10.2345723.46945×10−180.234623.46945×10−180.20.235533.46945×10−180.235580.0.30.2363655.20417×10−180.2364181.73472×10−180.40.2370835.20417×10−180.2371381.73472×10−180.50.2376861.73472×10−180.2377440.0.60.2381760.0.2382360.0.70.2385565.20417×10−180.2386180.0.80.2388263.46945×10−180.2388891.73472×10−180.90.2389871.73472×10−180.2390511.73472×10−181.00.2390415.20417×10−180.2391050.To study the solution behavior, we plot the present solution for the case ofλ=5,μ=0.05,β=0.53 when q=1.9, γ=0.9 and q=1.9, γ=1 in Figure 2. The solution for the case of λ=0.05,μ=0.5,β=0.6 when q=1.9, γ=0.9 and q=1.9, γ=1 is presented in Figure 4. It is clear that the fractional derivatives can change the solution behavior.Figure 4
GDTM solution:q=1.9 and γ=1 (line); q=1.9 and γ=0.9 (dashed).


## 6. Conclusions

Strongly nonlinear boundary value problem with multi-order fractional derivative that occurred in the chemical reaction theory has been successfully solved via the new algorithm based on GDTM. The constructed solution has been given in terms of convergent infinite series as seen in the provided theorem. The method was easy to apply, and the results have enough good accuracy as shown in the experimental results. The obtained solution was directly generated without any linearization or discretization of the domain. The given method is very powerful and can be easily applied for several kinds of fractional nonlinear boundary value problems in future.



---

*Source: 8655340-2022-06-18.xml*