Theoretical and numerical analysis for transmission dynamics of COVID-19 mathematical model involving Caputo–Fabrizio derivative

This manuscript is devoted to a study of the existence and uniqueness of solutions to a mathematical model addressing the transmission dynamics of the coronavirus-19 infectious disease (COVID-19). The mentioned model is considered with a nonsingular kernel type derivative given by Caputo–Fabrizo with fractional order. For the required results of the existence and uniqueness of solution to the proposed model, Picard’s iterative method is applied. Furthermore, to investigate approximate solutions to the proposed model, we utilize the Laplace transform and Adomian’s decomposition (LADM). Some graphical presentations are given for different fractional orders for various compartments of the model under consideration.


Introduction
COVID-19 is the nuisance word today that affected individuals all over the world as it marched forward across continents and spread throughout all countries. Coronavirus caused a total loss of motion of the development on the planet, which hit the economy and mended nature. It has made dread, despondency and uneasiness in individuals with social separating and the main treatment has become face cover, hand washing, and chemicals that are delivered on the basis of an everyday schedule. Governments stand weak in the face of startlingly high death rates particularly with the old and individuals who have reactions from interminable illnesses, the clinical community addressing patients in danger, made researchers search for an antibody and specialists develop scientific models to restrict the spread of this fatal infection.
The historical backdrop of coronaviruses started in 1965 [1]. It was found in the cultures of the fetal trachea organs obtained from the respiratory system of an adult suffering from a common cold. Recent research has revealed that coronary respiratory infections frequently emerge in the spring and winter more than in the autumn and summer seasons. COVID-19 added up to 35 of overall respiratory viral vigor through epidemics.
A mathematical modeling of disease transmission is intended to understand the dynamic spread of infection between individuals; see [2][3][4][5][6] and the references therein. Actually, it has been observed that infectious disease models have many advantages for forestalling and pointing out rising infectious diseases as well as controlling them like COVID-19. The spreading of infections was found to follow an exponential function and the growth rate differs from 2 to 5 days in some regions. According to the latest statistics to the date of July-01-2020, confirmed infections with the coronavirus have exceeded 10,609,665 worldwide, while the number of deaths has reached 514,449, and the number of people recovered has risen to 5,817,869, according to the Worldometer website that specializes in counting COVID-19 victims.
Infectious sicknesses represent a major threat to people, and additionally to the nation's economy. A careful comprehension of the behavior of diseases plays a noteworthy role in the decline of the infection in the population. Execution of an advantageous technique contra ailment transportation is another challenge. Various disease models were created in the recent literature that permit us to more readily control the spread of the infection. The vast majority of these models are established by ordinary differential equations; see [7][8][9][10][11][12][13]. Nonetheless, as of lately the job of fractional calculus that deals with fractional orders has shown up, as it has a noticeable job in the understanding of real-world phenomena, just as in demonstrating modeling of the exact description of hereditary qualities and memory [14][15][16][17]. As of lately, it has been seen that fractional differential equations can be employed to model worldwide problems with more subtlety; see [18,19].
Recently, new fractional operators were developed that give a precise description of memory and have succeeded with regards to modeling infectious ailments; see [20][21][22]. The worldwide issue of the spread of the ailment attracted the consideration of analysts from different fields, which prompted the rise of various propositions to examine and predict the development of the epidemic [23,24]. Our contribution relates to the consideration of the best-known class, which is the classification that presently shows up in clinical diaries [25][26][27][28][29][30][31][32][33][34][35][36]. This new work incorporates various theoretical and practical analyses on examining the dynamic conduct of the speed of spreading the coronavirus infection (COVID-19) disease and how to lessen the spread of infection in the public arena, and numerical simulations are likewise considered in this work.
As is well known, fractional calculus has a vast range of applicability in expounding convoluted dynamical systems with memory effects in different areas of engineering, biological sciences, and social sciences. Furthermore, the fractional order derivatives given by Caputo and Riemann-Liouville contain a singular type kernel. Fractional differential operators in fact are definite integrals which geometrically represent a complete spectrum of the functions or accumulation. The singular kernel some time creates difficulty during numerical analysis. This is because of its local singular kernel. To overcome this difficulty, Caputo and Fabrizo in 2016 introduced the concept of a nonlocal nonsingular kernel type derivative; see [37]. Recently many authors have proved that the mentioned derivative has interesting features in the descriptions of many processes and phenomena in the thermal sciences; see [38][39][40][41][42]. Keeping the importance of fractional derivatives, recently authors have investigated some models of COVID-19 from different aspects; see [43,44]. Very re-cently, authors [45] studied the behavior of COVID-19 transmission through new control strategies, involving all possible conditions of human-to-human transmission.
Mathematical models of infectious disease under fractional order derivatives provide a comprehensive description of the global and local dynamics. Furthermore, such a kind of models in which fractional calculus is involved more precisely describes the phenomena in the best way. Motivated by [37] and [45], we will study the following COVID-19 model with the Caputo-Fabrizio fractional derivative: where 0 ≤ t ≤ T < ∞ and CF D θ denotes the Caputo-Fabrizio fractional derivative of order 0 < θ ≤ 1. The details of the given model are described as follows: • The total population is classified into five compartments of individuals as follows: E(t) is for exposed (uninfected but surrounded by infection); I(t) for infected (with obvious clinical symptoms but not critical); C(t) for critically infected; H(t) for hospitalized individuals; D(t) for dead individuals due to COVID-19. • B represents a birth rate of exposed individuals. • δ is the rate of natural mortality. • α 1 is the rate of individuals transmission from exposed to infected compartment. • α 2 is the rate of critical cases of infected individuals. • α 3 is the rate of critical infected hospitalized. • α 4 is the rate of hospitalized individuals which not recovered and stay in critical case. • α 5 is the rate of death in critically infected individuals class. • α 6 is the rate of death in infected individuals class. • α 7 is the rate of infected people due to spreading infection from a dead body. • α 8 is the rate of infected individuals which hospitalized without passing in critical case. • α 9 is the rate of recovered which individuals hospitalized and get exposed again. • α 10 is the rate of recovered in infected individuals class due to powerful immunity and get exposed again. We first establish the existence theory of the model under the said derivative via a Picard type analysis of fixed point theory. Then on using LADM, we derive a semi-analytical solution to the problem under consideration. Here we remark that treating Caputo-Fabrizo type differential equations by LADM is very rare in the literature.
Our manuscript is arranged as follows. Some useful fundamentals are given in Sect. 2. Further theoretical results are given in Sect. 3. Numerical results are presented in Sect. 4. Finally, a brief conclusion is given in Sect. 5.

Preliminaries
In this section, we recall some useful fundamentals related to fractional calculus.
where N (γ ) is the normalization function which defined by N (γ ) = 2 2-γ and it satis-

Lemma 2.1 ([46])
The solution of the following system: is given by

Theoretical approach
In this section, we aim to present the existence and uniqueness result for a solution of the model (1.1)-(1.2) by using Picard's successive iterative approximation method [47]. For this purpose, let X = × × × × denote a Banach space with supremum norm Now, we rewrite the model (1.1) in the following form: such that and According to Lemma 2.1, the system (3.3) is equivalent to the following fractional integral equation: Proof Obviously, the solution of the model (1.1)-(1.2) is equivalent to the fractional integral equation (3.6). Consider (3.9) . Our target is a proof that the sequence {ϒ n (t)} converges to ϒ(t).

Numerical approach
Throughout this section, we introduce the series type solution of the proposed model (1.1)-(1.2), by utilizing the Laplace transform with the Adomian decomposition method [48]. For the convergence of such a method, we refer the reader to [49]. Applying the Laplace transform to both sides of (1.1), we get (4.1) Next, suppose the solution in the following series type: Moreover, by Adomian's polynomial we can decompose the nonlinear terms E(t)I(t) and E(t)D(t) as follows: where the Adomian polynomial A n (E, I) can be defined as In particular, we have A 0 (E, I) = E 0 (t)I 0 (t), A 1 (E, I) = E 1 (t)I 0 (t) + E 0 (t)I 1 (t). (4.5) Similarly, we can define the polynomial B n (E, D).

Numerical simulations and discussion
Here, we present the numerical simulation for the solution of the considered model (1.1)-(1.2) in the form of an infinite series as given in (4.8). Here the time is in days. The numerical values of the parameters utilized in the simulation are designated in Table 1 From Figs. 1-5, one can say that a huge population of exposed individuals becomes infected within a month. The decline in uninfected population is represented via different fractional order by taking the first ten terms. It is faster at smaller fractional order and slower at greater order. The decay occurs in the uninfected class under different fractional orders which is faster at lower fractional value by taking the first ten terms of the approximate solution in Fig. 1. As a result the exposed class will go up at a different rate due to the fractional order derivative as given in Fig. 2. This is because more people are exposed to infection. The infected class also goes on increasing. Here the growth is slow at lower fractional order as compared to higher order as in Fig. 3. Furthermore, the critically infected cases and hospitalization cases are also increasing as in Fig. 4. From Fig. 5, the death class leads to a fluctuation, maybe due to better care of infected people who recovered from the disease.

Figure 2
Graphical representation of the approximate solution I for the first ten terms at different fractional order Figure 3 Graphical representation of the approximate solution C for the first ten terms at different fractional order Figure 4 Graphical representation of the approximate solution H for the first ten terms at different fractional order Furthermore, we investigate the dynamical behavior on increasing the values of the three parameters α i (i = 1, 2, 3), corresponding to integer order.
In Figs. 6-10, we have presented the dynamical behavior at different values of the parameters α i (i = 1, 2, 3) corresponding to integer order. As we increase the values of α i Figure 5 Graphical representation of the approximate solution D for the first ten terms at different fractional order   1, 2, 3), the corresponding infection, hospitalization and chronic infection are increasing. Also the uncertain behavior has been observed on increasing the corresponding values of the aforesaid three parameters. Here we use a nonstandard finite difference scheme for the illustration of the dynamics at given values of the parameters.

Conclusion
First of all, it is necessary that corresponding to a real-world problem the model one built would exist. This question should be guaranteed and in this regard the fixed point approach is a powerful analysis which gives proper information about the existence of such model. On the other hand, using a nonsingular derivative of fractional order for real-world problems is a new field in the last few years. Such investigations have been proved to get significant information about the global dynamics of an infectious disease. Furthermore, treating such a type of model by the Laplace Adomian decomposition method is another best way to handle an approximate solution of such a type of problems. For the Caputo-Fabrizo case, this concept has been very rarely adopted. The mentioned techniques omit discretization of data and need no collocation to control the method. Therefore, with the help of Picard's iterative methods, we have successfully established a qualitative theory for a five compartment model of COVID-19 with Caputo-Fabrizo fractional order derivatives. Furthermore, some approximate analytical results have been developed via the Laplace Adomian decomposition method. The concerned solution has been presented via graphs for some numerical values. The fractional order derivative provides some more details of the transmission dynamics of the proposed model. In the future the concerned analysis can be extended to other mathematical models of infectious diseases.