International Journal of Mathematics and Mathematical Sciences
Volume 2017 (2017), Article ID 4804897, 19 pages
https://doi.org/10.1155/2017/4804897

Research Article

Optimal Control Techniques on a Mathematical Model for the Dynamics of Tungiasis in a Community

Jairos Kahuru,¹ Livingstone S. Luboobi,1,2 and Yaw Nkansah-Gyekye¹

¹School of Computational and Communication Science and Engineering, Nelson Mandela African Institution of Science and Technology, P.O. Box 447, Arusha, Tanzania
²Department of Mathematics, Makerere University, P.O. Box 7062, Kampala, Uganda

Correspondence should be addressed to Jairos Kahuru; shinzehk@nm-aist.ac.tz

Received 20 March 2017; Accepted 2 July 2017; Published 14 August 2017

Academic Editor: Nawab Hussain

Copyright © 2017 Jairos Kahuru et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Tungiasis is a permanent penetration of female sand flea “Tunga penetrans” into the epidermis of its host. It affects human beings and domestic and sylvatic animals. In this paper, we apply optimal control techniques to a Tungiasis controlled mathematical model to determine the optimal control strategy in order to minimize the number of infested humans, infested animals, and sand flea populations. In an attempt to reduce Tungiasis infestation in human population, the control strategies based on personal protection, personal treatment, educational campaign, environmental sanitation, and insecticidal treatments on the affected parts as well as on animal fur are considered. We prove the existence of optimal control problem, determine the necessary conditions for optimality, and then perform numerical simulations. The numerical results showed that the control strategy comprises all five control measures and that which involves the three control measures of insecticide control, insecticidal dusting on animal furs, and environmental hygiene has the significant impact on Tungiasis transmission. Therefore, fighting against Tungiasis infestation in endemic settings, multidimensional control process should be employed in order to achieve the maximum benefits.

1. Introduction

Tungiasis is a skin disease caused by the sand flea “Tunga penetrans”; the disease is endemic in some poor resource communities where various domestic and sylvatic animals act as reservoirs for this zoonosis [1]. The flea infestation is associated with poverty and occurs in many resource-poor communities in the Caribbean, South America, and sub-Saharan Africa [2]. Transmission of Tungiasis is strictly by infestation of humans and animal reservoirs by “Tunga penetrans” when they are in contact with sandy soil in which female fleas are present or when in contact with infested animal reservoirs as it is known that the animal reservoirs harbor the fleas [3]. Tungiasis results in significant morbidity, manifesting itself in a number of symptoms such as severe local inflammation, autoamputation of digits, deformation and loss of nails, formation of fissures and ulcers, gangrene, and walking difficulties [4]. Moreover it may result into secondary infection caused by transmission of blood-borne pathogens such as hepatitis B and C virus and possibly also HIV/AIDS when a single nonsterile instrument is used to remove the jiggers from different affected individuals [5].

Mathematical models have played a major role in increasing understanding of the underlying mechanisms which influence the spread of the diseases and provide guidelines as to how the spread can be controlled [6, 7]. Optimal control theory is a powerful mathematical tool which makes the decision involving complex dynamical systems. It is a standard method for solving dynamic optimization problems, when those problems are expressed in continuous time [8]. Optimal control theory was developed by the Russian mathematician Lev S. Pontryagin (1908–1988) and his coworkers with the formulation and proof of the Pontryagin Maximum Principle (Pontryagin et al., 1962). Optimal control is the process of determining control and state trajectories for a dynamic system over a period of time to minimize a performance index [9]. Optimal control problem is represented by a set of differential equations describing the paths of the control variables that minimize the cost functional and has been used successfully to make decisions involving biological or medical models [10]. The formulation of an optimal control problem requires a mathematical model of the system to be controlled, a specification of the performance index (cost function), and a specification of all boundary conditions on states and constraints to be satisfied by states and controls [11]. Pontryagin’s maximum (or minimum) principle of optimal control gives the fundamental necessary conditions for a controlled trajectory to be optimal [12]. The principle technique is to transform the constrained dynamic optimization problem into an unconstrained problem, by allowing each of the adjoint variables to correspond to each of the state variables accordingly and combining the results with the objective functional [13]. The resulting function is known as the Hamiltonian, which is used to solve a set of necessary conditions that an optimal control and corresponding state variables must satisfy. The necessary conditions are the optimality solutions and adjoint equations which form the optimality system. The optimality system consists of the state system and adjoint system with initial and transversal conditions together with characterization of optimal control.

To the best of our knowledge Tungiasis dynamical model with application of optimal control technique has not been done. Therefore we are going to refer to other infectious diseases with similar characteristics where the optimal control theory has been applied. Bonyah et al. [20] applied optimal control theory to a Buruli ulcer model that takes into account human, water bug, and fish populations as well as Microbacterium ulcerans in the environment. The control measures were applied on mass treatment, insecticide, and mass education to minimize the number of infected hosts, vectors, and infected fishes. The optimality system was determined and computed numerically for several scenarios. The results showed that the combination of all the control measures, mass treatment, insecticide, and mass education, is capable of helping reduce the number of infected humans, water bugs, small fishes, and Mycobacterium ulcerans in the environment. Isere et al. [21] developed the optimal control model that includes two time dependent control functions with one minimizing the contact between the susceptible and the bacteria and the others, the population of bacteria in water. The results from the numerical solutions showed that increasing the susceptible pool and the infected populations above some threshold values were responsible for reducing cholera epidemic and the difference between the growth rate and the loss rate of the bacteria played a huge role in the outbreak of the disease. Devipriya and Kalaivani [22] conducted the study on “Optimal Control of Multiple Transmission of Water-Borne Disease.” A controlled SIWR model was considered. The control measures represented an immune boosting and pathogen suppressing drugs. Their objective function was based on a combination of minimizing the number of infected individuals and the cost of the drugs dose. The numerical results have shown that both the vaccines resulted in minimizing the number of infected individuals and at the same time in a reduction of the budget related to the disease.

In this paper, the Tungiasis dynamical model with control measures is presented and a detailed qualitative optimal control model that minimizes the number of infested individuals (humans and animals) and sand fleas with minimal cost of implementing the control measures is developed. We establish the proof for existence of the optimal control and analyze the optimal control problem in order to determine the necessary conditions for optimality using the Pontryagin’s maximum principle (Pontryagins et al., 1962). We then determine numerically the optimality system for several scenarios. Our paper is arranged as follows. In Section 2, we formulate an optimal control model. In Section 3, we analyze the optimal control model by determining the conditions for existence of optimal control and the necessary conditions for optimality. In Section 4, we carry out numerical simulations and discussion of the results and Section 5 is the conclusion.

2. Formulation of Optimal Control Problem

We formulate an optimal control model for Tungiasis disease in order to derive five optimal control measures with minimal implementation cost to eradicate the disease after a defined period of time. We employ the control efforts in human, animal reservoirs, and adult flea populations and is the failure rate for the control efforts for . We let be the effort of controlling the flea infested soil environment with insecticides spraying, be the efforts of controlling the flea infested animal reservoir through dusting them with ant-flea compounds, be the efforts of controlling the transmission from flea infested environment to susceptible animals (this can be achieved by environmental hygiene and cementing the floors), be the efforts of controlling the transmission from flea infested animals to susceptible humans (this can be achieved by educating people not to live with animals in the same quarters or sharing common resting places), and be the efforts of controlling transmission from flea infested environment to susceptible humans (this can be achieved by environmental hygiene, cementing the floors, covering of feet with solid shoes, and application of plant based repellent (Zanzarin) within the time interval of ). Therefore we assume that the mortality rate of jigger fleas in the soil environment is increased by the factor , on-host spraying of infested animals will reduce the shedding rate of adult jigger fleas into the environment by a fraction , and the animal to animal effective contact rate is reduced at the same fraction because spraying insecticides on animal fur will reduce the transmission of infestation within animal population. The transmission rate from the soil environment to animal hosts is reduced by the factor , the factor reduces the transmission from severely infested animal reservoirs to susceptible humans, and the factor reduces the transmission from flea infested soil environment to susceptible humans. Here, we consider the model developed by Kahuru et al. [19] whereby we add a distinct epidemiological compartment which represents human beings under treatment and incorporate the five control measures , , , , and as defined above. In the submodel of human population, the total human population is subdivided into susceptible population mildly infested population , the severely infested population , and the human treatment class denoted by ; therefore we have . We assume that the humans are recruited into through birth by the adults at a rate . Individuals in class acquire infestation from the severely infested animal reservoirs and move to class at a rate and may also acquire infestation from the flea infested soil environment and move to class at a rate . may as well acquire infestation from the flea infested soil environment and progresses to class at a rate . Classes and seek treatment at the respective rates and and join class, and eventually the treated individuals revert back to join at a progression rate . Individuals in compartments , , and suffer a natural mortality rate and for the compartment they suffer a natural mortality at a rate and the disease induced mortality at a rate . In the submodel of animal reservoir population, the total animal reservoir population is subdivided into susceptible population mildly infested population and the severely infested population ; therefore we have . We assume that the animals are recruited into through birth by the adults at a rate . Individuals in class acquire infestation from the severely infested animal reservoirs and move to at a rate and also may acquire infestation from the flea infested environment and move to class at a rate . may as well acquire infestation from the flea infested soil environment and progresses to class at a rate . Individuals in compartments and suffer a natural mortality rate and for the compartment they suffer a natural mortality at a rate and a disease induced mortality at a rate . The submodel of environmental component consists of two compartments, a compartment of larvae denoted by and a compartment of adult sand fleas denoted by . The larvae population are recruited into through shedding of jigger eggs by and at a constant rate ; therefore we have the total contribution of and from infested humans and animal reservoir populations, respectively. The larvae in compartment mature into adult jigger fleas at a maturation rate and undergo a natural death at a rate . The adult jigger flea population are recruited into through maturation by larvae at a rate and from infested animal reservoirs who contributes the fleas into the soil environment at a rate . The adult fleas leave the compartment when they attack the hosts at a rate and when they undergo a natural death at a rate and the additional death due to insecticides control at a rate , therefore we have the flea total death rate of .

The variables and parameters that describe the flow rates between compartments are given, respectively, in Notations. The possible interactions between humans, animal reservoirs, and flea infested environment with control measures are presented by the model flow diagram in Figure 1 and the differential equations describing the model are given in system (1).

Figure 1: The flow chart showing the dynamics of Tungiasis with control measures incorporated.

2.1. Model Flow Chart with Control Measures

The dynamical model with submodels of humans, animal reservoirs, and flea infested environment that incorporates time dependent control measures is presented hereunder.

2.2. Equations of the Model with Control Measures

From compartmental flow chart in Figure 1 the nonlinear differential equations representing the controlled system of Tungiasis dynamical model are given by

with initial condition

where

2.3. Optimal Control Problem for Tungiasis Epidemic

In this section, we present the optimal control problem considering the performance index and the controlled model equations with initial state conditions whereby the goal is to find the optimal levels of the control measures needed to minimize the number of infested humans, animal reservoirs, and fleas as well as the cost of implementing the control strategies ( for ). The objective functional to be minimized is given by

The terms , , , , in the objective functional are the number of infested populations that need to be minimized. The terms represents the cost of off-host insecticides spraying, represents the cost of on-host spraying of infested domestic animals, represent the cost of implementing environmental hygiene, represents the cost of education campaign, and represents the cost of personal protection. , , , , , , , , , and are positive balancing coefficients (weights) which regularize the optimal control. Quadratic expressions of the controls are included to indicate nonlinear costs potentially arising at high intervention levels [23].

The optimal control problem is formulated to obtain the minimum number of infested populations , , , , under minimum cost. Therefore the objective function in (4) is minimized subject to the model equations in (1). We seek the optimal controls , , , , such that

3. Model Analysis

The basic framework of an optimal control is to prove the existence of the optimal control and then characterize the optimal control through optimality system [24]. Given the optimal control problem in (5) we prove the existence of optimal control problem using the approach by Fleming and Rishel [25] and by Lukes [26] and then characterizing it for optimality.

3.1. The Existence of Optimal Control Problem

To prove the existence of optimal control the following conditions should be satisfied:

(i)The set of controls and corresponding state variables are nonempty.(ii)The control set is convex and closed.(iii)The right-hand side of the state system is bounded by a linear function in the state and control variables.(iv)The integrand of the objective functional is convex.(v)The integrand of the objective functional is bounded below by , where and .

Given that, , , , , , , , ,  and   are the state variables of the controlled model system (1) and , , , , are the control variables. We state and prove Theorem 1 as follows.

Theorem 1. There exists an optimal control variables set such that

Proof. Using the results in Fleming and Rishel [25] and by Lukes [26], the control and the state variables are nonnegative values. In this minimization problem, the necessary convexity of the objective functional in , , , , is satisfied. The set of admissible Lebesgue measurable control variables is also convex and closed by the definition. The optimal system is bounded which determines the compactness needed for the existence of optimal control. In order to verify this argument we use the approach adopted by Sadiq et al. [27] and Namawejje et al. [28], whereby system (1) is put in the following form:

where

is the matrix

where

and denotes the derivative of with respect to time . System (4) is a nonlinear system with a bounded coefficient. We set

The second term on the right-hand side of , , satisfies

where the positive constant is independent of the state variables. Also we have , where . So, it follows that the function is uniformly Lipchitz continuous. From the definition of control variables and nonnegative initial conditions we can see that a solution of the system (1) exists [29].
The integrand in the objective functional (4) which is given by the following equation is convex in the control set .
We must verify the condition that there exist a constant and positive numbers and such that

The last condition is satisfied when , , and . This ends the proof.

3.2. Determination of the Necessary Conditions for Optimality

The necessary conditions include the optimality solutions and the adjoint equations that an optimal control must satisfy which come from Pontryagin’s maximum principle (Pontryagin’s et al., 1962). This principle converts systems (1) and (4) into a problem of minimizing pointwise Hamiltonian function , which is formed by allowing each of the adjoint variables to correspond to each of the state variables accordingly and combining the results with the objective functional [21]. The resulting equation is as given by

where

  for are the adjoint functions associated with the state equations in (1),  for is the right-hand side of the differential equations of th state variable in system (1).

The expanded form of Hamiltonian function in (11) is given by

The optimality equations are obtained when taking the partial derivative of the Hamiltonian function with respect to the control variables , respectively, and the time derivative of adjoint equation which is obtained by taking the negative partial derivative of with respect to the model state variables such that .

Theorem 2. There exist an optimal control set and their corresponding state solutions , , , , , , , , and that minimize , and therefore there exist adjoint functions such that

where and with transversality conditions, , and the control variables satisfy the following optimality conditions:

Proof. The differential equations for the adjoints are standard results from Pontryagin’s maximum principle (1962). Given the Hamiltonian function in (12) the adjoint equations can be easily computed by

Therefore, the adjoint system evaluated at optimal controls , , , , and and the corresponding model state variables , , , , , , , , is as given by

where and with transversality conditions, , and the characterization of optimal controls , , , , ; that is, the optimality equations are based on the conditions:

Subject to (17) the optimality condition given the Lebesgue measurable control set , where the control variables , , , , are measurable functions, is given by

If we set at the results are the same as in characterization in (14).
In order to satisfy the given bounds for the control functions (i.e., and ) the optimal control is restricted to . Therefore using the bounds for the controls , , , , we have the following.
The bounds and the impact notation for the control are, respectively, given as

Similar four-step arguments hold for optimal control schedules , , , and in the same way as in (19) based on the characterization in (14).

3.3. The Optimality System

The optimality system consists of the state system and adjoint system with initial and transversal conditions together with characterization of optimal control. Any optimal control pair must satisfy this optimality system as indicated in (20) and (22), respectively.

with initial conditions

where and with transversality conditions .

4. Numerical Simulation of the Optimal Control Model, Results, and Discussion

Sometimes it may not be possible to solve the optimality system analytically; instead numerical methods are used to approximate the solutions and display the results. The optimality system is a two-point boundary problem, because of the initial condition of the state system and the terminal condition of the adjoint system [30]. To solve the optimality system with initial conditions for the states and final time conditions for the adjoints, we use the Runge-Kutta fourth-order procedure which is more accurate and elaborative technique. A Runge-Kutta method is a multiple-step method, where the solution at time is obtained from a defined set of previous values and is the number of steps. This method is described in a book by Lenhart and Workman [8] and it is known as forward-backward sweep method. The process begins with an initial guess on the control variable and given initial conditions for states, we approximate solutions for state equations using Runge-Kutta forward sweep method. Given the state solutions from previous step and the final time conditions for adjoints, we approximate solutions for adjoint equations using Runge-Kutta backward sweep method. The value of control variables is updated by averaging the previous value and the new value arising from the control characterization. The process is repeated for forward numerical scheme and updating of the controls until successive values of all states, adjoints, and controls are sufficiently close or converge.

4.1. Numerical Simulation of Optimal Control Model

We conduct numerical simulation in order to investigate the effects of the control strategies on the transmission dynamics of Tungiasis. The simulations are performed using MALAB, and we set time in days. The estimated initial values of model state variables are , , , , , , , , and and for the adjoint system we have terminal conditions , where we set days. The cost coefficients corresponding to state variables are estimated to be , , , , and . The quadratic cost coefficient corresponding to control measures is as well estimated to be , , , , and . We use a set of parameter values whose sources are from literature and others are estimated as shown in Table 1.

Table 1: Parameter values for numerical simulation of optimal control model.

Parameter
Value/range
Source/references

Estimated

Estimated

0.0105 per day
Estimated

0.011 per day
Estimated

0.037 per day
Estimated

0.000045 per day
UNICEF. [14]

per day
Gaff et al. [15]; Radostits. (2001)

0.04 per day
Eisele et al. [16],

0.58 per day
Estimated

0.08 per day
Estimated

0.19 per day
Estimated

0.48 per day
Estimated

0.052 per day
Gaff et al. [15]

0.26 (0.091–0.9) per day
Allerson et al. [17]

0.00011 per day
TP, [18].

0.022 per day
Gaff et al. [15]

0.40 per day
Estimated

0.12 per day
Estimated

0.4
Estimated

0.6
Estimated

0.15
Estimated

0.15
Estimated

0.09
Estimated

Source: Kahuru et al. [19].

We then plot the graphs to show the effects of the control measures when implemented under different combination options. We first illustrate the situation when no optimal control strategy is implemented as shown in Figures 2(a)–2(d) and then we suggest seven control strategies with different combinations of control measures and compare their performance in order to determine the best option to control the disease for maximum benefit.

(a)

(b)

(c)

(d)

Figure 2: Simulations of model variables , , without intervention.

4.2. Discussion of the Results

In this section, we present the results of the numerical simulation of our optimal control problem by discussing the implications of implementing the seven optimal control strategies on the Tungiasis dynamical model. To observe the effects of the optimal control strategies, we plot the results from simulation of the uncontrolled system and that from the controlled system together in Figures 3–9. To compare the effects of these options, we plot the results together in Figures 10(a)–10(c). Thus we consider the following seven combination options.

(a)

(b)

(c)

(d)

Figure 3: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a)

(b)

(c)

(d)

Figure 4: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a)

(b)

(c)

(d)

Figure 5: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a)

(b)

(c)

(d)

Figure 6: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a)

(b)

(c)

(d)

Figure 7: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a)

(b)

(c)

(d)

Figure 8: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a)

(b)

(c)

(d)

Figure 9: Optimal solutions for model variables , , and the control profiles for , , , , with .

(a) The comparison of the control strategies on severely infested humans

(b) The comparison of the control strategies on severely infested animals

(c) The comparison of the control strategies on the adult flea population

Figure 10

4.2.1. When All the Control Strategies Are Not Implemented (i.e., )

Figures 2(a)–2(c) show the situation where no control strategy is implemented to the dynamical system whereby the model trajectories represented by a dotted lines for severely infested humans, severely infested animal reservoir, and flea populations remain unchanged and the control profiles in Figure 2(d) show that all control measures are at the lower bound (i.e., ).

4.2.2. Strategy  1: All Control Measures Are Implemented (i.e., )

Under strategy  1, all the control measures are used to optimize the objective functional . In Figure 3(d), the control measure is at upper bound at the beginning and after 40 days it gradually drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 45 days it gradually drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 60 days it drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 140 days it rapidly drops to the lower bound at the final time and the control measure starts at the upper bound at the beginning and remains there until it drops to the lower bound.

4.2.3. Strategy  2: The Control Measures (i.e., and )

Under strategy  2, the control measures are used to optimize the objective functional . In Figure 4(d), the control measure is at upper bound at the beginning and after 140 days it rapidly drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 175 days it rapidly drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 195 days it drops to the lower bound at the final time. The control measure starts at the upper bound at the beginning and remains there until it drops to the lower bound at the final time. The control measure is at the lower bound at the beginning and remains there till the final time.

4.2.4. Strategy  3: The Control Measures (i.e., and )

Under strategy  3, the control measures are used to optimize the objective functional . In Figure 5(d), the control measure is at upper bound at the beginning and after 140 days it rapidly drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 190 days it rapidly drops to the lower bound at the final time. The control measure starts at the upper bound at the beginning and remains there until it drops to the lower bound. The control measures and start at the lower bound at the beginning and remains there till the final time.

4.2.5. Strategy  4: The Control Measures (i.e., and )

Under strategy  4, the control measures are used to optimize the objective functional . In Figure 6(d), the control measure is at upper bound at the beginning and after 80 days it gradually drops to the lower bound at the final time. The control measure starts at the upper bound at the beginning and remains there until it drops to the lower bound. The control measures, , , and , start at the lower bound at the beginning and remain there till the final time.

4.2.6. Strategy  5: The Control Strategies (i.e., and )

Under strategy  5, the control measures are used to optimize the objective functional . In Figure 7(d), the control measure is at upper bound at the beginning and after 50 days it gradually drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 60 days it gradually drops to the lower bound at the final time. The control measure starts at the upper bound at the beginning and remains there until it drops to the lower bound at the final time. The control measures and start at the lower bound at the beginning and remains there till the final time.

4.2.7. Strategy  6: The Control Measures (i.e., and )

Under strategy  6, the control measures are used to optimize the objective functional . In Figure 8(d), the control measure is at upper bound at the beginning and after 140 days it rapidly drops to the lower bound at the final time. The control measures and start at the upper bound at the beginning and remain there until they drop to the lower bound. The control measures and start at the lower bound at the beginning and remain there till the final time.

4.2.8. Strategy  7: The Control Measures (i.e., and )

Under strategy  7, the control measures are used to optimize the objective functional . In Figure 9(d), the control measure is at the upper bound at the beginning and after 40 days it gradually drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 75 days it gradually drops to the lower bound at the final time. The control measure is at the upper bound at the beginning and after 170 days it rapidly drops to the lower bound at the final time. The control measure starts at the upper bound at the beginning and remains there until it drops to the lower bound at the final time and the control measure starts at the lower bound at the beginning and remains there till the final time.

Generally, it is observed in Figures 3(a)–3(c), 4(a)–4(c), 5(a)–5(c), 6(a)–6(c), 7(a)–7(c), 8(a)–8(c), and 9(a)–9(c) that the controlled trajectories represented by the solid lines decrease compared to uncontrolled trajectories represented by the dotted lines. This implies that the control strategies have the positive impacts which lead to a reduction in the number of infested humans and animals and also reduces the sand flea population in the soil environment.

4.2.9. The Comparison of the Control Strategies

To compare the performance of the control strategies under consideration, we plot the results on same graphs indicating the effects of the control strategies as indicated in Figures 10(a)–10(c).

From Figures 10(a)–10(c), it is observed that the control strategy involving the combination of all the five control measures , , , , has the significant impact on the reduction of disease transmission because it lowers the severely infested humans and animals and the sand flea populations to minimal levels compared to other control strategies. The red solid line as depicted in Figures 10(a), 10(b), and 10(c) represents the effect caused by the control strategy  1: with the combination of all five control measures. This is the best control strategy because it dominates both graphs as it decreases rapidly compared to other trajectories. Moreover we have observed that those control strategies, whose combinations involve the control measure based on insecticides applications to the premises , yield better results. These are control strategies  5, 6, and 7. But control strategy  5 denoted by the control measures performs better in the same way as the control strategy  1 denoted by the control measures

5. Conclusion

In this paper, the optimal control techniques have been applied on Tungiasis dynamical model with control strategies. We defined the control set including controlling the transmission of infestation from flea infested soil environment to human population, from the flea infested soil environment to animal population, and from flea infested animal to human population, controlling flea infested soil environment, and controlling the flea infested animal population. We proved the existence of optimal control problem and determined the necessary conditions for optimality using Pontryagin’s maximum principle which converts constrained optimization problem into unconstrained Hamiltonian function whereby optimality and adjoint equations are obtained. We, lastly, performed numerical simulations of the resulting control problem to investigate the effects of the control strategies under consideration and compare their performances. The numerical results showed that the control strategy that comprises all five control measures and that with control measures have the significant impact on reduced Tungiasis transmission and those control strategies involving insecticides control to the premises yielded better results, which implies that the insecticides application control is more effective than other individual control measures. In poor rural communities where resources are always scarce, we suggest that the combination option involving the controls of focal insecticides spraying, insecticidal dusting on animal furs, and environmental hygiene should be adopted, having observed from the comparison of all seven control strategies in Figures 10(a), 10(b), and 10(c) that there is no significant difference between this strategy and the strategy that involves the combination of the five control measures . Among others, poor housing conditions and the presence of domestic and sylvatic animals on the home compound are risk factors. Therefore, controlling of infested soils and animal reservoirs with insecticides control, ant-flea compounds or animal furs, environmental hygiene, and cementing the floors of houses may serve as a possible approach to control the epidemic and, thus, to fight against Tungiasis infestation in endemic settings multidimensional control process should be employed in order to achieve the maximum benefits.

Notations

The State Variables of the Model with Control Measures (Source: Kahuru et al. [19])

: | Number of humans in a susceptible class at time
: | Number of animals in a susceptible class at time
: | Number of humans in mildly infested class at time
: | Number of animals in mildly infested class at time
: | Number of humans in severely infested class at time
: | Number of animals in severely infested class at time
: | Number of humans who are receiving treatments at time
: | The density of fleas population in the environment at time
: | The density of larvae population in the environment at time
, : | Total human and animal reservoirs populations at time

The Parameters of the Model with Control Measures (Source: Kahuru et al. [19])

: | Maximal larval carrying capacity
: | Half saturation constant
: | Maturation (transition) rate from larvae to adult jigger fleas
, : | Disease induced death rates for humans and animal reservoirs, respectively
, , , : | Natural mortality rates for humans, animals, fleas, and larvae, respectively
: | The rate of removal of jigger fleas that leaves the soil to attack the hosts
: | Effective contact rate between environment and humans
: | Effective contact rate between environment and animal reservoirs
: | Effective contact rate between infested animals and susceptible humans
: | Effective contact rate between infested animals and susceptible animals
, : | Recruitment birth rates for humans and animal reservoirs, respectively
: | The rate of flea eggs deposit on the ground
: | Shedding rates for adult fleas into the environment
, : | The proportions of infestation for humans and animals, respectively
: | Progression rate of treated humans to susceptible class
, : | Progression rates from infested human classes to the treatment class.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

D. Pilger, S. Schwalfenberg, J. Heukelbach et al., “Controlling tungiasis in an impoverished community: an intervention study,” PLoS Neglected Tropical Diseases, vol. 2, no. 10, article e324, 2008.

J. Heukelbach, T. Wilcke, G. Harms, and H. Feldmeier, “Seasonal variation of tungiasis in an endemic community,” American Journal of Tropical Medicine and Hygiene, vol. 72, no. 2, pp. 145–149, 2005.

J. Heukelbach, A. M. L. Costa, T. Wilcke, N. Mencke, and H. Feldmeier, “The animal reservoir of Tunga penetrans in severely affected communities of north-east Brazil,” Medical and Veterinary Entomology, vol. 18, no. 4, pp. 329–335, 2004.

G. Collins, T. McLeod, N. I. Konfor, C. B. Lamnyam, L. Ngarka, and N. L. Njamnshi, “Tungiasis: a neglected health problem in rural cameroon,” International Journal of Collaborative Research on Internal Medicine and Public Health, vol. 1, no. 1, pp. 2–10, 2009.

P. Nordin, M. Thielecke, N. Ngomi, G. M. Mudanga, I. Krantz, and H. Feldmeier, “Treatment of tungiasis with a two-component dimeticone: a comparison between moistening the whole foot and directly targeting the embedded sand fleas,” Tropical Medicine and Health, vol. 45, no. 1, 2017.

M. Ozair, A. A. Lashari, I. H. Jung, and K. O. Okosun, “Stability analysis and optimal control of a vector-borne disease with nonlinear incidence,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 595487, 21 pages, 2012.

B. Seidu and O. D. Makinde, “Optimal control of HIV/AIDS in the workplace in the presence of careless individuals,” Computational and Mathematical Methods in Medicine, Article ID 831506, 19 pages, 2014.

S. M. Lenhart and J. T. Workman, Optimal Control Applied to Biological Models, CRC Press, 2007.

H. S. Rodrigues, M. T. T. Monteiro, and D. F. Torres, “Optimal control and numerical software: an overview,” https://arxiv.org/abs/1401.7279.

H. R. Joshi, S. Lenhart, M. Y. Li, and L. Wang, “Optimal control methods applied to disease models,” in Mathematical studies on human disease dynamics, vol. 410 of Contemp. Math., pp. 187–207, Amer. Math. Soc., Providence, RI, 2006.

P. D. Roberts and V. M. Becerra, “Optimal control of a class of discrete-continuous non-linear systems—decomposition and hierarchical structure,” Automatica. A Journal of IFAC, the International Federation of Automatic Control, vol. 37, no. 11, pp. 1757–1769, 2001.

H. Schättler and U. Ledzewicz, “The Pontryagin Maximum Principle: From Necessary Conditions to the Construction of an Optimal Solution,” in Geometric Optimal Control, vol. 38 of Interdisciplinary Applied Mathematics, pp. 83–194, Springer New York, New York, NY, 2012.

P. Ruhnau and C. Schnörr, “Optical stokes flow estimation: an imaging-based control approach,” PAMM, vol. 6, no. 1, pp. 863-864, 2006.

UNICEF, “State of The World’s Children 2015 Country Statistical Information,” 2015.

H. D. Gaff, D. M. Hartley, and N. P. Leahy, “An epidemiological model of rift valley fever,” Electronic Journal of Differential Equations, vol. 115, pp. 1–12, 2007.

M. Eisele, J. Heukelbach, E. Van Marck et al., “Investigations on the biology, epidemiology, pathology and control of Tunga penetrans in Brazil: I. Natural history of tungiasis in man,” Parasitology Research, vol. 90, no. 2, pp. 87–99, 2003.

M. Allerson, J. Deen, S. E. Detmer et al., “The impact of maternally derived immunity on influenza A virus transmission in neonatal pig populations,” Vaccine, vol. 31, no. 3, pp. 500–505, 2013.

Tanzania population 2016, “Tanzania Population 2017,” http://worldpopulationreview.com/countries/tanzania-population/.

J. Kahuru, L. Luboobi, and Y. Nkansah-Gyekye, “Modelling the dynamics of Tungiasis transmission in zoonotic areas,” Journal of Mathematical and Computational Science, vol. 7, no. 2, pp. 375–399, 2017.

E. Bonyah, I. Dontwi, and F. Nyabadza, “Optimal control applied to the spread of buruli uclcer disease,” American Journal of Computational and Applied Mathematics, vol. 4, no. 3, pp. 61–67, 2014.

A. O. Isere, J. E. Osemwenkhae, and D. Okuonghae, “Optimal control model for the outbreak of cholera in Nigeria,” African Journal of Mathematics and Computer Science Research, vol. 7, no. 2, pp. 24–30, 2014.

G. Devipriya and K. Kalaivani, “Optimal control of multiple transmission of water-borne diseases,” International Journal of Mathematics and Mathematical Sciences, vol. 2012, Article ID 421419, 2012.

R. L. M. Neilan, Optimal Control Applied to Population and Disease Models [Ph.D. Dissertation], University of Tennessee, Knoxville, Tennessee, Tenn, USA, 2009.

D. Kirschner, S. Lenhart, and S. Serbin, “Optimal control of the chemotherapy of HIV,” Journal of Mathematical Biology, vol. 35, no. 7, pp. 775–792, 1997.

W. H. Fleming and R. W. Rishel, Deterministic and Stochastic Optimal Control, vol. 1, Springer, New York, NY, USA, 1975.

D. L. Lukes, Differential Equations: Classical to Controlled, Mathematics in Science and Engineering, Academic Press, New York, NY, USA, 1982.

S. F. Sadiq, M. A. Khan, S. Islam, G. Zaman, H. Jung, and S. A. Khan, “Optimal control of an epidemic model of leptospirosis with nonlinear saturated incidences,” Annual Research & Review in Biology, vol. 4, no. 3, pp. 560–576, 2014.

H. Namawejje, L. S. Luboobi, D. Kuznetsov, and E. Wobudeya, “Modeling optimal control of rotavirus disease with different control strategies,” Journal of Mathematical and Computational Science, vol. 4, no. 5, p. 892, 2014.

G. Birkhoff and G. Rota, Ordinary Differential Equations, John Wiley & Sons, New York, NY, USA, 4th edition, 1989.

K. Blayneh, Y. Cao, and H.-D. Kwon, “Optimal control of vector-borne diseases: treatment and prevention,” Discrete and Continuous Dynamical Systems B, vol. 11, no. 3, pp. 587–611, 2009.