Effect of public awareness, behaviours and treatment on infection-age-structured of mathematical model for HIV/AIDS dynamics

Nurudeen Oluwasola Lasisi1

1Department of Mathematics and Statistics, Federal Polytechnic, Kaura Namoda, Nigeria

Mathematical Models in Engineering, Vol. 6, Issue 2, 2020, p. 103-121. https://doi.org/10.21595/mme.2020.21249
Received 19 December 2019; accepted 31 March 2020; published 30 June 2020

Copyright © 2020 Nurudeen Oluwasola Lasisi. 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.
Creative Commons License
Table of Contents Download PDF References
Cite this article
Views 36
Reads 18
Downloads 326
CrossRef Citations 0
Abstract.

Since the discovering of the HIV/AIDS epidemic, more people have been infected with the HIV/AIDS infection and several individuals have died of HIV/AIDS virus. In this research work, an infection-age-structured mathematical model for HIV/AIDS transmission dynamics is developed and investigated, taking into consideration the public awareness, Treatment and Sex behaviours. It is assumed that the infectious population is structured according to time and age of infection. We established the criteria for the existence and uniqueness of solution. An explicit formula for the effective reproduction number RE of the model is obtained. We showed that the disease-free-equilibrium (DFE) state is stable if RE<1. Constructing a Lyaponov function, the global stability of the DFE state of the system is obtained for RE<1. Endemic Equilibrium (EE) state was analyzed locally and globally with the method of linear and nonlinear Lyapunov function; it was found that EE is locally asymptotically stable if RE>1 and globally stable if RE>1 for X1<X2.

Highlights
  • An infection-age-structured model is used to study the effect of HIV public awareness program, with unaware HIV susceptible on HIV/AIDS dynamics.
  • The public awareness, Treatment and Sex behaviours are considered. The criteria for the existence and uniqueness of solution were obtained.
  • It was discovered that the model have two equilibrium; Disease Free Equilibrium (DFE) and Endemic Equilibrium (EE). An explicit formula for the effective reproduction number of the model is obtained.
  • We analysed the invariant region and shown the dynamics of model equations in the region Ω, The model equations was considered as been epidemiologically and mathematically well posed.
  • The results from stability analysis showed that DFE is locally asymptotically stable when R_E<1 and globally asymptotically stable when R_E<1.
  • Endemic Equilibrium was found to be locally asymptotically stable when R_E>1 and globally asymptotically stable when R_E>1 for X(1)<X(2).

Keywords: age-of-infection, equilibrium, mathematical model, HIV/AIDS, public awareness.

1. Introduction

Human Immunodeficiency Virus (HIV) Infection is a virus that destroys the immune system in an individual’s system which natural defence against illness. There are two types of HIV, which are HIV 1 and HIV 2, HIV 1 is the more fatal and deadly which is the cause of the Acquired Immunodeficiency Syndrome (AIDS) pandemic while, HIV 2 is also known to cause AIDS [1]. Symptoms of the virus are diarrheal, fever, constant cough, and skin problems usual infections, severe illnesses [2]. The virus infection spread from person to person including having unguarded sex with an infected person, sharing of the followings, needles, syringes as well as other items for injection drug use with an infected person. Blood transfusion or organ and tissue transplant can also transmit the virus [1]. The Virus does not spread through the following, skin to skin contact, shaking of hands, as well as kissing, through Air and water, drinking items, saliva, tears sharing of a toilet, towels and mosquitoes or insects biting [3].

There are stages of HIV infection which are Acute HIV infection is the first few weeks after infection. The second stage is called clinical latency is after the first month, HIV enters the clinical latency stage [3]. The third stage, AIDS which we called acquired immunodeficiency syndrome. Without treatment, HIV advance to AIDS, at this point, the immune system is too weak to fight off life-threatening disease and infection [1]. Whereas, untreated life expectancy with AIDS is about three years. HIV/AIDS destroys T cells called CD4 cell, these cells help immune system fight infections of individuals. More so, healthy adults have a CD4 count of 800 to 1,000 per cubic millimetres. If CD4 count falls below 200 per cubic millimetres, then it is likely to have progressed to AIDS [3]. AIDS weakens the immune system to the point where it can no longer fight off most diseases and infections. This makes individuals more vulnerable to a wide range of illnesses [3].

Meanwhile, more than 70 million individuals have been infected with the HIV virus and about 35 million individuals have died of HIV [1]. Therefore, Sub-Saharan Africa remains most severely affected, with nearly 1 in every 25 adults (4.2 %) living with HIV. Mathematical models have been used frequently to study both transmission and infectious dynamics of virus [4]. Meanwhile, there are numerous works on Transmission of HIV/AIDS Infection. In this research work, we considered infection-age-structured of mathematical model of HIV/AIDS dynamics incorporating public awareness, behaviours and treatment.

2. Model

We considered an infection-age-structured mathematical model for transmission dynamic of HIV/AIDS infection incorporating HIV public awareness, sexual behaviours of some individuals and treatment. The transmission dynamics of HIV was developed using ordinary and partial differential equations and integro-differential equation. The total population (N) is partitioned into five (5) compartments, namely; unaware susceptible individuals (S), aware susceptible individuals whose modify their sex behaviours (S1), aware susceptible individuals whose abstain and remain faithful to uninfected sex partners (S2), infected individuals class (I) and removed individuals class (R).

S t is generated by a recruitment rate of Λ and the natural death rate for all compartments is μ. As a result of counselling and awareness, a proportion of the susceptible leave the S class at the rate of ψ(0<ψ<1)which is the individuals abstained and remained faithful to uninfected sex partners and a fraction of 1-ψ is the rate at which individuals modify their sex behaviours. Also, ρ is the rate S of individuals knowing their HIV status. The expected reduction in unsafe sex behaviours by the S1 class as a result of awareness is accounted for by ε(0<ε<1). Thus, λ1 is the force of infection from S to I, where λ1=β1I(t) and β1 is the incidence rate of unaware susceptible (effective infection contact rate). λ2 is the force of infection from S1 to I, where λ2=β2I(t) and β2 is the incidence rate of aware susceptible. The infection spreads through direct contact between S and S1 with infected class (I). The disease induced death rate is α(τ) and received treatment at the rate of δ(τ), then moved to removed class (R), the R class whose are assumed in this study not be involved in transmission as a result of efficacy of ART and has additional disease induced death rate of υ, It is structured by the infection-age with the density function i(t,τ) following the idea of [3], where t is the parameter and τ is the infection-age. There is a maximum infection age T at which a number of the infected individuals form infected class leave the class via death, 0τT. Meanwhile, since the infection period of HIV is long, then, the population is not constant, it is varying. A schematic diagram for the model is shown in Fig. 1.

Fig. 1. Schematic diagram of HIV/AIDS infection

 Schematic diagram of HIV/AIDS infection

These assumptions lead to the following system of nonlinear integro-partial differential equations with non-local boundary conditions, which describes the dynamics of the transmission of the disease as:

(1)
d S d t = Π - μ S - λ 1 S - ψ ρ S - 1 - ψ ρ S ,
(2)
d S 1 d t = 1 - ψ ρ S - μ S 1 - ε λ 2 S 1 ,
(3)
d S 2 d t = ψ ρ S - μ S 2 ,
(4)
i ( t , τ ) t + i ( t , τ ) τ + ( μ + α ( τ ) + δ ( τ ) ) i ( t , τ ) = 0 ,
(5)
d R d t = 0 T δ i ( t , τ ) d τ - μ R - υ R ,

where:

(6)
λ 1 = 0 T β 1 i ( t , τ ) d τ ,
(7)
λ 2 = 0 T β 2 i ( t , τ ) d τ ,
(8)
i t , 0 = B t = λ 1 S + ε λ 2 S = β 1 S t I t + ε β 2 S t I t .

The total Population size N is given by:

(9)
N = S + S 1 + S 2 + I + R ,

where:

(10)
I t = 0 T i t , τ d τ ,
(11)
S 0 = S 0 ,           S 1 0 = S 10 ,             S 2 0 = S 20 ,             I 0 = I 0 ,             R 0 = R 0 ,             N 0 = N 0 ,
(12)
i 0 , τ = ϕ τ .

Integration of Eq. (4) over τ is setting as:

(13)
0 T i t , τ t + i t , τ τ   d τ + 0 T ( μ + α ( τ ) + δ ( τ ) ) i ( t , τ ) d τ = 0 ,
d I d t + 0 T i t , τ τ d τ + ( μ + α ( τ ) + δ ( τ ) ) 0 T i t , τ d τ = 0 ,
d I d t + 0 T i ( t , τ ) + ( μ + α ( τ ) + δ ( τ ) ) 0 T i t , τ d τ = 0 ,
d I d t + i ( t , τ ) 0 T + ( μ + α ( τ ) + δ ( τ ) ) I ( t ) = 0 ,
d I d t + i ( t , T ) - i ( t , 0 ) + ( μ + α ( τ ) + δ ( τ ) ) I ( t ) = 0 .

Substituting Eqs. (8), (10) into Eq. (13) yields:

(14)
d I d t + i ( t , T ) - β 1 S ( t ) I ( t ) - ε β 2 S ( t ) I ( t ) + ( μ + α ( τ ) + δ ( τ ) ) I ( t ) = 0 .

T is the maximum infection age and if τ=T then the infected person die of the disease i.e.:

(15)
i ( t , T ) = 0 .

Implies:

(16)
d I d t = β 1 S t I t + ε β 2 S t I t - μ + α τ + δ τ I t .

3. Basic properties of the model

The addition of Eqs. (1), (2), (3), (5) and (16) is given by:

(17)
d N d t = d S d t + d S 1 d t + d S 2 d t + d I d t + d R d t ,
d N d t = Λ - μ N - α τ I - υ R ,
(18)
d N d t Λ - μ N .

We apply theorem of differential inequality in [5] and separation of differential inequality of Eq. (18), we have:

(19)
d N Λ - μ N = d t ,
(20)
N t Λ μ - Λ - μ N 0 μ e - μ t .

Thus, as t, we have:

(21)
0 N t N 0 - Λ μ e - μ t + Λ μ .

Implies, as limt, N(t)<Λ/μ, where N(0)<Λ/μ.

However, if N(0)Λ/μ, Nt will decrease to Λ/μ. So Nt is a bounded function of time. We can say that Ω is bounded and at limiting equilibrium t, N(t)=Λ/μ. Besides any sum or difference of variables in Ω with positive initial values will remain in Ω or in a neighbourhood of Ω. Thus Ω is positively invariant and attracting with respect to the Eqs. (1)-(5) which are epidemiologically meaningful and mathematically well posed.

4. Existence and uniqueness of solution

We shall first and foremost establish the criteria for the existence and uniqueness of solution of the model Eqs. (4), (7) and (12). We define the derivative part of Eq. (4) as:

(22)
i ( t , τ ) t + i ( t , τ ) τ = l i m h 0 i ( t + h , τ + h ) - i ( t , τ ) h .

Let:

(23)
q h = i t 0 + h , τ 0 + h .

Then, the ordinary differential equation:

(24)
d q ( h ) d h + ( μ + α ( τ ) + δ ( τ ) ) q ( h ) = 0 ,

has the unique solution:

(25)
q h = q 0 e - 0 h μ + α s + δ s d s .

Using Eq. (23) in Eq. (25), we have:

(26)
i t 0 + h , τ 0 + h = i t 0 , τ 0 e - 0 h μ + α s + δ s d s ,

which gives the value of i at all point on the characteristics passing through the point i(t0,τ0), when i(t0,τ0)=(0,τ-t); h=t, we have that:

(27)
i t , τ = ϕ τ - t e - 0 h μ + α s + δ s d s ,             t τ .

Also, when i(t0,τ0)=(t-τ,0); h=t, we have:

(28)
i t , τ = B t - τ e - 0 h μ + α s + δ s d s ,           t > τ .

The Eqs. (27) and (28) give the solution of i(t,τ) in the positive quadrant t0; τ0.

Thus, problem Eqs. (4)-(12) has the following unique solution that exists for all time:

(29)
i ( t , τ ) = B t - τ π τ , t > τ , ϕ τ - t π τ , t < τ .

Eq. (29) is rewritten in the form:

(30)
i ( t , τ ) = B t - τ π τ , t > τ , ϕ ( τ - t ) π ( τ ) π ( τ - t ) , t < τ .

where:

(31)
π τ = e - 0 τ μ + α s + δ s d s = e - μ τ e - 0 h α s + δ s d s .

Eq. (31) denotes the survival probability, the probability for an individual to survive to age τ; thus π(τ)=0.

The expected life of an infected individual given by Eq. (26).

The function:

(32)
K τ = ϕ τ π τ ,           τ 0 , T .

Let, ϕ() is non-negative and belong to LT(0,τ), α() is non-negative and belong to L1loc([0,τ)), and

(33)
0 τ ( μ + α ( s ) + δ ( s ) ) d s = T ,             ϕ L 1 0 , τ ,               ϕ τ 0     a.e in     0 , τ .

We rewrite Eq. (8) in the form:

(34)
B t = φ 1 τ + φ 2 τ 0 T i t , τ d τ ,

where:

(35)
φ 1 = β 1 S ( t ) I ( t ) φ 2 = ε β 2 S ( t ) I ( t ) φ φ 1 , φ 2 .

Eq. (35) are also non-negative and belongs to LT(0,τ).

Substituting Eq. (30) into Eq. (34), we have for tτ:

(36)
B t = φ 1 τ + φ 2 τ 0 T B t - τ π τ d τ + 0 T ϕ t - τ π τ π τ - t d τ .

And for t>τ:

(37)
B t = 0 T φ 1 τ + φ 2 τ B t - τ π τ d τ .

Thus, Bt satisfies the following Volterra Integral equation of the second kind:

(38)
B t = F K t + 0 t K t - s B s d s .

With:

(39)
F t = 0 T φ τ π τ π τ - t ϕ τ - t d τ = 0 T φ τ + t π τ + t π τ ϕ τ d τ ,

where:

(40)
φ 1 τ + φ 2 τ φ τ ,
(41)
K t = φ t π t ,

where t0 and the function φ(φ1,φ2), π and ϕ are extended by zero outside the interval [0, T]. Eq. (38) is known as the renewal equation and also as the Lotka equation, we observed that the kernel Kt is an infectious function of newly infected individuals defined in Eq. (32). Eq. (38) is equivalent to Eqs. (4), (8), and (12); actually Eq. (38) is the main tool to investigate the existence of Eqs. (4), (8) and (12), the connection being provided by Eqs. (30), (39) and (41) respectively. The following Lemma states some properties of Eq. (38) on the basis of the assumptions Eq. (33).

Lemma 1: Let Eq. (33) be satisfied, then:

(42)
K t 0     a.e     K t = 0 ,           t > τ ,               K L 1 R + L T R + ,
(43)
F t 0 ,           K t = 0 ,             t > τ ,               F R + .

If moreover, according to [6]:

(44)
ϕ W 1,1 0 , τ ,           μ + α + δ ϕ L 1 0 , τ ,             F W 1 , T .

Proof:

Eq. (42) and the first part Eq. (43) are obvious. To prove that F(R+), we take t00 then we have:

(45)
F t = 0 T φ τ π τ π τ - t ϕ τ - t - ϕ τ - t 0 d τ + 0 T φ τ π τ π τ - t ϕ τ - t 0 d τ ,

where, ϕL1(R):

(46)
0 T φ ( τ ) π ( τ ) π ( τ - t ) ( ϕ ( τ - t ) - ϕ ( τ - t 0 ) ) d τ φ L T 0 T ϕ ( τ - t ) - ϕ ( τ - t 0 ) ) d τ .

As tt0, so that:

(47)
l i m t t 0 F t = 0 T φ τ π τ π τ - t ϕ τ - t 0 d τ = F t 0 .

Similarly, Eq. (44) implies FW1,T(R+). Thus, completes the proof. In the next theorem, we establish the existence and uniqueness of Eqs. (4), (8) and (12) through Eq. (38). Now we study Eqs. (4), (8) and (12) by considering Eqs. (38)-(41). Since the function F is continuous for the purpose of this section it is sufficient to study Eq. (38) in the class of continuous functions. Thus we introduce the following definition of a solution.

Definition 1. A Solution to Eq. (38) is a function B(R+,R) satisfying Eq. (38). First we have the following theorem, which is standard in the theory of Volterra equations. We provide a proof here for the sake of completeness.

Theorem 1. Let Eq. (33) be satisfied, then Eqs. (38)-(41) has a unique solution B(R+,R) such that B(t)0 for all t in addition if ϕ satisfied Eq. (44), then:

(48)
B W   l o c 1 , T R + ,             B 1 t = F k 1 t + K t B 0 + 0 t K t - s B 1 s d s .

Proof: first we assume that:

(49)
K L 1 ( R + ) = 0 T K ( s ) d s < 1 .

The solution of Eq. (38) is obtained using the standard Picard iteration defined, for t0 by:

(50)
B k + 1 t = F t + 0 t K s B k k - s d s .

And initialized by B0(t)=0 on IR+.

Therefore:

(51)
B 0 t = F t , B k + 1 = F t + 0 t K t - s B k s d s .

If we take any T>0, then by Eqs. (42)-(43) we obtained Bk([0,T]) and Bk(T)0, moreover:

(52)
B k + 1 ( t ) - B k ( t ) 0 t K t - s B k s - B k - 1 s d s ,
(53)
B k + 1 - B k C ( [ 0 , T ] ) K L 1 R + B k - B k - 1 C 0 , T .

Thus by Eq. (49) the sequence Bk(t) converges uniformly on [0,T] to a solution of Eq. (38) BC([0,T]) and B(t)0.

Concerning Uniqueness of this solution of the model, we set Bt and B-t to be the two solutions of Eq. (38), so that:

(54)
B - B ¯ C ( [ 0 , T ] ) K L 1 R + B - B ¯ C 0 , T .

This implies that by Eq. (49), we have B(t)=B¯(t).

In addition, if ϕ satisfied Eq. (44) then by Lemma 1 and Eq. (51), we have BkW1,T(R+) and setting Vk(t)=ddtBk(t) a.e we have VkLT(R+) and

(55)
V k + 1 t = F 1 t + K t F 0 + 0 t K t - s V k s d s .

This yield:

(56)
V k + 1 - V k L T ( R + ) K L 1 R + V k - V k - 1 L T R + .

Thus, again by Eq. (49) the sequence Vk converges in LT(R+) to V(t)=ddtB(t) are therefore Eq. (48) following from Eq. (55). In addition to Eq. (49), to make our argument valid we take α:

(57)
0 T e - α t K ( t ) d t < 1 .

Setting B¯=e-αtB(t), F¯=e-αtF(t), K¯=e-αtK(t), Eq. (38) is transformed into an equation to the follows:

(58)
B ¯ t t = F ¯ t + 0 t K ¯ t - s B ¯ s d s .

And since Eq. (58) satisfies Eq. (49), this is similar with Eq. (56). The next theorem allows us to state results for problem Eqs. (4), (8) and (12) through Eq. (30).

Theorem 2. Let Eq. (33) and Eq. (44) be satisfied, and we also assume that:

(59)
ϕ 0 = 0 T φ τ ϕ τ d τ .

And let i(t,τ) defined by Eq. (30), where Bt is the solution of Eq. (38)-(40), then:

(60)
i 0 , τ × R + ,             i t , τ 0 ,             ( μ + α ( ) + δ ( ) ) i ( t , ) L 1 ( 0 , τ ) t > 0 ,
(61)
i ( t , τ ) t + i ( t , τ ) τ   exist     a.e     in      0 , τ × 0 , T .

And problem Eqs. (4), (8) and (12) is satisfied. Moreover, according to [7], i(t,τ) is the only solution in the sense of Eqs. (60, 61).

Proof: the proof of Eqs. (60, 61) is quite straightforward and follows from the properties of Bt which is stated in Theorem 1. We only note the following inequality concerning the last part of Eq. (51):

(62)
0 τ ( μ + α ( τ ) + δ ( τ ) ) i ( t , τ ) d τ = 0 t τ ( μ + α ( τ ) + δ ( τ ) ) B ( t - τ ) π ( τ ) d τ
            + t τ τ μ + α τ + δ τ ϕ τ - t π τ π τ - t d τ
            M a x S [ 0 , t ] B ( t ) 0 t τ ( μ + α ( τ ) + δ ( τ ) ) π ( τ ) d τ
            + e 0 ( t τ ) - t ( μ + α ( s ) + δ ( s ) ) d s ϕ C ( 0 , τ ) t τ τ ( μ + α ( s ) + δ ( s ) ) d s
(63)
            M a x S 0 , t B t + e 0 t τ - t ( μ + α ( s ) + δ ( s ) ) d s ϕ C 0 , τ .

Eq. (63) provides the fact that Eq. (59) is intended to guarantee the continuity of i(t,τ) through the line τ=t; so that:

(64)
B 0 = 0 T φ τ ϕ τ d τ = ϕ 0 .

We note that the solution of problem Eqs. (4), (8) and (12) must be in the form of Eq. (30) with Bt satisfying Eqs. (38)-(40) such that Eqs. (30) and (44) are enough to provide a classical solution in the next following theorem.

Theorem 3. Let Eq. (33) be satisfied, then i(t,τ) defined by Eq. (30) has the following properties:

(65)
i t , C 0 , T ; L 1 0 , τ ,         i t , τ 0    a.e    in     0 , τ × R + ,
(66)
i ( t , ) L 1 e t φ L T ϕ L 1 ,
(67)
i ( t , τ )  is continuous for  τ<t  and satisfies Eq.8 for t>0,
(68)
l i m h 0 1 h i t + h , τ + h - i t , τ = - μ + α τ + δ τ i t , τ   a.e in   0 , τ × R + .

Proof: let us prove Eq. (66) first. From Eq. (39) we have:

(69)
F t φ L T ϕ L 1 ,           K t φ L T .

Then, from Eq. (38):

(70)
B t φ L T ϕ L 1 + φ L T 0 t B s d s .

Thus, by classical Gronwall’s inequality:

(71)
B t φ L T e t φ L T ϕ L 1 .

From this estimate, Eq. (30) yields:

(72)
i ( t , ) L 1 = 0 t B ( t - τ ) π ( τ ) d τ +   0 t ϕ ( τ ) π ( τ + t ) π ( τ ) d τ
            φ L T 0 t e t - τ φ L T d τ + 1 ϕ L 1 = e t φ L T ϕ L 1 .

Eq. (65) follows easily from Eq. (67). Now from a given ϕL1(0,τ), let ϕn be a sequence such that ϕn satisfy Eqs. (44) and (59):

(73)
l i m n T ϕ n - ϕ L 1 = 0 .

And let in be the solution of Eqs. (4), (8) and (12) corresponding to ϕn. Thus inC([0,T];L1(0,τ)) and by Eq. (66) and linearity, we have:

(74)
i n ( t , ) - i ( t , ) L 1 e t φ L T ϕ n - ϕ L 1 .

So that i is the limit of the sequence in in the Space C([0,T];L1(0,τ)) i.e Eq. (65) is true. Finally, Eqs. (67)-(68) are straightforward and that completes the proof.

This shows that even when the initial age of infection ϕ is not regular the solution i(t,τ) still has some regularity. We also note that the estimate Eq. (66) provides continuity of the solution iwith respect to initial age of infection ϕ. Hence, model Eqs. (4), (8) and (12) exist and are unique and well posed mathematically. In the norm of the space L1(0,τ): this illustrates the existence and uniqueness of model Eqs. (4), (8) and (12) which is in agreement with the biological meaning of the age of infection i(t,τ) for the period of disease.

In the next theorem, we used a similar approach to the work of [2] to prove the existence and uniqueness of solution for the model Eqs. (1), (2), (3) and (5).

Theorem 4. The exists of a unique solution of Eqs. (1), (2), (3) and (5) for t-t0a, S-S0b, S1-S10c,S2-S20d, I-I0e, R-R0f.

Proof:

(75)
f 1 S , S 1 , S 2 , I , R = Λ - μ S - β 1 S I - ψ ρ S - 1 - ψ ρ S ,               S 0 = S 0 ,
(76)
f 2 S , S 1 , S 2 , I , R = 1 - ψ ρ S - μ S 1 - ε β 2 S I ,             S 1 0 = S 10 ,
(77)
f 3 S , S 1 , S 2 , I , R = ψ ρ S - μ S 2 ,             S 2 0 = S 20 ,
(78)
f 4 S , S 1 , S 2 , I , R = δ I - μ R - υ R ,             R 0 = R 0 .

Then:

(79)
f 1 S = - ( μ + β 1 I + ρ ) ,             f 1 S 1 = 0 ,             f 1 S 2 = 0 ,             f 1 I = - ( β 1 S ) ,             f 1 R = 0 ,
(80)
f 2 S = 1 - ψ ρ - ε β 2 I ,             f 2 S 1 = - μ ,             f 2 S 2 = 0 ,             f 2 I = - ε β 2 S ,           f 2 R = 0 ,
(81)
f 3 S = ψ ρ ,             f 3 S 1 = 0 ,             f 3 S 2 = - μ ,             f 3 I = 0 ,             f 3 R = 0 ,
(82)
f 4 S = 0 ,             f 1 S 1 = 0 ,             f 1 S 2 = 0 ,             f 1 I = δ ,             f 1 R = - μ + υ .

Taking the absolute (module) values of Eqs. (79-82) show that fixj, where i= 1, 2, 3, 4 j=S,S1,S2,I,R are bounded. Hence, there is exists a unique solution of the model Eqs. (1-5).

5. Existence of equilibrium

Eqs. (1-5) are rewritten as follows:

(83)
d S d t = Λ - λ 1 S - A 1 S ,
(84)
d S 1 d t = 1 - ψ ρ S - μ S 1 - ε λ 2 S ,
(85)
d S 2 d t = ψ ρ S - μ S 2 ,
(86)
i ( t , τ ) t + i ( t , τ ) τ + A 2 i ( t , τ ) = 0 ,
(87)
d R d t = 0 T δ i ( t , τ ) d τ - A 3 R ,
(88)
λ 1 = 0 T β 1 i ( t , τ ) d τ ,             λ 2 = 0 T β 2 i ( t , τ ) d τ ,
(89)
B.C .    i t , 0 = B t = β 1 S t 0 T i t , τ d τ + ε β 2 S t 0 T i t , τ d τ ,
A 1 = μ + ρ ,             A 2 = μ + α τ + δ τ ,                 A 3 = μ + υ , .

At DFE, we let:

(90)
S t = x 0 ,             S 1 0 = y 0 ,           S 2 0 = u 0 ,             R 0 = r 0 ,             i t , 0 = z 0 ,             N t = n 0 .

And set dSdt+dS1dt+dS2dt+dIdt+dRdt=0, then we have the following results as the Disease Free Equilibrium (DFE):

(91)
E 0 = S 0 , S 1 0 , S 2 0 , I 0 , R 0 = Λ μ + ρ , 1 - ψ Λ ρ μ μ + ρ , ψ Λ ρ μ μ + ρ , 0.0 .

6. Effective reproduction number

In an infection age and age-structured models, the basic or effective reproduction number are often expressed as the sum of the infectivity of each infected compartment, see for example in [8, 9]. For a single infected compartment, the basic reproduction number is simply the product of the infection rate and the mean duration of the infection, the effective reproduction number of the model becomes:

(92)
R E = β 1 Λ μ + ρ 0 T π ( τ ) d τ + ε β 2 Λ μ + ρ 0 T π ( τ ) d τ .

This is the number of secondary cases the number of generated by individuals in the actively infected class and x0=(Λ/μ+ρ) represents the number of susceptible individuals in the absence of HIV/AIDS. The term π(τ)=e-0τ(μ+α(τ)+δ(τ))ds is the survival probability as a function of infection age τ in the actively infected class. The effective reproduction number obtained in Eq. (92), using this to determine the Local Stability of the Disease Free Equilibrium State.

7. Local stability of the disease free equilibrium (DFE) state

Theorem 5. The DFE state E0 is locally asymptotically stable if RE<1 and unstable if RE>1.

Proof: to show the local stability of DFE state given by Eq. (91), we have the following results:

(93)
z ¯ 0 = β 1 x 0 z ¯ 0 0 T e - λ τ π τ d τ + ε β 2 x 0 z ¯ 0 0 T e - λ τ π τ d τ .

Dividing both sides of Eq. (93) by z¯(0), we obtain:

(94)
1 = β 1 x 0 0 T e - λ τ π ( τ ) d τ + ε β 2 x 0 0 T e - λ τ π ( τ ) d τ ,

where x0 is defined as in Eq. (91).

Define a function G(λ) to be right hand side in Eq. (94). G(λ) is continuously differentiable function with limnG(λ)=0 this shown that G1(λ)<0 and therefore, G(λ) is a decreasing function. Hence, any real solution of Eq. (94) is negative if G(0)<1 and positive if G(0)>1. Thus, if G(0)>1, then DFE state is unstable.

To show that Eq. (94) has no complex solution with non-negative real part if G(0)<1. In fact, we set:

(95)
H τ = β 1 x 0 π τ ,             F τ = β 2 ε x 0 π τ .

Thus, we have:

(96)
G λ = 0 T e - λ τ H τ d τ + 0 T e - λ τ F τ d τ .

Suppose G(0)<0. Assume that λ=a1+b1i is a complex solution of Eq. (96) with a10 [7]:

(97)
G ( λ ) = 0 T e - λ τ H ( τ ) d τ + 0 T e - λ τ F ( τ ) d τ 0 T e - ( a 1 + i b 1 ) τ H ( τ ) d τ + 0 T e - ( a 1 + i b 1 ) τ F ( τ ) d τ
            0 T e - a 1 τ H τ d τ + 0 T e - a 1 τ F τ d τ = G ( a 1 ) G ( 0 ) < 1 .

It follows from Eq. (97) that Eq. (94) has solutions λ=a1+b1i only if a1<0. Thus, every solution of Eq. (94) must have a negative real part. We observed that RE=G(0). Therefore, the DFE state E0 is locally asymptotically stable if G(0)<1.

8. Global stability of the disease free equilibrium (DFE) state

We construct a Lyapunov function as outline in the work of [10] to study the global stability of the DFE state in this study.

Theorem 6. The DFE state E0 of the model Eqs. (1-5) is Globally Asymptotically Stable (GAS) in if RE1 while unstable if RE1.

Proof: Using the approach of [4], by constructing a suitable Lyapunov function as follows:

(98)
L ( t ) = x - x 0 - x 0 l n x x 0 + y - y 0 - y 0 l n y y 0 + u - u 0 - u 0 l n u u 0
            + 0 T r ( τ ) i ( t , τ ) π ( τ ) d τ .

Using Eq. (66) in Eq. (98) and at the DFE, we have the following results:

(99)
d L d t = - μ + ρ x ( x - x 0 ) 2 - μ y ( y - y 0 ) 2 - μ u ( u - u 0 ) 2 + i t , 0 R E - 1 ,
d L d t = - μ + ρ x ( x - x 0 ) 2 - μ y ( y - y 0 ) 2 - μ u ( u - u 0 ) 2 - i t , 0 1 - R E .

The equality dLdt=0 holds if and only if x=x0, y=y0, u=u0, i(t,0)=0. Thus, from the solution Eq. (29) for the model Eqs. (4) and (8) along the characteristics lines, it can be see that i(t,τ)=0 for all t>τ. Hence, i(t,τ)0 as t. Therefore, from the LaSalle invariant Principle, the Disease-Free Equilibrium (DFE) is globally asymptotically stable (GAS) if RE1.

9. Endemic equilibrium (EE) state

The endemic equilibrium state is referred to as the state in which HIV/AIDS still persist in the population. We let:

(100)
S = x * ,             S 1 ( t ) = y * ,             S 2 ( t ) = u * ,             R ( t ) = r * ,             i ( t , τ ) = i * ( τ ) ,
I t = z * ,             N t = n * ,             z * = 0 T i * τ d τ .

Let {x*,y*,u*,r*,i*(τ)} represents any arbitrary endemic equilibrium of the model Eq. (13). Thus, equilibrium satisfies the following equations:

(101)
Λ - β 1 x * z * - ( μ + ρ ) x * = 0 ,
(102)
( 1 - ψ ) ρ x * - μ y * - ε β 2 x * z * = 0 ,
(103)
ψ ρ x * - μ u * = 0 ,
(104)
i * ( τ ) τ + ( μ + α ( τ ) + δ ( τ ) ) i * ( τ ) = 0 ,
(105)
δ 0 T i * ( τ ) d τ - ( μ + υ ) r * = δ z * - ( μ + υ ) r * = 0 ,
(106)
i * 0 = β 1 x * z * + ε β 2 x * z * .

Solving Eq. (104) yields:

(107)
i * τ = C e - μ + α τ + δ τ τ .

From Eq. (107) we obtain:

(108)
i * τ = i * 0 e - μ + α τ + δ τ τ ,           i * 0 = C .

Eq. (108) is rewritten as:

(109)
i * τ = i * 0 π τ ,

where, π(τ) is defined as in Eq. (22).

Substituting Eq. (109) into Eq. (105) we get:

(110)
r * = δ ( μ + υ ) 0 T i * ( 0 ) π ( τ ) d τ .

Substituting Eq. (109) into Eq. (106), we have:

(111)
i * 0 = β 1 x * 0 T i * 0 π τ d τ + ε β 2 x * 0 T i * 0 π τ d τ .

Divide both side of Eq. (111) by i*(0), we get:

(112)
1 = β 1 x * 0 T π ( τ ) d τ + ε β 2 x * 0 T π ( τ ) d τ ,             x * = 1 β 1 π ¯ + ε β 2 π ¯ ,

where, π¯ is defined as in Eq. (26).

Substituting Eqs. (112) and (109) into Eq. (101) gives:

(113)
Λ - β 1 x * 0 T i * 0 π τ d τ - ( μ + ρ ) x * = 0 ,
Λ - β 1 x * i * ( 0 ) π ¯ - ( μ + ρ ) x * = 0 ,
i * 0 = Λ - μ + ρ x * β 1 x * π ¯ = Λ - μ + ρ 1 β 1 π ¯ + ε β 2 π ¯ β 1 1 β 1 π ¯ + ε β 2 π ¯ π ¯ ,
i * 0 = Λ β 1 π ¯ + ε β 2 π ¯ - μ + ρ β 1 π ¯ .

Substituting Eq. (112) into Eq. (110) gives:

(114)
r * = δ ( μ + υ ) 0 T Λ ( β 1 π ¯ + ε β 2 π ¯ ) - ( μ + ρ ) β 1 π ¯ π ( τ ) d τ ,
r * = δ ( μ + υ ) * Λ β 1 π ¯ + ε β 2 π ¯ - μ + ρ β 1 π ¯ π ¯ ,
(115)
r * = δ Λ ( β 1 π ¯ + ε β 2 π ¯ ) - δ ( μ + ρ ) β 1 ( μ + υ ) .

Substituting Eq. (112) into Eq. (103), we have:

(116)
u * = ψ ρ μ x * = ψ ρ μ ( β 1 π ¯ + ε β 2 π ¯ ) .

Substituting Eqs. (112), (109) into Eq. (102), we have:

1 - ψ ρ 1 β 1 π ¯ + ε β 2 π ¯ - μ y * - ε β 2 1 β 1 π ¯ + ε β 2 π ¯ 0 T i * 0 π τ = 0 ,
( 1 - ψ ) ρ β 1 π ¯ + ε β 2 π ¯ - μ y * - ε β 2 β 1 π ¯ + ε β 2 π ¯ * Λ β 1 π ¯ + ε β 2 π ¯ - μ + ρ β 1 = 0 ,
( 1 - ψ ) ρ β 1 π ¯ + ε β 2 π ¯ - μ y * - ε β 2 Λ β 1 π ¯ + ε β 2 π ¯ - μ + ρ β 1 β 1 π ¯ + ε β 2 π ¯ = 0 ,
(117)
y * = ( 1 - ψ ) ρ β 1 + ε β 2 ( μ + ρ ) - ε β 2 Λ ( β 1 π ¯ + ε β 2 π ¯ ) μ β 1 ( β 1 π ¯ + ε β 2 π ¯ ) .

Substituting Eq. (113) into Eq. (109) is setting as:

(118)
i * τ = Λ β 1 π ¯ + ε β 2 π ¯ - μ + ρ β 1 π ¯ * π ,
(119)
E * = ( S * , S 1 * , S 2 * , i * , R * )
            = 1 β 1 π ¯ + ε β 2 π ¯ , ( 1 - ψ ) ρ β 1 + ε β 2 ( μ + ρ ) - ε β 2 Λ ( β 1 π ¯ + ε β 2 π ¯ ) μ β 1 ( β 1 π ¯ + ε β 2 π ¯ ) ,             ψ ρ μ ( β 1 π ¯ + ε β 2 π ¯ ) , [ Λ ( β 1 π ¯ + ε β 2 π ¯ ) - ( μ + ρ ) ] * π β 1 π ¯ , δ Λ ( β 1 π ¯ + ε β 2 π ¯ ) - δ ( μ + ρ ) β 1 ( μ + υ ) .

10. Local stability of the endemic (EE) state

Theorem 7. The unique EE E*{x*,y*,u*,r*,i*(τ)} given by Eqs. (108), (115-118) with i*(0) given by Eq. (109) is locally asymptotically stable if RE> 1.

Proof: to show the local stability, we linearize the system Eq. (13) around the endemic equilibrium E* and we have the following results:

(120)
z ( 0 ) = i * ( 0 ) + ( β 1 x * + ε β 2 x * ) 0 t z ( 0 ) e - λ t π ( τ ) d τ
          + β 1 z * + ε β 2 z * - β 1 x * 0 t z 0 e - λ t π τ d τ q 1 + β 1 z * + β 1 0 t z 0 e - λ t π τ d τ
          + β 1 + ε β 2 - β 1 x * 0 t z 0 e - λ t π τ d τ 2 q 1 + β 1 z * + β 1 0 t z 0 e - λ t π τ d τ .

Divide through by z0, then we have:

(121)
1 = Q ( 0 ) + β 1 x * + ε β 2 x * + ( β 1 z * + ε β 2 z * ) ( - β 1 x * ) ( ( μ + ρ + λ ) + β 1 z * + β 1 0 t z ( 0 ) e - λ t π ( τ ) d τ )
      0 t e - λ t π τ d τ + β 1 + ε β 2 - β 1 x * μ + ρ + λ + β 1 z * + β 1 0 t z 0 e - λ t π τ d τ 0 t e - λ t π τ d τ 2 ,

where:

(122)
Q 0 = i * 0 z 0 .

Suppose that λ is real, then from the characteristics Eqs. (122) and (97) and the fact that RE>1, we obtained that λ<0, since for λ0, then right – hand side of Eq. (122) is non-negative while, the left-hand side is positive. Therefore, the endemic equilibrium E* is locally asymptotically stable, since the real root has the dominant real part.

11. Global stability of endemic equilibrium (EE) state

Theorem 8. If RE>1 then there exists a unique endemic equilibrium E*{x*,y*,u*,r*,i*(τ)} where x*, y*, u*, r* and i*(τ) satisfies the Eqs. (101-106).

Define L2(a,b)=a-b-blnab for a, b>0. We construct the following Lyapunov functional:

(123)
W 2 = W 3 + W 4 + W 5 + W 6 + W 7 .

With:

(124)
W 3 = x - x * - x * l n x x * ,               W 4 = y - y * - y * l n y y * ,               W 5 = u - u * - u * l n u u * , W 6 = r - r * - r * l n r r * ,                 W 7 = 0 T i τ L 2 i t , τ , i * τ d τ .

The following Lemmas evaluate the derivatives W3, W4, W5, W6, W7 along the solutions of Eqs. (1-5) respectively.

Lemma 2. Let W3 be defined as in Eq. (124), then:

W 3 1 = - ( μ + ρ ) ( x - x * ) 2 x - β 1 x - x * x x z - x * z * .

Proof: Direct differentiation of Eq. (124) is given by:

(125)
W 3 1 = x - x * x d x d t .

Using Eq. (1) in Eq. (125), we have:

(126)
W 3 1 = x - x * x Λ - μ x - λ 1 x - ψ ρ x - 1 - ψ ρ x .

From Eq. (101), we obtain:

(127)
Λ = β 1 x * z * + μ + ρ x * .

Substituting Eq. (127) into Eq. (126) we have:

(128)
W 3 1 = - ( μ + ρ ) ( x - x * ) 2 x + β 1 x - x * x x * z * - x z .

Lemma 3. Let W4 be defined as in Eq. (124), then:

W 4 1 = x y ε β 2 z * - z y - y * - μ 1 - y * x y x * y - y * .

Proof:

(129)
W 4 1 = y - y * y d y d t = y - y * y 1 - ψ ρ x - μ y - ε β 2 x z .

From Eq. (102), we obtain:

(130)
1 - ψ ρ = μ y * x * + ε β 2 z * ,
(131)
W 4 1 = - μ 1 - y * x y x * y - y * + x y ε β 2 z * - z y - y * .

Lemma 4. Let W5 be defined as in Eq. (124), then:

W 5 1 = ψ ρ y u u - u * - ψ ρ x * u * u - u * .

Proof: direct differentiation of Eq. (124) is given by:

(132)
W 5 1 = u - u * u d u d t = u - u * u ψ ρ y - ψ ρ x * u * u .
(133)
W 5 1 = ψ ρ y u u - u * - ψ ρ x * u * u - u * .

Lemma 5. Let W6 be defined as in Eq. (124), then:

(134)
W 6 = δ z r r - r * - δ z * r * r - r * .

Proof: Direct differentiation of Eq. (124) is given by:

(135)
W 6 = r - r * r d r d t = r - r * r δ z - μ + υ r .

From Eq. (105), we obtain:

(136)
W 6 = δ z r r - r * - δ z * r * r - r * .

We state the following lemma by [11] (Lemma 5.2) to aid us establish the proof of Lemma 6.

Lemma 6.

(137)
W * = 0 ϑ ( a ) L 2 ( i ( t , a ) , i * ( a ) ) d a .

then:

(138)
W * 1 = ϑ 0 L 2 i t , 0 , i * 0 + 0 ϑ 1 a - θ a ϑ a L 2 i t , a , i * a d a .

Lemma 7. Let W7 be defined as in Eq. (124), then:

W 7 1 = i * 0 1 - i t , 0 i * 0 + l n i t , 0 i * 0 + 2 0 T d i * τ d t 1 - i t , 0 i * 0 + l n i t , 0 i * 0 d τ .

Proof: applying the result in Lemma 6, we obtain the direct differentiation of Eq. (124) as:

(139)
W 7 = i * 0 L 2 i t , 0 , i * 0 + 0 T d i * τ d t - μ + α τ + δ τ i * τ L 2 i t , τ , i * τ d τ .

Using L2(a,b)=a-b-blnab for a, b>0 in Eq. (139) is given by:

W 7 1 = i * 0 i t , 0 - i * 0 - i * 0 l n i t , 0 i * 0
            + 0 T d i * τ d t - Φ τ i * τ i t , 0 - i * 0 - i * 0 l n i t , 0 i * 0 d τ ,

or:

(140)
W 7 1 = i * 0 1 - i t , 0 i * 0 + l n i t , 0 i * 0
            + 0 T d i * τ d t - Φ τ i * τ 1 - i t , 0 i * 0 + l n i t , 0 i * 0 d τ ,

where, Φ(τ)=μ+α(τ)+δ(τ).

Using Eq. (104) i.e (μ+α(τ)+δ(τ))=i*(τ)i*(τ)τ in Eq. (140), we have:

(141)
W 7 1 = i * 0 1 - i t , 0 i * 0 + l n i t , 0 i * 0 + 2 0 T d i * τ d t 1 - i t , 0 i * 0 + l n i t , 0 i * 0 d τ .

The addition of Eqs. (128), (131), (133), (136) and (141) take the for:

(142)
d W 2 ( t ) d t = d ( W 3 + W 4 + W 5 + W 6 + W 7 ) d t = X 1 - X 2 + X 3 ,
(143)
d W 2 ( t ) d t = X 1 - X 2 + X 3 ,

where:

(144)
X 1 = β 1 x - x * x x * z * - x z + x y ε β 2 z * - z y - y * + ψ ρ y u u - u * + δ z r r - r * ,
(145)
X 2 = ( μ + ρ ) ( x - x * ) 2 x + μ 1 - y * x y x * y - y * + ψ ρ x * u * u - u * + δ z * r * r - r * ,
(146)
X 3 = i * 0 1 - i t , 0 i * 0 + l n i t , 0 i * 0 + 2 0 T d i * τ d t 1 - i t , 0 i * 0 + l n i t , 0 i * 0 d τ ,

Since 1-a+lna0, a>0 implies from Eq. (146) that z0 with equality holding if and only if a=1. From Eq. (143), if X2>X1 then dW2dt will be negative definite, meaning that dW2dt<0. Also it follows that dW2dt=0 if and only if x=x*, y=y*, u=u*, r=r* and i(t,0)i*(0)=i(t,τ)i*(τ), τ0. Therefore, the largest compact invariant set in {x*,y*,u*,r*,i*(τ)Ω:dW2dt=0} is the singleton {E*} if RE>1, then by the Lassalle’s invariant principle, {E*} is globally asymptotically stable in Ω if X1<X2.

12. Conclusions

In this paper, an infection-age-structured mathematical model to study the effect of HIV public awareness program, with unaware HIV susceptible on HIV/AIDS dynamics was developed and analysed. We took into consideration the public awareness, Treatment and Sex behaviours. The criteria for the existence and uniqueness of solution were obtained. it was discovered that the model has two equilibrium states; Disease Free Equilibrium (DFE) and Endemic Equilibrium (EE). We analysed the invariant region and shown that the dynamics of model Eqs. (1-9) is in the region , the model equations was considered as been epidemiologically and mathematically well posed. An explicit formula for the effective reproduction number RE of the model is obtained. The results from stability analysis showed that DFE is locally asymptotically stable when RE<1 and globally asymptotically stable when RE<1. More so, the EE was found to be locally asymptotically stable when RE>1 and globally asymptotically stable when RE>1 for X1<X2.

References

  1. Udoo I. J. M. Modelling impact of public awareness and behaviour changes with treatment on HIV/AIDS dynamics in resource-poor settings. Journal of the Nigerian Society for Mathematics Biology (JNSMB), Vol. 1, 2018, p. 19-46. [Search CrossRef]
  2. Ashezua T. T. An Infection-age-structured Mathematical Model for Tuberculosis Disease Dynamics Incorporating Control Measures. Ph.D. Thesis, Federal University of Technology, Minna, Nigeria, 2016. [Publisher]
  3. Akinwande N. I. A mathematical model of the chaotic dynamics of the AIDS disease pandemic. Journal of Nigeria Mathematical Society (NMS), Vol. 24, Issue 1, 2005, p. 8-16. [Search CrossRef]
  4. Lasisi N. O., Akinwande N. I., Olayiwola R. O., et al. Mathematical model for Ebola virus infection in human with effectiveness of drug usage. Journal of Applied Sciences and Environmental Management, Vol. 22, Issue 7, 2018, p. 1089-1095. [Publisher]
  5. Birkhoff G., Rota G. C. Ordinary differential equations. 2nd edition, Oxford University Press, 1982. [Search CrossRef]
  6. Sebastian Anita Analysis and Control of Age-Dependent Population Dynamics. Kluwer Academic Publishers, Dordrecht/Boston/London, 2000. [Search CrossRef]
  7. Webb Glenn F. Theory of Nonlinear Age-Dependent Population Dynamics. A Series of Monographs, Textbook and Lecture Note in Pure and Applied Mathematics. Marcel Dekker, New York, 1985. [Search CrossRef]
  8. Shuai Z., Van Den Driessche P. Global stability of infectious disease models using Lyapunov functions. SIAM Journal on Applied Mathematics, Vol. 73, 2013, p. 1513-1532. [Publisher]
  9. Wang J., Zhang T. Mathematical analysis for an age-structured HIV infection model with saturation infection rate. Electronic Journal of Differential Equations, Vol. 33, 2016, p. 1-19. [Search CrossRef]
  10. Huang G., Liu X., Takeuchi Y. Lyapunov functions and global stability for age structured HIV infection model. Journal Applied Mathematics, Vol. 72, Issue 1, 2012, p. 25-38. [Search CrossRef]
  11. Brauer F., Shuai Z., Van Den Driessche P. Dynamics of an age of infection cholera model. Mathematics Biosciences and Engineering, Vol. 10, Issue 6, 2013, p. 1335-1349. [Search CrossRef]