Skip to main content

Theory and Modern Applications

New investigation of bats-hosts-reservoir-people coronavirus model and application to 2019-nCoV system

Abstract

According to the report presented by the World Health Organization, a new member of viruses, namely, coronavirus, shortly 2019-nCoV, which arised in Wuhan, China, on January 7, 2020, has been introduced to the literature. The main aim of this paper is investigating and finding the optimal values for better understanding the mathematical model of the transfer of 2019-nCoV from the reservoir to people. This model, named Bats-Hosts-Reservoir-People coronavirus (BHRPC) model, is based on bats as essential animal beings. By using a powerful numerical method we obtain simulations of its spreading under suitably chosen parameters. Whereas the obtained results show the effectiveness of the theoretical method considered for the governing system, the results also present much light on the dynamic behavior of the Bats-Hosts-Reservoir-People transmission network coronavirus model.

1 Introduction

Today the whole world has been witnessing and suffering from a big pandemic disease, the novel coronavirus pneumonia, named as “2019 novel coronavirus (2019-nCoV)” [1]. The first wave of outbreak happened in Wuhan, China, and then quickly spread into some other parts of China and even to other countries due to its high person-to-person infection rate. The evidence of human-to-human transmission of it was identified by the National Health Commission (NHC) of the People’s Republic of China [2]. Based on statistics, there are more than 80,000 infected confirmed patients in China and around 7000 reported cases out of China (by the date March 16, 2020). Until now, around 150 countries reported the confirmed infection cases. In other words, it has been a global serious infective disease. Several researchers compared 2019-nCoV with SARS (severe acute respiratory syndrome) epidemic, which was out-broken in China in 2003 and pointed that the basic reproduction number (R0) of 2019-nCoV is estimated to be even higher than SARS with its reproduction number 6.47 vs SARS’s 4.91, which indicates a high spread and infection of 2019-nCoV and also how severe and acute the disaster we are confronted.

On December 2019, the symptom of 2019-nCoV infected patients was identified as fever, cough, breathing difficulties, and some other. Due to the long incubation period and mild symptoms, the suspicious infected people need to be observed for around 14 days. To reduce population flow and restrict the spread of the virus, corresponding virus dissemination control policies and relevant actions are being carried out at different levels. On January 23, Wuhan was locked down by strict restrictions of transport, and soon some other provinces announced to lockdown as well [3]. Chinese people were suggested to stay at home and avoid gathering, assembling, celebrations, visiting, and so on to reduce the virus dissemination. People were required to wear respirators in public areas. Some researchers pointed out some effective ways to control the spread of infectious virus, including school closure, case isolation, household quarantine, internal travel restrictions, and border control, which were proved to be helpful for the delay or reduction of virus infections [46]. Almost all the economic activities in some countries have been paused, which caused countless damages to humans’ lives, development, and also a large amount of financial pressure on these countries. Globally, Australia and New Zealand firstly released regulations to ban travellers who had been to China in the past 14 days. From December 2019 up to now, 2019-nCoV virus suddenly out-broke in other countries in flood. Italy, Japan, South Korea, Spain, and other countries were in the clouds of the virus as well, and some of them started to lockdown cities, restrict transports, close schools, and so on. The whole world is in damage by the 2019-nCoV and in the combat against coronavirus. Some other researchers devoted themselves to the drug and vaccine developments, however, there is no way to effectively eliminate the virus in the human’s body.

A large number of researchers have studied the 2019-nCoV from a wide range of perspectives, including diverse infectious diseases, microbiology, virology, respiratory system, biochemistry molecular biology, immunology, public environmental occupational health, genetics heredity, veterinary sciences, environmental sciences ecology, and pathology. Most of them are conducted by the USA and China, followed by Saudi Arabia, South Korea, and Germany. According to Tian [7] and others, the differences of RBD between SARS-CoV and 2019-nCoV is of significance for the cross-reactivity of neutralizing antibodies, and a SARS-CoV-specific human monoclonal antibody CR3022 could bind potently with 2019-nCoV RBD (KD of 6.3 nM). As for the origins of 2019-nCoV virus, 2019-nCoV, Benvenuto et al. [8] held that 2019nCoV could be considered a coronavirus distinct from SARS virus, probably transmitted from bats or another host, where mutations conferred upon it the ability to infect humans, and they also proposed a preliminary evolutionary and molecular epidemiological analysis of this new virus, considering high genetic similarity between 2019-nCoV and the Severe Acute Respiratory Syndrome coronavirus (SARS-CoV), and leveraging existing immunological studies of SARS-CoV, Ahmed, Quadeer, and McKay [9] devoted themselves to seeking for gaining insights for vaccine design against 2019-nCoV.

With the increasing of virus spread and the ongoing related research both at domestic and international, one question still hinders human’s knowledge of the 2019-nCoV: What is the original source of such a virus and how can it transmit to human. In this paper, we intend to study a mathematical model called the Bats-Hosts-Reservoir-People coronavirus (BHRPC) model for the transfer of 2019-nCoV from the reservoir to people. By using a powerful numerical method we gain its spreading simulations under suitably chosen parameters. The obtained results show the effectiveness of the theoretical method considering the governing system and also present much light on the dynamic behavior of the Bats-Hosts-Reservoir-People transmission network coronavirus model. In this regards, more recently, some experts have investigated some important nonlinear models arising in real-world problems [1016].

One of such problems has been mathematically developed by Chen et al. for simulating the phase-based transmissibility of a novel coronavirus as 2019-nCoV [17] defined by

$$ \begin{gathered} \frac{du ( t )}{dt} = n_{p} - m_{p}u ( t ) - b_{p}u ( t ) \bigl[ y ( t ) + \kappa f ( t ) \bigr] - b_{w}u ( t )w ( t ), \\ \frac{dv ( t )}{dt} = b_{p}u ( t ) \bigl[ y ( t ) + \kappa f ( t ) \bigr] + b_{w}u ( t )w ( t ) - ( 1 - \delta_{p} ) \omega_{p}v ( t ) - \delta_{p}\omega '_{p}v ( t ) - m_{p}v ( t ), \\ \frac{dy ( t )}{dt} = ( 1 - \delta_{p} )\omega_{p}v ( t ) - ( \gamma_{p} + m_{p} )y ( t ), \\ \frac{df ( t )}{dt} = \delta_{p}\omega '_{p}v ( t ) - \bigl[ \gamma '_{p} + m_{p} \bigr]f ( t ), \\ \frac{dr ( t )}{dt} = \gamma_{p}y ( t ) + \gamma '_{p}f ( t ) - m_{p}r ( t ), \\ \frac{dw ( t )}{dt} = \varepsilon \bigl[ y ( t ) + cf ( t ) - w ( t ) \bigr], \end{gathered} $$
(1)

where \(n_{p}\), \(m_{p}\), \(b_{p}\), κ, \(b_{w}\), \(\delta_{p}\), ε, c are real nonzero constants. The initial conditions for this system are shortly given by

$$ \begin{gathered} u(0)=\beta_{1},\qquad v(0)=\beta_{2},\qquad y(0)=\beta_{3},\qquad f(0)= \beta_{4},\qquad r(0)=\beta_{5}, \\ \omega(0)=\beta_{6},\qquad \omega_{p}'(0)=\beta_{7}, \qquad w(0)=\beta_{10},\qquad \gamma_{p}'(0)=\beta_{8},\qquad \gamma(0)= \beta_{9}.\end{gathered} $$
(2)

Equation (1) is used to describe the phase-based transmissibility of a novel coronavirus from source to people. In Eq. (1), u is the susceptible people, v is used to symbolize exposed people, y is the symptomatic infected people, f is asymptomatic infected people, r is removed people (recovered and died people), \(n_{p}\) is the birth rate, \(m_{p}\) is the death rate of people, w is the reservoir (the seafood area), \(1 / \omega_{B}\) is the incubation period of bat infection, and \(1 / \gamma_{B}\) is the infectious period of bat infection [17]. Khan et al. [10] have investigated the endemic equilibria, stability, and global sensitivity of system (1). If \(r_{0} < 1\), then system (1) is locally asymptotically stable [10, 18], and the outbreak will fade away [18]. When \(r_{0} > 1\), the outbreak will occur [18], and it is not stable. The used data for system (1) are for Wuhan, China, [17, 19]. Moreover, mathematical analysis and applications of dengue fever outbreak and epidemiology in the sense of fractional have been investigated [20, 21]. Recently, a numerical scheme based on the Newton polynomial has been applied successfully to observe important properties of the spread of COVID-19 with new fractal-fractional operators [22]. Some close relationships of COVID-19 with HIV have been presented in [23]. Many applications of fractional- or integer-order mathematical models explaining more detailed informations about the real-world problems have been presented in a detailed manner [2450]. In this paper, we investigate the numerical distributions of 2019-nCoV according to time with the help of several approaching terms of VIM.

2 Some important properties of VIM

In 1999, VIM, one of the most powerful numerical methods, was firstly developed by He [5154] for numerical investigation and exceeds the difficulties of the perturbation or Adomian functions. Later, Wazwaz has applied VIM for investigating linear and nonlinear wave equations along with wave-like equations [55] and Laplace equation [56]. Moreover, many applications of VIM have been observed for various models [1417, 5160]. We consider differential equations of the form

$$ Lu + Nu = f ( x ), $$
(3)

where L and N are linear and nonlinear operators, respectively, [55], and \(f ( x )\) is a source inhomogeneous term. According to basic concepts of VIM presented by He, we construct the following iteration formula for Eq. (3) [5154]:

$$ u_{n + 1} ( x ) = u_{n} ( x ) + \int_{0}^{t} \lambda ( \tau ) \bigl[ L \widetilde{u}_{n} ( \tau ) + R\widetilde{u}_{n} ( \tau ) + N\widetilde{u}_{n} ( \tau ) - f ( \tau ) \bigr]\,\mathrm{d} \tau, $$
(4)

where the parameter λ is a general Lagrange multiplier, which can be optimally identified via the variational theory, the subscript n denotes the nth-order approximation, and \(\widetilde{u}_{n}\) is considered as a restricted variation, which means \(\delta \widetilde{u}_{n} = 0\). Clearly, the main steps of VIM first require the determination of the Lagrange multiplier λ, which needs to be optimally identified. Once λ is determined, the successive approximations \(u_{n + 1}\), \(n \ge 0\), of the solution u are obtained upon using a suitably selected function \(u_{0}\), which satisfies the boundary conditions. Then the solution is given by

$$ u = \lim_{n \to \infty} u_{n}. $$
(5)

3 Application of VIM to COVID-19 system

In this subsection, by using VIM we numerically investigate the Bats-Hosts-Reservoir-People coronavirus model. According to VIM iteration structure, we can write Eq. (1) in the following form:

$$ \begin{gathered} u_{k + 1} = u_{k} + \int_{0}^{t} \lambda \biggl[ \frac{du_{k}}{d\tau} - n_{p} + m_{p}u_{k} + b_{p}u_{k} [ y_{k} + \kappa f_{k} ] + b_{w}u_{k}w_{k} \biggr] \,d\tau, \\ v_{k + 1} = v_{k} + \int_{0}^{t} \lambda \biggl[ \frac{dv_{k}}{d\tau} - b_{p}u_{k} [ y_{k} + \kappa f_{k} ] - b_{w}u_{k}w_{k} \\ \phantom{v_{k + 1} =}{} + ( 1 - \delta_{p} )\omega_{p}v_{k} + \delta_{p}\omega '_{p}v_{k} - m_{p}v_{k} \biggr] \,d\tau, \\ y_{k + 1} = y_{k} + \int_{0}^{t} \lambda \biggl[ \frac{dy_{k}}{d\tau} - ( 1 - \delta_{p} )\omega_{p}v_{k} + ( \gamma_{p} + m_{p} )y_{k} \biggr] \,d \tau, \\ f_{k + 1} = f_{k} + \int_{0}^{t} \lambda \biggl[ \frac{df_{k}}{d\tau} - \delta_{p}\omega '_{p}v_{k} + \bigl[ \gamma '_{p} + m_{p} \bigr]f_{k} \biggr]\, d\tau, \\ r_{k + 1} = r_{k} + \int_{0}^{t} \lambda \biggl[ \frac{dr_{k}}{d\tau} - \gamma_{p}y_{k} - \gamma '_{p}f_{k} + m_{p}r_{k} \biggr] \,d\tau, \\ w_{k + 1} = w_{k} + \int_{0}^{t} \lambda \biggl[ \frac{dw_{k}}{d\tau} - \varepsilon y_{k} - \varepsilon cf_{k} + \varepsilon w_{k} \biggr]\, d\tau, \end{gathered} $$
(6)

where \(k = 0,1,2,3, \dots\). It produces the stationary condition

$$ \lambda = - 1. $$
(7)

Substituting Eq. (7) into Eq. (6), we find the following iteration equation:

$$ \begin{gathered} u_{k + 1} = u_{k} - \int_{0}^{t} \biggl[ \frac{du_{k}}{d\tau} - n_{p} + m_{p}u_{k} + b_{p}u_{k} [ y_{k} + \kappa f_{k} ] + b_{w}u_{k}w_{k} \biggr]\, d\tau, \\ v_{k + 1} = v_{k} - \int_{0}^{t} \biggl[ \frac{dv_{k}}{d\tau} - b_{p}u_{k} [ y_{k} + \kappa f_{k} ] - b_{w}u_{k}w_{k} \\ \phantom{v_{k + 1} =}{} + ( 1 - \delta_{p} )\omega_{p}v_{k} + \delta_{p}\omega '_{p}v_{k} - m_{p}v_{k} \biggr] \,d\tau, \\ y_{k + 1} = y_{k} - \int_{0}^{t} \biggl[ \frac{dy_{k}}{d\tau} - ( 1 - \delta_{p} )\omega_{p}v_{k} + ( \gamma_{p} + m_{p} )y_{k} \biggr] \,d\tau, \\ f_{k + 1} = f_{k} - \int_{0}^{t} \biggl[ \frac{df_{k}}{d\tau} - \delta_{p}\omega '_{p}v_{k} + \bigl[ \gamma '_{p} + m_{p} \bigr]f_{k} \biggr] \,d\tau, \\ r_{k + 1} = r_{k} - \int_{0}^{t} \biggl[ \frac{dr_{k}}{d\tau} - \gamma_{p}y_{k} - \gamma '_{p}f_{k} + m_{p}r_{k} \biggr] \,d\tau, \\ w_{k + 1} = w_{k} - \int_{0}^{t} \biggl[ \frac{dw_{k}}{d\tau} - \varepsilon y_{k} - \varepsilon cf_{k} + \varepsilon w_{k} \biggr] \,d\tau. \end{gathered} $$
(8)

With the help of some computational software algorithm, by considering Eq. (2) with initial values we get the first approximate components of the system for \(k=0\):

$$\begin{gathered} u_{1} = u_{0} - \int_{0}^{t} \biggl[ \frac{du_{0}}{d\tau} - n_{p} + m_{p}u_{0} + b_{p}u_{0} [ y_{0} + \kappa f_{0} ] + b_{w}u_{0}w_{0} \biggr] \,d\tau, \\ v_{1} = v_{0} - \int_{0}^{t} \biggl[ \frac{dv_{0}}{d\tau} - b_{p}u_{0} [ y_{0} + \kappa f_{0} ] - b_{w}u_{0}w_{0} \\ \phantom{v_{1} =}{}+ ( 1 - \delta_{p} )\omega_{p} ( 0 )v_{0} + \delta_{p}\omega '_{p} ( 0 )v_{0} - m_{p}v_{0} \biggr] \,d\tau, \\ y_{1} = y_{0} - \int_{0}^{t} \biggl[ \frac{dy_{0}}{d\tau} - ( 1 - \delta_{p} )\omega_{p} ( 0 )v_{0} + \bigl( \gamma_{p} ( 0 ) + m_{p} \bigr)y_{0} \biggr] \,d\tau, \\ f_{1} = f_{0} - \int_{0}^{t} \biggl[ \frac{df_{0}}{d\tau} - \delta_{p}\omega '_{p} ( 0 )v_{0} + \bigl[ \gamma '_{p} ( 0 ) + m_{p} \bigr]f_{0} \biggr] \,d\tau, \\ r_{1} = r_{0} - \int_{0}^{t} \biggl[ \frac{dr_{0}}{d\tau} - \gamma_{p} ( 0 )y_{0} - \gamma '_{p} ( 0 )f_{0} + m_{p}r_{0} \biggr] \,d\tau, \\ w_{1} = w_{0} - \int_{0}^{t} \biggl[ \frac{dw_{0}}{d\tau} - \varepsilon y_{0} - \varepsilon cf_{0} + \varepsilon w_{0} \biggr]\, d\tau. \end{gathered} $$

Thus we get the first approach of \(u_{n}\), \(v_{n}\), \(y_{n}\), \(f_{n}\), \(r_{n}\), \(w_{n}\) as follows:

$$\begin{aligned}& u_{1} ( t ) = \beta_{1} + \tau_{1}t, \\& v_{1} ( t ) = \beta_{2} + \tau_{2}t, \\& y_{1} ( t ) = \beta_{3} + \tau_{3}t, \\& f_{1} ( t ) = \beta_{4} + \tau_{4}t, \\& r_{1} ( t ) = \beta_{5} + \tau_{5}t, \\& w_{1} ( t ) = \beta_{10} + \tau_{6}t, \end{aligned}$$
(9)

where, for simplicity, we have taken

$$\begin{gathered} \tau_{1} = n_{p} - m_{p} \beta_{1} - b_{p}\beta_{1} \beta_{3} - \kappa b_{p}\beta_{1} \beta_{4} - b_{w}\beta_{1} \beta_{10}, \\ \tau_{2} = b_{p}\beta_{1} [ \beta_{3} + \kappa \beta_{4} ] + b_{w} \beta_{1}\beta_{10} - ( 1 - \delta_{p} ) \beta_{6}\beta_{2} - \delta_{p} \beta_{7}\beta_{2} + m_{p} \beta_{2}, \\ \tau_{3} = ( 1 - \delta_{p} )\beta_{6} \beta_{2} - ( \gamma_{0} + m_{p} ) \beta_{3}, \\ \tau_{4} = \delta_{p}\beta_{7} \beta_{2} - [ \beta_{8} + m_{p} ] \beta_{4}, \\ \tau_{5} = \gamma_{0}\beta_{3} + \gamma '_{k}\beta_{4} - m_{p}r_{0}, \\ \tau_{6} = \varepsilon \beta_{3} + \varepsilon c \beta_{4} - \varepsilon \beta_{10}. \end{gathered} $$

Now we obtain the second components of the variables \(u_{n}\), \(v_{n}\), \(y_{n}\), \(f_{n}\), \(r_{n}\), \(w_{n}\) for \(k=1\):

$$\begin{gathered} u_{2} = u_{1} - \int_{0}^{t} \biggl[ \frac{du_{1}}{d\tau} - n_{p} + m_{p}u_{1} + b_{p}u_{1} [ y_{1} + \kappa f_{1} ] + b_{w}u_{1}w_{1} \biggr]\, d\tau, \\ v_{2} = v_{1} - \int_{0}^{t} \biggl[ \frac{dv_{1}}{d\tau} - b_{p}u_{1} [ y_{1} + \kappa f_{1} ] - b_{w}u_{1}w_{1}+ ( 1 - \delta_{p} )\omega_{p}v_{1} + \delta_{p}\omega '_{p}v_{1} - m_{p}v_{1} \biggr] \,d\tau, \\ y_{2} = y_{1} - \int_{0}^{t} \biggl[ \frac{dy_{1}}{d\tau} - ( 1 - \delta_{p} )\omega_{p}v_{1} + ( \gamma_{p} + m_{p} )y_{1} \biggr]\, d\tau, \\ f_{2} = f_{1} - \int_{0}^{t} \biggl[ \frac{df_{1}}{d\tau} - \delta_{p}\omega '_{p}v_{1} + \bigl[ \gamma '_{p} + m_{p} \bigr]f_{1} \biggr] \,d\tau, \\ r_{2} = r_{1} - \int_{0}^{t} \biggl[ \frac{dr_{1}}{d\tau} - \gamma_{p}y_{1} - \gamma '_{p}f_{1} + m_{p}r_{1} \biggr] \,d\tau, \\ w_{2} = w_{1} - \int_{0}^{t} \biggl[ \frac{dw_{1}}{d\tau} - \varepsilon y_{1} - \varepsilon cf_{1} + \varepsilon w_{1} \biggr]\, d\tau, \end{gathered} $$

which gives

$$ \begin{gathered} u_{2} ( t ) = \beta_{1} + ( \tau_{1} + \tau_{7} )t + \frac{\tau_{8}}{2}t^{2} + \frac{\tau_{9}}{3}t^{3}, \\ v_{2} ( t ) = \beta_{2} + ( \tau_{2} + \tau_{10} )t + \frac{\tau_{11}}{2}t^{2} + \frac{\tau_{12}}{3}t^{3}, \\ y_{2} ( t ) = \beta_{3} + ( \tau_{3} + \tau_{13} )t + \frac{\tau_{14}}{2}t^{2}, \\ f_{2} ( t ) = \beta_{4} + ( \tau_{4} + \tau_{15} )t + \frac{\tau_{16}}{2}t^{2}, \\ r_{2} ( t ) = \beta_{5} + ( \tau_{5} + \tau_{17} )t + \frac{\tau_{18}}{2}t^{2}, \\ w_{2} ( t ) = \beta_{10} + ( \tau_{6} + \tau_{19} )t + \frac{\tau_{20}}{2}t^{2}, \end{gathered} $$
(10)

where

$$\begin{gathered} \tau_{7} = \tau + n_{p} - m_{p}\beta_{1} - b_{p}\beta_{1} \beta_{3} - b_{p}\beta_{1}\kappa \beta_{4} - \beta_{1}\beta_{10}b_{w}, \\ \tau_{8} = m_{p}\tau_{1} + b_{p}\tau_{1}\beta_{3} + \tau_{3}b_{p}\beta_{1} + \tau_{6}b_{w}\beta_{1} + b_{p} \tau_{1}\kappa \beta_{4} + \tau_{4}b_{p} \beta_{1}\kappa + \tau_{1}\beta_{10}b_{w}, \\ \tau_{9} = - \tau_{6}b_{w} \tau_{1}t^{2} - \tau_{3}b_{p} \tau_{1}t^{2} - \tau_{4}b_{p} \tau_{1}\kappa, \\ \tau_{10} = \tau_{2} - \beta_{2} ( \omega_{p} - \omega_{p}\delta_{p} ) + b_{p}\beta_{1}\beta_{3} + b_{p} \kappa \beta_{1}\beta_{4} + b_{w} \beta_{10}\beta_{1} - \beta_{2} \delta_{p}\omega '_{p} + m_{p} \beta_{2}, \\ \tau_{11} = - b_{p}\tau_{1} \beta_{3} - \tau_{3}b_{p}\beta_{1} - b_{p}\kappa \tau_{1}\beta_{4} - \tau_{4}b_{p}\kappa \beta_{1} - \beta_{10}b_{w}\tau_{1} - \tau_{6}b_{w}\beta_{1} \\ \phantom{\tau_{11} =}{}+ \tau_{2} ( \omega_{p} - \omega_{p}\delta_{p} ) + \tau_{2} \delta_{p}\omega '_{p} - m_{p} \tau_{2}, \\ \tau_{12} = \tau_{3}b_{p} \tau_{1} + \tau_{4}b_{p}\kappa \tau_{1} + \tau_{6}b_{w} \tau_{1}, \\ \tau_{13} = \tau_{3} - \beta_{2} ( - 1 + \delta_{p} )\omega_{p} - \beta_{3} ( \gamma_{p} + m_{p} ), \\ \tau_{14} = + \tau_{2} ( - 1 + \delta_{p} ) \omega_{p} + \tau_{3} ( \gamma_{p} + m_{p} ), \\ \tau_{15} = \tau_{4} + \beta_{2} \delta_{p}\omega '_{p} - \beta_{4} \bigl( \gamma '_{p} + m_{p} \bigr), \\ \tau_{16} = \tau_{4} \bigl( \gamma '_{p} + m_{p} \bigr) - \tau_{2}\delta_{p}\omega '_{p}, \\ \tau_{17} = \tau_{5} + \beta_{3} \gamma_{p} + \gamma '_{p} \beta_{4} - \beta_{5}m_{p}, \\ \tau_{18} = - \gamma '_{p} \tau_{4} + \tau_{5}m_{p} - \tau_{3}\gamma_{p}, \\ \tau_{19} = \tau_{6} + \varepsilon \beta_{3} + \varepsilon c\beta_{4} - \varepsilon \beta_{10}, \\ \tau_{20} = - \varepsilon c\tau_{4} - \varepsilon \tau_{3} + \varepsilon \tau_{6}. \end{gathered} $$

The remaining components of the iteration formula (5) can be found in the same manner using a similar algorithm via various computational schemes. In this work, we observe the spreading rate of the Bats-Hosts-Reservoir-People coronavirus model from reservoir to people under suitably chosen values of the parameters, reported by experts in Wuhan area of China. We obtain two-dimensional simulations of the second terms \(u_{2}\), \(v_{2}\), \(y_{2}\), \(f_{2}\), \(r_{2}\), \(w_{2}\) of \(u_{n}\), \(v_{n}\), \(y_{n}\), \(f_{n}\), \(r_{n}\), \(w_{n}\) as in Fig. 1.

Figure 1
figure 1

Two-dimensional surfaces of system (10)

4 Conclusions

In this paper, we have successfully applied VIM to numerical investigation of the 2019-nCoV model. This method is based on a series solution terms of iteration. Only in the second terms of iteration, we have obtained numerical results for system (10). Under suitably chosen values of the parameters, reported by WHO, we have plotted the numerical results. According to Figs. 1234, we observed that VIM produces similar distributions for susceptible people, which increase exponentially, when we compare the simulated data with Fig. 2 in [17]. Moreover, we can say that each class of system (1) has also simulated estimated behaviors from them. Furthermore, the spread of the 2019-nCoV with susceptible people \(u_{2}\) is faster than the others, such as \(v_{2}\), \(y_{2}\), \(f_{2}\), \(r_{2}\), \(w_{2}\). Finally, from the second terms of the proposed algorithm we can observe that susceptible people will affect more and more people from all over the world. As a future direction of this concept, application of powerful projected tools may be studied. These produce more comprehensive results on the mathematical system of 2019-nCoV.

Figure 2
figure 2

Two-dimensional surfaces of susceptible people and exposed people of system (10)

Figure 3
figure 3

Two-dimensional surfaces of symptomatic infected people and asymptomatic infected people of system (10)

Figure 4
figure 4

Two-dimensional surfaces of removed people (recovered or died people) and reservoir (the seafood area) of system (10)

References

  1. Jasper, F.-W.C., Kok, K.-H., Zheng, Z., et al.: Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from patients with acute respiratory disease in Wuhan, Hubei, China. Emerg. Microbes Infect. (2020)

  2. National Health Commission of the People’s Republic of China. Top expert: disease spread won’t be on scale of SARS 2020 [updated January 21; cited March 16, 2020]. http://en.nhc.gov.cn/2020-01/21/c_75991.htm

  3. South China Morning Post. China coronavirus: three cities join Wuhan in quarantine lockdown as Beijing tries to contain deadly outbreak 2020 [updated January 23; cited March 16, 2020]. https://www.scmp.com/news/china/society/article/3047278/wuhan-goes-shutdown-china-tries-contain-deadly-coronavirus

  4. Ferguson, N.M., Cummings, D.A., Fraser, C., Cajka, J.C., Cooley, P.C., Burke, D.S.: Strategies for mitigating an influenza pandemic. Nature 442(7101), 448–452 (2006)

    Google Scholar 

  5. Cowling, B.J., Ho, L.M., Leung, G.M.: Effectiveness of control measures during the SARS epidemic in Beijing: a comparison of the Rt curve and the epidemic curve. Epidemiol. Infect. 136(4), 562–566 (2008)

    Google Scholar 

  6. Koonin, L.M., Cetron, M.S.: School closure to reduce influenza transmission. Emerg. Infect. Dis. 15(1), 137–138 (2009). https://doi.org/10.3201/eid1501.081289

    Article  Google Scholar 

  7. Tian, X., Li, C., Huang, A., Xia, S., Lu, S., Shi, Z., Lu, L., Jiang, S., Yang, Z., Wu, Y., Ying, T.: Potent binding of 2019 novel coronavirus spike protein by a SARS coronavirus-specific human monoclonal antibody. Emerg. Microbes Infect. (2020). https://doi.org/10.1080/22221751.2020.1729069

    Article  Google Scholar 

  8. Benvenuto, D., Giovannetti, M., Ciccozzi, A., Spoto, S., Angeletti, S., Ciccozzi, M.: The 2019-new coronavirus epidemic: evidence for virus evolution. J. Med. Virol. (2020). https://doi.org/10.1002/jmv.25688

    Article  Google Scholar 

  9. Ahmed, S.F., Quadeer, A.A., McKay, M.R.: Preliminary identification of potential vaccine targets for the COVID-19 coronavirus (SARS-CoV-2) based on SARS-CoV immunological studies. Viruses (2020). https://doi.org/10.3390/v12030254

    Article  Google Scholar 

  10. Khan, M.A., Atangana, A.: Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative. Alex. Eng. J. (2020). https://doi.org/10.1016/j.aej.2020.02.033

    Article  Google Scholar 

  11. Is the world ready for the coronavirus? Editorial. The New York Times. 29 January 2020. Archived from the original on 30 January 2020

  12. China virus death toll rises to 41, more than 1,300 infected worldwide. CNBC. 24 January 2020. Archived from the original on 26 January 2020. Retrieved 26 January 2020. Retrieved 30 January 2020

  13. Gao, W., Veeresha, P., Prakasha, D.G., Baskonus, H.M., Yel, G.: New approach for the model describing the deathly disease in pregnant women using Mittag-Leffler function. Chaos Solitons Fractals 134, Article ID 109696 (2020)

    MathSciNet  Google Scholar 

  14. Owolabi, K.M., Atangana, A.: Mathematical analysis and computational experiments for an epidemic system with nonlocal and nonsingular derivative. Chaos Solitons Fractals 126, 41–49 (2019)

    MathSciNet  Google Scholar 

  15. Naik, P.A., Zu, J., Owolabi, K.M.: Modeling the mechanics of viral kinetics under immune control during primary infection of HIV-1 with treatment in fractional order. Phys. A, Stat. Mech. Appl. 13, Article ID 123816 (2019)

    Google Scholar 

  16. Gao, W., Ghanbari, B., Baskonus, H.M.: New numerical simulations for some real world problems with Atangana–Baleanu fractional derivative. Chaos Solitons Fractals 128, 34–43 (2019)

    MathSciNet  Google Scholar 

  17. Chen, T.M., Rui, J., Wang, Q.P., Zhao, Z.Y., Cui, J.A., Yin, L.: A mathematical model for simulating the phase-based transmissibility of a novel coronavirus. Infect. Dis. Poverty 9, Article ID 24 (2020). https://doi.org/10.1186/s40249-020-00640-3

    Article  Google Scholar 

  18. van den Driessche, P., Watmough, J.: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 180, 29–48 (2002)

    MathSciNet  MATH  Google Scholar 

  19. Li, Q., Guan, X., Wu, P., Wang, X., Zhou, L., Tong, Y., et al.: Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N. Engl. J. Med. 382, 1199–1207 (2020)

    Google Scholar 

  20. Qureshi, S., Atangana, A.: Mathematical analysis of dengue fever outbreak by novel fractional operators with field data. Phys. A, Stat. Mech. Appl. 526, Article ID 121127 (2019)

    MathSciNet  Google Scholar 

  21. Atangana, A.: Application of fractional calculus to epidemiology. In: Fractional Dynamics, pp. 174–190 (2015). https://doi.org/10.1515/9783110472097-011

    Chapter  Google Scholar 

  22. Atangana, A.: Modelling the spread of COVID-19 with new fractal-fractional operators: can the lockdown save mankind before vaccination? Chaos Solitons Fractals 136, Article ID 109860 (2020)

    MathSciNet  Google Scholar 

  23. Doungmo Goufo, E.F., Khan, Y., Chaudhry, Q.A.: HIV and shifting epicenters for COVID-19, an alert for some countries. Chaos Solitons Fractals 139, Article ID 110030 (2020)

    MathSciNet  Google Scholar 

  24. Ihan, E., Kıymaz, I.O.: A generalization of truncated M-fractional derivative and applications to fractional differential equations. Appl. Math. Nonlinear Sci. 5(1), 171–188 (2020)

    MathSciNet  Google Scholar 

  25. Baskonus, H.M.: Complex surfaces to the fractional \((2+1)\)-dimensional Boussinesq dynamical model with local M-derivative. Eur. Phys. J. Plus 134, Article ID 322 (2019)

    Google Scholar 

  26. Cattani, C.: A review on harmonic wavelets and their fractional extension. J. Adv. Eng. Comput. 2(4), 224–238 (2018)

    Google Scholar 

  27. Zhang, Y., Cattani, C., Yang, X.J.: Local fractional homotopy perturbation method for solving non-homogeneous heat conduction equations in fractal domains. Entropy 17, 6753–6764 (2015)

    MathSciNet  Google Scholar 

  28. Al-Ghafri, K.S., Rezazadeh, H.: Solitons and other solutions of \((3 + 1)\)-dimensional space–time fractional modified KdV–Zakharov–Kuznetsov equation. Appl. Math. Nonlinear Sci. 4(2), 289–304 (2019)

    MathSciNet  Google Scholar 

  29. Cattani, C., Pierro, G.: On the fractal geometry of DNA by the binary image analysis. Bull. Math. Biol. 75(9), 1544–1570 (2013)

    MathSciNet  MATH  Google Scholar 

  30. Yang, A.-M., Zhang, Y.-Z., Cattani, C., Xie, G.-N., Rashidi, M.M., Zhou, Y.-J., Yang, X.-J.: Application of local fractional series expansion method to solve Klein–Gordon equations on Cantor sets. Abstr. Appl. Anal. 2014, Article ID 372741 (2014)

    MathSciNet  MATH  Google Scholar 

  31. Cattani, C., Srivastava, H.M., Yang, X.-J.: Fractional dynamics. In: Fractional Dynamics, pp. 1–5. de Gruyter, Berlin (2015)

    Google Scholar 

  32. Gao, W., Yel, G., Baskonus, H.M., Cattani, C.: Complex solitons in the conformable \((2+1)\)-dimensional Ablowitz–Kaup–Newell–Segur equation. AIMS Math. 5(1), 507–521 (2020)

    Google Scholar 

  33. Durur, H., Senol, M., Kurt, A., Tasbozan, O.: Approximate solutions of the time fractional Kadomtsev–Petviashvili equation with conformable derivative. Erzincan Univ. J. Sci. Technol. 12(2), 796–806 (2019)

    Google Scholar 

  34. Yang, X.J., Baleanu, D., Lazarević, M.P., Cajić, M.S.: Fractal boundary value problems for integral and differential equations with local fractional operators. Therm. Sci. 19(3), 959–966 (2015)

    Google Scholar 

  35. Brzeziński, D.W.: Review of numerical methods for NumILPT with computational accuracy assessment for fractional calculus. Appl. Math. Nonlinear Sci. 3(2), 487–502 (2018)

    MathSciNet  Google Scholar 

  36. Yang, X.J., Gao, F.: A new technology for solving diffusion and heat equations. Therm. Sci. 21(1 Part A), 133–140 (2017)

    Google Scholar 

  37. Singh, J., Kumar, D., Hammouch, Z., Atangana, A.: A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. Appl. Math. Comput. 316, 504–515 (2018)

    MathSciNet  MATH  Google Scholar 

  38. Dokuyucu, M.A., Celik, E., Bulut, H., Baskonus, H.M.: Cancer treatment model with the Caputo–Fabrizio fractional derivative. Eur. Phys. J. Plus 133, Article ID 92 (2018)

    Google Scholar 

  39. Goswami, A., Singh, J., Kumar Sushila, D.: An efficient analytical approach for fractional equal width equations describing hydro-magnetic waves in cold plasma. Physica A 524, 563–575 (2019)

    MathSciNet  Google Scholar 

  40. Gao, W., Veeresha, P., Prakasha, D.G., Baskonus, H.M., Yel, G.: New numerical results for the time-fractional Phi-four equation using a novel analytical approach. Symmetry 12, Article ID 478 (2020)

    Google Scholar 

  41. Sulaiman, T.A., Bulut, H., Atas, S.S.: Optical solitons to the fractional Schrödinger–Hirota equation. Appl. Math. Nonlinear Sci. 4(2), 535–542 (2019)

    MathSciNet  Google Scholar 

  42. Durur, H., Taşbozan, O., Kurt, A., Şenol, M.: New wave solutions of time fractional Kadomtsev–Petviashvili equation arising in the evolution of nonlinear long waves of small amplitude. Erzincan Univ. J. Sci. Technol. 12(2), 807–815 (2019)

    Google Scholar 

  43. Dananea, J., Allalia, K., Hammouch, Z.: Mathematical analysis of a fractional differential model of HBV infection with antibody immune response. Chaos Solitons Fractals 136, Article ID 109787 (2020)

    MathSciNet  Google Scholar 

  44. Gao, W., Veeresha, P., Prakasha, D.G., Baskonus, H.M.: Novel dynamical structures of 2019-nCoV with nonlocal operator via powerful computational technique. Biology 9(5), Article ID 107 (2020)

    Google Scholar 

  45. Owolabi, K.M., Hammouch, Z.: Spatiotemporal patterns in the Belousov–Zhabotinskii reaction systems with Atangana–Baleanu fractional order derivative. Physica A 523, 1072–1090 (2019)

    MathSciNet  Google Scholar 

  46. Sulaiman, T.A., Yel, G., Bulut, H.: M-fractional solitons and periodic wave solutions to the Hirota–Maccari system. Mod. Phys. Lett. B 33, Article ID 1950052 (2019)

    MathSciNet  Google Scholar 

  47. Yokus, A., Gulbahar, S.: Numerical solutions with linearization techniques of the fractional Harry Dym equation. Appl. Math. Nonlinear Sci. 4(1), 35–42 (2019)

    MathSciNet  Google Scholar 

  48. Kumar, D., Singh, J., Baleanu, D.: Analysis of regularized long-wave equation associated with a new fractional operator with Mittag-Leffler type kernel. Physica A 492, 155–167 (2018)

    MathSciNet  Google Scholar 

  49. Gao, W., Veeresha, P., Baskonus, H.M., Prakasha, D.G., Kumar, P.: A new study of unreported cases of 2019-nCOV epidemic outbreaks. Chaos Solitons Fractals 138, Article ID 109929 (2020)

    MathSciNet  Google Scholar 

  50. Durur, H., Tasbozan, O., Kurt, A.: New analytical solutions of conformable time fractional bad and good modified Boussinesq equations. Appl. Math. Nonlinear Sci. 5(1), 447–454 (2020)

    MathSciNet  Google Scholar 

  51. He, J.H.: Approximate analytical solution for seepage flow with fractional derivatives in porous media. Comput. Methods Appl. Mech. Eng. 167, 57–68 (1998)

    MathSciNet  MATH  Google Scholar 

  52. He, J.H.: Variational iteration method: a kind of nonlinear analytical technique: some examples. Int. J. Non-Linear Mech. 34, 699–708 (1999)

    MATH  Google Scholar 

  53. He, J.H.: A new approach to nonlinear partial differential equations. Commun. Nonlinear Sci. Numer. Simul. 2(4), 203–205 (1997)

    MathSciNet  Google Scholar 

  54. He, J.H.: A variational iteration approach to nonlinear problems and its applications. Mech. Appl. 20(1), 30–31 (1998) (in Chinese)

    MathSciNet  Google Scholar 

  55. Wazwaz, A.M.: The variational iteration method: a reliable analytic tool for solving linear and nonlinear wave equations. Comput. Math. Appl. 54, 926–932 (2007)

    MathSciNet  MATH  Google Scholar 

  56. Wazwaz, A.M.: The variational iteration method for exact solutions of Laplace equation. Phys. Lett. A 363(4), 260–262 (2007)

    MathSciNet  MATH  Google Scholar 

  57. Momani, S., Abuasad, S.: Application of He’s variational iteration method to Helmholtz equation. Chaos Solitons Fractals 27, 1119–1123 (2006)

    MathSciNet  MATH  Google Scholar 

  58. Mousa, M.M., Kaltaye, A., Ragab, S.F.: Investigation of a transition from steady convection to chaos in porous media using piecewise variational iteration method. World Acad. Sci., Eng. Technol. 58, 1088–1097 (2009)

    Google Scholar 

  59. Odibat, Z.M., Momani, S.: Application of variational iteration method to nonlinear differential equations of fractional order. Int. J. Nonlinear Sci. Numer. Simul. 7(1), 27–34 (2006)

    MathSciNet  MATH  Google Scholar 

  60. Sweilam, N.H.: Fourth order integro-differential equations using variational iteration method. Comput. Math. Appl. 54(7–8), 1086–1091 (2007)

    MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The numerical calculations reported in this paper were fully/partially performed at Harran University High Performance Computing Center (Harran HPC resources).

Availability of data and materials

Data sharing is not applicable to this paper as no datasets were generated or analyzed during the current study.

Funding

There is not any funding.

Author information

Authors and Affiliations

Authors

Contributions

Authors contributed to the paper equally. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Haci Mehmet Baskonus.

Ethics declarations

Ethics approval and consent to participate

Authors declare that there is not any ethical approval.

Competing interests

Authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gao, W., Baskonus, H.M. & Shi, L. New investigation of bats-hosts-reservoir-people coronavirus model and application to 2019-nCoV system. Adv Differ Equ 2020, 391 (2020). https://doi.org/10.1186/s13662-020-02831-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13662-020-02831-6

Keywords