Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method

: In this study, we present the generalized form of the higher-order nonlinear fractional Bratu-type equation. In this generalization, we deal with a generalized fractional derivative, which is quite useful from an application point of view. Furthermore, some special cases of the generalized fractional Bratu equation are recognized and examined. To solve these nonlinear differential equations of fractional order, we employ the homotopy perturbation transform method. This work presents a useful computational method for solving these equations and advances our understanding of them. We also plot some numerical outcomes to show the efficiency of the obtained results.


Introduction
Numerous areas of pure and applied mathematics focus largely on fractional calculus, which works with derivatives and integrals of arbitrary order [1][2][3][4][5].Fractional differential equations have been used more frequently in applied mathematics and physics, as well as in the modeling and interpretation of various real-world issues.Biological science, physical science, chemical engineering, and other fields of research indicate nonlinear perspectives [6][7][8][9].
In mathematical modeling, the nonlinear fractional differential equations are quite helpful.However, it becomes more challenging to solve nonlinear systems of equations.It has been seen that there are many differential equations without an exact solution.Various innovative numerical and analytical approaches have been established to solve fractional differential equations, leading to numerical and analytical solutions to these types of problems.These techniques consist of Adomian decomposition method [10][11][12], homotopy analysis method [13][14][15]), generalized differential transform approach [16], homotopy perturbation method [17,18], modified Laplace decomposition method [19], homotopy perturbation transform method [20,21], and q-homotopy analysis transform method [22,23], etc.The homotopy perturbation method is the most effective for solving the differential equations of non-linear type because of its straightforward approach and fastest convergence rate.This technique was developed by the famous mathematician Mr. He [24][25][26].Its main benefit is that the solution determined in series form approaches its exact solution very quickly.
Many mathematicians have developed an appropriate definition and formula for fractional calculus.Major contributions to the theory of fractional calculus with singular kernels have been made by Riemann, Caputo, Liouville, Weyl, Kilbas, and others.Gorenglo, Miller-Ross, Yang-Abdel-Cattani, Atangana-Baleanu, Wiman, and many others encouraged the study of fractional integrals and derivatives having nonsingular kernels.Miller-Ross, Yang-Abdel-Cattani, Wiman, Atangana-Baleanu, and many others promoted the study of fractional integrals and derivatives with non-singular kernels.Riemann-Liouville and Caputo are the two fractional operators that are used most frequently.However, there are several problems that are not properly described by these operators.For fractional derivatives, Caputo-Fabrizio developed a new operator in 2015 with a nonsingular kernel [27].This nonsingular kernel makes the Caputo-Fabrizio fractional derivative stronger than the Caputo derivative.
The calculus operators of fractional type have been divided into three groups by Gomez and Atangana: weak, medium, and strong [28].They suggest that the Riemann-Liouville fractional derivative is significantly more useful in explaining physical difficulties than the Caputo version.However, the Mittag-Leffler memory pattern of the Atangana-Baleanu fractional operator provides a good clarification.Some authors have used fractional derivatives of Atangana-Baleanu-Caputo to characterize fractional differential equations [29].
A numerical slab model of a combustion problem is examined using the Bratu-type equation [30].Several additional issues have been solved using Bratu's problem, such as the Chandrasekhar model, the fuel ignition model, the thermal response model, and the framework for the electrospun nano-fiber creation method [31,32].To solve the Bratu type problem, Ghazanfari and Sephvandzadeh used the homotopy perturbation method [33], the modified variational iteration method [34], and the Adomian decomposition method [35].Recent studies have explored numerical approaches to addressing the fractional Bratu initial value problem.Investigations have been conducted using the Bézier curve [36], the fractional differential transform method [38], and a technique grounded in the CAS wavelet scheme [37] for solving fractional Bratu-type differential equations.
Applications for nonlinear differential equations, like the Bratu-type equation, can be found in many branches of science and engineering.It is used in heat transfer to investigate temperature distributions in materials, in combustion theory to understand flame propagation, and in chemical kinetics to replicate complex reaction kinetics.The modeling of complex phenomena like population dynamics, material diffusion, reactor neutron transport, electrochemical reactions, and pollutant spread is made possible by its involvement in mathematical biology, material science, nuclear physics, electrochemistry, and environmental science.In addition to having particular uses, the Bratu-type equation is a basic mathematical model that provides an understanding of the characteristics of nonlinear differential equations.Because of this, it is a flexible and essential instrument for scientific analysis and study.The "Nonlinear Local Fractional Bratu-Type Equation" is a delightful and theoretically complex differential equation that belongs to the field of fractional calculus.
As mathematicians, physicists, and engineers grapple with issues related to nonstandard calculus, nonlinear dynamics, and complex systems, they are entranced by ongoing research projects centered around the "Non-linear Local Fractional Bratu-type Equation".Having a thorough grasp of the characteristics and solutions of this equation could improve our understanding of a variety of real-world processes and possibly lead to useful applications in fields including chemical kinetics, biological modeling, and combustion theory.In the fields of science and engineering, the Bratu-type equation is quite useful.In all domains where standard differentiable calculus is relevant, this article presents some new insights into the concept of generalization.Other expansions of the differentiability concept can be found in [39][40][41].Very few researchers have studied the generalized form of a fractional Bratu-type equation, so finding the solution of a generalized fractional Bratu-type differential equation is a new thing in itself.Realizing how important it is, we turn our attention to examining a generalized fractional Bratu-type equation.
This study aims to define the generalized Bratu-type differential equation and expand the application of the homotopy perturbation transform method for identifying specific features of the generalized Bratu-type differential equation where δ is constant and E β (u) is the Mittag-Leffler function defined as , Γ indicates the gamma function and The structure of this paper is divided into seven sections: We provide a brief introduction and some fundamental terminology in Section 1.The basic homotopy perturbation transform approach is defined in Section 2. In Section 3, we discuss the existence and uniqueness of the generalized fractional Bratu-type differential equation.In Section 4, we describe the main findings of this work.In Section 5, we discuss the convergence analysis.In Section 6, we describe some special cases of generalized Bratu-type equations and provide a graphic representation of the solution.This paper is concluded in the last section.

Basic Definitions
This section provides the definitions and foundational results essential for the content presented in this paper.
A theory that deals with generalized fractional-order integration and differentiation operators is called generalized fractional calculus.Using an analogy with the concept of the Caputo and R-L fractional derivative, the analogous generalized fractional derivatives are defined as [42] C where 0 < κ < 1 is the order of derivative, f : [0, +∞) → R is the absolute continuous function with ḟ ∈ L 1 loc (0, +∞), ∆ l is known as the general kernel function, also for f ∈ L 1 loc (0, +∞), where 0 D −κ t represent the general form of the Riemann-Liouville fractional integral defined as Definition 1. Caputo derivative of fractional order κ is denoted by C a D κ t and defined by [43] C Definition 2. Caputo-Fabrizio derivative is denoted by CF a D κ t and defined by [44] CF where Definition 3. We have the Atangana-Baleanu derivative which is denoted by ABC a D κ t and defined by [45] ABC a where Definition 5.In accordance with the reference [47], we define the Sumudu transform for the Caputo derivative of order κ > 0 Definition 6.In accordance with the reference [47], we define the Sumudu transform for the Caputo-Fabrizio derivative of order κ > 0 where ∆(κ) = Normalization function with ∆(0) = 1, ∆(1) = 1.

Definition 7.
In accordance with the reference [47], we define the Sumudu transform for the Atangana-Baleanu fractional derivative of order κ > 0 where

Basic Homotopy Perturbation Transform Approach
This section covers the Homotopy Perturbation Transform Method (HPTM) procedure for solving fractional differential equations using generalized fractional operators.
The Homotopy Perturbation Technique (HPT), incorporating He's polynomials and a transform algorithm, is formulated as the homotopy perturbation transform method.To illustrate the fundamental workings of the homotopy perturbation transform method, we apply it to a nonlinear fractional-order differential equation of the form Here, u(t) represents the probability density function, D κ t denotes the generalized fractional derivative with order κ, L stands for the linear differential operator with an order less than that of D, and N represents the general nonlinear differential operator.
Apply the Sumudu transform on both sides of (9), as we know, Sumudu transform for generalized fractional derivatives [48], here, and ∆(κ) is the normalization function. Hence Now, applying the inverse Sumudu transform, we obtain where p n H n (u), here H n (u) is He's polynomial defined as , n = 0, 1, 2, 3, ... , and p ∈ (0, 1) is an embedding parameter. Then This equation represents the combination of the homotopy perturbation method, He's polynomial, and Sumudu transform.Now, by equating the coefficients corresponding to different powers of p, we obtain The remaining elements of u n (t), n ≥ 0 can be evaluated by applying a similar approach.As a result, the homotopy perturbation transforms method solutions are completely computed.
Thus, the approximate analytical solution of ( 9) is u n (t).

Existence and Uniqueness of the Solution
Theorem 1.The solution of the generalized Bratu-type fractional differential equation defined by (1) satisfies the Lipschitz condition and has the unique solution for 0 < b < 1, where b = ν|δ| Γ(1+κ) t κ is known as the Lipschitz constant.
Proof.Firstly, we write Equation (1) in the following form where The non-linear part G(t, u(t)) is Lipschitz continuous with the condition Lipschitz constant can be evaluated by making use of maximum norm ∥G∥ = max 0≤t≤1 we have Since the series is convergent where M is some constant. Thus, Now from (12), we obtain Thus, we can select the Lipschitz constant as ν = . Hence, the solution of Equation ( 1) exists.Now, suppose that Equation (1) has two different solutions u 1 and u 2 , so then This produces the following inequality, , where b = ν|δ| Γ(1+κ) t κ .Thus, ∥u 1 − u 2 ∥ ≤ 0 for 0 < b < 1 which shows that u 1 = u 2 .This proves that the solution of (1) is unique.

Main Result
The solution to the generalized Bratu-type equation using the homotopy perturbation transform method is given in this section.
Theorem 2. The solution of the generalized Bratu-type fractional differential equation is given by the series u = ∞ ∑ n=0 u n (t).
where δ is constant, E β (u) is the Mittag-Leffler function, and D 2nκ is a generalized differential operator with order 2n.
Proof.We have the generalized Bratu-type fractional differential equation Applying Sumudu transform, Here, where ∆(κ) is the normalization function.
Applying the inverse Sumudu transform Now, using the homotopy perturbation transform method, Equation ( 15) becomes where We have the Taylor's series expansion of E β (u), .
Also, we have, Now comparing the coefficient of the power of p in ( 16), we obtain The remaining elements of u n (t), n ≥ 0 can be evaluated by applying a similar approach through using ( 17) wherein the series solutions are completely determined.
We have The generalized differential operator reduces into three well-known differential operators, namely, the Atangana-Baleanu-Caputo, Caputo, and Caputo-Fabrizio fractional operators.Here, we discuss all three cases of a generalized Bratu-type fractional differential equation.
Proof.Let us consider how the sequence of partial sums of series ( 18) is {S n (t)} n≥0 defined by Now, first we have to show that {S n (t)} be a Cauchy sequence in Banach space Thus, S q (t) − S r (t) q,r→∞ → 0 as u 0 is bounded.This {S n (t)} is a Cauchy sequence in Banach space and hence convergent.Therefore, there exists u ∈ B such that We achieve the desired result.

Special Cases
We have the normal generalized fractional operator.However, it may be reduced into three well-known fractional operators: Atangana-Baleanu-Caputo, Caputo, and Caputo-Fabrizio.To derive the results, we separately discuss the cases for generalized operators.
If we take n = 1, then (1) becomes Moreover, consider the initial point u(0) = 0.This equation is known as a particular form of the generalized fractional Bratu equation.
We offer the solution to the generalized fractional Bratu equation in this section by using the process defined in Section 4.
Case 1.Consider generalized Bratu's initial value problem for the Caputo fractional derivative, with u(0) = 0. Applying Sumudu transform on (23), Applying the inverse Sumudu transform, Now, using the homotopy perturbation transform method and expanding u(t) as u Now, we have the Taylor's series expansion of E β (u), .

Conclusions
In this investigation, we have introduced and addressed the generalized nonlinear Bratu-type equation incorporating the generalized fractional derivatives operator.The solution is derived by applying an innovative analytical technique known as the homotopy perturbation transformation method.We have demonstrated the effectiveness of the generalized fractional derivatives operator and presented reduced results expressed in wellknown fractional operators, including Caputo, Caputo-Fabrizio, and Atangana-Baleanu.Furthermore, various special cases of the generalized fractional Bratu's equation have been identified and thoroughly examined.We also plot some numerical results to show the behavior of the obtained results.
Now, expanding u(t) into the series of power p as u(t) = ∞ ∑ n=0 p n u n (t) and considerNu(t) = ∞ ∑ n=0

Figure 1 .
Figure 1.Change in u for Case 1 with respect to time, for different values of κ.

Figure 2 .
Figure 2. Change in u for remark of Case 1 with respect to time, for different values of κ.
Now to plot the results by considering δ = 1.0, β = 0.7 for different values of, κ as illustrated in Figure3.

Figure 3 .
Figure 3. Change in u for Case 2 with respect to time, for different values of κ.

Figure 4 .
Figure 4. Change in u for remark of Case 2 with respect to time, for different values of κ.