Study of transmission dynamics of novel COVID-19 by using mathematical model

In this research work, we present a mathematical model for novel coronavirus-19 infectious disease which consists of three different compartments: susceptible, infected, and recovered under convex incident rate involving immigration rate. We first derive the formulation of the model. Also, we give some qualitative aspects for the model including existence of equilibriums and its stability results by using various tools of nonlinear analysis. Then, by means of the nonstandard finite difference scheme (NSFD), we simulate the results for the data of Wuhan city against two different sets of values of immigration parameter. By means of simulation, we show how protection, exposure, death, and cure rates affect the susceptible, infected, and recovered population with the passage of time involving immigration. On the basis of simulation, we observe the dynamical behavior due to immigration of susceptible and infected classes or one of these two.


Introduction
Recently, the whole world has been suffering due to a novel coronavirus pandemic. It was named novel coronavirus infectious disease  which was claimed to outbreak first in Wuhan, central China (see [1]). Novel coronavirus-19 is a new chain of corona group of viruses that had not been identified in human history before December 2019. For the first time COVID-19 was found in Wuhan, China in December 2019 and has spread to various urban areas in China as well as round about 196 different countries of the world. It has since been declared an outbreak by World Health Organization (WHO). According to the data reported by WHO (World Health Organization), on 11 June 2020, the reported laboratory confirmed that the number of affected humans reached more than 7.5 million including more than 0.425 million death cases recorded. Some researchers have also claimed that there are other sources of this corona virus including dogs, pangolin, etc. As per recoded data, the death rate is different in different countries. Currently the highest death rate has been observed in Europe, USA, and Iran. The number of confirmed cases has been growing on a very fast track on a daily basis, and it has been declared a worldwide pandemic disease.
On 31 December 2019, the WHO reported a novel corona virus (2019-nCoV) in Wuhan City, Hubei Province of China in humans, see [2,3]. It was named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by the International Committee on Taxonomy of Viruses on 11 February 2020 (for details, we refer to [4][5][6][7][8][9][10]). Firstly, this outbreak was identified in Wuhan with most early cases being reported in the city and later spread to other countries at an alarming rate and became a lethal disease. There are different schools of thought behind the origin of COVID-19: some say that it might be of bat origin (see [11]), some say that it might be related to a seafood market exposure (see [12]). International travel of any form has been a potential reason for the fast spread of COVID-19 [2,[12][13][14]. So, immigration has a severe impact on the severity of spreading of COVID-19. It has been stated (fact) that the origin of novel COVID-19 is the transmission from animal to human as many infected cases claimed that they had been to a local fish and wild animal market in Wuhan in November [15]. Soon, some researchers confirmed that the transmission also happened person to person (see [2]). In the present situation this pandemic has produced a very harmful effect on health, economics, and social life of the whole globe. In the whole world researchers, policy makers, and doctors are struggling to control this serious pandemic so that the lives of maximum people may be secured. They observed this disease from their own point of view. Also it is a fact that most people infected with novel COVID-19 will experience mild to moderate respiratory illness. Common symptoms are fever, tiredness, dry cough, and throat infection. Some people may experience aches and pains, nasal congestion, runny nose, sore throat, diarrhoea, etc.
Since mathematical models are powerful tools to understand the dynamics of real world phenomena, particulary the transmission of an infectious disease, in literature large numbers of mathematical models of infectious disease have been studied, we refer to a few of them [16][17][18][19][20][21][22]. Also the area of modeling has been extended recently to noninteger order and nonlocal derivatives of fractional order [23,24]. By using mathematical models for understanding the transmission dynamics of a disease can help the researchers to make future prediction and to adopt some precautionary measures to save maximum population from being lost. Also the mentioned tools help to make strategies to control or eliminate the disease from society. Same as the case of current novel COVID-19, which has been studied from different aspects in the last few months (for details, see [25][26][27][28][29]). Therefore, motivated by the aforesaid discussion, we observed that immigration has major roles in spreading the current disease in our society. It has been observed that due to immigration, this disease has spread in the whole globe within two to three months. Therefore in this work we construct a modified SIR type model involving immigration rate to investigate the transmission dynamics of the aforementioned disease. For numerical simulation, we use the nonstandard finite difference (NSFD) scheme. The concerned method is an efficient and powerful method to find numerical solutions to many nonlinear problems. Therefore various researchers have used this method for the numerical simulation of many bathetical models (for instance, see [17,30] and the references therein).

Model formulation
This part of the paper is devoted to constructing the mathematical model for our proposed problem. We take here three compartments: susceptible S(t), infected I(t), and recovered R(t). We construct the required model under convex incidence rate which is assumed to be a convex function with respect to the infective class due to host population. The benefit of using convex incidence rate is that it corresponds to an increased rate of infection  Table 1 The physical interpretation of the parameter Parameters The physical interpretation The recruitment rate , where both C, γ are positive constants. This is an interesting example for nonlinear incidence rate already used by some authors [17,31,32]. The flow chart of the model is shown in Fig. 1. The dynamics of the population are described by the following differential equations: The parameters involved in model (1) are described as in Table 1.
First, for the equilibrium of model (1), we consider its existence. Corresponding to some values of parameters, there exists a disease-free equilibrium for system (1) denoted by E 0 = (a/μ, 0, 0). To compute the nonnegative equilibrium, we have To find the basic reproduction number R 0 , let x = (S(t), I(t)) in model (1). Then for the disease-free equilibrium Jacobian of F is and Jacobian of V to deduce the disease-free equilibrium is given by Hence, for model (1), by simple calculation, we have Hence the basic reproduction number (reproductive rate) R 0 is From (2), we clearly observe that (i) There are no positive equilibria of model (6) if R 0 ≤ 1; (ii) A unique positive equilibrium also known as endemic equilibrium The endemic equilibrium is given by , The value of Ω is given as Next, we will elaborate on the characteristics of these equilibria and a global mathematical analysis of system (1).

Dynamical behavior of the model
To elaborate the dynamic of system (1), we have the following lemma. Then (4) implies that Hence for (5) we present the general solution as which completes our conclusion. Now we reduce system (1), because obviously the limit set of (1) on plane S(t) + I(t) + R(t) = a μ has a limit set: We have the following theorem with regards to the nonexistence of cyclical shells in system (6), which shows the nonexistence of cyclical shells of system (1) by Lemma 1.

Theorem 1
There do not exist nontrivial periodic orbits corresponding to system (6).

D I(t), R(t) = 1 + γ I(t) CI(t) .
Then we have By adding all equations of (7), we have Hence which proves the conclusion of the theorem.
To study S 0 disease-free equilibrium and its properties, and also the endemic equilibrium S * , we recall (6) with One can obtain from (6) with p = β μ + α , Note Keep in mind that (0, 0) may be obtained from system (9). In fact the disease-free equilibrium S 0 of system (1) and (x * , y * ) of system (9) is the unique positive equilibrium, which is, in fact, the endemic equilibrium S * of system (1) under the condition n -B < 0 with x * = B-n q-1 and y * = px * . At first glance, we investigate for (0, 0) the stability and topological type trivial equilibrium. At the point (0, 0), the Jacobian matrix of system (9) is given by The dynamic of system (9) is equivalent to (11). If Bn = 0, then there exists a small neighborhood N 0 of (0, 0).
From (11) (0, 0) is a saddle node. The next results is important.

Theorem 2
The trivial equilibrium point of system (1) possesses the following properties: (i) As a result the system has a hyperbolic saddle if n < B.
(ii) As a result the system has a saddle node if n = B.
(iii) As a result the system has a stable hyperbolic node if n > B.
Proof When n -B < 0, we study the topological type of endemic equilibrium (x * , y * ) and stability. From (9) at (x * , y * ), we have the Jacobian matrix Thus det(M 1 ) has not a unique sign due to Relation (12) tells that S1 > 0 yields det(M1) > 0 and (x * , y * ) is a node (focus or center). Also, for the stability of (x * , y * ), one can find the given results.

Theorem 3
The equilibrium (x * , y * ) of system (9) is locally stable in a unique way, and also it has a stable node if n -B < 0.
The following theorem summarizes the results for the stability of the original system (1) in terms of the basic reproduction number. (2) we define R 0 .

Numerical results and discussion
We present numerical simulation for system (6) with the used values. We take two different sets of values for immigration parameters involved in model (6) and real data of Wuhan city of China to simulate the results.
According to the NSFD scheme, the first equation of our considered model (6) may be expressed as which is decomposed in the NSFD scheme as follows: Like (14), we can decompose model (6) in the NSFD scheme and write the whole system as follows: S j+1 = S j + h a -CI j (t)S j (t) 1 + γ I j (t) -μS j (t) + αR j (t) , The birth rate of infection 0.1243 [28] μ Natural death 0.002 [28] δ Death due to corona 0.05 [28] b The immigration rate rate 0.0205 [29] β Corona infection recovery rate 0.09871 [27] C Infection rate 0.580 [28] γ Rate at which recovered individuals lose immunity 0.0003 [29] α Rate of recovery 0.854302 [27]   During the first thirty days in the presence of excessive rate of immigration the susceptible population is decreasing as shown in Fig. 2. When the immigration rate is high, the decline in the population of uninfected (susceptible) people is observed, because they are exposed to infection. Hence the higher the immigration rate, the faster the growth rate of the infected population and vice versa. As a result, more deaths will occur along with the recovery from the disease. Therefore the growth in the recovery class is also different against different immigration. The concerned dynamics are presented by Figs. 3 and 4, respectively.
Further we present by graphs in Figs. 5-6 the dynamical behavior of the transmission dynamics corresponding to the second set of values of immigration rate assumed as [0.0099, 0.0064, 0.0042, 0.0011] as Case II.
We see that the immigration slightly reduced during the thirty days from 10 March to 10 April. The decline in the susceptible population at different rate is shown by Fig. 5, while the corresponding dynamics of the infectious and recovered classes are presented via Figs. 6 and 7, respectively. As the immigration rate is decreasing, the susceptibility is decreasing at rapid speed, and consequently the infection rate is going down. The recov- ered population is also growing with faster speed when immigration rate is low, because the chance of catching infection is decreasing.

Conclusion
A mathematical model addressing the current novel COVID-19 under three compartments, susceptible, infected, and recovered, has been studied. By nonlinear analysis the existence of global and local stability analysis has been demonstrated. On using nonstandard finite difference numerical method, we have simulated the results by using the real data of Wuhan city during the last sixty days from 10 February 2020 to 10 April 2020. Our model has been simulated for the fixed values of the parameters except immigration rate. In the first set of data we have simulated the model against the highest values of immigration rate. We observed that due to this the infection has been rapidly transmitted from person to person during the first thirty days in the mentioned place. During this time more deaths occurred and the recorded rate also increased accordingly. After reducing the concerned immigration rate properly, the dynamics was greatly affected and the infection rate started decreasing, and the recovery rate also increased with different rate because the rate of immigration was different. With smaller immigration rate, the rate of spread of infection is slow as compared at higher order and vice versa. We concluded that minimizing the immigration during this outbreak can cause the increase in protection rate. In other words avoiding unnecessary immigration of people will greatly help in reducing or controlling this disease. Therefore this model is an indication for further study in this area.