Caputo Fractional Derivative for Analysis of COVID-19 and HIV/AIDS Transmission

Cheneke, Kumama Regassa

doi:https://doi.org/10.1155/2023/6371148

Abstract and Applied Analysis

On this page

Abstract Introduction Model Formulation Results and Discussion Conclusion Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2023 | Article ID 6371148 | https://doi.org/10.1155/2023/6371148

Caputo Fractional Derivative for Analysis of COVID-19 and HIV/AIDS Transmission

Kumama Regassa Cheneke^1,2

Academic Editor: Abd Allah Hyder

Received24 Jun 2022

Revised11 Aug 2023

Accepted30 Aug 2023

Published29 Sept 2023

Abstract

In this study, Caputo fractional derivative model of HIV and COVID-19 infections is analyzed. Moreover, the well-posedness of a model is verified to depict that the developed model is mathematically meaningful and biologically acceptable. Particularly, Mittag Leffler function is incorporated to show that total population size is bounded whereas fixed point theory is applied to show the existence and uniqueness of solution of the constructed Caputo fractional derivative model of HIV and COVID-19 infections. The study depicts that as the order of fractional derivative increase the size of the infected variable decrease as time increase. Additionally, memory effects correspond to order of derivative in the reduction of a number of populations infected both with HIV and COVID-19 infections. Numerical simulations are performed using MATLAB platform.

1. Introduction

COVID-19 wreaked havoc on both human life and the immense economic development. A severe economic crisis and the struggle for survival gripped the entire planet. Even the most economically developed nation lost faith and was unable to supply the urgently required resources asked by doctors at medical facilities [1]. The dynamic transmission of COVID-19 crises and attack on the upper respiratory organs cause the breathing system to quickly narrow, which has a significant negative impact [2]. Additionally, the prevalence of the virus has decreased overall as a result of the use of masks and the COVID-19 vaccine that reinforced through activation of memory effects by public health interventions [3]. Recently, in the worldwide it is registered that 692,576,573 cases, 6,903,976 deaths, and 664,687,106 recovered of COVID-19 [4]. On the other hand, the transmission of human immunodeficiency virus transmission through act of sexual practices, blood transfusion, and possible exposure to the virus is a burden in the worldwide [5]. Hence, different mathematical models have been developed to address the impact of diseases in the whole world [6–8]. For instance in engineering problems, the application of a regularized ψ–Hilfer fractional derivative is described by Jajarmi et al. [9] whereas Caputo fractional derivative is applied in the analysis of accelerated mass-spring system by Defterli et al. [10]. Recently, fractional derivative get a great attention due to indication of order of fractional derivative impact on memory effects to control the transmission and prevalence of infection in the population [11]. Caputo fractional derivative is applied in controlling the transmission dynamics of Nipah virus [12]. Also, some studies supported with mathematical modeling show that the occurrence of COVID-19 has significant impact on the HIV carriers [13, 14]. Moreover, the availability of the infections within the community inspires further investigation to consider the impact of memory effects through public health education intervention [15–18]. Also, the Caputo fractional derivative is incorporated to investigate banking data competition [19]. Recently, fractional derivative models are fitted to real data to describe the behavior of the real-world phenomena [20–26]. Therefore, in this study the memory effect is incorporated and a model with Caputo fractional derivative is developed and analyzed. The remaining portion of the paper is divided into the following sections: Section 2, which introduces preliminary ideas and the formulation of a mathematical model; Section 3, which analyzes the fractional model mathematically; Section 4, which uses numerical simulation; Section 5, which discusses the results; and Section 6, which contains the paper’s conclusion.

2. Mathematical Preliminary and Model Formulation

2.1. Mathematical Preliminary

In this section, revise basic definitions and theorems in fractional calculus that supports the main results of the current study.

Definition 2.1. [27] The order fractional integral of Riemann–Liouville of a function with is given by

Definition 2.2. [28] Let be an element of the space of absolutely continuous function on and . The left and right Caputo fractional derivative of a function with order is given by

Definition 2.3. [29] The Laplace transform of the left Caputo fractional derivative is defined as follows:

Definition 2.4. [29] The Laplace transform of the function is given bywhere is the two-parameter Mittage–Leffler function with parameters .

Definition 2.5. [29] Let be Mittage–Leffler function. Then

Lemma 1. (Generalized Mean Value Theorem) Suppose that and for , then we havewith , for all .

Moreover, if , then there is a neighborhood N of such that . Also, if , then there is a neighborhood N of such that .

2.2. Model Formulation

To construct our current model, we have modified the coinfection model of HIV and cholera virus developed in [30]. We formulate model of HIV and COVID-19 coinfection by dividing the total population under consideration as classes consists of (i) susceptible individuals (S). They are infection free individuals with possibility of acquiring HIV from HIV infected individuals through sexual contact and acquiring COVID-19 from COVID-19 infected individuals only through effective direct contacts; (ii) COVID-19 infected individuals (C). They are individuals infected with COVID-19 and capable of transmitting infection to susceptible population through direct contact; (iii) COVID-19 recovered individuals (R). They are individuals recovered from COVID-19 infection with temporary immunity; (iv) HIV/AIDS infected individuals (H). They are individual infected with only HIV/AIDS and capable of transmitting HIV to susceptible population through sexual contacts; (v) coinfected individuals (). They are individuals infected with HIV and COVID-19. They can transmit disease with effective contacts with susceptible population; (vi) COVID-19 recovered and HIV/AIDS individuals (). They are HIV/AIDS individuals recovered from COVID-19 with temporary immunity. The following assumptions are also stated:(i)For this study a total population size is not constant;(ii)Susceptible persons recruited at the rate ;(iii)A naturally death rate is ;(iv)Only COVID-19 infected individuals die as a result of infection at the rate ;(v)Only HIV/AIDS infected individuals die at the rate ;(vi)Death rate of coinfected individuals die due to infection at the rate ;(vii)COVID-19 transmission rate due to direct contact is ;(viii)Transmission rate of HIV/AIDS due to unsafe exposure to virus is ;(ix)Recovery rate of only COVID-19 infected with temporary immunity is ;(x)Recovery rate of coinfected with temporary immunity is .(xi) is immunity loss rate of COVID-19 recovered individuals ;(xii) is COVID-19 infection vaccination rate for only HIV infected individuals;(xiii) is COVID-19 infection vaccination rate for only susceptible individuals;(xiv) is immunity loss rate of COVID-19 infection recovered individuals .

The total size of population under study is denoted by and defined as where is the size of susceptible population at time , size of COVID-19 population at time , size of only COVID-19 recovered population at time , size of HIV/AIDS population at time , size of coinfected population at time , and size of COVID-19 recovered HIV population at time.

The aforementioned assumptions and supporting flow diagram described in Figure 1 lead to a model with Caputo fractional derivative as given by subsequent equations.with initial value conditions for the problem .

3. Fractional Model Analysis

3.1. Invariant Region

Theorem 1. The solutions of the model (1) are invariant in the region proper subset of six-dimensional space over nonnegative real numbers such that

Proof. To show the boundedness of solution we employ the method employed by Qian et al. [31]. So, adding corresponding terms on left and right of qualities in model (1) we obtainImplies,Moreover, for computation purpose, we set the preceding expression to the formAlso, from the preceding expression, we obtainMoreover, applying the Laplace transform on both sides of the preceding equation yieldsAccording to Podlubny [29], applying Laplace transform definition in Caputo fractional derivative sense such that , the preceding equation reduces to the formAlso, for , he preceding equation reduces toAlso, by summation property the preceding equation reduces toMoreover, arranging the terms, the preceding equation can be written as follows:Further, solving for , the foregoing equation reduces toAlso, applying inverse Laplace transform on both sides of the preceding equation, we obtainAgain applying property of inverse Laplace transform, the preceding equation givesAlso, using the relationship between inverse Laplace transform and Mittag Leffler function given by Qian et al. [31], , the preceding equation reduces toImplies,Also, the definition of Mittag Leffler function, , the preceding equation reduces toImplies,Hence, taking into consideration the Equation (3), we obtain:Hence, as time gets larger and larger the total population size is bounded between 0 and .

3.2. Positivity Property

Theorem 2. The solution of the developed fractional model (1) is positive for all time in the invariant region .

Proof. To show positivity of solutions, we employ the techniques applied by Baleanu et al. [12]. In similar fashion, the trajectory of solution of solution of model (1) along only one state-axis, where other state variables vanishes, givesSimilar to the procedures done for the boundedness computed in Section 3.1, the solution of the foregoing equation is solved, respectively, as follows:Therefore, the preceding computed results show that the solution is positively invariant along the axis of state variables. Moreover, since the solution of the fractional model (1) is positive in the plane plane, let such that and . On this plane,Moreover, by Caputo fractional derivative mean value theorem sated and applied by Musa et al. [32], we haveTherefore, we get which contradicts our earlier assumption for . Hence, the state variable is nonnegative for all time . Similarly, other state variables are nonnegative for all time . Therefore, any solution of fractional derivative is nonnegative for all time .

3.3. Existence and Uniqueness of Solution

In this section, as presented by Ahmad et al. [28], we show uniqueness and existence of the solution. However, before we prove the existence and uniqueness of solution define the kernel functions from fractional derivative model (1) as follows:

Moreover, let .

Then model (1) can be written as follows:

In the preceding expression, the condition is to be taken component wise. Problem (32) which is equivalent to model (1), can be described by the integral of

Next we shall analyze model (1) through the integral representation given above. For that case, let denotes Banach space of all continuous functions that maps from (0, ) to endowed with the normwhere . Note that, all belong to . Moreover, we define the operator by

Hence, the operator is well-defined due to obvious continuity of .

Theorem 4. Let , the function defined above satisfiesfor some .

Proof. The first component of gives,Let , then the preceding equation can be reduced towhere,Therefore,In similar fashion the remaining can be shown. Consequently, we can conclude thatwhere .

Theorem 5. Let the result of preceding theorem holds . If , then there exists a unique solution of model (1) on () which is uniformly Lyapunov stable.

Proof. The function is obviously continuous on its stated domain. Thus, the existence of the solution follows from the works done by Ahmad et al. [28].
For uniqueness, we apply Banach contraction mapping on operator defined above. Next, we show is both a self-map and a contraction. By definition, . Let us now define and a closed convex set . Thus, for self-map property it is enough to show that . So let , thenIt follows that the operator and is indeed a self-map. Next, we prove that is a contraction. Let and satisfy the abbreviated dynamical system. The using the result of theorem above, we haveHence, if then is contraction mapping and by Banach contraction mapping principle, has a unique fixed point on which is a solution of model (1). Moreover, the uniformly Lyapunov stability of solutions follows as stated by Ahmad et al. [28].

3.4. Disease-Free Equilibrium

The disease-free equilibrium, of fractional model (1) is a steady state point where the disease extinct. Setting the state variables for disease equal to zero, at steady state model (1) gives a disease-free equilibrium as follows:

That is, we set model (1) equal to zero and solve as given below

Moreover, at vanishing point of , and , we obtain

Therefore, the disease-free equilibrium of fractional model is computed by applying simple mathematical calculation as follows:

3.5. Endemic Equilibrium

The endemic equilibrium, , of fractional derivative model (1) is a point where the diseases persist in the population. Hence, we can write the endemic equilibrium as follows:where , and are obtained by solving the subsequent equation

Moreover, due to trick of the problem, we obtain the endemic equilibrium from numerical simulation.

3.6. Basic Reproduction Number

The basic reproduction number is the indicator of disease status at any time of investigation and defined as the average number of infected persons flow into a fully susceptible population due to generation of infected individuals by typical infectious individual. In mathematical biology of epidemiology, the next-generation matrix method is widely applied to systematically obtain basic reproduction number from a model with more compartments. Let consider model (1) and construct the subsequent Jacobian matrices at disease-free equilibrium as follows:

Moreover, the next-generation matrix constructed from the preceding matrices and is given bywhere .

The eigenvalues of the preceding matrices are computed as follows:

We know that the basic reproduction number is the spectral radius of next-generation matrix. Thus, the basic reproduction is computed as given below

4. Numerical Simulations

4.1. Numerical Schemes for Caputo Fractional Derivative

The numerical solution of Caputo fractional derivative is performed using FDE12 [33]. Moreover, the initial population sizes are . Moreover, parameter value is given in Table 1.

5. Results and Discussion

In Figure 2 the size of susceptible population decrease as time increase which shows the higher activation of human memory toward vaccination the less the population are susceptible. Moreover, the higher the order of fractional derivative the higher the mobilization of population toward vaccination. Figure 3 describes the dynamics of population size subject to the transmission dynamics of novel COVID-19 and memory effects through public health education. It is noticed that the intervention to activate memory of population significantly contributes to control the pandemic COVID-19. Moreover, the simulation shows that less memory effect contributes less in controlling of the infection whereas high involvement in the activation of human memory results in better way to control the infection.

In Figure 4, dynamics of only HIV-infected population size is simulated. It is observed that the memory effect works well to control the transmission dynamics of HIV infection. Moreover, the simulation depicts there is comprise of population in the long run. Thus, effectively apply the memory effects before HIV infects many people. In Figure 5, the simulation of coinfected population size is depicted. It is shown that the number of individuals infected with both infections is high at the beginning and the memory effect is less effective. On the other hand, the memory effects work well for reducing coinfected individuals as time increase. In Figure 6, the dynamics of only HIV infected population with COVID-19 recovered status is simulated. It is observed that memory effect activation works continually for controlling COVID-19 for HIV infected population. Hence, a continuum memory effect helps HIV-infected population for recovery from pandemic COVID-19.

In Figure 7, the dynamics of population size is simulated. It is shown that the impact of COVID-19 can be controlled with effective memory effect whereas the HIV infection persists with population. In Figure 8, we observe that the more individuals are susceptible the more they get infected with COVID-19 infection whereas less the number of susceptible individuals the less the number of individuals infected with COVID-19 infection. In Figure 9, we observe that the more individuals are susceptible to HIV the more they get attacked with HIV infection whereas the less the number of susceptible individuals toward HIV the less the number of individuals exposed to the HIV infection. Figure 10 depicts the transferring of individuals to R(t) with immunity reduces the number of susceptible individuals.

Figure 11 shows that the number of coinfected individuals decreases or increases along susceptible population size. Figure 12 describes that the number of HIV infection recovered individuals increases as susceptible decrease, but lately decrease as the number of susceptible individuals decreases. Figure 13 describes the dynamics of , and . The behavior of and alike as they decreases from the beginning whereas increase initially and decrease lately. Figure 14 describes the dynamics of , and . The behaviors of and are similarly get decrease from the beginning whereas increase initially and decrease lately. Figure 15 describes the dynamics of , and . The behaviors of and are similarly get increases from the beginning and decreases lately whereas decreases from the beginning.

6. Conclusion

In this study, fractional derivative is incorporated and the model of HIV and COVID-19 dynamics is developed and analyzed. Based on the simulations of HIV and COVID-19 dynamics, infections can be controlled for basic reproduction number is less than unity. The study depicts that HIV/AIDS and COVID-19 coinfection may increase, through direct interaction of coinfected individuals. Also, fractional derivative analysis shows increasing memory of population toward COVID-19 and HIV infections will increase the ability of population to control the infections. In the future study, one can extend our model to a new generalized fractional derivative and use real data to describe the biological phenomena. Numerical simulations are performed using MATLAB platform.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declares that there is no conflict of interest.

References

K. Hattaf, A. A. Mohsen, J. Harraq, and N. Achtaich, “Modeling the dynamics of COVID-19 with carrier effect and environmental contamination,” International Journal of Modeling, Simulation, and Scientific, vol. 12, no. 3, Article ID 2150048, 2021.
View at: Publisher Site | Google Scholar
K. Hattaf, M. I. El Karimi, A. A. Mohsen, Z. Hajhouji, M. El Younoussi, and N. Yousfi, “Mathematical modeling and analysis of the dynamics of RNA viruses in presence of immunity and treatment: a case study of SARS-CoV-2,” Vaccines, vol. 11, no. 2, Article ID 201, 2023.
View at: Google Scholar
K. Hattaf and N. Yousfi, “Dynamics of SARS-CoV-2 infection model with two modes of transmission and immune response,” Mathematical Biosciences, vol. 17, no. 5, pp. 5326–5340, 2020.
View at: Google Scholar
“02: 21 GMT,” August 2023, https://www.worldometers.info/ corona virusvirus/.
View at: Google Scholar
A. A. Mohsen and K. Hattaf, “Dynamics of a generalized fractional epidemic model of COVID-19 with carrier effect,” Advances in Systems Science and Applications, vol. 22, no. 3, pp. 36–48, 2022.
View at: Google Scholar
K. Hattaf and H. Dutta, Mathematical Modeling and Analysis of Infectious Diseases, vol. 302, Springer, Berlin, 2020.
K. Hattaf and H. Dutta, “Modeling the dynamics of viral infections in presence of latently infected cells,” Chaos, Solitons & Fractals, vol. 136, Article ID 109916, 2020.
View at: Publisher Site | Google Scholar
D. Baleanu, A. Jajarmi, H. Mohammadi, and S. Rezapour, “A new study on the mathematical modelling of human liver with Caputo–Fabrizio fractional derivative,” Chaos, Solitons & Fractals, vol. 134, Article ID 109705, 2020.
View at: Google Scholar
A. Jajarmi, D. Baleanu, S. S. Sajjadi, and J. J. Nieto, “Analysis and some applications of a regularized ψ–Hilfer fractional derivative,” Journal of Computational and Applied, vol. 415, Article ID 114476, 2022.
View at: Google Scholar
O. Defterli, D. Baleanu, A. Jajarmi, S. S. Sajjadi, N. Alshaikh, and J. Asad, “Fractional treatment: an accelerated mass-spring system,” Romanian Reports in Physics, vol. 74, no. 4, Article ID 122, 2022.
View at: Google Scholar
K. Hattaf, “A new generalized definition of fractional derivative with non-singular kernel,” Computation, vol. 8, no. 2, Article ID 49, 2020.
View at: Google Scholar
D. Baleanu, P. Shekari, L. Torkzadeh, H. Ranjbar, and A. Jajarmi, “Stability analysis and system properties of Nipah virus transmission: a fractional calculus case study,” Chaos, Solitons & Fractals, vol. 166, Article ID 112990, 2023.
View at: Google Scholar
B. L. Jewell, E. Medium, J. Stover et al., “Potential effects of disruption to HIV programmes in sub-Saharan Africa caused by COVID-19: results from multiple mathematical models.,” The Lancet HIV, vol. 7, no. 9, pp. e629–e640, 2020.
View at: Google Scholar
E. F. D. Goufo, Y. Khan, and Q. A. Chaudhry, “HIV and shifting epicenters for COVID-19, an alert for some countries,” Chaos, Solitons & Fractals, vol. 139, Article ID 110030, 2020.
View at: Google Scholar
S. Ahmad, S. Owyed, A. H. Abdel-Aty, E. E. Mahmoud, K. Shah, and H. Alrabaiah, “Mathematical analysis of COVID-19 via new mathematical model,” Chaos, Solitons & Fractals, vol. 143, Article ID 110585, 2021.
View at: Google Scholar
L. P. James, J. A. Salomon, C. O. Buckee, and N. A. Menzies, “The use and misuse of mathematical modeling for infectious disease policymaking: lessons for the COVID-19 pandemic,” Medical Decision Making, vol. 41, no. 4, pp. 379–385, 2021.
View at: Publisher Site | Google Scholar
Fatmawati, E. Yuliani, C. Alfiniyah, M. L. Juga, and C. W. Chukwu, “On the modeling of COVID-19 transmission dynamics with two strains: insight through caputo fractional derivative,” Fractal and Fractional, vol. 6, no. 7, Article ID 346, 2022.
View at: Publisher Site | Google Scholar
K. R. Cheneke, K. P. Rao, and G. K. Edessa, “Application of a new generalized fractional derivative and rank of control measures on cholera transmission dynamics,” International Journal of Mathematics and Mathematical Sciences, vol. 2021, Article ID 2104051, 9 pages, 2021.
View at: Publisher Site | Google Scholar
M. H. DarAssi, M. A. Khan, Fatmawati, and M. S. Alqarni, “Analysis of the competition system using parameterized fractional differential equations: application to real data,” Symmetry, vol. 15, no. 2, Article ID 542, 2023.
View at: Google Scholar
P. A. Naik, M. Yavuz, S. Qureshi, J. Zu, and S. Townley, “Modeling and analysis of COVID-19 epidemics with treatment in fractional derivatives using real data from Pakistan,” The European Physical Journal Plus, vol. 135, pp. 1–42, 2020.
View at: Google Scholar
S. Qureshi and A. Atangana, “Fractal–fractional differentiation for the modeling and mathematical analysis of nonlinear diarrhea transmission dynamics under the use of real data,” Chaos, Solitons & Fractals, vol. 136, Article ID 109812, 2020.
View at: Google Scholar
Sania Qureshi and Abdullahi Yusuf, “Fractional derivatives applied to MSEIR problems: comparative study with real world data,” The European Physical Journal, vol. 134, no. 4, Article ID 171, 2019.
View at: Publisher Site | Google Scholar
A. Atangana and M. A. Khan, “Modeling and analysis of competition model of bank data with fractal–fractional Caputo–Fabrizio operator,” Alexandria Engineering Journal, vol. 59, no. 4, pp. 1985–1998, 2020.
View at: Google Scholar
I. Area, H. Batarfi, J. Losada, J. J. Nieto, W. Shammakh, and Á. Torres, “On a fractional order Ebola epidemic model,” Advances in Difference, vol. 2015, no. 1, pp. 1–12, 2015.
View at: Google Scholar
D. Baleanu, F. A. Ghassabzade, J. J. Nieto, and A. Jajarmi, “On a new and generalized fractional model for a real cholera outbreak,” Alexandria Engineering Journal, vol. 61, no. 11, pp. 9175–9186, 2022.
View at: Google Scholar
C. J. Silva, C. Cruz, D. F. Torres et al., “Optimal control of the COVID-19 pandemic: controlled sanitary deconfinement in Portugal,” Scientific Reports, vol. 11, no. 1, Article ID 3451, 2021.
View at: Publisher Site | Google Scholar
K. Logeswari, C. Ravichandran, and K. S. Nisar, “Mathematical model for spreading of COVID-19 virus with the Mittag–Leffler kernel,” Numerical Methods for Partial Differential Equations, vol. 26, Article ID 33362342, 2020.
View at: Publisher Site | Google Scholar
S. Ahmad, A. Ullah, Q. M. Al-Mdallal, H. Khan, K. Shah, and A. Khan, “Fractional order mathematical modeling of COVID-19 transmission,” Chaos, Solitons & Fractals, vol. 139, Article ID 110256, 2020.
View at: Google Scholar
I. Podlubny, Fractional Differential Equations, Academic Press, New York, 2009.
K. R. Cheneke, K. P. Rao, and G. K. Edessa, “Modeling and analysis of HIV and cholera direct transmission with optimal control,” Discrete Dynamics in Nature and Society, vol. 2022, 16 pages, 2022.
View at: Publisher Site | Google Scholar
D. Qian, C. Li, R. P. Agarwal, and P. J. Wong, “Stability analysis of fractional differential system with Riemann–Liouville derivative,” Mathematical and Computer Modelling, vol. 52, no. 5-6, pp. 862–874, 2010.
View at: Google Scholar
S. S. Musa, S. Qureshi, S. Zhao, A. Yusuf, U. T. Mustapha, and D. He, “Mathematical modeling of COVID-19 epidemic with effect of awareness programs,” Infectious Disease Modeling, vol. 6, pp. 448–460, 2021.
View at: Google Scholar
S. Rosa and D. F. Torres, “Numerical fractional optimal control of respiratory syncytial virus infection in Octave/MATLAB,” Mathematics, vol. 11, no. 6, Article ID 1511, 2023.
View at: Google Scholar

Copyright

Copyright © 2023 Kumama Regassa Cheneke. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

900

Downloads

517

Citations