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.

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]:where , and are the Péclet number, Damköhler number, and adiabatic temperature rise, respectively. The relative boundary conditions are given by

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:where and are the Caputo derivative operators along with fractional orders , subject to the boundary conditions. (1) can be considered as a particular form of (3) by fixing the orders and . 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 [2127] 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) [2831].

This paper is organized as follows. In Section 2, 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 the th derivative of the analytic function is given bywhere -times, and denotes the Caputo fractional differential operator of order given by

Here is the integer-order differential operator of order and is the Riemann–Liouville integral operator of order with , which is given by

The generalized differential inverse transform of is defined aswhich practically can be approximated by the following finite series:

Because the initial conditions are represented as integer-order derivatives, the GDTM defines the transformation of the initial conditions as follows:where is the order of considered FDEs.

Putting (4) in (7) yields

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 [32]. Table 1 lists some of the essential features of GDTM derived from equations (4) and (5).

3. Solution Procedure

For solving (2) and (3), we proceed with the following algorithm steps.(1)Choose proper value of that satisfies .(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:where , and is the generalized differential transform of which is given by(3)The boundary conditions given in (2) for are transformed by employing (9) as follows:where is the initial condition. We can defineThe conditions in (2) for are transformed by employing (7) as follows:(4)Equations (11) and (13) are utilized to find up to any -terms. Then, by using (15), the value of is evaluated. Also, the th order of approximation is

For simplicity, we can generate the solution for by assuming . Then, applying (11) for , we have

The condition in (15) gives

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 , , has a radius of convergence if and only if the fractional one has a radius of .

Theorem 1. The fractional power series:where the coefficients are mentioned in (11), has a positive radius of convergence.

Proof. From equation (8), we havewhereandNow, consider the seriesfor which , , andfor and for .
Then, defineNow, we have the function of two variables:which is analytic in the plane with the characteristics and . Since is an analytic function in a neighborhood of the point of the -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 Section 3 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

Firstly, we start to generate the results for standard fractional derivative . By fixing , 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 , and . Figure 1 is in good agreement with Figure 2 given in [5].

Now, we explore the impact of the fractional derivative on the solution of the model. In Figure 3, we show the convergence of obtained missing condition with increasing when and It is clear that the value of starts to be fixed when .

Tables 4 and 5 introduce the approximate solutions and the residual errors of problems (2) and (3) for different values of , and . The residual error indicator demonstrates that the results are accurate for at least .

To study the solution behavior, we plot the present solution for the case of when , and , in Figure 2. The solution for the case of when , and , is presented in Figure 4. It is clear that the fractional derivatives can change the solution behavior.

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.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


Marina Murillo-Arcila was supported by MCIN/AEI/10.13039/501100011033, Project no. PID2019-105011GBI00, and by Generalitat Valenciana, Project no. PROMETEU/2021/070.