The analysis of a time delay fractional COVID-19 model via Caputo type fractional derivative.

Novel coronavirus (COVID-19), a global threat whose source is not correctly yet known, was firstly recognised in the city of Wuhan, China, in December 2019. Now, this disease has been spread out to many countries in all over the world. In this paper, we solved a time delay fractional COVID-19 SEIR epidemic model via Caputo fractional derivatives using a predictor-corrector method. We provided numerical simulations to show the nature of the diseases for different classes. We derived existence of unique global solutions to the given time delay fractional differential equations (DFDEs) under a mild Lipschitz condition using properties of a weighted norm, Mittag-Leffler functions and the Banach fixed point theorem. For the graphical simulations, we used real numerical data based on a case study of Wuhan, China, to show the nature of the projected model with respect to time variable. We performed various plots for different values of time delay and fractional order. We observed that the proposed scheme is highly emphatic and easy to implementation for the system of DFDEs.

INTRODUCTION

Novel coronavirus (COVID‐19), a global threat whose source is not correctly yet known, was firstly recognised in the city of Wuhan, China, in December 2019. , Now, this disease has been spread out to many countries in all over the world. From the date of its origin, it grows exponentially in mankind and infected more than 20,254,685 with 738,930 deaths and 13,118,618 recoveries on 11 August throughout the globe. For the preventions to this disease, social control measures have been extended up by increased public awareness such as through social or physical distancing measures, good hygiene and not walking out in the open environment. Billions of people have been infected by this virus, and in some particular cases, the virus has also been recognised in animals. The symptoms of the coronavirus are fever, sneeze, a runny nose, fatigue, dry cough, bilateral lung infiltration to severely ill and breathing problems, in which fever and dry cough are two most common symptoms. So many newly reported cases have been associated to the travel history from an epidemic region or a meeting history to people from the region in the early outbreak phases. Because a person who was infected by COVID‐19 can travel to any other region and spread the virus. This coronavirus outbreak has received considerable global attention when, since 31 Dec 2019, the Wuhan Municipal Health Commission reported 27 cases of viral pneumonia, with 7 critically ill cases. Many researchers and doctors have suggested their specialisation to the study of this virus. COVID‐19 has been studied from different points of view, in many fields of study, including epidemiology, microbiology, environmental and occupational health, virology, veterinary science, economics, sociology and media studies. In reaction to early outbreaks of the virus, the United States, Korea and China come out as the leading countries in COVID‐19 research. Other research groups focused on the improvement of healthcare systems and virus prevention propounding that the way to restraining virus transmission is the reduction of human contact and social distance by dint of the enclose of public transport, school closures, the adjournment of common activities, etc. To date, there is no specific vaccine or clinically established treatment for coronavirus. Nevertheless, according to the suggestions of different prudent infectious disease specialists, several countries have started using chloroquine combined with azithromycin as alternative drugs. But still there are some countries that successfully controlled the pandemic. These include Iceland, New Zealand, Fiji, Vatican City, Tanzania, Montenegro, Seychelles, Papua New Guinea and Mauritius. On 28 February 2020, New Zealand reported its first coronavirus case and then enforced the one of the strictest lockdowns in the world. The United States, Brazil, India and Peru are the countries which are most infected by this virus. Particularly, the United States is one of the most infected country of this virus with 5,251,446 cases, 166,192 deaths and 2,715,934 recoveries on 11 August. In Brazil, 3,057,470 cases with 101,857 deaths and 2,163,812 recoveries have been recorded. In India, there are 2,267,153 cases with 45,353 deaths and 1,581,640 recoveries founded on August 11. Because of this pandemic, lots of social programmes and activities have been cancelled or postponed. T‐20 cricket world cup 2020 in Australia, Summer Olympics which was planned to be held in Tokyo, has been postponed. One of the most popular cricket tournament, Indian Premier League (IPL), has been shifted from India to UAE.

Mathematical models estimate the number of cases in best and worst case scenarios when disease dynamics are unclear and can also be used to estimate the effects of preventive (healthcare) measures. A lot of research articles have been come in literature to analyse the effects of COVID‐19 via mathematical modelling. , , , , In Khan et al, authors formulated a new mathematical model for the dynamics of COVID‐19 with quarantine and isolation. They studied the dynamics of the COVID‐19 using fractal–fractional Atangana–Baleanu derivative. A study on the mathematical modelling of the impact of nonpharmaceutical interventions on the dynamics of COVID‐19 with optimal control analysis is given in Ullah and Khan. A study on solution of a COVID‐19 model via new generalised fractional derivative (FD) is done by Erturk and Kumar. They used a modified predictor–corrector scheme to perform the numerical simulations with the real data of Wuhan, China. The applications of fractional calculus have been received in various branches of Science and Engineering. , , , , Unlike an ordinary derivative (OD) operator, an FD operator is non‐local in nature. Due to this non‐local behaviour of the FD operator, it can formulate processes having memory and hereditary properties. Recently, an attention to the delay fractional differential equations (DFDEs) has considered cause of their applications in the mathematical modelling of real‐world problems. The delay differential equation is a differential equation in which the derivative of the function at any time depends on the solution at previous time. It is well known that an ordinary delay differential equation has a unique local solution under some Lipschitz conditions; furthermore, using continuation property, one can derive global solutions as well, but in the noninteger, the existence of unique solutions (local and global) is more intricate because of the noninteger order feature of the equation which implies history dependence of the solutions; hence, among others, the continuation property is not applicable.

The aim of this article is to study the time delay fractional COVID‐19 epidemic model using a real numerical data of a case study of Wuhan, China, from the literature to show the nature of the given model. We also performed some analysis to show the role of time delay parameter by the help of Caputo FD with predictor–corrector algorithm. The paper is formulated as follows. In Section 2, we recall some important definitions of the FDs. Section 3 is devoted for the description of the ordinary differential equations (ODE) model following the fractional order model. Existence and uniqueness analysis of the problem are performed in Section 4. Solution of the projected model is done in Section 5. Simulation results are performed in Section 6. A conclusion completes the paper.

PRELIMINARIES

Here, we remind some basic definitions and properties.

(Podlubny The Riemann–Liouville (R‐L) definition of noninteger order integral of order of a function is defined by

(Podlubny The R‐L definition of noninteger order derivative of order of a function is defined by where and [ is the integer part of

(Podlubny The Caputo definition of noninteger order derivative of order of a function is defined by where and [ is the integer part of

(Podlubny The Mittag–Leffler type function with one parameter is defined as follows

MODEL DESCRIPTION

Description of the ODE model

There are so many mathematical models have been introduced in the literature to study the outbreaks of COVID‐19. So many researchers have analysed various type of models to study the dynamics of COVID‐19 with the case study of different particular countries. In this section, we analyse the time delay ordinary model studied by Cakan. to introduce the COVID‐19 epidemic. The considered model consists of four compartments with individuals of susceptible , exposed , infectious and recovered . The author presented and derived the projected model in the sense of ordinary time delay differential equations as follows:

In the given ordinary model (4), b is the birth rate, is the contact rate of susceptible to infected persons and natural death rate denoted by d. and are the progression rate of exposed humans into infectious population and rate of exposed to removed, respectively. The brief description of all other parameters with the numerical values is given in Table 1. The function presents the available opportunities level by healthcare system to public who are infected at time t. It can be observed that all hospital facilities are almost consumed away when tends to zero and all hospital facilities (opportunities) can be used fully when with respect to the time variable t. is a time delay corresponding the latent period of the COVID‐19.

TABLE 1

Parameter values for simulations

Parameter	Description	Value/range	Reference
b	Birth rate	3,210	Estimated
β	Contact rate susceptible to infected	0.62 × 10−8	Yang and Wang 24
d	Natural death rate	3.57 × 10−5	Khan et al 11
γ	Rate of exposed to infected	0.143	Yang and Wang 24
δ	Rate of exposed to removed	0.006	Fitted
c	Level of available opportunities by health care systems	[0, 1]	Yang and Wang 24
α 1	Natural recovery rate of the infectious class	0.0005	Fitted
α 2	Recovery rate of the infectious class	0.0667	Yang and Wang 24
μ 1	Minimum disease‐induced death rate	0.01	Yang and Wang 24
μ 1 + μ 2	Maximum disease‐induced death rate	0.02	Fitted

Description of the Caputo fractional model

The theory of the time DFDEs is a well known phenomena of fractional calculus. Now, we generalise the above ODE model in the Caputo FD sense. Here, we are including the death equation also. In this generalisation, we replace the OD operator by the Caputo fractional operator So the generalisation of the given ordinary time delay differential equation system into the time DFDEs system is as follows: where denotes the Caputo FD operator.

By means of the next generation matrix method, here, we find the basic reproduction number for the fractional model (5) as follows:

Let us assume Then, the model (5) can be written as where and is the unique disease‐free equilibrium point of the model (6).

Then, Jacobian matrices at the disease‐free equilibrium point of and by looked to the FDs with respect to I and S are obtained as and respectively. So F and which are the new infection terms and the remaining transfer terms of the model, respectively, are determined as and Also, the characteristic polynomial of is It follows from Van den Driessche and Watmough that the basic reproduction number of the model (5) is the spectral radius of the next generation matrix In that case,

The equivalent compact form of the above system (5) is as follows: with the initial conditions taken as and

MATHEMATICAL ANALYSIS OF THE FRACTIONAL MODEL

Existence and uniqueness analysis

In this section, we give the existence of unique solution for the projected fractional time delay COVID‐19 model by the help of the consequences of fixed point theory. In this regard, many results have been given in the literature, and here, we are following the procedure proposed by Cong and Tuan. We show the analysis for , and for other equations of the system (8), it will be similar. Let us consider the fractional time delay equations with the initial condition where is continuous.

( be the n‐dimensional Euclidean space defined with a norm ‖.‖)

(Cong and Tuan The function ( space of continuous functions with the sup norm ‖.‖ is a solution of the initial value problem (IVP) (Equations (9), (10)) on the interval [−  if it is a solution of the delay integral equation with the initial condition

Note: In this particular way to show the existence of unique global solutions, we don't need to require Lipschitz property of with respect to the time delay variable of , but only the Lipschitz property of with respect to the nondelay (second) variable.

(Existence and uniqueness of global solutions) Assume that is continuous and agree with the following Lipschitz condition with respect to the nondelay variable: there exists a non‐negative continuous function such that

Then, the IVP (Equations (9), (10)) has a unique global solution on the interval [−

According to Lemma (1), Equation (9) with Equation (10) is equivalent to the initial value problem (11)–(12). First, we take the case 0 <  In this case, Equation (11) has the form: For this integral equation, by Tisdell [ Theorem 6.4, p. 310], there exists a unique solution on the interval [0,  Identify that solution by and put

Then, is the unique solution of the problem (11)–(12) on [−

Now, in another case when we break the interval [0,  into [0,  where and 0 ≤  On the interval [−  using the same arguments as above, we can derive a unique solution of the IVP (Equations (11), (12)), which is identified by We will prove the existence of the unique solution on the interval [−  by induction. Let us assume that the problem (11)–(12) has a unique solution on the interval [−  for some 1 ≤ . We denote that solution by On [ we define an operator as follows: Let be a positive constant satisfying On the space we define a new metric by where is the Mittag–Leffler function which is defined in Definition 4. Then, the space equipped the metric is complete. We will show that the operator is contractive on . Indeed, for any and any , we have

This implies that for all Therefore, for all By a statement of the Banach fixed point theorem, there exists a unique fixed point of in Put

Then, is the unique solution of the problem (11)–(12) on [−

Finally, on the interval [ we derive an operator by Let be a positive constant satisfying On the space we define a new metric by and repeating arguments as above, we can show that the operator has a unique fixed point on [ Define a function

It is evident that is the unique solution of the problem (11)–(12) on the interval [−

SOLUTION OF THE PROJECTED MODEL USING PREDICTOR– CORRECTOR ALGORITHM

We know that the numerical methods or techniques used for solving ordinary differential equations cannot be used directly to solve noninteger order differential equations. Lots of numerical methods have been used to solve DFDEs arising in biology. In this section, we find the solution of the projected model by the help of Adams–Bashforth–Moulton predictor–corrector scheme described in Bhalekar and Daftardar‐Gejji for solving DFDEs.

Let us consider the fractional delay differential equation

Consider a uniform grid , where n and are integers such that and Let and note that

Suppose we have already calculated the approximations and we want to calculate using the volterra integral equation equivalent to the Equations (20a) and (20b)

We use approximations for in Equation (23). Further, the integral in Equation (23) is evaluated using product trapezoidal quadrature formula. The corrector formula is thus where The unknown term present on both sides of Equation (24) and due to non‐linearity of Equation (24) cannot be solved explicitly for So we replace the term on the right hand side by an approximation called predictor. Product rectangle rule is used in Equation (24) to evaluate predictor term where So from the above calculations, the corrector formulas for all five equations of system (8) are where Similarly, the predictor terms are where

SIMULATION RESULTS

To perform numerical simulations, we use parameter values based on a case study of Wuhan, China, cited from literature and with some assumptions, summarize in Table 1.

We also use the following initial conditions , , , and for Table 1. The function is taken as Initial population is considered as 89,995,587, and birth rate is

Figure 1A exemplifies the behaviour of achieved results by given solution procedure for at different values of fractional order . In Figure 1B, we analyse the exposed individuals for different fractional order with respect to time t. Figure 1C, D and E analyses the nature of infectious, recovered and deaths classes at and 1, respectively. Figure 1A–E examines the nature of all given classes at time delay Similarly, the graphical simulations from Figure 2A–E are for time delay for parameter values given in Table 1. In Figure 2A, we show the nature of susceptible individuals, in Figure 2B, exposed classes; in Figure 2C, infectious individuals; in Figure 2D, recovered and in Figure 2E, the deaths, with respect to time variable t for parameter values used from Table 1. In this same manner, we observed the plots for time delay in Figure 3A–E. All graphs are computed using Mathematica software. For the all above simulations, we used time delay and 10, respectively, and time step was 0.05. From the cited figures, we can observe that the given model exceedingly depends on the order and gestures more degree of flexibility. Moreover, the fractional method gives more interesting results than the integer order model and permit to better examine the obtained results. We observed that time delay variable plays an important role in this COVID‐19 dynamics. The biological meaning of time delay in COVID‐19 is cleared from the above simulations.

FIGURE 1

Nature of the given classes at time delay= 2 for different fractional order values

FIGURE 2

Nature of the given classes for time delay= 6 at different fractional order values

FIGURE 3

Nature of the given classes for time delay= 10 at different fractional order values

In Figure 4A–C, we study the nature of basic reproductive number at different fractional order values and for different time delay We can see that the maximum value of basic reproductive number for the given time range is near about 1.038. Numerically, we have calculated that when tends to one, tends to 1.032, and when tends to zero, then tends to 3.880. It is clear that the recovery rate and death rate are directly related to whether healthcare capacity is exceeded or not. For the all above graphical simulations, we can observe that the given model works well to study the dynamics of reproductive number respect to the time t. We have clearly observed that the parameter c plays a very important role in the simulations and time delay can show the effect of healthcare opportunities.

FIGURE 4

Nature of basic reproductive number for different time delay

CONCLUSIONS

In this paper, we solved a time delay fractional COVID‐19 compartmental model via Caputo FDs applying predictor–corrector algorithm. After preliminaries and the model description, the existence and uniqueness analysis of the given time delay fractional system are established by the applications of fixed point theory (particularly, a mild Lipschitz condition using properties of a weighted norm, Mittag–Leffler functions and the Banach fixed point theorem). All necessary graphical simulations are done to specified the nature of the achieved solutions in Caputo noninteger order derivative sense. We analysed the role of time delay in the coronavirus epidemic by the graphical simulations for different values. We also presented the plotting of the of basic reproductive number at different fractional order values and various values of time delay. The projected scheme is strong and highly credible in finding the solution to delay fractional models of biological, physical and medical importance.

CONFLICT OF INTEREST

This work does not have any conflicts of interest.

AUTHORS CONTRIBUTION

Conceptualization, formal analysis, investigation, methodology, resources, visualization and writing—original draft: Pushpendra Kumar. Conceptualization, investigation, software and writing—review & editing: Vedat Suat Erturk.