1. Introduction
Mathematical models can play an important role in analyzing and studying different epidemiological issues. From epidemiological models, one can predict not only the dynamics of a disease but also the long-term behavior derived from the dynamical nature of early invasions, or, in fact, the influence of immunization on disease transmission. The Susceptible-Infectious-Recovered model, popularly known as the SIR model, is a fundamental epidemiological model used to analyze and predict the spread of infectious diseases within a population [1]. This model partitions the populace into three compartments: people who are not yet infected but are in danger of contamination, people who are already infected with the disease and spread the sickness, and recovered people who have either endured contamination and acquired resistance or have capitulated to the illness.
Though the use of epidemiological models has skyrocketed with the growth of computing, there are a number of drawbacks when applying mathematical models. The SIR model assumes that individuals move through these compartments based on a set of differential equations capturing the dynamics of disease transmission. While the classic SIR model provides valuable insights into the general trends of an epidemic, it may oversimplify the complex spatial and temporal variations that often characterize real-world scenarios. To address this limitation, the concept of patch transfer has been introduced. It is undeniable that the inclusion of population dispersal in an epidemiological model makes it more realistic. Patch transfer extends the traditional SIR model by incorporating spatial heterogeneity into the analysis. Instead of treating the entire population as a homogeneous unit, the region under consideration is divided into patches, each representing a distinct subpopulation with its own set of infection dynamics. This approach allows for a more nuanced understanding of how infectious diseases spread across different geographical areas or social groups.
Incorporating patch transfer into the SIR model involves accounting for the movement of individuals between patches to reflect the mobility patterns observed in real-world populations [2,3,4]. This addition enables researchers to explore how factors such as traveling, migration, and social interactions contribute to the spatial transmission of infectious diseases. By combining the classic SIR model with patch transfer, researchers can gain a more accurate and realistic representation of the dynamics of disease spread in diverse and interconnected populations. This enhanced model is particularly valuable for addressing public health challenges in regions with complex geographical and social structures, providing a more comprehensive foundation for designing effective intervention strategies and mitigating the impact of infectious diseases. A two-patch epidemic system refers to a mathematical model which is used to study the spread and dynamics of infectious diseases among two connected populations. In this system, two populations are connected by migration or movement; the spread of an infectious disease in one population can impact the other population. A two-patch epidemic system acts as a useful tool for studying the transmission dynamics of infectious diseases in various contexts, including the global spread of pandemics, the transmission of diseases between animal and human populations, and the impact of vaccination programs on disease spread [5,6]. Let us consider a real-life scenario of the cross-border spread of infectious diseases. Imagine that an infectious disease outbreak occurs in two neighboring countries, Country A and Country B. These countries have different population densities, healthcare infrastructures, and social behavioral patterns. People frequently travel between these countries for work, tourism, and other reasons. Features of both the countries are as follows:
Country A:
Higher population density
Advanced healthcare facilities
Urban areas with dense social interactions
Country B:
Lower population density
Limited healthcare facilities
Rural areas with less frequent social interactions
The infectious disease can spread within each country, with transmission rates influenced by population density [7] and healthcare accessibility. People move between the two countries, facilitating the cross-border spread of the disease. Individuals infected in Country A may travel to Country B, introducing the disease to new areas. Conversely, individuals from Country B may travel to Country A, contributing to the spread of the disease in more densely populated regions. The model helps public health officials understand how the movement of people between different regions influences the overall spread of the disease. Each nation can develop strategies in order to manage and contain the epidemic while considering the distinct qualities of both urban and rural areas, as well as the effects of cross-border exchanges.
On the other hand, fractional calculus is a branch of mathematics that deals with the generalization of derivatives and integrals to non-integer orders [8]. In recent years, it has been increasingly applied to epidemic models in order to capture the complex dynamics of infectious diseases [9,10]. The use of fractional calculus in epidemic systems has shown several effects, some of which are highlighted below:
Increased accuracy: the use of fractional calculus in epidemic models can provide a more accurate representation of the dynamics of infectious diseases. This is because it can capture long-term memory effects, which are not accounted for in traditional integer-order differential equations.
Increased complexity: fractional calculus introduces a new level of complexity to epidemic models, making them more challenging to analyze. This complexity can lead to the emergence of new behavioral patterns that are not present in traditional models.
Improved control strategies: fractional calculus can help design more effective control strategies for infectious diseases. This is because it can capture the impact of interventions such as vaccination campaigns and social distancing measures on the long-term dynamics of the disease [11].
New insights into disease dynamics: the use of fractional calculus in epidemic models can provide new insights into the dynamics of infectious diseases. For example, it can help in order to explain why some diseases exhibit periodic outbreaks or why they show complex oscillatory behavior.
In real-life examples, the patch transfer concept in the SIR model helps capture the dynamics of an infectious disease spreading across borders while considering the distinct features of different regions and the movement of individuals between them. Using mathematical models, researchers can better understand how diseases spread and how interventions such as vaccination programs or travel restrictions can impact disease transmission. Additionally, the two-patch epidemic system can help policymakers to make informed decisions about how to allocate resources and respond to disease outbreaks in different populations. Overall, the two-patch epidemic system is an important tool for understanding the spread and dynamics of infectious diseases in interconnected populations and has the potential to inform public health interventions and policies aimed at controlling and preventing the spread of infectious diseases.
Today, fractional calculus is becoming increasingly popular for model simulation and for studying the complex nature of disease transmission. The use of fractional calculus in epidemic models can provide a more accurate and nuanced understanding of the dynamics of infectious diseases. It can help design more effective control strategies and offer new insights into the behavior of infectious diseases. Considering all of the above facts, we have constructed our model in the Caputo differential framework. In this context, the works of Das and Samanta [9,12,13] on fractional order dynamics should be noted. Several academics have recently made substantial additions to the numerous SIR models in both integer and fractional order systems.
There are many examples of infectious diseases in which people’s movement has played a key role in the spread of the virus [14,15]. Consider the example of SARS; it had spread over almost all of China and several other regions of the world solely as a result of human movement. Moreover, the world has recently witnessed the pandemic situation caused by the COVID-19 coronavirus. Each and every country imposed strict restrictions on non-pharmaceutical interventions and people’s migration, whether within states or between countries. Numerous articles on COVID-19 have already been published, demonstrating how control measures and public awareness have reduced the disease load [16,17,18]. Epidemic models intended to explain the dynamic behavior of disease propagation between patches have been presented and analyzed by Takeuchi et al. [19], Wang and Wang [20], Wang and Mulone [21], and Wang and Zhao [22,23]. Other articles have been published in which researchers focused on epidemiological models in a two-patch environment under optimal treatment and vaccination [4,24]. Multi-city epidemic models were suggested by Wan and Cui [25] and Arino and Vanden Driessche [26] in order to investigate the dynamics of infectious illnesses. Another transport-related pandemic model with entrance and exit checks was examined by Liu et al. [27]. In his study, Chen proposed an epidemic model with adaptive dispersal rates in order to observe how proper adaptation helps reduce disease load in connected groups or patches [28]. The impact of population dispersal among ‘n’ patches on the transmission of a disease was examined by Arino and Driessche [26]. According to their numerical calculations, a disease that becomes endemic will not become extinct in two patches. It becomes extinct in two patches when people disperse within a certain range. Their simulation demonstrates how population dispersal might help in the eradication of endemic illnesses in isolated patches if the basic reproduction numbers of those patches are not very high and the contact rates between the two patches are not high either. Chen et al. [29] presented an SIR model that includes illnesses due to transport. Despite the publications mentioned above, further study is required on epidemic models involving population dispersal. In the current study, we have examined an epidemic model that accounts for population dispersal across different locations. To keep things simple, we have merely taken into account the population distribution between the two cities. It is considered that susceptible and recovered people are free to move between these two regions, while there are strict restrictions on the people who are diagnosed to be infected. The main aim of our work is to observe how the infection level changes in each of the patches when people disperse at the optimal rate, which is considered to be time-dependent. Moreover, we have not restricted ourselves to the deterministic approach, and have also analyzed the fractional-order system.
We have organized this article into several sections. Model formations in both the integer-order system and fractional system are presented in Section 2, followed by the establishment of the positivity and boundedness of the solutions of the integer-order system. Equilibrium analysis, sensitivity analysis, and stability analysis of the model are studied in Section 4, Section 5 and Section 6, respectively. Next, we have presented our numerical simulations for the system (integer and fractional) in Section 7 and Section 8, respectively. In Section 9, an optimal control problem has been proposed for both fractional and traditional integer-order systems, while the corresponding numerical studies of control problems are analyzed in Section 10 and Section 11, respectively. Finally, we have presented our conclusions in Section 12.
2. Model Formulation
This article deals with a compartmental SIRS model in a two-patch environment where the population in each patch is infected with a disease. The susceptible and recovered populations are allowed to disperse between the patches, while there are strict restrictions for the infected population, who cannot leave their respective patches. It is assumed that the infected individuals are advised to live in their own patches. The total population in Patch-i for is subdivided into the three following compartments in each patch: the susceptible population , infected population , and recovered population (see Figure 1).
It is considered that the interaction between individuals is a mass action type and the population is homogeneously mixed in the environment; hence, the parameters and represent the disease transmission rates from susceptible to infected people in Patch-1 and Patch-2, respectively. The recruitment rate in Patch-1 is denoted by , whereas for Patch-2 it is . The natural mortality rates are parameterized by and , respectively, and are considered in each compartment of each patch. The parameters and denote the disease-induced death rates of Patch-1 and Patch-2, respectively, which are considered in the infected class of each patch only. The infected people gain immunity at a rate in Patch-i (for ) either by natural immunity or by clinical treatment and then move to recovered class; however, as it is considered that the recovery is only temporary, a portion of the recovered people becomes susceptible again in Patch-i with a rate for . The term m represents the dispersal speed of the population between two patches. The parameters and are the probabilities of dispersion of the susceptible population from Patch-2 to Patch-1 and from Patch-1 to Patch-2, respectively, while and are the probabilities of dispersion of the recovered population from Patch-2 to Patch-1 and from Patch-1 to Patch-2, respectively. A diagram in Figure 1 is provided to show the epidemic model in a two-patch environment. It should be noted that in the absence of dispersal, i.e., when , the population in each patch can evolve independently. Thus, the SIRS epidemic model with population-dispersal dynamics is described as follows:
(1)
Next, we consider the above SIRS model under the Caputo fractional order system framework. The main reasons behind our consideration are: (1) observing the impact of the long-term memory effect on the transmission dynamics of the SIRS model; (2) looking for better numerical results; and (3) studying the effect of fractional order to achieve a disease-free environment. When dealing with real-world problems, the Caputo derivative is very useful because it allows traditional initial and boundary conditions to be included in the formulation of the problem; moreover, the derivative of a constant is zero, which is not the case with the Riemann–Liouville fractional derivative. For this reason, we intend to consider the Caputo derivative.
(2)
Here, is the Caputo fractional derivative operator of order . The system (2) is well balanced in the time dimension. For simplicity, we have omitted the power in each parameter containing . The system we intend to study is mentioned below:
(3)
Here, are the initial stages of system (3).
For an absolutely continuous function ), is defined as
where, is the Gamma function, , and n is a natural number. In particular, for ,For , the model (2) is transformed into the model (1).
3. Positivity and Boundedness of System (1)
This section shows the biological relevance of the proposed deterministic system by proving the positivity and boundedness of the system variables with time.
Any solution of system (1) in is positive for .
The right-hand side of the model (1) is continuous and locally Lipschitzian on the domain , which implies that a solution of the system including the initial conditions exists uniquely on , where . The second equation of system (1) with the initial conditions provides
Similarly, we can show that for all . Now, we want to show that . If this does not hold, then such that and . Then, we need to show that . If it does not hold, then such that and . From the last equation of (1), we have
which contradicts So, , while from the third equation we obtain which contradicts Thus, . From the previous steps, it follows that .Next, we want to prove that . If this statement is not true, then such that and . Thus, we need to prove that . If this is not true, then such that and . From the fourth equation of (1), we have
which contradicts So, . Lastly, from the first equation we have which contradicts Thus, . From the above-mentioned steps, we obtain , where , hence, the theorem. □Any solution of system (1) initiating from is bounded with time.
Let us consider the overall population as . Then the time derivative along the solution trajectories for model (1) is obtained as follows:
where and . Now, using the notion of differential inequality for , we can derive the following result: , where is the total population size at the initial time. Thus, for any .The solutions of the system converge in the region □
4. Equilibrium Analysis
In this section, we first derive the basic reproduction numbers in both the presence and absence of population dispersal, then analyze the endemic equilibrium states. In a susceptible environment, the number of newly infected individuals generated from a single infected person is represented by ; in this work, we have followed the process developed by van den Driessche and Watmough [30] to obtain . Let us denote for and for .
Consider the situation when there is no population dispersal between the patches . For , systems (1) and (2) possess a disease-free equilibrium (DFE) where . The basic reproduction numbers for Patch-1 and Patch-2 in this case are found to be and , respectively.
The Patch-2 infection-free equilibrium in absence of population dispersal is noted as , where , while the Patch-1 infection-free equilibrium is , where . This shows that feasibility of and occurs according to and , respectively.
The endemic equilibrium point in absence of population dispersal is , where and . Thus, leads to the feasibility of and , respectively.
Next, we consider the situation in when people can disperse towards Patch-2 from Patch-1 but cannot move back to Patch-1 (i.e., ). In this case, systems (1) and (2) possess a disease-free equilibrium (DFE) , where . Here, we obtain the basic reproduction number as .
The Patch-2 infection-free equilibrium when is noted as , where , while Patch-1 infection-free equilibrium is , where . This shows that the feasibility of and occurs according to and , respectively.
Let us denote, and .
The endemic equilibrium point here is , where , , and .
Here, when and when . Thus, a disease can invade in Patch-2 even if the corresponding reproduction number lies below the unit value ().
In the next situation, people can disperse towards Patch-1 from Patch-2 but cannot move back to Patch-2 (i.e., ). In this case, systems (1) and (2) possess a disease-free equilibrium (DFE) where . We obtain the basic reproduction number as .
The Patch-2 infection-free equilibrium when is noted as , where , while the Patch-1 infection-free equilibrium is , where . This shows that the feasibility of and occurs according to and , respectively.
Let us denote, and .
The endemic equilibrium point here is , where , , and .
Here, when and when . Thus, a disease can invade in Patch-1 even if the corresponding reproduction number lies below the unit value ().
Lastly, we consider the situation where people from both patches can disperse (i.e., ). In this case, systems (1) and (2) possess a disease-free equilibrium (DFE) where .
: To obtain in the presence of population dispersal, we can consider . Then, we have where
Here, and contain the compartment containing the new infection term and the other terms, respectively; thus, at the disease-free equilibrium we have
The spectral radius of the next generation matrix is , and is provided by
(4)
The endemic equilibrium point is , where , , .
Let us consider . Then, we have and .
As , we have when and when .
When the population disperses among two patches, several different scenarios may occur: (IV.a), when Patch-1 shows SIRS epidemic dynamics but Patch-2 contains only the susceptible population; (IV.b), when Patch-2 shows SIRS epidemic dynamics but Patch-1 contains only susceptible population; (IV.c), when Patch-1 shows a SIRS epidemic dynamics but Patch-2 is infection-free; and (IV.d), when Patch-1 shows a SIRS epidemic dynamics but Patch-2 is infection-free.
Patch-2 infection-free equilibrium , where, and .
Here, . Thus, only when .
Here. . Thus, only when .
Altogether, we have come up with equilibrium points for different situations that may occur due to population dispersal, before becoming infected or after recovery. It is observed here that even if the basic reproduction number of one patch becomes less than 1, the whole system may become infected. Moreover, it has been noted prominently from the model (2) that almost all parameters contain as a power which implies that the reproduction number clearly depends on for the fractional-order system.
5. Sensitivity Analysis
The basic reproduction number depends on system parameters such as the disease transmission rates , dispersal speed , probabilities of susceptibles to disperse between patches , etc. However, not all of the parameters can be controlled, as they may depend only on environmental changes. A few parameters have more impact on disease propagation compared to others; thus, in this section the effect of on is analyzed.
Now, , where for and . Let us denote . Then, we have
Now, the normalized forward sensitivity index for a parameter say , is provided as [31]:
Then, the normalized forward sensitivity index for the parameters , and is:It is observed that and maintain monotonically increasing relation with and , respectively which is evident as the disease transmission rates of each patch mainly decide the invasion of infection on that patch. Moreover, the population dispersal speed can also be considered as a contributory factor in disease propagation in the patches. Then it will be the movement probability of the population that will decide the infection fatality. For example, an increased dispersal speed will help to spread the infection in Patch-1 only when holds. On the other hand, if is satisfied, then an increasing dispersal speed will make Patch-2 infected. Moreover, if more people start to disperse towards Patch-1, then the chance of people becoming infected there will increase only and Patch-2 will remain safe. This holds the other way around also. So from the analysis, it is observed that regulation on the population dispersal speed will help to control the prevalence to a certain extent during a disease outbreak.
6. Stability Analysis
The local stability conditions for the equilibrium points of systems (1) and (2) are discussed in this section. The Routh–Hurwitz criterion helps to derive the stability conditions of the system (1), whereas the stability conditions for system (2) follow from the results stated in Theorem A2 of Appendix B. According to this criterion, an equilibrium point is said to be locally asymptotically stable (LAS) if its corresponding Jacobian matrix provides eigenvalues with negative real parts. Let for and let for .
Consider the case when the population at each patch is not allowed to disperse from one patch to another (). Here, Patch-1 and Patch-2 act as independent regions.
In the absence of population dispersal, the disease-free equilibrium of the proposed system is locally asymptotically stable (LAS) for for Patch-1 and for Patch-2.
The Jacobian matrix of systems (1) and (2) in the absence of population dispersal is provided as follows:
This provides the eigenvalues , which are always negative; the other two eigenvalues are and . Thus, and will be negative only when , . Moreover, the arguments of all eigenvalues will be if , . Hence, is locally asymptotically stable for both systems when for Patch-1 and for Patch-2. □
In the absence of population dispersal, the endemic equilibrium of the proposed system is locally asymptotically stable (LAS) for for Patch-1 and for Patch-2.
At , the Jacobian matrix is
The eigenvalues are roots of the following characteristic equation:
which means that
Thus, the eigenvalues have negative real parts if and , and the arguments of all the eigenvalues lie in the second or third quadrant, in this case as . Consequently, it can be stated that is locally asymptotically stable in both system (1) and system (2) when and hold. □
Similarly, we can analyze the stability of the system around and using the nature of eigenvalues of the corresponding Jacobian matrices.
Consider the case when the population from Patch-1 can disperse towards Patch-2 but cannot move back to Patch-1 again ().
If people are not allowed to disperse to Patch-1, then the disease-free equilibrium of systems (1) and (2) is locally asymptotically stable (LAS) for for Patch-1 and for Patch-2, i.e., .
The Jacobian matrix of systems (1) and (2) when people can not disperse to Patch-1 is provided as follows:
From the Jacobian matrix, we have the two eigenvalues and . The rest of the eigenvalues are , and , which are always negative. Thus, and when , , respectively. For all eigenvalues, here; thus, is locally asymptotically stable for both systems when and hold simultaneously, i.e., . □The endemic equilibrium of the proposed systems (1) and (2), if it exists, is locally asymptotically stable (LAS).
The Jacobian matrix at endemic equilibrium point is provided as follows:
The characteristic equation corresponding to is whereHere, the characteristic equation provides all of the eigenvalues with negative real parts for feasible and . In addition, the arguments of all the eigenvalues lie in the second or third quadrant in this case, i.e., . Thus, it can be stated that , whenever it exists, is locally asymptotically stable. □
Similarly, we can analyze the stability of the system around and using the nature of eigenvalues of the corresponding Jacobian matrices.
Consider the case when the population from Patch-2 can disperse towards Patch-1 but cannot move back to Patch-2 again ().
If people are not allowed to disperse to Patch-2, then the disease-free equilibrium of systems (1) and (2) is locally asymptotically stable (LAS) for for Patch-1 and for Patch-2, i.e., .
The Jacobian matrix of systems (1) and (2) when people can not disperse to Patch-2 is provided as follows:
Here, we have the two eigenvalues and . The rest of the eigenvalues are , and , which are always negative. The arguments of all the eigenvalues of the characteristic equation lie in the second or third quadrant, i.e., in this case. Thus, is locally asymptotically stable when and hold simultaneously for both system (1) and system (2), i.e., . □
The endemic equilibrium of the proposed systems (1) and (2), if it exists, is locally asymptotically stable (LAS).
The Jacobian matrix at endemic equilibrium point is provided as follows:
The characteristic equation corresponding to is where Here, the characteristic equation provides all of the eigenvalues with negative real parts for feasible and . Thus, all of the eigenvalues have negative real parts and the arguments of all the eigenvalues lie in the second or third quadrant. Therefore, , and it can be stated that , whenever it exists, is locally asymptotically stable. □Similarly, we can analyze the stability of the system around and using the nature of eigenvalues of the corresponding Jacobian matrices.
Consider the case when the population is allowed to disperse between the patches. Here, susceptible and recovered people can disperse among Patch-1 and Patch-2 ().
In the presence of population dispersal, the disease-free equilibrium of systems (1) and (2) is locally asymptotically stable (LAS) for and , i.e., .
The Jacobian matrix of systems (1) and (2) in the presence of population dispersal is provided as follows:
From the Jacobian matrix we have the two eigenvalues and . Thus, and hold when and , respectively. The other eigenvalues are roots of the equation
Here, all the roots will only have negative real parts, which leads to the conclusion that the arguments of all the eigenvalues should be greater than . Consequently, is locally asymptotically stable when and , i.e., . □
The endemic equilibrium of the proposed system, if it exists, is locally asymptotically stable.
The Jacobian matrix at endemic equilibrium point is provided as follows:
The characteristic equation corresponding to is whereHere, the characteristic equation provides all of the eigenvalues with negative real parts for feasible and . Thus, it can be stated that , whenever it exists, is locally asymptotically stable if and , and the arguments of all the eigenvalues lie in the second or third quadrant (i.e., ). □
The equilibrium point of the proposed system is locally asymptotically stable for .
The Jacobian matrix at endemic equilibrium point is provided as follows:
From the Jacobian matrix, we have one eigenvalue, ; thus, holds when . The other eigenvalues are roots of the characteristic equation corresponding to , which is whereHere, the characteristic equation provides all of the eigenvalues with negative real parts for a feasible equilibrium point. Thus, it can be stated that is locally asymptotically stable if , and the arguments of all the eigenvalues lie in the second or third quadrant (i.e., ). □
The equilibrium point of the proposed systems (1) and (2) is locally asymptotically stable for .
The Jacobian matrix at endemic equilibrium point is provided as follows:
From the Jacobian matrix, we have one eigenvalue, ; thus, holds when . The other eigenvalues are the roots of the characteristic equation corresponding to , which is where
Here, the characteristic equation provides all of the eigenvalues with negative real parts for a feasible equilibrium point. Thus, it can be stated that is locally asymptotically stable if , and the arguments of all the eigenvalues lie in the second or third quadrant (i.e., ). □
7. Numerical Simulation of System (1) without Implementation of Control Strategy
In this section, we validate the analytical results through numerical simulation. The scenarios presented here can help to visualize the dynamical behavior of the system in both the presence and absence of population dispersal. The parametric values used in this section are listed below in Table 1.
First, we discuss the situation when the population does not disperse between patches; Figure 2 depicts this scenario. Choosing , it is observed that the system converges to an infection-free state for both the patches, we obtain a stable , and the basic reproduction number becomes (see Figure 2a). Now, if starts to increase, then becomes a saddle point and we can observe that the infection invades Patch-1. In Figure 2b, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-2 and infected Patch-1 . On the other hand, if is increased instead of , then we can observe that the infection invades Patch-2. In Figure 2c, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-1 and infected Patch-2 . Now, if both and increase, then we have an infected system even if people do not disperse between the two patches. Figure 2d depicts the scenario at , where both Patch-1 and Patch-2 are infected with the disease and we obtain a stable with .
Next, we discuss the case when the population from Patch-1 can disperse towards Patch-2 but cannot move back to Patch-1 again (); Figure 3 depicts this scenario. Choosing , it is observed that the system converges to an infection-free state for both patches, we obtain a stable , and the basic reproduction number becomes (see Figure 3a). From this situation, if starts to increase then becomes a saddle point and we can observe that the infection invades in Patch-1. In Figure 3b, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-2 and infected Patch-1 . On the other hand, if is increased instead of , then we can observe that the infection invades in Patch-2, while Patch-1 remains infection-free. In Figure 3c, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-1 and infected Patch-2 . If both and increase, then we have an infected system even if people do not disperse between the two patches. Figure 3d depicts the scenario at , where both Patch-1 and Patch-2 are infected with the disease and we obtain a stable with .
Next, we discuss the case when the population from Patch-2 can disperse towards Patch-1 but cannot move back to Patch-2 again (); Figure 4 depicts this scenario. Choosing , it is observed that the system converges to an infection-free state for both the patches, we obtain a stable , and the basic reproduction number becomes (see Figure 4a). Now, increasing the value of leads to instability as becomes saddle, and we can observe that the infection invades in Patch-1. In Figure 4b, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-2 and infected Patch-1 . On the other hand, if is increased instead of , then we can observe that the infection invades Patch-2 while Patch-1 remains infection-free. In Figure 4c, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-1 and infected Patch-2 . Now, if both and increase, then we have an infected system even if people do not disperse between the two patches. Figure 4d depicts the scenario at , where both Patch-1 and Patch-2 are infected with the disease and we obtain a stable with .
Figure 5, shows the scenario when the population disperses between Patch-1 and Patch-2 (). Choosing , it is observed that the system converges to an infection-free state for both patches, we obtain a stable , and the basic reproduction number becomes (see Figure 5a). Now, this stable behavior of switches with increasing value of as the equilibrium becomes saddle, and we obtain a scenario where the infection invades in Patch-1. In Figure 5b, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-2 and infected Patch-1 . On the other hand, if is increased instead of , then we can observe that the infection invades Patch-2 while Patch-1 remains infection-free. In Figure 5c, we obtain the basic reproduction number as for and obtain a stable steady state containing an infection-free Patch-1 and infected Patch-2 . Now, if both and increase, then we have an infected system even if people do not disperse between the two patches. Figure 5d depicts the scenario at , where both Patch-1 and Patch-2 are infected with the disease and we obtain a stable with . Moreover, from Figure 2d and Figure 5d it can be observed that although the infected population count has increased slightly in Patch-1 in the presence of population dispersal, it has significantly decreased in Patch-2. Moreover, the highest number of recovered individuals is obtained in Patch-1 when people can disperse between the patches.
In Figure 6, we have analyzed the sensitivity of some of the system parameters in disease transmission, and have calculated corresponding sensitivity indices for the parametric values in Table 1 with . The fact that the spread of the disease has a greater ascendancy for than can be observed in Figure 6a. Moreover, an increase in the disease transmission rates actually proliferates the illness in a system, as the chances of spreading the infection become higher. On the other hand, for the chosen parametric values there is an inclination in the basic reproduction number for Patch-1 and a declination for Patch-2 with the increase in the population dispersal speed . In addition, if more people move towards Patch-1 , then the prevalence of the illness in that patch is expanded, while in the other patch, the spread is somewhat more controlled. Furthermore, the inversely proportional relationship of with indicates that the epidemic situation in Patch-1 can be controlled by increasing people’s movement towards Patch-2. However, the infectivity in Patch-2 will increase in this case. Figure 6b contains the sensitivity indices for the parametric values in Table 1, which are calculated as follows: and , , . On the other hand, Figure 6c indicates that the spread of disease has an ascendency for increasing .
In Figure 7, we have plotted the count of infected people in Patch-1 and Patch-2 for different movement scenarios. We mainly considered three scenarios: (i) people move to Patch-2 in larger numbers ; (ii) people move to Patch-1 with larger probabilities ; and (iii) the probabilities are same for both patches . Figure 7a shows that the count of infected people in Patch-1 is lowest when both susceptible and recovered people move to Patch-2, and that this situation leads to an inclination in the infectivity there (see Figure 7b). In fact, higher population dispersal probabilities (either of susceptible, recovered, or both) towards Patch-1 suppress the infection level in Patch-2 even for smaller dispersal speeds.
8. Numerical Simulation of System (2) without Implementation of Control Strategy
The numerical simulation of the system (2) is described in this section. We have used the FDE12 MatLab function, which is based on the Adams–Bashforth–Moulton algorithm introduced by Roberto Garrappa ([32]), and used a predictor-corrector method.
Consider the situation when there is no population dispersal between the patches (m = 0)
To analyze the dynamical behavior of all people, the values of the parameters in Table 1 are employed. The initial populations of both patches are taken in the following manner: . Here, we are more concerned about the disease-free equilibrium. The population of the two patches has been considered as 100,000 people per patch, and we have taken t = 1 day as the time unit and . At first, we have considered . From the direct calculation, we have ; therefore, the disease-free equilibrium is stable (see Figure 8).
Consider the situation when there is population dispersal between the patches ()
In this case, we fix and consider ; then, we have . The time series in Figure 9 shows the stability of the disease-free equilibrium.
9. Optimal Control Problem
The deterministic system (1) is further proposed with an effective control strategy to reduce the infection load. People’s movement between the two patches plays an important role in disease transmission. During an outbreak, people try to move back to their homes as soon as possible, and this behavior induces the transmission rate only. Furthermore, if people travel from a disease-prone area, it increases the chance of spreading the virus only. Imposing strict restrictions sometimes becomes necessary to control the continuous upsurge in transmission; thus, the dispersal speed of the population is considered as the control policy in this work, as it depends on the eruption of the infection and changes with time. Here, we study the characteristics of the applied time-dependent control policy mathematically. As the public authorities usually take immediate actions and impose strict restrictions only when the numbers of reported cases start to exceed the safe limit, instead of constant values, time-dependent dispersal rate function is considered here with . By implementing this control policy, the overall chances of becoming infected with the virus strain are lessened. Here, 1 indicates that people’s dispersal works with its fullest intensity, while 0 denotes full disconnection between the patches.
Our main focus is to determine the optimal control intervention with minimum implemented cost; hence, the region for the control intervention is provided as
where is the final time up to which the control policies are executed and is a measurable and bounded function.Deduction of Total Cost Which Needs to Be Minimized Due to Dispersal
The total cost incurred to manage people’s emigration (or immigration) from (or to) the patches is provided by . The term denotes the cost incurred in lost manpower because of the infected population in the patches. This is known as the opportunity loss, which includes productivity loss due to the infection. The term represents the cost associated with managing the dispersal of population among the patches, which includes spreading information about the disease via different programs, informing the population about the importance of maintaining the proposed rules and regulations, etc. As this includes the cost of associated efforts towards convincing people, it is relatively higher. In this work, the nonlinearity term for the mitigation strategy is chosen up to the second order [33]. This work analyzes how the optimal control strategy representing people’s movement reduces the overall count of infected people in the system.
The following control problem is considered based on previous discussions, along with the cost functional J to be minimized [34]:
(5)
subject to the model system(6)
with initial conditions and We have already considered for and for . The functional J denotes the total incurred cost as stated, and the integrand denotes the cost at time t. The positive parameters are weight constants balancing the units of the integrand [33,35]. The optimal control interventions exist in , and mainly minimize the cost functional J.The optimal control intervention in Ψ of the control systems (5) and (6) exists such that .
The proof is presented in Appendix A. □
For the optimal control and corresponding optimal states of system (6), there exist adjoint variables satisfying the canonical equations:
(7)
with transversality conditions for . The corresponding optimal control which minimizes in Ψ is provided as follows:(8)
The proof is presented in Appendix A. □
Now, the model system in the Caputo fractional system with the implemented control intervention is as follows:
(9)
with initial conditions , andFor the optimal control and corresponding optimal states of system (9), there exist adjoint variables satisfying the canonical equations:
(10)
where is the α-ordered Riemann–Liouville derivative with initial point 0 and transversality conditions for . The corresponding optimal control which minimizes in Ψ is provided as follows:(11)
The proof is similar to the proof of Theorem 14 mentioned in the Appendix A.2 of Appendix A. □
10. Numerical Simulation of the Optimal Control Problem (6)
In the proposed system (6), we have considered the population dispersal speed as the control intervention in order to analyze its impact on reducing the infection level. The control intervention is assumed to be time-independent, as it changes according to the disease prevalence, fatality, and governmental rules and restrictions. Through employing a forward–backward sweep approach, numerical simulations have been carried out to demonstrate the effects of different control strategies on the behavior of the system [36,37]. The parametric values are chosen from Table 1 along with , and the positive weight constants are taken as , and . Additionally, it is presumed that the control techniques will be used consistently for two months, or days. The initial population size is chosen as follows: , , , , , .
The dynamics of model (6) when the dispersal speed is chosen as a constant value instead of a time-dependent control measure are shown in Figure 10. In this scenario, the population size is , , , , , at .
If the control strategy is properly implemented through people’s dispersal (or migration) strategy, the burden of the infection may be managed more effectively. Hence, we have emphasized the fact that susceptible and recovered people move or disperse between Patch-1 and Patch-2 according to disease severity, allowing the increased disease transmission to be curbed. The population trajectories in the presence of the control strategy are shown in Figure 11, and at , the population is represented as , , , , , . In this situation, the count of the infected population increases in Patch-1 and declines in Patch-2. Surprisingly, we find the highest number of recovery cases in both patches. Higher recovery results in a slightly lower susceptible population count in the patches than in the scenario where dispersal speed is not used as a control approach. Figure 12b reveals the optimal graphs of the control intervention, where can be seen that its intensity remains at the highest level for almost 10 days, then starts to diminish in the last few days.
The significance of an implemented control measure is determined by how cost-effective it is. Figure 12 and Figure 13 illustrate the impact of the time-dependent control measure on the cost design analysis (Figure 12a) and the infected population count of both patches (Figure 13). Without the control measure, the cost is due to the impact of the disease load on the production loss. When the control measure is not implemented, the infected population in Patch-1 is high enough, as is that in Patch-2; in fact, it is at its maximum count, which results in higher opportunity loss and leads to a rise in the economic burden. It is shown in Figure 13a that although there is a small inclination in in the presence of optimal dispersal, we have the fewest number of infected people in Patch-2 in the presence of this control policy, and this situation makes the optimal cost much lower. This finding indicates that implementing the control strategy results in a lower cost profile than not implementing the intervention. Thus, it may be said that the control policy is financially viable.
10.1. Effect of the Weight Constant on the Optimal Control Policy
When determining the weight that needs to be applied to a particular control, weight constants are essential. Here, we vary the weight constant in order to identify its effectiveness. Figure 14c illustrates that placing higher weight on migration results in an upsurge in overall cost. Moreover, increasing decreases the infection level in Patch-1 but spreads the infectivity in Patch-2 (see Figure 14a,b). It can be seen in Figure 14d that if the related weight ) is larger, then more effort on the part of the control is needed to lower the cost and total number of infected people.
10.2. Effect of Behavioral Changes and People’s Awareness on Disease Transmission When
In this work, we have mainly focused on the dynamical behavior of an SIRS model in which the population disperses between two patches. In Figure 15, we have tried to show how the infection level is affected when people’s behavioral changes, mainly the population’s awareness towards preventive measures or precautionary measures, are taken into consideration in the presence of optimal dispersal. It can be observed that the count of the infected population in both patches is reduced when people start to pay attention to restrictions and precautions. Figure 15c shows that if people adopt proper behavioral changes, then a comparatively lower effort in terms of the control policy is required to reduce the count of infected individuals in the system.
11. Numerical Simulation of the Optimal Control Problem (9)
In this section, we discuss the effect of optimal control on the disease transmission dynamics. We can set the dispersal rate (m) under the control such that the cost function is minimized. Here, we have assumed that , and . All other parametric values are taken from Table 1. We performed the simulation for . We have assumed t = 1 day as a time unit and .
We have developed a fractional order optimal control problem using an iterative approach (Euler’s forward and backward) in the MatLab interface [32,38]. The procedure is summarized below. A two-point boundary value problem with a set of fractional-order differential equations defines the optimality of the system (9). The initial value system (9) is an initial value problem and the adjoint system is a boundary value problem. The state system is solved using the forward iteration approach, while the co-state system is solved using the backward iteration method using the Matlab code presented below.
The state system (9) is solved using the following iterative scheme:
where , , and is the length of the time step; moreover, is the value of at iteration. The last term of each of the above system of equations stands for memory. The adjoint system (10) is solved using the backward iteration method with terminal conditions with the following iterative scheme: where for . The system below updates the optimal control:We have developed MatLab code using the above algorithm and chose throughout the numerical simulations. The number of infected people in both patches can be gradually reduced by controlling . The optimal value of the migration rate quickly approaches 1 as the order of the derivative increases. The effect of the other parameters has already been mentioned in Section 10. In this section, we have only observed the impact of the order of derivatives in the control system. Figure 16 depicts the time series of the S, I, and R variables for . Figure 17 provides a comparison of the time series of the system variables along with the control variable and cost function. It is apparent that the optimal stage can be archived quickly for a lower value of .
12. Conclusions
When it comes to population dispersal, infectious diseases tend to be a cause for concern. In densely populated areas, respiratory diseases such as flu or COVID-19 have a tendency to be highly infectious, as they can spread very quickly. There are several other aspects that causes high transmission. It includes mode of transmission, period of infectivity, efficacy of implemented preventative measures, etc., which can be taken into account. In order to control the spread of disease during population movement or traveling, public health is essential. The migration of people, which has increased with time and growing socioeconomic conditions, expanding tourism, etc., acts as a contributing factor. This movement alleviates the spread of an infectious disease; thus, analyzing an epidemic model in the presence of population dispersal is of utmost importance.
In this work, a compartmental SIRS model has been proposed in a two-patch environment, where it is assumed that only susceptible and recovered populations are allowed to disperse between the patches. A deterministic model and a fractional-order system are analyzed here in order to observe the detailed dynamical scenario. Four situations are taken into consideration while analyzing the equilibrium points of the system: (i) no dispersal between the patches; (ii) people can travel from Patch-1 to Patch-2 but cannot move back to Patch-1; (iii) people can travel from Patch-2 to Patch-1 but cannot move back to Patch-2; and (iv) people can freely disperse between the two patches. In these cases, we respectively reach a disease-free equilibrium state, an infection-free Patch-2 and infected Patch-1 state, an infection-free Patch-1 and infected Patch-2 state, and an infected state in which both patches are infected. The numerical scenario shows that the disease transmission rates in both patches play significant roles in the stability of the system, as adjusting these parameters decides whether the disease will invade the patch or not. It is notable that although the infected population count increases slightly in Patch-1 in the presence of population dispersal, it significantly decreases in Patch-2. For the chosen parametric values, the results depict that there is an inclination in the basic reproduction number for Patch-1 and a declination for Patch-2 with the increase in the population dispersal speed . In the later part of the paper, a corresponding optimal control problem is formulated in which the dispersal speed of the population is considered to be time-dependent instead of a constant value. The numerical scenario reveals that the count of the infected population increases in Patch-1 and declines in Patch-2 in the presence of optimal dispersal, and that the highest number of recovery cases has been obtained in both patches.
In the model simulation, fractional order has a noticeable impact. It is noteworthy that when the derivative order decreases, the recovery rate increases. Both patches have a comparable rise in the infected population, although at a lower derivative rate. In the control situation, the order of derivatives is less, and the optimal state is reached faster. These effects can be quite useful when trying to fit actual data into real-life scenarios.
While the model proposed here has shown some interesting dynamical behavior, there are many remaining items that can be implemented later as an extension of this work. In this work, we have only dealt with the deterministic and fractional approaches; environmental stochasticity has not been incorporated here. It would be interesting to explore whether Gaussian noise plays a significant role in this dispersed epidemic model. Secondly, we have only considered scenarios in which susceptible and recovered people are allowed to move between two patches; any asymptomatic or pre-symptomatic stage has not been considered here. It is quite evident that people belonging to these stages shall also move between the patches, and so the future researchers can possibly observe the impact of this movement on the dynamics of the system. Moreover, a certain amount of time is required for a person to be infected with a certain type of virus. As it is not an instantaneous process, consideration of a delay parameter can make the model more realistic. The present work can be studied further by incorporating the facts mentioned above in order to observe various realistic dynamical behaviors.
All the authors have participate equally in all the aspects of this paper: conceptualization, methodology, investigation, formal analysis, writing—original draft preparation, writing—review and editing. All authors have read and agreed to the published version of the manuscript.
The data used to support the findings of the study are available within the article.
The authors are grateful to the anonymous referees, for their careful reading, valuable comments, and helpful suggestions, which have helped to significantly improve the presentation of this work.
The authors state that they do not have any conflicts of interest in respect of this study.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Figure 2. Trajectory profiles of the populations in Patch-1 and Patch-2 around (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.], and (d) [Forumla omitted. See PDF.] in the absence of dispersal.
Figure 3. Trajectory profiles of populations in Patch-1 and Patch-2 around (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.], and (d) [Forumla omitted. See PDF.] when people disperse towards Patch-2 from Patch-1 but do not move back to Patch-1.
Figure 4. Trajectory profiles of populations in Patch-1 and Patch-2 around (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.], and (d) [Forumla omitted. See PDF.] when people disperse towards Patch-1 from Patch-2 but do not move back to Patch-2.
Figure 5. Trajectory profiles of populations in Patch-1 and Patch-2 around (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.], and (d) [Forumla omitted. See PDF.] when people disperse between Patch-1 and Patch-2.
Figure 6. (a) Different plots of the basic reproduction number [Forumla omitted. See PDF.] with the variation of [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.], (b) sensitivity index of the parameters [Forumla omitted. See PDF.], and [Forumla omitted. See PDF.] for [Forumla omitted. See PDF.], and (c) plots of the basic reproduction number [Forumla omitted. See PDF.] with the variation of [Forumla omitted. See PDF.].
Figure 7. Count of infected population in (a) Patch-1 [Forumla omitted. See PDF.] and (b) Patch-2 [Forumla omitted. See PDF.] for different dispersal probabilities.
Figure 8. Time series of different populations of model (2) for different values of [Forumla omitted. See PDF.] when [Forumla omitted. See PDF.].
Figure 9. Trajectory profiles of (a) susceptible, (b) infected, and (c) recovered populations around [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], and [Forumla omitted. See PDF.] when people disperse between Patch-1 and Patch-2 for different values of [Forumla omitted. See PDF.] when [Forumla omitted. See PDF.].
Figure 10. Profiles of populations in the absence of time-dependent control intervention.
Figure 11. Population profiles in the presence of time-dependent control intervention.
Figure 12. (a) Cost distribution with time-dependent control policy and (b) profile of optimal control [Forumla omitted. See PDF.].
Figure 13. Profiles of infected population in (a) Patch-1 and (b) Patch-2 under time-dependent control policy.
Figure 14. Impact of the weight constant [Forumla omitted. See PDF.] on (a) the infected population in Patch-1 [Forumla omitted. See PDF.], (b) the infected population in Patch-2 [Forumla omitted. See PDF.], (c) the optimal cost [Forumla omitted. See PDF.], and (d) the profile of the optimal control [Forumla omitted. See PDF.].
Figure 15. Impact of behavioral changes on disease transmission. Profiles of infected population on (a) Patch-1 [Forumla omitted. See PDF.], (b) Patch-2 [Forumla omitted. See PDF.], and (c) the optimal control [Forumla omitted. See PDF.] are shown. For purple ([Image omitted. Please see PDF.]) colored curves, disease transmission rates are [Forumla omitted. See PDF.], while for green ([Image omitted. Please see PDF.]) colored curves disease transmission rates are [Forumla omitted. See PDF.] for [Forumla omitted. See PDF.].
Figure 16. Time series of (a) susceptible, (b) infected and (c) recovered populations profile of model (9) without control, with yellow ([Image omitted. Please see PDF.]) colored curves for [Forumla omitted. See PDF.] and red ([Image omitted. Please see PDF.]) colored curves for [Forumla omitted. See PDF.].
Figure 17. Comparison of time series of different population profiles (a) susceptible, (b) infected and recovered , (c) control function and cost function of model (9) with optimal control for [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.].
Parameter values used for numerical simulations of both inter-ordered and fractional-ordered systems.
Parametric Values | |||||||
---|---|---|---|---|---|---|---|
| | | | | | | |
10 | | | | | | | 12 |
| | | | m | | | - |
| | | | | | | - |
Appendix A
Appendix A.1. Existence of Optimal Control Functions
Here, the existence of optimal control interventions which minimize the cost functional within a finite time interval are provided.
In model (
Here, system (
Here,
Appendix A.2. Characterization of Control Interventions
Pontryagin’s Principle helps to find the control interventions of an optimal system [
For system (
From the optimality conditions
Appendix A.3. Optimal System
The optimal system is stated below along with the optimal control policy
Appendix B
Appendix B.1. Some Results on Fractional Order Systems
([
([
([
([
References
1. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. A; 1927; 115, pp. 700-721.
2. Wang, W. Epidemic models with population dispersal. Mathematics for Life Science and Medicine; Springer: Berlin/Heidelberg, Germany, 2007; pp. 67-95.
3. Bichara, D.; Kang, Y.; Castillo-Chavez, C.; Horan, R.; Perrings, C. SIS and SIR epidemic models under virtual dispersal. Bull. Math. Biol.; 2015; 77, pp. 2004-2034. [DOI: https://dx.doi.org/10.1007/s11538-015-0113-5]
4. Jana, S.; Haldar, P.; Kar, T. Optimal control and stability analysis of an epidemic model with population dispersal. Chaos Solitons Fractals; 2016; 83, pp. 67-81. [DOI: https://dx.doi.org/10.1016/j.chaos.2015.11.018]
5. Hu, L.; Wang, S.; Zheng, T.; Hu, H.; Kang, Y.; Nie, L.F.; Teng, Z. The Effects of Migration and Limited Medical Resources of the Transmission of SARS-CoV-2 Model with Two Patches. Bull. Math. Biol.; 2022; 84, 55. [DOI: https://dx.doi.org/10.1007/s11538-022-01010-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35377056]
6. Jin, Y.; Wang, W. The effect of population dispersal on the spread of a disease. J. Math. Anal. Appl.; 2005; 308, pp. 343-364. [DOI: https://dx.doi.org/10.1016/j.jmaa.2005.01.034]
7. Abhishek, V.; Srivastava, V. SIS Epidemic Spreading under Multi-layer Population Dispersal in Patchy Environments. IEEE Trans. Control. Netw. Syst.; 2023; pp. 1-12. [DOI: https://dx.doi.org/10.1109/TCNS.2023.3272839]
8. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999.
9. Das, M.; Samanta, G.P. Optimal Control of Fractional Order COVID-19 Epidemic Spreading in Japan and India 2020. Biophys. Rev. Lett.; 2020; 15, pp. 207-236. [DOI: https://dx.doi.org/10.1142/S179304802050006X]
10. Guo, T. The Necessary Conditions of Fractional Optimal Control in the Sense of Caputo. J. Optim. Theory Appl.; 2013; 156, pp. 115-126. [DOI: https://dx.doi.org/10.1007/s10957-012-0233-0]
11. Agarwal, O. A general formulation and solution scheme for fractional optimal control problems. Nonlinear Dyn.; 2004; 38, pp. 323-337. [DOI: https://dx.doi.org/10.1007/s11071-004-3764-6]
12. Das, M.; Samanta, G.; De la Sen, M. A Fractional Ordered COVID-19 Model Incorporating Comorbidity and Vaccination. Mathematics; 2021; 9, 2806. [DOI: https://dx.doi.org/10.3390/math9212806]
13. Das, M.; Samanta, G.; De la Sen, M. A Fractional Order Model to Study the Effectiveness of Government Measures and Public Behaviours in COVID-19 Pandemic. Mathematics; 2022; 10, 3020. [DOI: https://dx.doi.org/10.3390/math10163020]
14. Prosper, O.; Ruktanonchai, N.; Martcheva, M. Assessing the role of spatial heterogeneity and human movement in malaria dynamics and control. J. Theor. Biol.; 2012; 303, pp. 1-14. [DOI: https://dx.doi.org/10.1016/j.jtbi.2012.02.010] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22525434]
15. Yan, A.W.; Black, A.J.; McCaw, J.M.; Rebuli, N.; Ross, J.V.; Swan, A.J.; Hickson, R.I. The distribution of the time taken for an epidemic to spread between two communities. Math. Biosci.; 2018; 303, pp. 139-147. [DOI: https://dx.doi.org/10.1016/j.mbs.2018.07.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30089576]
16. Saha, S.; Samanta, G.P.; Nieto, J.J. Epidemic model of COVID-19 outbreak by inducing behavioural response in population. Nonlinear Dyn.; 2020; 102, pp. 455-487. [DOI: https://dx.doi.org/10.1007/s11071-020-05896-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32863581]
17. Saha, S.; Samanta, G.; Nieto, J.J. Impact of optimal vaccination and social distancing on COVID-19 pandemic. Math. Comput. Simul.; 2022; 200, pp. 285-314. [DOI: https://dx.doi.org/10.1016/j.matcom.2022.04.025] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35531464]
18. Saha, S.; Samanta, G.P. Analysis of a COVID-19 Model Implementing Social Distancingas an Optimal Control Strategy. Integrated Science of Global Epidemics. Integrated Science; Rezaei, N. Springer: Berlin/Heidelberg, Germany, 2023; Volume 14, pp. 211-258. [DOI: https://dx.doi.org/10.1007/978-3-031-17778-1_10]
19. Takeuchi, Y.; Liu, X.; Cui, J. Global dynamics of SIS models with transport-related infection. J. Math. Anal. Appl.; 2007; 329, pp. 1460-1471. [DOI: https://dx.doi.org/10.1016/j.jmaa.2006.07.057]
20. Wang, L.; Yang, W. Global dynamics of a two-patch SIS model with infection during transport. Appl. Math. Comput.; 2011; 217, pp. 8458-8467. [DOI: https://dx.doi.org/10.1016/j.amc.2011.03.045]
21. Wang, W.; Mulone, G. Threshold of disease transmission in a patch environment. J. Math. Anal. Appl.; 2003; 285, pp. 321-335. [DOI: https://dx.doi.org/10.1016/S0022-247X(03)00428-1]
22. Wang, W.; Zhao, X.Q. An epidemic model in a patchy environment. Math. Biosci.; 2004; 190, pp. 97-112. [DOI: https://dx.doi.org/10.1016/j.mbs.2002.11.001]
23. Wang, W.; Zhao, X.Q. An age-structured epidemic model in a patchy environment. SIAM J. Appl. Math.; 2005; 65, pp. 1597-1614. [DOI: https://dx.doi.org/10.1137/S0036139903431245]
24. Nandi, S.K.; Jana, S.; Mandal, M.; Kar, T. Mathematical analysis of an epidemic system in the presence of optimal control and population dispersal. Biophys. Rev. Lett.; 2018; 13, pp. 1-17. [DOI: https://dx.doi.org/10.1142/S1793048018500017]
25. Wan, H. An SEIS epidemic model with transport-related infection. J. Theor. Biol.; 2007; 247, pp. 507-524. [DOI: https://dx.doi.org/10.1016/j.jtbi.2007.03.032]
26. Arino, J.; van den Driessche, P. A multi-city epidemic model. Math. Popul. Stud.; 2003; 10, pp. 175-193. [DOI: https://dx.doi.org/10.1080/08898480306720]
27. Liu, X.; Chen, X.; Takeuchi, Y. Dynamics of an SIQS epidemic model with transport-related infection and exit–entry screenings. J. Theor. Biol.; 2011; 285, pp. 25-35. [DOI: https://dx.doi.org/10.1016/j.jtbi.2011.06.025]
28. Cheng, C.Y. Adaptive dispersal effect on the spread of a disease in a patchy environment. Appl. Math. Model.; 2017; 47, pp. 17-30. [DOI: https://dx.doi.org/10.1016/j.apm.2017.03.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32287941]
29. Chen, Y.; Yan, M.; Xiang, Z. Transmission dynamics of a two-city SIR epidemic model with transport-related infections. J. Appl. Math.; 2014; 2014, 764278. [DOI: https://dx.doi.org/10.1155/2014/764278]
30. van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci.; 2002; 180, pp. 29-48. [DOI: https://dx.doi.org/10.1016/S0025-5564(02)00108-6]
31. Arriola, L.; Hyman, J. Lecture Notes, Forward and Adjoint Sensitivity Analysis: With Applications in Dynamical Systems, Linear Algebra and Optimisation. Math. Theor. Biol. Inst.; 2005.
32. Garrappa, R. On linear stability of predictor-corrector algorithms for fractional differential equations. Int. J. Comput. Math.; 2010; 87, pp. 2281-2290. [DOI: https://dx.doi.org/10.1080/00207160802624331]
33. Kassa, S.M.; Ouhinou, A. The impact of self-protective measures in the optimal interventions for controlling infectious diseases of human population. J. Math. Biol.; 2015; 70, pp. 213-236. [DOI: https://dx.doi.org/10.1007/s00285-014-0761-3]
34. Castilho, C. Optimal control of an epidemic through educational campaigns. Electron. J. Differ. Equations (EJDE); 2006; 2006, pp. 1-11.
35. Gaff, H.; Schaefer, E. Optimal control applied to vaccination and treatment strategies for various epidemiological models. Math. Biosci. Eng.; 2009; 6, pp. 469-492. [DOI: https://dx.doi.org/10.3934/mbe.2009.6.469]
36. Yusuf, T.T.; Benyah, F. Optimal control of vaccination and treatment for an SIR epidemiological model. World J. Model. Simul.; 2012; 8, pp. 194-204.
37. Kirk, D.E. Optimal Control Theory: An Introduction; Dover Books on Electrical Engineering Dover Publications: New York, NY, USA, 2012.
38. Kheiri, H.; Jafari, M. Optimal control of a fractional-order model for the HIV/AIDS epidemic. Int. J. Biomath.; 2018; 11, 1850086. [DOI: https://dx.doi.org/10.1142/S1793524518500869]
39. Coddington, A.; Levinson, N. Theory of Ordinary Differential Equations; International series in pure and applied mathematics; Tata McGraw-Hill Companies: New York, NY, USA, 1955.
40. Fleming, W.; Rishel, R. Deterministic and Stochastic Optimal Control; Springer: New York, NY, USA, 1975; Volume 1.
41. Pontryagin, L.S. Mathematical Theory of Optimal Processes; CRC Press: Boca Raton, FL, USA, 1987.
42. Liang, S.; Wu, R.; Chen, L. Laplace transform of fractional order differential equations. Electron. J. Differ. Equ.; 2015; 2015, pp. 1-15.
43. Delavari, H.; Baleanu, D.; Sadati, J. Stability analysis of Caputo fractional-order non linear system revisited. Nonlinear Dyn.; 2012; 67, pp. 2433-2439. [DOI: https://dx.doi.org/10.1007/s11071-011-0157-5]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Migration or dispersal of population plays an important role in disease transmission during an outbreak. In this work, we have proposed an SIRS compartmental epidemic model in order to analyze the system dynamics in a two-patch environment. Both the deterministic and fractional order systems have been considered in order to observe the impact of population dispersal. The following analysis has shown that we can have an infected system even if the basic reproduction number in one patch becomes less than unity. Moreover, higher dispersal towards a patch controls the infection level in the other patch to a greater extent. In the optimal control problem (both integer order and fractional), it is assumed that people’s dispersal rate will depend on the disease prevalence, and as such will be treated as a time-dependent control intervention. The numerical results reveal that there is a higher amount of recovery cases in both patches in the presence of optimal dispersal (both integer order and fractional). Not only that, implementation of people’s awareness reduces the infection level significantly even if people disperse at a comparatively higher rate. In a fractional system, it is observed that there will be a higher amount of recovery cases if the order of derivative is less than unity. The effect of fractional order is omnipotent in achieving a stable situation.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details



1 MS2 Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, ON N2L 3C5, Canada;
2 Department of Mathematics, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India;