Modelling and Simulation of Complex Biological SystemsView this Special Issue
Global Stability for Fractional Diffusion Equations in Biological Systems
This paper proposes a new method of construction of Lyapunov functionals for the dynamical systems described by fractional differential equations and fractional partial differential equations. The proposed method is rigorously presented. Furthermore, the method is applied to establish the global stability of some fractional biological models with and without diffusion.
In recent years, fractional differential equations (FDEs) are used to describe the temporal dynamics of various systems in many fields. These equations are the generalization of the classical ordinary differential equations (ODEs). However, fractional partial differential equations (FPDEs) are the generalization of the partial differential equations (PDEs) which can be an effective tool to describe the spatiotemporal dynamics of several phenomena with memory or have hereditary properties.
The construction of Lyapunov functionals to prove the global stability of fractional dynamic systems has attracted the attention of some authors. Aguila-Camacho et al.  established a new lemma for fractional derivative in Caputo sense with order . They used this lemma to demonstrate the stability of some fractional-order systems by mean of quadratic Lyapunov functionals. Duarte-Mermoud et al.  extended the lemma of  to a vector of differentiable functions, and they used Lyapunov functionals containing general quadratic forms in order to analyze the stability of fractional-order model reference adaptive control (FOMRAC) schemes. Vargas-De-León  extended the Volterra-type Lyapunov function to fractional-order epidemic systems via an inequality to estimate the Caputo fractional derivative of this function. On the other hand, a study in  has been devoted to establish the global stability for some diffusion equations in biology by means of Lyapunov functionals.
The methods mentioned above are applied for particular Lyapunov functionals such as quadratic or Volterra-type. Likewise, the work in  is especially applicable for the models formulated by PDEs. Therefore, the main goal of this study is to develop a new mathematical method to construct the Lyapunov functionals for FDEs and FPDEs based on those of ODEs. To do this, the next section deals with the description of the method and the last section is devoted to the application of our method to investigate the global stability of some mathematical models in epidemiology as well as in virology.
2. Description of the Method
Consider the following FDE:where is the fractional derivative in the Caputo sense of order , the state variable is a non-negative vector of concentrations , …, , and : is a function. It is obvious that if , then (1) becomes the following ordinary differential equation:
Let be a bounded domain in with smooth boundary and with . Assume that is a steady state of (1). Then, is also the steady state of the following fractional diffusion system with homogeneous Neumann boundary condition:where represents the Laplacian operator and denotes the outward normal derivative on the boundary .
Let be a function defined on some domain in and is a solution of (1). Further, we suppose that the range of is contained in the domain of andwhose equality holds if .
We observe that the right-hand side of the above inequality is given by the scalar product of the gradient of the function and the vector field . Hence, the right-hand side is defined without the fact that is a solution of (1), which is very important for the construction of Lyapunov functionals.
Let be a solution of (3). Denote
The fractional time derivative of along the positive solution of (3) satisfies
By applying Green’s formula, we find
According to on , we have
Additionally, we assume that the function satisfies the following condition:
From the above, it is not difficult to obtain the following result.
Theorem 1. Let be a Lyapunov functional for the ordinary differential equation (2).(i)If satisfies the condition (4), then is also a Lyapunov functional for fractional differential equation (1).(ii)If satisfies the conditions (4) and (10), then the function defined by (5) is a Lyapunov functional for fractional diffusion system (3).In the literature, several researchers constructed the Lyapunov functional in the following form:
Corollary 2. If is a Lyapunov functional for ordinary differential equation (2) of the form given in (11), then is a Lyapunov functional for fractional differential equation (1). Moreover, the function defined by (5) is a Lyapunov functional for fractional diffusion system (3).
Proof. We haveBy applying Lemma 3.1 in , we getThen, satisfies the condition (4). It follows from Theorem 1 (i) that is also a Lyapunov functional for fractional differential equation (1).
On the other hand, we havewhich implies that satisfies the condition (10). According to Theorem 1 (ii), we deduce that the function given by (4) is also a Lyapunov functional for fractional diffusion system (3). This completes the proof.
Remark 3. The method described above can be used to prove the stability of many fractional systems with and without diffusion. It is very important to recall that the steady state is stable if there exists a Lyapunov functional satisfying . Moreover, if for all , then is asymptotically stable. Additionally, according to , if and the largest invariant set in is the singleton , then is asymptotically stable. This means that the solution of the system starting from any initial conditions converges to .
This section focuses on the application of the method described in the above section in order to establish the global stability of some fractional diffusion biological models by constructing Lyapunov functionals from those of the corresponding systems which are formulated by ODEs.Example 1. Consider the SIR epidemic model described by the following nonlinear system of FDEs:where , , and are the populations of susceptible, infected, and recovered individuals, respectively. The parameters , , , and are, respectively, the recruitment rate, the natural death rate, the death rate due to disease, and the recovery rate. The incidence function of system (15) is described by Hattaf–Yousfi functional response  of the form , where the non-negative constants , , measure the saturation, inhibitory, or psychological effects, and the positive constant is the infection rate. This functional response covers the most famous forms existing in the literature such as the classical bilinear incidence, the saturated incidence, the Beddington–DeAngelis functional response , the Crowley–Martin functional response , and the specific functional response introduced in . Further, the fractional models proposed in [10, 11] are particular cases of model (15); it suffices to take and for  and for .
Due to the great mobility of individuals inside or outside a country or region, we consider the following fractional model:where and . The parameters and are the diffusion coefficients for the susceptible and infected populations, respectively. Also, we consider model (17) with homogeneous Neumann boundary conditions:and initial conditions:
For , model (16) becomes the following nonlinear system of ODEs:
Obviously, model (20) has always one disease-free equilibrium . By a simple computation, the basic reproduction number is given by
When , model (20) has another equilibrium named endemic equilibrium , wherewith .
System (20) is a special case of the mathematical model presented in . Thus, the disease-free equilibrium is globally asymptotically stable when . However, becomes unstable and the endemic equilibrium is globally asymptotically stable if .
Additionally, we have
Now, we construct the Lyapunov functional for fractional diffusion model (17) at as follows:
In this case, we have
For the global stability of the endemic equilibrium , we consider the following function:which is a Lyapunov functional for ODE model (20) at . On the other hand, we have
It is not hard to show that
Consequently, we have the following:(i)When , the disease-free equilibrium of (17)–(19) is globally asymptotically stable.(ii)When , the endemic equilibrium of (17)–(19) is globally asymptotically stable.Example 2. Consider the following fractional viral infection model:where healthy cells () are produced at rate , die at rate , and become infected by contact with virus at rate and by contact with infected cells at rate . Infected cells () die at rate . Free viral particles () are released from infected cells at rate and decay at rate . Note that when and , we obtain the basic model which was used to study some viral infections such as HIV, HBV, and HCV [13–16].
Now, we extend the above fractional model (33) by taking into account the mobility of the virus as well as the cells. Then, model (33) becomeswhere , , and denote the concentrations of healthy cells, infected cells, and free viral particles at position and at time , respectively. The parameters , , and represent, respectively, the diffusion coefficients of healthy cells, infected cells, and free viral particles.
The cells and virus are biological quantities which should be non-negative, and they also should not cross the boundary of the domain. Then, we consider system (34) with the homogeneous Neumann boundary conditions:and initial conditions:
The corresponding ordinary system of the above models is given by
This system admits a unique infection-free equilibrium , where . Further, the basic reproduction number is given by
This basic reproduction number is the sum of two basic reproduction numbers due to virus-to-cell and cell-to-cell transmissions. If , ODE model (37) has another biological steady state called the infection chronic which is labeled by , where
Calculating the time derivative of along the solution of (37), we get
If , then with equality if and only if , and . Therefore, the infection-free equilibrium is globally asymptotically stable when . Sinceit follows from Theorem 1 (i) that is also a Lyapunov functional for FDE model (33) at .
Next, we construct a Lyapunov functional for fractional diffusion model (34) at as follows:
Since has the form given in (11), we conclude, by applying Corollary 2, that is a Lyapunov functional for FDE model (33) and is a Lyapunov functional for fractional diffusion systems (34)–(36) at when .
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
M. A. Duarte-Mermoud, N. Aguila-Camacho, J. A. Gallegos, and R. Castro-Linares, “Using general quadratic Lyapunov functions to prove Lyapunov uniform stability for fractional order systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 22, no. 1-3, pp. 650–659, 2015.View at: Publisher Site | Google Scholar