Research Article  Open Access
Olumuyiwa James Peter, "Transmission Dynamics of Fractional Order Brucellosis Model Using Caputo–Fabrizio Operator", International Journal of Differential Equations, vol. 2020, Article ID 2791380, 11 pages, 2020. https://doi.org/10.1155/2020/2791380
Transmission Dynamics of Fractional Order Brucellosis Model Using Caputo–Fabrizio Operator
Abstract
In this paper, a noninteger order Brucellosis model is developed by employing the Caputo–Fabrizio noninteger order operator. The approach of noninteger order calculus is quite new for such a biological phenomenon. We establish the existence, uniqueness, and stability conditions for the model via the fixedpoint theory. The initial approachable approximate solutions are derived for the proposed model by applying the iterative Laplace transform technique. Finally, numerical simulations are conducted for the analytical results to visualize the effect of various parameters that govern the dynamics of infection, and the results are presented using plots.
1. Introduction
Brucellosis is a contagious zoonotic infection that is attributable to diverse species of Brucella. Although the disease is well contained in Australia and the UK, the annual global incidence of Brucellosis is estimated above 500 000 [1, 2]. Brucellosis is endemic in the Caribbean, South America, Mediterranean Europe, Africa, Asia, and in the Middle East and poses considerable economic consequences in some of the regions despite the enforcement of firm veterinary hygiene interventions [3]. Typically, animal Brucellosis is attributed to either direct or indirect contacts with the Brucellosis bacteria [4]. While direct contact is from contact with the infected secretions, droppings, or vertical transmission from mothertooffspring during delivery, indirect contact occurs through exposure to contaminated environment or objects. It may also occur through respiratory or digestive tracks especially during feeding [5]. Human Brucellosis, on the contrary, is triggered by contamination from infectious animals or their products. Orchitis, epididymitis, placentitis, and abortions are the commonest consequences of Brucellosis [6]. The primary hosts for Brucellosis are sheep, goats, cattle, buffalo, swine, horses, dogs, and man [1]. Control strategies for Brucella within domestic animals involve detection, inoculation, and elimination of the affected livestock [7]. However, rapid control of Brucellosis is achieved by the combination of the interventions.
The quantification of intricate and complex biological systems in terms of mathematical structures has posed a serious challenge to contemporary scientists worldwide. A good number of models with soliton solutions that are designed in terms of integerorder derivatives had been extensively studied in the last two decades [8–12]. Some of these models in the literature include Tzizeicatype equation [13], Fokas–Lenells equation [14], complex Ginzburg–Landau equation [15], differentialdifference equation [16], nonlinear oscillator with discontinuities [17], and Bernstein wavelet equations [18]. The aforementioned models play significant roles in advancing suitable control procedures in a biological system [19].
Scientific and technical models of Brucellosis dynamics can be used to interpret the primary outcomes and condense the discoveries to the best advantage of protection at various stages; still, a number of those models in the form of nonfractional derivative were thoroughly investigated [20–23]. Researchers in mathematics have been engaged in the last decade in the formulation of mathematical models of Brucellosis to investigate the stability and dynamical behaviours of the disease [24, 25]. The primary model of these kinds, which examine the indirect transmission mechanism of Brucellosis with the influence of immunology, was jointly investigated by Zhang et al. [26]. In [1], Li et al. proposed a model that quantified the spread of Brucellosis from sheep to man and among sheep. Unlike in [1], Brucellosis models that described the transmission of Brucellosis in the cattle population and that characterized stillbirths in dairy cattle are designed and analyzed by Tumwiine and Robert [27] and Nie et al. [6], respectively.
As time goes by, advancing mathematical models using noninteger order became a significant area of study since the evolutionrelated realities and evidence are represented more effectively in terms of noninteger calculus [28, 29]. The elaborate investigations of such models play an exemplary role through advancing implementations of various reallife problem. A large number of logical and numerical procedures and schemes are proposed to define arbitrary and real solutions with noninteger operators for differential equations have been playing an exceptional task in the advancement of implementations of numerous problems. A good number of logical and numerical procedures and schemes are proposed to figure out indefinite and actual solutions of the differential equations with noninteger operators [30–35]. However, up till now, there is no work which has been designed to analyze the dynamics of Brucellosis with the fractionalorder calculus.
There are four compartments into which the Brucellosis model in [26] is separated, denoted as susceptible dairy cows S(t), the exposed dairy cow E(t), infectious dairy cows I(t), and the Brucella contaminated environment B(t). The sum of the population at any time t is represented by N = S + E + I + B. In [26], Zhang et al. proposed a Brucellosis model and analyzed the qualitative characteristics of the suggested model which is denoted bywhere denotes the yearly recruitment rate of dairy cows, the yearly birth rate of cow is p, the annual culling rate is d, the treatment outcome rate is , the rate of yearly elimination for the positive dairy cows is , the yearly number of Brucella is n, the yearly rate of natural death of Brucella is , the rate of sterilization in a disinfection is c, and the yearly number of disinfections is h, and is the indirect transmission rate and is the pathogen densitydependent component.
Motivated by the previously mentioned literature, we have applied Caputo–Fabrizio fractional order , where to explore dynamics of the Brucellosis model discussed in Zhang et al. [26]:with initial conditions
The typical model in [26] is derived for . As the Caputo–Fabrizio fractional order quantifies the previously mentioned phenomenon accurately. The numerical results are obtained by employing the Laplace transform method. Arbitrary values are allocated to the initial conditions and parameters to validate our results.
The remainder of this script is arranged as follows. In Section 2, some fundamental remarks and definitions associated with noninteger calculus are stated. In Section 3, existence and uniqueness conditions are derived by using the fixedpoint theory. In Section 4, stability results are obtained via iterative Laplace transform technique. In Section 5, the proposed method is employed to noninteger order Brucellosis model and the results of the simulations are illustrated graphically. Section 6 is devoted to the discussions of numerical results. Finally, the conclusion of the work is stated in Section 7.
2. Preliminaries
Definition 1. (see [36]). Let then the Caputo fractional operator is defined aswhere is the normalization function which depends on , such that
Definition 2. (see [37]). The Caputo–Fabrizio fractional integral operator of order is defined as
Definition 3. (see [36]). The Laplace transformation for the Caputo fractional operator of order and is expressed asGenerally, we have
3. Existence and Uniqueness
This section deals with the existence and uniqueness of the fractional Brucellosis model using fixedpoint theory [38–45]. By applying the fractional integral operator (2) into (4), we have
For clarity, we consider the following kernels:
Theorem 1. The kernels , , , and in (8) satisfy the Lipschitz condition, provided the following inequality holds:
Proof. Let and for the kernel given , and for the kernel , and for the kernel , and and for the kernel be the respective functions associated with the following:By applying the properties of the norm on (10), respectively, we havewhere are bounded functions. Subsequently, we havewhere
We construct the following recursive formula:Furthermore, by applying triangular inequality, we have, , , and
The kernels , , , and which satisfy the Lipschitz condition are therefore obtained:This proves the result.
Theorem 2. The SEIB Brucellosis model involving the fractional operator in (2) has a solution (existence of the solution).
Proof. Following the result in (15) and with respect to the recursive expression, we get the following:Hence, (16) exists, additionally, we show that (16) is the solution of (2) with the following assumption:where the remaining terms of the solution are , and . Hence, we haveWe make use of norm and Lipschitz condition to obtainUsing in (19), we haveTherefore,In the same way, as , we obtain the following.
Using the same approach, in (21), we haveSimilarly, from (21), there exists a solution of (2).
Theorem 3. The SEIB Brucellosis model involving the fractional operator in (2) has a unique solution.
Proof. Proof We adopt another solution (2) say , and ; then, we haveBy applying the norm on (23),By considering Theorems 1 and 2, we obtain the following results:The solution functions in (19) satisfy the following inequalities:Going by the last equation, we conclude that
4. Stability
This section deals with the application of iterative Laplace transformation technique on the fractional Brucellosis model, obtaining the stability condition for the approximate solution.
4.1. Iterative Laplace Transformation Technique
We consider the Brucellosis model in (2) with the initial conditions in (3). Applying the Laplace transformation to both sides of equation (2), we have
By rearranging the following, we have
The inverse Laplace transform of (29) gives
We obtain the following infinite series solution for the SEIB Brucellosis model in (2):
5. Description
In this section, we analyze the numerical simulations of the Caputo–Fabrizio Brucellosis model (2) with the following initial conditions and for specific fractionalorder value The parameters used and their values are given in Table 1.

6. Numerical Results and Discussion
Figures 1–4 display the numerical results for the susceptible cows S(t), exposed cows E(t), infectious cows I(t), and the Brucella contaminated environment B(t) obtained for different values of ( = 1, 0.9, 0.8, 0.7) with different initial conditions and by utilizing the iterative Laplace transform technique of noninteger order Brucellosis model which confirms that the iterative Laplace transform method is promising and can predict the behaviour of the variables accurately within the region of consideration. The integer order has a little impact on the dynamics of the Brucellosis model. This negligible effect has been taken care of in our proposed fractionalorder model and our simulations have shown that a change in the values of the order influences the transmission dynamics of the Brucellosis model. Figure 1 shows that susceptible cows S(t) increase steadily in the first six months when = 0.8 but falls afterwards. However, the population of the susceptible cows falls continuously at other various fractional values of . The plot in Figure 2 of exposed cows E(t) demonstrates that the population of the exposed cow rises rapidly in the first three months when = 0.8 but drops afterwards. The population of the exposed cows also falls continuously for other fractional values of in Figure 2. The populations of the infectious cows I(t) and the Brucella contaminated environment B(t) depicted in Figures 3 and 4, respectively, also follow similar pattern and behaviour as in Figures 1 and 2. It can be inferred that, as the fractionalorder rises from = 0.5 to = 1, the population of cows infected with Brucellosis asymptotically tends to zero.
Figures 5–8 shows different behaviours of the population of susceptible cows S(t), exposed cows E(t), infectious cows I(t), and the Brucella contaminated environment B(t), where the required simulations have been carried out for t = 12 months while we set the value of the fractionalorder parameter . Furthermore, it is observed that the numerical outcome of the simulation is a function of the noninteger derivative . This shows that the special noninteger operator such as the Caputo–Fabrizio operator is efficient and plays significant roles in inaccurate predictions. Besides, the hybrid feature of Caputo–Fabrizio is strong to capture the complex dynamics and offer reliable predictions.
7. Conclusion
In this paper, we have analyzed the Caputo–Fabrizio noninteger order Brucellosis model. Using fixedpoint theorem, the steady results for the existence and uniqueness of solutions and the stability of the proposed noninteger order Brucellosis model have been derived. The iterative Laplace transform technique is applied to conduct the simulations by using a set of values whose sources are from the literature for the model parameters. The effects of various values for the order of the Caputo–Fabrizio noninteger order derivative Brucellosis model on the transmission dynamics of the disease are revealed by the results of the simulations. The results of the simulations, therefore, confirm the validity and efficiency of the Caputo–Fabrizio noninteger derivative in quantifying dynamics of Brucellosis disease and other related kind of problems.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The author declares that they have no conflicts of interest.
References
 M. Li, G. Sun, G. Sun et al., “Transmission dynamics and control for a Brucellosis model in Hinggan League of Inner Mongolia, China,” Mathematical Biosciences and Engineering, vol. 11, no. 5, pp. 1115–1137, 2014. View at: Publisher Site  Google Scholar
 T. Akhvlediani, C. T. Bautista, N. Garuchava et al., “Epidemiological and clinical features of Brucellosis in the country of Georgia,” PLoS One, vol. 12, no. 1, Article ID e0170376, 2017. View at: Publisher Site  Google Scholar
 B. MariaLaura, F. Vincent, and D. O’Callaghan, “Brucellosis: a worldwide zoonosis,” Current Opinionin Microbiology, vol. 4, no. 1, pp. 58–64, 2001. View at: Publisher Site  Google Scholar
 B. G. Mantur and S. K. Amarnath, “Brucellosis in India—a review,” Journal of Biosciences, vol. 33, no. 4, pp. 539–547, 2008. View at: Publisher Site  Google Scholar
 A. F. Brouwer, M. H. Weir, M. C. Eisenberg, R. Meza, and J. N. S. Eisenberg, “Doseresponse relationships for environmentally mediated infectious disease transmission models,” PLoS Computational Biology, vol. 13, no. 9, Article ID e1005765, 2017. View at: Publisher Site  Google Scholar
 J. Nie, G.Q. Sun, X.D. Sun et al., “Modeling the transmission dynamics of dairy cattle Brucellosis in Jilin Province, China,” Journal of Biological Systems, vol. 22, no. 4, pp. 533–554, 2014. View at: Publisher Site  Google Scholar
 M.T. Li, G.Q. Sun, W.Y. Zhang, and Z. Jin, “Modelbased evaluation of strategies to control Brucellosis in China,” International Journal of Environmental Research and Public Health, vol. 14, no. 3, p. 295, 2017. View at: Publisher Site  Google Scholar
 B. Ghanbari and N. Raza, “An analytical method for soliton solutions of perturbed Schrödinger’s equation with quadraticcubic nonlinearity,” Modern Physics Letters B, vol. 33, no. 3, Article ID 1950018, 2019. View at: Publisher Site  Google Scholar
 B. Ghanbari, “Abundant soliton solutions for the HirotaMaccari equation via the generalized exponential rational function method,” Modern Physics Letters B, vol. 33, no. 9, Article ID 1950106, 2019. View at: Publisher Site  Google Scholar
 B. Ghanbari, A. Yusuf, M. Inc, and D. Baleanu, “The new exact solitary wave solutions and stability analysis for the (2 + 1) $(2 + 1)$ dimensional ZakharovKuznetsov equation,” Advances in Difference Equations, vol. 2019, no. 1, p. 49, 2019. View at: Publisher Site  Google Scholar
 B. Ghanbari and M. Inc, “A new generalized exponential rational function method to find exact special solutions for the resonance nonlinear Schrödinger equation,” The European Physical Journal Plus, vol. 133, no. 4, pp. 133–142, 2018. View at: Publisher Site  Google Scholar
 B. Ghanbari, M. S. Osman, and D. Baleanu, “Generalized exponential rational function method for extended ZakharovKuzetsov equation with conformable derivative,” Modern Physics Letters A, vol. 34, no. 20, Article ID 1950155, 2019. View at: Publisher Site  Google Scholar
 B. Ghanbari, I. Mustafa, M. Inc, and L. Rada, “Solitary wave solutions to the tzitzéica type equations obtained by a new efficient approach,” Journal of Applied Analysis & Computation, vol. 9, no. 2, pp. 568–589, 2019. View at: Publisher Site  Google Scholar
 M. S. Osman and B. Ghanbari, “New optical solitary wave solutions of FokasLenells equation in presence of perturbation terms by a novel approach,” Optik, vol. 175, pp. 328–333, 2018. View at: Publisher Site  Google Scholar
 M. S. Osman, B. Ghanbari, and J. A. T. Machado, “New complex waves in nonlinear optics based on complex GinzburgLandau equation with kerr law nonlinearity,” The European Physical Journal Plus, vol. 134, no. 1, p. 20, 2019. View at: Publisher Site  Google Scholar
 A. Yildirima and T. Ozis, “Solution of singular ivps of LaneEmden type by homotopy perturbation method,” Physics Letters A, vol. 369, no. 12, pp. 70–76, 2007. View at: Publisher Site  Google Scholar
 J.H. He, “The homotopy perturbation method nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 15, no. 1, pp. 87–292, 2004. View at: Publisher Site  Google Scholar
 E. Mahsa, A. Saeid, and B. Esmail, “Numerical solution of fractional partial differential equations with normalized Bernsteinwavelet method,” Applications and Applied Mathematics: An International Journal, vol. 14, no. 2, pp. 890–909, 2019. View at: Google Scholar
 S. S. Amaj and K. Sooppy Nisar, “Transmission dynamics of fractional order typhoid fever model using CaputoFabrizio operator,” Chaos, Solitons and Fractals, vol. 128, pp. 355–365, 2019. View at: Publisher Site  Google Scholar
 X.X. Zhan, C. Liu, G. Zhou et al., “Coupling dynamics of epidemic spreading and information diffusion on complex networks,” Applied Mathematics and Computation, vol. 332, pp. 437–448, 2018. View at: Publisher Site  Google Scholar
 M.T. Li, G.Q. Sun, Y.F. Wu, J. Zhang, and Z. Jin, “Transmission dynamics of a multigroup Brucellosis model with mixed cross infection in public farm,” Applied Mathematics and Computation, vol. 237, no. 1, pp. 382–394, 2014. View at: Publisher Site  Google Scholar
 Q. Hou, X. Sun, J. Zhang, Y. Liu, Y. Wang, and Z. Jin, “Modeling the transmission dynamics of sheep Brucellosis in inner Mongolia autonomous region, China,” Mathematical Biosciences, vol. 242, no. 1, pp. 51–58, 2013. View at: Publisher Site  Google Scholar
 G.Q. Sun and Z.K. Zhang, “Global stability for a sheep Brucellosis model with immigration,” Applied Mathematics and Computation, vol. 246, pp. 336–345, 2014. View at: Publisher Site  Google Scholar
 G.G. Jorge and R. Naulin, “Analysis of a model of bovine Brucellosis using singular perturbations,” Journal of Mathematical Biology, vol. 33, no. 2, pp. 211–223, 1994. View at: Publisher Site  Google Scholar
 B. E. Aïnseba, C. Benosman, and P. Magal, “A model for ovine Brucellosis incorporating direct and indirect transmission,” Journal of Biological Dynamics, vol. 4, no. 1, pp. 2–11, 2010. View at: Publisher Site  Google Scholar
 W. Zhang, J. Zhang, Y.P. Wu, and L. Li, “Dynamical analysis of the SEIB model for Brucellosis transmission to the dairy cows with immunological threshold,” Complexity, vol. 2019, Article ID 6526589, 13 pages, 2019. View at: Publisher Site  Google Scholar
 J. Tumwiine and G. Robert, “A mathematical model for treatment of bovine Brucellosis in cattle population,” Journal of Mathematical Modeling, vol. 5, no. 2, pp. 137–152, 2017. View at: Google Scholar
 A. A. Ayoade, M. O. Ibrahim, and O. Odetunde, “Analytical solution of a fractional order cholera model,” Nigerian Journal of Scientific Research, vol. 17, no. 2, pp. 158–164, 2018. View at: Google Scholar
 I. Ameen, “The effect of vaccination and treatment of measles disease described by a fractional order model,” World J of Modeling and Simulation, vol. 14, no. 1, pp. 30–38, 2018. View at: Google Scholar
 B. R. Sontakke and A. Shaikh, “Approximate solutions of time fractional Kawahara and modified Kawahara equations by fractional complex transform,” Communications in Numerical Analysis, vol. 2016, no. 2, pp. 218–229, 2016. View at: Publisher Site  Google Scholar
 B. Ghanbari and D. Kumar, “Numerical solution of predatorprey model with beddingtondeangelis functional response and fractional derivatives with MittagLeffler kernel,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 29, no. 6, Article ID 063103, 2019. View at: Publisher Site  Google Scholar
 B. Ghanbari and J. F. GómezAguilar, “Modeling the dynamics of nutrientphytoplanktonzooplankton system with variableorder fractional derivatives,” Chaos, Solitons & Fractals, vol. 116, pp. 114–120, 2018. View at: Publisher Site  Google Scholar
 B. R. Sontakke, A. S. Shaikh, and K. S. Nisar, “Approximate solutions of a generalized HirotaSatscima coupled KDV and a coupled MKDV systems with time fractional derivatives,” Malaysian Journal of Mathematical Sciences, vol. 12, no. 2, pp. 173–193, 2018. View at: Google Scholar
 A. Shaikh, A. Tassaddiq, K. S. Nisar, and D. Baleanu, “Analysis of differential equations involving CaputoFabrizio fractional operator and its applications to reactiondiffusion equations,” Advances in Difference Equations, vol. 2019, no. 1, p. 78, 2019. View at: Publisher Site  Google Scholar
 B. Ghanbari, D. Baleanu, and M. Al Qurashi, “New exact solutions of the generalized BenjaminBonaMahony equation,” Symmetry, vol. 11, no. 1, pp. 73–85, 2018. View at: Publisher Site  Google Scholar
 M. Caputo and M. Fabrizio, “A new definition of fractional derivative without singular kernel,” Progress in Fractional Differentiation and Applications, vol. 1, no. 2, pp. 73–85, 2015. View at: Google Scholar
 T. Bashiri, S. M. Vaezpour, and J. J. Nieto, “Approximating solution of fabriziocaputo volterra's model for population growth in a closed system by homotopy analysis method,” Journal of Function Spaces, vol. 2018, Article ID 3152502, 10 pages, 2018. View at: Publisher Site  Google Scholar
 A. Atangana and B. T. Alkahtani, “Analysis of nonhomogenous heat model with new trend of derivative with fractional order,” Chaos Solitons Fractals, vol. 89, pp. 566–571, 2016. View at: Publisher Site  Google Scholar
 A. Atangana and B. Alqahtani, “Analysis of the Keller–Segel model with a fractional derivative without singular kernel,” Entropy, vol. 17, no. 12, pp. 439–4453, 2015. View at: Publisher Site  Google Scholar
 D. Kumar, J. Singh, D. Baleanu, and Sushila, “Analysis of regularized longwave equation associated with a new fractional operator with MittagLeffler type kernel,” Physica A: Statistical Mechanics and Its Applications, vol. 492, pp. 155–167, 2018. View at: Publisher Site  Google Scholar
 M. A. Khan, Z. Hammouch, and D. Baleanu, “Modeling the dynamics of hepatitis E via the caputo–fabrizio derivative,” Mathematical Modelling of Natural Phenomena, vol. 14, no. 3, p. 311, 2019. View at: Publisher Site  Google Scholar
 U. Saif, A. K. Muhammad, F. Muhammad, H. Zakia, and B. Dumitru, “A fractional model for the dynamics of tuberculosis infection using CaputoFabrizio derivative,” Discrete & Continuous Dynamical Systems, vol. 13, no. 3, pp. 975–993, 2020. View at: Publisher Site  Google Scholar
 S. Jagdev, K. Devendra, H. Zakia, and A. Abdon, “A fractional epidemiological model for computer viruses pertaining to a new fractional derivative,” Applied Mathematics and Computation, vol. 316, pp. 504–515, 2018. View at: Publisher Site  Google Scholar
 S. Qureshi and A. Yusuf, “Modeling chickenpox disease with fractional derivatives: from caputo to atanganabaleanu,” Chaos, Solitons & Fractals, vol. 122, pp. 111–118, 2019. View at: Publisher Site  Google Scholar
 S. Qureshi and A. Yusuf, “Fractional derivatives applied to MSEIR problems: comparative study with real world data,” The European Physical Journal Plus, vol. 134, no. 4, p. 171, 2019. View at: Publisher Site  Google Scholar
 J. Nie, G.Q. Sun, X.D. Sun et al., “Modeling the transmission dynamics of dairy cattle Brucellosis in Jilin Province, China,” Journal of Biological Systems, vol. 22, no. 4, pp. 533–554, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Olumuyiwa James Peter. 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.