Global attractivity of the equilibria of the diffusive SIR and SEIR epidemic models with multiple parallel infectious stages and nonlinear incidence mechanism
Introduction
In this paper, we consider the following two diffusive epidemic models with the nonlinear incidence mechanism : and where ; stand for, respectively, the densities at location and time of susceptible individuals, exposed individuals and infected individuals, and the positive constants , and are their respective diffusion coefficients. The formula suggests that the susceptible individuals carrying capability are linearly increasing, where the positive constant and account for the rates of birth and death of the susceptible individuals, respectively. The positive constant measures the infection force. The constant is the rate for an infected individual to enter the th infective compartment and is the rate of infected individuals of th compartment leaving this compartment. In addition, is the sum of the recovery rate and the death rate in system (1.1), while it is the sum of the rates of death, recovered and the amelioration of infection individuals in system (1.2) (and hence ). The positive constant is the rate that exposed hosts leave -compartment . In (1.1), (1.2), is a bounded smooth domain of , the homogeneous Neumann boundary conditions imply that there is no population flux across the boundary .
By the biological implication, it is natural to assume that the initial data are nonnegative on . To ensure the local existence and uniqueness of solutions of (1.2), we additionally assume that are Hölder continuous on , and on when . Thus, the standard existence theory of parabolic equations shows that for such initial data, (1.2) admits a classical solution . It further follows from the parabolic comparison principle that and for all and . Furthermore, if for all , then for all and . A similar assumption is imposed for the initial data of the system (1.1), and the same conclusion also holds for (1.1).
Nowadays, people have recognized that individual motility are significant factors in the study of the spread of infectious diseases, and many reaction–diffusion epidemic models have been developed to investigate the dynamics of disease transmissions. These models are derived from the ODE compartmental epidemic models by introducing random diffusion terms to describe the movement of individuals; see [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], to list a few.
The corresponding ODE systems of (1.1), (1.2) with the standard bilinear incidence mechanism were studied in [16]. Korobeinikov [16] found that the levels of contagiousness and the lengths of the infectious periods did not always the same by studying the Hepatitis B. Then they proposed SIR and SEIR systems with several parallel infectious stages (there has alternative infectious pathways and noninteracting infective subclasses ) and the bilinear incidence mechanism . This kind of incidence mechanism was first introduced by Kermack and McKendrick [17] to describe the plague epidemic in Bombay. Since then, such a bilinear infection mechanism has been used in many models, for instance, [4], [18], [19], [20], [21]. However, as mentioned by [22], [23], [24], this kind of infection mechanism has some shortcomings and may require modifications in certain situations to better model the disease transmission. Hence, many different nonlinear infection mechanisms have been proposed. Among them, the nonlinear incidence form () has been used to describe the transmission of infectious disease, for example, in [22], [24], [25], [26], [27], [28]. Especially, the nonlinear incidence mechanism () was investigated in this article.
For the corresponding ODE systems of (1.1), (1.2) with the nonlinear incidence mechanism is replaced by , Korobeinikov [16] analysed the global stability of the disease-free equilibrium and endemic equilibrium by the Lyapunov method. His results implies that the classic compartment models of infectious diseases are staunch and that the diversity of the infection level and the lengths of the infectious period do not affect the basic properties of these models. For more studies on multistage or multigroup models of infectious diseases, one may refer to [29], [30], [31], [32], [33] and the references therein.
In general, the global stability of systems (1.1), (1.2) is determined by the basic reproduction number. Next, we define the basic reproduction number [34] for systems (1.1), (1.2) as follows:
The main purpose of this paper is to discuss the global stability of the disease-free equilibrium (DFE) and the endemic equilibrium (EE) of systems (1.1), (1.2). For any , it can be easily shown that the system (1.1) has two equilibria: a DFE , and a unique EE with such that , where provided that .
For the system (1.1), we can state the following result on the global stability of the disease-free equilibrium and the endemic equilibrium.
Theorem 1.1 The following assertions hold for (1.1). If , then the disease-free equilibrium is globally attractive. If , then the unique endemic equilibrium is globally attractive.
We now turn to the system (1.2). For any , the system (1.2) has DFE . Next, we need to determine the existence of endemic equilibrium. The EE of (1.2) satisfies , , .
For any , we have . Inserting this into the second equation, we can get with and . It follows from and that is a positive constant. Adding the first and second equation, it is easily seen that (1.2) has a unique positive solution with provided that .
For the system (1.2), we have the following result on the global stability of the disease-free equilibrium and the endemic equilibrium.
Theorem 1.2 The following assertions hold for the system (1.2). If , then the disease-free equilibrium is globally attractive. If , then the endemic equilibrium is globally attractive.
To prove Theorem 1.1, Theorem 1.2, as in [16], we shall employ the Lyapunov function approach. However, it is worth mentioning that the construction of the Lyapunov function in the case of is rather different from that for the case of ; in the former case, we just need to modify the ones used in [16], but in the latter case we have to use some new types of Lyapunov function. Moreover, in order to derive the desired results, we need to establish the uniform bounds of solutions to (1.1), (1.2).
The rest of this paper is devoted to the uniform bounds of solutions to (1.1), (1.2) and the proof of the global attractivity of the disease-free equilibrium and the endemic equilibrium of (1.1), (1.2).
Section snippets
Proof of Theorems 1.1 and 1.2
In this section, we will present the proof of Theorem 1.1, Theorem 1.2.
Conclusion
In this paper, we consider the reaction–diffusion models (1.1), (1.2) with multiple parallel infectious stages and a more general nonlinear incidence function. The main goal of this paper is to study the effect of the general nonlinear incidence mechanism, the diversity infectious pathways and the population dispersal on the global dynamics of the solutions of the two models. The related problems of SIR and SEIR epidemic models with one infection stage in [37], [38], and the ODE model with
Acknowledgements
We would like to thank the handling editor and the two anonymous referees for their careful reading and valuable suggestions which improved the manuscript. We are also grateful to Prof. Rui Peng for the initial discussion.
References (38)
- et al.
Dynamics and asymptotic profiles of steady states of an epidemic model in advective environments
J. Differential Equations
(2017) - et al.
Long-time dynamics of an SIRS reaction–diffusion epidemic model
J. Math. Anal. Appl.
(2019) - et al.
Analysis on a diffusive SIS epidemic model with logistic source
Z. Angew. Math. Phys.
(2017) - et al.
Vary total population enhances disease persistence: qualitative analysis on a diffusive SIS epidemic model
J. Differential Equations
(2017) Asymptotic profiles of the positive steady state for an SIS epidemic reaction–diffusion model part I
J. Differential Equations
(2009)- et al.
Global stability of the steady states of an SIS epidemic reaction–diffusion model
Nonlinear Anal.
(2009) - et al.
Asymptotic profile of the positive steady state for an SIS epidemic reaction–diffusion model: Effects of epidemic risk and population movement
Phys. D.
(2013) - et al.
A spatial SEIRS reaction–diffusion model in heterogeneous environment
J. Differential Equations
(2019) - et al.
Asymptotic profiles of steady states for a diffusive SIS epidemic model with mass action infection mechanism
J. Differential Equations
(2016) - et al.
Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission
Math. Biosci.
(2002)
A semilinear parabolic system arising in the theory of superconductivity
J. Differential Equations
Global stability of equilibria of a diffusive SEIR epidemic model with nonlinear incidence
Appl. Math. Lett.
Asymptotic profiles of the steady states for an SIS epidemic reaction–diffusion model
Discrete Contin. Dyn. Syst.
Asymptotic profiles of the endemic equilibrium of a reaction–diffusion-advection SIS epidemic model with saturated incidence rate
Discrete Contin. Dyn. Syst. Ser. B
Dynamics of an SIS epidemic reaction–diffusion model
Proc. Roy. Soc. Edinburgh Sect. A
Concentration profile of the endemic equilibria of a reaction–diffusion-advection SIS epidemic model
Calc. Var. Partial Differential Equations
Dynamics and asymptotic profiles of endemic equilibrium for two frequency-dependent SIS epidemic models with cross-diffusion
European J. Appl. Math.
A reaction–diffusion SIS epidemic model in a time-periodic environment
Nonlinearity
Analysis on a diffusive SIS epidemic system with linear source and frequency-dependent incidence function in a heterogeneous environment
Math. Biosci. Eng.
Cited by (2)
Threshold dynamics of stochastic H7N9 model with Markov switching and hybrid strategy
2024, Journal of the Franklin InstituteFixed-Time Synergetic Control Based on SEAIQR Model for COVID-19 Epidemic
2024, IEEE Transactions on Circuits and Systems II: Express Briefs
- 1
C. Lei was partially supported by NSF of China (No. 11801232), the NSF of Jiangsu Province (No. BK20180999) and the Foundation of Jiangsu Normal University (No. 17XLR008).