A Cubic Spline Collocation Method to Solve a Nonlinear Space-Fractional Fisher’s Equation and Its Stability Examination

: This article seeks to show a general framework of the cubic polynomial spline functions for developing a computational technique to solve the space-fractional Fisher’s equation. The presented approach is demonstrated to be conditionally stable using the von Neumann technique. A numerical illustration is given to demonstrate the proposed algorithm’s effectiveness. The novelty of the present work lies in the fact that the results suggest that the presented technique is accurate and convenient in solving such problems.


Introduction
Fractional calculus is a generalization of classical calculus that deals with fractional order differentiation and integration operations. Everyone who has studied elementary calculus is familiar with the differentiation operator, D = d/dx. Additionally, the n th derivative of a suitable function, D n f (x) = d n f (x)/d n x, is properly defined as long as n is a positive integer. In 1695, L'Hopital enquired to Leibniz about the meaning of D n f if n was fractional. Many prominent mathematicians, including Euler, Laplace, Fourier, Abel, Liouville, Riemann, and Laurent, have since studied time-fractional calculus. However, it was not until 1884 that the theory of generalized operators progressed to the point where it could be used as a starting point for modern mathematicians. Fractional calculus has been utilized in thermodynamics, viscoelasticity, bioengineering, aerodynamics, control theory, electromagnetics, finance, chemistry, and signal processing in recent decades [1][2][3]. Fractional calculus theory dates back to Leibniz in the sixteenth century, and many different types of fractional operators have been created and developed virtually as far as the classical theory has evolved. The Riemann-Liouville, Caputo, conformable, and Riesz operator approaches are among these [4][5][6][7][8][9]. Fractional differential equations (FDEs) appear in engineering processes, a number of physical phenomena, financial products, and biological systems, such as non-exponential relaxation patterns and anomalous diffusion [10]. Several numerical methods have been developed to obtain an approximate solution of FDEs. These methods include the Adomian decomposition method [11], finite difference method [12], variational iteration method [13], spectral methods [14], and homotopy perturbation method [15,16].
Fisher [17] introduced the nonlinear fractional Fisher equation as a mathematical model to describe the kinetic advancing rate of an advantageous gene. The fractional Fisher's equation is also used in autocatalytic chemical reactions, nuclear reactor theory, neurophysiology, flame propagation, and branching Brownian motion processes [18][19][20]. In the open literature, many works describe solving the nonlinear fractional Fisher equation, see [20][21][22][23]. In our article, we propose a cubic polynomial spline-based technique for obtaining approximation solutions for the space-fractional Fisher's equation in the following form [24]: subject to the conditions: and where α is a parameter in the Caputo sense that describes the fractional order of the space derivative, and 1 < α ≤ 2; if we take α = 2, Equation (1) reduces to the classical nonlinear Fisher's equation. Several scholars have studied and solved the space-fractional nonlinear Fisher's equation using various numerical techniques, for example, the generalized differential transform method (GDTM) and the variational iteration method (VIM) [25], the radial basis functions method [26], the quadratic polynomial spline-based method (QPSM) [24], and the Legendre spectral collocation method [27]. The above-mentioned articles show how to solve the space-fractional Fisher's equation in which the exact solution is unknown, but, as far as we are aware, no papers are using cubic polynomial spline for space-fractional Fisher's equation; hence, that is the motivation for our study.
Definition 1. The Caputo space-fractional derivative operator of order α > 0 is given by [25]: where Γ( ) is the gamma function. The right and left Caputo fractional derivatives of order α are given, respectively, by [28]: The numerical solutions of fractional differential equations using numerical methods with spline functions and their use is an important area of research in numerical analysis, see [29][30][31][32][33][34]. Such problems appear very often in physics, engineering, medicine, biology, and other sciences. Spline functions can be integrated and differentiated due to being piecewise polynomials and, since they have a basis with small support, many of the integrals that occur in numerical methods are zero. Numerical methods with spline functions in obtaining the numerical solution of the differential equations lead to band matrices with a low computational cost. Cubic splines were used instead of using linear or quadratic splines, because they produce smooth approximations while maintaining a low computational cost. [35,36]:

Definition 2. Polynomial Spline Function
Let a = x 0 < · · · < x n−1 < x n = b be a subdivision of the interval [a, b] and m ∈ N. A function P : [a, b] → R is called a spline of degree m for this subdivision if P is (m − 1) times continuously differentiable on [a, b] and if the restriction of P to each subinterval [x i , x i+1 ] for i = 0, 1, · · · , n − 1 reduces to a polynomial of degree at most m. Thus, the polynomial spline function has the form: such that each polynomial spline segment is given by: The remainder of this article is structured as follows. In Section 2, a method is proposed that is based on the cubic spline polynomial method. In Section 3, the von Neumann method is used to discuss stability theoretically. In Section 4, we use a numerical example to demonstrate the suggested method's efficiency and accuracy. The conclusion is presented in Section 5.

Derivation of the Method
To approximate u(x, t) by a cubic polynomial spline-based method over finite elements, let the region M = [a, b] × [0, T] be discretized by a set of points M ij , which are the vertices of a grid of points N for each i = 0, 1, . . . , N and t j = jk, k = ∆t for j = 0, 1, . . . , R.

Remark 1: For Equation
and can be obtained by expanding this equation in terms of u x i , t j and its derivatives in the Taylor series, as follows: Then, the last expression T * j i can be written as: . From the expression (21) of the local truncation error, our scheme (17) The fractional Fisher's Equation (1) can now be defined as follows in terms of S j i : Using the finite difference method, we obtain: This allows us to express S j i−1 , S j i , and S j i+1 as follows: where σ j i = 1 − W j i . Substituting (22) into (17) gives us the following system: where: . System (23) consists of N − 1 equations in the unknowns W i , i = 0, 1, . . . , N. Two additional equations are required to solve this system. The boundary conditions in (2) are used to derive these equations, which can be represented as: Writing (23) and (24) in matrix form gives: where: , and, where Φ * and Φ are (N + 1) × (N + 1) matrices and W j and r j are (N + 1) vectors. Now, the initial condition, u(x, t 0 ) = g(x) for a ≤ x ≤ b, implies that W 0 i for i = 0, 1, . . . , N. These values can be utilized in (25) to obtain the value of W 1 i for i = 0, 1, . . . , N. If the technique is reapplied once all the approximations W 1 i are determined, the values of W 2 i , W 3 i , . . . , W R i can be derived in the same way. The nonlinear term in the system is linearized using the technique in (25). For example, if j = 1, we approximate σ 1 i by σ 1 * i computed from w 0 i , which obtains an approximation to w 1 i . For j = m, we approximate σ m i by σ m * i computed from w m−1 i and obtain an approximation to w m i .

Stability Analysis
The von Neumann technique will be used to investigate our system's stability (23). To do this, we must linearize the nonlinear term u(x, t)(1 − u(x, t)) of Fisher's Equation (1), by assuming that the corresponding quantities σ j i−1 , σ j i , and σ j i+1 are equal to a local constant λ * in (23). The key part of the von Neumann analysis is to assume a solution of the form [39]: where q = √ −1, h is the element size, ϕ is the mode number, and ξ j is the amplification factor at time level j. As j increases (more time steps are computed), adding the latter W j i expression to the system (23) gives the following characteristic equation: where: After simple basic calculations, (27) becomes: where φ = ϕh. After using Euler's formula exp(qφ) = cos φ + q sinφ, the last equation (28) can be simplified as: Equation (29) can be rewritten in the forms: or The quantity (1 − cos φ) is positive or equal to zero but 2β cos φ + γ is surely positive if we choose γ > 0 and β > 0 such that γ > 2β. If we choose k, β, and γ small enough to make λ * (2βk cos(φ) + γk) → 0 , then last (31) is close to: For stability, we must have |ξ| ≤ 1 (otherwise ξ j in (26) would expand in an unbounded manner). This condition is satisfied for γ > 0 and β > 0 such that γ > 2β. Finally, we can say that our numerical scheme is conditionally stable for γ > 0, β > 0, and γ > 2β, if k, β, and γ values chosen are small enough.

Numerical Results
In this part, we employ the suggested method for solving the nonlinear space-fractional Fisher's equation [17,18] with α = 1.5, f (x, t) = −x 2 , and u(x, 0) = x and subject to conditions: The solution of the space-fractional Fisher's equation, like that of most fractional partial differential equations, is unknown, so we will compare it with the known computational approaches such as GDTM and VIM. We used k = 0.0005 and N = 80 in all our numerical results. Tables 1-3 compare our method, which was improved in Section 2, to other existing methods.

Conclusions
The nonlinear space-fractional Fisher's equation is numerically treated in this article using a collocation method with cubic polynomial spline functions. The proposed method is demonstrated to be conditionally stable using the von Neumann stability technique. Compared to other useful techniques, such as VIM, GDTM, and QPSM, the numerical results illustrate that our proposed method retains efficiency, capability, and good versatility. The findings of our study are hoped to pique the attention of academics interested in using the cubic polynomial spline to numerically solve nonlinear fractional partial differential equations of the same type. It is proven that our scheme's local truncation error is to O(h α ), 1 < α ≤ 2. Furthermore, the novelty of the present work lies in the fact that the results illustrate that the proposed technique is convenient, accurate, and very efficient in solving such problems.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.