Prospective Analysis of Time-Fractional Emden–Fowler Model Using Elzaki Transform Homotopy Perturbation Method

: The present study presents a combination of two famous analytical techniques for the analytical solutions of linear and nonlinear time-fractional Emden–Fowler models. We combine the Elzaki transform (ET) and the homotopy perturbation method (HPM) for the development of the Elzaki transform homotopy perturbation method (ET-HPM). In this paper, we demonstrate that the Elzaki transform (ET) simplifies fractional differential problems by transforming them into algebraic formulas within the transform space. On the other hand, the HPM has the ability to discretize the nonlinear terms in fractional problems. The fractional orders are considered in the Caputo sense. The main purpose of this strategy is to use an alternative approach that has never been employed in the time-fractional Emden–Fowler model. This strategy does not require any variable or hypothesis constraints that ruin the physical nature of the actual problem. The derived series yields a convergent series using the Taylor series formula. The analytical data and visual illustrations for several kinds of fractional orders validate the effectiveness of the suggested scheme. The significant results demonstrate that our recommended strategy is quick and simple to use on fractional problems.


Introduction
In the twenty-first century, there has been a significant amount of interest in fractional calculus (FC) and its diverse applications in mathematical science, astronomy, and biological sciences.It is widely used in different applications of sciences and technology, including device control concepts, computer networks, statistical computation, optical science, electrical chemistry, signal analysis, and chemical substances.The fractional derivative has a global scope instead of being limited to a local area; therefore, it is extremely beneficial, and unlimited variations are anticipated [1,2].It is capable of taking impacts into account more precisely.This nonlocality is extremely useful for analyzing physical reactions involving memory effects.These systems are challenging to investigate using traditional calculus.Consequently, this topic has been focused on in various fascinating studies during the last few decades.FC is a widely recognized mechanism in a variety of scientific and technological domains, and it has been explained by fractional differential equations [3,4].Several other phenomena, such as electromagnetic, hydrodynamic, thermal, acoustic, and electrodynamic phenomena, can also be effectively simulated using these equations.To date, numerous researchers have shown that it can handle multiple challenges, particularly in mathematics and physical science.Indeed, the reliability of fractional operators has demonstrated their suitability for simulating the local instability of components in time or space, which has a critical function in many circumstances; however, it cannot be understood using normal mathematical approaches [5,6].Many scientists have studied the approximate solutions of these fractional differential problems and showed that Fractal Fract.2024, 8, 363 2 of 15 these problems appear substantially more complex than integer-order variants for dealing with accurate results [7][8][9][10].
The time-fractional Emden-Fowler equation is employed in modeling complex systems across various scientific and engineering disciplines.It simulates the propagation of waves in media with memory, where waves exhibit non-standard diffusion characteristics as a result of the fractional time derivative.In addition, this fractional model is being investigated in the context of nonlinear dynamics and pattern formation in a variety of physical systems, including those exhibiting fractal and chaotic dynamics.The astrophysical researchers Homer Lane and Robert Emden first examined the Lane-Emden model, in which they explored the heat distribution of a sphere-shaped ball of oxygen reacting through molecular interactions with other molecules under the traditional rules of thermodynamics [11].Most of the singular problems described by Lane-Emden models have been applied in various applications of applied sciences [12].The Emden-Fowler model is a form of differential equation that is used in astrophysics and computational astrophysics.Due to singularity actions at the point (y = 0), the solution of the Emden-Fowler problem and its initial value problems at the singularity are computationally difficult.Chowdhury and Hashim [13] investigated the use of the HPM to find approximate solutions for second-order ordinary differential problems of the generalized Emden-Fowler form.Mall and Chakraverty [14] introduced a neural network scheme to solve singular initial value problems of Emden-Fowler equations.Syam [15] proposed a scheme based on series results to derive the analytical outcomes of higher-order Emden-Fowler equations [16].
The homotopy perturbation method (HPM) is one of the particularly prominent semianalytical approaches that combine homotopy and perturbation strategies because there are no strict constraints on the selection of its linear operator, and its series solution probably frequently fails to converge [17].He and El-Dib [18] argued that their proposed approach depends on the traditional Taylor series, which does not need a higher-level approximation.The strategy of the HPM has been applied to various problems, such as algebraic problems [19], nonlinear diffusion equations [20], initial and boundary value problems [21], nonlinear wave problems [22], harmonic oscillation [23], and other challenges [24][25][26].Tarig Elzaki [27] developed the scheme of Elzaki transformation to obtain the results of some ordinary and partial differential equations in a time analysis more easily.Later, Aggarwal et al. [28] established the definitions of duality relations for numerous advanced integral transformations.They enable the analysis of a rapidly converging series that leads to precise solutions.This approach is regarded as highly effective in resolving differential problems including both integer and noninteger orders, which may be linear or nonlinear and homogeneous or nonhomogeneous.
In the current study, we develop a strategy called the Elzaki transform homotopy perturbation method (ET-HPM) for an analytical view of the time-fractional Emden-Fowler model.The Elzaki transform simplifies the handling of fractional operators, enhances convergence properties, provides straightforward operational rules, and improves the treatment of initial conditions.It offers the direct treatment of fractional differential equations despite the further transformation required.This can substantially speed up the solution process and decrease potential errors caused by several conversions.These features make the Elzaki transform an effective tool for solving fractional differential equations and provide an alternative approach to conventional methods.By utilizing the HPM, we derive He's polynomials with an easy computational process.Thus, outcomes like iterative series that lead toward the precise results of fractional problems are achieved.Our proposed scheme is independent of any hypotheses, constraints, or restrictions on variables that may ruin the actual problem.This scheme has more direct applications to the fractional problem than other schemes discussed in the literature.The remaining work is designed as follows: We present the highlights of FC and features of the ET in Section 2. Section 3 provides an overview of the proposed strategy, along with certain convergence theorems.In Section 6, we implement the ET-HPM formulation to derive the outcomes of the time-fractional Emden-Fowler model.Section 7 displays the graphical structures of the obtained results with respect to various fractional orders.At last, we provide our conclusion in Section 8.

Definitions and Fundamental Principles
This section presents basic definitions concerning FC and the ET with their properties.Definition 1.The expression for the fractional-order Riemann-Liouville integral is [29] J α ϑ(φ) = 1 The fractional-order Caputo's derivative is expressed as [29] Let A be a set of functions expressed as [27,30] A = ϑ(t) : ∃M, where M represents the constant of a finite term.On the other hand, q 1 and q 2 can be either finite or infinite.Thus, the integral formulation of the Elzaki transform is expressed as where σ is the transform parameter of is the inverse of the ET.
Definition 4. The ET of a fractional operator in Caputo form is described as [31]

Statements:
The ET has the following features in differential cases [32]: The ET for the nth derivative is defined as

Formulation of ET-HPM
This section presents the development of the ET-HPM strategy for the analytical treatment of the time-fractional Emden-Fowler model.The Elzaki transform is useful for handling initial constraints and transforming differential equations into algebraic form.This scheme expands the solution into a power series of the embedding parameter p and solves the transformed equations iteratively.The greatest advantage of this scheme is that it is useful for solving fractional differential equations, which are often challenging for traditional methods in comparison to the traditional HPM.The process of this technique begins with the consideration of a nonlinear fractional differential problem as follows: subject to the condition where ϑ represents the function of time φ, L and N show linear and nonlinear operators, and g(ψ, φ) is a known term.Applying the ET to Equation (1), we obtain Using the propositions of the ET, we obtain Applying the inverse ET, we obtain ( Equation ( 3) is called the formulation of the ET-HPM of Equation ( 1), where Consider the following solution of Equation (1): where p ∈ [0, 1] is a small homotopy parameter, and ϑ 0 (ϑ, φ) is the initial condition.The nonlinear component of the homotopy polynomial is considered as follows: which can be calculated as By using Equations ( 4) and ( 5), Equation (3) becomes With the comparison of p on both sides, we obtain In the end, our analytical solution acts in series, as follows: Now, we present a theorem that explains the conditions of the convergence for the series solutions (6).

Convergence Analysis
This section presents a study of the convergence theorem and states its proof for the analytical results of the considered fractional model under the suggested scheme.We show that the obtained results in the form of a series converge to the precise results of the problem very rapidly.
Theorem 1 (The Banach fixed-point theorem).Consider X a Banach space to the nonlinear mapping T : X → X , and also suppose that It is stated that T represents a singular fixed point, and the sequence produced by the ET-HPM is considered to converge to ϑ n+1 = Tϑ n with an arbitrary ϑ 0 ∈ X, and therefore, The aforementioned theorem serves as a prerequisite for the subsequent analysis, which is explicable through the utilization of the Banach fixed-point theorem.

The Uniqueness Theorem
Let the analytical results of a fractional model obtained by the ET-HPM be unique whenever 0 < γ < 1, i.e., Proof.Since Equation ( 14) is the derived solution of the ET-HPM, we have where L and N agree with the Lipschitz constraints.Applying the ET, we obtain Using the properties of the ET, we have Now, applying the inverse ET to Equation ( 15), we obtain Consider that there are two unique results, ϑ(ψ, φ) and ω(ψ, φ); we have .
From the triangle inequality, we can deduce that By applying the convolution theorem, Given that L and N agree with the Lipschitz constraints, we can conclude that L is a bounded operator with the property that |L(ϑ) − L(ω)|≤ µ|ϑ − ω|, where µ is a constant.Additionally, N is defined as |N(ϑ) − N(ω)|≤ ε|ϑ − ω| for any positive value of ε.Thus, Equation ( 16) becomes By applying the mean value theorem for integrals to Equation (17), we can establish that M is the maximum value of M = max(φ − τ) n for φ in the interval [0, φ].Consequently, Equation ( 17) can be expressed as Fractal Fract.2024, 8, 363 This means that ϑ = ω whenever γ < 1 and γ ∈ (0, 1).Hence, the solution is unique.□

Numerical Applications
Here, we implement the proposed scheme to derive the approximate results of the time-fractional Emden-Fowler model.We consider two examples: linear and nonlinear models.We note that the resultant series converges to the precise responses just after a couple of cycles, which shows the accuracy and authenticity of the proposed scheme.We adopt Mathematica 11 in our computational processes and graphical illustrations.

Problem 2
Next, assume the following homogeneous nonlinear time-fractional Emden-Fowler model: subject to the condition Employing the ET with Equation ( 25), we obtain 1 Taking the inverse ET, we obtain By adopting the HPM in Equation ( 27), we can obtain After analyzing identical parts of p in Equation ( 28), we obtain The following possible outcomes are obtained: . . .Hence, the analytical solution of Equation ( 25) yields In the case of α = 1, the precise result of Equation ( 25) yields

Numerical Findings and Analysis
This section reveals the physical interpretation of linear and nonlinear time-fractional Emden-Fowler problems.Our aim is to present the physical behavior of these problems with different fractional orders.In our first problem, we consider a linear time-fractional Emden-Fowler problem with a sufficient condition.
Figure 1a displays the physical behavior of the ET-HPM solution for a fractional order α = 1, whereas Figure 1b shows the physical behavior of the precise results including 0 ≤ ψ ≤ 1 and 0 ≤ φ ≤ 2 for Problem 1.The 3-D graphical visuals show that both graphs are close and reveal the significance of our proposed scheme.Figure 1c demonstrates the contour plot, which examines the correlation across three variables on a single graph and the expression of ψ and φ that gives the ideal results of ϑ(ψ, φ).The contour lines relate multiple combinations of the variables that generate similar parameters of ϑ(ψ, φ) at 0 ≤ ψ ≤ 1 and 0 ≤ φ ≤ 1. Figure 1d depicts the error distribution of ϑ(ψ, φ) at 0 ≤ ψ ≤ 1 and φ = 0.1 with α = 0.6 (red), α = 0.8 (blue), and α = 1 (green) and the precise result (the dotted line).Table 1 shows the error distribution between the ET-HPM and precise results.The ET-HPM results are obtained at different fractional orders, such as α = 0.4, 0.6, 0.8, and 1.By increasing the parameter of the fractional order, our computed results become closer to the precise results.
1.By increasing the parameter of the fractional order, our computed results become closer to the precise results.Similarly, we display the physical behavior of the ET-HPM solution for a fractional order α = 1 in Figure 2a, whereas Figure 2b shows precise results including 0 ≤ ψ ≤ 0.05 and 0 ≤ φ ≤ 0.01 for Problem 2. The 3-D graphical visuals show that both graphs are close and reveal the significance of our proposed scheme.Figure 2c demonstrates the contour plot, which examines the correlation across three variables on a single graph and the expression of ψ and φ that gives the ideal results of ϑ(ψ, φ).The contour lines relate multiple combinations of the variables that generate similar parameters of ϑ(ψ, φ) at 0 ≤ ψ ≤ 1 and 0 ≤ φ ≤ 1. Figure 2d depicts the error distribution of ϑ(ψ, φ) at 0 ≤ ψ ≤ 0.1 and φ = 0.01 with α = 0.6 (red), α = 0.8 (blue), and α = 1 (green) and the precise result (the dotted line).Table 2 represents the error distribution between the ET-HPM and precise results.The ET-HPM results are derived at various fractional orders, such as α = 0.4, 0.6, 0.8, and 1.By extending the fractional-order value, we obtain results that are close to the precise results.
Similarly, we display the physical behavior of the ET-HPM solution for a fractional order α = 1 in Figure 2a, whereas Figure 2b shows precise results including 0 ≤ ψ ≤ 0.05 and 0 ≤ φ ≤ 0.01 for Problem 2. The 3-D graphical visuals show that both graphs are close and reveal the significance of our proposed scheme.Figure 2c demonstrates the contour plot, which examines the correlation across three variables on a single graph and the expression of ψ and φ that gives the ideal results of ϑ(ψ, φ).The contour lines relate multiple combinations of the variables that generate similar parameters of ϑ(ψ, φ) at 0 ≤ ψ ≤ 1 and 0 ≤ φ ≤ 1. Figure 2d depicts the error distribution of ϑ(ψ, φ) at 0 ≤ ψ ≤ 0.1 and φ = 0.01 with α = 0.6 (red), α = 0.8 (blue), and α = 1 (green) and the precise result (the dotted line).Table 2 represents the error distribution between the ET-HPM and precise results.The ET-HPM results are derived at various fractional orders, such as α = 0.4, 0.6, 0.8, and 1.By extending the fractional-order value, we obtain results that are close to the precise results.Contour graphics are an effective visualization tool for representing three-dimensional data in two dimensions.The contour graphs demonstrate the result ϑ(ψ, φ) of the timefractional Emden-Fowler equation across a domain in the space ψ and time φ planes.These graphical simulations are helpful in understanding the dynamic behavior of the obtained results.The contour lines depict the constant parameters of the results ϑ(ψ, φ), providing a visual representation of the solution's evolution over space and time.Only three iterations are considered for analytical computation results.It is evident that the accuracy of the results can be significantly enhanced by incorporating additional parameters, and the errors will come closer to zero.By extending the value of α, the nonlinearity impact is influenced, Contour graphics are an effective visualization tool for representing three-dimensional data in two dimensions.The contour graphs demonstrate the result ϑ(ψ, φ) of the timefractional Emden-Fowler equation across a domain in the space ψ and time φ planes.These graphical simulations are helpful in understanding the dynamic behavior of the obtained results.The contour lines depict the constant parameters of the results ϑ(ψ, φ), providing a visual representation of the solution's evolution over space and time.Only three iterations are considered for analytical computation results.It is evident that the accuracy of the results can be significantly enhanced by incorporating additional parameters, and the errors will come closer to zero.By extending the value of α, the nonlinearity impact is influenced, even though the wave amplitude decreases.We have observed that the ET-HPM is entirely designed for handling both linear and nonlinear time-fractional models.

Conclusions
In this paper, we successfully obtained the analytical solutions of linear and nonlinear time-fractional Emden-Fowler models via the ET-HPM strategy.Two numerical problems were adequately examined traditionally using this proposed approach.The fractional orders are introduced in Caputo form.Our proposed scheme shows excellence in obtaining iterative results that are very close to the precise solution.By simplifying the fractional problem in the transform domain, the ET-HPM can provide more valuable convergence aspects of the solution series.The suggested approach generates convergence series systems with simple identifiable variables without the need for linearization, interference, or restricting constraints.The analytical and visualized findings obtained by the recommended approach are considerably more appealing and effective in determining the solutions of time-fractional models.This approach has significant promise for solving fractional differential problems.We will consider expanding this scheme for real-life problems with larger dimensions to other fractal and fractional phenomena encountered in engineering and science trends in future work.

Table 1 .
A comparison of the exact and ET-HPM results of Problem 1 in Section 6.1 at different fractional orders.

Table 1 .
A comparison of the exact and ET-HPM results of Problem 1 in Section 6.1 at different fractional orders.

Table 2 .
A comparison of the exact and ET-HPM results of Problem 2 in Section 6.2 at different fractional orders.