Abstract
This research work investigates some theoretical and semi-analytical results for the mathematical model of tuberculosis disease via derivative due to Caputo and Fabrizio. The concerned derivative involves exponential kernel and very recently it has been adapted for various applied problems. The required results are established by using some fixed point approach of Krasnoselskii and Banach. Further, by the use of iterative tools of Adomian decomposition and Laplace, the semi-analytical results are studied. Some graphical results are given with discussion.
Similar content being viewed by others
1 Introduction
One of the most important diseases faced by human beings ever is tuberculosis (TB). TB is a spreadable, airborne bacterial infection which is caused by mycobacterium tuberculosis. This bacterium normally affects the lungs (pulmonary tuberculosis). This bacterium may also affect several other systems like kidneys, the brain, the lymphatic system, the central nervous system, spinal cord, etc. The presence of TB disease has been found since ancient times in various civilizations like Egypt, China, Roman, etc. (see [1]). One-third of the world population at the present time is infected with TB, and the number of infectious individuals increases at a rate of one per second [2]. The aforementioned disease was among the top ten death causes around the globe in the year 2015, where about 10.4 million individuals were infected from it. In the same year, about 1.8 million infectious individuals lost their lives from various diseases including 0.4 million with HIV and TB. 60% of the tuberculosis cases around the globe were concentrated in the six countries (Pakistan, India, Nigeria, China, Indonesia, and South Africa) [3]. Dye [4] gave some information that the major cause of death worldwide, in particularly in Sub-Saharan Africa, is due to TB and HIV. Further, HIV epidemic is a serious threat in many countries of the world. It is a clear evidence that worldwide children are protected from the early infection of the disease by vaccination like Bacillus Calmatte–Guerine(BCG) [5]. Therefore, detection and treatment of latent TB by modern therapy have been used recently to prevent the breakdown of rate of spread of the disease as only the members of the infectious class can transmit the disease to others.
Worldwide various procedures and methods have been used to understand the cause and control of these diseases in society. One of the powerful tools is the mathematical modeling through which we understand the dynamics of various disease transmission and suggest procedure how to control them in society. The mentioned area originated during 1927 for the first time. A variety of models have been developed and studied so far (refer to [6–10]). In this regard a five-compartment model for TB has also been constructed in [11] as follows:
In the above model the whole populace is categorized into five classes: the immunized class \(\mathscr{M}\), the susceptible class \(\mathscr{S}\), the infected latently class \(\mathscr{L}\), the infectious class \(\mathscr{I}\), and the recovered class \(\mathscr{R}\). Parameters of the model under consideration are explained as follows: The constant of recruitment is represented by the symbol ρ, θ denotes the immunized portion at birth, α represents the rate of weaning off the vaccine, the natural death rate is denoted by the symbol μ, β represents the tuberculosis contraction rate, the successful treatment of infectious latent is denoted by the letter σ, the symbol τ is the rate of breakdown of latent TB into infectious TB, the successful cure of infectious TB individuals and the death resulting from the disease are respectively denoted by the symbols γ, δ.
Usually an integer order derivatives do not explore the dynamics of real world problems related to biology and physics well. In order to overcome this deficiency, fractional calculus has been given attention for the last few decades. Also we know that fractional calculus is increasingly used by mathematicians for mathematical modeling. Derivatives and integrals of noninteger order may be defined by a number of ways. Some well-known definitions are those given by Riemann and Liouville [12], Caputo, etc. (see [13]). The mentioned derivatives involve kernel of singular type. Frequently these two definitions have been increasingly used since fractional differential operator is in fact a definite integral operator for whom the definition of kernel is not unique or not regular. Further due to high degree of freedom in derivative of arbitrary order, researchers have given much attention to studying applied problems under these concepts. In this sense very recently some authors replaced the singular kernel by some nonlocal nonsingular and produced new definitions. Hence in 2015 Caputo and Fabrizio replaced the singular kernel by exponential kernel in the usual Caputo derivative and called it Caputo–Fabrizio fractional derivative (abbreviated as CFFD); for details, see [14–16]. Therefore various researchers investigated different problems of applied nature under this concept. In various cases the mentioned derivative has produced significant results as compared to other forms of derivatives (see [17–19]). The Caputo–Fabrizio derivative omits singular kernel by exponential kernel and hence makes the concerned differential operator nonlocal. Conventional fractional derivatives contain singular kernel which sometimes causes difficulty in explanation of some characteristics of various materials. To overcome this, Caputo and Fabrizio introduced a new definition of fractional integral and derivative which involves exponential kernel instead of singular one. Various studies can be found in the literature that have focused on the Caputo–Fabrizio fractional order derivatives, see, for instance, [20–24]. Further, in various papers related to thermal sciences, the mentioned derivative has been proved powerful as compared to other type, we refer to [25–30]. Keeping these points in mind and the nonsingular nature of the proposed derivative, we investigate the considered model for existence and analytical results.
Now the question is how to treat the problems involving derivative of fractional orders. For this need researchers have successfully updated the usual tools and methods to handle differential equations of fractional order (FODEs). Usual perturbation techniques and decomposition methods were greatly utilized to deal with ordinary FODEs. Also, for the mentioned problems, Adomian decomposition coupled with some integral transforms has been used very well (see [31–34]). On the other hand, since the FODEs involving new type derivatives are very rarely used, very frequently authors [35] established some algorithms to handle such FODEs containing CFFD.
Hence we investigate the model given in (1) under CFFD as follows:
We study model (2) subject to the biologically feasible initial conditions
Initially we establish some conditions about the existence of solutions for model (2) by using some fixed point results like Banach and Krasnoselskii. After that, by using the considered tool of “Laplace Adomian decomposition method (LADM)” for \(\eta \in (0,1]\), we compute semi-analytical results. Finally, the approximate results are presented via graphs.
2 Preliminaries
In this section of the manuscript, we present some fundamental definitions as follows.
Definition 1
([15])
Let \(\varphi \in \mathcal{H}^{1}(a,b)\), \(b>a\), \(\eta \in (0,1)\), then the CFFD is given as
where \(\mathcal{K(\eta )}\) in (3) is the normalization function with \(\mathcal{K}(1)=\mathcal{K}(0)=1\). If the function failed to exist in \(\mathcal{H}^{1}(a,b)\), then the above derivative can be reformulated as
Definition 2
([36])
Let \(\eta \in (0,1]\), then the integral of fractional order η of the function φ is
Definition 3
The Laplace transform of CFFD \({}_{0}^{CF}D_{t}^{\eta }\) of \(M(t) \) is given as
2.1 Equilibrium points and the basic reproduction number
Before proceeding further, we consider it advantageous to find the equilibrium points and the basic reproduction number of the model under consideration. There are two types of possible equilibrium points of the model. The first one is the point where there is no disease in the community, i.e., the disease-free equilibrium point. This is obtained by setting the right-hand side of each equation in the model to zero along with \(\mathscr{L}=\mathscr{I}=\mathscr{R}=0\). Solving the system then gives \(\mathscr{M}^{0}=\frac{\theta \rho }{\alpha +\mu }\) and \(\mathscr{S}^{0}=\frac{\alpha +\mu -\mu \theta }{\alpha +\mu }\). Thus the disease-free equilibrium point of the model under investigation is given by \(\mathscr{E}^{0}= (\frac{\theta \rho }{\alpha +\mu }, \frac{\alpha +\mu -\mu \theta }{\alpha +\mu },0,0,0 )\).
To find the basic reproduction number, we consider only the infectious classes of the model. Let \(\mathbf{{\mathit{{V}}}= (\mathscr{L},\mathscr{I} )^{T}}\), by the help of the given model one can write . The Jacobian matrices of \(\mathscr{F}\) and \(\mathscr{V}\) are given by and . The inverse matrix of \(\mathcal{V}\) is given by . Hence the next generation matrix \(\mathcal{F}\mathcal{V}^{-1}\) is calculated as
The spectral radius of the next generation matrix (4) gives the threshold quantity \(R_{0}\) [38]. Thus
Three-dimensional plots of the basic reproduction number \(R_{0}\) versus different parameters in the model under consideration are depicted in Fig. 1. This quantity plays the key role in stability analysis and in finding conditions for the said purpose.
3 Existence and uniqueness results for tuberculosis disease model of fractional order
In the following we derive existence results related to our model (2) exploiting the so-called fixed point theorem due to Banach. To proceed further, we first of all define the functions given below:
where
So our problem becomes
where \(M(0)=\mathbf{N}_{1} \), \(S(0)=S_{0} \), \(L(0)=L_{0} \), \(I(0)=I_{0} \), \(R(0)=R_{0}\). Application of \({}_{0}^{CF}I^{\eta }\) on both sides of (2) gives the following system of integral equations:
where \(G=\frac{(1-\eta )}{ \mathcal{K}(\eta )}\), \(\overline{G}= \frac{\eta }{\mathcal{K}(\eta )} \). Further, we will use \(f_{0i}=f_{i}(0, M(0), S(0), L(0), I(0), R(0))\), \(i=1,2,3,4,5\). Using the initial conditions, we have
Here, we denote Banach space by \(X=C([0, T]\times \mathscr{R}^{5}, \mathscr{R})\) under the norm
where \(T>0\) such that \(0\leq t\leq T<\infty \).
Theorem 1
(Krasnoselskii fixed point theorem)
Let X be a Banach space and D be a closed and convex subset of X, then there exist two operators A, B for which the following hold:
-
1.
The sum \(Ax+By\) belongs to D;
-
2.
The operator A is a contraction, while the operator B is continuous and compact;
-
3.
∃ at least one solution z in a way that \(Az+Bz=z\) holds.
Let us assume
Therefore, system (8) reduces to
For further analysis, we suppose that the following assumptions hold:
- (\(\mathcal{H}_{1}\)):
-
There exists \(\mathcal{K}_{F}>0\) for which
$$ \bigl\vert F(t,W)-F(t,\overline{W}) \bigr\vert \leq \mathcal{K}_{F} \vert W-\overline{W} \vert . $$ - (\(\mathcal{H}_{2}\)):
-
There exist \(C_{F}>0\) and \(M_{F}>0\) such that
$$ \bigl\vert F(t,W) \bigr\vert \leq C_{F} \vert W \vert +M_{F}. $$
Theorem 2
In the view of Theorem 1, problem (9) has at least one solution provided \(G \mathcal{K}_{F}\) is less than unity.
Proof
To prove the theorem, we define a compact and closed set D such that \(D={W \in X: \|W\|\leq r}\). Next, we define operators A and B as follows:
To verify that the operator A in (10) is a contraction, we assume \(W,\overline{W} \in X\), so that
from which we have
which clearly indicates that the operator A is a contraction.
Now we show that the operator B is compact and continuous. Consider
Taking max of (11), we have
This implies that B is bounded in (12). Let us assume that in the domain of t we have \(t_{1}< t_{2}\). One may write
If \(t_{2}\) approaches \(t_{1}\), then the right-hand side of (13) goes to zero. Consequently, \(t_{2} \rightarrow t_{1}\), which leads to
It follows that B is equicontinuous and, consequently, B is compact continuous. This implies that B is a completely continuous operator. Hence, all the conditions of Theorem 1 are satisfied. One may immediately conclude that model (2) has at least one solution. □
Theorem 3
There is a unique solution of the model under consideration (2) if the functions \(f_{1}\), \(i=1,2,3,4\), are continuous and \(\overline{G} K_{F}(1+T)<1\).
Proof
To prove the theorem, we define an operator \(P:X\rightarrow X\) such that
Let \(W, \overline{W} \in X\), we may write
It follows that
Consequently, model (2) under investigation has a unique solution. □
4 Construction of general algorithm for the required solution of the considered model
To derive the series type solution for the considered problem, we take \(\mathcal{K}(\eta )=1\) and apply the Laplace transform on both sides of (2). We construct the following algorithm:
After rearranging terms in (15), we have
Using the initial conditions of system (2), one has
Let the solution we compute be in the form of an infinite series as follows:
and decompose the nonlinear term SI in terms of Adomian polynomial as follows:
where \(A_{n}=\frac{1}{\Gamma {(n+1)}}\frac{d^{n}}{d\lambda ^{n}} [ ( \sum_{k=0}^{n} \lambda ^{k} S_{k} ) (\sum_{k=0}^{n} \lambda ^{k} I_{k} ) ]|_{\lambda =0}\)
In view of these values, model (8) becomes
Now, comparing terms on both sides of (18), we get the following series of problems.
Case 1. When \(n=0 \), we have
Evaluating the inverse Laplace transform, we get
Case 2. When \(n=1 \), we have
Evaluating the inverse Laplace transform, we get
Case 3. When \(n=2 \), we have
Evaluating the inverse Laplace transform, we get
Case 4. When \(n=3 \), we have
and so on. The other terms may similarly be computed.
Evaluating the inverse Laplace transform, we get
and so on. In this way, next terms of the series solution may be computed. Therefore, we get the required solution as follows:
Theorem 4
Let \(\mathscr{X}\) be the Banach space and \(\mathbf{T} : \mathscr{X}\rightarrow \mathscr{X}\) be a contractive nonlinear operator such that, for all \(W, {\bar{W}} \in \mathscr{X}\), \(\|\mathbf{T}(W)-\mathbf{T}({\bar{W}}) \|_{\mathscr{X}}\leq \kappa \|W-{\bar{w}}\|_{\mathscr{X}}\), \(0<\kappa <1\). Using the Banach contraction principle, T has a unique point W such that \(\mathbf{T}{W}={W}\), where \(W=(x,y,z)\). By applying LADM, the series given in (26) can be written as
and let \(W_{0}=W_{0}\in B_{\varepsilon }(W)\), where \(B_{\varepsilon }(W)={\bar{\mathbf{w}}\in \mathscr{X}:\| \bar{\mathbf{w}}-W\|_{\mathscr{X}}<\varepsilon }\), then one has
-
(i)
\(W_{n} \in B_{r}(W)\);
-
(ii)
\(\lim_{n\rightarrow \infty }W_{n}=W\).
Proof
The proof of the above theorem can be derived in a way similar to that in [39]. □
5 Numerical results and discussion
In this part of the paper, we present numerical results along with illustration regarding the approximate solution of the model under discussion. We take the approximate values for the parameters as given in Table 1. In light of these values, we get the series solution as follows:
Now we plot the solution up to five terms as given in (28) in Figures 1–5, corresponding to different fractional orders.
It can be observed from Fig. 1 that the immunized population decreases with different fractional orders at different ratio. In the same way the susceptible population is increasing as shown in Fig. 2. The infected and the latently infected population is also increasing as shown in Figs. 3 and 4, respectively. Because the susceptible population is converted to infected or latently infected. Proper cure is applied, then the recovered population will increase as shown in Fig. 5. The process of increase or decrease is initially fastest at lower fractional order and some time after the process is reversed, and the greater is the fractional order the faster is the increasing or decreasing process of population of the respective compartments. It means that fractional order derivatives can express the behavior more globally. The recovered population gradually increases and converges to equilibrium state as time progresses, as shown in Fig. 6.
6 Conclusion
We have investigated a biological model of TB under Caputo–Fabrizio fractional derivative. We have also established some sufficient results regarding the existence as well as the uniqueness of solution for the considered problem with the help of fixed point theorems. Further we have used a hybrid type method to compute series type solutions for the proposed model. To the best of our knowledge, the said techniques were very rarely used to handle the analytical solutions of FODEs involving nonsingular derivative of Caputo–Fabrizio type in the past. Further the numerical results have been displayed via graphs indicating that the established technique can be used to handle solution of those FODEs involving CFFD efficiently. Further, the mentioned method can be utilized to investigate more nonlinear problems of FODEs involving CFFD.
Availability of data and materials
Data sharing is not applicable to this paper.
References
Morse, D., Brothwell, D.R., Ucko, P.J.: Tuberculosis in ancient Egypt. Am. Rev. Respir. Dis. 90(4), 524–541 (1964)
Aparicio, J.P., Capurro, A.F., Castillo-Chavez, C.: Transmission and dynamics of tuberculosis on generalized households. J. Theor. Biol. 206, 327–341 (2000)
Floyd, K., Glaziou, P., Zumla, A., Raviglione, M.: The global tuberculosis epidemic and progress in care, prevention, and research: an overview in year 3 of the end TB era. Lancet Respir. Med. 6(4), 299–314 (2018)
Dye, C.: Global epidemiology of tuberculosis. Lancet 367(9514), 938–940 (2006)
Colditz, G.A., Brewer, T.F., Berkey, C.S., Wilson, M.E., Burdick, E., Fineberg, H.V., Mosteller, F.: Efficacy of BCG vaccine in the prevention of tuberculosis: meta-analysis of the published literature. JAMA 271(9), 698–702 (1994)
Arqub, O.A., El-Ajou, A.: Solution of the fractional epidemic model by homotopy analysis method. J. King Saud Univ., Sci. 25(1), 73–81 (2013)
Rafei, M., Ganji, D.D., Daniali, H.: Solution of the epidemic model by homotopy perturbation method. Appl. Math. Comput. 187(2), 1056–1062 (2007)
Zhao, S., Xu, Z., Lu, Y.: A mathematical model of hepatitis B virus transmission and its application for vaccination strategy in China. Int. J. Epidemiol. 29(4), 744–752 (2000)
Haq, F., Shah, K., Khan, A., Shahzad, M., Rahman, G.: Numerical solution of fractional order epidemic model of a vector born disease by Laplace Adomian decomposition method. Punjab Univ. J. Math. 49(2), 13–22 (2017)
Ullah, I., Ahmad, S., Al-Mdallal, Q., Khan, Z.A., Khan, H., Khan, A.: Stability analysis of a dynamical model of tuberculosis with incomplete treatment. Adv. Differ. Equ. 2020(1), 499 (2020)
Enagi, A.I., Ibrahim, M.O., Akinwande, N.I., Bawa, M., Wachin, A.A.: A mathematical model of tuberculosis control incorporating vaccination, latency and infectious treatments (case study of Nigeria). Int. J. Math. Comput. Sci. 12(2), 97 (2017)
Podlubny, I.: Fractional Differential Equations: Mathematics in Science and Engineering. Academic Press, New York (1999)
Kilbas, A.A., Srivastava, H.M., Trujillo, J.J.: Theory and Applications of Fractional Differential Equations. North-Holland, Amsterdam (2006)
Caputo, M., Fabrizio, M.: A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 1(2), 73–85 (2015)
Caputo, M., Fabrizio, M.: Application of new time and spatial fractional derivatives with exponential kernels. Prog. Fract. Differ. Appl. 2, 1–11 (2016)
Baleanu, D., Jajarmi, A., Mohammadi, H., Rezapour, S.: A new study on the mathematical modelling of human liver with Caputo–Fabrizio fractional derivative. Chaos Solitons Fractals 134, 109705 (2020)
El-Saka, H.A.A.: The fractional-order SIS epidemic model with variable population size. J. Egypt. Math. Soc. 22(1), 50–54 (2014)
Toledo-Hernandez, R., Rico-Ramirez, V., Iglesias-Silva, G.A., Diwekar, U.M.: A fractional calculus approach to the dynamic optimization of biological reactive systems. Part I: fractional models for biological reactions. Chem. Eng. Sci. 117, 217–228 (2014)
Wang, Z., Yang, D., Ma, T., Sun, N.: Stability analysis for nonlinear fractional-order systems based on comparison principle. Nonlinear Dyn. 75(1–2), 387–402 (2014)
Baleanu, D., Rezapour, S., Saberpour, Z.: On fractional integro-differential inclusions via the extended fractional Caputo–Fabrizio derivation. Bound. Value Probl. 2019(1), 79 (2019)
Baleanu, D., Etemad, S., Rezapour, S.: A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions. Bound. Value Probl. 2020(1), 1 (2020)
Aydogan, M.S., Baleanu, D., Mousalou, A., Rezapour, S.: On high order fractional integro-differential equations including the Caputo–Fabrizio derivative. Bound. Value Probl. 2018(1), 90 (2018)
Baleanu, D., Mousalou, A., Rezapour, S.: On the existence of solutions for some infinite coefficient-symmetric Caputo–Fabrizio fractional integro-differential equations. Bound. Value Probl. 2017(1), 1 (2017)
Ullah, I., Ahmad, S., Rahman, M., Arfan, M.: Investigation of fractional order Tuberculosis (TB) model via caputo derivative. Chaos Solitons Fractalss (2020, in press)
Algahtani, O.J.J.: Comparing the Atangana–Baleanu and Caputo–Fabrizio derivative with fractional order: Allen Cahn model. Chaos Solitons Fractals 89, 552–559 (2016)
Saad, K.M.: Comparing the Caputo, Caputo–Fabrizio and Atangana–Baleanu derivative with fractional order: fractional cubic isothermal auto-catalytic chemical system. Eur. Phys. J. Plus 133(3), 1–12 (2018)
Goufo, E.F.D.: Application of the Caputo–Fabrizio fractional derivative without singular kernel to Korteweg–de Vries–Burgers equation. Math. Model. Anal. 21(2), 188–198 (2016)
Ahmad, B., Alsaedi, A., Nazemi, S.Z., Rezapour, S.: Some existence theorems for fractional integro-differential equations and inclusions with initial and non-separated boundary conditions. Bound. Value Probl. 2014(1), 249 (2014)
Rezapour, S., Samei, M.E.: On the existence of solutions for a multi-singular pointwise defined fractional q-integro-differential equation. Bound. Value Probl. 2020(1), 1 (2020)
Baleanu, D., Rezapour, S., Mohammadi, H.: Some existence results on nonlinear fractional differential equations. Philos. Trans. R. Soc. A, Math. Phys. Eng. Sci. 371(1990), 20120144 (2013)
Haq, F., Shah, K., Rahman, G., Shahzad, M.: Numerical solution of fractional order smoking model via Laplace Adomian decomposition method. Alex. Eng. J. 57(2), 1061–1069 (2018)
Ali, A., Shah, K., Khan, R.A.: Numerical treatment for traveling wave solutions of fractional Whitham–Broer–Kaup equations. Alex. Eng. J. 57(3), 1991–1998 (2018)
Kiymaz, O.: An algorithm for solving initial value problems using Laplace Adomian decomposition method. Appl. Math. Sci. 3(29–32), 1453–1459 (2009)
Khuri, S.A.: A Laplace decomposition algorithm applied to a class of nonlinear differential equations. J. Appl. Math. 1(4), 141–155 (2001)
Shaikh, A., Tassaddiq, A., Nisar, K.S., Baleanu, D.: Analysis of differential equations involving Caputo–Fabrizio fractional operator and its applications to reaction-diffusion equations. Adv. Differ. Equ. 2019(1), 178 (2019)
Losada, J., Nieto, J.J.: Properties of a new fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 1, 87–92 (2015)
Khan, S.A., et al.: Existence theory and numerical solutions to smoking model under Caputo–Fabrizio fractional derivative. Chaos, Interdiscip. J. Nonlinear Sci. 29(1), 013128 (2019)
Zhao, X.Q.: The theory of basic reproduction ratios. In: Dynamical Systems in Population Biology, pp. 285–315. Springer, Cham (2017)
Shah, K., Khalil, H., Khan, R.A.: Analytical solutions of fractional order diffusion equations by natural transform method. Iran. J. Sci. Technol., Trans. A, Sci. 42(3), 1479–1490 (2018)
Acknowledgements
Not applicable.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
All authors contributed equally in writing of this manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ahmad, S., Ullah, R. & Baleanu, D. Mathematical analysis of tuberculosis control model using nonsingular kernel type Caputo derivative. Adv Differ Equ 2021, 26 (2021). https://doi.org/10.1186/s13662-020-03191-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13662-020-03191-x