key: cord-0786227-8wov472l authors: Cao, Qi; Liu, Yuying; Yang, Wensheng title: Global dynamics of a diffusive SIR epidemic model with saturated incidence rate and discontinuous treatments date: 2022-03-18 journal: Int J Dyn Control DOI: 10.1007/s40435-022-00935-3 sha: 184e5322cd50bca1e21250177bbf740ee0e4ba20 doc_id: 786227 cord_uid: 8wov472l In this paper, we study a diffusive SIR epidemic model with saturated incidence rate and discontinuous treatments under Neumann boundary conditions. Firstly, the existence and boundedness of the solution of the system are addressed. Then, on the basis of the differential inclusions theory, we analysis the existence of endemic equilibrium. Furthermore, by constructing different appropriate Lyapunov functions, we investigate the global asymptotic stability of the disease free equilibrium(DFE) and the endemic equilibrium(EE), respectively. Additionally, numerical simulations are given to confirm the correctness of theorem. Finally, we give a brief conclusion and discussion in the end of the paper. Nowadays, The SARS-Cov-2 pandemic pose a great threat to human health, which has introduced an evident research boom into biophysical and mathematical modeling of infection expansions. In order to develop control strategies to prevent disease epidemics and reduce the number of infections, many scholars have used mathematical models of reactiondiffusion equations to study the infection mechanisms of infectious diseases, for example in [1] [2] [3] [4] [5] [6] . Moreover, reaction -diffusion mechanisms have been successfully applied to many patterning phenomena in predator-prey system [7] [8] [9] [10] [11] [12] . It should be noted that most of the above models have continuous terms. But in practice, when an infectious disease occurs in the host population, some treatment measures need to be considered, which usually was described by some discontinuous (or non-smooth) control functions when constructing the mathematical model. Moreover, the discontinuous control strategy has been also extensively studied For instance, discontinuous harvesting on fishery was investigated in [13] . Li et al. [14] studied the the dynamic behaviour of a computer worm virus system with discontinuous control strategy. In [15] , Guo et al. considered the impact of discontinuous treatments to SIR epidemic system. Zhang and Zhao considered a predator-prey model with discontinuous harvesting policy with spatial diffusion in [16] . However, there are few results in the public literature on the effects of discontinuous treatments in SIR epidemic diffusive system. Recently, in [17] , Li et al. studied the global dynamics of a diffusive SIR epidemic system with linear incidence rate and discontinuous term. Other relate work can be found in [18] [19] [20] . Inspired by [17] and the above discussions. In this article, we consider a diffusive SIR epidemic model with nonlinear saturated incidence rate and discontinuous treatments. First, we denote the densities of susceptible and infected individuals at position x and time t by S(x, t) and I (x, t), respectively. Besides, a saturated incidence rate function g(I ) = β I 1+m I is considered, which was first proposed by Capasso and Serio in [21] , subsequently, many scholars have done relevant research about this type of incidence rate (for example, [22] [23] [24] ). Where m > 0 is the saturation coefficient, β > 0 is the rate of disease transmission. This incidence rate seems to be more realistic in some cases, as the number of effective contacts between infected and susceptible individuals may be saturated at high levels of infection due to congestion in infected individuals or protective measures in susceptible individuals. Therefore, we focus on the following epidemic reaction-diffusion model under Neumann boundary conditions where S(x, t), I (x, t) and m are described as before. is a bounded open set in R n , n is the outward unit normal vector of the boundary ∂ . The homogeneous Neumann boundary conditions indicate that the epidemic system is self-inclusion and zero population flux across the boundary. d 1 > 0 and d 2 > 0 are the diffusion rates of susceptible and infected individuals, respectively. And the positive number A is the recruitment rate of susceptible individuals. The positive coefficients μ and μ 1 represents natural mortality and mortality caused by diseases, respectively. Function h(I ) represents the treatment strategy which λ is positive coefficient. We give some properties of function h(I ) as follows. (H1) h(I ) is continuous except in a cluster of countable isolated points ρ k , where h(ρ + k ) and h(ρ − k ) represents the right and left limits, respectively, with h(ρ + k ) > h(ρ − k ). Besides, h(I ) has a finite discontinuous points in any compact interval [0, +∞). Due to the function h(I ) is not continuous, here, we applying differential inclusion theory instead of theories and methods in ordinary differential equations, So the system (1.1) can be written as the following differential inclusion When saturation coefficient m = 0, the nonlinear saturated incidence rate become the linear incidence rate, which is corresponding to the system in literature [17] . To the best of the authors' knowledge, there are few analysis results about diffusive SIR epidemic model with saturated incidence rate and discontinuous term in the open literature. The aim of this paper is to investigate whether we can conclude some different dynamic behaviour of the system (1.1) owing to the saturated incidence rate, and whether saturation coefficient m will affects the global stability of the endemic equilibrium of system (1.1). The rest of this paper is arranged as follows. Firstly, we mainly investigate the existence and some properties of the solution of the diffusive epidemic system (1.1) in Sect. 2. Secondly, we discuss about the existence of the endemic equilibrium of the system in Sect. 3. Thirdly, by constructing different suitable Lyapunov functions, we prove the global asymptotic stability of the disease free equilibrium(DFE) and endemic equilibrium(EE) respectively in Sect. 4. Finally, we give a brief conclusion in Sect. 5. In this section, we are concerned with the existence and properties of the solution of the system (1.1). Firstly, we assume that the initial values S(x, 0), I (x, 0) of system (1.1) satisfy the following condition. ( (2.1) Proof Firstly, to prove the boundedness of solutions (S(x, t), I (x, t)), we give the invariant rectangle Obviously, S 0 (x) and I 0 (x) are closed for the rectangle . The vector field of the system (1.1) is given all points on the rectangle point inside. Secondly, through the above analysis, we have the following conclusions: (i) On the left side of the first quadrant invariant rectangle with S = S 1 , I 1 < I < I 2 , by the definition of S 1 , it satisfies the following estimate: (ii) On the right side of the the first quadrant invariant rectangle with S = S 2 , I 1 < I < I 2 , by the definition of S 2 , it satisfies the following estimate: On the bottom side of the the first quadrant invariant rectangle with I = I 1 , S 1 < S < S 2 , by the definition of I 1 , it satisfies the following estimate: (iv) On the left side of the the first quadrant invariant rectangle with I = I 2 , S 1 < S < S 2 , by the definition of I 2 , it satisfies the following estimate: Finally, by virtue of the definition on [20] , in view of the above discussions, we can conclude that is the invariant rectangle of the vector field (2.3). Thus, we can choose M 1 = min{S 1 , I 1 } and M 2 = max{S 2 , I 2 }, which completes the proof. Now, we give the following definition. Firstly, we denote (2.9) Obviously, we can obtain that the map ϕ 1 (S, I ) is bounded. By assumptions (H1), (H2) and (H3), with combining the above discussions, it is easy to know ϕ 2 (S, I ) is an upper semi-continuous bounded set-valued mapping with nonempty compact convex values. Next, we denote Let X be a real Banach space, U = (S, I ) be the solution of system (1.1) with initial value U (0) = (S 0 , I 0 ) > 0, and we defined : D( ) X → X is the infinitesimal generator of a C 0 -semigroup of linear contractions I (t). Furthermore, we defined : [0, T ] × X → X be a function which measurable in t and Lipschitz continuous in X , uniformly with respect to t ∈ [0, T ]. Then we have the following conclusion: (i) If U 0 ∈ X , there exists a unique weak solution of system (2.9). (ii) If X is a Hilbert space, is a self-adjoint and dissipative on X with U (0) ∈ D( ), and we can obtained which shows that the weak solution of system (2.9) is actually a strong solution. [29] , if the condition T max < +∞, then But by Lemma 2.1, where the invariant is an L ∞ a priori bound for the solution (S, I ) of the system (1.1), it is a contradiction which implies T max = ∞, i.e., so for all (x, t) ∈ ×[0, +∞), the solution of system (1.1) exists and bounded. Which completes the proof. In this section, we are concerned about the existence of disease free equilibrium(DFE) and endemic equilibrium(EE) of the system (1.1). We using the analysis method in [17] . Firstly, when h(I ) = 0, the DFE E 0 (S * 0 , Through a simple transformation, we define g(I ) as follows (3.2) We are divided into the following three steps to prove. Step 1. We prove that system Step 3. We prove thatĪ is the unique positive solution of system (3.1). Let I 1 =Ī is a solution of (3.1), and I 2 = I 1 is another positive solution of (3.1), then, there exists γ 1 ∈ co[h(I 1 )] and γ 2 ∈ co[h(I 2 )], so we have From the monotonicity of h(I ), it implies that H = γ 1 −γ 2 I 1 −I 2 ≥ 0. However, after subtraction of the two equations of (3.3), we obtain which is a contradiction. The proof of the lemma is completed. A direct result of Lemma 3.1 is the following theorem of endemic equilibrium. when Aβ > μμ 1 + μλ holds, I * is decreasing with the saturation coefficient m, it means when m → +∞, the infective individuals may cannot persist. In this section, we discussed the global stability of disease free equilibrium and endemic equilibrium in the invariant rectangle , respectively. Firstly, the global stability of DFE E 0 is discussed as follows. So from the above discussions, by calculating ∇V (S, I ) · w, we can obtain that (4.4) Owing to the homogeneous Neumann boundary condition, we obtain Furthermore, when Aβ ≤ μμ 1 , we have Therefore, we can obtain that Thus, the disease free equilibrium E 0 is stable, and when (S, By the Lasalle invariance principle [31] , E 0 is globally asymptotically stable for system (1.1). The proof is completed. Aβ > μμ 1 + μλ hold, then EE E * = (S * , I * ) of system (1.1) is globally asymptotically stable. Proof We define is a smooth function, based on assumptions (H1) and (H2), there exist a function γ 2 ∈ co[h(I )], we have By calculating ∇V (S, I ) · w, we can obtain (4.10) Owing to the homogeneous Neumann boundary condition, we can get (4.11) (4.14) As a result, combined with (4.11), we have Thus, the endemic equilibrium E * is stable, and when (S, I ) = (S * , I * ), d dt V 2 (t) = 0. So the singleton E * is the maximum compact invariant set in = {(S, I )| d dt V 2 (t)| = 0}. By the Lasalle invariance principle [31] , E * is globally asymptotically stable for system (1.1). The proof is completed. In this section, we show numerical simulations regarding our model to illustrate and support the theoretical results of the previous sections. We illustrate the DFE is globally asymptotically stable if the basic regeneration number R 0 is less than unity in Fig.1 Then we obtain R 0 = 0.667 < 1, which implies that DFE E 0 = (4, 0) is globally asymptotically stable by Theorem 4.1 where disease will extinct. And for Fig.2 , we show that the EE is globally asymptotically stable, for Then we obtain Aβ > μμ 1 + μλ, which implies the EE is globally asymptotically stable by Theorem 4.2 where disease spread in the human world. In this paper, we investigate the dynamic of a diffusive SIR epidemic model under discontinuous treatments. Due to the discontinuous term, the existence of strong solution or weak solution of system(1.1) is proved under the framework of differential inclusion. Based on the differential inclusions theory, we analysis the existence of endemic equilibrium of the system (1.1). Moreover, by constructing different suitable Lyapunov functions, we investigate the global asymptotic stability of the disease free equilibrium(DFE) and the endemic equilibrium(EE), respectively. Compared to literature [17] , a different incidence rate g(I ) = β I 1+m I is considered in the system (1.1). When saturated incidence rate m = 0, the nonlinear saturated incidence rate become the linear incidence rate, which is corresponding to the system in [17] . When saturated incidence rate m → +∞, we conclude that the infective individuals cannot persist in Remark 3.1. This result seems to be consistent with realistic intuition : the more behavioral changes of susceptible individuals or the inhibition of crowding effect of infected individuals, the better disease control. When saturated incidence rate 0 < m < +∞, from Theorem 4.2, it is shown that saturated incidence rate m may affect the global stability of the endemic equilibrium. Furthermore, if linear incidence rate is considered and the discontinuous treatments term does not exist in the system(1.1), that is m = 0 and λ = 0, then from Theorem 4.1 and Theorem 4.2 we can easily define the basic reproduction number R 0 = Aβ μμ 1 . When R 0 < 1, the DFE is globally asymptotically stable, and when R 0 > 1, the EE is globally asymptotically stable. However, in the system (1.1), due to nonlinear incidence rate and the discontinuous term, it follows from Fig.1 , we can obtain that when R 0 < 1, the DFE is also globally asymptotically stable, and when R 0 > 1 + λ μ 1 , it follows from Fig.2 , the EE is also globally asymptotically stable. But, when 1 < R 0 < 1 + λ μ 1 , whether the endemic equilibrium is globally asymptotically stable is still unknown, which will be considered in our future work. A nonlinear epidemiological model considering asymptotic and quarantine classes for SARS CoV-2 virus Modelling the role of optimal social distancing on disease prevalence of COVID-19 epidemic Mathematical model for COVID-19 management in crowded settlements and high-activity areas Global asymptotic stability of an SIR epidemic model with distributed time delay On the nonlinear stability of an epidemic SEIR reaction-diffusion model Dynamics and asymptotical profiles of an age-structured viral infection model with spatial diffusion Pattern dynamics of a reaction-diffusion predator-prey system with both refuge and harvesting Dynamic behaviour of a reactiondiffusion predator-prey model with both refuge and harvesting Stationary patterns induced by self-and cross-diffusion in a Beddington-DeAngelis predator-prey model Consequences of refuge and diffusion in a spatiotemporal predator-prey model Existence of complex patterns in the Beddington-DeAngelis predator-prey model Coexistence states of a Lotka Volterra cooperative system with cross diffusion Impact of discontinuous harvesting on fishery dynamics in a stock-effort fishing model Dynamics of a discontinuous computer worm system Impact of discontinuous treatments on disease dynamics in an SIR epidemic model Global stability of a diffusive predatorprey model with discontinuous harvesting policy Global dynamics of a controlled discontinuous diffusive SIR epidemic system Optimal control of harvest and bifurcation of a prey-predator model with stage structure Dynamical analysis for a fractionalorder prey-predator model with Holling III type functional response and discontinuous harvest Turing patterns in the Lengyel-Epstein system for the CIMA reaction A generalization of the Kermack-McKendrick deterministic epidemic model Analysis on a diffusive SIS epidemic model with saturated incidence rate and linear source in a heterogeneous environment Qualitative analysis on an SIRS reaction diffusion epidemic model with saturation infection mechanism Asymptotic profile of endemic equilibrium to a diffusive epidemic model with saturated incidence rate A priori L ∞ estimates for solutions of a class of reaction-diffusion systems A new two-grid method for expanded mixed finite element solution of nonlinear reaction diffusion equations Invariant sets for weakly coupled parabolic and elliptic systems On p Laplacian differential inclusions-Global existence, compactness properties and asymptotic behavior Global existence and boundedness in reaction-diffusion systems Stability and stabilization of discontinuous systems and nonsmooth Lyapunov functions Stability theory for ordinary differential equations Acknowledgements This work is supported by the Natural Science Foundation of China (11672074), and the Natural Science Foundation of Fujian Province (2018J01655).Author Contributions All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.Funding The research was supported by the Natural Science Foundation of China (11672074), and the Natural Science Foundation of Fujian Province (2018J01655). The authors declare that there is no conflict of interests.