Next Article in Journal
Second-Order Dynamic Equations with Noncanonical Operator: Oscillatory Behavior
Next Article in Special Issue
Order of Convergence, Extensions of Newton–Simpson Method for Solving Nonlinear Equations and Their Dynamics
Previous Article in Journal
Frac-Vector: Better Category Representation
Previous Article in Special Issue
Stability and Bifurcation Analysis of Fifth-Order Nonlinear Fractional Difference Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectral Collocation Approach via Normalized Shifted Jacobi Polynomials for the Nonlinear Lane-Emden Equation with Fractal-Fractional Derivative

by
Youssri Hassan Youssri
1,* and
Ahmed Gamal Atta
2
1
Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt
2
Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Cairo 11341, Egypt
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(2), 133; https://doi.org/10.3390/fractalfract7020133
Submission received: 6 January 2023 / Revised: 24 January 2023 / Accepted: 28 January 2023 / Published: 31 January 2023
(This article belongs to the Special Issue Applications of Iterative Methods in Solving Nonlinear Equations)

Abstract

:
Herein, we adduce, analyze, and come up with spectral collocation procedures to iron out a specific class of nonlinear singular Lane–Emden (LE) equations with generalized Caputo derivatives that appear in the study of astronomical objects. The offered solution is approximated as a truncated series of the normalized shifted Jacobi polynomials under the assumption that the exact solution is an element in L 2 . The spectral collocation method is used as a solver to obtain the unknown expansion coefficients. The Jacobi roots are used as collocation nodes. Our solutions can easily be a generalization of the solutions of the classical LE equation, by obtaining a numerical solution based on new parameters, by fixing these parameters to the classical case, we obtain the solution of the classical equation. We provide a meticulous convergence analysis and demonstrate rapid convergence of the truncation error concerning the number of retained modes. Numerical examples show the effectiveness and applicability of the method. The primary benefits of the suggested approach are that we significantly reduce the complexity of the underlying differential equation by solving a nonlinear system of algebraic equations that can be done quickly and accurately using Newton’s method and vanishing initial guesses.

1. Introduction

The Lane–Emden equation appears in astrophysics and is a dimensionless form of Poisson’s equation for the gravitational potential of a Newtonian self-gravitating, spherically symmetric, polytropic fluid. This equation was first named in 1870 after astrophysicists Jonathan Homer Lane and Robert Emden [1]. The classical form of the LE equations is
θ ( ξ ) + 2 ξ 1 θ ( ξ ) + θ n ( ξ ) = 0 ,
subject to the standard boundary conditions
θ ( 0 ) = 1 , θ ( 0 ) = 0 ,
where ξ is an independent variable representing the dimensionless radius, θ ( ξ ) is the density function, and thus the pressure ρ = ρ c θ n is the central density, and the number n denotes the polytropic index. This equation has no closed-form solution in general; for the three cases n = 0 , 1 , and 5, we have the following closed-form solutions
θ ( ξ ) = θ 0 ( ξ ) = 1 ξ 2 6 , θ ( ξ ) = θ 1 ( ξ ) = sin ξ ξ , θ ( ξ ) = θ 5 ( ξ ) = 1 1 + 1 3 ξ 2 .
The LE equation describes the dimensionless density distribution in an isothermal gas sphere and plays an important role in galactic dynamics and in the theory of stellar structure and evolution [2]. The LE equation was further investigated with more general nonlinear terms of the form f ( ξ , θ ) in both classical and fractional cases, and recently, many numerical and analytical methods have been applied to handle this LE equation, for instance, the high-order wavelets method [3], the Bernoulli wavelets method [4], the Monte Carlo method [5], the Taylor series method [6], the Morlet wavelets neural network method [7], the Bessel polynomials method [8], the neuro-evolution approach [9], the power series method [10], the implicit method [11], the composite Chebyshev finite difference method [12], the computational approximate method [13], the collocation method [14], and the tau and Galerkin methods [15,16]. Within our study, we use the typical shifted Jacobi collocation method; the collocation method is very important in solving many types of differential problems [17,18,19,20].
The study of the calculus of nonintegral order generalizes the concept of differentiation and integration in the classical sense; it takes into consideration the memory effect of the phenomena, which extensively helps to manipulate real-life phenomena that cannot be precisely characterized by the classical differentiation definition. For the purpose of establishing a formula for the fractional order, Riemann used a generalization of a Taylor series and invented an arbitrary complementary function. This idea motivated Caputo, in 1967, to precisely define the most significant definition of the fractional derivatives, which was named Caputo fractional derivative [21]. Later, some seminal textbooks discussed in-depth these definitions in the theory of fractional differential equations, for instance, [22,23].
Recently, a more general helpful generalization of fractional integral operators was introduced and used to describe the derivatives with more nonlocality properties by introducing a fractional-order integral of a given function related to another function [24,25,26]. This definition and its promising properties allowed researchers to study the availability of spectral methods to handle the LE equation in a more general form by replacing the term θ with the fractional derivative in the fractal–fractional sense and replacing the term θ n with a nonlinear quadratic function in θ .
Orthogonal Jacobi polynomials [27,28,29] have many lucrative fixtures that make them very important in the numerical solution of different kinds of differential problems, primarily via spectral methods. The most important features of Jacobi polynomials are orthogonality, exponential accuracy, and the existence of two parameters that may affect the variety of the approximate solutions, which make these polynomials appropriate for solving diverse problems. In our work, we use the orthonormal Jacobi polynomials, which help us construct a spectral collocation algorithm to handle the nonlinear fractal–fractional LE equation.
It is worth reporting that, within this study, we handled a more general form of the LE equation with a generalized fractional derivative, which helps the interpretation of the solutions with a wide temporal capture of the phenomena of astronomical objects. As far as we know, this is the first time in the literature that the orthonormal shifted Jacobi polynomials are used as basis functions to solve spectrally the LE equation with generalized a fractal–fractional derivative with a detailed study of the truncation error. The main advantages of the proposed methods are that we drastically convert the underlying problem to that of solving a nonlinear system of algebraic equations that can be easily solved via Newton’s method with vanishing initial guesses, with stunted computational time; for recent advances in the field of numerical treatment of related models, the interested readers are referred to [30,31,32].
In bullet points, the main contributions of this study are as follows:
  • We suggest orthonormal Jacobi polynomials as the basis of the solution.
  • We build and prove all derivatives needed within the algorithm.
  • We construct and implement a collocation scheme to handle the nonlinear LE equation with a generalized fractional derivative.
  • We study in detail the truncation error of the method.
  • We perform some numerical examples with comparisons, when possible, with other existing methods.
This manuscript is organized as follows: in Section 2, we report all needed definitions, and state and prove very essential lemmas and theorems. Section 3 is devoted to the structure of the spectral collocation algorithm for handling the generalized fractal–fractional LE equation. In Section 4, we give an upper estimate of the truncation error. In Section 5, we perform some numerical experiments with comparisons to test and validate the method; some concluding remarks are reported in Section 6.

2. Preliminaries

In this section, the essential definitions of the generalized Caputo fractional operators and some relevant properties of the orthonormal normalized Jacobi polynomials are reported, which are subsequently of important use.

2.1. Generalized Caputo Type Fractional Derivative

Definition 1
([24]). Let ρ be a positive constant, the generalized fractional integral of a continuous function φ ( z ) of order α > 0 is defined by
I α , ρ φ ( z ) = ρ 1 α Γ ( α ) 0 z n ρ 1 ( z ρ n ρ ) α 1 φ ( n ) d n .
Definition 2
([24]). Let ρ > 0 , s 1 < α < s , and  s = α ; the generalized Caputo fractional derivative of order α of a function φ ( z ) C s [ 0 , 1 ] is defined as
D α , ρ φ ( z ) = ρ α s + 1 Γ ( s α ) 0 z n ρ 1 ( z ρ n ρ ) s α 1 n 1 ρ d d n s φ ( n ) d n .
The generalized Caputo fractional derivative satisfies:
D α , ρ const . = 0 ,
if α ( s 1 , s ) , > s 1 , then
D α , ρ z ρ = ρ α Γ ( + 1 ) Γ ( + 1 α ) z ρ ( α ) , N 0 , α , 0 , N 0 , < α .
It should be noted here that, in Definition 2, if we set ρ = 1 , we directly obtain the usual Caputo derivative.

2.2. An Account of Shifted Orthonormal Normalized Jacobi Polynomials

The shifted normalized Jacobi polynomials on the interval [ 0 , 1 ] associated with the real parameters ( a , b > 1 , ) are a sequence of polynomials { ϕ m ( a , b ) ( 2 t 1 ) : m = 0 , 1 , 2 , } that can be defined as
ϕ m ( a , b ) ( 2 t 1 ) = m ! Γ ( a + 1 ) Γ ( a + m + 1 ) P m ( a , b ) ( 2 t 1 ) ,
where P m ( a , b ) ( 2 t 1 ) are the shifted classical Jacobi polynomials.
Suppose that
ψ m ( a , b ) ( t ) = h m ϕ m ( a , b ) ( 2 t 1 ) ,
where
h m = ( a + b + 2 m + 1 ) Γ ( a + m + 1 ) Γ ( a + b + m + 1 ) m ! Γ ( a + 1 ) 2 Γ ( b + m + 1 ) .
The orthonormality relation of ψ m ( a , b ) ( t ) is given by
0 1 ψ n ( a , b ) ( t ) ψ m ( a , b ) ( t ) ω ( t ) d t = δ n , m ,
where ω ( t ) = ( 1 t ) a t b and
δ n , m = 1 , if n = m , 0 , if n m .
The power-form representation of ψ m ( a , b ) ( t ) can be represented as
ψ j ( a , b ) ( t ) = r = 0 j ( 1 ) j + r ( a + b + 2 j + 1 ) Γ ( a + b + j + r + 1 ) h j r ! Γ ( a + 1 ) ( j r ) ! Γ ( b + r + 1 ) t r .
Moreover, the inversion formula of ψ m ( a , b ) ( t ) is
t p = r = 0 p B r , p ψ r ( a , b ) ( t ) ,
where
B r , p = ( a + b + 2 r + 1 ) ( a + 1 ) r Γ ( b + p + 1 ) ( p r + 1 ) r Γ ( a + b + r + 1 ) h r r ! Γ ( b + r + 1 ) Γ ( a + b + p + r + 2 ) .
Lemma 1.
For all nonnegative integers i and j, the following linearization formula is valid [33]
ψ i ( a , a + 1 ) ( t ) ψ j ( a , a + 1 ) ( t ) = p = i j i + j h i h j h p χ p , a , i , j ψ p ( a , a + 1 ) ( t ) ,
where
χ p , a , i , j = 2 2 a + 1 i ! j ! Γ ( a + 1 ) π Γ a + 3 2 Γ ( i + 2 a + 2 ) Γ ( j + 2 a + 2 ) × Γ 1 2 ( i + j p + 3 ) + a Γ 1 2 ( i j + p + 3 ) + a Γ 1 2 ( i + j + p + 3 ) + a Γ 1 2 ( i + j + p + 4 a + 4 ) Γ 1 2 ( i + j p + 2 ) Γ 1 2 ( i j + p + 2 ) Γ 1 2 ( i + j + p + 2 ) Γ 1 2 ( i + j + p + 3 ) + a , i f ( i + j p ) e v e n Γ 1 2 ( i + j p ) + a + 1 Γ 1 2 ( i j + p + 2 ) + a Γ 1 2 ( i + j + p + 2 ) + a Γ 1 2 ( i + j + p + 4 a + 5 ) Γ 1 2 ( i + j p + 1 ) Γ 1 2 ( i j + p + 1 ) Γ 1 2 ( i + j + p + 1 ) Γ 1 2 ( i + j + p + 4 ) + a , o t h e r w i s e .
The generalized hypergeometric function is defined by
r F s p 1 , p 2 , , p r q 1 , q 2 , , q s t = n = 0 ( p 1 ) n ( p 2 ) n ( p r ) n ( q 1 ) n ( q 2 ) n ( q s ) n t n n ! .
Theorem 1
([34]). The qth derivative of ϕ n ( a , b ) ( t ) is given explicitly by
D q ϕ n ( a , b ) ( t ) = i = 0 n q λ i , n , q ( a , b ) ϕ i ( a , b ) ( t ) ,
where
λ i , n , q ( a , b ) = 2 q n ! Γ ( a + b + i + 1 ) ( a + b + n + 1 ) q ( a + i + q + 1 ) i + n q ( a + b + n + q + 1 ) i i ! ( i + n q ) ! Γ ( a + b + 2 i + 1 ) ( a + i + 1 ) n i × 3 F 2 i n + q , a + b + i + n + q + 1 , a + i + 1 a + i + q + 1 , a + b + 2 i + 2 1 .
Corollary 1.
The first derivative of ψ i ( a , b ) ( t ) is given explicitly by
d ψ i ( a , b ) ( t ) d t = s = 0 i 1 θ s , i ( a , b ) ψ s ( a , b ) ( t ) ,
where
θ s , i ( a , b ) = i ! h i Γ ( 1 + a + b + s ) ( a + b + i + 1 ) ( a + s + 2 ) s + i 1 ( a + b + i + 2 ) s s ! h s ( s + i 1 ) ! Γ ( 1 + a + b + 2 s ) ( a + s + 1 ) i s × 3 F 2 a + s + 1 , s i + 1 , a + b + s + i + 2 a + s + 2 , a + b + 2 s + 2 1 .
Corollary 2.
The second derivative of ψ i ( a , b ) ( t ) is given explicitly by
d 2 ψ i ( a , b ) ( t ) d t 2 = s = 0 i 2 μ s , i ( a , b ) ψ s ( a , b ) ( t ) ,
where
μ s , i ( a , b ) = i ! h i Γ ( a + b + s + 1 ) ( a + b + i + 1 ) 2 ( a + s + 3 ) s + i 2 ( a + b + i + 3 ) s s ! h s ( s + i 2 ) ! Γ ( a + b + 2 s + 1 ) ( a + s + 1 ) i s × 3 F 2 a + s + 1 , s i + 2 , a + b + s + i + 3 a + s + 3 , a + b + 2 s + 2 1 .
Proof.
The proof of Corollaries 1 and 2 are a direct result from using the definition of ψ i ( a , b ) ( t ) along with Theorem 1.    □
Corollary 3.
The following relation is valid for ρ > 0 and 0 < β < 1
D β , ρ t m = ρ β Γ m ρ + 1 Γ m ρ β + 1 t m β ρ , m 1 .
Proof.
Putting φ ( t ) = t m and s = 1 in Definition 2, one has
D β , ρ t m = m ρ β Γ ( 1 β ) 0 t ( t ρ n ρ ) β n m 1 d n .
Integrating the right-hand side of the previous equation, we get
D β , ρ t m = m ρ β 1 Γ ( 1 β ) β 1 β , m ρ t m β ρ ,
where β ( . , . ) is the well-known beta function. Hence, the last equation can be simplified and written alternatively in the following form
D β , ρ t m = ρ β Γ m ρ + 1 Γ m ρ β + 1 t m β ρ .
This completes the proof of this corollary.    □
Theorem 2.
The following formula holds for 0 < β < 1 ,
D β , ρ ψ i ( a , b ) ( t ) = t ρ β p = 0 i r = p i A r , i B p , r ψ p ( a , b ) ( t ) ( B 0 , 0 A 0 , i ) ψ 0 ( a , b ) ( t ) .
Proof.
Using Corollary 3 along with the power form of ψ i ( a , b ) ( t ) , one has
D β , ρ ψ i ( a , b ) ( t ) = p = β i A p , i t p ρ β ,
where
A p , i = ρ β ( 1 ) i + p ( a + b + 2 i + 1 ) Γ p ρ + 1 Γ ( a + b + i + p + 1 ) p ! Γ ( a + 1 ) h i ( i p ) ! Γ ( b + p + 1 ) Γ p ρ β + 1 .
Now, t p ρ β can be written with the aid of the inversion formula (7) as
t p ρ β = t ρ β r = 0 p B r , p ψ r ( a , b ) ( t ) .
Inserting Equation (22) into Equation (20) yields
D β , ρ ψ i ( a , b ) ( t ) = t ρ β p = β i r = 0 p A p , i B r , p ψ r ( a , b ) ( t ) .
After expanding and rearranging the right-hand side of the last equation, we get the desired result. This completes the proof of the theorem.    □

3. Collocation Approach for the Nonlinear Generalized Fractional LE Equation

Consider the following nonlinear generalized fractional LE equation:
u t t ( t ) + 2 t D β , ρ u ( t ) + a 0 u 2 ( t ) + a 1 u ( t ) = f ( t ) , t ( 0 , 1 ) ,
subject to the initial conditions
u ( 0 ) = y 0 , u ( 0 ) = y 1 ,
where 0 < β 1 , a 0 , a 1 , y 0 , y 1 are known constants, and  f ( t ) is a known continuous source term.
As we know, the set { ψ i ( a , b ) ( t ) , i : 0 , , } forms an orthonormal basis function in the space function L ω ( t ) 2 ( 0 , 1 ) . Consequently, any function u ( t ) L ω ( t ) 2 ( 0 , 1 ) can be written as
u ( t ) = i = 0 c i ψ i ( a , b ) ( t ) ,
and approximated as
u ( t ) u N ( t ) = i = 0 N c i ψ i ( a , b ) ( t ) .
Now, we present our technique when b = a + 1 for the following two cases:
  • The case in which 0 < β < 1 and ρ > 0 .
  • The case in which β = 1 and ρ = 1 .

3.1. The Case in Which 0 < β < 1 and ρ > 0

The application of Lemma 1 and Corollary 2 along with Theorem 2 enables us to write the residual R ( t ) of Equation (24) as
R ( t ) = d 2 u N ( t ) d t 2 + 2 t D β , ρ u N ( t ) + a 0 u N 2 ( t ) + a 1 u N ( t ) f ( t ) = i = 0 N s = 0 i 2 c i μ s , i ( a , a + 1 ) ψ s ( a , a + 1 ) ( t ) + 2 t ρ β 1 i = 0 N p = 0 i r = p i c i A r , i B p , r ψ p ( a , a + 1 ) ( t ) ( B 0 , 0 A 0 , i ) ψ 0 ( a , a + 1 ) ( t ) + a 0 i , j = 0 N p = i j i + j c i h i h j h p χ p , a , i , j ψ p ( a , a + 1 ) ( t ) + a 1 i = 0 N c i ψ i ( a , a + 1 ) ( t ) f ( t ) .
Multiplying the last equation by t 1 + ρ β , we get
R ¯ ( t ) = t 1 + ρ β i = 0 N s = 0 i 2 c i μ s , i ( a , a + 1 ) ψ s ( a , a + 1 ) ( t ) + 2 i = 0 N p = 0 i r = p i c i A r , i B p , r ψ p ( a , a + 1 ) ( t ) ( B 0 , 0 A 0 , i ) ψ 0 ( a , a + 1 ) ( t ) + a 0 t 1 + ρ β i , j = 0 N p = i j i + j c i h i h j h p χ p , a , i , j ψ p ( a , a + 1 ) ( t ) + a 1 t 1 + ρ β i = 0 N c i ψ i ( a , a + 1 ) ( t ) t 1 + ρ β f ( t ) ,
Now, equations t , t 1 + ρ β , and t 1 + ρ β f ( t ) may be approximated as
t = m = 0 1 a m ψ m ( a , a + 1 ) ( t ) , t 1 + ρ β = n = 0 M b n ψ n ( a , a + 1 ) ( t ) , t 1 + ρ β f ( t ) = l = 0 M g l ψ l ( a , a + 1 ) ( t ) ,
where M is an arbitrary positive number and
a m = 0 1 t ψ m ( a , a + 1 ) ( t ) ω ( t ) d t , m = 0 , 1 , b n = 0 1 t 1 + ρ β ψ n ( a , a + 1 ) ( t ) ω ( t ) d t , g l = 0 1 t 1 + ρ β f ( t ) ψ l ( a , a + 1 ) ( t ) ω ( t ) d t .
Inserting Equation (30) into Equation (29), R ¯ ( t ) can be obtained. Hence, the application of the collocation method [35,36,37] enables us to get the following ( N + 1 ) nonlinear algebraic system of equations in the unknown expansion coefficients c i
R ¯ ( t i ) = 0 , i = 1 , 2 , , N 1 , u N ( 0 ) = y 0 , u N ( 0 ) = y 1 ,
where { t i : i = 1 , 2 , , N 1 } are the first ( N 1 ) distinct roots of ψ N + 1 ( a , a + 1 ) ( t ) . Hence, the system (32) can be solved with the aid of the well-known Newton iterative method.

3.2. The Case in Which β = 1 and ρ = 1

Using similar steps as those given in the previous case along with Lemma 1, Corollary 1, and Corollary 2, we get the ( N + 1 ) nonlinear algebraic system of equations in the unknown expansion coefficients c i that can be solved using Newton’s iterative method.
Remark 1.
Algorithm 1 shows all the steps required to obtain the numerical solution of Equation (24) governed by the conditions (25).
  Algorithm 1: Coding algorithm for the proposed scheme.
  Input β , a , ρ , a 0 , a 1 , y 0 , y 1 , N , M and f ( t ) .
  Step 1. Assume an approximate solution u N ( t ) as in (27).
  Step 2. Compute R ¯ ( t ) as in (29).
  Step 3. Apply the collocation method to obtain the system in (32).
  Step 4. Use FindRoot command with initial guess { c i = 10 i , i : 0 , 1 , . . . , N } ,
                  to solve the system in (32) to get c i .
  Output u N ( t )

4. Error Bound

We first define the following error norms:
L 2 [ 0 , 1 ] = { u : u 2 = 0 1 u 2 t b ( 1 t ) a d t 1 2 < } ,
L [ 0 , 1 ] = { u : u = max t [ 0 , 1 ] | u | < } .
Consider the following space functions
V N ( a , b ) = span { ψ i ( a , b ) ( t ) : i = 0 , 1 , , N }
and assume that u N ( t ) V N ( a , b ) is the best approximation of u ( t ) ; then, by the definition of the best approximation, we have
| | u ( t ) u N ( t ) | | | | u ( t ) v N ( t ) | | , v N ( t ) V N ( a , b ) .
It turns out that the previous inequality is also true if v N ( t ) denotes the interpolating polynomial for u ( t ) at points t i , where t i are the roots of ψ i ( a , b ) ( t ) . Then, by similar procedures as in [38]
u ( t ) v N ( t ) = d N + 1 u N ( η ) d t N + 1 ( N + 1 ) ! i = 0 N ( t t i ) ,
where η [ 0 , 1 ] and hence, one has
| | u ( t ) v N ( t ) | | max t [ 0 , 1 ] d N + 1 u N ( η ) d t N + 1 | | i = 0 N ( t t i ) | | ( N + 1 ) ! .
Since u ( t ) is a smooth function on [ 0 , 1 ] , then there exist a constant k, such that
max t [ 0 , 1 ] d N + 1 u N ( t ) d t N + 1 k
To minimize the factor | | i = 0 N ( t t i ) | | , let us use the one-to-one mapping t = 1 2 ( z + 1 ) between the intervals [ 1 , 1 ] and [ 0 , 1 ] to deduce that
min t i [ 0 , 1 ] max t [ 0 , 1 ] i = 0 N ( t t i ) = min z i [ 1 , 1 ] max z [ 1 , 1 ] i = 0 N 1 2 ( z z i ) = 1 2 N + 1 min z i [ 1 , 1 ] max z [ 1 , 1 ] i = 0 N ( z z i ) = 1 2 N + 1 min z i [ 1 , 1 ] max z [ 1 , 1 ] ψ N + 1 ( a , a + 1 ) ( z ) y N ( a , a + 1 ) ,
where y N ( a , a + 1 ) = 2 N ( N + 1 ) h N + 1 Γ ( a + 1 ) Γ ( 2 a + 2 N + 2 ) Γ ( a + N + 2 ) Γ ( 2 a + N + 2 ) is the leading coefficient of ψ N + 1 ( a , a + 1 ) ( z ) and z i are the roots of ψ N + 1 ( a , a + 1 ) ( z ) .
It is known that
max z [ 1 , 1 ] ψ N + 1 ( a , a + 1 ) ( z ) = ψ N + 1 ( a , a + 1 ) ( 1 ) = h N + 1
Therefore, inequality (37) along with Equations (38) and (39) enable us to get the following desired result
| | u ( t ) u N ( t ) | | k h N + 1 2 N + 1 y N ( a , a + 1 ) ( N + 1 ) ! = k Γ ( a + N + 2 ) Γ ( 2 a + N + 2 ) 2 ( N + 1 ) ( N + 1 ) ! Γ ( a + 1 ) Γ ( 2 a + 2 N + 2 ) .
Hence, an upper bound of the absolute error is obtained for the approximate and exact solutions.

5. Illustrative Examples and Comparisons

Before providing the following examples, we would like to mention that all codes were written and debugged using Mathematica 11 on an HP Z420 workstation, with an Intel (R) Xeon(R) CPU E5-1620 3.6 GHz processor, 16 GB RAM DDR3, and 512 GB storage.
Example 1
([39]). Consider the following nonlinear generalized fractional LE equation
u t t ( t ) + 2 t D β , ρ u ( t ) + u ( t ) = 0 , t ( 0 , 1 ) ,
subject to the initial conditions
u ( 0 ) = 1 , u ( 0 ) = 0 ,
where u ( t ) = sin ( t ) t is the exact solution when β = 1 . Figure 1 shows that the approximate solutions have smaller variations for values of ρ and β near the value ρ = β = 1 when N = 12 , M = 20 , and a = 2 .  Table 1 presents a comparison of the L error between our method at ρ = β = a = 1 , M = 20 and the methods in [39,40,41,42,43].
Example 2.
Consider the following nonlinear generalized fractional LE equation
u t t ( t ) + 2 t D β , ρ u ( t ) + u 2 ( t ) + u ( t ) = f ( t ) , t ( 0 , 1 ) ,
subject to the initial conditions
u ( 0 ) = 1 , u ( 0 ) = β ,
where f ( t ) is chosen such that the exact solution is u ( t ) = e β t . Table 2 shows the L 2 and L errors at different values of ρ, β, a, and N when M = 20 . Figure 2 illustrates the absolute error (left) and approximate solution (right) at ρ = 1 , β = 0.8 , a = 1 , N = 12 , and M = 20 .
Example 3.
Consider the following nonlinear generalized fractional LE equation
u t t ( t ) + 2 t D β , ρ u ( t ) + u 2 ( t ) + u ( t ) = t 2 2 ρ β Γ 1 + 4 ρ Γ β + 4 ρ + 1 t 1 ρ β + t 6 + t 2 + 12 , t ( 0 , 1 ) ,
subject to the initial conditions
u ( 0 ) = 0 , u ( 0 ) = 0 ,
where u ( t ) = t 4 is the exact solution. Figure 3 illustrates the absolute errors at different values of β when ρ = 1.01 , a = 2 , N = 4 , and M = 20 . Table 3 shows the L 2 and L errors at different values of ρ, β, and N at a = 1 and M = 20 . Moreover, we have reported the CPU running times in Table 3.
Example 4.
Consider the following nonlinear generalized fractional LE equation
u t t ( t ) + 2 t D β , ρ u ( t ) + 2 u 2 ( t ) + 3 u ( t ) = f ( t ) , t ( 0 , 1 ) ,
subject to the initial conditions
u ( 0 ) = 0 , u ( 0 ) = 0 ,
where f ( t ) is chosen such that the exact solution is u ( t ) = l n 1 + t 2 . Figure 4 shows the maximum absolute errors when β = 0.95 at different values of ρ, a, and N when M = 20 , while Figure 5 illustrates the absolute errors at different values of β, ρ, and a when N = 20 and M = 20 . It can be seen that the approximate solution is quite near to the precise one.
Example 5.
Consider the following nonlinear generalized fractional LE equation
u t t ( t ) + 2 t D β , ρ u ( t ) + u 2 ( t ) + 3 u ( t ) = 0 , t ( 0 , 1 ) ,
subject to the initial conditions
u ( 0 ) = 0 , u ( 0 ) = 0 .
Since the exact solution is not available, we define the following error norm
E = max t [ 0 , 1 ] t u N t t ( t ) + 2 D β , ρ u N ( t ) + t u N 2 ( t ) + 3 t u N ( t ) .
We applied our technique with when N = 8 and M = 20 . The values of E at various values of ρ, β, and a are listed in Table 4. Moreover, we report the CPU running times in that table.

6. Closing Remarks

Within this research work, we presented and analyzed an accurate collocation solver for a specific nonlinear LE equation with a generalized fractal–fractional Caputo derivative. We also discussed the truncation error of the suggested approximate orthonormal Jacobi solution. Some numerical results and comparisons were exhibited to check and verify the validity and the accuracy of the proposed algorithm. We believe that the offered scheme can be extended to more general models in different disciplines in engineering, mathematics, and physics. As an expected future work, we aim to employ the developed theoretical results in this paper along with suitable spectral methods to treat numerically some other types of generalized fractional differential equations; see, e.g., [44,45].

Author Contributions

Conceptualization, Y.H.Y.; methodology, Y.H.Y. and A.G.A.; software, Y.H.Y. and A.G.A.; formal analysis, A.G.A.; data curation, Y.H.Y. and A.G.A.; writing—original draft, A.G.A.; writing—review and editing, Y.H.Y.; supervision, Y.H.Y.; project administration, Y.H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The authors received no funding for this study.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lane, H.J. On the theoretical temperature of the sun, under the hypothesis of a gaseous mass maintaining its volume by its internal heat, and depending on the laws of gases as known to terrestrial experiment. Am. J. Sci. 1870, 2, 57–74. [Google Scholar] [CrossRef] [Green Version]
  2. Liu, F. Polytropic gas spheres: An approximate analytic solution of the Lane-Emden equation. MNRAS 1996, 281, 1197–1205. [Google Scholar] [CrossRef] [Green Version]
  3. Singh, M.; Singh, K. An efficient technique based on higher order Haar wavelet method for Lane–Emden equations. Math. Comput. Simul. 2023, 206, 21–39. [Google Scholar]
  4. Kumbinarasaiah, S.; Manohara, G.; Hariharan, G. Bernoulli wavelets functional matrix technique for a system of nonlinear singular Lane Emden equations. Math. Comput. Simul. 2023, 204, 133–165. [Google Scholar]
  5. El-Essawy, S.H.; Nouh, M.I.; Soliman, A.A.; Rahman, H.I.; Abd-Elmougod, G.A. Monte Carlo simulation of Lane–Emden type equations arising in astrophysics. Astron. Comput. 2023, 42, 100665. [Google Scholar] [CrossRef]
  6. He, J.H.; Ji, F.Y. Taylor series solution for Lane–Emden equation. J. Math. Chem. 2019, 57, 1932–1934. [Google Scholar] [CrossRef]
  7. Sabir, Z.; Wahab, H.A.; Umar, M.; Sakar, M.G.; Raja, M.A. Novel design of Morlet wavelet neural network for solving second order Lane–Emden equation. Math. Comput. Simul. 2020, 172, 1–14. [Google Scholar] [CrossRef]
  8. Izadi, M.; Srivastava, H.M.; Adel, W. An Effective Approximation Algorithm for Second-Order Singular Functional Differential Equations. Axioms 2022, 11, 133. [Google Scholar] [CrossRef]
  9. Sabir, Z.; Raja, M.A.Z.; Botmart, T.; Weera, W. A neuro-evolution heuristic using active-set techniques to solve a novel nonlinear singular prediction differential model. Fractal Fract. 2022, 6, 29. [Google Scholar] [CrossRef]
  10. Awonusika, R.O. Analytical Solutions of a Class of Fractional Lane–Emden Equation: A Power Series Method. Int. J. Appl. Comput. Math. 2022, 8, 155. [Google Scholar] [CrossRef]
  11. Momoniat, E.; Harley, C. Approximate implicit solution of a Lane-Emden equation. New Astron. 2006, 11, 520–526. [Google Scholar] [CrossRef]
  12. Aydinlik, S.; Kiris, A.; Roul, P. An effective approach based on Smooth Composite Chebyshev Finite Difference Method and its applications to Bratu-type and higher order Lane–Emden problems. Math. Comput. Simul. 2022, 202, 193–205. [Google Scholar] [CrossRef]
  13. Bisheh-Niasar, M. A computational method for solving the Lane-Emden initial value problems. Comput. Methods Differ. Equ. 2020, 8, 673–684. [Google Scholar]
  14. Adel, W.; Sabir, Z.; Rezazadeh, H.; Aldurayhim, A. Application of a Novel Collocation Approach for Simulating a Class of Nonlinear Third-Order Lane–Emden Model. Math. Probl. Eng. 2022, 2022, 5717924. [Google Scholar] [CrossRef]
  15. Abd-Elhameed, W.; Ahmed, H.M. Tau and Galerkin operational matrices of derivatives for treating singular and Emden–Fowler third-order-type equations. Int. J. Mod. Phys. C 2022, 33, 2250061. [Google Scholar] [CrossRef]
  16. Abdelhakem, M.; Youssri, Y.H. Two spectral Legendre’s derivative algorithms for Lane-Emden, Bratu equations, and singular perturbed problems. Appl. Numer. Math. 2021, 169, 243–255. [Google Scholar] [CrossRef]
  17. Luo, W.H.; Li, C.; Huang, T.Z.; Gu, X.M.; Wu, G.C. A high-order accurate numerical scheme for the Caputo derivative with applications to fractional diffusion problems. Numer. Funct. Anal. Optim. 2018, 39, 600–622. [Google Scholar] [CrossRef]
  18. Bu, S. A collocation methods based on the quadratic quadrature technique for fractional differential equations. AIMS Math. 2022, 7, 804–820. [Google Scholar] [CrossRef]
  19. Derakhshan, M. The stability analysis and numerical simulation based on Sinc Legendre collocation method for solving a fractional epidemiological model of the Ebola virus. Partial Differ. Equ. Appl. Math. 2021, 3, 100037. [Google Scholar] [CrossRef]
  20. Dadkhah, E.; Shiri, B.; Ghaffarzadeh, H.; Baleanu, D. Visco-elastic dampers in structural buildings and numerical solution with spline collocation methods. J. Appl. Math. Comput. 2020, 63, 29–57. [Google Scholar] [CrossRef]
  21. Caputo, M. Linear models of dissipation whose Q is almost frequency independent—II. Geophys. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  22. Podlubny, I. Fractional differential equations. Math. Sci. Eng. 1999, 198, 41–119. [Google Scholar]
  23. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  24. Odibat, Z.; Baleanu, D. Numerical simulation of initial value problems with generalized Caputo-type fractional derivatives. Appl. Numer. Math. 2020, 156, 94–105. [Google Scholar] [CrossRef]
  25. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef]
  26. Youssri, Y.H. Orthonormal ultraspherical operational matrix algorithm for fractal–fractional Riccati equation with generalized Caputo derivative. Fractal Fract. 2021, 5, 100. [Google Scholar] [CrossRef]
  27. Shen, J.; Tang, T.; Wang, L. Spectral Methods: Algorithms, Analysis and Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011; Volume 41. [Google Scholar]
  28. Ashry, H.; Abd-Elhameed, W.M.; Moatimid, G.M.; Youssri, Y.H. Robust Shifted Jacobi-Galerkin Method for Solving Linear Hyperbolic Telegraph Type Equation. Palestine J. Math. 2022, 11, 504–518. [Google Scholar]
  29. Youssri, Y.H.; Hafez, R.M. Exponential Jacobi spectral method for hyperbolic partial differential equations. Math. Sci. 2019, 13, 347–354. [Google Scholar] [CrossRef] [Green Version]
  30. Marasi, H.; Sharifi, N.; Piri, H. Modified differential transform method for singular Lane-Emden equations in integer and fractional order. TWMS J. Appl. Eng. Math. 2015, 5, 124–131. [Google Scholar]
  31. Zhao, Y.L.; Zhu, P.Y.; Gu, X.M.; Zhao, X.L.; Jian, H.Y. A preconditioning technique for all-at-once system from the nonlinear tempered fractional diffusion equation. J. Sci. Comput. 2020, 83, 10. [Google Scholar] [CrossRef] [Green Version]
  32. Abd-Elhameed, W.M.; Alsuyuti, M.M. Numerical Treatment of Multi-Term Fractional Differential Equations via New Kind of Generalized Chebyshev Polynomials. Fractal Fract. 2023, 7, 74. [Google Scholar] [CrossRef]
  33. Abd-Elhameed, W.M.; Ali, A.; Youssri, Y.H. Newfangled Linearization Formula of Certain Nonsymmetric Jacobi Polynomials: Numerical Treatment of Nonlinear Fisher’s Equation. J. Funct. Spaces 2023, 2023, 9833404. [Google Scholar]
  34. Abd-Elhameed, W.M.; Doha, E.H.; Youssri, Y.H. Efficient spectral-Petrov-Galerkin methods for third-and fifth-order differential equations using general parameters generalized Jacobi polynomials. Quaest. Math. 2013, 36, 15–38. [Google Scholar] [CrossRef]
  35. Atta, A.G.; Moatimid, G.M.; Youssri, Y.H. Generalized Fibonacci operational collocation approach for fractional initial value problems. Int. J. Appl. Comput. Math. 2019, 5, 1–11. [Google Scholar] [CrossRef]
  36. Hafez, R.; Youssri, Y. Shifted Gegenbauer–Gauss collocation method for solving fractional neutral functional-differential equations with proportional delays. Kragujev. J. Math. 2022, 46, 981–996. [Google Scholar] [CrossRef]
  37. Atta, A.G.; Youssri, Y.H. Advanced shifted first-kind Chebyshev collocation approach for solving the nonlinear time-fractional partial integro-differential equation with a weakly singular kernel. Comput. Appl. Math. 2022, 41, 1–19. [Google Scholar] [CrossRef]
  38. Bhrawy, A.H.; Zaky, M.A. A method based on the Jacobi tau approximation for solving multi-term time-space fractional partial differential equations. J. Comput. Phys. 2015, 281, 876–895. [Google Scholar] [CrossRef]
  39. Aghazadeh, N.; Ahmadnezhad, G.; Mohammadi, A. Chebyshev-quasilinearization method for solving fractional singular nonlinear Lane-Emden equations. Commun. Math. 2022, 30, 201–228. [Google Scholar]
  40. Wazwaz, A.M. Adomian decomposition method for a reliable treatment of the Emden–Fowler equation. Appl. Math. Comput. 2005, 161, 543–560. [Google Scholar] [CrossRef]
  41. Shiralashetti, S.C.; Deshi, A.B.; Desai, P.B.M. Haar wavelet collocation method for the numerical solution of singular initial value problems. Ain Shams Eng. J. 2016, 7, 663–670. [Google Scholar] [CrossRef] [Green Version]
  42. Khodier, A.M.M.; Hassan, A.Y. One-dimensional adaptive grid generation. Internat. J. Math. Math. Sci. 1997, 20, 577–584. [Google Scholar] [CrossRef]
  43. Saeed, U. Haar Adomian method for the solution of fractional nonlinear Lane-Emden type equations arising in astrophysics. Taiwan. J. Math. 2017, 21, 1175–1192. [Google Scholar] [CrossRef]
  44. Gu, X.M.; Sun, H.W.; Zhao, Y.L.; Zheng, X. An implicit difference scheme for time-fractional diffusion equations with a time-invariant type variable order. Appl. Math. Lett. 2021, 120, 107270. [Google Scholar] [CrossRef]
  45. Gu, X.M.; Huang, T.Z.; Zhao, Y.L.; Lyu, P.; Carpentieri, B. A fast implicit difference scheme for solving the generalized time–space fractional diffusion equations with variable coefficients. Numer. Methods Partial Differ. Equ. 2021, 37, 1136–1162. [Google Scholar] [CrossRef]
Figure 1. Different solutions of Example 1.
Figure 1. Different solutions of Example 1.
Fractalfract 07 00133 g001
Figure 2. The absolute error (left) and approximate solution (right) of Example 2.
Figure 2. The absolute error (left) and approximate solution (right) of Example 2.
Fractalfract 07 00133 g002
Figure 3. The absolute errors of Example 3.
Figure 3. The absolute errors of Example 3.
Fractalfract 07 00133 g003
Figure 4. The maximum absolute errors of Example 4.
Figure 4. The maximum absolute errors of Example 4.
Fractalfract 07 00133 g004
Figure 5. The absolute errors of Example 4.
Figure 5. The absolute errors of Example 4.
Fractalfract 07 00133 g005
Table 1. Comparison of L error of Example 1.
Table 1. Comparison of L error of Example 1.
N = 64
E A D M [40] E H W C M [41] E H W A G M [42] E H W C A M [43] E C W C Q M [39]Our method at N = 12
8.0733 × 10 3 3.2854 × 10 7 9.1857 × 10 8 9.7345 × 10 6 7.2153 × 10 12 3.33067 × 10 16
Table 2. L 2 and L errors of Example 2.
Table 2. L 2 and L errors of Example 2.
a N ρ = 1.2 β = 0.1 ρ = 0.9 β = 0.5 ρ = 0.7 β = 0.9
L 2 Error L Error L 2 Error L Error L 2 Error L Error
42.63416 × 10 10 1.53419  × 10 9 3.56771 × 10 7 4.94507 × 10 6 5.18941 × 10 6 7.99767 × 10 5
181.61411 × 10 15 6.39142 × 10 14 7.61436 × 10 14 3.64833 × 10 13 9.68827 × 10 12 1.42884 × 10 10
121.76776 × 10 15 1.01641 × 10 13 5.56363 × 10 16 3.31957 × 10 14 2.09734 × 10 16 8.46545 × 10 15
41.69006 × 10 10 1.42379 × 10 9 2.37063 × 10 7 5.65026 × 10 6 3.3588 × 10 6 9.01696 × 10 5
281.20331 × 10 17 1.94289 × 10 16 7.22686 × 10 14 6.96721 × 10 13 9.06465 × 10 12 2.20747 × 10 10
121.19946 × 10 17 1.76942 × 10 16 1.57298 × 10 17 1.94289 × 10 16 2.76152 × 10 17 5.82867 × 10 16
49.71126 × 10 11 1.75633 × 10 11 1.41421 × 10 7 5.98682 × 10 6 1.99222 × 10 6 9.46116 × 10 5
385.09106 × 10 18 1.11022 × 10 16 5.46102 × 10 14 1.07617 × 10 12 6.80349 × 10 12 3.10428 × 10 10
121.77111 × 10 17 2.89699 × 10 16 2.21721 × 10 17 3.88578 × 10 16 4.23896 × 10 18 6.38378 × 10 16
Table 3. L 2 and L errors of Example 3.
Table 3. L 2 and L errors of Example 3.
N ρ = 0.8 β = 0.85 ρ = 0.5 , β = 0.95
L 2 ErrorCPU Time L ErrorCPU Time L 2 ErrorCPU Time L ErrorCPU Time
408.972.78666 × 10 14 8.92307.8474.32987 × 10 15 7.831
63.86623 × 10 18 8.4413.19189 × 10 15 8.0513.50141 × 10 18 9.9838.88178 × 10 16 9.437
Table 4. Residual error of Example 5.
Table 4. Residual error of Example 5.
a ρ = 0.5 , β = 0.5 ρ = 0.7 , β = 0.7 ρ = 0.9 β = 0.9
ECPU TimeECPU TimeECPU Time
14.14113  × 10 24 15.1325.24048  × 10 25 10.8723.75029  × 10 11 16.706
32.99225 × 10 17 15.8812.20543 × 10 20 16.6937.55597 × 10 25 16.069
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Youssri, Y.H.; Atta, A.G. Spectral Collocation Approach via Normalized Shifted Jacobi Polynomials for the Nonlinear Lane-Emden Equation with Fractal-Fractional Derivative. Fractal Fract. 2023, 7, 133. https://doi.org/10.3390/fractalfract7020133

AMA Style

Youssri YH, Atta AG. Spectral Collocation Approach via Normalized Shifted Jacobi Polynomials for the Nonlinear Lane-Emden Equation with Fractal-Fractional Derivative. Fractal and Fractional. 2023; 7(2):133. https://doi.org/10.3390/fractalfract7020133

Chicago/Turabian Style

Youssri, Youssri Hassan, and Ahmed Gamal Atta. 2023. "Spectral Collocation Approach via Normalized Shifted Jacobi Polynomials for the Nonlinear Lane-Emden Equation with Fractal-Fractional Derivative" Fractal and Fractional 7, no. 2: 133. https://doi.org/10.3390/fractalfract7020133

APA Style

Youssri, Y. H., & Atta, A. G. (2023). Spectral Collocation Approach via Normalized Shifted Jacobi Polynomials for the Nonlinear Lane-Emden Equation with Fractal-Fractional Derivative. Fractal and Fractional, 7(2), 133. https://doi.org/10.3390/fractalfract7020133

Article Metrics

Back to TopTop