Abstract
In this paper, we consider a fractional-order mathematical system comprising four different compartments for the recent pandemic of SARS-CoV-2 with regard to global and singular kernels of Caputo fractional operator. The SARS-CoV-2 fractional mathematical model is analyzed for series-type solution by Laplace–Adomian decomposition techniques (LADM) and homotopy perturbation method (HPM). The whole quantity of each compartment is divided into small parts, and then the sum of these all parts is written as a series solution for each agent of the system, while the nonlinear part is decomposed using the Adomian polynomial. The model is also checked for approximate solution by HPM through a comparison of the parameter power, , for each equation. The numerical simulation for both methods is provided in different fractional orders along with comparison with each other as well as with natural order 1.
1. Introduction
The novel coronavirus (SARS-CoV-2), considered the most dangerous virus of this decade, belong to family of severe acute respiratory syndrome (SARS) [1, 2]. Therefore, this new virus is related to the viruses associated with the syndrome. It has become a new novel strain of the SARS family, which was recognised in humans before [3, 4]. SARS-CoV-2 affected not only humans but also several animals. This virus has been transmitted from human to human, and it occurs similarly in animals. But many times, it has become a mystery that how the animals have been affected by it in certain security. Infected humans and different species of various animals are also recognised as an active cause of the spreading of the virus [5]. In the past, some similar viruses like the Middle East respiratory syndrome coronavirus (MERS-CoV) were spread from camels to the human population, and for SARS-CoV-1, the civet cats were recognised as the source of spreading into humans [6].
Mathematical models concerning infections have vastly been used since the last century to study the dynamics and transmission of various pandemics and epidemics and to apply valuable techniques for their control or minimization. Scholars investigating pandemics in the various areas of sciences are working to control these epidemics or reduce these negative impacts to a stable situation [7–9]. They also give the concept of the globalization of the ODEs that have the natural-order differentiation providing more degrees of freedom at any order. The equations which contain the FO differential equation , with FO , may be studied in [10–12]. Mostly, in epidemiological problems, FDEs refer to models with memory effects [12]. Next, the fractional-order derivative has the term of integration providing the knowledge of the past spreading for an infection. We can expect the behavior of the transmission based on the previous results and studies. The hereditary and historical characteristics point to the past transmission of diseases, which is very beneficial for making predictions. Therefore, these characteristics can be tested by the application of noninteger order derivatives, and it impacts the transmission of an epidemic [10–12].
Modern calculus goods production arbitrary order may be used in different fields of clinical and physical sciences, such as goods production by engineering, control theory, economics, financing, and infectious disease conditions. The extensive huge study of FODEs in modelling global phenomena is because of more attracting characteristics which are not explained in the natural order derivatives. The natural order differential equations are local in nature, while FDEs are concerned with nonlocality, which provides more globalization of their dynamics. Usage of integer operator is hot area and recently caught the researchers interest, whereas the noninteger order operators have been studied and used intensively. Infectious problems for the endemic are the most realistic area for the researchers as a research gap and are applied to test them in recent times. Moreover, the analysis for the mathematical models referring to the real-globe situation is made using the theory of stability, existence of solution, and with the optimal problem as in [13–15].
We proposed a new vaccinated SARS-CoV-2 epidemic model that has four quantities including the susceptible class , the acutely infected class , chronically infected class , and the recovered class , which takes the following form:
The parameters used in system (1) are described in Table 1:
The analysis of model (1) under fractional-order derivative with regard to the Caputo operator is given aswith general initial approximation , , , and .
In the seventeenth century, many researchers Euler, L'Hôpital, Fourier, Abels, Riemann–Liouville, etc., made fundamental contributions in this field of modern calculus and were known as the pioneers of fractional calculus. Besides them, many other scientists have made significant findings and discovered some fractional models as seen in [10–12]. Most of the basic properties needed for real phenomena such as memory, globality, and hereditary are involved in various fractional operators, while the integer-order differential operator have no such properties; therefore, modern calculus gives more realistic result. Modern calculus refers to the appliances of biological, physical, and engineering as in [11, 16–18]. Other properties of fractional operators like nonlocality, globality, singularity, and nonsingularity also attract the interest of many researchers. These properties are more applicable to most real-world problems.
It should be kept in mind that fractional operators do have not a unique definition and are formulated by different formulae. Most of the operators have the definite integral with singular and nonsingular kernels. These kernels can also be found in the various fractional integration formulae. These kernels are mostly taken from the Cauchy integral formula. These integral kernels present in the fractional operator may be analyzed by various techniques. Some researchers have used Laplace–Adomian decomposition techniques for both linear and nonlinear fractional differential models. For the terms of noninearity, Adomian polynomial is applied, which decomposes that term into a small one (see [19, 20]). Similarly, the homotopy perturbation techniques are also applied for a series solution, which can be obtained by perturbed by a small factor and then comparing the power of parameter on both sides of the equation. For the iterative numerical solution of integer-order models, mostly Runge–Kutta techniques (RK4 and RK2) were used. They can be used for fractional-order equations by involving the fractional terms. In this article, we will apply Laplace transformation along with Adomian decomposition techniques by converting the while quantity into small ones. We will also apply the homotopy perturbation techniques (HPM) for the series solution of the said problem [21, 22].
Different global real-life situations may be represented by noninteger differential equations such as physics problems, problems of control theory, chemistry problems (polymerization, rheology, and electronics), physical biology, heat, aerodynamics, infectious diseases, electro-statics, electrical circuits, and blood flow [11, 23–25]. In the last few years, the interesting description of the existence of solutions of by changing the boundary conditions counting integer order phenomena, nonlocal, nonsingular kernel, periodic/antiperiodic, and multipoints has been investigated and has provided good supposition, as can be seen in [26–28]. Different researchers have tested the solution of the coupled models by differential equations of noninteger orders which gives the key role in applied fields of mathematics. This is the establishment for the model developed from biochemistry, biology, ecology, and other areas of engineering and physical sciences such as in [29].
2. Basic Definition
In this section, some definitions are written from [19, 20, 30].
Definition 1. Consider to be a mapping; then, the Riemann–Liouville noninteger order integral of order is as follows:where the right integral must exist.
Definition 2. Let be an operator, then the Caputo noninteger order derivative can be defined asand the R.H.S of the integration must exist and . If , then someone has
Lemma 1. From the noninteger order differential equation, the following is satisfied:
Definition 3. With regard to the Caputo operator, the Laplace transformation can be written as
Definition 4. On application of the homotopy perturbation techniques to an equation having linear and nonlinear classes, we may write a homotopy .where and is the embidding parameter.
3. Series Solution of Problem (2) via LADM
This section is devoted to the analysis of general-typed series techniques for the considered model (2) along with some starting conditions. On the application of the Laplace transformation to the considered problem (2), we get the following:
Applying the initial approximation, problem (9) becomes
Taking the infinite series solution for asthe nonlinear term can be decomposed as
Substitute equations (11) and (12) into equation (10), and on comparing the terms on each side, we obtain
Upon using the inverse Laplace transform in expression (14), we get
In the same way, the remaining terms can be obtained, and the terms given in equation (14) are as follows:
4. Approximate Solution for Problem (2) via HPM
Furthermore, we apply the HPM to obtain the approximate solution for the proposed problem (2), according to [20, 21] as follows:
By applying in equation: (16), a model of arbitrary differential equations can be obtained:
Solutions for the equation (17) are straightforward. If we use in equation (16), we will obtain the same system as (2). The infinite series solution can be write in the following form as
So, by comparing each term , in equation (18), we get the original model. Plugging equation (18) in to equation (16) and on comparison of each terms with the same power of , we get
Similarly,
Similarly, the high terms can be obtained, and the required terms are given in the part above. Hence, we obtain the same high terms as we obtained using LADM. Both the methods are applied as strong techniques for nonlinear fractional-order equations.
5. Discussion along with Numerical Simulation for the Proposed System (2)
Next, we compute the numerical solution for the considered system (2), by assigning different values for the parameter given in Table 1 used in problem (2).
The first four terms of the considered model (2), by using values from the table are given as follows:
The solution for the first few terms are the following:
Now, evaluating equation (23) with , we get
Similarly, for , the approximate solution for (23) isand for , the approximate solution for (23) is
For the verification of our semianalytical solution by both methods, we provide numerical simulations for problem (2). Our simulation pertains to the qualitative point analysis, and the parameters are considered through a biological feasibility approach. We take parameter numerical value and then different initial class sizes for each of the compartment, namely, susceptible , acutely infectious , chronic carries, , and recovered , from Table 2. Figure 1(a) shows a a quick decline in the starting in the quantity of susceptible class at at different fractional orders. The occurred decline took much more time at low nonnatural order and very slow at high noninteger order. We observe that, as the arbitrary order values increased, the curve of simulation goes, converging to the order 1. Besides this, after few days, the said individuals show low growth and then show converging towards stability to the equilibrium point. Figure 1(b) represents the acute infectious class, increases for some beginning days at various fractional orders. Beside that, the curve declines and stay stable through all arbitrary orders and converges to an integer order. In Figure 1(c), one may see the chronical infectious cases. decreases at the starting of few days and then moves up to the maximum value of 1.5 at different noninteger order as no treatment and precautionary measures are done. But after keeping the precautions, the infectious case declines to 0.9 and then their after show stability and convergency. From Figure 1(d) we see that the recovery class goes up to 2.7 at the starting at various fractional orders and by keeping precautions and treatment the recovered class also stabilized.
(a)
(b)
(c)
(d)
Next, we have provided the comparison of different agents of the considered model by both the methods of LADM and HPM, as given from Figures 2(a)–3(b). Figure 4 represents the comparison of various agents at .
(a)
(b)
(a)
(b)
6. Conclusions
The current investigation is the development of the four compartmental fractional-order SARS-CoV-2 model with regard to the Caputo fractional-order operative of Caputo having a singular kernel. The analysis for the series type solution of the proposed problem has been successfully achieved by two methods, one is the Laplace transform along with the Adomian polynomial (LADM) for a nonlinear term and the other one is the homotopy perturbation method (HPM). The semianalytical type solution has been obtained by both the methods which are comparable with each other. The numerical simulation for a few terms has been plotted using the available data given in Table 1 for four different values of . We also compare few fractional-order values with that of integer 1, and as increasing the fractional-order value, we achieved the behavior of order 1. A complete spectrum for all compartments has established, and we can use other fractional values between 0 and 1. As a result, we say that fractional-order analysis provides better results than those of the integer-order problems.
Data Availability
Data are available upon the request by email and according to the type of collaboration.
Conflicts of Interest
The author declares no conflicts of interest.
Acknowledgments
The author expresses his appreciation to the Deanship of Scientific Research at King Saud University for funding this work through research group no. RG-1437-017.