Research Article | | Peer-Reviewed

Mathematical Modelling and Optimal Control of Tuberculosis Transmission Dynamics in Werder District, Somali Regional State, Ethiopia

Received: 18 January 2025     Accepted: 13 June 2025     Published: 16 January 2026
Views:       Downloads:
Abstract

Tuberculosis (TB) is an infectious disease caused by mycobacterium tuberculosis bacteria (tubercle bacilli) and spreads through the air from an infectious person to a susceptible person. This research aims to investigate the transmission dynamics of tuberculosis (TB) using a mathematical model with optimal control and to estimate the reproduction number in the Werder district of the Somali regional state. To achieve the objective of the study, we employed different mathematical tools and available cumulative case data from the Werder district of the Somali regional state health officials. The model is developed using a system of nonlinear ordinary differential equations and basic reproduction number is determine to perform analysis for both local and global stability of the equilibrium points of the model. In addition to this, sensitivity analysis is studied based on the reproduction number to reveal influential parameters of transmission dynamics of tuberculosis. Then the optimal control strategy was found by minimizing the number of exposed and infected individual populations with tuberculosis taking into account the cost of implementation. The existence of optimal control and characterization was established with the help of Pontryagin's Maximum Principle and we conducted different numerical simulation cases to see agreements with analytical results. Moreover, the model was calibrated with real data collected from the Werder district of the Somali regional state health office from 2016 – 2022. The simulation results also clearly shown that the tuberculosis can be minimized by applying TB prevention programs to susceptible individuals and preventing the failure of treatment in active tuberculosis individuals through continuous supervision and helping patients during treatment period with optimal implementation cost. Furthermore, the optimal control model also predicted that tuberculosis is under control in Werder district.

Published in Science Development (Volume 7, Issue 1)
DOI 10.11648/j.scidev.20260701.11
Page(s) 1-22
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2026. Published by Science Publishing Group

Keywords

Tuberculosis (TB), Differential Equations, Basic Reproduction Number, Stability, Optimal Control

References
[1] Adebisi, B., Yusuff. A and Adebayo. M., 2020. Burden of tuberculosis and challenges facing its Eradication in West Africa.
[2] Agyeman. D., Akosua. A and Richard. O., 2017. Tuberculosis-an overview. Journal of Pub Health and Emergency.
[3] Alene. A., Kefyalew. Aand Alachew. A., 2019. Mapping tuberculosis treatment outcomesin Ethiopia. BMC Infectious diseases.
[4] Alebachew, Z., 2011. Ethiopian National TB Prevalence survey 2010‐2011 PreliminaryResult. Available from:
[5] Annabel. B., Anna. D and Hannah. M., 2021. Global tuberculosis report 2021. Geneva: World Health Organization.
[6] Athithan, S., and Mini Ghosh. M., 2013. Mathematical modelling of TB with the effects ofCase Detection and treatment. International Journal of Dynamics and Control.
[7] Banerjee. K and Sandip. S., 2014. Mathematical modeling: models, analysis andApplications. CRC Press.
[8] Castillo. C and Baojun. S., 2004. Dynamical models of tuberculosis and their applications. Mathematical Biosciences and Engineering.
[9] CDC, 2022. Tuberculosis - Disease Surveillance Epidemiology Program
[10] CSA, 2007. Population and Housing Census Report: Somali Region. Available from:
[11] Daniel. M., 2006. The history of tuberculosis. Respiratory medicine, 100(11): 1862-1870.
[12] Egonmwan. A and Okuonghae. D., 2019. Analysis of a mathematical model forTuberculosis with Diagnosis. Journal of Applied Mathematics and Computing.
[13] Elshafie. M., 2017. Review of the access to tuberculosis care in Doolo Zone in SomaliRegion State of Ethiopia.
[14] Fatmawati. J and Ahmadin. A., 2018. Mathematical modelling of Tuberculosis in aLogistically growing human population with optimal control. World Journal ofModelling and Simulation.
[15] Fleming. W., Wendell. H., and Raymond. W., 2012. Deterministic and stochastic optimalcontrol. Springer Science and Business Media.
[16] FMOH., 2015. Annual TBL Bulletin 2015.
[17] FMOH., 2020. Health Sector Strategic Plan III, 2015 - 2020.
[18] Gaff. A., Holly. J and Elsa. S., 2009. Optimal control applied to vaccination and treatment.
[19] Strategies for various epidemiological models. Mathematical Biosciences and Engineering.
[20] Gugsa. N., Abebe. M and Desalew. T., 2018. Knowledge, Attitude And Practice OnPrevention And Control Of Tuberculosis Among Bilida Kebele Residents, Mana Woreda, Jimma Zone, Oromia Region, South West Ethiopia.
[21] Herzog., 1998. History of tuberculosis. Respiration.
[22] Kelemu. M., Ashenafi. B and Peter. J., 2019. Modeling the Effects of Vaccination andTreatment On Tuberculosis Transmission Dynamics. Journal of Applied Mathematics.
[23] Kereyu, D., & Demie, S., 2021. Transmission dynamics model of Tuberculosis with optimal control strategies in Haramaya district, Ethiopia. Advances in Difference Equations, 2021(1).
[24] Kokomo, E. and Emvudu, Y. (2016). Mathematical modeling: Proposal of a generalmethodology and application of this methodology to epidemiology. Academic Journal of Applied Mathematical Sciences, 2(11): 147–156.
[25] Liu. J. and Zhang. T., 2011. Global stability for a tuberculosis model. Mathematical and Computer Modelling.
[26] Julie. H and Marias. J., 2019. Multidrug-resistant tuberculosis infection and disease in children: a review of new and repurposed drugs. National center for biotechnology information.
[27] Malik. M., Monica. L and Aldila. D., 2018 Mathematical modeling and numerical Simulation of Tuberculosis spread with diabetes effect. IOP Conf. Series: Journal of Physics.
[28] MOFED., 2022. Federal Government of Ethiopia Ministry of Finance and EconomicDevelopment Public Finance Management Assessment Somali RegionalGovernment.. Available from:
[29] Nkamba. A., Leontine. S and Nkague. M., 2019. Mathematical model to assess vaccinationand Effective contact rate impact in the spread of tuberculosis. Journal of Biological Dynamics.
[30] Raviglione. M., Mario. A and Giorgia. S., 2016. Tuberculosis 2015: burden, challengesand Strategy for control and elimination. Infectious disease reports.
[31] Revelle. L., Charles. S., Walter. R and Feldmann. F., 1967. Mathematical models for the economic Allocation of tuberculosis control activities in developing nations. American review of Respiratory disease.
[32] Reves. R and Sahil. A., 2016. Ethiopia moves toward tuberculosis elimination, success Requires Higher investment. Washington DC: Center for Strategic and International Studies.
[33] Ricks, M., Philip. M and Richard. O., 2009. Estimating the burden of tuberculosis among foreign-born persons Acquired prior to entering the US.
[34] Rosa. S and Torres. D., 2018. Optimal control and sensitivity analysis of a fractional order TB Model.
[35] Shirima. S., Theresia. J., Ismail. I and Kuznetsov. D., 2020. Modeling the Transmission Dynamics of Bovine Tuberculosis. International Journal of Mathematics andMathematical Sciences.
[36] Silva. H and Torres. M., 2017. Optimal control of a tuberculosis model with state and Control Delays. Math. Biosci.
[37] Silva. J and Torres. M., 2015. Optimal control of tuberculosis: a review. Dynamics, Games and Science. Springer, Cham.
[38] SRS Bureau of Health., 2016. Somali Adjusted Woreda overview July final 2016.
[39] Tahir. M., Ali. Sh and Gul. Z., 2019. Prevention strategy for super infection mathematical model Tuberculosis and HIV associated with AIDS. Cogent Mathematics] and Statistics.
[40] Turner. S., Richard. D and Graham. H., 2015. Cough and the transmission of tuberculosis. The Journal of infectious diseases.
[41] Ullah. S., Zhien. M., Liu. L., Yang. Y., Jianquan. L., 2019. Modeling and analysis of Tuberculosis (TB) in Khyber Pakhtunkhwa, Pakistan. Mathematics and Computers in Simulation.
[42] Waaler. H., Geser. A and Andersen. S., 1962. The Use of Mathematical Models in the Study of The Epidemiology of Tuberculosis. American Journal of Public Healthand the Nation Health.
[43] Wallis. K and Robert. S., 2016. Mathematical models of tuberculosis reactivation and relapse. Frontiers in Microbiology.
[44] WHO, 2020. Global Tuberculosis report 2020 Geneva. WHO Report.
[45] WHO, 2022. Global tuberculosis report 2021 Geneva. WHO Report.
[46] WHO, 2021. Media Centre-Visual Impairment and Blindness Fact Sheet.
[47] WHO, 2020. Guidelines for treatment of drug-susceptible tuberculosis and patient care.
[48] Woldeyohannes. D., Sisay, S., Mengistu, B., Kassa, H., 2015. Directly observed treatment short course (DOTS) for treatment of new tuberculosis cases in Somali Regional State, Eastern Ethiopia: ten years retrospective study. BMC Res Notes.
[49] Yang. D., Yali. M and Liu, L., 2010. Global stability of two models with incomplete treatment for tuberculosis. Chaos, Solitons and Fractals, Elsevier.
[50] Yohanes. F and Temesgen. T. 2021. Post exposure vaccination for the dynamics of Tuberculosis with fast and slow latent period: A Mathematical model analysis.
[51] Merkin, David R. Introduction to the Theory of Stability. Vol. 24. Springer Science & Business Media, 2012.
[52] La Salle, Joseph P. The stability of dynamical systems. Society for Industrial and Applied Mathematics, 1976.
[53] Van den Driessche, Pauline, and James Watmough. “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission.” Mathematical biosciences 180.1-2 (2002): 29-48.
[54] Wang, Yi, and Jinde Cao. “Global dynamics of multi-group SEI animal disease models with indirect transmission.” Chaos, Solitons & Fractals 69 (2014): 81-89.
[55] Wang, Yi, and Jinde Cao. “Global stability of general cholera models with nonlinear incidence and removal rates.” Journal of the Franklin Institute 352.6 (2015): 2464-2485.
[56] Gaff, H. and Schaefer, E., 2009. Optimal control applied to vaccination and treatment strategies for various epidemiological models. Math. Biosci. Eng, 6(3), pp. 469-492.
[57] Coddington, E. A. and Levinson, N., Tata McGraw-Hill Education; 1955. Theory of ordinary differential equations.
[58] Pontryagin, L. S. and Boltyanskii, V. G., 1962. R. V. Gamkrelidze, and E. F. Mishchenko,“The Mathematical Theory of Optimal Processes,” translated by K. N. Trirogoff, Interscience. New York.
[59] Chong, E. K., Lu, W. S. and Żak, S. H., 2023. An introduction to optimization. John Wiley & Sons.
[60] Lenhart, S. and Workman, J. T., 2007. Optimal control applied to biological models. CRC press.
[61] Rodrigues, H. S. F., 2012. Optimal control and numerical optimization applied to epidemiological models (Doctoral dissertation, Universidade de Aveiro. (Portugal).
[62] Parvaneh. E., Baba, A., Abdullahi. I., Rabiu. A 2020. Analysis of tuberculosis model withSaturated incidence rate and optimal control. Statistical Mechanicsand its Applications.
Cite This Article
  • APA Style

    Taddese, T. A., Silase, T. M. G., Tessema, T. M., Adem, S. D. (2026). Mathematical Modelling and Optimal Control of Tuberculosis Transmission Dynamics in Werder District, Somali Regional State, Ethiopia. Science Development, 7(1), 1-22. https://doi.org/10.11648/j.scidev.20260701.11

    Copy | Download

    ACS Style

    Taddese, T. A.; Silase, T. M. G.; Tessema, T. M.; Adem, S. D. Mathematical Modelling and Optimal Control of Tuberculosis Transmission Dynamics in Werder District, Somali Regional State, Ethiopia. Sci. Dev. 2026, 7(1), 1-22. doi: 10.11648/j.scidev.20260701.11

    Copy | Download

    AMA Style

    Taddese TA, Silase TMG, Tessema TM, Adem SD. Mathematical Modelling and Optimal Control of Tuberculosis Transmission Dynamics in Werder District, Somali Regional State, Ethiopia. Sci Dev. 2026;7(1):1-22. doi: 10.11648/j.scidev.20260701.11

    Copy | Download

  • @article{10.11648/j.scidev.20260701.11,
      author = {Tasew Ayele Taddese and Tsega Mhretab Gebre Silase and Tsegaw Muche Tessema and Sadam Dawud Adem},
      title = {Mathematical Modelling and Optimal Control of Tuberculosis Transmission Dynamics in Werder District, Somali Regional State, Ethiopia},
      journal = {Science Development},
      volume = {7},
      number = {1},
      pages = {1-22},
      doi = {10.11648/j.scidev.20260701.11},
      url = {https://doi.org/10.11648/j.scidev.20260701.11},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.scidev.20260701.11},
      abstract = {Tuberculosis (TB) is an infectious disease caused by mycobacterium tuberculosis bacteria (tubercle bacilli) and spreads through the air from an infectious person to a susceptible person. This research aims to investigate the transmission dynamics of tuberculosis (TB) using a mathematical model with optimal control and to estimate the reproduction number in the Werder district of the Somali regional state. To achieve the objective of the study, we employed different mathematical tools and available cumulative case data from the Werder district of the Somali regional state health officials. The model is developed using a system of nonlinear ordinary differential equations and basic reproduction number is determine to perform analysis for both local and global stability of the equilibrium points of the model. In addition to this, sensitivity analysis is studied based on the reproduction number to reveal influential parameters of transmission dynamics of tuberculosis. Then the optimal control strategy was found by minimizing the number of exposed and infected individual populations with tuberculosis taking into account the cost of implementation. The existence of optimal control and characterization was established with the help of Pontryagin's Maximum Principle and we conducted different numerical simulation cases to see agreements with analytical results. Moreover, the model was calibrated with real data collected from the Werder district of the Somali regional state health office from 2016 – 2022. The simulation results also clearly shown that the tuberculosis can be minimized by applying TB prevention programs to susceptible individuals and preventing the failure of treatment in active tuberculosis individuals through continuous supervision and helping patients during treatment period with optimal implementation cost. Furthermore, the optimal control model also predicted that tuberculosis is under control in Werder district.},
     year = {2026}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Mathematical Modelling and Optimal Control of Tuberculosis Transmission Dynamics in Werder District, Somali Regional State, Ethiopia
    AU  - Tasew Ayele Taddese
    AU  - Tsega Mhretab Gebre Silase
    AU  - Tsegaw Muche Tessema
    AU  - Sadam Dawud Adem
    Y1  - 2026/01/16
    PY  - 2026
    N1  - https://doi.org/10.11648/j.scidev.20260701.11
    DO  - 10.11648/j.scidev.20260701.11
    T2  - Science Development
    JF  - Science Development
    JO  - Science Development
    SP  - 1
    EP  - 22
    PB  - Science Publishing Group
    SN  - 2994-7154
    UR  - https://doi.org/10.11648/j.scidev.20260701.11
    AB  - Tuberculosis (TB) is an infectious disease caused by mycobacterium tuberculosis bacteria (tubercle bacilli) and spreads through the air from an infectious person to a susceptible person. This research aims to investigate the transmission dynamics of tuberculosis (TB) using a mathematical model with optimal control and to estimate the reproduction number in the Werder district of the Somali regional state. To achieve the objective of the study, we employed different mathematical tools and available cumulative case data from the Werder district of the Somali regional state health officials. The model is developed using a system of nonlinear ordinary differential equations and basic reproduction number is determine to perform analysis for both local and global stability of the equilibrium points of the model. In addition to this, sensitivity analysis is studied based on the reproduction number to reveal influential parameters of transmission dynamics of tuberculosis. Then the optimal control strategy was found by minimizing the number of exposed and infected individual populations with tuberculosis taking into account the cost of implementation. The existence of optimal control and characterization was established with the help of Pontryagin's Maximum Principle and we conducted different numerical simulation cases to see agreements with analytical results. Moreover, the model was calibrated with real data collected from the Werder district of the Somali regional state health office from 2016 – 2022. The simulation results also clearly shown that the tuberculosis can be minimized by applying TB prevention programs to susceptible individuals and preventing the failure of treatment in active tuberculosis individuals through continuous supervision and helping patients during treatment period with optimal implementation cost. Furthermore, the optimal control model also predicted that tuberculosis is under control in Werder district.
    VL  - 7
    IS  - 1
    ER  - 

    Copy | Download

Author Information
  • Department of Mathematics, Kebri Dehar University, Kebri Dehar, Ethiopia

  • Department of Mathematics, Kebri Dehar University, Kebri Dehar, Ethiopia

  • Department of Mathematics, Kebri Dehar University, Kebri Dehar, Ethiopia

  • Department of Mathematics, Jigjiga University, Jigjiga, Ethiopia

  • Sections