key: cord-266090-f40v4039
authors: Gao, Wei; Baskonus, Haci Mehmet; Shi, Li
title: New investigation of bats-hosts-reservoir-people coronavirus model and application to 2019-nCoV system
date: 2020-08-03
journal: Adv Differ Equ
DOI: 10.1186/s13662-020-02831-6
sha: 
doc_id: 266090
cord_uid: f40v4039

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.

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 [4] [5] [6] . 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 [10] [11] [12] [13] [14] [15] [16] .

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

where n p , m p , b p , κ, b w , δ p , ε, c are real nonzero constants. The initial conditions for this system are shortly given by

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/ω B is the incubation period of bat infection, and 1/γ 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 . In this paper, we investigate the numerical distributions of 2019-nCoV according to time with the help of several approaching terms of VIM.

In 1999, VIM, one of the most powerful numerical methods, was firstly developed by He [51] [52] [53] [54] 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 [14] [15] [16] [17] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] . We consider differential equations of the form

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) [51] [52] [53] [54] :

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 u n is considered as a restricted variation, which means δ 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 ≥ 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

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:

where k = 0, 1, 2, 3, . . . . It produces the stationary condition λ = -1.

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

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:

Thus we get the first approach of u n , v n , y n , f n , r n , w n as follows:

f 1 (t) = β 4 + τ 4 t, r 1 (t) = β 5 + τ 5 t,

where, for simplicity, we have taken

Now we obtain the second components of the variables u n , v n , y n , f n , r n , w n for k = 1:

y 2 (t) = β 3 + (τ 3 + τ 13 )t + τ 14 2 t 2 , f 2 (t) = β 4 + (τ 4 + τ 15 )t + τ 16 2 t 2 , r 2 (t) = β 5 + (τ 5 + τ 17 )t + τ 18 2 t 2 , w 2 (t) = β 10 + (τ 6 + τ 19 )t + τ 20 2 Figure 1 Two-dimensional surfaces of system (10) where

τ 14 = +τ 2 (-1 + δ p )ω p + τ 3 (γ p + m p ), 10 ,

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 .

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 suit- Figure 2 Two-dimensional surfaces of susceptible people and exposed people of system (10) Figure 3 Two-dimensional surfaces of symptomatic infected people and asymptomatic infected people of system (10) Figure 4 Two-dimensional surfaces of removed people (recovered or died people) and reservoir (the seafood area) of system (10) ably chosen values of the parameters, reported by WHO, we have plotted the numerical results. According to Figs. 1, 2, 3, 4, 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.

Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from patients with acute respiratory disease in Wuhan

National Health Commission of the People's Republic of China. Top expert: disease spread won't be on scale of SARS 2020

South China Morning Post. China coronavirus: three cities join Wuhan in quarantine lockdown as Beijing tries to contain deadly outbreak 2020

Strategies for mitigating an influenza pandemic

Effectiveness of control measures during the SARS epidemic in Beijing: a comparison of the Rt curve and the epidemic curve

School closure to reduce influenza transmission

Potent binding of 2019 novel coronavirus spike protein by a SARS coronavirus-specific human monoclonal antibody

The 2019-new coronavirus epidemic: evidence for virus evolution

Preliminary identification of potential vaccine targets for the COVID-19 coronavirus (SARS-CoV-2) based on SARS-CoV immunological studies

Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative

Is the world ready for the coronavirus? Editorial. The New York Times

China virus death toll rises to 41, more than 1,300 infected worldwide

New approach for the model describing the deathly disease in pregnant women using Mittag-Leffler function

Mathematical analysis and computational experiments for an epidemic system with nonlocal and nonsingular derivative

Modeling the mechanics of viral kinetics under immune control during primary infection of HIV-1 with treatment in fractional order

New numerical simulations for some real world problems with Atangana-Baleanu fractional derivative

A mathematical model for simulating the phase-based transmissibility of a novel coronavirus

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Mathematical analysis of dengue fever outbreak by novel fractional operators with field data

Application of fractional calculus to epidemiology

Modelling the spread of COVID-19 with new fractal-fractional operators: can the lockdown save mankind before vaccination?

HIV and shifting epicenters for COVID-19, an alert for some countries

A generalization of truncated M-fractional derivative and applications to fractional differential equations

Complex surfaces to the fractional (2 + 1)-dimensional Boussinesq dynamical model with local M-derivative

A review on harmonic wavelets and their fractional extension

Local fractional homotopy perturbation method for solving non-homogeneous heat conduction equations in fractal domains

Solitons and other solutions of (3 + 1)-dimensional space-time fractional modified KdV-Zakharov-Kuznetsov equation

On the fractal geometry of DNA by the binary image analysis

Application of local fractional series expansion method to solve Klein-Gordon equations on Cantor sets

Fractional dynamics

Complex solitons in the conformable (2 + 1)-dimensional Ablowitz-Kaup-Newell-Segur equation

Approximate solutions of the time fractional Kadomtsev-Petviashvili equation with conformable derivative. Erzincan Univ

Fractal boundary value problems for integral and differential equations with local fractional operators

Review of numerical methods for NumILPT with computational accuracy assessment for fractional calculus

A new technology for solving diffusion and heat equations

A fractional epidemiological model for computer viruses pertaining to a new fractional derivative

Cancer treatment model with the Caputo-Fabrizio fractional derivative

An efficient analytical approach for fractional equal width equations describing hydro-magnetic waves in cold plasma

New numerical results for the time-fractional Phi-four equation using a novel analytical approach

Optical solitons to the fractional Schrödinger-Hirota equation

New wave solutions of time fractional Kadomtsev-Petviashvili equation arising in the evolution of nonlinear long waves of small amplitude. Erzincan Univ

Mathematical analysis of a fractional differential model of HBV infection with antibody immune response

Novel dynamical structures of 2019-nCoV with nonlocal operator via powerful computational technique

Spatiotemporal patterns in the Belousov-Zhabotinskii reaction systems with Atangana-Baleanu fractional order derivative

M-fractional solitons and periodic wave solutions to the Hirota-Maccari system

Numerical solutions with linearization techniques of the fractional Harry Dym equation

Analysis of regularized long-wave equation associated with a new fractional operator with Mittag-Leffler type kernel

A new study of unreported cases of 2019-nCOV epidemic outbreaks

New analytical solutions of conformable time fractional bad and good modified Boussinesq equations

Approximate analytical solution for seepage flow with fractional derivatives in porous media

Variational iteration method: a kind of nonlinear analytical technique: some examples

A new approach to nonlinear partial differential equations

A variational iteration approach to nonlinear problems and its applications

The variational iteration method: a reliable analytic tool for solving linear and nonlinear wave equations

The variational iteration method for exact solutions of Laplace equation

Application of He's variational iteration method to Helmholtz equation

Investigation of a transition from steady convection to chaos in porous media using piecewise variational iteration method

Application of variational iteration method to nonlinear differential equations of fractional order

Fourth order integro-differential equations using variational iteration method