Loading [MathJax]/jax/element/mml/optable/MathOperators.js
Research article Special Issues

Stability and bifurcation analysis of a multi-delay model for mosaic disease transmission

  • A mathematical model is developed for analysis of the spread of mosaic disease in plants, which account for incubation period and latency that are represented by time delays. Feasibility and stability of different equilibria are studied analytically and numerically. Conditions that determine the type of behavior exhibited by the system are found in terms of various parameters. We have derived the basic reproduction number and identify the conditions resulting in eradication of the disease, as well as those that lead to the emergence of stable oscillations in the population of infected plants, as a result of Hopf bifurcation of the endemic equilibrium. Numerical simulations are performed to verify the analytical results and also to illustrate different dynamical regimes that can be observed in the system. In this research, the stabilizing role of both the time delay has been established i.e. when delay time is large, disease will persist if the infection rate is higher. The results obtained here are useful for plant disease management.

    Citation: Fahad Al Basir, Konstantin B. Blyuss, Ezio Venturino. Stability and bifurcation analysis of a multi-delay model for mosaic disease transmission[J]. AIMS Mathematics, 2023, 8(10): 24545-24567. doi: 10.3934/math.20231252

    Related Papers:

    [1] Xin-You Meng, Miao-Miao Lu . Stability and bifurcation of a delayed prey-predator eco-epidemiological model with the impact of media. AIMS Mathematics, 2023, 8(7): 17038-17066. doi: 10.3934/math.2023870
    [2] Sahabuddin Sarwardi, Hasanur Mollah, Aeshah A. Raezah, Fahad Al Basir . Direction and stability of Hopf bifurcation in an eco-epidemic model with disease in prey and predator gestation delay using Crowley-Martin functional response. AIMS Mathematics, 2024, 9(10): 27930-27954. doi: 10.3934/math.20241356
    [3] Yougang Wang, Anwar Zeb, Ranjit Kumar Upadhyay, A Pratap . A delayed synthetic drug transmission model with two stages of addiction and Holling Type-II functional response. AIMS Mathematics, 2021, 6(1): 1-22. doi: 10.3934/math.2021001
    [4] San-Xing Wu, Xin-You Meng . Dynamics of a delayed predator-prey system with fear effect, herd behavior and disease in the susceptible prey. AIMS Mathematics, 2021, 6(4): 3654-3685. doi: 10.3934/math.2021218
    [5] Heping Jiang . Complex dynamics induced by harvesting rate and delay in a diffusive Leslie-Gower predator-prey model. AIMS Mathematics, 2023, 8(9): 20718-20730. doi: 10.3934/math.20231056
    [6] Muhammad Rifqy Adha Nurdiansyah, Kasbawati, Syamsuddin Toaha . Stability analysis and numerical simulation of rabies spread model with delay effects. AIMS Mathematics, 2024, 9(2): 3399-3425. doi: 10.3934/math.2024167
    [7] Ting Gao, Xinyou Meng . Stability and Hopf bifurcation of a delayed diffusive phytoplankton-zooplankton-fish model with refuge and two functional responses. AIMS Mathematics, 2023, 8(4): 8867-8901. doi: 10.3934/math.2023445
    [8] Qinghui Liu, Xin Zhang . Chaos detection in predator-prey dynamics with delayed interactions and Ivlev-type functional response. AIMS Mathematics, 2024, 9(9): 24555-24575. doi: 10.3934/math.20241196
    [9] Haowen Gong, Huijun Xiang, Yifei Wang, Huaijin Gao, Xinzhu Meng . Strategy evolution of a novel cooperative game model induced by reward feedback and a time delay. AIMS Mathematics, 2024, 9(11): 33161-33184. doi: 10.3934/math.20241583
    [10] Changjin Xu, Chaouki Aouiti, Zixin Liu, Qiwen Qin, Lingyun Yao . Bifurcation control strategy for a fractional-order delayed financial crises contagions model. AIMS Mathematics, 2022, 7(2): 2102-2122. doi: 10.3934/math.2022120
  • A mathematical model is developed for analysis of the spread of mosaic disease in plants, which account for incubation period and latency that are represented by time delays. Feasibility and stability of different equilibria are studied analytically and numerically. Conditions that determine the type of behavior exhibited by the system are found in terms of various parameters. We have derived the basic reproduction number and identify the conditions resulting in eradication of the disease, as well as those that lead to the emergence of stable oscillations in the population of infected plants, as a result of Hopf bifurcation of the endemic equilibrium. Numerical simulations are performed to verify the analytical results and also to illustrate different dynamical regimes that can be observed in the system. In this research, the stabilizing role of both the time delay has been established i.e. when delay time is large, disease will persist if the infection rate is higher. The results obtained here are useful for plant disease management.



    In this manuscript we discuss the existence of multiple non-negative solutions of the boundary value problem (BVP)

    {(p(t)u(t))+f(t,u(t))=0, t0,αu(0)βu(0)=H[u], u(+)=0, (1.1)

    where p,f are continuous functions on their domains with p>0, α>0, β0 and R>0 exists such that f(t,v)0 for (t,v)[0,R]×[0,). Further, H is a suitable functional on the space C[0,R]. No sign condition is assumed on f for t>R.

    The functional formulation of the boundary conditions (BCs) in (1.1) covers, as special cases, the interesting setting of (not necessarily linear) multi-point and integral BCs. A wide literature has been devoted on this topic; for the case of the bounded intervals we refer the reader to the recent paper [11] and references therein. Nonlocal BCs have also received attention in the case of unbounded intervals, one example is the paper [22] where the case β=0 and H[u]=u(ξ) (ξ(0,R]) has been investigated by means of a fixed point theorem of Leggett–Williams type. The case β=0 and H[u]=m2i=1ηiu(ξi) (ξi(0,R]) has been studied, for instance, in [16] for a differential equation with a different linear part in (1.1), by means of the well-known Krasnosel'skiĭ-Guo theorem on cone-compressions and cone-expansions. The approach in [16,22] is to rewrite the BVP into an integral equation of Hammerstein type on the half-line, and this approach is possible due to the linearity of the BCs considered. Our methodology to solve the BVP (1.1) is different, and fits within the line of the papers [7,8,10,20,21]. It consists in considering two auxiliary BVPs separately, the first one on the compact interval [0,R], where the functional BC acts and f is nonnegative, and the second one on the half-line [R,), where f is allowed to change its sign. Unlike the above cited articles, in which the problem of gluing the solutions is solved with some continuity arguments and an analysis in the phase space, here both the auxiliary problems have the same slope condition in the junction point (namely the condition u(R)=0), which simplifies the arguments. This kind of decomposition is some sort of an analogue of one employed by Boucherif and Precup [2] utilized for equations with nonlocal initial conditions, where the associated nonlinear integral operator is decomposed into two parts, one of Fredholm-type (that takes into account the functional conditions) and another one of Volterra-type.

    We make the following assumptions on the terms that occur in (1.1).

    p:[0,+)(0,+) is continuous, with

    P=+R1p(t)dt<+. (1.2)

    f:[0,+)×[0,+)R is continuous, with f(t,v)0 for (t,v)[0,R]×[0,+), f(t,0)=0 for t0.

    ● There exist two continuous functions b1,b2:[R,+)R, with b20, and two nondecreasing C1-functions F1,F2:[0,+)[0,+), with Fj(0)=0, Fj(v)>0 for v>0, j=1,2, such that

    b1(t)F1(v)f(t,v)b2(t)F2(v),for all (t,v)[R,+)×[0,+). (1.3)

    ● For j=1,2

    lim supv0+Fj(v)v<+. (1.4)

    ● Let b1,b+1 be the negative and the positive part of b1, i.e. b1(t)=max{0,b1(t)}, b+1(t)=max{0,b1(t)}. Then

    B1=+Rb1(t)dt<+, (1.5)
    +R1p(t)tRb+1(s)dsdt=+. (1.6)

    Notice that (1.2), (1.6) imply

    Rb+1(t)dt=+,

    and, in particular, the function b1 cannot be negative in a neighbourhood of infinity.

    If (1.6) is not satisfied, then (1.1) may have no solution, since it may happen that no positive solution of the equation in (1.1), satisfying the functional boundary conditions, tends to zero as t tends to infinity. This happens, for instance, if f(t,u)0 for all (t,u)[R,+)×R (i.e., b+1=b2=0) and u(R)0, in which case all positive solutions are nondecreasing. However, as we will show in Section 3 (see also Theorem 4.1), when the condition (1.6) is not satisfied, then our approach leads to sufficient conditions for the existence of a bounded non-negative solution of the differential equation in (1.1), namely a solution of the problem

    {(p(t)u(t))+f(t,u(t))=0, t0,αu(0)βu(0)=H[u], u0 and bounded on [0,+). (1.7)

    The problem of the existence and multiplicity of the solutions for the equation in (1.1), which are non-negative in the interval [0,R] and satisfy the functional BCs and the additional assumptions at u(R)=0, is considered in Section 2 and is solved by means of the classical fixed point index for compact maps. A BVP on [R,+) is examined in Section 3, where we deal with the existence of positive global solutions which have zero initial slope and are bounded or tend to zero at infinity. This second problem is solved by using a fixed point theorem for operators defined in a Frechét space, by a Schauder's linearization device, see [5,Theorem 1.3], and does not require the explicit form of the fixed point operator, but only some a-priori bounds. These estimates are obtained using some properties of principal and nonprincipal solutions of auxiliary second-order linear equations, see [15,Chapter 11] and [9]. Finally, the existence and multiplicity of solutions for the BVP (1.1) and (1.7) is obtained in Section 4, thanks to the fact that the problem in [R,+) has at least a solution for every initial value u(R) sufficiently small. A couple of examples completes the paper, illustrating our two main results.

    In this Section we investigate the existence of multiple positive solutions of the BVP

    {(p(t)u(t))+f(t,u(t))=0, t[0,R],αu(0)βu(0)=H[u], u(R)=0. (2.1)

    First of all we recall some results regarding the linear BVP

    {(p(t)u(t))=0,t[0,R],αu(0)βu(0)=0, u(R)=0. (2.2)

    It is known, see for example [19], that the Green's function k for the BVP (2.2) is given by

    k(t,s):=1α{βp(0)+αs01p(μ)dμ,st,βp(0)+αt01p(μ)dμ, st.

    and satisfies the inequality (see [19,Lemma 2.1])

    c(t)Φ(s)k(t,s)Φ(s), (t,s)[0,R]2,

    where

    Φ(s):=βαp(0)+s01p(μ)dμ, and c(t):=βp(0)+αt01p(μ)dμβp(0)+αR01p(μ)dμ. (2.3)

    Note that the constant function 1α solves the BVP

    {(p(t)u(t))=0,t[0,R],αu(0)βu(0)=1, u(R)=0.

    We associate to the BVP (2.1) the perturbed Hammerstein integral equation

    u(t)=Fu(t)+H[u]α:=Tu(t), (2.4)

    where

    Fu(t):=R0k(t,s)f(s,u(s))ds.

    We seek fixed points of the operator T in a suitable cone of the space of continuous functions C[0,R], endowed with the usual norm w:=max{|w(t)|,t[0,R]}.

    We recall that a cone K in a Banach space X is a closed convex set such that λxK for xK and λ0 and K(K)={0}. In the following Proposition we recall the main properties of the classical fixed point index for compact maps, for more details see [1,12]. In what follows the closure and the boundary of subsets of a cone K are understood to be relative to K.

    Proposition 2.1. Let X be a real Banach space and let KX be a cone. Let D be an open bounded set of X with 0DK and ¯DKK, where DK=DK. Assume that T:¯DKK is a compact operator such that xTx for xDK. Then the fixed point index iK(T,DK) has the following properties:

    (i) If there exists eK{0} such that xTx+λe for all xDK and all λ>0, then iK(T,DK)=0.

    (iii) If Txλx for all xDK and all λ>1, then iK(T,DK)=1.

    (iv) Let D1 be open bounded in X such that ¯D1KDK. If iK(T,DK)=1 and iK(T,D1K)=0, then T has a fixed point in DK¯D1K. The same holds if iK(T,DK)=0 and iK(T,D1K)=1.

    The assumptions above allow us to work in the cone

    K:={uC[0,R]:u0,mint[a,b]u(t)cu}, (2.5)

    a type of cone firstly used by Krasnosel'skiĭ, see [17], and D. Guo, see e.g. [12]. In (2.5) [a,b] is a suitable subinterval of [0,R] and c:=mint[a,b]c(t), with c(t) given by (2.3). We have freedom of choice of [a,b], with the restriction a>0 when β=0. Note also that the constant function equal to r0 (that we denote, with abuse of notation r) belongs to K, so K{0}.

    Regarding the functional H we assume that

    H:K[0,+) is continuous and maps bounded sets in bounded sets.

    With these ingredients it is routine to show that T leaves K invariant and is compact. We make use of the following open bounded set (relative to K)

    Kρ:={uK:u<ρ}.

    We now employ some local upper and lower estimates for the functional H, in the spirit of [13,14]. We begin with a condition which implies that the index is 1.

    Lemma 2.2. Assume that

    (I1ρ) there exists ρ>0, such that the following algebraic inequality holds:

    1m¯fρ+1α¯Hρ<ρ, (2.6)

    where

    ¯fρ:=max(t,u)[0,R]×[0,ρ]f(t,u), ¯Hρ:=supuKρH[u]  and 1m:=supt[0,R]R0k(t,s)ds.

    Then iK(T,Kρ) is 1.

    Proof. Note that if uKρ then we have 0u(t)ρ for every t[0,R]. We prove that μuTu for every μ1 and uKρ. In fact, if this does not happen, there exist μ1 and uKρ such that, for every t[0,R], we have

    μu(t)=Tu(t)=Fu(t)+1αH[u].

    Then we obtain, for t[0,R],

    μu(t)R0k(t,s)¯fρds+1α¯Hρ1m¯fρ+1α¯Hρ. (2.7)

    Taking the supremum for t[0,R] in (2.7) and using the inequality (2.6) we obtain μρ<ρ, a contradiction that proves the result.

    Now we give a condition which implies that the index is 0 on the set Kρ.

    Lemma 2.3. Assume that

    (I0ρ) there exists ρ>0 such that the following algebraic inequality holds:

    1Mf_ρ+1αH_ρ[u]>ρ, (2.8)

    where

    f_ρ:=min(t,u)[a,b]×[cρ,ρ]f(t,u), H_ρ:=infuKρH[u]  and 1M:=inft[a,b]bak(t,s)ds.

    Then iK(T,Kρ) is 0.

    Proof. Note that the constant function 1 belongs to K. We prove that uTu+λ1 for every uKρ and for every λ0. If this is false, there exist uKρ and λ0 such that u=Tu+λ1. Then we have, for t[a,b],

    ρu(t)=Fu(t)+1αH[u]+λ1Fu(t)+1αH[u]bak(t,s)f(s,u(s))ds+1αH[u]bak(t,s)f_ρds+1αH_ρ[u]1Mf_ρ+1αH_ρ[u]. (2.9)

    Using the inequality (2.8) in (2.9) we obtain ρ>ρ, a contradiction that proves the result.

    In view of the Lemmas above, we may state our result regarding the existence of one or more nontrivial solutions. Here, for brevity, we provide sufficient conditions for the existence of one, two or three solutions. It is possible to obtain more solutions, by adding more conditions of the same type, see for example [18].

    Theorem 2.4. The BVP (2.1) has at least one non-negative solution u1, with ρ1<u1(R)<ρ2 if either of the following conditions holds.

    (S1) There exist ρ1,ρ2(0,+) with ρ1<ρ2 such that (I0ρ1) and (I1ρ2) hold.

    (S2) There exist ρ1,ρ2(0,+) with ρ1<ρ2 such that (I1ρ1) and (I0ρ2) hold.

    The BVP (2.1) has at least two non-negative solutions u1 and u2, with ρ1<u1(R)<ρ2<u2(R)<ρ3, if either of the following conditions holds.

    (S3) There exist ρ1,ρ2,ρ3(0,+) with ρ1<ρ2<ρ3 such that (I0ρ1),(I1ρ2)and(I0ρ3) hold.

    (S4) There exist ρ1,ρ2,ρ3(0,+) with ρ1<ρ2<ρ3 such that (I1ρ1),(I0ρ2)and(I1ρ3) hold.

    The BVP (2.1) has at least three non-negative solutions u1, u2 and u3, with ρ1<u1(R)<ρ2<u2(R)<ρ3<u3(R)<ρ4 if either of the following conditions holds.

    (S5) There exist ρ1,ρ2,ρ3,ρ4(0,+) with ρ1<ρ2<ρ3<ρ4 such that (I0ρ1),(I1ρ2),(I0ρ3) and (I1ρ4) hold.

    (S6) There exist ρ1,ρ2,ρ3,ρ4(0,+) with ρ1<ρ2<ρ3<ρ4 such that (I1ρ1),(I0ρ2),(I1ρ3) and (I0ρ4) hold.

    Proof. Assume condition (S1) holds, then, by Lemma 2.3, we have iK(T,Kρ1)=0 and, by Lemma 2.2, iK(T,Kρ2)=1. By Proposition 2.1 we obtain a solution u1 for the integral equation (2.4) in Kρ2¯Kρ1. Furthermore note that

    tk(t,s):=1α{0, s<t,1p(t), s>t,

    and therefore u10 in [0,R], which, in turn, implies u1(R)=u1.

    Assume now that condition (S3) holds, then we obtain in addition that iK(T,Kρ3)=0. By Proposition 2.1 we obtain the existence of a second solution u2 of the integral equation (2.4) in Kρ3¯Kρ2. A similar argument as above holds for the monotonicity of u2.

    The remaining cases are dealt with in a similar way.

    In this section we state sufficient conditions for the existence of solutions of the following BVP

    (p(t)u(t))+f(t,u(t))=0, tR, (3.1)
    u(R)=0,u(R)=u0,u(t)>0,limtu(t)=0, (3.2)

    where u0>0 is a given constant. The BVP (3.1), (3.2) involves both initial and asymptotic conditions, and a global condition (i.e., the positivity on the whole half-line). The continuability at infinity of solutions of (3.1) is not a simple problem, see for example [4]. For instance, the Emden–Fowler equation

    (p(t)u(t))+g(t)|u(t)|βsgn(u(t))=0, (3.3)

    if β>1 and g is allowed to take negative values, has solutions which tend to infinity in finite time, see [3,4]. Moreover, again in the superlinear case β>1, if g is non-negative with isolated zeros, then (3.3) may have solutions which change sign infinitely many times in the left neighborhood of some ˉt>R, and so these solutions are not continuable to infinity, see [6]. Further, even if global solutions exist (that is, solutions which are defined in the whole half-line [R,+)), their positivity is not guaranteed in general. Indeed, (3.1) may exhibit the coexistence of nonoscillatory and oscillatory solutions; further, nonoscillatory solutions may have an arbitrary large number of zeros.

    The problem (3.1), (3.2) has been consider in [9] for nonlinear equations with p-laplacian operator and nonlinear term f(t,u(t))=b(t)F(u(t)). We address the reader to such a paper for a complete discussion on the issues related to the BVP (3.1), (3.2) and for a review of the existing literature on related problems. The approach used in [9] to solve the BVP was based on a combination of the Schauder's (half)-linearization device, a fixed point result in the Fréchet space of continuous functions on [R,+), and comparison results for principal and nonprincipal solutions of suitable auxiliary half-linear equations, which allow to find good upper and lower bounds for the solutions of the (half)-linearized problem. The same approach, with minor modifications, allow us to treat also the present case with a general nonlinearity f(t,u(t)), under the assumptions (1.3), (1.4). In the following Proposition we recall the fixed-point result [9,Theorem 1], based on [5,Theorem 1.3], in the form suitable for the present problem.

    Proposition 3.1. Let J=[t0,). Consider the BVP

    {(p(t)u)+f(t,u)=0,tJ,uS, (3.4)

    where f is a continuous function on J×R and S is a subset of C1(J,R). Let g be a continuous function on J×R2 such that

    g(t,c,c)=f(t,c)for all(t,c)J×R,

    and assume that there exist a closed convex subset Ω of C1(J,R) and a bounded closed subset S1 of SΩ which make the problem

    {(p(t)y)+g(t,y,q)=0,tJ,yS1 (3.5)

    uniquely solvable for all qΩ. Then the BVP (3.4) has at least one solution in Ω.

    In view of this result, no topological properties of the fixed-point operator are needed to be checked, since they are a direct consequence of a-priori bounds for the solutions of the "linearized" problem (3.5).

    We state here the existence results in the form which will be used in the next section, addressing to [9,Theorem 2] for the general result, in case of a factored nonlinearity. For reader's convenience we provide a short proof, focusing only on those points which require some adjustments due to the more general nonlinearity. We point out that the present results are obtained by using the Euler equation

    t2y+nty+(n12)2y=0,n>1, (3.6)

    as a Sturm majorant of the linearized auxiliary equation (see also [9,Corollary 3]), but any other linear equation having a positive decreasing solution can be used as a majorant equation, obtaining different conditions. The first result states sufficient conditions for the existence of a global positive solution of (3.1), bounded on [R,+).

    Theorem 3.2. Assume that (1.2)–(1.5) hold, and let

    Mj(d)=supv(0,d]Fj(v)v,j=1,2

    where d>0 is a fixed constant. If

    B1log2M1(d)P, (3.7)

    and

    p(t)tn,b2(t)(n1)24M2(d)tn2,for alltR, (3.8)

    for some n>1, then for every u0(0,d/2] the equation (3.1) has a solution u, satisfying u(R)=u0,u(R)=0,0<u(t)2u0 for tR.

    Proof. The result follows from [9,Theorem 2 and Corollary 3], with some technical adjustments due to the actual general form of the nonlinearity. Indeed, it is sufficient to observe that, for every continuous function q:[R,)(0,d] fixed, the equations

    (p(t)w(t))M1(d)b1(t)w(t)=0 (3.9)

    is a Sturm minorant of the linearized equation

    (p(t)u(t))+f(t,q(t))q(t)u(t)=0, (3.10)

    and (3.6) is a Sturm majorant, due to (3.8) (see for instance [15]). Since (3.6) is nonoscillatory, and has the solution y(t)=t(n1)/2, which is positive decreasing on [R,+), from [9,Lemma 3] the solution of (3.10) satisfying the initial conditions u(R)=u0, u(R)=0 exists and is positive on [R,+), since it satisfies u(t)w0(t) for all tR, where w0 is the principal solution of (3.9). Further, double integration of (3.10) on [R,t], t>R, gives

    U(t)u0+M1(d)B1tRU(s)p(s)ds, where U(t)=maxs[R,t]u(s),

    and the Gronwall lemma, together with (3.7), gives the upper bound u(t)U(t)2u0. Thus, put

    S={qC1[t0,):q(R)=u0, q(R)=0, q(t)>0 for tR},Ω={qC1[1,):q(R)=u0,q(R)=0,w0(t)q(t)2u0}.

    We have SΩ=Ω and the set S1=¯T(Ω), where T is the operator which maps every qΩ into the unique solution of (3.10) satisfying the initial conditions u(R)=u0, u(R)=0, satisfies S1SΩ=Ω and is bounded in C1[R,+) (for the detailed proof see [9,Theorem 2]). Then Proposition 3.1 can be applied, and the existence of a solution of (3.1) in the set SΩ is proved.

    Remark 3.3. Clearly, the result in Theorem 3.2 holds also if we allow a different upper bound for the solution. More precisely, if we look for a solution satisfying 0<u(t)ku0 with k>1, then it is sufficient to put logk instead of log2 in (3.7), and we get the existence of global positive solutions of the Cauchy problem associated with (3.1), for every u0 such that 0<ku0d.

    Remark 3.4. Since M2(d) is nondecreasing, condition (3.7) can be seen as an upper bound for the values of u0 for which (3.1) has a global bounded solution. For instance, if F1(v)=vβ, β>1, then M1(d)=dβ1 and (3.7) can be read as

    2u0(log2PB1)1β1ifB10,

    while, if F1(v)=v, then M1(d)=1 for all d>0, and (3.1), (3.2) has solution for all u0>0, provided (3.8) is satisfied.

    If in the statement of Theorem 3.2 we assume in addition the condition (1.6), we get the existence of a solution of the BVP (3.1), (3.2).

    Theorem 3.5. Assume that (1.2)–(1.6) hold, and that d>0 exists such that (3.7) and (3.8) are satisfied for some n>1. Then, for every u0(0,d/2], the BVP (3.1), (3.2) has at least a solution u, satisfying

    0<u(t)2u0fort[R,), u(t)<0for larget.

    Proof. The proof is analogous to the second part of the proof of [9,Theorem 2], with obvious modifications due to the more general form of the nonlinear term here considered. In particular notice that (3.1) with conditions (1.3) gives the inequality

    (p(t)u(t))F1(u(t))+b1(t)0

    for every positive solution u of (3.1) and all tR. Thus the arguments in the proof of [9,Theorem 2] apply also to the present case.

    We conclude this Section pointing out that the case b1(t)0 for tR is included in the previous results, and in this case Theorems 3.2, 3.5 have a more simple statement. Indeed, B1=0, and therefore (1.5) and (3.7) are trivially satisfied. Further, every solution of (3.1) is nonincreasing on the whole half-line.

    Combining Theorems 2.4 and 3.2 or 2.4 and 3.5, we obtain an existence result for one or more solutions of the BVP (3.1) and (3.2), respectively. We limit ourself to state results for the existence of one or two solutions, for sake of simplicity. Clearly, as pointed out in Section 2, adding more conditions, with similar arguments it is possible to obtain sufficient conditions for the existence of three solutions (see Theorem 2.4) or more solutions.

    Theorem 4.1. Suppose that (1.2)–(1.5) are satisfied, and that there exist ρ1,ρ2(0,+), with ρ1<ρ2, such that either (S1) or (S2) holds. If n>1 exists, such that (3.7), (3.8) are satisfied with d=2ρ2, then the BVP (1.7) has at least one solution u1, with u1(t)2ρ2 for all t0.

    If, in addition, there exists ρ3(0,+), with ρ2<ρ3, such that either (S3) or (S4) holds, and (3.7), (3.8) are satisfied with d=2ρ3, then the BVP (1.7) has an additional solution u2, with u2(t)2ρ3.

    Notice that, from Theorem 2.4, the solutions u1,u2 satisfy ρ1<u1(R)<ρ2<u2(R)<ρ3 and therefore they are distinct solutions. Further, since solutions on [0,R] are increasing, then u1(t)<ρ2,u2(t)<ρ3, for all t[0,R].

    In case also assumption (1.6) is satisfied, from the above Theorem we obtain sufficient conditions for the existence of solutions of the BVP (1.1).

    Theorem 4.2. Suppose that (1.2)–(1.6) are satisfied, and that there exist ρ1,ρ2(0,+), with ρ1<ρ2, such that either (S1) or (S2) holds. If n>1 exists, such that (3.7), (3.8) are satisfied with d=2ρ2, then the BVP (1.1) has at least one solution u1, with u1(t)2ρ2 for all t0.

    If, in addition, there exists ρ3(0,+), with ρ2<ρ3, such that either (S3) or (S4) holds, and (3.7), (3.8) are satisfied with d=2ρ3, then the BVP (1.1) has an additional solution u2, with u2(t)2ρ3.

    Remark 4.3. If we have that f(t,v)0 for (t,v)[0,+)×[0,+), then every solution of (3.1) satisfying u(R)=0 is nonincreasing for tR. Thus our method provides one (or more) solutions having a global maximum at t=R, as solutions on [0,R] are nondecreasing. Clearly, if f changes sign in some t1>R, f for t\in [R, t_1] , then t = R is only a point of local maximum (not necessarily global) for the solutions.

    We conclude this section with the following examples that illustrate our results.

    Example 4.4. Let us consider the BVP

    \begin{equation} \begin{cases} (p(t)u'(t))'+b(t)u^2(t) = 0,\ t\geq 0,\\ u(0) = H[u],\ u(+\infty) = 0, \end{cases} \end{equation} (4.1)

    where

    p(t) = \begin{cases} 1, &0\leq t\leq 1,\\ t^{n},& t > 1, \end{cases} \quad b(t) = \begin{cases} 1, &0\leq t\leq 1,\\ \sin^+(\frac{\pi}{2}t)- \dfrac{\mu}{t^2}\sin^-(\frac{\pi}{2}t), & t > 1, \end{cases}

    for some n > 1 and \mu > 0 , and

    H[u] = \frac{1}{2}\sqrt{\int_0^1u(t)\, dt}.

    The definition of H leads to the natural choice [0, R] = [0, 1] . By direct computation one has m = 2 . The choice of [a, b] = [1/2, 1] leads to c = 1/2 and M = 4 . Furthermore note that \underline{f}_{\rho} = \frac{1}{4}\rho^2 , \overline{f}_{\rho} = \rho^2 and

    \frac{1}{2}\sqrt{\rho}\geq H[u]\geq \frac{1}{2}\sqrt{\int_{1/2}^1 u(t)\, dt}\geq \frac{1}{4}\sqrt{\rho},\; {for\; every}\; u\in \partial K_\rho.

    Observe that the inequalities

    \begin{align*} \frac{1}{4}\underline{f}_{\rho_1}+\underline{H}_{\rho_1}[u]&\geq \frac{1}{4}\bigl(\frac{1}{4}{\rho_1}^2+\sqrt{\rho_1}\bigr) > \rho_1,\\ \frac{1}{2} \overline{f}_{\rho_2} + \overline{H}_{\rho_2}&\leq \frac{1}{2} \bigl({\rho_2}^2 +\sqrt{\rho_2}\bigr) < \rho_2,\\ \frac{1}{4}\underline{f}_{\rho_3}+\underline{H}_{\rho_3}[u]&\geq \frac{1}{4}\bigl(\frac{1}{4}{\rho_3}^2+\sqrt{\rho_3}\bigr) > \rho_3. \end{align*}

    are satisfied for \rho_1 = 1/20 , \rho_2 = 3/4 and \rho_3 = 15 . Thus ( S_3 ) holds for \rho_1 = 1/20 , \rho_2 = 3/4 and \rho_3 = 15 . Note that (1.2)–(1.5) are satisfied, with P = 1/(n-1) , F_1(v) = F_2(v) = v^2 , b_1(t) = b(t), \, b_2(t) = 1 , and B_1^- < \mu holds. Straightforward calculations show that also (1.6) is satisfied. Since M_1(d) = M_2(d) = d , applying Theorem 4.2 we get the following result:

    If (n-1)^2\geq 6 then (4.1) has at least a positive solution u_1 , with 0 \leq u_1(t)\leq 3/2 , for every \mu \leq 2(n-1) \log 2/3 .

    If (n-1)^2\geq 120 then (4.1) has at least two distinct positive solutions u_1, \, u_2 , with 0 \leq u_1(t)\leq 3/2 , 0 \leq u_2(t)\leq 30 , for every \mu \leq (n-1) \log 2/30 .

    Example 4.5. Let us consider the BVP

    \begin{equation} \begin{cases} (p(t)u'(t))'+b(t)u^2(t) = 0,\ t\geq 0,\\ u(0) = H[u],\ u'(1) = 0, \, u(t)\geq 0 \;{ and\; bounded \;on }\; [0, +\infty), \end{cases} \end{equation} (4.2)

    where this time we take

    p(t) = \begin{cases} 1, &0\leq t\leq 1,\\ t^{5},& t > 1, \end{cases} \quad b(t) = \begin{cases} 1, &0\leq t\leq 1,\\ 2-t, & 1 < t\leq 3,\\ -e^{-(t-3)}, &t > 3, \end{cases}

    and we keep H[u] = \frac{1}{2}\sqrt{\int_0^1 u(t)\, dt} , as in Example 4.3. With the natural choice [0, R] = [0, 1] , following the same computations as in Example 4.3, we get that ( S_1 ) holds for \rho_1 = 1/20 , \rho_2 = 3/4 . Furthermore (1.2)–(1.5) are satisfied, with P = 1/4 , F_1(v) = F_2(v) = v^2 , b_1(t) = b(t), \, b_2(t) = 1 , and B_1^- = 3/2 , while (1.6) does not hold. Since M_1(d) = M_2(d) = d , and (3.7), (3.8) are satisfied for all d\leq \min\{\frac83 \ln 2, 4\} = \frac83 \ln 2 > \frac 32 = 2 \rho_2 , applying Theorem 4.1 we get that (4.2) has at least a solution u_1 , with 0 \leq u_1(t)\leq 3/2 for all t \geq 0 .

    Note that u_1 does not tend to zero as t\to +\infty . Indeed, as u_1'(1) = 0 and b(t) > 0 in [1, 2) , then u_1 is decreasing on [1, 2] . Thus, for all t \in [1, 2] it holds

    t^5 u_1'(t) = -\int_1^t (2-s) u_1^2(s) \, ds\geq -u_1^2(1)\int_1^2 (2-s) \, ds = -u_1^2(1)/2.

    Integration of the above inequality on [1, t], t \leq 2 gives

    u_1(t)\geq u_1(1)-\frac{u_1^2(1)}{2}\int_1^t \frac{1}{s^5}\, ds \geq u_1(1)\left(1-\frac{u_1(1)}{8}\right).

    Thus we have u_1'(2)\geq -u_1^2(1)/64 and u_1(2)\geq u_1(1)\left(1-u_1(1)/8\right) > 0 if u_1(1) < 8 . Now, for t \geq 2 , it holds (t^5 u_1'(t))' = -b(t) u_1^2(t)\geq 0, thus t^5 u_1'(t) \geq 32 u_1'(2)\geq -u_1^2(1)/2 , and integrating this inequality on [2, t], t\geq2 we obtain

    \begin{align*} u_1(t) &\geq u_1(2) - \frac{u_1^2(1)}{2} \int_2^t \frac{1}{s^5}\, ds \geq u_1(1)\left( 1-\frac{u_1(1)}{8}-\frac{u_1(1)}{2} \int_2^\infty \frac{1}{s^5}\, ds\right)\\ & = u_1(1)\left( 1-\frac{17}{128} \, u_1(1)\right). \end{align*}

    Thus, since 1/20 < u_1(1) < 3/4 , from the above inequality we have

    u_1(t)\geq \frac{1}{20}\Big(1-\frac{51}{512}\Big), \quad \forall t \geq 2,

    and therefore u_1 has no zero limit as t \to +\infty .

    The authors wish to thank the Referee for the careful reading of the paper and the constructive comments. G. Infante and S. Matucci were partially supported by G.N.A.M.P.A. - INdAM (Italy).

    The authors declare no conflict of interest in this paper.



    [1] D. S. A. Narayana, K. S. Shankarappa, M. R. Govindappa, H. A. Prameela, M. R. G. Rao, K. T. Rangaswamy, Natural occurrence of Jatropha mosaic virus disease in India, Curr. Sci., 91 (2006), 584–586.
    [2] O. J. Alabi, P. L. Kumar, R. A. Naidu, Cassava mosaic disease: A curse to food security in sub-Saharan Africa, APSnet Features, 2011, 1–17.
    [3] J. E. Duffus, Whitefly transmission of plant viruses, In: Current topics in vector research, New York: Springer, 1987, 73–91. https://doi.org/10.1007/978-1-4612-4712-8_3
    [4] S. Q. Gao, J. Qu, N. H. Chua, J. Ye, A new strain of Indian cassava mosaic virus causes a mosaic disease in the biodiesel crop Jatropha curcas, Arch. Virol., 155 (2010), 607–612. https://doi.org/10.1007/s00705-010-0625-0 doi: 10.1007/s00705-010-0625-0
    [5] J. P. Tewari, H. D. Dwivedi, P. Madhvi, S. K. Srivastava, Incidence of a mosaic disease in Jatropha curcas L. from eastern Uttar Pradesh, Curr. Sci., 93 (2007), 1048–1049.
    [6] B. D. Kashina, M. D. Alegbejo, O. O. Banwo, S. L. Nielsen, M. Nicolaisen, Molecular identification of a new begomovirus associated with mosaic disease of Jatropha curcas L. in Nigeria, Arch. Virol., 158 (2013), 511–514. https://doi.org/10.1007/s00705-012-1512-7 doi: 10.1007/s00705-012-1512-7
    [7] H. Czosnek, A. Hariton-Shalev, I. Sobol, R. Gorovits, M. Ghanim, The incredible journey of begomoviruses in their whitefly vector, Viruses, 9 (2017), 273. https://doi.org/10.3390/v9100273 doi: 10.3390/v9100273
    [8] A. S. Costa, Whitefly-transmitted plant diseases, Annu. Rev. Phytopathol., 14 (1976), 429–449. https://doi.org/10.1146/annurev.py.14.090176.002241 doi: 10.1146/annurev.py.14.090176.002241
    [9] M. Leclerc, T. Doré, C. A. Gilligan, P. Lucas, J. A. N. Filipe, Estimating the delay between host infection and disease (incubation period) and assessing its significance to the epidemiology of plant diseases, PLoS One, 9 (2014), e86568. https://doi.org/10.1371/journal.pone.0086568 doi: 10.1371/journal.pone.0086568
    [10] A. van Maanen, X. M. Xu, Modelling plant disease epidemics, Eur. J. Plant Path., 109 (2003), 669–682. https://doi.org/10.1023/A:1026018005613 doi: 10.1023/A:1026018005613
    [11] B. Buonomo, M. Cerasuolo, The effect of time delay in plant-pathogen interactions with host demography, Math. Biosci. Eng., 12 (2015), 473–490. https://doi.org/10.3934/mbe.2015.12.473 doi: 10.3934/mbe.2015.12.473
    [12] T. Erneux, Applied delay differential equations, New York: Springer, 2009. https://doi.org/10.1007/978-0-387-74372-1
    [13] Y. Kuang, Delay differential equations with applications in population dynamics, New York: Academic Press, 1993.
    [14] J. E. Duffus, Whitefly transmission of plant viruses, In: Current topics in vector research, New York: Springer, (1987) 73–91.
    [15] J. P. Legg, African cassava mosaic disease, In: Encyclopedia of virology, Third Edition, Elsevier, 2008, 30–36.
    [16] F. A. Basir, S. Adhurya, M. Banerjee, E. Venturino, S. Ray, Modelling the effect of incubation and latent periods on the dynamics of vector-borne plant viral diseases, Bull. Math. Biol., 82 (2020), 94. https://doi.org/10.1007/s11538-020-00767-2 doi: 10.1007/s11538-020-00767-2
    [17] J. E. Van der Plank, Plant diseases: Epidemics and control, New York: Academic Press, 2013.
    [18] M. J. Jeger, J. Holt, F. van den Bosch, L.V. Madden, Epidemiology of insect-transmitted plant viruses: Modelling disease dynamics and control interventions, Physiol. Entomol., 29 (2004), 291–304. https://doi.org/10.1111/j.0307-6962.2004.00394.x doi: 10.1111/j.0307-6962.2004.00394.x
    [19] M. Jackson, B. M. Chen-Charpentier, A model of biological control of plant virus propagation with delays, J. Comput. Appl. Math., 330 (2018), 855–865. https://doi.org/10.1016/j.cam.2017.01.005 doi: 10.1016/j.cam.2017.01.005
    [20] M. Jackson, B. M. Chen-Charpentier, Modeling plant virus propagation with delays, J. Comput. Appl. Math., 309 (2017), 611–621. https://doi.org/10.1016/j.cam.2016.04.024 doi: 10.1016/j.cam.2016.04.024
    [21] G. Neofytou, Y. N. Kyrychko, K. B. Blyuss, Time-delayed model of immune response in plants, J.Theor. Biol., 389 (2016), 28–39. https://doi.org/10.1016/j.jtbi.2015.10.020 doi: 10.1016/j.jtbi.2015.10.020
    [22] K. Cooke, P. van den Driessche, X. Zou, Interaction of maturation delay and nonlinear birth in population and epidemic models, J. Math. Biol., 39 (1999), 332–352. https://doi.org/10.1007/s002850050194 doi: 10.1007/s002850050194
    [23] G. A. Ngwa, On the population dynamics of the malaria vector, Bull. Math. Biol., 68 (2006), 2161–2189. https://doi.org/10.1007/s11538-006-9104-x doi: 10.1007/s11538-006-9104-x
    [24] M. Martcheva, O. Prosper, Unstable dynamics of vector-borne diseases: Modeling through delay-differential equations, In: Dynamic models of infectious diseases, New York: Springer, 2013, 43–75. https://doi.org/10.1007/978-1-4614-3961-5_2
    [25] M. Banerjee, Y. Takeuchi, Maturation delay for the predators can enhance stable coexistence for a class of prey-predator models, J. Theor. Biol., 412 (2017), 154–171. https://doi.org/10.1016/j.jtbi.2016.10.016 doi: 10.1016/j.jtbi.2016.10.016
    [26] F. A. Basir, Y. N. Kyrychko, K. B. Blyuss, S. Ray, Effects of vector maturation time on the dynamics of cassava mosaic disease, Bull. Math. Biol., 83 (2021), 87. https://doi.org/10.1007/s11538-021-00921-4 doi: 10.1007/s11538-021-00921-4
    [27] F. A. Basir, K. B. Blyuss, S. Ray, Modelling the effects of awareness-based interventions to control the mosaic disease of Jatropha curcas, Ecol. Complex, 36 (2018), 92–100. https://doi.org/10.1016/j.ecocom.2018.07.004 doi: 10.1016/j.ecocom.2018.07.004
    [28] N. Rakshit, F. A. Basir, A. Banerjee, S. Ray, Dynamics of plant mosaic disease propagation and the usefulness of roguing as an alternative biological control, Ecol. Complex., 38 (2019), 15–23. https://doi.org/10.1016/j.ecocom.2019.01.001 doi: 10.1016/j.ecocom.2019.01.001
    [29] K. B. Blyuss, Control of mosaic disease using microbial biostimulants: Insights from mathematical modelling. Ricerche mat., 69 (2020), 437–455. https://doi.org/10.1007/s11587-020-00508-6 doi: 10.1007/s11587-020-00508-6
    [30] F. A. Basir, S. Ray, Impact of farming awareness based roguing, insecticide spraying and optimal control on the dynamics of mosaic disease, Ricerche mat., 69, (2020), 393–412. https://doi.org/10.1007/s11587-020-00522-8 doi: 10.1007/s11587-020-00522-8
    [31] E. Venturino, P. K. Roy, F. A. Basir, A. Datta, A model for the control of the mosaic virus disease in Jatropha curcas plantations, Energ. Ecol. Environ., 1 (2016), 360–369. https://doi.org/10.1007/s40974-016-0033-8 doi: 10.1007/s40974-016-0033-8
    [32] F. A. Basir, E. Venturino, P. K. Roy, Effects of awareness program for controlling on mosaic disease in Jatropha curcas plantations, Math. Method Appl. Sci., 40 (2017), 2441–2453. https://doi.org/10.1002/mma.4149 doi: 10.1002/mma.4149
    [33] J. Holt, M. J. Jeger, J. M. Thresh, G. W. Otim-Nape, An epidemiological model incorporating vector population dynamics applied to African cassava mosaic virus disease, J. Appl. Ecol., 34 (1997), 793–806. https://doi.org/10.2307/2404924 doi: 10.2307/2404924
    [34] S. R. Chant, Studies on the transmission of cassava mosaic virus by Bemisia spp. (Aleyrodidae), Ann. Appl. Biol., 46 (1958), 210–215. https://doi.org/10.1111/j.1744-7348.1958.tb02198.x doi: 10.1111/j.1744-7348.1958.tb02198.x
    [35] H. D. Thurston, Tropical plant diseases, St. Paul: APS press, 1998.
    [36] D. Fargette, M. Jeger, C. Fauquet, L. D. Fishpool, Analysis of temporal disease progress of African cassava mosaic Virus, Phytopathology, 84 (1994), 191–198.
    [37] C. Fauquet, D. Fargette, African cassava mosaic virus: Etiology, epidemiology, and control, Plant Dis., 74 (1990), 404–411.
    [38] J. K. Hale, Theory of functional differential equations, New York: Springer-Verlag, 1977. https://doi.org/10.1007/978-1-4612-9892-2
    [39] B. Anderson, J. Jackson, M. Sitharam, Descartes' rule of signs revisited, Amer. Math. Mon., 105 (1998), 447–451. https://doi.org/10.2307/3109807 doi: 10.2307/3109807
    [40] S. A. Levin, Descartes rule of signs-how hard can it be? http://sepwww.stanford.edu/oldsep/stew/descartes.pdf
    [41] J. Dushoff, W. Huang, C. Castillo-Chavez, Backwards bifurcations and catastrophe in simple models of fatal diseases, J. Math. Biol., 36 (1998), 227–248. https://doi.org/10.1007/s002850050099 doi: 10.1007/s002850050099
    [42] C. Castillo-Chavez, B. Song, Dynamical models of tuberculosis and their applications, Math. Biosci. Eng., 1 (2004), 361–404. https://doi.org/10.3934/mbe.2004.1.361 doi: 10.3934/mbe.2004.1.361
    [43] B. Buonomo, D. Lacitignola, On the backward bifurcation of a vaccination model with nonlinear incidence, Nonlinear Anal.-Model., 16 (2011), 30–46. https://doi.org/10.15388/NA.16.1.14113 doi: 10.15388/NA.16.1.14113
    [44] A. K. Misra, P. K. Tiwari, P. Chandra, Modeling the control of algal bloom in a lake by applying some external efforts with time delay, Differ. Equ. Dyn. Syst., 29 (2021), 539–568. https://doi.org/10.1007/s12591-017-0383-5 doi: 10.1007/s12591-017-0383-5
    [45] J. M. Heffernan, R. J. Smith, L. M. Wahl, Perspectives on the basic reproductive ratio, J. R. Soc. Interface, 2 (2005), 281–293. https://doi.org/10.1098/rsif.2005.0042 doi: 10.1098/rsif.2005.0042
  • This article has been cited by:

    1. Nabeela Anwar, Aqsa Ghaffar, Muhammad Asif Zahoor Raja, Iftikhar Ahmad, Muhammad Shoaib, Adiqa Kausar Kiani, Stochastic supervised neuro-architecture design for analyzing vector-borne plant virus epidemics with latency and incubation effects, 2024, 139, 2190-5444, 10.1140/epjp/s13360-024-05883-9
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1504) PDF downloads(106) Cited by(1)

Figures and Tables

Figures(7)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog