Modelling and Simulating the Novel Coronavirus with Implications of Asymptomatic Carriers

The World Health Organization declared that the total number of conﬁrmed cases tested positive for SARS-CoV-2, aﬀecting 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 eﬀective 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 eﬀorts 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.


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 "SARS-CoV-2" [1]. Few weeks later, COVID-19 spread widely in China and other countries, making the World Health Organization (WHO) to declare a global pandemic in March 2020 [2]. ere 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, population-level 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 COVID-19. ey 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 Covid-19 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. ey found that the large asymptomatic population is associated with a high risk of Covid-19 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].
is model was developed based on studies of COVID-19 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.
is article begins with a section that describes the proposed epidemic model and gives an idea about each compartment of the model. e second section presents a discussion of the well-posedness and equilibria of the proposed model. e stability of the disease-free equilibrium, the use of next-generation operator to get the basic reproductive number R 0 , and the proof of the global stability of the disease-free equilibrium are presented in the third section. e results of the numerical simulation and their discussion are given in the fourth section.

Proposed Epidemic Model
As per the recent report from the WHO, a significant number of COVID-19 patients have very mild symptoms.
is 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.
e susceptible state is a person without the disease at time t. e exposed state refers to those near the infected person and has higher risk to get infected. e 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.
e model equations are listed below (system 1):  Table 1 shows the definition of symbols in the proposed model.

Well-Posedness and Equilibria
e 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: where and P is a C 1 function mapping R 8 into itself, defined by  International Journal of Differential Equations 3 By the fundamental theory of functional differential equations, system (1) has a unique solution (S(t), E(t),

L(t), I(t), A(t), H(t), R(t), D(t))
with respect to the initial data Z 0 such that If we put S + E + L + I + A + R + H + D ≤ N, then the following theorems hold.
Proof. From the first equation of system (1), we have Hence, S(t) is nonnegative for all t ≥ 0. From the other equations of system (1), we set erefore, From the seventh equation of system (1), we can easily deduce the positivity of R(t) for all t ≥ 0.
Hence, the positivity of solution has been proved. □

Stability of Disease-Free Equilibrium
is section is devoted to studying the stability of the diseases-free equilibrium U 0 of system (1).
e equilibrium of the model is obtained by setting dS (t)/dt � dE(t)/dt � dL(t)/dt � dI(t)/dt � dA(t)/dt � dH (t)/dt � dR(t)/dt � dD(t)/dt � 0, and the given system has a disease-free equilibrium given by U 0 � (N, 0, 0, 0, 0, 0, 0, 0). In the proposed model, we have multiple types of infected individuals, and then to investigate the stability of the disease-free equilibrium, we will define the basic reproductive number (R 0 ) using the next-generation operator [10], who define R 0 as the spectral radius of the nextgeneration operator.
Let Y be the vector of infected classes, such as infectious, exposed, carrier, etc., and X be the vector of uninfected classes, such as susceptible, recovered, etc. Y) is the vector of new infection rates (flows from X to Y) and V(X, Y) is the vector of all others rates (not a new infection). ese rates include flows from Y to X (for instance, recovery rates), flows within Y, and flows leaving the system (for instance, death rates). For each compartment, inflow in V is negative and outflow in V is positive. e formula of the basic reproductive number (R 0 ) is defined by the maximum eigenvalue of FV − 1 , where F � (zF/zY) (X * ,0) and V � (zV/zY) (X * ,0) , and FV − 1 is called the next-generation matrix. where Multiplying F and V − 1 together and calculating the spectral radius, we have International Journal of Differential Equations

Global Stability Condition for Disease-Free Equilibrium.
In this section, we present two conditions, and if they are satisfied, they allow us to say that the disease-free equilibrium is globally asymptotically stable. First, system (1) must be written in the following form: (H 1 ) for dX/dt � W(X, 0), X * is globally asymptotic stable (GAS), ( where Q � (zG/zY) U 0 is an Metzler matrix (M-matrix, the off-diagonal elements of Q are nonnegative). If the system satisfies these two conditions, then the following lemma holds.

Proof. Let X � (S, R) and Y � (E, L, I, A, H, D). We will have
At the point (X, 0), W(X, 0) � (0), X * � (N, 0) is globally asymptotically stable for dX/dt � W(X, 0). On the other hand, In addition, this follows that Q � (zG/zY) U 0 , 6 International Journal of Differential Equations e off-diagonal elements of Q are nonnegative, i.e., an M-matrix and Since S ≤ N, 1 − (S/N) ≥ 0, and thus G(X, Y) ≥ 0. So, the conditions (H 1 ) and (H 2 ) of Lemma 1 are satisfied, and then this completes the global asymptotic stability of the disease-free equilibrium of the proposed system (1) for R 0 < 1.

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) e mean incubation period was 5.2 days (95% confidence interval) [9]. erefore, μ 1 � 0.1923. (b) ere is a mean 5-day delay from symptom onset to detection/hospitalization of a case. e 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 μ 3 � 0.1724.
(c) e recovery rate c 1 � 0.15 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 μ 1 � μ 2 and μ 4 � 0.08. (e) Symptomatic and asymptomatic infection rates α 1 and α 2 are defined as follows: 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 90% 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 International Journal of Differential Equations asymptomatic infected cases to 6, making the healthcare system not overloaded at all and only around 5 deaths.
e simulation results show that after 150 days, with the presence of asymptomatic infection the basic reproduction number R 0 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 (90% of population respect of social distancing) even with the presence of asymptomatic infection, the basic reproduction number becomes R 0 � 0.68 after 150 days and the healthcare system will not be overloaded. e paper [12] proposed a SEIAR model to simulate the epidemic dynamics of "SARS-COV-2" 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.  International Journal of Differential Equations

Conclusions
e world continues to pay a heavy price as the COVID-19 pandemic spreads. is 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 disease-free equilibrium for the "SARS-CoV-2" using the proposed model. e 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.). e 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
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.