Next Article in Journal
New Classifier Ensemble and Fuzzy Community Detection Methods Using POP Choquet-like Integrals
Previous Article in Journal
Evolution Analysis of Strain Waves for the Fractal Nonlinear Propagation Equation of Longitudinal Waves in a Rod
Previous Article in Special Issue
Lyapunov-Type Inequalities for Systems of Riemann-Liouville Fractional Differential Equations with Multi-Point Coupled Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Stability, Global Stability, and Simulations in a Fractional Discrete Glycolysis Reaction–Diffusion Model

1
Department of Mathematics, Faculty of Science, Al Zaytoonah University of Jordan, Amman 11733, Jordan
2
Laboratory of Dynamical Systems and Control, University of Oum EL-Bouaghi, Oum El Bouaghi 04000, Algeria
3
Department of Mathematics, Faculty of Science, The Hashemite University, Zarqa 13133, Jordan
4
Data Science and Artificial Intelligence Department, Zarqa University, Zarqa 13110, Jordan
5
Department of Mathematics, Faculty of Science and Technology, Irbid National University, Irbid 21110, Jordan
6
Department of Mathematics and Computer Science, University of Oum EL-Bouaghi, Oum El Bouaghi 04000, Algeria
*
Author to whom correspondence should be addressed.
Submission received: 29 May 2023 / Revised: 2 July 2023 / Accepted: 24 July 2023 / Published: 29 July 2023
(This article belongs to the Special Issue Stability Analysis for Fractional-Order Equations)

Abstract

:
In the last few years, reaction–diffusion models associated with discrete fractional calculus have risen in prominence in scientific fields, not just due to the requirement for numerical simulation but also due to the described biological phenomena. This work investigates a discrete equivalent of the fractional reaction–diffusion glycolysis model. The discrete fractional calculus tool is introduced to the discrete modeling of diffusion problems in the Caputo-like delta sense, and a fractional discretization diffusion model is described. The local stability of the equilibrium points in the proposed discrete system is examined. We additionally investigate the global stability of the equilibrium point by developing a Lyapunov function. Furthermore, this study indicates that the L1 finite difference scheme and the second-order central difference scheme can successfully preserve the characteristics of the associated continuous system. Finally, an equivalent summation representing the model’s numerical formula is shown. The diffusion concentration is further investigated for different fractional orders, and examples with simulations are presented to corroborate the theoretical findings.

1. Introduction

Reaction–diffusion equations have recently gained popularity due to their applicability to a wide range of biological and chemical processes as well as their unique characteristics in terms of the propagation speed and the new oscillations not seen in standard two-species reaction–diffusion systems. The biochemical process that converts glucose to pyruvate is known as glycolysis, which is present in the vast majority of living creatures today [1]. Throughout this process, two molecules of pyruvate are generally created for every molecule of glucose. The prevalence of glycolysis in a wide range of animals on Earth suggests that it is a primitive metabolic route [2]. Glycolysis is a series of reactions that divide pyruvates (three-carbon molecules) into two components in order to obtain energy from glucose. The same reaction strategy underpins all glycolysis models. The model differs due to differences in the mechanism for the main enzyme reaction. Because such systems have spatial variables, they are widely utilized to analyze certain unusual behaviors, which can include self-replicating spikes, self-excitation, and spatiotemporal chaos. The need to comprehend dynamics is emphasized for this model in order to actualize a broader range of real-world scenarios. Numerous significant attempts have been made to explore the behavior of the reaction–diffusion glycolysis model. The application of a linear control rule to synchronize the reaction–diffusion glycolysis model is discussed in [3]. In [4], the numerical solution of a general glycolysis reaction–diffusion system was studied. In addition, the stability of the Selkov–Schnakenberg reaction–diffusion system was addressed in [5].
In modern-world experiences, fractional-order nonlinear equations are often employed to simulate a variety of physical phenomena [6,7,8,9]. Due to its numerous applications in the fields of chemistry, biology, physics, technology, thermal technology, electricity, mechanics, data processing, electrohydraulics, population modeling, machine learning, and control, the concept of fractional calculus continues to attract the attention of scientists [10,11,12,13]. Due to its use in and relevance to many fields of research and engineering, fractional reaction–diffusion models have gained attention in recent years. In this field, for example, in [14], a pair fractional reaction–diffusion model was explored. In fractional reaction–diffusion systems, they examined nonlinear oscillations and the stability domain. In [15], a solution to a specific fractional reaction-diffusion problem was proposed. In addition, several breakthroughs have been made in the research on the dynamics of the fractional reaction diffusion glycolysis model, for example, see [16,17,18]. In the works described earlier, the fractional component is dedicated to the time partial derivative; nevertheless, in some circumstances, the standard Laplacian cannot imitate the diffusion, and the mean square movement of the diffusing components does not progress linearly over time. Following numerous findings in biological systems, emphasis has been placed on fractional subdiffusion, which is defined as the mean square movement of a population propagating as a sublinear power law in time. It is now widely known that if particles are confined for arbitrarily extended periods of time, subdiffusion occurs (see [19,20,21]).
Discrete reaction–diffusion systems have received less attention than continuous reaction–diffusion systems. Discrete reaction–diffusion systems are formed of cells or lattices that contain chemical processes and are linked together via mass transfer across their membranes. Discrete reaction diffusion systems are more tightly linked to biological systems than continuous systems, and in reaction to a small disturbance, they may exhibit extremely intriguing and diverse patterns in a two-dimensional space [22,23,24]. However, fractional discrete reaction diffusion equations have received little attention [25]. In [26], a fractional discrete diffusion equation was proposed. The chaotic behavior of a variable-order fractional diffusion equation on discontinuous time scales was examined in [27]. In [28], the global dynamics for a class of discrete fractional epidemic models with reaction–diffusion were also investigated. Apparently, we have a gap in our understanding of the dynamics of such systems.
In recent decades, both integer and noninteger fractional orders of the well-known glycolysis reaction diffusion model have received extensive study. The majority of this work has focused on the dynamics and behavior of solutions ([29,30,31,32]). However, to the authors’ knowledge, this is the first investigation of the dynamics and stability of the discrete fractional reaction–diffusion glycolysis model. This prompted us to study the stability of the equilibrium state in a discrete fractional-order reaction–diffusion glycolysis model. A discrete fractional glycolisis reaction diffusion model in the Caputo-like delta sense is described in this study, employing a second-order central difference scheme and an L1 finite difference scheme. Both the local and global stability were explored, and simulations were conducted to test the model’s applicability.
The following is a summary of the paper’s content. Section 2 provides an introduction to fractional discrete calculus, as well as the fundamental concepts associated with the delta Caputo h-discrete difference operator. The mathematical model under examination is described in Section 3. The analysis of the local stability of the equilibrium state for the presence and absence of diffusion is highlighted in Section 4. The global stability of the equilibrium state, which relies on the fractional orders of the investigated model, is discussed in Section 5. Section 6 includes numerical simulations that back up the findings. Section 7 outlines a number of results.

2. Preliminaries

This section provides the necessary terminology and stability theory for the fractional discrete calculus.
Definition 1
([33]). Given that z : N R , the forward difference operator Δ is stated by
Δ z ( r ) = z ( r + 1 ) z ( r ) ; r N .
Also, the recursive identification of the operator Δ n , n = 1 , 2 , 3 , is Δ n z ( r ) = Δ ( Δ n 1 z ( r ) ) , r N .
In particular, the second-order difference operator of the function z ( t ) is provided by
Δ 2 z ( r ) = z ( r + 2 ) 2 z ( r + 1 ) + z ( r ) .
Theorem 1
([33]). We have the summation by parts formulae for the following two functions: z ; y : R R and a ; b N ; a < b ;
j = a b 1 z ( j ) Δ y ( j ) = z ( j ) y ( j ) | a b j = a b 1 y ( j + 1 ) Δ z ( j ) ,
j = a b 1 z ( j + 1 ) Δ y ( j ) = z ( j ) y ( j ) | a b j = a b 1 y ( j ) Δ z ( j ) .
Definition 2
([34,35]). Let z ( N ) a R . For a given ϑ > 0 , the ϑ -th order ℏ-sum is given by
Δ a ϑ z ( t ) = Γ ( ϑ ) t ϑ s = a ( t σ ( s ) ) ( ϑ 1 ) z ( s ) , σ ( s ) = ( s + 1 ) , t ( N ) a + ϑ ,
the initial value is a R , and the ℏ-falling factorial function is specified by
t ( ϑ ) = ϑ Γ ( t + 1 ) Γ ( t + 1 ϑ ) ,
while
( N ) a + ϑ = { a + ( 1 ϑ ) , a + ( 2 ϑ ) , } .
Definition 3
([36,37]). Considering a function z ( t ) given on ( N ) a and for a particular ϑ > 0 , the Caputo ℏ-difference operator is represented by
C Δ a ϑ z ( t ) = Δ a ( n ϑ ) Δ n z ( t ) ,
where Δ n z ( t ) = z ( t + ) z ( t ) .
Lemma 1
([34]). The following inequality holds
C Δ a ϑ z 2 ( t ) 2 z ( t + ϑ ) C Δ a ϑ z ( t ) , t ( N ) a + ϑ ,
where 0 < ϑ 1 .
Let us consider the nonlinear fractional-order difference system:
C Δ a ϑ z ( t ) = ψ ( t + ϑ , z ( t + ϑ ) ) , t ( N ) a + ϑ .
Theorem 2
([38]). Assume that the equilibrium point of (7) is z . z is asymptotically stable if all of the eigenvalues of ψ ( z ) are found in S ϑ , where
S ϑ = w C : | A r g ( w ) | > ϑ π 2 or | w | > 2 ϑ ϑ cos ϑ A r g ( w ) ϑ .
Theorem 3
([34]). Assume that the equilibrium point of the system (21) is z = 0 . If there is a positively definite and decreasing scalar function, the equilibrium point is asymptotically stable, such that C Δ a ϑ V ( t , z ( t ) ) 0 .

3. The Fractional Discrete Glycolysis Reaction–Diffusion System

The model under consideration is approximated employing two commonly used approaches in this section. This is, to our knowledge, the first discrete model in the literature. The glycolysis reaction–diffusion system, as previously stated, was provided as a model. Refs. [4,5] cite it as follows:
z t = d 1 Δ z + b w z + z 2 w , z Ω , t > 0 , w t = d 2 Δ w + a b w z 2 w , z Ω , t > 0 , z = w = 0 , z Ω , t > 0 , u ( z , 0 ) = z 0 ( z ) > 0 , w ( z , 0 ) = w 0 ( z ) > 0 , z Ω .
In this case, the variables z ( z , t ) and w ( z , t ) denote the chemical concentrations, d 1 and d 2 are the diffusion coefficients, a is the dimensionless input flow, and b is the dimensionless constant rate for the low activity state. Ω is a bounded domain in R n with sufficiently smooth boundaries. For a two-cell system, this model has some stability findings (see [39,40]). The Sel’klov model, which has attracted a lot of interest recently ([41,42,43]), is known to exist for b = 0 .
Researchers have investigated time-fractional systems extensively; hence, the following fractional-time Glycolysis reaction-diffusion system was presented:
0 C D t δ z d 1 Δ z = b w z + z 2 w , z Ω , t > 0 0 C D t δ w d 2 Δ w = a b w z 2 w , z Ω , t > 0 ,
where Δ = i = 1 n 2 z i 2 , 0 < δ 1 is the fractional order, and 0 C D t δ describes the Caputo fractional derivative. d 1 , d 2 , and σ are strictly positive constants with the same initial conditions and Neumann boundary conditions.
Relying on the model (10) and the method of discretization employed in [26], provided that z [ 0 , L ] , we obtain z i + 1 = z i + k , i = 0 , , m , and by applying (1), 2 z ( z , t ) z 2 and 2 w ( z , t ) z 2 may be approximated as
2 z ( z , t ) z 2 z i + 1 ( t ) 2 z i ( t ) + z i 1 ( t ) k 2 , 2 w ( z , t ) z 2 w i + 1 ( t ) 2 w i ( t ) + w i 1 ( t ) k 2 .
By defining the second-order difference operator of z i and w i , we obtain
2 z ( z , t ) z 2 Δ 2 z i 1 ( t ) k 2 , 2 w ( z , t ) z 2 Δ 2 w i 1 ( t ) k 2 .
Thus, we introduce the following discrete-time reaction–diffusion fractional glycolysis system:
C Δ t 0 ϑ z i ( t ) = d 1 k 2 Δ 2 z i 1 ( t + ϑ ) + b w i ( t + ϑ ) z i ( t + ϑ ) + z i 2 ( t + ϑ ) w i ( t + ϑ ) , C Δ t 0 ϑ w i ( t ) = d 2 k 2 Δ 2 w i 1 ( t + ϑ ) + a b w i ( t + ϑ ) z i 2 ( t + ϑ ) w i ( t + ϑ ) ,
with the periodic boundary conditions
z 0 ( t ) = z m ( t ) , z 1 ( t ) = z m + 1 ( t ) , w 0 ( t ) = w m ( t ) , w 1 ( t ) = w m + 1 ( t ) ,
and the initial condition
z i ( t 0 ) = ϕ 1 ( z i ) 0 , w i ( t 0 ) = ϕ 2 ( z i ) 0 .

4. Local Stability

To analyze the asymptotic stability of the previously stated discrete fractional glycolysis model, we consider the single equilibrium produced by this system:
d 1 k 2 Δ 2 z + b w z + z 2 w = 0 , d 2 k 2 Δ 2 w + a b w z 2 w = 0 .
Considering Theorem 1, the unique equilibrium point of system (11) is provided by
( z , w ) = a , a a 2 + b .

4.1. Local Stability of the Free Diffusions System

In this section, we provide sufficient conditions for the system’s local asymptotic stability.
C Δ t 0 ϑ z ( t ) = b w ( t + ϑ ) w ( t + ϑ ) + u 2 ( t + ϑ ) w ( t + ϑ ) , C Δ t 0 ϑ w ( t ) = a b z ( t + ϑ ) z 2 ( t + ϑ ) w ( t + ϑ ) .
Theorem 4.
The conditions that follow must be satisfied for system (15) to be locally asymptotically stable at the positive steady state ( z , w ) :
  • If a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 4 a 2 + b 2 and a 2 b 2 a 2 + b 2 < a 2 + b 2 .
  • If a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 < 4 a 2 + b 2 and a 2 b 2 a 2 + b 2 a 2 + b 2 .
Proof. 
The characteristic equation for the eigenvalues is discovered around this stable state using linear stability analysis:
J ( z , w ) = ψ z ψ w Ψ z Ψ w = 2 a 2 a 2 + b 2 1 a 2 + b 2 2 a 2 a 2 + b 2 ( a 2 + b 2 ) ,
with
ψ ( z , w ) = b w ( t + ϑ ) z ( t + ϑ ) + z 2 ( t + ϑ ) w ( t + ϑ ) ,
and
Ψ ( z , w ) = a b w ( t + ϑ ) z 2 ( t + ϑ ) w ( t + ϑ ) .
From the Jacobian matrix, we can see that
tr ( J ) = 2 a 2 a 2 + b 2 ( a 2 + b 2 ) 1 , det ( J ) = a 2 + b 2 ,
and J ( z , w ) has the characteristic equation given by
Λ 2 tr ( J ) Λ + det ( J ) = 0 ,
or its discriminant is
N Λ = tr 2 ( J ) 4 det ( J ) = a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 4 a 2 + b 2 .
As a consequence, we may conclude the following according to the sign of (21).
  • If N Λ > 0 , it is evident that det ( J ) > 0 . Thus, the sign of tr ( J ) determines the eigenvalues’ negativity, and the eigenvalues Λ 1 and Λ 2 are real and may be expressed as
    Λ 1 = tr ( J ) N Λ 2 , Λ 2 = tr ( J ) + N Λ 2 .
    -
    If tr ( J ) < 0 , we have
    Λ 1 = tr ( J ) N Λ 2 < 0 .
    Hence, A r g ( Λ 1 ) = π . The two eigenvalues are real; therefore, A r g ( Λ 1 ) = A r g ( Λ 2 ) = π . As a result, ( z , w ) is asymptotically stable according to Theorem 2.
    -
    If tr ( J ) > 0 , then we have
    Λ 2 = tr ( J ) + N Λ 2 > 0 .
    Therefore, A r g ( Λ 2 ) = 0 , and based on Theorem 2, system (15) is unstable.
  • If N Λ < 0 , then
    Λ 1 = tr ( J ) i N Λ 2 , Λ 2 = tr ( J ) + i N Λ 2 .
    We may examine the solutions relying on the sign of tr ( J ) .
    -
    If tr ( J ) < 0 or tr ( J ) > 0 , consequently, system (15) is asymptotically stable following the same case that was previously analyzed.
    -
    If tr ( J ) = 0 , then
    A r g i N Λ 2 = A r g i N Λ 2 = π 2 ,
    and system (15) is asymptotically stable.
  • If N Λ = 0 , it is impossible for tr ( J ) to equal zero, since det ( J ) > 0 . The sign of the eigenvalues and that of tr ( J ) are the same. As a consequence, ( z , w ) is unstable if tr ( J ) > 0 , and it is asymptotically stable for all ϑ ( 0 , 1 ] , if tr ( J ) < 0 .
The proof is complete. □

4.2. Local Stability of the Diffusion System

The steady state ( z , w ) may be stable in the presence of diffusion under specific parameter conditions. We follow the same procedure as in [44], beginning by examining the eigenvalues of this equation:
Δ 2 z i 1 ( t + ϑ ) + Λ i z i ( t + ϑ ) = 0 ,
with the periodic boundary conditions:
z 0 ( t ) = z m ( t ) , z 1 ( t ) = z m + 1 ( t ) .
We obtain
C Δ t 0 ϑ z i ( t ) = d 1 k 2 Λ i z i ( t + ϑ ) + b w i ( t + ϑ ) z i ( t + ϑ ) + z i 2 ( t + ϑ ) w i ( t + ϑ ) , C Δ t 0 ϑ w i ( t ) = d 2 k 2 Λ i w i ( t + ϑ ) + a b w i ( t + ϑ ) z i 2 ( t + ϑ ) w i ( t + ϑ ) .
When linearizing system (28) about the steady state ( z , w ) , we obtain the following.
J i = d 1 k 2 Λ i + 2 a 2 a 2 + b 2 1 a 2 + b 2 2 a 2 a 2 + b 2 d 2 k 2 Λ i ( a 2 + b 2 ) .
Then, the following result is obtained.
Theorem 5.
We suppose that
a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 > 4 a 2 + b 2 and a 2 b 2 a 2 + b 2 < a 2 + b 2 ;
system (11) is asymptotically stable at the steady state ( z , w ) , if the following hold:
  • If d 1 < d 2 and d 1 k 2 Λ 1 2 a 2 a 2 + b 2 1 .
  • If d 1 > d 2 and d 1 k 2 Λ 1 2 a 2 a 2 + b 2 1 ,
    and in addition, the eigenvalues
    μ j ( Λ i ) = t r ( J i ) + ̲ t r ( J i ) 2 4 d e r t ( J i ) 2 , j = 1 , 2 ,
    satisfy A r g ( μ j ( Λ i ) ) > ϑ π 2 .
Proof. 
We linearize the system to investigate its local asymptotic stability. By applying basic linear operator theory and taking the system’s fractional nature into account, we may argue that ( u , v ) is asymptotically stable if the eigenvalues of the linearized system satisfy the conditions of Theorem 2. We obtain
J i = 6 d 1 k 2 Λ i + 2 a 2 a 2 + b 2 1 a 2 + b 2 2 a 2 a 2 + b 2 6 d 2 k 2 Λ i ( a 2 + b 2 ) = J μ ( Λ i ) I ,
which has the eigenvalue equation
μ 2 ( Λ i ) tr ( J i ) μ ( Λ i ) + det ( J i ) = 0 ,
where
tr ( J i ) = d 1 k 2 + d 2 k 2 Λ i + a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) ,
and
det ( J i ) = d 1 k 2 d 2 k 2 Λ i 2 + d 1 k 2 ( a 2 + b 2 ) d 2 k 2 a 2 b 2 a 2 + b 2 Λ i + a 2 + b 2 ,
and its discriminant is
N i = tr 2 ( J i ) 4 det ( J i ) = d 1 k 2 + d 2 k 2 Λ i + a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 4 d 1 k 2 d 2 k 2 Λ i 2 + d 1 k 2 ( a 2 + b 2 ) d 2 k 2 a 2 b 2 a 2 + b 2 Λ i + a 2 + b 2 , = d 1 k 2 d 2 k 2 2 Λ i 2 2 d 1 k 2 d 2 k 2 a 2 b 2 a 2 + b 2 + ( a 2 + b 2 ) Λ i + Δ Λ .
The stability of ( z , w ) depends on the sign of Δ 1 i . The discriminant of Δ 1 i in relation to Λ i is
N Λ i = 4 d 1 k 2 d 2 k 2 2 a 2 b 2 a 2 + b 2 + ( a 2 + b 2 ) 2 4 d 1 k 2 d 2 k 2 2 a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 4 a 2 + b 2 = 32 a d 1 k 2 d 2 k 2 2 .
Clearly, N Λ i > 0 , and because d 1 d 2 , we distinguish two cases,
  • If d 1 < d 2 , then a 2 b 2 a 2 + b 2 ( a 2 + b 2 ) 2 > 4 a 2 + b 2 , and the two solutions of the equation N Λ i = 0 are both negative. Thus, N Λ i > 0 , and the roots of (31) are
    μ 1 ( Λ i ) = tr ( J i ) + tr ( J i ) 2 4 det ( J i ) 2 , μ 2 ( Λ i ) = tr ( J i ) tr ( J i ) 2 4 det ( J i ) 2 .
    It should be noted that the solutions are real, and μ ( Λ i ) 1 < 0 . In addition, if d 1 k 2 Λ 1 2 a 2 a 2 + b 2 1 , then μ ( Λ i ) 2 < 0 . This leads to
    | A r g ( μ 1 ( Λ i ) ) | = | A r g ( μ 2 ( Λ i ) 2 ) | = π ,
    which ensures ( z , w ) is asymptotic stability.
  • If d 1 > d 2 , we have d 1 k 2 Λ 1 2 a 2 a 2 + b 2 1 . This returns us to the previous scenario Again, for d 1 k 2 Λ 1 3 a 5 2 5 1 + a 5 2 , det ( J i ) > 0 ; hence, λ 1 and λ 2 are negative and must meet the conditions of Theorem 2.

5. Global Stability

In this section, we show how the constant steady-state solution is globally asymptotically stable.
Theorem 6.
If
b w + 2 w z 2 0 ,
then system (11) is globally asymptotically stable.
Proof. 
To demonstrate this result, we employ the identical Lyapunov function as in [28]. We consider the following function for this purpose:
l ( z ) = z 1 ln ( z ) ,
where this function has a strict global minimum, i.e., l ( 1 ) = 0 . Let the following be a nonnegative function:
L ( t ) = L 1 ( t ) + L 2 ( t ) ,
where
L 1 ( t ) = i = 1 m z l z i z , L 2 ( t ) = i = 1 m w l w i w .
First, we calculate L 1 ( t )
C Δ t 0 ϑ L 1 ( t ) = i = 1 m C Δ t 0 γ z l z i ( t ) z , i = 1 m 1 z z i ( t + ϑ ) C Δ t 0 γ z i ( t ) , i = 1 m 1 z z i ( t + ϑ ) ( d 1 k 2 Δ 2 z i 1 ( t + ϑ ) + b w i ( t + ϑ ) z i ( t + ϑ ) + z i 2 ( t + ϑ ) w i ( t + ϑ ) ) , i = 1 m 1 z z i ( t + ϑ ) ( b w i ( t + ϑ ) b w z i ( t + ϑ ) + z + z i 2 ( t + ϑ ) w i ( t + ϑ ) z 2 w ) + d 1 k 2 1 z z i ( t + ϑ ) ( z i + 1 ( t + ϑ ) 2 z i ( t + ϑ ) + z i 1 ( t + ϑ ) ) .
Using (6), we obtain
C Δ t 0 ϑ L 1 ( t ) i = 1 m b w 1 z z i ( t + ϑ ) 1 w i ( t + ϑ ) w z 2 w 1 z z i ( t + ϑ ) 1 z i 2 ( t + ϑ ) w i ( t + ϑ ) z 2 w + z 1 z z i ( t + ϑ ) 1 z i ( t + ϑ ) z + i = 1 m d 1 k 2 ( z i + 1 ( t + γ ) 2 z i ( t + γ ) + z i 1 ( t + γ ) ) d 1 k 2 z z i + 1 ( t + ϑ ) z i ( t + ϑ ) + z i 1 ( t + ϑ ) z i ( t + ϑ ) 2 , i = 1 m b w l z z i ( t + ϑ ) + l w i ( t + ϑ ) w l w i ( t + ϑ ) z w z i ( t + ϑ ) + z 2 w l z z i ( t + ϑ ) l w i ( t + ϑ ) z i ( t + ϑ ) z w + l z 2 ( t + ϑ ) w ( t + ϑ ) z 2 w z l z z i ( t + ϑ ) + l z i ( t + ϑ ) z + d 1 k 2 z m + 1 ( t + ϑ ) 2 z m ( t + ϑ ) + z 0 ( t + h γ ) z 1 ( t + ϑ ) z d 1 k 2 z i + 1 ( t + ϑ ) z i ( t + ϑ ) 2 + z i 1 ( t + ϑ ) z i ( t + ϑ ) , i = 1 m ( b w + z 2 w z ) l z z i ( t + ϑ ) b w l z i ( t + ϑ ) z z l w i ( t + ϑ ) w + z 2 w l z 2 ( t + ϑ ) w ( t + ϑ ) z 2 v z 2 w l w i ( t + ϑ ) z i ( t + ϑ ) z w + z l w i ( t + ϑ ) z w z i ( t + ϑ ) z d 1 k 2 z i + 1 ( t + ϑ ) z i ( t + ϑ ) 2 + z i 1 ( t + ϑ ) z i ( t + ϑ ) .
Now, we have
C Δ t 0 ϑ L 2 ( t ) = i = 1 m C Δ t 0 ϑ w l w i ( t ) w , i = 1 m 1 w w i ( t + ϑ ) C Δ t 0 ϑ w i ( t ) , i = 1 m 1 w w i ( t + ϑ ) d 2 k 2 Δ 2 w i 1 ( t + ϑ ) + a b w i ( t + ϑ ) z i 2 ( t + ϑ ) w i ( t + ϑ ) , i = 1 m 1 w w i ( t + ϑ ) b w + z 2 w b w i ( t + ϑ ) z i 2 ( t + ϑ ) w i ( t + ϑ ) , + d 2 k 2 i = 1 m 1 w w i ( t + ϑ ) ( z i + 1 ( t ) 2 z i ( t ) + z i 1 ( t ) ) + o ( 1 ) , i = 1 m b w 1 w w i ( t + ϑ ) 1 w i ( t + ϑ ) w + w z 2 1 w w i ( t + h γ ) 1 w i ( t + ϑ ) z i 2 ( t + ϑ ) w u 2 + i = 1 m d 2 k 2 w i + 1 ( t + ϑ ) 2 w i ( t + ϑ ) + w i 1 ( t + ϑ ) d 2 k 2 w w i + 1 ( t + ϑ ) w i ( t + ϑ ) + w i 1 ( t + ϑ ) w i ( t + ϑ ) 2 , i = 1 m b w l w w i ( t + ϑ ) + l w i ( t + ϑ ) w + w z 2 l w w i ( t + ϑ ) + l z 2 z i 2 ( t + ϑ ) l w i ( t + ϑ ) z i 2 ( t + ϑ ) w z 2 + d 2 k 2 w m + 1 ( t + ϑ ) 2 w m ( t + ϑ ) + w 0 ( t + ϑ ) w 1 ( t + ϑ ) w d 2 k 2 w i + 1 ( t + ϑ ) w i ( t + ϑ ) 2 + w i 1 ( t + ϑ ) w i ( t + ϑ ) , i = 1 m ( b w w z 2 ) l w w i ( t + ϑ ) b w l w i ( t + ϑ ) w + w z 2 l z i 2 ( t + ϑ ) z 2 w z 2 l w i ( t + ϑ ) z i 2 ( t + ϑ ) w z 2 w d 2 k 2 w i + 1 ( t + ϑ ) w i ( t + ϑ ) 2 + w i 1 ( t + ϑ ) w i ( t + ϑ ) .
According to [28], we have
z i + 1 ( t + ϑ ) z i ( t + ϑ ) 2 + z i 1 ( t + ϑ ) z i ( t + ϑ ) 0 , w i + 1 ( t + ϑ ) w i ( t + ϑ ) 2 + w i 1 ( t + ϑ ) w i ( t + ϑ ) 0 .
We conclude that
C Δ t 0 ϑ L ( t ) = C Δ t 0 ϑ L 1 ( t ) + C Δ t 0 ϑ L 2 ( t ) , i = 1 m ( b w + z 2 w z ) l z z i ( t + ϑ ) + ( b w + 2 w z 2 ) l z i ( t + ϑ ) z ( b w + z ) l w i ( t + ϑ ) w ( b w + w z 2 ) l w w i ( t + ϑ ) z 2 w l w i ( t + ϑ ) z i ( t + ϑ ) z w z d 1 k 2 z i + 1 ( t + ϑ ) z i ( t + ϑ ) 2 + z i 1 ( t + ϑ ) z i ( t + ϑ ) w d 2 k 2 w i + 1 ( t + ϑ ) w i ( t + ϑ ) 2 + w i 1 ( t + ϑ ) w i ( t + ϑ ) .
Hence, since b w + 2 w z 2 0 , then,
C Δ t 0 ϑ L ( t ) 0
and C Δ t 0 ϑ L ( t ) = 0 if and only if ( z i , w i ) = ( 0 , 0 ) . ( z , w ) is therefore globally asymptotically stable under condition (36) according to Theorem 3. □

6. Simulations

Here, we demonstrate several simulations that illustrate certain theoretical aspects of the stability of the discrete fractional glycolysis reaction–diffusion system. Matlab was used to carry out all of the simulations in this section.
By changing the characteristics and order of the system, we are able to observe how it behaves. When we use the numerical solution shown below, system (11) looks like:
z i ( n ) = ϕ 1 ( z i ) + ϑ Γ ( ϑ ) r = 1 n Γ ( n r + ϑ ) Γ ( n r + 1 ) [ d 1 z i + 1 ( ( r 1 ) ) 2 z i ( ( r 1 ) ) + z i 1 ( ( r 1 ) ) k 2 + b z i ( ( r 1 ) ) + z i 2 ( ( r 1 ) ) w i ( ( r 1 ) ) ] , w i ( n ) = ϕ 2 ( z i ) + ϑ Γ ( ϑ ) r = 1 n Γ ( n r + ϑ ) Γ ( n r + 1 ) [ d 2 w i + 1 ( ( r 1 ) ) 2 w i ( ( r 1 ) ) + w i 1 ( ( r 1 ) ) k 2 + a b z i ( ( r 1 ) ) w i ( ( r 1 ) ) z i 2 ( ( r 1 ) ) w i ( ( r 1 ) ) ] , 1 i m , n > 0 .
Example 1.
To demonstrate our point, we consider the following parameter values: ( d 1 , d 2 , a , b ) = ( 1 , 2 , 0.1 , 0.01 ) , N = 100 , = 0.1 , t [ 0 , 10 ] , and z [ 0 , 20 ] and the boundary conditions ( z 0 ( t ) , w 0 ( t ) ) = ( 3 , 1 ) , ( z 1 ( t ) , w 1 ( t ) ) = ( 3 , 1 ) . Also, we provide the following appropriate initial conditions:
ϕ 1 ( z i ) = 3.5 + cos ( π z i ) 5 , ϕ 2 ( z i ) = 2 + cos ( π z i ) 5 .
We can see that all of the solutions to our model ultimately reach the same unique positive equilibrium ( z , w ) = ( 0.1 , 4 ) . The unique equilibrium is hence asymptotically stable. For a, which denotes the dimensionless input flow, and b, which denotes the dimensionless constant rate for the low activity state, the conditions in Theorem 6 are clearly demonstrated, and this numerical finding is consistent with our earlier theoretical conclusions.
The results stated previously for various fractional orders are supported by Figure 1, Figure 2 and Figure 3, which describe the dynamic behavior of the system (11) and show that as the fractional order ϑ approaches 1, the model’s solution becomes more stable as it rapidly converges to the unique equilibrium.

7. Conclusions

In this study, a discrete equivalent of the fractional reaction–diffusion glycolysis model was investigated using the second-order central difference and L1 finite difference schemes. This unique approach’s capacity to keep the system confined and positive should be emphasized. Moreover, the local stability of the discrete system was explored. The results were consistent with those of the continuous system that serves as a continuous analog. This suggests that the discrete system generated by the considered schemes can preserve as continuous system characteristics like the positivity, boundedness, and global stability of the equilibrium point, which are first evaluated as a way to assess the stability of discrete fractional reaction diffusion systems using Lyapunov functions. Additionally, computations and simulations show the effectiveness and viability of the suggested system. In addition, the Lyapunov functional and the linearization approach may be utilized to solve the stability issue in discrete fractional reaction–diffusion models. More importantly, the results of this study might be easily extended to a variety of reaction–diffusion discrete fractional spatiotemporal systems, as well as other dynamical problems including chaos and synchronization control.

Author Contributions

Formal analysis, T.H. and O.A.; Investigation, A.H.; Software, A.H.; Super-vision, A.O.; Validation, Y.A.A.-K.; Visualization, A.A.-H.; Writing—original draft, A.H.; Writing—review and editing, A.O. All authors haveread and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mojikon, F.D.; Kasimin, M.E.; Molujin, A.M.; Gansau, J.A.; Jawan, R. Probiotication of Nutritious Fruit and Vegetable Juices: An Alternative to Dairy-Based Probiotic Functional Products. Nutrients 2022, 14, 3457. [Google Scholar] [CrossRef] [PubMed]
  2. Romano, A.H.; Conway, T. Evolution of carbohydrate metabolic pathways. Res. Microbiol. 1996, 147, 448. [Google Scholar] [CrossRef] [PubMed]
  3. Ouannas, A.; Batiha, I.M.; Bekiros, S.; Liu, J.; Jahanshahi, H.; Aly, A.A.; Alghtani, A.H. Synchronization of the glycolysis reaction-diffusion model via linear control law. Entropy 2021, 23, 1516. [Google Scholar] [CrossRef] [PubMed]
  4. Ahmed, N.; Ss, T.; Imran, M.; Rafiq, M.; Rehman, M.A.; Younis, M. Numerical analysis of auto-catalytic glycolysis model. AIP Adv. 2019, 9, 085213. [Google Scholar] [CrossRef] [Green Version]
  5. Al Noufaey, K.S. Stability analysis for Selkov-Schnakenberg reaction-diffusion system. Open Math. 2021, 19, 46–62. [Google Scholar] [CrossRef]
  6. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: Hoboken, NJ, USA, 1993. [Google Scholar]
  7. Oldham, K.; Spanier, J. The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order; Elsevier: Amsterdam, The Netherlands, 1974. [Google Scholar]
  8. Dababneh, A.; Djenina, N.; Ouannas, A.; Grassi, G.; Batiha, I.M.; Jebril, I.H. A new incommensurate fractional-order discrete COVID-19 model with vaccinated individuals compartment. Fractal Fract. 2022, 6, 456. [Google Scholar] [CrossRef]
  9. Saadeh, R.; Abbes, A.; Al-Husban, A.; Ouannas, A.; Grassi, G. The Fractional Discrete Predator–Prey Model: Chaos, Control and Synchronization. Fractal Fract. 2023, 7, 120. [Google Scholar] [CrossRef]
  10. Hilfer, R. (Ed.) Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  11. Silva, M.F.; Machado, J.T. Fractional order PDϑ joint control of legged robots. J. Vib. Control 2006, 12, 1483–1501. [Google Scholar] [CrossRef]
  12. Duarte, F.B.; Machado, J.T. Chaotic phenomena and fractional-order dynamics in the trajectory control of redundant manipulators. Nonlinear Dyn. 2002, 29, 315–342. [Google Scholar] [CrossRef]
  13. Tenreiro Machado, J.A.; Silva, M.F.; Barbosa, R.S.; Jesus, I.S.; Reis, C.M.; Marcos, M.G.; Galhano, A.F. Some applications of fractional calculus in engineering. Math. Probl. Eng. 2010, 2010, 639801. [Google Scholar] [CrossRef] [Green Version]
  14. Gafiychuk, V.; Datsko, B.; Meleshko, V. Mathematical modeling of time fractional reaction–diffusion systems. J. Comput. Appl. Math. 2008, 220, 215–225. [Google Scholar] [CrossRef] [Green Version]
  15. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Further solutions of fractional reaction–diffusion equations in terms of the H-function. J. Comput. Appl. Math. 2011, 235, 1311–1316. [Google Scholar] [CrossRef]
  16. Che, H.; Wang, Y.L.; Li, Z.Y. Novel patterns in a class of fractional reaction–diffusion models with the Riesz fractional derivative. Math. Comput. Simul. 2022, 202, 149–163. [Google Scholar] [CrossRef]
  17. Belmahi, N.; Shawagfeh, N. A new mathematical model for the glycolysis phenomenon involving Caputo fractional derivative: Well posedness, stability and bifurcation. Solitons Fractals 2021, 142, 110520. [Google Scholar] [CrossRef]
  18. Abdel-Aty, A.H.; Khater, M.M.; Baleanu, D.; Abo-Dahab, S.M.; Bouslimi, J.; Omri, M. Oblique explicit wave solutions of the fractional biological population (BP) and equal width (EW) models. Adv. Differ. Equ. 2020, 2020, 552. [Google Scholar] [CrossRef]
  19. Angstmann, C.N.; Henry, B.I. Time Fractional Fisher–KPP and Fitzhugh–Nagumo Equations. Entropy 2020, 22, 1035. [Google Scholar] [CrossRef]
  20. Abad, E.; Yuste, S.B.; Lindenberg, K. Reaction-subdiffusion and reaction-superdiffusion equations for evanescent particles performing continuous-time random walks. Phys. Rev. E 2010, 81, 031115. [Google Scholar] [CrossRef] [Green Version]
  21. Angstmann, C.N.; Donnelly, I.C.; Henry, B.I. Continuous time random walks with reactions forcing and trapping. Math. Model. Nat. Phenomena 2013, 8, 17–27. [Google Scholar] [CrossRef] [Green Version]
  22. Klein, C.T.; Mayer, B. A model for pattern formation in gap-junction coupled cells. J. Theor. Biol. 1997, 186, 107–115. [Google Scholar] [CrossRef]
  23. Maselko, J. Mosaic Pattern Formations in Multicellular Chemical Systems. J. Phys. Chem. 1995, 99, 2949–2952. [Google Scholar] [CrossRef]
  24. Anakira, N.; Hioual, A.; Ouannas, A.; Oussaeif, T.E.; Batiha, I.M. Global Asymptotic Stability for Discrete-Time SEI Reaction-Diffusion Model. In Proceedings of the International Conference on Mathematics and Computations, Zarqa, Jordan, 11–13 May 2022; Springer Nature: Singapore, 2022; pp. 345–357. [Google Scholar]
  25. Almatroud, O.A.; Hioual, A.; Ouannas, A.; Grassi, G. On Fractional-Order Discrete-Time Reaction Diffusion Systems. Mathematics 2023, 11, 2447. [Google Scholar] [CrossRef]
  26. Wu, G.C.; Baleanu, D.; Zeng, S.D.; Deng, Z.G. Discrete fractional diffusion equation. Nonlinear Dyn. 2015, 80, 281. [Google Scholar] [CrossRef]
  27. Mesdoui, F.; Ouannas, A.; Shawagfeh, N.; Grassi, G.; Pham, V.T. Synchronization methods for the Degn-Harrison reaction-diffusion systems. IEEE Access 2020, 8, 91829–91836. [Google Scholar] [CrossRef]
  28. Zhenzhen, L.; Yongguang, Y.; Weiyi, X.; Guojian, R.; Xiangyun, M. Global dynamics for a class of discrete fractional epidemic model with reaction-diffusion. arXiv 2022, arXiv:2208.06548. [Google Scholar]
  29. Bhargava, S.C. On the higgins model of glycolysis. Bull. Math. Biol. 1980, 42, 829–836. [Google Scholar] [CrossRef] [PubMed]
  30. Peng, R. Qualitative analysis of steady states to the Sel’kov model. J. Differ. Equ. 2007, 241, 386–398. [Google Scholar] [CrossRef] [Green Version]
  31. Lee, K.J.; Swinney, H.L. Replicating spots in reaction-diffusion systems. Int. J. Bifurc. Chaos 1997, 7, 1149–1158. [Google Scholar] [CrossRef]
  32. You, Y. Asymptotical dynamics of Selkov equations. Discret. Contin. Dyn. Syst.-S 2008, 2, 193–219. [Google Scholar] [CrossRef]
  33. Kelley, W.G.; Peterson, A.C. Difference Equations: An Introduction with Applications; Academic Press: Cambridge, MA, USA, 2001. [Google Scholar]
  34. Baleanu, D.; Wu, G.C.; Bai, Y.R.; Chen, F.L. Stability analysis of Caputo–like discrete fractional systems. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 520–530. [Google Scholar] [CrossRef]
  35. Hioual, A.; Ouannas, A.; Oussaeif, T.E.; Grassi, G.; Batiha, I.M.; Momani, S. On Variable-Order Fractional Discrete Neural Networks: Solvability and Stability. Fractal Fract. 2022, 6, 119. [Google Scholar] [CrossRef]
  36. Abdeljawad, T. On Riemann and Caputo fractional differences. Comput. Math. Appl. 2011, 62, 1602–1611. [Google Scholar] [CrossRef] [Green Version]
  37. Abdeljawad, T. Different type kernel h—Fractional differences and their fractional h—Sums. Solitons Fractals 2018, 116, 146–156. [Google Scholar] [CrossRef]
  38. Cermák, J.; Nechvátal, L. On a problem of linearized stability for fractional difference equations. Nonlinear Dyn. 2021, 104, 1253–1267. [Google Scholar] [CrossRef]
  39. Ashkenazi, M.; Othmer, H.G. Spatial patterns in coupled biochemical oscillators. J. Math. Biol. 1977, 5, 305–350. [Google Scholar] [CrossRef]
  40. Tyson, J.; Kauffman, S. Control of mitosis by a continuous biochemical oscillation: Synchronization; spatially inhomogeneous oscillations. J. Math. Biol. 1975, 1, 289–310. [Google Scholar] [CrossRef]
  41. Davidson, F.A.; Rynne, B.P. A priori bounds and global existence of solutions of the steady-state Sel’kov model. Proc. R. Soc. Edinb. Sect. A Math. 2000, 130, 507–516. [Google Scholar] [CrossRef]
  42. Furter, J.E.; Eilbeck, J.C. Analysis of bifurcations in reaction–diffusion systems with no-flux boundary conditions: The Sel’kov model. Proc. R. Soc. Edinb. Sect. A Math. 1995, 125, 413–438. [Google Scholar] [CrossRef]
  43. Peng, R.; Wang, M.; Yang, M. Positive steady-state solutions of the Sel’kov model. Math. Comput. Model. 2006, 44, 945–951. [Google Scholar] [CrossRef]
  44. Xu, L.; Zhao, L.J.; Chang, Z.X.; Feng, J.T.; Zhang, G. Turing instability and pattern formation in a semi-discrete Brusselator model. Mod. Phys. Lett. B 2013, 27, 1350006. [Google Scholar] [CrossRef]
Figure 1. The state trajectories of r z i ( t ) and w i ( t ) for ϑ = 0.1 .
Figure 1. The state trajectories of r z i ( t ) and w i ( t ) for ϑ = 0.1 .
Fractalfract 07 00587 g001
Figure 2. The dynamic behaviors of z i ( t ) and w i ( t ) for N = 10 , ( a , b , d 1 , d 2 ) = ( 1 , 2 , 0.1 , 0.01 ) and ϑ = 0.1 .
Figure 2. The dynamic behaviors of z i ( t ) and w i ( t ) for N = 10 , ( a , b , d 1 , d 2 ) = ( 1 , 2 , 0.1 , 0.01 ) and ϑ = 0.1 .
Fractalfract 07 00587 g002
Figure 3. The dynamic behaviors of z i ( t ) and w i ( t ) for N = 10 , ( a , b , d 1 , d 2 ) = ( 1 , 2 , 0.1 , 0.01 ) and ϑ = 0.9 .
Figure 3. The dynamic behaviors of z i ( t ) and w i ( t ) for N = 10 , ( a , b , d 1 , d 2 ) = ( 1 , 2 , 0.1 , 0.01 ) and ϑ = 0.9 .
Fractalfract 07 00587 g003
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

Hamadneh, T.; Hioual, A.; Alsayyed, O.; AL-Khassawneh, Y.A.; Al-Husban, A.; Ouannas, A. Local Stability, Global Stability, and Simulations in a Fractional Discrete Glycolysis Reaction–Diffusion Model. Fractal Fract. 2023, 7, 587. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7080587

AMA Style

Hamadneh T, Hioual A, Alsayyed O, AL-Khassawneh YA, Al-Husban A, Ouannas A. Local Stability, Global Stability, and Simulations in a Fractional Discrete Glycolysis Reaction–Diffusion Model. Fractal and Fractional. 2023; 7(8):587. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7080587

Chicago/Turabian Style

Hamadneh, Tareq, Amel Hioual, Omar Alsayyed, Yazan Alaya AL-Khassawneh, Abdallah Al-Husban, and Adel Ouannas. 2023. "Local Stability, Global Stability, and Simulations in a Fractional Discrete Glycolysis Reaction–Diffusion Model" Fractal and Fractional 7, no. 8: 587. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7080587

Article Metrics

Back to TopTop