Next Article in Journal
On the Stability Domain of a Class of Linear Systems of Fractional Order
Previous Article in Journal
Fourth-Order Numerical Solutions for a Fuzzy Time-Fractional Convection–Diffusion Equation under Caputo Generalized Hukuhara Derivative
Previous Article in Special Issue
Investigation of Novel Piecewise Fractional Mathematical Model for COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Modeling of COVID-19 Transmission Using a Fractional Order Derivative

by
Badr S. Alkahtani
Department of Mathematics, College of Science, King Saud University, P.O. Box 1142, Riyadh 11989, Saudi Arabia
Fractal Fract. 2023, 7(1), 46; https://doi.org/10.3390/fractalfract7010046
Submission received: 20 October 2022 / Revised: 12 December 2022 / Accepted: 15 December 2022 / Published: 30 December 2022
(This article belongs to the Special Issue Fractal Approaches in Materials: Structure and Mechanics)

Abstract

:
In this article, the mathematical model of COVID-19 is analyzed in the sense of a fractional order Caputo operator with the consideration of an asymptomatic class. The suggested model is comprised of four compartments. The results from fixed point theory are used to theoretically analyze the existence and uniqueness of solution of the model in fractional perspective. For the numerical approximation of the suggested problem, a numerical iterative scheme is used, which is based on the Newton polynomial interpolation. For the efficiency and applicability of the suggested technique with a fractional Caputo operator, we simulate the results for various fractional orders.

1. Introduction

Wuhan, China, was the first to be infected with the new virus (2019-nCoV), which is extremely pathogenic and transmissible. This unique illness has spread over the world, causing severe acute respiratory syndrome.The COVID-19 epidemic has attracted much attention, since the first case of COVID-19 was detected. Unfortunately, the globe has reached a gloomy coronavirus milestone as of 20 December 2021, with 5.36 million fatalities recorded and confirmed cases of about 275.6 million. The total deaths have far outpaced that of the other two coronaviruses (severe acute respiratory syndrome coronavirus, SARS-CoV, and Middle East respiratory syndrome coronavirus, MERS-CoV), according to [1,2]. The worldwide epidemic is ongoing, wreaking havoc on public health and the economy throughout the world. Despite the country’s collective efforts, China’s COVID-19 pandemic has not been effectively contained. The influence of international imports and population movement (particularly among subclinical patients) on contagious control, on other hand, cannot be overlooked. It is critical to build a model for the contagion as a potent instrument to study the mechanisms of contagious disease control and transmission, in order to control contagious diseases and limit their incidence. When COVID-19 was first discovered, researchers employed epizootic data or dynamic models, as well as the effective times of COVID-19’s reproduction to predict the epidemic’s peak timing and size both locally and internationally [3,4]. Researchers also created forecast models on the impact of resuming work on the progression of the disease in Hubei Province, if not the entire country [5]. Asamoah et al. studied the global stability and cost-effectiveness of COVID-19, with the consideration of the environmental impact using Ghana’s data [6].
The early research lacked adequate raw data. As a result, the majority of these researched forecasts of the pandemic scenario differed from the actual reality.
More crucially, these investigations ignored the novel coronavirus’s high conveyance capability throughout the incubation phase, subclinical infection, and the effect of population movement on the epidemic conveyance. After extensive study on the new coronavirus pneumonia, it was oberved that social distancing had a strong impact on decreasing the number of infections. For this, several works have been presented; for instance, Qian et al. worked on COVID-19 and social distancing [7], Mwalili studied an SEIR COVID-19 model by incorporating social distancing and environment [8], Elgazzar presented the controling of COVID-19 through social distancing and awareness [9].
Most academics included the conveyance characteristics of individuals with a latent period or subclinical diseases into new coronavirus pneumonia models in later investigations [10]. Mandal et al. [11] developed a model that included a quarantine class and government action. This study found that the most important component in attaining disease management was minimizing the exposure to exposed and vulnerable individuals. Khan and Atangana et al. [12] devised a new mathematical model for COVID-19 dynamics in the presence of quarantine and isolation. Some researchers have proposed mathematical models to investigate COVID-19’s proliferation and conveyance in the community, particularly the role of asymptomatic infected persons. The authors of [13] suggested the COVID-19 spread using a constituent mathematical model in order to concentrate on super-spreaders’ contagiousness. However, only the infectiousness of exposed people was included in this model, not the infectiousness of asymptomatic people with illnesses. Indeed, research employing data from early Chinese reports paired with Bayesian inference methodology found that asymptomatic illnesses increased the epidemic’s spread [14]. Furthermore, a great number of mathematical models or studies focusing on COVID-19 and other relevant subjects have been established [15]. Yet, it is uncommon to find a model that takes into account the influence of both asymptomatic infected people’s illness features and the population’s mobility on COVID-19 conveyance.
Kang et al. [16] provided a simplified mathematical model, which could be utilized to observe the influence of the mobility of the population and subclinical infected persons on the development of COVID-19. As opposed to the other models described above, they took into account the following assumptions. To begin with, because symptomatic infected people have apparent ailments, they will be segregated and treated as soon as they are discovered. This indicates that the rate of viral infection in susceptible persons from symptomatic diseases is relatively low. As a result, the infection rate from symptomatic illnesses to vulnerable persons was overlooked. Second, since May, the pandemic situation in the provinces has mostly stabilized. This also suggests that in China, both the medical issues and the treatment methods have remained relatively steady. As a result, our model incorporated the constant cure rate. However, when people return to regular life, the chance of an epidemic re-emerging cannot be overlooked due to the large population, particularly asymptomatic persons. They then concentrated on the influence of population movement on the epidemic’s progression. Furthermore, recurrence incidences in people who have been cured are quite rare. As a result, they believed that there was no transfer of recovery to the vulnerable population.
In assessing the dynamics of an infectious disease, epidemic models of fractional order are more informative and reliable than typical integer order models, and they are being used more frequently [17,18]. For some diseases, fractional order models exhibit a substantially superior match to the real data. Many fractional operators were suggested in [19,20]; these fractional operators’ applications were presented in [21,22,23]. The Caputo operator has been used by several researchers to analyze a variety of real world phenomena. For instance, Saifullah et al. investigated a nonlinear wave model in [24], Khan et al. studied a four dimensional dynamical system with Caputo’s operator [25]. Similarly, Alqahtani et al. analyzed the bioethanol production model under a generalized nonlocal operator in the Caputo sense [26]. Furthermore, in [27], for instance, they investigated equations of fractional diffusion and their analysis. The authors of [28] employed a novel approach for numerically solving fractional order differential equations. In [29], the authors studied the nonlinear wave equation with solitary/shock solutions. In [30], the author investigated coronavirus modeling, simulations, and potential control using a mathematical model. The authors recently investigated the fractional derivative analysis of the coronavirus model [31]. The authors of [32] investigated a mathematical model to study COVID-19 and its control analysis. Several authors have conducted work on fractional differential equations [33,34,35].
We reformulated the model [16] in this paper using the Caputo derivative with a nonlinear incidence rate and fixed input and fixed treatment rates. Fractional order models provide a deeper grasp of the epidemic and provide additional insights. The remainder of the paper is organized as follows: Section 2 presents the basic preliminaries, and Section 3 presents the model formulation in terms of the parameter estimates and curve fitting. Section 4 gives the model derivation. The model’s analysis is presented in Section 5, whereas numerical simulations are shown in Section 6. Section 7 has brief concluding remarks.

2. Preliminaries

Definition 1. 
We consider U ( t ) to be a continuous function; then, the definition of a Caputo operator of order α, where m 1 < α 1 is
0 C D t α U ( t ) = 1 Γ ( m α ) 0 t ( t τ ) m α 1 [ U ( τ ) ] d τ .
Definition 2. 
We consider U ( t ) to be a continuous function; then, the definition of the Riemann–Liouville integration with respect to t is
0 R I t α U ( t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 U ( ø ) d τ , α > 0 ,
with a converging integral.

3. Classical Integer Order Model Formulation

Here, we examine the SAIR mathematical model, which models COVID-19’s dynamic conveyance in mobile populations. Susceptible individuals S ( t ) , asymptomatic people A ( t ) , symptomatic individuals I ( t ) , and recovered individuals R ( t ) make up the population size N ( t ) . The following is a system of ordinary differential equations [16].
d S ( t ) d ( t ) = Π A ( t ) S ( t ) N ( t ) υ S ( t ) , d A ( t ) d ( t ) = A ( t ) S ( t ) N ( t ) υ A ( t ) q 1 A ( t ) q 2 A ( t ) , d I ( t ) d ( t ) = q 1 A ( t ) υ I ( t ) κ I ( t ) ν I ( t ) d R ( t ) d ( t ) = q 2 A ( t ) + ν I ( t ) υ R ( t ) .
The starting state was S ( 0 ) > 0 , A ( 0 ) 0 , I ( 0 ) 0 , and R ( 0 ) 0 . The comprehensive input rate was represented by the parameter Π > 0 , and υ > 0 represented the natural death rate. The death rate owing to disease was denoted by κ > 0 . Vertical conveyance was not taken into account in this model; therefore, all babies were at risk. Asymptomatic infections, according to clinical practice, are not in need of therapy; yet, this population is a main infection source. The rate of asymptomatic people who have recovered becoming symptomatic was determined by the parameter q 1 . ⋏ represented the transmission rate between asymptomatic and susceptible people. The rate at which asymptomatic people recover was determined by parameter q 2 > 0 . The rate at which symptomatic people recover was ν > 0 . In clinical terms, the recurrence in treated patients is uncommon. As a result, from recovered people, no conveyance to the vulnerable class was envisaged in this model. The authors in [16] studied the local and global stability based on the disease-free as well as the endemic equilibrium points by constructing an appropriate Lyapunov function under appropriate parameter conditions. Moreover, some important results have been presented and verified by numerical simulations.

4. Model Derivation in the Caputo Operator

We used a Caputo fractional derivative to reformulate the COVID-19 model (1) to study the memory effects and learn more about the epidemic. The following is a system of Caputo fractional differential equations. Moreover, to obtain the same dimension on both sides of the coefficient, we included one auxiliary parameter Λ [36,37,38]
1 Λ 1 α C D t α S ( t ) = Π A ( t ) S ( t ) N ( t ) υ S ( t ) , 1 Λ 1 α C D t α A ( t ) = A ( t ) S ( t ) N ( t ) υ A ( t ) q 1 A ( t ) q 2 A ( t ) , 1 Λ 1 α C D t α I ( t ) = q 1 A ( t ) υ I ( t ) κ I ( t ) ν I ( t ) 1 Λ 1 α C D t α R ( t ) = q 2 A ( t ) + ν I ( t ) υ R ( t ) ,
with S ( 0 ) > 0 , A ( 0 ) 0 , I ( 0 ) 0 , and R ( 0 ) 0 .

5. Analysis of the Model

5.1. Existence Results

Here, the existence and uniqueness of the solution of Model (2) are presented using some fixed point results. For this purpose, Model (2) can be written as
1 Λ 1 α C D t α S ( t ) = G 1 ( t , S ) , 1 Λ 1 α C D t α A ( t ) = G 2 ( t , A ) , 1 Λ 1 α C D t α I ( t ) = G 3 ( t , I ) , 1 Λ 1 α C D t α R ( t ) = G 4 ( t , R ) .
System (3) may be expressed as
1 Λ 1 α C D t α F ( t ) = H ( t , F ( t ) ) F ( 0 ) = F 0 ,
where
F ( t ) = S ( t ) , A ( t ) , I ( t ) , R ( t ) . H ( t , F ( t ) ) = G 1 , G 2 , G 3 , G 4 .
Now, to define an orbitally complete metric space, let ( H , d ) represent the metric space and T : H H ; if x 0 H , then the orbit of x 0 is the set
O ( x 0 ) = { T n x 0 : n = 0 , 1 , 2 , 3 } ,
where T n is the nth iteration of T , and D ( x 0 ) is the diameter of O ( x 0 ) . If all Cauchy sequences belonging to O ( x ) converge in H for some x H , then ( H , d ) is known as T -orbitally complete metric space.
Theorem 1 
([39]). Let ( H , d ) represent T -orbitally complete metric space, T : H H , and θ : H N ; if ν > 0 and x 0 H , with 0 < D < , such as
d ( T θ ( x ) ( p ) , T θ ( x ) ( q ) ) e ν d ( p , q ) ,
then T has a unique fixed point.
Let H = C × C × C , where C [ 0 , T ] , R represents space containing continuous mappings, and d ( p , q ) = sup t [ 0 , T ] | p q | is the metric on H .
The existence and uniqueness of the solution of the problem (4) are presented and proved in the following theorem.
Theorem 2. 
Let T : H H be defined by T ( F ( t ) ) = F 0 + 1 α Γ ( α ) 0 t ( t τ ) α 1 H ( τ , F ( τ ) ) d τ and
| H ( t , F ( t ) ) H ( t , K ( t ) ) | Γ ( α + 1 ) α T ^ α e ν | | F | | K | |
| H ( t , F ( t ) ) H ( t , K ( t ) ) | Γ ( α + 1 ) α T ^ α e ν | | F | + | K | | .
Then, there must be a unique solution to problem (4).
Proof. 
By applying the integral on Equation (4), we obtain
F ( t ) = F 0 + 1 α Γ ( α ) 0 t ( t τ ) α 1 H ( τ , F ( t ) ) d τ = T F .
To prove that there exists a unique fixed point of T ,
| T F T K | = 1 α Γ ( α ) 0 t ( t τ ) α 1 ( H ( τ , F ) ) d τ 1 α Γ ( α ) 0 t ( t τ ) α 1 | H ( τ , F ) H ( τ , K ) ) | d τ Γ ( α + 1 ) e ν Γ ( α ) T ^ α 0 t ( t τ ) α 1 | | F | | K | | d τ = α e ν T ^ α sup t [ 0 , T ] | | F | | K | | 0 t ( t τ ) α 1 d τ e ν sup t [ 0 , T ] | | F | | K | | .
Moreover,
| T F | + | T K | = 1 α Γ ( α ) 0 t ( t τ ) α 1 ( H ( τ , F ) ) d τ | + | 1 Γ ( α ) 0 t ( t τ ) α 1 ( H ( τ , K ) ) d τ 1 α Γ ( α ) 0 t ( t τ ) α 1 ( | H ( τ , F ) | + | H ( τ , K ) | ) d τ Γ ( α + 1 ) e ν Γ ( α ) T ^ α 0 t ( t τ ) α 1 | | F | + | K | | d τ = α e ν T ^ α sup t [ 0 , T ] | | F | + | K | | 0 t ( t τ ) α 1 d τ e ν sup t [ 0 , T ] | | F | + | K | | sup t [ 0 , T ] | | F | + | K | | .
Now,
d ( T 2 F , T 2 K ) = sup t [ 0 , T ] | T 2 F T 2 K | = sup t [ 0 , T ] | T F T K | × sup t [ 0 , T ] | T F + T K | sup t [ 0 , T ] | T F T K | × sup t [ 0 , T ] ( | T F | + | T K | ) e ν sup t [ 0 , T ] | | F | | K | | × sup sup t [ 0 , T ] | | F | + | K | | = e ν sup t [ 0 , T ] | | F | | K | | e ν sup t [ 0 , T ] | F K | = e ν d ( F , K ) .
If we take θ : H N , such that θ ( F ) = 2 for each F H , then all conditions of Theorem 1 hold; therefore, there exists a fixed point of T , which is unique, and consequently, problem (4) has a unique solution. □

5.2. Ulam–Hyers Stability

Here, we present the U-H and the generalized U-H stability [40,41] of the suggested system. Let us assume ε , with the following inequality
| C D p Ʊ ˜ ( t ) X ( t , Ʊ ˜ ( t ) ) | ε , t J , ε = max ( ε i ) T , i = 1 , 2 , 3 , 4 .
Definition 3. 
System (2) is UH stable, if there exist U X > 0 for all ε > 0 ; then the solution of Ʊ ˜ Y holds for (6); there is a unique solution Ʊ Y for Equation (4), such that
| Ʊ ˜ ( t ) Ʊ ( t ) | U X ε , t J , U X = max ( U X j ) T .
Definition 4. 
System (2) is generalized UH stable, if there exists a continuous function Φ : R + R + and Φ ( 0 ) = 0 , so that for all solutions Ʊ ˜ Y of (6), there is a unique solution Ʊ Y for (4), with the following
| Ʊ ˜ ( t ) Ʊ ( t ) | Φ X ε , t J , Φ X = max ( Φ X j ) T .
Remark 1. 
A function Ʊ ˜ Y satisfies (6), if and only if there exists a function Φ Y with the following properties:
(I
| Φ ( t ) | ε , Φ = max ( Φ j ) , t J .
(II
C D p Ʊ ˜ ( t ) = X ( t , Ʊ ˜ ( t ) ) + Φ ( t ) , t J .
Lemma 1. 
If Ʊ ˜ Y holds for Equation (6), then Ʊ ˜ also holds for the following
| Ʊ ˜ ( t ) Ʊ ˜ 0 ( t ) 1 α Γ ( α ) 0 t ( t ϑ ) α 1 X ( ϑ , Ʊ ˜ ( ϑ ) ) d ϑ | ε .
Proof. 
Using (II), we have
C D p Ʊ ˜ ( t ) = X ( t , Ʊ ˜ ( t ) ) + Φ ( t ) ,
and along with Lemma 1, we obtain
Ʊ ˜ ( t ) = Ʊ ˜ 0 ( t ) + 1 α Γ ( α ) 0 t ( t ϑ ) α 1 X ( ϑ , Ʊ ˜ ( ϑ ) ) d ϑ + 1 α Γ ( α ) 0 t ( t ϑ ) α 1 Φ ( ϑ ) d ϑ .
Next, using (I) gives
| Ʊ ˜ ( t ) Ʊ ˜ 0 ( t ) 1 α Γ ( α ) 0 t ( t ϑ ) α 1 X ( ϑ , Ʊ ˜ ( ϑ ) ) d ϑ | 1 α Γ ( α ) 0 t ( t ϑ ) α 1 | Φ ( ϑ ) | d ϑ ε .
Therefore, the proof is finished. □
Theorem 3. 
For all Ʊ Y , X : J × R 4 R with the assumption ( C 1 ) holds, and 1 p L X > 0 . Equation (4) is equal to Equation (2) and is UH stable and, consequently, generalized UH stable.
Proof. 
Suppose that Ʊ , Ʊ ˜ Y is a unique solution of (4); therefore, for all ε > 0 , t J , along with Lemma 1, we have
| Ʊ ˜ ( t ) Ʊ ( t ) | = max t J | Ʊ ˜ ( t ) Ʊ 0 1 α Γ ( α ) 0 t ( t ϑ ) α 1 X ( ϑ , Ʊ ( ϑ ) ) d ϑ | max t J | Ʊ ˜ ( t ) Ʊ ˜ 0 1 α Γ ( α ) 0 t ( t ϑ ) α 1 X ( ϑ , Ʊ ( ϑ ) ) d ϑ | + max t J 1 α Γ ( α ) 0 t ( t ϑ ) α 1 | X ( ϑ , Ʊ ˜ ( ϑ ) ) X ( ϑ , Ʊ ( ϑ ) ) | d ϑ | Ʊ ( t ) Ʊ ˜ 0 1 α Γ ( α ) 0 t ( t ϑ ) α 1 X ( ϑ , Ʊ ˜ ( ϑ ) ) d ϑ | + L X Γ ( α ) 0 t ( t ϑ ) α 1 | Ʊ ˜ ( ϑ ) Ʊ ( ϑ ) | d ϑ α ε + α L X | Ʊ ˜ ( t ) Ʊ ( t ) | ,
which gives
Ʊ ˜ Ʊ U X ε .
From (10), we may write
U X = α 1 α L X .
Hence, equating Φ X ( ε ) = U X ε , so that Φ X ( 0 ) = 0 , we conclude that the solution of (2) is stable for both the UH and the generalized UH. □

5.3. Numerical Algorithm

Next, we establish a numerical scheme for the proposed model (2). We consider the first equation of System (3)
1 Λ 1 α C D t α S ( t ) = G 1 ( t , S ) .
Equation (12) can also be rearranged in the following way:
S ( t ) S ( 0 ) = 1 α Γ ( α ) 0 t G 1 ( τ , S ( τ ) ) ( t τ ) α 1 d τ .
We can write the following at point t ı + 1 = ( ı + 1 ) Δ t :
S ( t ı + 1 ) S ( 0 ) = 1 α Γ ( α ) 0 t ı + 1 G 1 ( τ , S ( τ ) ) ( t ı + 1 τ ) α 1 d τ .
As a result,
S ( t ı + 1 ) = S ( 0 ) + 1 α Γ ( α ) = 2 ı t t + 1 G 1 ( τ , S ( τ ) ) ( t ı + 1 τ ) α 1 d τ .
By substituting Equation (15) for the Newton polynomial, we obtain
S ı + 1 = S 0 + 1 α Γ ( α ) = 2 ı t t + 1 G 1 ( t 2 , S 2 ) + G 1 ( t 1 , S 1 ) G 1 ( t 2 , S 2 ) Δ t ( τ t 2 ) + G 1 ( t , S ) 2 G 1 ( t 1 , S 1 ) ( G 1 ( t 2 , S 2 ) ) 2 ( Δ t ) 2 × ( τ t 2 ) ( τ t 1 ) × ( t ı + 1 τ ) α 1 d τ .
As a result, the preceding equation may be rearranged in the following manner:
S ı + 1 = S 0 + 1 α Γ ( α ) = 2 ı t t + 1 G 1 ( t 2 , S 2 ) ( t ı + 1 τ ) α 1 d τ + t t + 1 G 1 ( t 1 , S 1 ) G 1 ( t 2 , S 2 ) Δ t ( τ t 2 ) ( t ı + 1 τ ) α 1 d τ + t t + 1 G 1 ( t , S ) 2 G 1 ( t 1 , S 1 ) ( G 1 ( t 2 , S 2 ) ) 2 ( Δ t ) 2 ( τ t 2 ) ( τ t 1 ) ( t ı + 1 τ ) α 1 d τ .
Consequently,
S ı + 1 = S 0 + 1 α Γ ( α ) = 2 ı G 1 ( t 2 , S 2 ) t t + 1 ( t ı + 1 τ ) α 1 d τ + 1 Γ ( α ) = 2 ı G 1 ( t 1 , S 1 ) G 1 ( t 2 , S 2 ) Δ t × t t + 1 ( τ t 2 ) ( t ı + 1 τ ) α 1 d τ + 1 Γ ( α ) = 2 ı G 1 ( t , S ) 2 G 1 ( t 1 , S 1 ) + ( G 1 ( t 2 , S 2 ) ) 2 ( Δ t ) 2 × t t + 1 ( τ t 2 ) ( τ t 1 ) ( t ı + 1 τ ) α 1 d τ .
In Equation (18), we may compute the aforementioned integrals as follows:
t t + 1 ( t ı + 1 τ ) α 1 d τ = ( Δ t ) α α ( ı + 1 ) α ( ı ) α , t t + 1 ( τ t 2 ) ( t ı + 1 τ ) α 1 d τ = ( Δ t ) α + 1 α ( α + 1 ) × ( ı + 1 ) α ( ı + 3 + 2 α ) ( ı ) α ( ı + 3 + 3 α ) ,
t t + 1 ( τ t 2 ) ( τ t 1 ) ( t ı + 1 τ ) α 1 d τ = ( Δ t ) α + 2 α ( α + 1 ) ( α + 2 ) × ( ı + 1 ) α 2 ( ı ) 2 + ( 3 α + 10 ) ( ı ) + 2 α 2 + 9 α + 12 ( ı ) α 2 ( ı ) 2 + ( 5 α + 10 ) ( ı ) + 6 α 2 + 18 α + 12 .
We may see the following strategy, if we insert these values into Equation (18):
S ı + 1 = S 0 + ( Δ t ) α Γ ( α + 1 ) = 2 ı G 1 ( t 2 , S 2 ) × ( ı + 1 ) α ( ı ) α + ( Δ t ) α Γ ( α + 2 ) = 2 ı G 1 ( t 1 , S 1 ) G 1 ( t 2 , S 2 ) × ( ı + 1 ) α ( ı + 3 + 2 α ) ( ı ) α ( ı + 3 + 3 α ) × + ( Δ t ) α 2 Γ ( α + 3 ) = 2 ı G 1 ( t , S ) 2 G 1 ( t 1 , S 1 ) + ( G 1 ( t 2 , S 2 ) ) × ( ı + 1 ) α 2 ( ı ) 2 + ( 3 α + 10 ) ( ı ) + 2 α 2 + 9 α + 12 ( ı ) α 2 ( ı ) 2 + ( 5 α + 10 ) ( ı ) + 6 α 2 + 18 α + 12 .
Similarly from the second, third, and fourth equations of system (12), we can write
A ı + 1 = A 0 + ( Δ t ) α Γ ( α + 1 ) = 2 ı G 2 t 2 , A 2 ( ı + 1 ) α ( ı ) α + ( Δ t ) α Γ ( α + 2 ) = 2 ı G 2 t 1 , A 1 G 2 t 2 , A 2 × ( ı + 1 ) α ( ı + 3 + 2 α ) ( ı ) α ( ı + 3 + 3 α ) + ( Δ t ) α 2 Γ ( α + 3 ) = 2 ı G 2 t , A 2 G 2 t 1 , A 1 + G 2 t 2 , A 2 × ( ı + 1 ) α 2 ( ı ) 2 + ( 3 α + 10 ) ( ı ) + 2 α 2 + 9 α + 12 ( ı ) α 2 ( ı ) 2 + ( 5 α + 10 ) ( ı ) + 6 α 2 + 18 α + 12 ,
and
I ı + 1 = I 0 + ( Δ t ) α Γ ( α + 1 ) = 2 ı G 3 t 2 , I 2 ( ı + 1 ) α ( ı ) α + ( Δ t ) α Γ ( α + 2 ) = 2 ı G 3 t 1 , I 1 G 3 t 2 , I 2 × ( ı + 1 ) α ( ı + 3 + 2 α ) ( ı ) α ( ı + 3 + 3 α ) + ( Δ t ) α 2 Γ ( α + 3 ) = 2 ı G 3 t , I 2 G 3 t 1 , I 1 + G 3 t 2 , I 2 × ( ı + 1 ) α 2 ( ı ) 2 + ( 3 α + 10 ) ( ı ) + 2 α 2 + 9 α + 12 ( ı ) α 2 ( ı ) 2 + ( 5 α + 10 ) ( ı ) + 6 α 2 + 18 α + 12 .
Moreover,
R ı + 1 = R 0 + ( Δ t ) α Γ ( α + 1 ) = 2 ı G 4 t 2 , R 2 ( ı + 1 ) α ( ı ) α + ( Δ t ) α Γ ( α + 2 ) = 2 ı G 4 t 1 , R 1 G 4 t 2 , R 2 × ( ı + 1 ) α ( ı + 3 + 2 α ) ( ı ) α ( ı + 3 + 3 α ) + ( Δ t ) α 2 Γ ( α + 3 ) = 2 ı G 4 t , R 2 G 4 t 1 , R 1 + G 4 t 2 , R 2 × ( ı + 1 ) α 2 ( ı ) 2 + ( 3 α + 10 ) ( ı ) + 2 α 2 + 9 α + 12 ( ı ) α 2 ( ı ) 2 + ( 5 α + 10 ) ( ı ) + 6 α 2 + 18 α + 12 .
Hence (20), (21), (22), and (23) are the required numerical solutions of the proposed Model (2).

6. Simulation and Discussion

The purpose of this section is to provide the simulations of the results obtained in the above section via the Newton polynomial interpolation scheme. The values of the parameters present in the suggested system were considered in two different sets, as presented in Table 1. The parameters in set 1 were used for Figure 1, Figure 2 and Figure 3, while the parameter values of set 2 were used for Figure 4, Figure 5 and Figure 6. The initial values were considered in three different sets. The initial conditions in set 1 were S = 8000 , A = 1400 , I = 2400 , and R = 2400 ; in set 2, they were S = 10000 , A = 200 , I = 1200 , and R = 1200 , and in set 3, the initial conditions were considered as S = 60,000, A = 800 , I = 1800 , and R = 1800 . In the figures, the initial values of set 1 were used in Figure 1 and Figure 4, that of set 2 were used in Figure 2 and Figure 5, and finally the values of set 3 were used in Figure 3 and Figure 6.
In Figure 1a, Figure 2a, Figure 3a, Figure 4a, Figure 5a, Figure 6a, the population behavior of the susceptible population is presented. We see in Figure 1a, Figure 2a, Figure 3a that when the value of ⋏ was large, there were oscillations in the behavior, which moved towards stability when t = 60 . Further, Figure 1b, Figure 2b, Figure 3b, Figure 4b, Figure 5b, Figure 6b show the dynamics of the asymptomatic individuals, where it can be seen that the asymptomatic population increased as the time passed, while the individuals in this class became constant as t = 70 at lower fractional orders. Similarly, Figure 1c, Figure 2c, Figure 3c, Figure 4c, Figure 5c, Figure 6c and Figure 1d, Figure 2d, Figure 3d, Figure 4d, Figure 5d, Figure 6d show the population dynamics in the infected and recovered individuals, respectively. From Figure 1c and Figure 2c, we observed that the infected population increased at the beginning, which showed a decrease as the time passed, while a fast decrease was observed in the infected population in Figure 3c, when the fractional order was 0.96 . Further, the recovered population kept increasing and then became stable at t = 60 , it was also observed that the recovered population was very large when compared to the those infected with the virus.
Figure 4, Figure 5 and Figure 6 are projected to show the dynamics in various state variables where the basic reproductive number was less than zero or where the disease died out. For this purpose, we considered the parameter = 0.3 . Here, the fractional orders were considered to be ( b l u e , 0.99 ) , ( r e d , 0.98 ) , ( b l a c k , 0.97 ) , and ( g r e e n , 0.96 ) . From all of the following, we see that the state variables became stable at lower fractional orders as compared to the high orders.

7. Conclusions

In this article, we studied the fractional order COVID-19 ( SAIR ) model in the sense of a Caputo operator. This model had a nonlinear incidence rate, constant input rate, and constant treatment rate. The existence and uniqueness of the associated solution was studied through the tools from fixed point theory. Numerically, the solution of the model was approximated with the Newton Polynomial interpolation scheme. Simulations of the results were presented, where it was observed that when the value of ⋏ was large, there were oscillations in the behaviors of the state variables, which moved toward stability faster at lower fractional orders. Similarly, the recovered population was observed to be increasing as time passed and then became stable at t = 60 ; we also observed that the recovered population was very large as compared to those who were infected with the virus.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The author extend his appreciation to research supporting project number (RSPD2023R526), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Poutanen, S.M.; Low, D.E.; Henry, B.; Finkelstein, S.; Rose, D.; Green, K.; Tellier, R.; Draker, R.; Adachi, D.; Ayers, M.; et al. Identification of severe acute respiratory syndrome in Canada. N. Engl. J. Med. 2003, 348, 1995–2005. [Google Scholar] [CrossRef] [PubMed]
  2. De Groot, R.J.; Baker, S.C.; Baric, R.S.; Brown, C.S.; Drosten, C.; Enjuanes, L.; Fouchier, R.A.M.; Galiano, M.; Gorbalenya, A.E.; Memish, Z.A.; et al. Commentary: Middle east respiratory syndrome coronavirus (mers-cov): Announcement of the coronavirus study group. J. Virol. 2013, 87, 7790–7792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Read, J.M.; Bridgen, J.R.; Cummings, D.A.; Ho, A.; Jewell, C.P. Novel coronavirus 2019-nCoV: Early estimation of epidemiological parameters and epidemic predictions. medRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  4. Shen, M.; Peng, Z.; Xiao, Y.; Zhang, L. Modeling the epidemic trend of the 2019 novel coronavirus outbreak in China. Innovation 2020, 1, 100048. [Google Scholar] [CrossRef]
  5. Yu, Z.; Zhang, G.; Liu, Q.; Lv, Q. The outbreak assessment and prediction of 2019-nCoV based on time-varying SIR model. J. Univ. Electron. Sci. Technol. China 2020, 49, 357–361. [Google Scholar]
  6. Asamoah, J.K.K.; Owusu, M.A.; Jin, Z.; Oduro, F.T.; Abidemi, A.; Gyasi, E.O. Global stability and cost-effectiveness analysis of COVID-19 considering the impact of the environment: Using data from Ghana. Chaos Solitons Fractals 2020, 140, 110103. [Google Scholar] [CrossRef]
  7. Qian, M.; Jiang, J. COVID-19 and social distancing. J. Public Health 2022, 30, 259–261. [Google Scholar] [CrossRef]
  8. Mwalili, S.; Kimathi, M.; Ojiambo, V.; Gathungu, D.; Mbogo, R. SEIR model for COVID-19 dynamics incorporating the environment and social distancing. BMC Res. Notes 2020, 13, 352. [Google Scholar] [CrossRef]
  9. Elgazzar, A.S. Simple mathematical models for controlling COVID-19 transmission through social distancing and community awareness. Z. Naturforschung C 2021, 76, 393–400. [Google Scholar] [CrossRef]
  10. Gao, W.; Veeresha, P.; Prakasha, D.G.; Baskonus, H.M. Novel dynamic structures of 2019-nCoV with nonlocal operator via powerful computational technique. Biology 2020, 9, 107. [Google Scholar] [CrossRef]
  11. Mandal, M.; Jana, S.; Nandi, S.K.; Khatua, A.; Adak, S.; Kar, T.K. A model based study on the dynamics of COVID-19: Prediction and control. Chaos Solitons Fractals 2020, 136, 109889. [Google Scholar] [CrossRef] [PubMed]
  12. Khan, M.A.; Atangana, A.; Alzahrani, E. The dynamics of COVID-19 with quarantined and isolation. Adv. Differ. Equ. 2020, 2020, 425. [Google Scholar] [CrossRef] [PubMed]
  13. Ndaïrou, F.; Area, I.; Nieto, J.J.; Torres, D.F.M. Mathematical modeling of COVID-19 transmission dynamics with a case study of Wuhan. Chaos Solitons Fractals 2020, 135, 109846. [Google Scholar] [CrossRef] [PubMed]
  14. Li, R.; Pei, S.; Chen, B.; Song, Y.; Zhang, T.; Yang, W.; Shaman, J. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science 2020, 368, 489–493. [Google Scholar] [CrossRef] [Green Version]
  15. Cai, J.; Sun, W.; Huang, J.; Gamber, M.; Wu, J.; He, G. Indirect virus transmission in cluster of COVID-19 cases, Wenzhou, China, 2020. Emerg. Infect. Dis. 2020, 26, 1343. [Google Scholar] [CrossRef]
  16. Kang, X.; Hu, Y.; Liu, Z.; Sarwar, S. Forecast and evaluation of asymptomatic COVID-19 patients spreading in China. Results Phys. 2022, 34, 105195. [Google Scholar] [CrossRef]
  17. Xu, C.; Saifullah, S.; Ali, A. Theoretical and numerical aspects of Rubella disease model involving fractal fractional exponential decay kernel. Results Phys. 2022, 34, 105287. [Google Scholar] [CrossRef]
  18. Haidong, Q.; ur Rahman, M.; Arfan, M. Fractional model of smoking with relapse and harmonic mean type incidence rate under Caputo operator. J. Appl. Math. Comput. 2022, 1–18. [Google Scholar] [CrossRef]
  19. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. arXiv 2016, arXiv:1602.03408. [Google Scholar] [CrossRef] [Green Version]
  20. Podlubny, I. An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. Math. Sci. Eng. 1999, 198, 340. [Google Scholar]
  21. Zhang, L.; ur Rahman, M.; Arfan, M.; Ali, A. Investigation of mathematical model of transmission co-infection TB in HIV community with a non-singular kernel. Results Phys. 2021, 28, 104559. [Google Scholar] [CrossRef]
  22. Shen, W.-Y.; Chu, Y.-M.; ur Rahman, M.; Mahariq, I.; Zeb, A. Mathematical analysis of HBV and HCV co-infection model under nonsingular fractional order derivative. Results Phys. 2021, 28, 104582. [Google Scholar] [CrossRef]
  23. Xu, C.; ur Rahman, M.; Baleanu, D. On fractional-order symmetric oscillator with offset-boosting control. Nonlinear Anal. Model. Control 2022, 27, 994–1008. [Google Scholar] [CrossRef]
  24. Saifullah, S.; Ali, A.; Khan, Z.A. Analysis of nonlinear time-fractional Klein-Gordon equation with power law kernel. AIMS Math. 2022, 7, 5275–5290. [Google Scholar] [CrossRef]
  25. Khan, Z.A.; Khan, J.; Saifullah, S.; Ali, A. Dynamics of Hidden Attractors in Four-Dimensional Dynamical Systems with Power Law. J. Funct. Spaces 2022, 2022, 3675076. [Google Scholar] [CrossRef]
  26. Alqahtani, R.T.; Ahmad, S.; Akgül, A. Dynamical analysis of bio-ethanol production model under generalized nonlocal operator in Caputo sense. Mathematics 2021, 9, 2370. [Google Scholar] [CrossRef]
  27. Jain, S. Numerical analysis for the fractional diffusion and fractional Buckmaster equation by the two-step Laplace Adam-Bashforth method. Eur. Phys. J. Plus 2018, 133, 19. [Google Scholar] [CrossRef]
  28. Atangana, A.; Jain, S. A new numerical approximation of the fractal ordinary differential equation. Eur. Phys. J. Plus 2018, 133, 37. [Google Scholar] [CrossRef]
  29. Saifullah, S.; Ali, A.; Irfan, M.; Shah, K. Time-fractional Klein–Gordon equation with solitary/shock waves solutions. Math. Probl. Eng. 2021, 2021, 6858592. [Google Scholar] [CrossRef]
  30. Araz, S.İ. Analysis of a COVID-19 model: Optimal control, stability and simulations. Alex. Eng. J. 2021, 60, 647–658. [Google Scholar] [CrossRef]
  31. Awais, M.; Alshammari, F.S.; Ullah, S.; Khan, M.A.; Islam, S. Modeling and simulation of the novel coronavirus in Caputo derivative. Results Phys. 2020, 19, 103588. [Google Scholar] [CrossRef] [PubMed]
  32. Ullah, S.; Khan, M.A. Modeling the impact of non-pharmaceutical interventions on the dynamics of novel coronavirus with optimal control analysis with a case study. Chaos Solitons Fractals 2020, 139, 110075. [Google Scholar] [CrossRef] [PubMed]
  33. Li, B.; Liang, H.; He, Q. Multiple and generic bifurcation analysis of a discrete Hindmarsh-Rose model. Chaos Solitons Fractals 2021, 146, 110856. [Google Scholar] [CrossRef]
  34. Li, B.; Liang, H.; Shi, L.; He, Q. Complex dynamics of Kopel model with nonsymmetric response between oligopolists. Chaos Solitons Fractals 2022, 156, 111860. [Google Scholar] [CrossRef]
  35. Eskandari, Z.; Avazzadeh, Z.; Ghaziani, R.K.; Li, B. Dynamics and bifurcations of a discrete-time Lotka–Volterra model using nonstandard finite difference discretization method. Math. Methods Appl. Sci. 2022. [Google Scholar] [CrossRef]
  36. Baleanu, D.; Abadi, M.H.; Jajarmi, A.; Vahid, K.Z.; Nieto, J.J. A new comparative study on the general fractional model of COVID-19 with isolation and quarantine effects. Alex. Eng. J. 2022, 61, 4779–4791. [Google Scholar] [CrossRef]
  37. Zhang, L.; Addai, E.; Ackora-Prah, J.; Arthur, Y.D.; Asamoah, J.K.K. Fractional-Order Ebola-Malaria Coinfection Model with a Focus on Detection and Treatment Rate. Comput. Math. Methods Med. 2022, 2022, 6502598. [Google Scholar] [CrossRef]
  38. Gómez-Aguilar, J.F.; Rosales-García, J.J.; Bernal-Alvarado, J.J.; Córdova-Fraga, T.; Guzmán-Cabrera, R. Fractional mechanical oscillators. Rev. Mex. Física 2012, 58, 348–352. [Google Scholar]
  39. Afshari, H.; Hosseinpour, H.; Marasi, H.R. Application of some new contractions for existence and uniqueness of differential equations involving Caputo–Fabrizio derivative. Adv. Differ. Equ. 2021, 2021, 321. [Google Scholar] [CrossRef]
  40. Ulam, S.M. A Collection of Mathematical Problems; Interscience Publishers: New York, NY, USA, 1960; p. 29. [Google Scholar]
  41. Ulam, S.M. Problem in Modren Mathematics; Courier Corporation: Mineola, NY, USA, 2004. [Google Scholar]
Figure 1. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Figure 1. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Fractalfract 07 00046 g001
Figure 2. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Figure 2. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Fractalfract 07 00046 g002
Figure 3. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Figure 3. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Fractalfract 07 00046 g003
Figure 4. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Figure 4. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Fractalfract 07 00046 g004
Figure 5. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Figure 5. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Fractalfract 07 00046 g005
Figure 6. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Figure 6. The dynamics of the state variables’ SAIR in model (2) with different fractional orders α vs. time t with initial conditions.
Fractalfract 07 00046 g006
Table 1. The parameters and their values for model (2).
Table 1. The parameters and their values for model (2).
ParametersSet 1Set 2
Π 26002600
0.60.3
υ 0.0650.065
q 1 0.30.3
q 2 0.10.1
κ 0.30.3
ν 0.10.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alkahtani, B.S. Mathematical Modeling of COVID-19 Transmission Using a Fractional Order Derivative. Fractal Fract. 2023, 7, 46. https://doi.org/10.3390/fractalfract7010046

AMA Style

Alkahtani BS. Mathematical Modeling of COVID-19 Transmission Using a Fractional Order Derivative. Fractal and Fractional. 2023; 7(1):46. https://doi.org/10.3390/fractalfract7010046

Chicago/Turabian Style

Alkahtani, Badr S. 2023. "Mathematical Modeling of COVID-19 Transmission Using a Fractional Order Derivative" Fractal and Fractional 7, no. 1: 46. https://doi.org/10.3390/fractalfract7010046

APA Style

Alkahtani, B. S. (2023). Mathematical Modeling of COVID-19 Transmission Using a Fractional Order Derivative. Fractal and Fractional, 7(1), 46. https://doi.org/10.3390/fractalfract7010046

Article Metrics

Back to TopTop