Global dynamics for a class of discrete fractional epidemic model with reaction-diffusion

Lu Zhenzhen,Yu Yongguang,Xu Weiyi,Ren Guojian,Meng Xiangyun

DOI: https://doi.org/10.48550/arXiv.2208.06548

2022-08-13

Abstract:In recent years, discrete fractional epidemic models with reaction-diffusion have become increasingly popular in the literature, not only for its necessity of numerical simulation, but also for its defined physical processes. In this paper, by second order central difference scheme and L1 nonstandard finite difference scheme, a discrete counterpart of time-fractional reaction-diffusion epidemic model with generalized incidence rate is considered. More importantly, the main idea in choosing an nonstandard finite difference scheme is to obtain unconditionally positivity in the proposed system, which leads to the proposal of the discrete epidemic model with time delay. Furthermore, the global properties of the proposed discrete system are studied, including the global boundedness of positive solutions, the existence and the global stability of equilibrium points, which are consistent with the corresponding continuous systems. Meanwhile, it shows that L1 nonstandard finite difference scheme and second order central difference scheme can keep the properties of the corresponding continuous system well. It is worth noting that, different from discrete epidemic model with integer-order, the memory Lyapunov function is constructed in this paper, which depends on the previous historical information of the proposed system. This is consistent with the non-local property of Caputo fractional derivatives. Finally, numerical results are given to verify the theoretical results.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to construct and analyze a discrete fractional - order epidemic model with reaction - diffusion to better understand and predict the spread dynamics of infectious diseases. Specifically, the authors focus on the dynamic behavior of the epidemic model under the influence of time - fractional - order derivatives and spatial diffusion.

### Main Problem Analysis

1. **Model Selection and Construction**:
- Traditional integer - order differential equation models cannot fully describe the spread characteristics of some infectious diseases, especially when these diseases depend on the current state and historical state.
- Therefore, the authors introduce the Caputo fractional - order derivative to consider historical dependence and combine reaction - diffusion terms to simulate the spatial spread effect.

2. **Choice of Numerical Method**:
- To discretize the time - fractional - order derivative, the authors adopt the L1 non - standard finite difference scheme to ensure the unconditional positivity and stability of the system.
- In the spatial direction, the second - order central difference scheme is used to handle the diffusion term.

3. **Study of System Dynamic Properties**:
- The authors study the global properties of the discrete system, including the global boundedness of positive solutions, the existence and global stability of equilibrium points.
- A memory Lyapunov function depending on the system's historical information is constructed to ensure the system's stability.

4. **α - Robustness**:
- By constructing a Lyapunov function, the authors prove the α - robustness of the L1 scheme when the fractional - order approaches the integer - order, that is, as the fractional - order parameter α approaches 1, the stability of the discrete system is consistent with that of the integer - order system.

### Formula Representation

- The Caputo fractional - order derivative is defined as:
\[
{}_0^C D_t^\alpha f(t)=\frac{1}{\Gamma(n - \alpha)}\int_0^t\frac{f^{(n)}(s)}{(t - s)^{\alpha - n + 1}}ds
\]
where \(n - 1<\alpha <n\), and \(\Gamma(\cdot)\) is the Gamma function.

- The L1 discrete scheme is defined as:
\[
{}_0^C D_t^\alpha f(t_n)\approx(\Delta t)^{-\alpha}\frac{1}{\Gamma(2 - \alpha)}\left[f(t_n)-b_{n - 1}^{(1 - \alpha)}f(t_0)-\sum_{j = 1}^{n - 1}(b_{j - 1}^{(1 - \alpha)}-b_j^{(1 - \alpha)})f(t_{n - j})\right]
\]
where \(b_j^{(1 - \alpha)}=(j + 1)^{1 - \alpha}-j^{1 - \alpha}\).

### Summary

The main contribution of this paper is to propose a new discrete fractional - order reaction - diffusion epidemic model and prove the global stability and α - robustness of this model through strict mathematical analysis. This provides a powerful theoretical tool for understanding the spatio - temporal spread of complex infectious diseases.