key: cord-0989545-zej34z06 authors: Tong, Zhao-Wei; Lv, Yu-Pei; Din, Rahim Ud; Mahariq, Ibrahim; Rahmat, Gul title: Global Transmission Dynamic of SIR Model in The Time of SARS-CoV-2 date: 2021-05-13 journal: Results Phys DOI: 10.1016/j.rinp.2021.104253 sha: 19d05a52acc67ebb95550bdaebc54dc3448767b7 doc_id: 989545 cord_uid: zej34z06 This current work studies a new mathematical model for SARS-CoV-2. We show how immigration, protection, death rate, exposure, cure rate and interaction of infected people with healthy people affect the population. Our model is SIR model, which has three classes including susceptible, infected and recovered respectively. Here, we find the basic reproduction number and local stability through jacobean matrix. Lyapunvo function theory is used to calculate the global stability for the problem under investigation. Also a nonstandard finite difference sachem (NSFDS) is used to simulate the results. It was reported in December 2019, that a new kind of virus named Corona has affected badly the Wuhan city in China. The concerned disease was later named COVID-19 and announced a pandemic throughout the globe by WHO (see [1, 5, 6, 7] ). Firstly, the said virus and it's resultant outbreak hit the Wuhan city and later it affects almost the whole world. It took the lives of hundreds of thousands lives throughout the world. It is hard to find a single view about the origin of the origin of the said virus, for example, it is due to seafood market; migration of people from people one place to other places by transmission from animals to human or may be it is due to human to human interaction. Currently it has almost devastated everything around the world. Social life, health, economy, educational almost every field of the life have been affected badly. Researchers of the field of health, policy makers of the countries and health field are puzzled how to tackle with this deadly outbreak. They all have their own point of views observing the situations. They are trying hard to at least minimize the number of deaths due to this outbreak. People infected with this pandemic experiences a mild respiratory problems. Fever, dry cough, throat infection, and tiredness are the symptoms of this disease. People may have these symptoms; nasal infection, aches, and sore throat. Investigation of dynamics of real world phenomenon through mathematical models is well known area. Particularly, the transmission of diseases which are very important to understand. A number of mathematical models of infectious diseases have been analyzed in the field of mathematics we refer few as [8, 10, 11, 14] . The area of modeling can recently be widened to non-integer order and nonlocal derivatives of fractional order see [14, 15] . Mathematical models are used to understand the Dynamics of transmission of a specific diseases and it can be helpful and to present a palm to save people from lost. Following this technique and with the help of mentioned tools will help to make a strategy to have control on the spread of disease or to eliminate it thoroughly. The recently discovered virus and its resultant outbreak has been studied from different point of views from last few months (see [14, 11, 16, 2] ). Inspired from the above discussion, it is the migration which has played a major role in the spreading of the virus in the societies. It took a short period of almost there months that the said virus spread in almost every part of the world. In this section of manuscript , we formulate our new model for (SARS-CoV-2) in the form (1). We split the whole population N(t) into three classes Susceptible S (t), Infected I(t) and Recovered R(t) class to form system (1) as The dynamic of system (1) is presented in the following flow chart: The physical interpretation for the model (1) are present in Table 1 . The sum of all equations of model (1), we get Where represent From (2), we get With lim t→∞ supN(t) = N 0 if and only if lim t→∞ sup I(t) = 0. The 1 st equation of the system (1) shows We consider the existence of equilibrium for the system (1). Disease free equilibrium exists for some values of the variables in (1) denoted by E 0 = (S 0 , 0, 0). The reproductive number R 0 provides the transmission dynamics information about the disease. From R 0 , we look how the disease id spread in population and we can control it from this. The method of finding R 0 is below. Let X(t) = (S , I), then form system (1), one has We get We get the expression for R 0 is To computes the basic reproduction number we obtained R 0 = 3.13 from the parameters used in Table 2 . This show that the COVID-19 that occurred in Pakistan is not well controlled by govt of Pakistan. We have the following theorem on the basis of (4). (1); (1) called the endemic equilibrium. We reduced our model (1) to obtained the result for local stability as Subject to initial condition Here we established the coming result. Proof. At E 0 , we find Jacobian Matrix as The auxiliary equation of J 0 is given by We have (6) show that the "Routh-Hurtwriz criteria is satisfied" as k 1 > 0, k 2 > 0, k 3 > 0 and which show the system (1) is locally asymptotically stable at E 0 . (1) is locally asymptotically stable, otherwise unstable. Proof. The jacobian matrix of the model (5), at the endemic equilibrium E * is After some operations on matrix J 1 , we get We calculated det(J 1 ) and trace of J 1 is given below traJ 1 = −(bI * (t) + 2δ + µ + e) < 0. and det J 1 = (δ + bI * (t))(µ + δ + e) + δbS * (t). > 0. The determinant of J 1 > 0.The "endemic equilibrium" of the system (5) at E * (t) has negative real part. Thus, with condition R 0 > 1, one gets that the endemic equilibrium E * of system (5) is locally asymptotically stable. In this section of our work, we study global stability of model (1) . For global stability, we constructed "Lyapunov function" for diseases free and endemic equilibria. Theorem 5.1. The disease free equilibrium of the system (1) is globally asymptotically stable if R 0 < 1. Proof. To prove this required result, we construct a Lyapunov function as following : such that α 1 , α 2 > 0 are constant may be computed later. With respect to time t taking derivative of (7) with , we have dΨ dt = α 1 (β + a − δS (t) − bS (t)I(t)) + α 2 (bS (t)I(t) + c − (δ + µ + e)I(t)). Let assume α 1 = α 2 = 1, we get finally Hence "globally asymptotically stable" for system (1) with R 0 < 1 has reached. Proof. On construction of Lyapunov function as Take the derivative (8) w.r.t t along with system (5), Putting the values from (1) yields After some rearrangement, one has Thus dΘ dt < 0, "The endemic equilibrium" E * (t) of the model (1) is " Globally asymptotically stable", which show that R 0 > 1. In this part of our manuscript, we study simulation for model (1) by using the values of Table 2 . From first January upto next 300 days, we simulate the considered model (1) for Pakistan a. Here, using NSFD [16, 18] to rewrite the system in difference equations form. Hence from the first equation Thank to NSFD, (9) yields Like (10), Using NSFD scheme, we decomposed the system (1) as follow By using NSFD scheme, developed in (11) . Using the values from 2, we plotted the system (1). Using real data , we testify our system (1) for the real data of Pakistan with used in Pakistan at different values of a, c. During the 1st January 2021 in the present of immigration rate the susceptible population dynamic is given in Figure 1 . When the rate of immigration is high the decline in the host population of susceptible is decreasing. This yields that they are exposing to infection and hence higher the immigration rate faster the growth rate of infected papulation and vice versa. From the recovery SARS-CoV-2 more death will occurs. In the recovery compartment the growth is different under various immigrant rates. The concerned dynamical behavior is provided Figure 2 and Figure 3 . The recovered papulation is also growing with faster speed when immigration rate is low. This is due to the people getting ride from infection. As the recovery rate is high in Pakistan as majority of population is consisted on young people between 20 and 45 years. Hence bringing change in certain parameters of the considered model, significant change can be observed in the dynamics of respective compartments. Further the adopted numerical scheme is good to use to simulate mathematical model. To overcome the pandemic the migration should be decreased for the sake of saving humanity. Also the immigration of exposed population to infected community increase the infection. Isolation of infected one is the best option to secured the healthy community. It is necessary to judge the spread and model with various parameters for proper supervision. The proper treatment of this pandemic is to keep infected away from healthy people. High internal defense system aids to get healthy soon while the low internal defense system need more attention. This is model is just an indication to see the interaction of different compartments and its dynamics by varing the values of certain parameters. Since precautionary measures are not strict in well populated country like Pakistan which cause faster growth of infection in the second wave of COVID-19. " Competing interest In this work no competing interest exists. Funding This work has been supported by the Huzhou city emergency science and technology plan project Availability of Data This is not applicable in this research work. All authors played their role equally." Naming the coronavirus disease (COVID-19) and the virus that causes it The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health-The latest 2019 novel coronavirus outbreak in Wuhan, China Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China Estimating the unreported number of novel coronavirus (2019-nCoV) cases in China in the first half of January 2020, a data-driven Modelling analysis of the early outbreak Statistics based predictions of coronavirus 2019-nCoV spreading in mainland China Breaking down of healthcare system: Mathematical modelling for controlling the novel coronavirus (2019-nCoV) outbreak in Wuhan Contribution to the theory of periodic reactions On the Volterra and other nonlinear models of interacting populations A pneumonia outbreak associated with a new coronavirus of probable bat origin Early transmission dynamics in Wuhan, China, of novel coronavirus infected pneumonia Periodicity in epidemiological models,Pneumonia of Unknown Etiology in Wuhan, China: Potential for International Spread Via Commercial Air Travel Pakistan Covid-19 Health Advisory platform by minnistry Modelling strategies for controlling SARS out breaks Coronavirus (COVID-19) Mortality Rate Now casting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China, a modelling study Host vector dynamics of pine wilt disease model with convex incidence rate Introduction to non-standard finite-difference modelling On nonlinear classical and fractional order dynamical system addressing COVID-19 Study of global dynamics of COVID-19 via a new mathematical model Mohammed Al-Shomrani, and Abioye Abioye Idowu. A new mathematical model of COVID-19 using real data from Pakistan Nonlinear equations with global differential and integral operators: existence, uniqueness with application to epidemiology Facemasks simple but powerful weapons to protect against COVID-19 spread: Can they have sides effects