Research Article  Open Access
Ghassane Benrhmach, Khalil Namir, Jamal Bouyaghroumni, "Modelling and Simulating the Novel Coronavirus with Implications of Asymptomatic Carriers", International Journal of Differential Equations, vol. 2020, Article ID 5487147, 9 pages, 2020. https://doi.org/10.1155/2020/5487147
Modelling and Simulating the Novel Coronavirus with Implications of Asymptomatic Carriers
Abstract
The World Health Organization declared that the total number of confirmed cases tested positive for SARS‐CoV‐2, affecting 210 countries, exceeded 3 million on 29 April 2020, with more than 207,973 deaths. In order to end the global COVID‐19 pandemic, public authorities have put in place multiple strategies like testing, contact tracing, and social distancing. Predictive mathematical models for epidemics are fundamental to understand the development of the epidemic and to plan effective control strategies. Some hosts may carry SARS‐CoV‐2 and transmit it to others, yet display no symptoms themselves. We propose applying a model (SELIAHRD) taking in consideration the number of asymptomatic infected people. The SELIAHRD model consists of eight stages: Susceptible, Exposed, Latent, Symptomatic Infected, Asymptomatic Infected, Hospitalized, Recovered, and Dead. The asymptomatic carriers contribute to the spread of disease, but go largely undetected and can therefore undermine efforts to control transmission. The simulation of possible scenarios of the implementation of social distancing shows that if we rigorously follow the social distancing rule then the healthcare system will not be overloaded.
1. Introduction
By the end of 2019, the city of Wuhan, China, and the rest of the world have seen the apparition of a deadly virus, the most recently discovered coronavirus, named “SARSCoV2” [1]. Few weeks later, COVID19 spread widely in China and other countries, making the World Health Organization (WHO) to declare a global pandemic in March 2020 [2]. There have been 3,090,445 confirmed cases and 217,769 confirmed deaths in 212 countries around the globe until 30 April 2020 according to the World Health Organization.
In order to stop the rapid spread of the coronavirus, many countries have introduced widespread physical distancing measures, populationlevel movement restrictions, and health education knowledge [3, 4]. However, one of the biggest problems is the existence of asymptomatic infection (having very mild symptoms) among the infected people, and unaware of their contagious ability increases the transmission rate [5]. Recently, and to keep up with the biologists, many mathematical models have been developed to enable a more effective response to this hazard. In [6], Peng et al. have numerically studied a SEIR model that includes intrinsic impact of hidden exposed infected population to describe the epidemic of COVID19. They have estimated the time of the outbreak of this epidemic in several cities in China. However, they did not incorporate the transmission paths and speed of such an epidemic. In [7], M. Peirlinck et al. have introduce the effect of asymptomatic infection on the outbreak of Covid19 by proposing a SEIIR model and focusing in decomposing the infectious population into a symptomatic and an asymptomatic group and assuming a similar transmission rate β for both groups. They found that the large asymptomatic population is associated with a high risk of Covid19 spread. On the other hand, the epidemic progresses across countries, and hospital staff are faced with a massive influx of patients and a large number of deaths, causing a huge problem in resource management. Motivated by these reasons, we propose a SELIAHRD model including eight states (Susceptible, Exposed, Latent, Symptomatic, Asymptomatic, Hospitalized, Recovered, and Death) taking into consideration the total population as susceptible, age, and population structure is not considered [7, 8]. This model was developed based on studies of COVID19 and the H1N1 influenza epidemic model proposed by Lopez et al. [9] however, in this model, they suppose the existence of vaccinated people and used Generalized regression neural network for simulation without taking in consideration the social distancing scenarios.
This article begins with a section that describes the proposed epidemic model and gives an idea about each compartment of the model. The second section presents a discussion of the wellposedness and equilibria of the proposed model. The stability of the diseasefree equilibrium, the use of nextgeneration operator to get the basic reproductive number , and the proof of the global stability of the diseasefree equilibrium are presented in the third section. The results of the numerical simulation and their discussion are given in the fourth section.
2. Proposed Epidemic Model
As per the recent report from the WHO, a significant number of COVID19 patients have very mild symptoms. This means the existence of asymptomatic infections, which increase the transmission rate dramatically in some countries.
To take into consideration the effect of undetectable infected persons, we propose the SELIAHRD model including eight compartments shown in Figure 1.
The susceptible state is a person without the disease at time . The exposed state refers to those near the infected person and has higher risk to get infected. The exposed persons who got infection but not yet infectious are called latent people. From the latent compartment, we can classify the people as symptomatic and asymptomatic people; the symptomatic people will be admitted to the hospital. Asymptomatic and hospitalized people have the probability to go either to the recovered compartment or to the death compartment.
The model equations are listed below (system 1):
Table 1 shows the definition of symbols in the proposed model.

3. WellPosedness and Equilibria
The existence, the positivity, and the boundedness of solutions of the proposed system (1) need to be proved to ensure that the model has a mathematical and biological meaning.
System (1) can be rewritten as follows:whereand is a function mapping into itself, defined by
By the fundamental theory of functional differential equations, system (1) has a unique solution with respect to the initial data such that
If we put , then the following theorems hold.
Theorem 1. The following biological feasible region of system (1) is positively invariant and attracting.
Theorem 2. Let and the initial conditions satisfy , then the solution of system (1) are positive for all .
Proof. From the first equation of system (1), we haveHence, is nonnegative for all . From the other equations of system (1), we setTherefore, are all nonnegative for all .
From the seventh equation of system (1), we can easily deduce the positivity of for all .
Hence, the positivity of solution has been proved.
4. Stability of DiseaseFree Equilibrium
This section is devoted to studying the stability of the diseasesfree equilibrium of system (1).
The equilibrium of the model is obtained by setting , and the given system has a diseasefree equilibrium given by . In the proposed model, we have multiple types of infected individuals, and then to investigate the stability of the diseasefree equilibrium, we will define the basic reproductive number () using the nextgeneration operator [10], who define as the spectral radius of the nextgeneration operator.
Theorem 3. If , then DFE (diseasefree equilibrium) is locally asymptotically stable (LAS). If , then the DFE is unstable.
Let be the vector of infected classes, such as infectious, exposed, carrier, etc., and be the vector of uninfected classes, such as susceptible, recovered, etc.
Let denote the diseasefree equilibrium, that is, , where is the vector of new infection rates (flows from to ) and is the vector of all others rates (not a new infection). These rates include flows from to (for instance, recovery rates), flows within , and flows leaving the system (for instance, death rates). For each compartment, inflow in is negative and outflow in is positive.The formula of the basic reproductive number () is defined by the maximum eigenvalue of , where and , and is called the nextgeneration matrix.where
Multiplying and together and calculating the spectral radius, we have
4.1. Global Stability Condition for DiseaseFree Equilibrium
In this section, we present two conditions, and if they are satisfied, they allow us to say that the diseasefree equilibrium is globally asymptotically stable. First, system (1) must be written in the following form:
() for is globally asymptotic stable (GAS), () for for , where is an Metzler matrix (Mmatrix, the offdiagonal elements of are nonnegative). If the system satisfies these two conditions, then the following lemma holds.
Lemma 1. The diseasefree equilibrium of system (2) is globally asymptotically stable (GAS) provided that and assumptions () and () are satisfied.
Now, we state the following theorem.
Theorem 4. The diseasefree equilibrium of system (1) is globally asymptotically stable if .
Proof. Let and . We will haveAt the point is globally asymptotically stable for . On the other hand,In addition, this follows that ,The offdiagonal elements of are nonnegative, i.e., an Mmatrix and Since and thus . So, the conditions () and () of Lemma 1 are satisfied, and then this completes the global asymptotic stability of the diseasefree equilibrium of the proposed system (1) for .
5. Numerical Results and Discussion
In this section, we present the simulation of our SELIAHRD model in two different scenarios. In order to adjust the social distance, we multiply the transmission rate by a social distancing factor in order to obtain the two scenarios.(a)The mean incubation period was 5.2 days ( confidence interval) [9]. Therefore, .(b)There is a mean 5day delay from symptom onset to detection/hospitalization of a case. The duration from illness onset to the first medical visit for the 45 patients with illness onset before January 1 was estimated to have a mean of 5.8 days [11]. So, we set .(c)The recovery rate is determined by the average duration of recovery D of infection. After this period, they enter the recovery phase.(d)Since there were no data on the proportion of asymptomatic infection of the virus, we simulated and .(e)Symptomatic and asymptomatic infection rates are defined as follows: (symptomatic infected)/susceptible. (asymptomatic infected)/susceptible.
In this first simulation, we assume that there is no social distancing and a total population of 1.5 million. In this scenario, Figures 2 and 3 show that the model predicts in the next 85 days the total symptomatic infected person will reach a peak of 31000 and the healthcare system will be overload (only 15,869 of them get treatment). As a result, the model predicts that the peak of deaths would be 15,665.
In the second simulation, we assume that of population obey the rule of social distancing.
As can be seen in Figures 4 and 5, the peak of symptomatic infected cases drops significantly to 50 and asymptomatic infected cases to 6, making the healthcare system not overloaded at all and only around 5 deaths.
The simulation results show that after 150 days, with the presence of asymptomatic infection the basic reproduction number of coronavirus was 3.25 if there is no social distancing and the healthcare system will be overloaded. On the other hand, if the lockdown is fully operational ( of population respect of social distancing) even with the presence of asymptomatic infection, the basic reproduction number becomes after 150 days and the healthcare system will not be overloaded. The paper [12] proposed a SEIAR model to simulate the epidemic dynamics of “SARSCOV2” and demonstrate the efficiency of quarantine and government measure but didn’t take in consideration the velocity of transition between different stages and the effect of lockdown on the healthcare system.
6. Conclusions
The world continues to pay a heavy price as the COVID19 pandemic spreads. This paper proposes the mathematical model SELIAHRD (Susceptible, Exposed, Latent, Symptomatic, Asymptomatic, Hospitalized, Recovered, and Death), with two kinds of infected people, that is, symptomatic and asymptomatic. We assume that asymptomatic people will not be hospitalized and have the probability to go either recovered or dead. We calculate the basic reproduction number R_{0} and dynamics of diseasefree equilibrium for the “SARSCoV2” using the proposed model.
The SELIAHRD model will provide in the presence of asymptomatic infected the number of hospitalized and of deaths for the public authorities to develop the prevention policies (cities lockdown, social distancing, testing, contact tracing, etc.). The simulation shows the importance of the social distancing in making the number of infected and deaths decrease significantly even with the presence of asymptomatic infection and not overloading the healthcare system.
Data Availability
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.
References
 World Health Organization, Pneumonia of Unknown Cause – China, World Health Organization, Geneva, Switzerland, 2020.
 World Health Organization, WHO DirectorGeneral’s Opening Remarks at the Media Briefing on COVID19—11 March 2020, World Health Organization, Geneva, Switzerland, 2020.
 A. R. Tuite, D. N. Fisman, and A. L. Greer, “Mathematical modelling of COVID19 transmission and mitigation strategies in the population of Ontario, Canada,” Canadian Medical Association Journal, vol. 192, no. 19, pp. E497–E505, 2020. View at: Publisher Site  Google Scholar
 Z. Wu and J. M. Mcgoogan, “A Characteristics of and important lessons from the coronavirus disease 2019 (COVID19) outbreak in China,” JAMA, vol. 323, 2020. View at: Publisher Site  Google Scholar
 Y. Bai, L. Yao, T. Wei et al., “Presumed asymptomatic carrier transmission of COVID19,” JAMA, vol. 323, no. 14, p. 1406, 2020. View at: Publisher Site  Google Scholar
 L. Peng, W. Yang, D. Zhang, C. Zhuge, and L. Hong, “Epidemic analysis of COVID19 in China by dynamical modeling,” 2002, http://arxiv.org/abs/2002. View at: Google Scholar
 M. Peirlincka, K. Linkaa, F. Sahli Costabalb et al., Visualizing the Invisible: The Effect of Asymptomatic Transmission on the Outbreak Dynamics of COVID19, Cold Spring Harbor Laboratory Press, New York, NY, USA, 2020.
 B. H. Sze and H. H. Ying, “On the role of asymptomatic infection in transmission dynamics of infectious diseases,” Bulletin of Mathematical Biology, vol. 70, no. 1, pp. 134–155, 2008. View at: Publisher Site  Google Scholar
 D. Lopez, G. Manogaran, and J. Mohan, “Modelling the H1N1 influenza using mathematical and neural network approaches,” Biomedical Research, vol. 28, no. 8, 2017. View at: Google Scholar
 P. Van Den Driessche and J. Watmough, “Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, no. 12, pp. 29–48, 2002. View at: Publisher Site  Google Scholar
 Q. Li, X. Guan, P. Wu et al., “Early transmission dynamics in Wuhan, China, of novel coronavirusinfected pneumonia,” The New England Journal of Medicine, vol. 382, pp. 1–9, 2020. View at: Google Scholar
 L. Pribylova and V. Hajnova, “SEIAR model with asymptomatic cohort and consequences to efficiency of quarantine government measures in COVID19 epidemic,” 2020, https://arxiv.org/abs/2004.02601. View at: Google Scholar
Copyright
Copyright © 2020 Ghassane Benrhmach et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.