key: cord-0982816-q597thb0
authors: Alkahtani, Badr Saad T.; Jain, Sonal
title: Numerical analysis of COVID-19 model with constant fractional order and variable fractal dimension
date: 2020-12-10
journal: Results Phys
DOI: 10.1016/j.rinp.2020.103673
sha: 44adbe0e134b00661f346bc3d442132b6fae2490
doc_id: 982816
cord_uid: q597thb0

This work has considered a mathematical model describing the spread of COVID-19 in a given population. The model comprised 5 systems of equations that take into account different classes describing the impact of COVID-19 in a given population. The time differential operator was replaced with three different types of nonlocal operators. These operators are defined as the convolution of variable order fractal differential operators with different kernels including power law, exponential decay law, and Mittag-Leffler functions. We presented the well-poseness of the models for different differential operators that were presented in detail. A novel numerical scheme was used to solve numerically the system and numerical simulations were provided.

Differential operators with non-local characters have been noticed to be able to replicate several complexities occurring in nature in the last decades. One can mention among which fractional differential and integral operators with power law, exponential decay and Mittag-Leffler kernels on the other hand differential operators defined as convolution of power law, exponential decay and Mittag-Leffler kernels with fractal differential operator, these cases are called fractal-fractional operators [1, 2, 5, 10] . Finally fractal-fractional differential operators with variable fractal orders. These last one are considered to be adequate to modelling complex real world problems, for example real world problems with anomalous patterns could be replicated using variable order differential operators. These operators have been successfully applied in many academics disciplines, however much attention have not been devoted to epidemiologic models [3, 18, 15, 4, 12] . They could be very useful tools to modelling epidemiologic problems as the spread of some infectious diseases are anomalous therefore neither the classical differentiation nor fractional with constant can be applied in these cases. Very recently, the world have been surprised with an outbreak of a fatal disease called COVID-19, which was first observed in Wuhan, China, December 2019 [9, 7, 8] . From this time to 07 May 2020 the disease has infected 3836183, the total number of deaths from COVID-19 is about 265 364 of course this number is for those registered or declared by each nation around the globe. However, 1307 608 have been recovered. The exponential spread of this disease leads humans no choice than to undertake serious researchers activities in all field of science. In applied mathematics, many new mathematical models have been suggested, some including fractional differential and fractal fractional operators. In this paper, we aim to revert the model suggested by Shafiq and Atangana using the fractal-fractional with variable fractal-order [4] . One of the major concern of differential operators with variable orders is perhaps their solvability as analytical methods cannot be used to provide their exact solutions. Thus, numerical methods are adequate to providing approximate solutions to such models. In the last decades, many numerical methods have been provided all with their advantages and limitations. In the case of nonlinear equations, it is known that the Adams-Bashforth is powerful numerical scheme to provide approximate solutions of nonlinear equations [13, 11, 16, 17, 19] . However, the case with fractional differential has some limitation thus, Atangana and Toufit [20] suggested an alternative efficient numerical scheme to be used, and the method has been used in various nonlinear equations arising in many fields of science, engineering and technology. In this paper, we aim at using such scheme to solve the model of COVID-19 suggested by Altaf and Atangana [14] , where the time derivative is reverted to fractal-fractional with fractal variable order. The paper is organized as follow, we start with the model description, then, the wellposeness of the model is presented in the case of fractal-fractional with variable order fractal and finally numerical analysis for different cases.

Definition 0.1. A discontinuous media can be described by fractal

The more generalized version is given as:

Definition 0.2. [6] If u(t) is continuous is an closed interval [a,b], then the fractal integral of u with order ϑ is defined as:

Definition 0.3. [14] Let f be a differentiable function. Let ϑ be a constant fractional order, such that 0 < ϑ ≤ 1. Let ξ(t) be continuous function ξ(t) > 0. then a fractional derivative of f with order ϑ and fractal dimension ξ(t) is defined as

The above definition is with power law kernel. With exponential kernel we have

with Mittag-Leffer Law, we have

Corollary 0.1. [14] If ξ(t) is differentiable then

Definition 0.4. [14] The fractional integral operators associate to the suggested differential operators are given as: For power law we have

For exponential decay law, we have

For Mittag-Leffler law we have

The total population of people is denoted by N p which is classified further into five subgroups such as S p , E p , I p , A p and R p which represent respectively, the susceptible, exposed, infected (symptomatic), asymptotically infected, and the recovered or the removed people. Considering that the 2019-nCoV can be imported in short time to the seafood market with enough source of virus and thus, without loss of generality, ignoring of the interaction among bats and hosts, then model can be represent to the system below:

with the initial conditions (2.

2)

The total dynamics of the people can be obtained by adding the first five equations of the model 2.1, given by

The feasible reason of the model 2.1 is given by

The birth and natural death rate of the people is given by the parameters π p and μ p respectively. The susceptible people S p will be infected through sufficient contacts with the infected people I p through the term given by η p S p I p where the η p is the disease transmission coefficient [14] . The transmission among the asymptotically infected people with health people could take place at form ψη p S p A p , where ψ the transmissibility multiple of A p to that I p and ψ ∈ [0,1], when psi = 0, no transmissibility multiple will exists and hence vanish, and if ψ = 1, then the same will take place like I p infection. The parameter θ p is the proportion of asymptomatic infection. The parameters ω p and ρ p respectively represent the transmission rate after completing the incubation period and becomes infected, joining the class I p and A p . The people in the symptomatic class I p and asymptomatic class A p joining these class R p with the removal or recovery rate respectively by τ p and τ ap . The class M which is denoted be the reservoir or the seafood place or market. The susceptible people infected after the interaction with M, given by η w MS p , where η w is the disease transmission coefficient from M to S p . The parameters ϱ p and ϖ p of the infected symptomatic and asymptotically infected respectively contributing the virus into the seafood market M. The removing rate of the virus from the seafood market M is given by the rate π.

So we can write this model in the fractal fractional differentiation sense as follows:

with the initial conditions

(2.6)

we now present that for each case of S p(0) , E p (0), I p (0), A p (0), R p (0) and M(0) are positive then the solution are also positives. We shall start with power law case: To do this, we define the following norm,

we suppose that the all the classes constituting function of time have the same sign, therefore any product of the classes is positive. 

and finally

, (3.13) we consider the case with exponential decay law

(3.20)

and

Now we consider the case with the generalized Mittag-Leffler Kernel with same hypothesis, we have

(3.24)

Then following the same continue

(3.25)

(3.26)

(3.27)

(3.28)

(3.29)

Recently Atangana and Shafiq [4] suggested an alternative numerical scheme for solving fractal fractional differential equation. So in this section we can solve this method by using that numerical scheme as follows:

Let us consider the following cauchy problem

at the given point t m+1 , here we chose m = 0,1,2,…, then the above equation can be written as

For simplicity we consider

within [t j , t j+1 ], we can apply the Lagrange polynomial interpolation, Thus within this interval, we approximate

Then replacing above in original equation, we obtain

As given in the Atangana and Toufik [20] , the above equation can be reformulated as after integration

replacing F(t j , y(t j )) and F(t j− 1 ), y(t j− 1 ) by their respective values, we obtain

So for applying the procedure in (4.10), our model taking the shape below: 

(4.13)

(4.14)

We next consider the corresponding Cauchy with exponential decay kernel. The corresponding Volterra type is given by

So at t = t m+1 , we have

at t m = t, we have

Therefore

Now using the procedure by Adams Bashforth, we obtain

so our model can be written in exponential law form as

so we can solve the (2.5) equation numerically with constant fractional order and variable fractal dimension (see [14] ) by using above solution and get

The existence and uniqueness of this solution is given in the paper Atangana and Shafiq.

Now we consider the case of Mittag Leffler Kernel. Using the corresponding integral operator we convert the above equation as

So we consider the following general nonlinear equation

at t m+1 we have

For simplicity, we put

we approximate g(τ, y(τ)) within the interval [t j , t j+1 ]

replacing g(τ, y(τ))⋍q j (τ) by its value and integrating, we obtain

(4.44) 

Mathematical models with non local variable orders operators have been known to replicate sometime accurately as they are able to include into mathematical model real representation of complex patterns observed in nature. However, mathematical models depicted with such different and integral operators cannot be solved analytically due to highly non linearity of the operators. To have a solution of such model only numerical approximations are useful. In this paper a system of nonlinear ordinary differential equations with non local variable order operators were considered. The system describe the spread of COVID-19 in a given population. Numerical schemes based on Lagrange polynomial was used to derive numerical solutions for three cases, numerical simulations were depicted for different values of frac.

The first and second authors have both agreed to use the mathematical models and the methodology. Both authors did numerical schemes. The first author did numerical simulations. Both authors wrote the first original draft.

This article does not contain any studies with human or animal subjects.

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

A new numerical approximation of the fractal ordinary differential equation

Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system

The role of power decay, exponential decay and Mittag-Leffler functions waiting time distributions: application of cancer spread

Differential and integral operators with constant fractional order and variable fractional dimension

Models of fluid owing in non-conventional media: New numerical analysis

Investigation on fractional and fractal derivative relaxation-oscillation models

A mathematical model for simulating the transmission of Wuhan novel Coronavirus

Novel coronavirus: where we are and what we know

Short term outcome and risk factors for adverse clinical outcomes in adults with severe acute respiratory syndrome (SARS)

Chaos and multiple attractors in a fractal-fractional Shinrikis oscillator model

New Two Step Laplace Adam-Bashforth Method for Integer an Non integer Order Partial Differential Equations, Numerical Methods for Partial Differential Equations

Analytical investigation of third grade nanofluidic flow over a riga plate using Cattaneo-Christov model

Numerical analysis for the fractional diffusion and fractional Buckmaster's equation by two step Adam-Bashforth method

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

A mathematical model of treatment and vaccination interventions of pneumococcal pneumonia infection dynamics

Mathematical analysis and numerical simulation of patterns in fractional and classical reaction-diffusion systems

Robust and adaptive techniques for numerical simulation of nonlinear partial differential equations of fractional order

Analysis of Lassa hemorrhagic Fever Model with non-local and non-singular fractional derivatives

Influence of single-and multiwall carbon nanotubes on magnetohydrodynamic stagnation point nanofluid flow over variable thicker surface with concave and convex effects

New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models. A

The authors extend their appreciation to the Deputyship for Research & Innovation, "Ministry of Education" in Saudi Arabia for funding this research work through the project number IFKSURG-1437-017.

(4.39)Now replacingf (t j , y(t j )) (4.40)(4.42) so our model can be written in exponential law form asso we can solve the (2.5) equation numerically with constant fractional order and variable fractal dimension by using above solution and get(4.47)(4.48)(4.49)

In this section, using the obtained numerical solutions we present in this section numerical simulation for various fractional order and variable fractal order. The numerical simulations are depicted in Fig. 1-7 .