1
|
Butt AIK, Ahmad W, Rabbani HG, Rafiq M, Ahmad S, Ahmad N, Malik S. Exploring optimal control strategies in a nonlinear fractional bi-susceptible model for Covid-19 dynamics using Atangana-Baleanu derivative. Sci Rep 2024; 14:31617. [PMID: 39738098 PMCID: PMC11686349 DOI: 10.1038/s41598-024-80218-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2024] [Accepted: 11/15/2024] [Indexed: 01/01/2025] Open
Abstract
In this article, a nonlinear fractional bi-susceptible [Formula: see text] model is developed to mathematically study the deadly Coronavirus disease (Covid-19), employing the Atangana-Baleanu derivative in Caputo sense (ABC). A more profound comprehension of the system's intricate dynamics using fractional-order derivative is explored as the primary focus of constructing this model. The fundamental properties such as positivity and boundedness, of an epidemic model have been proven, ensuring that the model accurately reflects the realistic behavior of disease spread within a population. The asymptotic stabilities of the dynamical system at its two main equilibrium states are determined by the essential conditions imposed on the threshold parameter. The analytical results acquired are validated and the significance of the ABC fractional derivative is highlighted by employing a recently proposed Toufik-Atangana numerical technique. A quantitative analysis of the model is conducted by adjusting vaccination and hospitalization rates using constant control techniques. It is suggested by numerical experiments that the Covid-19 pandemic elimination can be expedited by adopting both control measures with appropriate awareness. The model parameters with the highest sensitivity are identified by performing a sensitivity analysis. An optimal control problem is formulated, accompanied by the corresponding Pontryagin-type optimality conditions, aiming to ascertain the most efficient time-dependent controls for susceptible and infected individuals. The effectiveness and efficiency of optimally designed control strategies are showcased through numerical simulations conducted before and after the optimization process. These simulations illustrate the effectiveness of these control strategies in mitigating both financial expenses and infection rates. The novelty of the current study is attributed to the application of the structure-preserving Toufik-Atangana numerical scheme, utilized in a backward-in-time manner, to comprehensively analyze the optimally designed model. Overall, the study's merit is found in its comprehensive approach to modeling, analysis, and control of the Covid-19 pandemic, incorporating advanced mathematical techniques and practical implications for disease management.
Collapse
Affiliation(s)
- Azhar Iqbal Kashif Butt
- Department of Mathematics and Statistics, College of Science, King Faisal University, 31982, Al-Ahsa, Saudi Arabia.
| | - Waheed Ahmad
- Department of Mathematics, GC University, Lahore, Pakistan.
| | - Hafiz Ghulam Rabbani
- Abdus Salam School of Mathematical Sciences, GC University, 68-B, New Muslim Town, Lahore, Pakistan
| | - Muhammad Rafiq
- Department of Mathematics, Namal University, 30KM Talagang Road, Mianwali, 42250, Pakistan
| | - Shehbaz Ahmad
- Abdus Salam School of Mathematical Sciences, GC University, 68-B, New Muslim Town, Lahore, Pakistan
| | - Naeed Ahmad
- Department of Mathematics, Government Murray Graduate College, Sialkot, Pakistan
| | - Saira Malik
- Department of Mathematics, GC University, Lahore, Pakistan
| |
Collapse
|
2
|
He JJ, Aljohani A, Mustafa S, Shokri A, Khalsaraei MM, Mukalazi H. Stability analysis of a nonlinear malaria transmission epidemic model using an effective numerical scheme. Sci Rep 2024; 14:17413. [PMID: 39075079 PMCID: PMC11286927 DOI: 10.1038/s41598-024-66503-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2024] [Accepted: 07/02/2024] [Indexed: 07/31/2024] Open
Abstract
Malaria is a fever condition that results from Plasmodium parasites, which are transferred to humans by the attacks of infected female Anopheles mosquitos. The deterministic compartmental model was examined using stability theory of differential equations. The reproduction number was obtained to be asymptotically stable conditions for the disease-free, and the endemic equilibria were determined. More so, the qualitatively evaluated model incorporates time-dependent variable controls which was aimed at reducing the proliferation of malaria disease. The reproduction number R o was determined to be an asymptotically stable condition for disease free and endemic equilibria. In this paper, we used various schemes such as Runge-Kutta order 4 (RK-4) and non-standard finite difference (NSFD). All of the schemes produce different results, but the most appropriate scheme is NSFD. This is true for all step sizes. Various criteria are used in the NSFD scheme to assess the local and global stability of disease-free and endemic equilibrium points. The Routh-Hurwitz condition is used to validate the local stability and Lyapunov stability theorem is used to prove the global asymptotic stability. Global asymptotic stability is proven for the disease-free equilibrium whenR 0 ≤ 1 . The endemic equilibrium is investigated for stability whenR 0 ≥ 1 . All of the aforementioned schemes and their effects are also numerically demonstrated. The comparative analysis demonstrates that NSFD is superior in every way for the analysis of deterministic epidemic models. The theoretical effects and numerical simulations provided in this text may be used to predict the spread of infectious diseases.
Collapse
Affiliation(s)
- Jian Jun He
- School of Humanities and Law, Gannan University of Science and Technology, Ganzhou, 341000, Jiangxi, People's Republic of China
| | - Abeer Aljohani
- Department of Computer Science, Applied College, Taibah University, Medina, 42353, Saudi Arabia
| | - Shahbaz Mustafa
- Department of Mathematics, Institute of Numerical Sciences, Gomal University, Dera Ismail Khan, 29050, KPK, Pakistan
| | - Ali Shokri
- Department of Mathematics, Faculty of Sciences, University of Maragheh, Maragheh, 83111-55181, Iran
| | | | - Herbert Mukalazi
- Department of Mathematics and Statistics, Kyambogo University, Kampala, Uganda.
| |
Collapse
|
3
|
Yasin MW, Ahmed N, Saeed J, Baber MZ, Ali SM, Akgül A, Muhammad S, Hassani MK, Ali M. Numerical study of diffusive fish farm system under time noise. Sci Rep 2024; 14:14711. [PMID: 38926460 PMCID: PMC11208429 DOI: 10.1038/s41598-024-62304-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 05/15/2024] [Indexed: 06/28/2024] Open
Abstract
In the current study, the fish farm model perturbed with time white noise is numerically examined. This model contains fish and mussel populations with external food supplied. The main aim of this work is to develop time-efficient numerical schemes for such models that preserve the dynamical properties. The stochastic backward Euler (SBE) and stochastic Implicit finite difference (SIFD) schemes are designed for the computational results. In the mean square sense, both schemes are consistent with the underlying model and schemes are von Neumann stable. The underlying model has various equilibria points and all these points are successfully gained by the SIFD scheme. The SIFD scheme showed positive and convergent behavior for the given values of the parameter. As the underlying model is a population model and its solution can attain minimum value zero, so a solution that can attain value less than zero is not biologically possible. So, the numerical solution obtained by the stochastic backward Euler is negative and divergent solution and it is not a biological phenomenon that is useless in such dynamical systems. The graphical behaviors of the system show that external nutrient supply is the important factor that controls the dynamics of the given model. The three-dimensional results are drawn for the various choices of the parameters.
Collapse
Affiliation(s)
- Muhammad Waqas Yasin
- Department of Mathematics and Statistics, The University of Lahore, Lahore, Pakistan
- Department of Mathematics, University of Narowal, Narowal, Pakistan
| | - Nauman Ahmed
- Department of Mathematics and Statistics, The University of Lahore, Lahore, Pakistan.
- Department of Computer Science and Mathematics, Lebanese American University, Beirut, Lebanon.
| | - Jawaria Saeed
- Department of Mathematics and Statistics, The University of Lahore, Lahore, Pakistan
| | | | - Syed Mansoor Ali
- Department of Physics and Astronomy, College of Science, King Saud University, P.O. BOX 2455, 11451, Riyadh, Saudi Arabia
| | - Ali Akgül
- Department of Mathematics, Art and Science Faculty, Siirt University, 56100, Siirt, Turkey.
- Department of Computer Science and Mathematics, Lebanese American University, Beirut, Lebanon.
| | - Shah Muhammad
- Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, 11451, Riyadh, Saudi Arabia
| | | | - Mubasher Ali
- School of Engineering and Digital Arts, University of Kent, Canterbury, Kent, UK
| |
Collapse
|
4
|
Nian F, Li F. Epidemic propagation risk study with effective fractal dimension. RISK ANALYSIS : AN OFFICIAL PUBLICATION OF THE SOCIETY FOR RISK ANALYSIS 2024; 44:612-622. [PMID: 37435956 DOI: 10.1111/risa.14192] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 06/11/2023] [Accepted: 06/15/2023] [Indexed: 07/13/2023]
Abstract
In this article, the risk of epidemic transmission on complex networks is studied from the perspective of effective fractal dimension. First, we introduce the method of calculating the effective fractal dimensionD B ${D}_B$ of the network by taking a scale-free network as an example. Second, we propose the construction method of administrative fractal network and calculate theD B ${D}_B$ . using the classical susceptible exposed infectious removed (SEIR) infectious disease model, we simulate the virus propagation process on the administrative fractal network. The results show that the larger theD B ${D}_B$ is, the higher the risk of virus transmission is. Later, we proposed five parameters P, M, B, F, and D, where P denotes population mobility, M denotes geographical distance, B denotes GDP, F denotesD B ${D}_B$ , and D denotes population density. The new epidemic growth index formulaI = ( P + ( 1 - M ) + B ) ( F + D ) $I = {( {P + ( {1 - M} ) + B} )}^{( {F + D} )}$ was obtained by combining these five parameters, and the validity of I in epidemic transmission risk assessment was demonstrated by parameter sensitivity analysis and reliability analysis. Finally, we also confirmed the reliability of the SEIR dynamic transmission model in simulating early COVID-19 transmission trends and the ability of timely quarantine measures to effectively control the spread of the epidemic.
Collapse
Affiliation(s)
- Fuzhong Nian
- School of Computer and Communication, Lanzhou University of Technology, Lanzhou, P.R. China
| | - Fangfang Li
- School of Computer and Communication, Lanzhou University of Technology, Lanzhou, P.R. China
| |
Collapse
|
5
|
Butt AIK, Imran M, Aslam J, Batool S, Batool S. Computational analysis of control of hepatitis B virus disease through vaccination and treatment strategies. PLoS One 2023; 18:e0288024. [PMID: 37883381 PMCID: PMC10602336 DOI: 10.1371/journal.pone.0288024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2023] [Accepted: 06/16/2023] [Indexed: 10/28/2023] Open
Abstract
Hepatitis B disease is an infection caused by a virus that severely damages the liver. The disease can be both acute and chronic. In this article, we design a new nonlinear SVEICHR model to study dynamics of Hepatitis B Virus (HBV) disease. The aim is to carry out a comprehensive mathematical and computational analysis by exploiting preventive measures of vaccination and hospitalization for disease control. Mathematical properties of proposed model such as boundedness, positivity, and existence and uniqueness of the solutions are proved. We also determine the disease free and endemic equilibrium points. To analyze dynamics of HBV disease, we compute a biologically important quantity known as the reproduction number R0 by using next generation method. We also investigate the stability at both of the equilibrium points. To control the spread of disease due to HBV, two feasible optimal control strategies with three different cases are presented. For this, optimal control problem is constructed and Pontryagin maximum principle is applied with a goal to put down the disease in the population. At the end, we present and discuss effective solutions obtained through a MATLAB code.
Collapse
Affiliation(s)
- Azhar Iqbal Kashif Butt
- Department of Mathematics and Statistics, College of Science, King Faisal University, Al-Ahsa, Saudi Arabia
- Department of Mathematics, GC University, Lahore, Pakistan
| | - Muhammad Imran
- Department of Mathematics, GC University, Lahore, Pakistan
- Tandy School of Computer Science, The University of Tulsa, Tulsa, OK, United States of America
| | - Javeria Aslam
- Department of Mathematics, GC University, Lahore, Pakistan
| | - Saira Batool
- Government Associate College (W) Kamar Mashani, Mianwali, Pakistan
| | - Saira Batool
- Department of Mathematics, GC University, Lahore, Pakistan
| |
Collapse
|
6
|
Khan IU, Mustafa S, Shokri A, Li S, Akgül A, Bariq A. The stability analysis of a nonlinear mathematical model for typhoid fever disease. Sci Rep 2023; 13:15284. [PMID: 37714901 PMCID: PMC10504385 DOI: 10.1038/s41598-023-42244-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 09/07/2023] [Indexed: 09/17/2023] Open
Abstract
Typhoid fever is a contagious disease that is generally caused by bacteria known as Salmonella typhi. This disease spreads through manure contamination of food or water and infects unprotected people. In this work, our focus is to numerically examine the dynamical behavior of a typhoid fever nonlinear mathematical model. To achieve our objective, we utilize a conditionally stable Runge-Kutta scheme of order 4 (RK-4) and an unconditionally stable non-standard finite difference (NSFD) scheme to better understand the dynamical behavior of the continuous model. The primary advantage of using the NSFD scheme to solve differential equations is its capacity to discretize the continuous model while upholding crucial dynamical properties like the solutions convergence to equilibria and its positivity for all finite step sizes. Additionally, the NSFD scheme does not only address the deficiencies of the RK-4 scheme, but also provides results that are consistent with the continuous system's solutions. Our numerical results demonstrate that RK-4 scheme is dynamically reliable only for lower step size and, consequently cannot exactly retain the important features of the original continuous model. The NSFD scheme, on the other hand, is a strong and efficient method that presents an accurate portrayal of the original model. The purpose of developing the NSFD scheme for differential equations is to make sure that it is dynamically consistent, which means to discretize the continuous model while keeping significant dynamical properties including the convergence of equilibria and positivity of solutions for all step sizes. The numerical simulation also indicates that all the dynamical characteristics of the continuous model are conserved by discrete NSFD scheme. The theoretical and numerical results in the current work can be engaged as a useful tool for tracking the occurrence of typhoid fever disease.
Collapse
Affiliation(s)
- Ihsan Ullah Khan
- Department of Mathematics, Institute of Numerical Sciences, Gomal University, Dera Ismail Khan, 29050, KPK, Pakistan
| | - Shahbaz Mustafa
- Department of Mathematics, Institute of Numerical Sciences, Gomal University, Dera Ismail Khan, 29050, KPK, Pakistan
| | - Ali Shokri
- Department of Mathematics, Faculty of Science, University of Maragheh, Maragheh, 83111-55181, Iran
| | - Shuo Li
- School of Mathematics and Data Sciences, Changji University, Changji, 831100, Xinjiang, People's Republic of China
| | - Ali Akgül
- Department of Computer Science and Mathematics, Lebanese American University, Beirut, 5053, Lebanon
- Mathematics Research Center, Department of Mathematics, Near East University, Near East Boulevard, 99138, Nicosia, Mersin, Turkey
- Department of Mathematics, Art and Science Faculty, Siirt University, 56100, Siirt, Turkey
| | - Abdul Bariq
- Department of Mathematics, Laghman University, Mehtarlam City, Laghman, 2701, Afghanistan.
| |
Collapse
|
7
|
Implementation of computationally efficient numerical approach to analyze a Covid-19 pandemic model. ALEXANDRIA ENGINEERING JOURNAL 2023; 69:341-362. [PMCID: PMC9918432 DOI: 10.1016/j.aej.2023.01.052] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 01/05/2023] [Accepted: 01/23/2023] [Indexed: 09/17/2023]
Abstract
Corona virus disease (Covid-19) which has caused frustration in the human community remains the concern of the globe as every government struggles to defeat the pandemic. To deal with the situation, we have extensively studied a deadly Covid-19 model to provide a deep insight into the disease dynamics. A mathematical analysis of the model utilizing preventive measures is performed with the aim to reduce the disease burden. Some comprehensive mathematical techniques are employed to demonstrate several essential properties of solutions. To start with, we proved the existence and uniqueness of solutions. Equilibrium points are stated both in the absence and presence of the pandemic. Biologically important quantity known as threshold parameter is computed to handle the future disease dynamics and analyzed for its sensitivity. We proved the stability of the proposed model at equilibrium points by employing necessary conditions on threshold parameter. A reliable and competitive numerical analysis is conducted to observe the effectiveness of implemented strategies and to verify obtained analytical results. The most sensitive parameters are determined through sensitivity analysis. An important feature of this study is to employ Non-Standard Finite Difference (NSFD) numerical scheme to solve the system instead of other standard methods like Runge–Kutta method of order 4 (RK4). Finally, several numerical simulations are performed to validate our former theoretical analysis. Numerical results exhibiting dynamical behavior of Covid-19 system under the influence of involved parameters suggest that both the implemented strategies, especially quarantine of exposed individuals, are effective for the substantial reduction in the diseased population and to achieve the herd immunity.
Collapse
|
8
|
Theoretical Analysis of a COVID-19 CF-Fractional Model to Optimally Control the Spread of Pandemic. Symmetry (Basel) 2023. [DOI: 10.3390/sym15020380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
In this manuscript, we formulate a mathematical model of the deadly COVID-19 pandemic to understand the dynamic behavior of COVID-19. For the dynamic study, a new SEIAPHR fractional model was purposed in which infectious individuals were divided into three sub-compartments. The purpose is to construct a more reliable and realistic model for a complete mathematical and computational analysis and design of different control strategies for the proposed Caputo–Fabrizio fractional model. We prove the existence and uniqueness of solutions by employing well-known theorems of fractional calculus and functional analyses. The positivity and boundedness of the solutions are proved using the fractional-order properties of the Laplace transformation. The basic reproduction number for the model is computed using a next-generation technique to handle the future dynamics of the pandemic. The local–global stability of the model was also investigated at each equilibrium point. We propose basic fixed controls through manipulation of quarantine rates and formulate an optimal control problem to find the best controls (quarantine rates) employed on infected, asymptomatic, and “superspreader” humans, respectively, to restrict the spread of the disease. For the numerical solution of the fractional model, a computationally efficient Adams–Bashforth method is presented. A fractional-order optimal control problem and the associated optimality conditions of Pontryagin maximum principle are discussed in order to optimally reduce the number of infected, asymptomatic, and superspreader humans. The obtained numerical results are discussed and shown through graphs.
Collapse
|
9
|
Numerical analysis of Atangana-Baleanu fractional model to understand the propagation of a novel corona virus pandemic. ALEXANDRIA ENGINEERING JOURNAL 2022; 61:7007-7027. [PMCID: PMC8692131 DOI: 10.1016/j.aej.2021.12.042] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Revised: 12/14/2021] [Accepted: 12/16/2021] [Indexed: 06/16/2023]
Abstract
In this manuscript, we formulated a new nonlinear SEIQR fractional order pandemic model for the Corona virus disease (COVID-19) with Atangana-Baleanu derivative. Two main equilibrium points F0∗,F1∗ of the proposed model are stated. Threshold parameter R0 for the model using next generation technique is computed to investigate the future dynamics of the disease. The existence and uniqueness of solution is proved using a fixed point theorem. For the numerical solution of fractional model, we implemented a newly proposed Toufik-Atangana numerical scheme to validate the importance of arbitrary order derivative ρ and our obtained theoretical results. It is worth mentioning that fractional order derivative provides much deeper information about the complex dynamics of Corona model. Results obtained through the proposed scheme are dynamically consistent and good in agreement with the analytical results. To draw our conclusions, we explore a complete quantitative analysis of the given model for different quarantine levels. It is claimed through numerical simulations that pandemic could be eradicated faster if a human community selfishly adopts mandatory quarantine measures at various coverage levels with proper awareness. Finally, we have executed the joint variability of all classes to understand the effectiveness of quarantine policy on human population.
Collapse
|
10
|
Masandawa L, Mirau SS, Mbalawata IS, Paul JN, Kreppel K, Msamba OM. Modeling nosocomial infection of COVID-19 transmission dynamics. RESULTS IN PHYSICS 2022; 37:105503. [PMID: 35469342 PMCID: PMC9021122 DOI: 10.1016/j.rinp.2022.105503] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Revised: 03/03/2022] [Accepted: 04/08/2022] [Indexed: 06/14/2023]
Abstract
COVID-19 epidemic has posed an unprecedented threat to global public health. The disease has alarmed the healthcare system with the harm of nosocomial infection. Nosocomial spread of COVID-19 has been discovered and reported globally in different healthcare facilities. Asymptomatic patients and super-spreaders are sough to be among of the source of these infections. Thus, this study contributes to the subject by formulating a S E I H R mathematical model to gain the insight into nosocomial infection for COVID-19 transmission dynamics. The role of personal protective equipment θ is studied in the proposed model. Benefiting the next generation matrix method,R 0 was computed. Routh-Hurwitz criterion and stable Metzler matrix theory revealed that COVID-19-free equilibrium point is locally and globally asymptotically stable wheneverR 0 < 1 . Lyapunov function depicted that the endemic equilibrium point is globally asymptotically stable whenR 0 > 1 . Further, the dynamics behavior ofR 0 was explored when varying θ . In the absence of θ , the value ofR 0 was 8.4584 which implies the expansion of the disease. When θ is introduced in the model,R 0 was 0.4229, indicating the decrease of the disease in the community. Numerical solutions were simulated by using Runge-Kutta fourth-order method. Global sensitivity analysis is performed to present the most significant parameter. The numerical results illustrated mathematically that personal protective equipment can minimizes nosocomial infections of COVID-19.
Collapse
Affiliation(s)
- Lemjini Masandawa
- School of Computational and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology, P.O. Box 447, Arusha, Tanzania
| | - Silas Steven Mirau
- School of Computational and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology, P.O. Box 447, Arusha, Tanzania
| | - Isambi Sailon Mbalawata
- African Institute for Mathematical Sciences, NEI Globla Secretariat, Rue KG590 ST, Kigali, Rwanda
| | - James Nicodemus Paul
- School of Computational and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology, P.O. Box 447, Arusha, Tanzania
| | - Katharina Kreppel
- School of Computational and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology, P.O. Box 447, Arusha, Tanzania
| | - Oscar M Msamba
- Arusha Technical College, P.O. Box 296, Arusha, Tanzania
| |
Collapse
|