arxiv: v1 [math.ds] 14 Dec 2017
|
|
- Lorin Gilbert
- 5 years ago
- Views:
Transcription
1 Global properties of an age-structured virus model with saturated antibody immune response, multi-target cells and general incidence rate Ángel Cervantes-Pérez, Eric Ávila-Vales arxiv: v1 [math.ds] 14 Dec 217 Facultad de Matemáticas, Universidad Autónoma de Yucatán, Anillo Periférico Norte, Tablaje 13615, C.P , Mérida, Yucatán, Mexico Abstract Some viruses, such as human immunodeficiency virus, can infect several types of cell populations. The age of infection can also affect the dynamics of infected cells and production of viral particles. In this work, we study a virus model with infectionage and different types of target cells which takes into account the saturation effect in antibody immune response and a general non-linear infection rate. We construct suitable Lyapunov functionals to show that the global dynamics of the model is completely determined by two critical values: the basic reproduction number of virus and the reproductive number of antibody response. 1 Introduction In recent times, several mathematical models have been proposed in order to try to understand the mechanism of virus infections. These models often describe the changes through time in the concentration of infected and uninfected target cells and viral particles in the blood of an infected individual [1 16]. Modelling the effect of antibody immune response in the neutralization of virus is a very important topic for research since it can provide useful insights into the dynamics of the infection and offer suggestions for clinical treatment. The production of new virions by an infected cell does not occur at a constant rate, and the death rate of infected cells can vary with the time the cell has been infected. Thus, the incorporation of an age structure in models allows us to have a more realistic picture of produced viral particles and the mortality of infected cells; see, for example, the models studied in [1 6]. A recent work by Duan and Yuan [3] introduced a model with age-structure and antibody immune response, described by the system iθ, t θ dt t βt tv t λ dt t 1 + αv t iθ, t + δθiθ, t t dv t P θiθ, t cv t qatv t dat katv t h + At bat, 1
2 with the initial and boundary conditions i, t βt tv t 1 + αv t, T T, iθ, i θ, V V, A A. In this model, T, i, V, A denote the concentration of uninfected susceptible host cells, infected host cells, free virus particles, and antibody responses released from B cells, respectively. The variable θ is the infection age, i.e., the time that has elapsed since a virion has infected a cell, while iθ, t is the density of infected T cells with infection age θ at time t. The functions P θ and δθ are respectively the time-since-infection structured virion production rate and the death rate of productively infected cells. The authors in [3] showed that the global dynamics of this model is completely determined by the basic reproduction number R and the viral reproductive number R 1. Most mathematical models for human immunodeficiency virus HIV focus on the infection of CD4 + T cells. However, macrophages and dendritic cells are also susceptible to be infected with HIV [4, 5]. Therefore, we need to incorporate in our virus dynamic models a multi-group component to study virus infection in different populations of cells. Viral infection models dealing with the interaction of virus with more than one class of target cells have been studied in [7 13] and references cited therein. When modelling the rate of infection of susceptible target cells by contact with virus, the bilinear incidence term βt V is used frequently in literature [1, 2, 14]. However, there is certain debate that this rate is insufficient to realistically describe the infection process. Some authors have incorporated in their models several types of non-linear incidence functions, including saturated Holling type II incidence βt V/1 + αv [3, 5, 15], Beddington-DeAngelis functional response βt V/1+aT +bv [6, 7], or even more general functions, e.g., of the forms ct fv [16] or ht, V [4, 8, 9]. In this paper, we propose a viral infection model that includes age-structure, multitarget cells, and a general non-linear rate of viral infection. We establish the global properties of its equilibria using semigroup methods, uniform persistence and the construction of Lyapunov functionals. The rest of this paper is organized as follows: in Section 2 we present the equations and assumptions for our model, we calculate the basic reproduction number R and show boundedness and non-negativity of solutions. In Section 3 we rewrite the system as a semilinear Cauchy problem, we prove that the semiflow generated by the system has a global compact attractor and that it is uniformly persistent when R > 1. In Section 4 we determine the steady states of the system and deal with their local stability. In Section 5 we construct several Lyapunov functionals to determine the global stability of equilibria. In Section 6 we present a special case in which the model can be reduced to a system of delay differential equations. In Section 7 we provide some numerical simulations which illustrate the different kinds of global behaviour the model can have. Finally, in Section 8 we present some concluding remarks. 2 Description of the model We consider a within-host viral dynamics model for HIV with multiple target cell populations, which takes into account the saturation effects of antibody immune response and incorporates a general incidence function h j T j, V. The model is given by the following 2
3 system of differential equations: i j θ, t θ with initial conditions dt j t + i jθ, t t dv t dat λ j T j t h j T j t, V t δ j θi j θ, t P j θi j θ, t cv t qatv t katv t h + At bat 2.1 and boundary conditions T j T j, i j θ, i jθ, V V, A A 2.2 i j, t h j T j t, V t. 2.3 This model considers n classes of target cells, denoted by the subscript j 1,..., n. For each class, T j represents the population of uninfected target cells, while i j θ, t denotes the population of infected cells with infection age θ at time t. We denote by V the concentration of free virus particles, and by A the number of antibody responses released from B cells. We assume that uninfected cells in the j-th class are recruited at a rate λ j and die at a rate. Uninfected cells become infected from contact with the virus according to the incidence function h j T j t, V t. Infected cells in the j-th class die at a rate δ j θ, which depends on the infection age θ. The function P j θ represents the viral production rate of an infected cell from the j-th class with infection age θ, while the parameter c is the viral clearance rate. The virus is removed by the immune system at a rate qa. The antibody responses are generated by contact with the virus at a rate katv t/h + At and vanish at rate b. All parameters are assumed positive. Furthermore, we assume that the functions δ j and P j satisfy δ j, P j L +,, R \ { L } and δ j θ δ min > for all θ, j 1,..., n. Let P max max {ess sup P j j 1,..., n}, where ess sup P j denotes the essential supremum of P j θ for θ,. We make the following hypotheses on the incidence functions h j, j 1,..., n: H1 h j T j, V V h j T j, V with h j, h j C 2 R 2 + R + and hj, V h j T j, for all T j, V. H2 h j T j, V is strictly increasing with respect to both variables when T j, V >. H3 h j T j, V for all T j, V and h j T j, V as V. 2.1 Basic reproduction number and boundedness of solutions For each infected cell in the j-th class, the probability of still being infected after θ time units is given by σ j θ e θ δ jτ dτ, 3
4 so the total number of virions produced by an infected cell from the j-th class in its entire life span is N j P j θσ j θ. Hence the corresponding basic reproduction number for system 2.1 when the j-th class is the unique class of target cells is given by R j N j c h j λj,, j 1,..., n. Therefore, we can define the basic reproduction number R of virus for the model as R R j 1 c h j λj N j,, 2.4 which represents the expected number of viruses that one virion produces in a fully uninfected cell population. The second equation of 2.1 is a linear transport equation with decay, so it can be solved with the initial conditions 2.2 and boundary conditions 2.3 by the method of integration along the characteristic lines t θ constant. The solution is given by { hj Tj t θ, V t θ σ j θ if θ < t, i j θ, t i j θ t σ jθ σ j θ t if θ t, 2.5 where σ j θ e θ δ jτ dτ. We will now prove that the model is biologically well-posed by showing the nonnegativity and boundedness properties of the system. Theorem 2.1. The solutions of system 2.1 with non-negative initial conditions 2.2 and boundary conditions 2.3 remain non-negative for t and are positively bounded. Proof. Let T 1, i 1,..., T n, i n, V, A be a solution of 2.1 with non-negative initial conditions. Suppose that T j t loses its non-negativity for some j and let t inf{t T j t < }. By continuity of solutions, we have T j t and dt j t. But the first equation of 2.1 implies that dt j t λ j >, a contradiction. Therefore, no such t exists and thus, T j t > for all t. Next, we show that V t > and i j θ, t for all t, θ. Define t 1 min { inf{t V t < }, inf{t i 1, t / L 1 +, },..., inf{t i n, t / L 1 +, } }. Suppose first that t 1 inf{t i j, t / L 1 +, } for some j. Then 2.5 implies that i j θ, t 1 for all θ, which is a contradiction. Hence, t 1 inf{t V t < }. At the time instant t 1, we have dv t 1. However, since i j, t 1 > for all j, the third equation of 2.1 implies that dv t 1 P j θi j θ, t 1 >, which leads to a contradiction. Thus, V t > and i j θ, t for all t, θ. Finally, from the last equation of the model we have At Ae t kv τ h+aτ b dτ. 4
5 Therefore, the solution remains non-negative on its maximal interval of existence. Next, we consider the boundedness of solutions. From the first equation of 2.1, we have dt j λ j T j for j 1,..., n. This implies that Let Gt lim sup T j t λ j. t T j t + i j θ, t be the total number of cells uninfected and infected at time t. Then dgt dtj t i j θ, t + t [ λ j T j t h j T j t, V t + [ λ j T j t λ j ] δ j θi j θ, t [ T j t + δ min ] i j θ, t, since δ j θ δ min. Let d min{d 1,..., d n, δ min }. Then and thus, dgt λ j d lim sup t n where M λ j /d. so Now, from the equation for dv dv t when t is large enough. Thus, [ T j t + i jθ, t δ j θi j θ, t θ ] i j θ, t λ j dgt, i j θ, t lim sup Gt M, t in 2.1, we have P max dv t lim sup t P j θi j θ, t cv t i j θ, t cv t, P max M cv t V t P maxm. c ] 5
6 Finally, from the last equation in 2.1, we have for large t, so dat kv t bat kp maxm c lim sup t which completes the proof of the theorem. At kp maxm, bc bat Let Ω { T1, i 1,..., T n, i n, V, A T j λ j, T j + n i j θ [ R + L 1 +,, R ] R 2 + } M, V P maxm, A kp maxm. c bc Then Ω is a positively invariant set for system 2.1. From now on, we will always assume that the initial value T 1, i 1,..., T n, i n, V, A is in Ω. 3 Integrated semigroup formulation We will now reformulate system 2.1 with the initial condition 2.2 as a semilinear Cauchy problem. In order to take into account also the boundary condition 2.3, we need to expand the state space as follows. Let M R R L 1,, R, P R {} L 1,, R, N R {} W 1,1,, R n, X M R R, where W 1,1 denotes a Sobolev space. Consider the linear operator A : DomA X X given by T 1 d 1 T 1 i 1 i 1 i 1 δ 1θi 1.. A T n d n T n, i n i n i n δ n θi n V cv A ba with DomA n i1 P R R. Then n DomA N R R : X i1 i1 6
7 and we can consider the non-linear map F : X X given by T 1 λ 1 h 1 T 1 t, V t h1 T 1 t, V t i 1 L1.. F T n λ n h n T n t, V t. hn T n t, V t i n L1 V n P j θi j θ, t qatv t A Let katv t h+at T ut T 1 t,,..., T i 1, t n t,, V t, At, i n, t P + R + {} L,, R n, X + P + R + R +. Thus, we can rewrite system 2.1 with the boundary and initial conditions as the following abstract Cauchy problem: i1 dut Aut + F ut for t, with u x X By applying the results given in [17] and [18], we can conclude the following theorem. Theorem 3.1. System 2.1 generates a unique continuous semiflow {Ut} t on X + that is asymptotically smooth and bounded dissipative. Furthermore, the semiflow {Ut} t has a global compact attractor in X +, which attracts the bounded sets of X Uniform persistence We will now establish the uniform persistence of system 2.1, using arguments that are highly motivated by those in [6] and [3]. Define M T1,, i 1,..., T n,, i n, V, A X + V + i j θ >, and M X + \ M. For j 1,..., n, let ϑ j sup{θ, P j θ > }. Note that ϑ j can possibly be +. Lemma 3.2. M is positively invariant under the semiflow {Ut} t generated by system 2.1. Moreover, the infection-free equilibrium E λ 1 λ d 1,, 2 λ d 2,,..., n d n,,, is globally asymptotically stable restricted to M. Proof. Let y M and T 1 t,, i 1, t,..., T n t,, i n, t, V t, At Uty. Then we have i j θ, t θ + i jθ, t t δ j θi j θ, t 7
8 dv t i j, t h j T j t, V t, V V. P j θi j θ, t cv t qatv t i j θ, i jθ, Since T j t λ j for large t, then i j θ, t ĩ j θ, t and V t solution of the system Ṽ t, where ĩ, Ṽ is the ĩ j θ, t θ + ĩ j θ, t δ j θĩ j θ, t t dṽ t P j θĩ j θ, t cṽ t ĩ j, t h j λj, Ṽ t, Ṽ V. ĩ j θ, i jθ 3.2 From the first and third equation of 3.2, we have { λj hj d ĩ j θ, t j, Ṽ t θ σ j θ if θ < t, i j θ t σ jθ σ j θ t if θ t and by substituting this expression in the second equation of 3.2, we obtain where dṽ t t h j λj, Ṽ t θ P j θσ j θ + F j t t P j θi σ j θ jθ t σ j θ t. 3.3 F j t cṽ t, 3.4 For each j, we will prove that F j t for all t. In fact, if t ϑ j, then P j θ for θ > t ϑ j, so F j t t P j θi j θ t. If t < ϑ j, then F j t t P j θi jθ t ϑj P max i jθ t + t ϑ j P j θi jθ t Since y M, then ϑ j t i j θ t V + n i j θ, and also P jθ for θ > ϑ j, so the two terms in the right-hand side of the aforementioned inequality are equal to zero. Thus, Fj t. Accordingly, 3.4 is an autonomous equation and has a unique solution Ṽ t. dat Hence V t Ṽ t and bat for large t, which implies that V t and At approach as t. Now, it follows from 3.3 that ĩ j θ, t for θ t. For t < θ, we have ĩ j θ, t L 1 σ j θ i jθ t σ j θ t e δ mint i j L 1. L 1 Therefore, i j θ, t as t. Since i j θ, t ĩ j θ, t, then Uty approaches E as t and the proof is complete. 8
9 Theorem 3.3. If R > 1, then the semiflow {Ut} t is uniformly persistent with respect to the pair M, M, i.e., there exists ɛ > such that for each y M, lim inf t d Uty, M ɛ. Furthermore, there exists a compact subset A M that is a global attractor for {Ut} t in M. Proof. Since E is globally asymptotically stable in M, by Theorem 4.2 in [19] we only need to investigate the behaviour of the solutions starting in M in some neighbourhood of E. Then, we will show that W s {E } M, where W s {E } { y X + lim t Uty E }. Assume by contradiction that there exists y W s {E } M. It follows that there exists t > such that V t + n i j θ, t >. Using the same argument as in the proof of Lemma 3.6i in [2], we have that V t > for t and i j θ, t > for any θ, t. Define the function γ j x x P j θe θ x δ jτ dτ. Note that γ j is bounded and satisfies γ j x δ jxγ j x P j x for all x. Consider the function which satisfies dφt Φt cv t γ j θi j θ, t + V t, N j h j T j t, V t cv t qat c 1. Since y W s {E }, we have that T j t λ j, V t, and At as t. Since R > 1, then lim t N j h j T j t, V t cv t qat c 1 N j c h j λj, 1, so Φt is non-decreasing for sufficiently large t. Thus, there exists t 1 > such that Φt Φt 1 for all t > t 1. Since Φt 1 >, this prevents the function V t, i j, t from converging to, L 1 as t, which contradicts that T j t λ j. This completes the proof. 4 Steady states and local stability In this section, we will study the existence and local properties of the steady states for system 2.1. In order to prove the existence of infected equilibria, we first need to prove the following results. 9
10 Lemma 4.1. For j 1,..., n, T j λ j and V, the equation λ j T j h j T j, V 4.1 has a unique solution given by T j f j V, where f j : R +, λ j ] is a decreasing function with f j λ j. Proof. Let ϕt j, V λ j T j h j T j, V. If V, it is clear that the only solution of 4.1 is T j λ j. If V >, then assumption H2 implies that ϕt j, V is strictly decreasing with respect to T j. Since ϕ, V λ j > and ϕ λ j, V < λ j λ j, then ϕt j, V has a solution T j f j V in the interval, λ j, and it is unique due to monotonicity of ϕ. Thus, the function f j is well defined and positive. Moreover, differentiation of 4.1 gives + h j T j T j, V dt j h j T j, V dv, so f jv dt j dv hj T j, V / + h j T j, V, T j which is negative for V > by assumption H2, Therefore, f j is decreasing. Lemma 4.2. The function ψ j : R + R + defined by ψ j V h j f j V, V is strictly decreasing for V, and ψ j h j λj,. Proof. Let V 1 < V 2. Then f j V 2 < f j V 1 by Lemma 4.1. By H2 this implies that h j f j V 2, V 2 < h j f j V 1, V 2, so h j f j V 2, V 2 < h j f j V 1, V 2. By H3 we also have that h j is non-increasing with respect to the second variable, so h j f j V 1, V 2 h j f j V 1, V 1. Therefore, by transitivity, h j f j V 2, V 2 < h j f j V 1, V 1, i.e., ψ j V 2 < ψ j V 1. For V, we have λj h j f j, V h j d ψ j h j f j, lim lim j, V λj h j, V V V V hence the result. h j λj,, We can now define the viral reproduction number R of model 2.1, which is given by R 1 c N j h j fj bh k, bh k 4.2 and can be interpreted as the average number of antibodies that are activated by the introduction of a single virion within the host, under the condition that there is no previous antibody response. Since bh k >, Lemma 4.2 implies that h j fj bh k, bh k < hj f j, h j λj,, and by using the expression in 2.4 for the basic reproduction number R we can see that R < R. These two numbers, R and R, determine whether system 2.1 has one, two, or three steady states, as asserted in the following theorem. 1
11 Theorem 4.3. System 2.1 always has an infection-free steady state E T 1,, T 2,,..., T n,,,, where T j λ j/. In addition, when R > 1, system 2.1 has a unique immune-free infected steady state E T 1, i 1 θ,..., T n, i nθ, V,, where T j f j V, i jθ h j fj V, V σ j θ, and V is the unique positive solution of N j h j fj V, V c. 4.3 Also, if R > 1, then system 2.1 has a unique antibody-immune infected steady state Ê ˆT1, î 1 θ,..., ˆTn, î n θ, ˆV, Â, where ˆT j f j ˆV, î j θ h j fj b k h + Â, b k h +  b σ j θ, ˆV h + Â, k and  is the unique positive solution of N j h j fj b k h + Â, b k h +  c qâ. 4.4 Proof. It is clear that for 2.1 the infection-free steady state E always exists. To obtain the immune-free infected steady state E, we have to assume that A and V. We obtain the system λ j T j h j T j, V di j θ δ j θi jθ P j θi jθ cv. 4.5 By Lemma 4.1, the solution to the first equation of this system is T j f jv, while the second equation together with the boundary condition 2.3 yields i j θ h jt j, V σ j θ h j f j V, V σ j θ. Substituting in the third equation of 4.5, we get P j θσ j θh j f j V, V cv and since V, we can divide the above equation by V and write it as F V with F V N j h j f j V, V c. By Lemma 4.2, the function V F V is continuous and strictly decreasing for V, and since lim V h jf j V, V by H3, then lim F V V N j lim h jf j V, V c c <. V 11
12 This implies that F V has no positive solutions if F, while it has exactly one positive solution if F >. Since h j f j, ψ j h j λj,, then F h j λj N j, c h j λj N j, 1 1, R so the steady state E exists and is positive if and only if R > 1. For the antibody-immune infected steady state Ê, we assume  and ˆV, so we obtain the system λ j ˆTj h j ˆT j, ˆV dî j θ δ j θî j θ P j θî j θ c ˆV qâ ˆV k ˆV h +  b. 4.6 The last equation of this system implies ˆV b k h+â. From the first and second equations, we have ˆT j f j ˆV and î j θ h j ˆT j, ˆV σ j θ h j fj b k h + Â, b k h +  σ j θ. Substituting in the third equation of 4.6, we get P j θσ j θh j fj b k h + Â, b k h +  c ˆV qâ ˆV and since ˆV, this is equivalent to GÂ, where GA N j h j fj b k h + A, b k h + A c qa. By Lemma 4.2 we know that h j f j V, V decreases with respect to V, so A GA is a continuous and decreasing function for A. Since V b k h + A tends to infinity as A and lim h jf j V, V, then V lim GA A N j lim h jf j V, V + lim c qa. V A This implies that GA has no positive solutions if G, while it has exactly one positive solution if G >. We have G N j h j fj bh k, bh k c N j h j fj bh k, bh k so the steady state Ê exists and is positive if and only if R > R, The above theorem implies that system 2.1 has only one steady state when R 1, two steady states when R 1 < R and three when 1 < R. We will now analyse their local stability by means of the characteristic equation of the system. 12
13 Theorem 4.4. The infection-free steady state E of system 2.1 is locally asymptotically stable if R < 1 and it is unstable if R > 1. Proof. Using the expression in 2.5 for i j θ, t, we have P j θi j θ, t t t P j θh j Tj t θ, V t θ σ j θ + k j θh j Tj t θ, V t θ + g j t, where k j θ P j θσ j θ and g j t t P j θi j θ t σ jθ σ j θ t 2.1 as the following system: t P j θi σ j θ jθ t σ j θ t. Thus, we can rewrite dt j t dv t dat λ j T j t h j T j t, V t [ t k j θh j Tj t θ, V t θ ] + g j t cv t qatv t katv t h + At bat. The linearization of this system at E is dt j t T j t h j λj, V t dv t hj λj, t k j θv t θ cv t dat bat, and the characteristic equation is n s + hj λj s + c, k j θe sθ s + b. This equation has the negative roots s b and s, j 1,..., n. The rest of its roots are the solutions of hj λj φs : s + c, k j θe sθ. 4.7 Suppose that 4.7 has a root s with non-negative real part. Then c s + c hj λj, k j θe sθ hj λj, k j θ e s θ, 13
14 but so k j θ e s θ c k j θ hj λj, N j cr. P j θσ j θ N j, This implies that 1 R. Therefore, when R < 1 all solutions of 4.7 have negative real part and thus, in such case, E is locally asymptotically stable. Otherwise, if R > 1 we have φ c hj λj, k j θ c1 R <, while φs as s. Then φs has at least one positive real root. Therefore, E is unstable when R > 1. We will now study the stability of the immune-free infected steady state E. For that, we define the reproductive number of antibody response R AN as R AN kv bh. We will also need the following result, which shows the relation between R AN and the viral reproduction number R. Lemma 4.5. R < 1 R AN < 1, R 1 R AN 1, and R > 1 R AN > 1. Proof. Suppose that R AN < 1. Then V < bh k. Recall that V satisfies F V N j h j fj V, V c. Since F is strictly decreasing, the condition V < bh k that is, N j h j fj bh k, bh k c <, cr 1 <. implies that F bh k < F V, i.e., Thus, R < 1. The proof that R AN 1 R 1 and that R AN > 1 R > 1 is similar, from which the result follows. Theorem 4.6. The immune-free infected steady state E is unstable when R > 1. 14
15 Proof. The linearization of system 2.1 at E is dt j t + h j Tj, V T j t h j T j T j, V V t dv t t hj Tj, V k j θt j t θ T j t hj + T j, V k j θv t θ cv t dat kv h b At. Thus, one of the roots of the characteristic equation is s 1 kv h b. If R > 1, then R AN kv bh > 1, which implies that s 1 is a positive root. Hence the theorem. 5 Global stability In order to establish the global stability of the equilibrium E when R 1, we will make use of LaSalle s invariance principle and a Lyapunov functional, similar to those used in [4] for a multi-target cell model with general incidence. Theorem 5.1. If R 1, then the infection-free steady state E is globally asymptotically stable. Proof. Assume R 1. For j 1,..., n, define the function γ j x x P j θe θ x δ jτ dτ. 5.1 We will use the Lyapunov functional W W T 1, i 1,,, T n, i n,, V, A given by W Tj N j T j Tj h j Tj lim, V Tj V h j η j, V dη j + γ j θi j θ, t + V. Then W is non-negative and W E. The time derivative of W along the solutions of 2.1 is dw h j Tj N j 1 lim, V λ j T j h j T j, V V h j T j, V ij θ, t γ j θ + δ j θi j θ, t θ + P j θi j θ, t cv qav h j Tj N j 1 lim, V V h j T j, V T j T j h j Tj N j 1 lim, V h j T j, V V h j T j, V 15
16 cv qav. Since T j λ j Thus, we have dw γ j θ di j θ, t h j T 1 lim j,v V h j T j,v γj θδ j θ P j θ i j θ, t h j Tj, V for all V, so 1 lim,v V h j T j,v. Tj, then h jt j, V h j Tj Tj T j and also qav, so N j 1 lim V h j T j, V h j T j, V h j T j, V γj θδ j θ P j θ i j θ, t cv N j h j T j, V + γ j θi j θ, t θ θ h j Tj N j h j T j, V lim, V V h j T j, V i j θ, tγ jθ γj θδ j θ P θ i j θ, t cv, γ j θ di j θ, t where we have used integration by parts to expand γ j θ di j θ, t. By 5.1 we know that γ j satisfies γ j x δ jxγ j x P j x, so i j θ, tγ jθ γj θδ j θ P j θ i j θ, t, and since γ j N j, then N j h j T j, V γ j θi j θ, t θ. Thus, dw N j h j T j, V + γ j θi j θ, t h j Tj N j h j T j, V lim, V V h j T j, V cv θ θ h j Tj N j h j T j, V lim, V V h j T j, V γ j θi j θ, t cv θ cv N j h j T j, V h j Tj lim, V cv V h j T j, V 1 γ j θi j θ, t. θ Since V and h j T j, V is non-increasing with respect to V, then We also have h j T j, V h j T j, h j T j, V h j T j, lim h j V V T j,. h j Tj lim, V V h j T j, V lim V h j T j,v h jt j, V h j T j,v h jt j, V 16 h j λj, h j T j,
17 so N j h j T j, V cv h j Tj lim, V V h j T j, V 1 N j c h jt j, V N j c R 1. h j h j λj, 1 λj, h j T j, 1 Hence we have dw and the equality holds if and only if T j Tj, i j, V. Thus, the largest positively invariant subset of the state space where dw is the set S { T 1, i 1,..., T n, i n, V, A Ω T j T j, i j, V }. By LaSalle s invariance principle, this implies that all solutions with positive initial conditions approach S as t. On the other hand, any solution of 2.1 contained in S satisfies da ba and thus, At as t, that is, the solution approaches E. Therefore, we conclude that E is globally asymptotically stable in the state space Ω. We will now use a Lyapunov functional that is constructed as a combination of those used in [3] and [4] to prove the global stability of E. Theorem 5.2. If R < 1 < R, then the immune-free infected steady state E is globally asymptotically stable. Proof. Consider the Lyapunov functional W 1 + V H Tj N j T j Tj T j h j T j, V h j ξ j, V dξ j + γ j θi jθh i j θ, t i j θ V V + qh k A, 5.2 where Hx x 1 ln x and γ j is defined by 5.1. Then this functional is well-defined based on the uniform persistence of the system. The function H is non-negative and equals zero only at 1, so it is clear that W 1 has a global minimum at the equilibrium E. The derivative of W 1 along the solutions of 2.1 is dw 1 N j 1 h jtj, V λj h j T j, V T j h j T j, V + γ j θ 1 i j θ ij θ, t i j θ, t t + 1 V V + qh kav k h + A ba. P j θi j θ, t cv qav 17
18 By 4.3 we have 1 V V cv cv V n N jh j Tj, V V V. Also, notice that λ j Tj + h jtj, V and γ j θ 1 i j θ ij θ, t i j θ, t t γ j θi jθ t H i j θ, t i j θ γ j θi i j θ, t θ i j θ, t jθh i j θ + H i θ j θ γ j i i j, t jh i j γ j θi i j θ, t jθh i j θ θ N j h j Tj, V h j T j, V H h j Tj, V γ j θi i j θ, t jθh i j θ P j θi i j θ, t jθh i j θ, γ θi jθ + γ j θ di θ θ P j θi jθh i j θ, t i j θ since γ θi j θ + γ jθ di θ P j θi j θ, γ j N j, and i j, t h j T j, V. Then we obtain dw 1 N j 1 h jtj, V dj h j T j, V Tj T j + h j Tj, V h j T j, V N [ + N j h j Tj, V h j T j, V H h j Tj, V γ j θi i j θ, t jθh i j θ θ ] P j θi i j θ, t jθh i j θ + P j θi j θ, t V P j θi j θ, t V + N j h j Tj, V V N j Tj T j + N j h j Tj, V 1 h jtj, V h j T j, V ln h jt j, V h j T j, V + ln h jt j, V + + qav P j θi jθ h h + A 1 N j h j Tj, V V qav + qav + qhav h + A qhba k γ j θi i j θ, t jθh i j θ θ h jtj, V h j T j, V h jt j, V h j Tj, V + h jt j, V h j T j, V + h jt j, V h j Tj, V h j T j, V + 1 V V V i j θ, t V i j θ ln i jθ, t i j θ + qa k kv bh. 18
19 By H2, we know that h jtj,v h j T j,v < 1 when T j < T j and h jtj,v h j T j,v > 1 when T j > T j, so Tj T j 1 h jtj,v h j T j,v for all T j. It is clear that γ j θi j θh ij θ,t θ i j θ. We also have qav h h+a 1 since h h + A, and qa k kv bh qabh k R AN 1 since the hypothesis R < 1 < R implies that R AN < 1. Thus, dw 1 N j h j Tj, V h jtj, V h j T j, V h jt j, V h j Tj, V + h jt j, V h j T j, V + h jt j, V h j Tj, V ln h jt j, V h j T j, V + ln h jtj, V h j T j, V + 1 V V P j θi V i j θ, t jθ V i j θ 1 ln i jθ, t i j θ. [ ] N j h j Tj, V hj T j, V V H h j T j, V H V P j θi V i j θ, t jθh V i j θ. By H2 and H3, we have that V V V implies V h jt j,v h j T j,v 1, while V V implies 1 h jt j,v h j T j,v V V. Since Hx x 1 ln x is decreasing for x, 1 and increasing for x 1,, then H H V V for all V. Also, hj T j,v h j T j,v non-negativity of H implies that P j θi j θh V i j θ,t V i j θ. Thus, the above argument shows that dw 1 and that the equality holds if and only if T j Tj, i jθ, t i j θ, V V, and A. Therefore, by LaSalle s invariance principle we conclude that E is globally asymptotically stable. Similarly, inspired in [3] and [4], we make a modification to the Lyapunov functional 5.2 in order to prove that Ê is globally asymptotically stable whenever it exists. Theorem 5.3. If R > 1, then the antibody-immune infected steady state Ê is globally asymptotically stable. Proof. Consider the Lyapunov functional W 2 N j T j ˆT j + ˆV H Tj ˆT j h j ˆT j, ˆV h j ξ j, ˆV dξ j + VˆV q A h + ττ  dτ, k  τ + i j θ, t γ j θî j θh î j θ A h+ττ  τ where H and γ j are the same as for 5.2. It can be seen that the term q k dτ is positive for A  and equals zero for A Â, so W 2 has a global minimum at Ê. We calculate the derivative of W 2 along the solutions of 2.1 as follows: dw 2 N j 1 h j ˆT j, ˆV λj h j T j, ˆV T j h j T j, V 19
20 + + 1 ˆV V γ j θ 1 îjθ i j θ, t + q kv k h + AA Â h + A b. i j θ, t t P j θi j θ, t cv qav Following the same steps as in the proof of Theorem 5.2, we obtain dw 2 N j ˆT j T j 1 h j ˆT j, ˆV i j θ, t h j T j, ˆV γ j θî j θh î j θ θ + N j h j ˆT j, ˆV h j ˆT j, ˆV h j T j, ˆV h jt j, V h j ˆT j, ˆV + h jt j, V h j T j, ˆV + h jt j, V h j ˆT j, ˆV ln h jt j, V h j T j, ˆV + ln h j ˆT j, ˆV h j T j, ˆV + 1 VˆV + P j θî j θ ˆV i j θ, t V î j θ ln i jθ, t î j θ qav + qa ˆV + q kv k h + AA Â h + A b. We have N j ˆT j T j ˆV bh + Â/k, we also have so dw 2 1 h j ˆT j, ˆV h j T j, ˆV qav + qa ˆV + q k h + AA Â kv h + A b and γ j θi j θh ij θ,t i j θ θ. qav + qa ˆV + qv A Â qb qb h + AA + h + AÂ k [ k ] qb A Â h + A q ˆV k qb k A Â2, Since N j h j ˆT j, ˆV h j ˆT j, ˆV h j T j, ˆV h jt j, V h j ˆT j, ˆV + h jt j, V h j T j, ˆV + h jt j, V h j ˆT j, ˆV ln h jt j, V h j T j, ˆV + ln h j ˆT j, ˆV h j T j, ˆV + 1 VˆV ˆV ij θ, t P j θî j θ V î j θ 1 ln i jθ, t î j θ [ ] N j h j ˆT j, ˆV h j T j, V H h j T j, ˆV H VˆV ˆV ij θ, t P j θî j θh V î j θ. Thus, dw 2 and the equality holds if and only if T j ˆT j, i j θ, t î j θ, V ˆV, and A Â. Therefore, by LaSalle s invariance principle we conclude that Ê is globally asymptotically stable. 2
21 6 Related models In this section, we consider a special case in which the model 2.1 can be reduced to a system of ordinary differential equations with one time delay. For that, we will use a particular form of incidence rate, which is given by the saturated incidence function h j T j, V β jt j tv t 1+α j V t. We will assume that it takes time τ for virus to enter into the target cell and that there is an intracellular delay ω that describes the time required for an infected cell to produce virus, so the death rate of infected cells and viral production rate are given by the functions { if θ < τ, δ j θ 6.1 if θ τ and P j θ δ j { if θ < ω, P j if θ ω. 6.2 We also assume that the initial condition i j θ, i j θ satisfies lim θ i j θ, which means that the initial number of infected cells tends to as the infection age tends to infinity. If I j t i j θ, t is the number of infected cells in the j-th class at time t, then but by 2.5 we have di j t i j θ, t i j θ, t t ij θ, t θ θ θ + δ j θi j θ, t δ j θi j θ, t, so lim i jθ, t lim i σ j θ θ θ jθ t σ j θ t lim θ i jθ te δ j t, di j t i j, t β jt j tv t 1 + α j V t δ j θi j θ, t τ i j θ, t τ δ j i j θ, t. Using the function 6.2 we obtain P j θi j θ, t P j P j ω ω i j θ, t i j θ ω, t ωe δ jθω. Thus, with these assumptions, the model can be reformulated equivalently as the following 21
22 system: dt j t di j t dv t dat λ j T j t β jt j tv t 1 + α j V t β jt j tv t τ 1 + α j V t i j θ, t δj i j θ, t τ i j θ ω, t ωe δjθω cv t qatv t P j ω katv t h + At bat. 6.3 In this model we have two intracellular time delays, τ and ω. If we consider the case when the death rate of infected cells is constant, i.e., δ j θ δj for all θ, then the above system can be further simplified to the following system of delay differential equations: dt j t di j t dv t dat λ j T j t β jt j tv t 1 + α j V t β jt j tv t 1 + α j V t δ j I j t Pj e δ j ω I j t ω cv t qatv t katv t h + At bat, 6.4 which contains no integral terms and only one time delay. In this case, we have σ j θ e δ j θ and N j Pj e δ j ω /δj, so the basic reproduction numbers of virus and antibodies for 6.4 can be calculated as R 1 c Pj e δ j ω λ j β j δj, R 1 c 7 Numerical examples δ j kpj e δ j ω λ j β j kdj + α j + β j bh. We will now perform some numerical simulations to illustrate the analytical results obtained for the dynamics of model Model with one class of target cells We first simulate the case n 1 for system 2.1 with the particular forms of δ 1 θ and P 1 θ used in Section 6 to obtain the simplified system 6.4. Hence, we are considering a model with intracellular delay ω in viral production, constant death rate δ for infected cells, and saturated incidence rate. We use the parameters shown in Table 1, which are similar to those used in [7, 21, 22] for models of HIV infection with only one population of target cells. We consider three different values for the infection rate β, given in virions/day, in order to obtain the three possible scenarios for the dynamics of the model. For each case, we show the graphs of solutions with respect to time for several initial conditions, which are given in cells/ml for T t and It, in virions/ml for V t, and in µg for At. 22
23 Parameter Value Units λ 1 46 cells ml 1 day 1 d 1.46 day 1 α 1.5 δ1.1 day 1 P virions day 1 cell 1 ω.5 days c.25 day 1 q.3 µg 1 day 1 k.15 virion 1 day 1 h.2 b 2.9 day 1 Table 1: Parameters used for numerical simulations of system 6.4 with one class of target cells. Figure 1: Solutions of system 6.4 when R < R < 1 and E is globally asymptotically stable. 23
24 Figure 2: Solutions of system 6.4 when R < 1 < R and E is globally asymptotically stable. Figure 3: Solutions of system 6.4 when 1 < R < R and stable. Ê is globally asymptotically 24
25 We first consider the case when β 1 1 8, obtaining R.454 and R.1547, which are less than 1. According to Theorem 5.1, all solutions approach the uninfected equilibrium E as t, as we can see in Figure 1. For β 5 1 8, we have R , R.927, so according to Theorem 5.2 all solutions converge to the immune-free infected equilibrium E, as we can see in Figure 2. Lastly, if we take β 5 1 7, then we have R , R , so according to Theorem 5.3 all solutions converge to the immune infected equilibrium Ê, as shown in Figure Model with two classes of target cells We will now assume that there are two different populations of target cells: CD4 + T cells denoted by the subscript j 1 and macrophages j 2. For the CD4 + T cells, we consider a viral production kernel P 1 θ similar to that used in [2], which is given by { if θ < θ 1, P 1 θ 1 e rθ θ 1, if θ θ 1, P 1 where P virions day 1 cell 1, while r 1 is a saturation parameter and θ 1.25 days is the age at which reverse transcription is completed. We use a death rate for infected cells δ 1 θ of the form 6.1 with d 1.46 day 1, δ1 1.5 day 1 and a saturated infection rate h 1 T 1, V β 1T 1 V 1+α 1 V with α 1.5, β virions day 1. For the macrophages, we use a constant death rate δ 2 θ 1/14.1 day 1 and we assume that viral particles at produced at a rate P 2 θ θ [23]. We assume that the infection rate is given by the function h 2 T 2, V 1.19T 2V [24]. The rest 1 6 +V of the parameters used are given in Table 2 and are based on [2, 7, 23]. Figure 4 shows the solutions for each variable of system 2.1 with the above parameters, using the initial conditions T 1 5, T 2 1 6, i 1 θ, i 2 θ, e 1θ, V.1, A.1. In this case, the basic reproduction numbers are R for the population of CD4 + T cells and R for macrophages, so the total basic reproduction number of virus is R R 1 + R From this, we can see that the relative contribution of T cells to infection is considerably stronger than that of macrophages. The viral reproduction number with antibody response is R Since R > 1 and R > 1, all solutions converge to a positive steady state, which corresponds to the chronic phase of HIV, as can clearly be seen in the simulations. 8 Conclusions In this paper, we have studied an age-structured within-host viral dynamic model that includes multiple populations of target cells. Our model incorporates saturated antibody immune response and a general non-linear incidence rate, so it can be viewed as a generalization of several models that have been studied previously in the literature. This kind of models can be used to evaluate the relative contribution of viral production from different compartments of target cells, such as CD4 + cells, macrophages and dendritic cells, and this can help us improve our understanding of the dynamics of HIV infection. We have extended the results published in [3], where the authors studied a similar model with a saturated incidence that did not include multiple classes of target cells. We proved that the global dynamics of the model is completely determined by the basic 25
26 Parameter Value Units λ cells ml 1 day 1 d 1.1 day 1 λ cells ml 1 day 1 d 2.24 day 1 c 23 day 1 q.3 µg 1 day 1 k.15 virion 1 day 1 h.2 b 2.9 day 1 Table 2: Parameters used for numerical simulations of system 2.1 with two classes of target cells. Figure 4: Time evolution of the solution of system 2.1 with two classes of target cells. 26
27 reproduction number and the reproductive number of antibody response. If R is less than one, the only steady state of the model is the infection-free equilibrium and the infection is predicted to die out. When R < 1 < R, an infected equilibrium E appears which is globally asymptotically stable, so infection becomes chronic but with no persistent antibody immune response. Lastly, if R > 1, the antibody-immune infected equilibrium Ê appears and all solutions converge to it; this case corresponds to the scenario when the infection becomes chronic with a persistent antibody immune response. References 1. Nelson, P. W., Gilchrist, M. A., Coombs, D., Hyman, J. M. & Perelson, A. S. An age-structured model of HIV infection that allows for variations in the production rate of viral particles and the death rate of productively infected cells. Math. Biosci. Eng 1, Rong, L., Feng, Z. & Perelson, A. S. Mathematical analysis of age-structured HIV-1 dynamics with combination antiretroviral therapy. SIAM Journal on Applied Mathematics 67, Duan, X. & Yuan, S. Global dynamics of an age-structured virus model with saturation effects. Mathematical Methods in the Applied Sciences 4, Wang, X., Lou, Y. & Song, X. Age-Structured Within-Host HIV Dynamics with Multiple Target Cells. Studies in Applied Mathematics 138, Wang, S., Wu, J. & Rong, L. A note on the global properties of an age-structured viral dynamic model with multiple target cell populations. Mathematical biosciences and engineering: MBE 14, Yang, Y., Ruan, S. & Xiao, D. Global stability of an age-structured virus dynamics model with Beddington DeAngelis infection function. Math. Biosci. Eng 12, Tian, X. & Wang, J. Stability analysis for viral infection model with multitarget cells, Beddington-DeAngelis functional response, and humoral immunity. Discrete Dynamics in Nature and Society Wang, J., Tian, X. & Wang, X. Stability analysis for delayed viral infection model with multitarget cells and general incidence rate. International Journal of Biomathematics 9, Wang, J., Lang, J. & Li, F. Constructing Lyapunov functionals for a delayed viral infection model with multitarget cells, nonlinear incidence rate, state-dependent removal rate. J Nonlinear Sci Appl 9, Elaiw, A. M. Global properties of a class of HIV models. Nonlinear Analysis: Real World Applications 11, Elaiw, A. M. Global properties of a class of virus infection models with multitarget cells. Nonlinear Dynamics 69, Elaiw, A. M. & Azoz, S. A. Global properties of a class of HIV infection models with Beddington DeAngelis functional response. Mathematical Methods in the Applied Sciences 36, Wang, X., Chen, Y., Liu, S. & Song, X. A class of delayed virus dynamics models with multiple target cells. Computational and Applied Mathematics 32,
28 14. Culshaw, R. V. & Ruan, S. A delay-differential equation model of HIV infection of CD4+ T-cells. Mathematical biosciences 165, Li, D. & Ma, W. Asymptotic properties of a HIV-1 infection model with time delay. Journal of Mathematical Analysis and Applications 335, Georgescu, P. & Hsieh, Y.-H. Global stability for a virus dynamics model with nonlinear incidence of infection and removal. SIAM Journal on Applied Mathematics 67, Hale, J. K. Asymptotic behavior of dissipative systems 25 American Mathematical Soc., Magal, P. & Thieme, H. R. Eventual compactness for semiflows generated by nonlinear age-structured models. Commun. Pure Appl. Anal 3, Hale, J. K. & Waltman, P. Persistence in infinite-dimensional systems. SIAM Journal on Mathematical Analysis 2, Demasse, R. D. & Ducrot, A. An age-structured within-host model for multistrain malaria infections. SIAM Journal on Applied Mathematics 73, Pawelek, K. A., Liu, S., Pahlevani, F. & Rong, L. A model of HIV-1 infection with two time delays: mathematical analysis and comparison with patient data. Mathematical Biosciences 235, Wang, J., Pang, J., Kuniya, T. & Enatsu, Y. Global threshold dynamics in a fivedimensional virus model with cell-mediated, humoral immune responses and distributed delays. Applied Mathematics and Computation 241, Duffin, R. P. & Tullis, R. H. Mathematical models of the complete course of HIV infection and AIDS. Computational and Mathematical Methods in Medicine 4, Wasserstein-Robbins, F. A mathematical model of HIV infection: simulating T4, T8, macrophages, antibody, and virus via specific anti-hiv response in the presence of adaptation and tropism. Bulletin of mathematical biology 72,
A Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and Cure of Infected Cells in Eclipse Stage
Applied Mathematical Sciences, Vol. 1, 216, no. 43, 2121-213 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.216.63128 A Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and
More informationDYNAMICS OF EVOLUTIONARY COMPETITION BETWEEN BUDDING AND LYTIC VIRAL RELEASE STRATEGIES. Xiulan Lai and Xingfu Zou. (Communicated by Yang Kuang)
MATHEMATICA IOSCIENCES doi:1.3934/mbe.214.11.191 AND ENGINEERING Volume 11, Number 5, October 214 pp. 191 1113 DYNAMICS OF EVOUTIONARY COMPETITION ETWEEN UDDING AND YTIC VIRA REEASE STRATEGIES Xiulan ai
More informationGlobal Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response
Global Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response Gang Huang 1,2, Wanbiao Ma 2, Yasuhiro Takeuchi 1 1,Graduate School of Science and Technology, Shizuoka University,
More informationPermanence Implies the Existence of Interior Periodic Solutions for FDEs
International Journal of Qualitative Theory of Differential Equations and Applications Vol. 2, No. 1 (2008), pp. 125 137 Permanence Implies the Existence of Interior Periodic Solutions for FDEs Xiao-Qiang
More informationA MATHEMATICAL MODEL SIMULATING THE EFFECT OF VACCINE INDUCED RESPONSES ON HIV-1 INFECTION JOHN GUARDIOLA, GIUSEPPE IZZO, AND ANTONIA VECCHIO
Dynamic Systems and Applications 16 (2007) 407-424 A MATHEMATICAL MODEL SIMULATING THE EFFECT OF VACCINE INDUCED RESPONSES ON HIV-1 INFECTION JOHN GUARDIOLA, GIUSEPPE IZZO, AND ANTONIA VECCHIO Istituto
More informationA Stochastic Viral Infection Model with General Functional Response
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 9, 435-445 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.664 A Stochastic Viral Infection Model with General Functional Response
More informationGLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 304, pp. 1 8. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationGLOBAL DYNAMICS OF A MATHEMATICAL MODEL OF TUBERCULOSIS
CANADIAN APPIED MATHEMATICS QUARTERY Volume 13, Number 4, Winter 2005 GOBA DYNAMICS OF A MATHEMATICA MODE OF TUBERCUOSIS HONGBIN GUO ABSTRACT. Mathematical analysis is carried out for a mathematical model
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationDynamics of a diffusive viral model with Beddington-DeAngelis incidence rate and CTL immune response
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 (2017, 5753 5762 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Dynamics of a diffusive viral
More informationGlobal Analysis of a HCV Model with CTL, Antibody Responses and Therapy
Applied Mathematical Sciences Vol 9 205 no 8 3997-4008 HIKARI Ltd wwwm-hikaricom http://dxdoiorg/02988/ams20554334 Global Analysis of a HCV Model with CTL Antibody Responses and Therapy Adil Meskaf Department
More informationResearch Article On the Stability Property of the Infection-Free Equilibrium of a Viral Infection Model
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume, Article ID 644, 9 pages doi:.55//644 Research Article On the Stability Property of the Infection-Free Equilibrium of a Viral
More informationNonlinear Dynamical Systems Lecture - 01
Nonlinear Dynamical Systems Lecture - 01 Alexandre Nolasco de Carvalho August 08, 2017 Presentation Course contents Aims and purpose of the course Bibliography Motivation To explain what is a dynamical
More informationarxiv: v2 [q-bio.pe] 3 Oct 2018
Journal of Mathematical Biology manuscript No. (will be inserted by the editor Global stability properties of renewal epidemic models Michael T. Meehan Daniel G. Cocks Johannes Müller Emma S. McBryde arxiv:177.3489v2
More informationGlobal Stability of SEIRS Models in Epidemiology
Global Stability of SRS Models in pidemiology M. Y. Li, J. S. Muldowney, and P. van den Driessche Department of Mathematics and Statistics Mississippi State University, Mississippi State, MS 39762 Department
More informationATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.
ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS Emerson A. M. de Abreu Alexandre N. Carvalho Abstract Under fairly general conditions one can prove that
More informationDiscrete Population Models with Asymptotically Constant or Periodic Solutions
International Journal of Difference Equations ISSN 0973-6069, Volume 6, Number 2, pp. 143 152 (2011) http://campus.mst.edu/ijde Discrete Population Models with Asymptotically Constant or Periodic Solutions
More informationDelay SIR Model with Nonlinear Incident Rate and Varying Total Population
Delay SIR Model with Nonlinear Incident Rate Varying Total Population Rujira Ouncharoen, Salinthip Daengkongkho, Thongchai Dumrongpokaphan, Yongwimon Lenbury Abstract Recently, models describing the behavior
More informationA NONLOCAL REACTION-DIFFUSION POPULATION MODEL WITH STAGE STRUCTURE
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 11, Number 3, Fall 23 A NONLOCAL REACTION-DIFFUSION POPULATION MODEL WITH STAGE STRUCTURE Dedicated to Professor Paul Waltman on the occasion of his retirement
More informationEquilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains
Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School
More informationGLOBAL ATTRACTIVITY IN A CLASS OF NONMONOTONE REACTION-DIFFUSION EQUATIONS WITH TIME DELAY
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 1, Spring 2009 GLOBAL ATTRACTIVITY IN A CLASS OF NONMONOTONE REACTION-DIFFUSION EQUATIONS WITH TIME DELAY XIAO-QIANG ZHAO ABSTRACT. The global attractivity
More informationGLOBAL DYNAMICS OF A TIME-DELAYED DENGUE TRANSMISSION MODEL
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number 1, Spring 212 GLOBAL DYNAMICS OF A TIME-DELAYED DENGUE TRANSMISSION MODEL Dedicated to Herb Freedman on the occasion of his 7th birthday ZHEN WANG
More informationA delayed SIR model with general nonlinear incidence rate
Liu Advances in Difference Equations 2015) 2015:329 DOI 10.1186/s13662-015-0619-z R E S E A R C H Open Access A delayed SIR model with general nonlinear incidence rate Luju Liu * * Correspondence: lujuliu@126.com
More informationGlobal compact attractors and their tripartition under persistence
Global compact attractors and their tripartition under persistence Horst R. Thieme (joint work with Hal L. Smith) School of Mathematical and Statistical Science Arizona State University GCOE, September
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationPERMANENCE IN LOGISTIC AND LOTKA-VOLTERRA SYSTEMS WITH DISPERSAL AND TIME DELAYS
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 60, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) PERMANENCE
More informationResearch Article Global Stability of Humoral Immunity HIV Infection Models with Chronically Infected Cells and Discrete Delays
Discrete Dynamics in Nature and Society olume 215, Article ID 37968, 25 pages http://dx.doi.org/1.1155/215/37968 Research Article Global Stability of Humoral Immunity HI Infection Models with Chronically
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationCompetitive Exclusion in a Discrete-time, Size-structured Chemostat Model
Competitive Exclusion in a Discrete-time, Size-structured Chemostat Model Hal L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 1804, USA E-mail: halsmith@asu.edu Xiao-Qiang Zhao
More informationGlobal Analysis for an HIV Infection Model with CTL Immune Response and Infected Cells in Eclipse Phase
Article Global Analysis for an HIV Infection Model with CTL Immune Response and Infected Cells in Eclipse Phase Karam Allali 1 ID Jaouad Danane 1 ID and Yang Kuang 2 * ID 1 Laboratory of Mathematics and
More informationGlobal Analysis of an SEIRS Model with Saturating Contact Rate 1
Applied Mathematical Sciences, Vol. 6, 2012, no. 80, 3991-4003 Global Analysis of an SEIRS Model with Saturating Contact Rate 1 Shulin Sun a, Cuihua Guo b, and Chengmin Li a a School of Mathematics and
More informationEpidemic Model of HIV/AIDS Transmission Dynamics with Different Latent Stages Based on Treatment
American Journal of Applied Mathematics 6; 4(5): -4 http://www.sciencepublishinggroup.com/j/ajam doi:.648/j.ajam.645.4 ISSN: -4 (Print); ISSN: -6X (Online) Epidemic Model of HIV/AIDS Transmission Dynamics
More informationGlobal dynamics of a predator-prey system with Holling type II functional response
4 Nonlinear Analysis: Modelling and Control, 011, Vol. 16, No., 4 53 Global dynamics of a predator-prey system with Holling type II functional response Xiaohong Tian, Rui Xu Institute of Applied Mathematics
More informationLOCAL AND GLOBAL STABILITY OF IMPULSIVE PEST MANAGEMENT MODEL WITH BIOLOGICAL HYBRID CONTROL
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 11 No. II (August, 2017), pp. 129-141 LOCAL AND GLOBAL STABILITY OF IMPULSIVE PEST MANAGEMENT MODEL WITH BIOLOGICAL HYBRID CONTROL
More informationChapter 3 Pullback and Forward Attractors of Nonautonomous Difference Equations
Chapter 3 Pullback and Forward Attractors of Nonautonomous Difference Equations Peter Kloeden and Thomas Lorenz Abstract In 1998 at the ICDEA Poznan the first author talked about pullback attractors of
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationA REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES KING-YEUNG LAM
More informationGlobal attractors, stability, and population persistence
Global attractors, stability, and population persistence Horst R. Thieme (partly joint work with Hal L. Smith*) School of Mathematical and Statistical Sciences Arizona State University partially supported
More informationVariational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized
More informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationMathematical models on Malaria with multiple strains of pathogens
Mathematical models on Malaria with multiple strains of pathogens Yanyu Xiao Department of Mathematics University of Miami CTW: From Within Host Dynamics to the Epidemiology of Infectious Disease MBI,
More informationStability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay Wassim M. Haddad and VijaySekhar Chellaboina School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationMathematical Analysis of HIV/AIDS Prophylaxis Treatment Model
Applied Mathematical Sciences, Vol. 12, 2018, no. 18, 893-902 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.8689 Mathematical Analysis of HIV/AIDS Prophylaxis Treatment Model F. K. Tireito,
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationDynamical systems. January 10, 2016
Dynamical systems January 10, 2016 1 Introduction In this text the term dynamical system means nothing else than a system of ordinary differential equations. Why is this term used here? The aim is to emphasize
More informationGlobal Stability of a Computer Virus Model with Cure and Vertical Transmission
International Journal of Research Studies in Computer Science and Engineering (IJRSCSE) Volume 3, Issue 1, January 016, PP 16-4 ISSN 349-4840 (Print) & ISSN 349-4859 (Online) www.arcjournals.org Global
More informationStability Analysis of a SIS Epidemic Model with Standard Incidence
tability Analysis of a I Epidemic Model with tandard Incidence Cruz Vargas-De-León Received 19 April 2011; Accepted 19 Octuber 2011 leoncruz82@yahoo.com.mx Abstract In this paper, we study the global properties
More informationCan multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species
Can multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species Xingfu Zou Department of Applied Mathematics University of Western Ontario London, Ontario, Canada (Joint work
More informationStochastic Viral Dynamics with Beddington-DeAngelis Functional Response
Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response Junyi Tu, Yuncheng You University of South Florida, USA you@mail.usf.edu IMA Workshop in Memory of George R. Sell June 016 Outline
More informationTopic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis
Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of
More informationAsymptotic behavior of Ginzburg-Landau equations of superfluidity
Communications to SIMAI Congress, ISSN 1827-9015, Vol. 3 (2009) 200 (12pp) DOI: 10.1685/CSC09200 Asymptotic behavior of Ginzburg-Landau equations of superfluidity Alessia Berti 1, Valeria Berti 2, Ivana
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Integrating Mathematical and Statistical Models Recap of mathematical models Models and data Statistical models and sources of
More information6. Age structure. for a, t IR +, subject to the boundary condition. (6.3) p(0; t) = and to the initial condition
6. Age structure In this section we introduce a dependence of the force of infection upon the chronological age of individuals participating in the epidemic. Age has been recognized as an important factor
More informationThreshold solutions and sharp transitions for nonautonomous parabolic equations on R N
Threshold solutions and sharp transitions for nonautonomous parabolic equations on R N P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract This paper is devoted to
More informationSYNCHRONIZATION OF NONAUTONOMOUS DYNAMICAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 003003, No. 39, pp. 1 10. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu login: ftp SYNCHRONIZATION OF
More informationSpread of Viral Plaque
Spread of Viral Plaque Don Jones, Gergely Rost, Hal Smith, Horst Thieme A R I Z O N A S T A T E U N I V E R S I T Y Everything Disperses to Miami () Spread of Viral Plaque December 212 1 / 2 Outline 1
More informationResearch Article Frequent Oscillatory Behavior of Delay Partial Difference Equations with Positive and Negative Coefficients
Hindawi Publishing Corporation Advances in Difference Equations Volume 2010, Article ID 606149, 15 pages doi:10.1155/2010/606149 Research Article Frequent Oscillatory Behavior of Delay Partial Difference
More informationPREPUBLICACIONES DEL DEPARTAMENTO DE MATEMÁTICA APLICADA UNIVERSIDAD COMPLUTENSE DE MADRID MA-UCM
PREPUBLICACIONES DEL DEPARTAMENTO DE MATEMÁTICA APLICADA UNIVERSIDAD COMPLUTENSE DE MADRID MA-UCM 2009-13 Extremal equilibria for monotone semigroups in ordered spaces with application to evolutionary
More informationON THE HYPERBOLIC RELAXATION OF THE CAHN-HILLIARD EQUATION IN 3-D: APPROXIMATION AND LONG TIME BEHAVIOUR
ON THE HYPERBOLIC RELAXATION OF THE CAHN-HILLIARD EQUATION IN 3-D: APPROXIMATION AND LONG TIME BEHAVIOUR ANTONIO SEGATTI Abstract. In this paper we consider the hyperbolic relaxation of the Cahn-Hilliard
More informationExistence of homoclinic solutions for Duffing type differential equation with deviating argument
2014 9 «28 «3 Sept. 2014 Communication on Applied Mathematics and Computation Vol.28 No.3 DOI 10.3969/j.issn.1006-6330.2014.03.007 Existence of homoclinic solutions for Duffing type differential equation
More informationESTIMATES ON THE NUMBER OF LIMIT CYCLES OF A GENERALIZED ABEL EQUATION
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 ESTIMATES ON THE NUMBER OF LIMIT CYCLES OF A GENERALIZED ABEL EQUATION NAEEM M.H. ALKOUMI
More informationSTUDY OF THE DYNAMICAL MODEL OF HIV
STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application
More informationExistence and uniqueness of solutions for a diffusion model of host parasite dynamics
J. Math. Anal. Appl. 279 (23) 463 474 www.elsevier.com/locate/jmaa Existence and uniqueness of solutions for a diffusion model of host parasite dynamics Michel Langlais a and Fabio Augusto Milner b,,1
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationDynamical behavior of a pest management model with impulsive effect and nonlinear incidence rate*
Volume 30, N. 2, pp. 381 398, 2011 Copyright 2011 SBMAC ISSN 0101-8205 www.scielo.br/cam Dynamical behavior of a pest management model with impulsive effect and nonlinear incidence rate* XIA WANG 1, ZHEN
More informationPersistence theory applied to Keen s model a link between mathematical biology and mathematical economics
Persistence theory applied to Keen s model a link between mathematical biology and mathematical economics Jianhong Wu and Xiang-Sheng Wang Mprime Centre for Disease Modelling York University, Toronto Persistence
More informationStability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates
Published in Applied Mathematics and Computation 218 (2012 5327-5336 DOI: 10.1016/j.amc.2011.11.016 Stability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates Yoichi Enatsu
More informationA Note on Some Properties of Local Random Attractors
Northeast. Math. J. 24(2)(2008), 163 172 A Note on Some Properties of Local Random Attractors LIU Zhen-xin ( ) (School of Mathematics, Jilin University, Changchun, 130012) SU Meng-long ( ) (Mathematics
More informationOBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 27 (27, No. 6, pp. 2. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu OBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationERRATA: Probabilistic Techniques in Analysis
ERRATA: Probabilistic Techniques in Analysis ERRATA 1 Updated April 25, 26 Page 3, line 13. A 1,..., A n are independent if P(A i1 A ij ) = P(A 1 ) P(A ij ) for every subset {i 1,..., i j } of {1,...,
More informationNontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
More informationarxiv: v3 [math.ds] 22 Feb 2012
Stability of interconnected impulsive systems with and without time-delays using Lyapunov methods arxiv:1011.2865v3 [math.ds] 22 Feb 2012 Sergey Dashkovskiy a, Michael Kosmykov b, Andrii Mironchenko b,
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
More informationImpulsive differential equations with applications to infectious diseases
Impulsive differential equations with applications to infectious diseases Rachelle Miron Thesis submitted to the Faculty of Graduate and Postdoctoral Studies in partial fulfillment of the requirements
More informationA Quasi-Linear Parabolic Partial Differential Equation with Accretive Property
ONLINE ISSN 8-749 : Volume 3, Issue, 433-438 A Quasi-Linear Parabolic Partial Differential Equation with Accretive Property Aminu U. Bawa *, Micheal O. Egwurube and Murtala M. Ahmad 3 Department of Computer
More informationConservation law equations : problem set
Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions
More informationA Model for Hepatitis B with Chronological and Infection Ages
Applied Mathematical Sciences, Vol. 7, 213, no. 12, 5977-5993 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.213.3742 A Model for Hepatitis B with Chronological and Infection Ages Kouakep Tchaptchie
More informationResearch Article A Delayed Epidemic Model with Pulse Vaccination
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 2008, Article ID 746951, 12 pages doi:10.1155/2008/746951 Research Article A Delayed Epidemic Model with Pulse Vaccination
More informationDynamics of an HIV pathogenesis model with CTL immune response and two saturated rates
ISSN 746-7233 England UK World Journal of Modelling and Simulation Vol. (24) No. 3 pp. 25-223 Dynamics of an HIV pathogenesis model with CTL immune response and two saturated rates Youssef Tait Khalid
More informationMyopic Models of Population Dynamics on Infinite Networks
Myopic Models of Population Dynamics on Infinite Networks Robert Carlson Department of Mathematics University of Colorado at Colorado Springs rcarlson@uccs.edu June 30, 2014 Outline Reaction-diffusion
More informationGlobal unbounded solutions of the Fujita equation in the intermediate range
Global unbounded solutions of the Fujita equation in the intermediate range Peter Poláčik School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA Eiji Yanagida Department of Mathematics,
More informationGlobal Analysis of a Mathematical Model of HCV Transmission among Injecting Drug Users and the Impact of Vaccination
Applied Mathematical Sciences, Vol. 8, 2014, no. 128, 6379-6388 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.48625 Global Analysis of a Mathematical Model of HCV Transmission among
More informationSingular boundary value problems
Department of Mathematics, Faculty of Science, Palacký University, 17. listopadu 12, 771 46 Olomouc, Czech Republic, e-mail: irena.rachunkova@upol.cz Malá Morávka, May 2012 Overview of our common research
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationMonotone variational inequalities, generalized equilibrium problems and fixed point methods
Wang Fixed Point Theory and Applications 2014, 2014:236 R E S E A R C H Open Access Monotone variational inequalities, generalized equilibrium problems and fixed point methods Shenghua Wang * * Correspondence:
More informationViral evolution model with several time scales
Viral evolution model with several time scales AA Archibasov Samara National Research University, 4 Moskovskoe Shosse, 4486, Samara, Russia Abstract In this paper a viral evolution model with specific
More informationA convergence result for an Outer Approximation Scheme
A convergence result for an Outer Approximation Scheme R. S. Burachik Engenharia de Sistemas e Computação, COPPE-UFRJ, CP 68511, Rio de Janeiro, RJ, CEP 21941-972, Brazil regi@cos.ufrj.br J. O. Lopes Departamento
More informationLIFE SPAN OF BLOW-UP SOLUTIONS FOR HIGHER-ORDER SEMILINEAR PARABOLIC EQUATIONS
Electronic Journal of Differential Equations, Vol. 21(21), No. 17, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu LIFE SPAN OF BLOW-UP
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationOn CTL Response against Mycobacterium tuberculosis
Applied Mathematical Sciences, Vol. 8, 2014, no. 48, 2383-2389 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.43150 On CTL Response against Mycobacterium tuberculosis Eduardo Ibargüen-Mondragón
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More information