Next Article in Journal
Approximation of the Constant in a Markov-Type Inequality on a Simplex Using Meta-Heuristics
Next Article in Special Issue
Nontrivial Solutions of Systems of Perturbed Hammerstein Integral Equations with Functional Terms
Previous Article in Journal
A Fuzzy Evaluation Decision Model for the Ratio Operating Performance Index
Previous Article in Special Issue
The Approximate Analytic Solution of the Time-Fractional Black-Scholes Equation with a European Option Based on the Katugampola Fractional Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Permanence for Nonautonomous Differential Systems with Delays in the Linear and Nonlinear Terms

Departamento de Matemática and CMAFCIO, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal
Mathematics 2021, 9(3), 263; https://doi.org/10.3390/math9030263
Submission received: 21 October 2020 / Revised: 21 January 2021 / Accepted: 25 January 2021 / Published: 28 January 2021
(This article belongs to the Special Issue Functional Differential Equations and Applications 2020)

Abstract

:
In this paper, we obtain sufficient conditions for the persistence and permanence of a family of nonautonomous systems of delay differential equations. This family includes structured models from mathematical biology, with either discrete or distributed delays in both the linear and nonlinear terms, and where typically the nonlinear terms are nonmonotone. Applications to systems inspired by mathematical biology models are given.
MSC:
34K12; 34K25; 34K20; 92D25

1. Introduction

In this paper, we investigate the persistence and permanence for a class of multidimensional nonautonomous delay differential equations (DDEs), which includes some structured models used in population dynamics, epidemiology, and other fields.
We start by setting the abstract framework for the DDEs which we deal with in the next sections. For τ 0 , consider the Banach space C : = C ( [ τ , 0 ] ; R n ) with the norm ϕ = max θ [ τ , 0 ] | ϕ ( θ ) | , where | · | is a fixed norm in R n . We shall consider DDEs written in the abstract form:
x ( t ) = L ( t ) x t + f ( t , x t ) , t t 0 ,
where x t C denotes the segment of a solution x ( t ) given by x t ( θ ) = x ( t + θ ) , τ θ 0 , the operator L ( t ) : C R n is linear bounded, for t t 0 , and the nonlinearities are given by a continuous function f : [ t 0 , ) × C [ 0 , ) n . For simplicity, we set t 0 = 0 .
Recently, there has been a renewed interest in questions of persistence and permanence for DDEs. A number of methods has been proposed to tackle different situations, depending on whether the equations are autonomous or not, scalar or multi-dimensional, monotone or nonmonotone. See [1,2,3,4,5,6,7,8,9,10,11] and references therein for an explanation of the models and motivation from real world applications.
Here, the investigation concerning permanence in [5,6] is pursued. In [5] only cooperative systems were considered, whereas in [6], sufficient conditions for the permanence of systems:
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n a i j ( t ) x j ( t ) + k = 1 m i β i k ( t ) h i k ( t , x i ( t τ i k ( t ) ) ) , i = 1 , , n , t 0 ,
were established. Clearly, family (2) is a particular case of (1). In this paper, the more general framework of systems (1) with distributed delays in both L and f is considered. Although many works only consider constant or time-varying discrete delays, it is apparent that in many contexts distributed delays are more appropriate, for instance, the maturation delay for a species is not exactly the same for all individuals, or the incubation time for a disease is within a time frame (e.g., between 3 and 14 days for Covid-19), rather than a fixed value. These situations are better portrayed by distributed delays. In fact, the use of integro-differential equations with distributed (possibly infinite) delays in predator-prey models goes back to the works of Volterra in population dynamics [12], where typically the “memory functions” appear as L 1 integral kernels. PDEs (partial differential equations) models with discrete or distributed delays have also been largely analyzed. More recently, mathematicians have proposed fractional calculus to account for memory effects, in epidemic models as well as other fields – moreover, both time-delays and fractional-order derivatives have been used simultaneously, for different types of memory [13,14].
The criteria for permanence in [2,4,6] and other works demand that all coefficients are bounded. More recently, some authors have relaxed this restriction [8,9,15,16], though they still imposed some boundedness requirements. Here, the boundedness of all coefficients in (1) will not be a priori assumed. We also emphasize that typically the nonlinearites f i ( t , ϕ ) in (1) are not monotone in the second variable – which is the case of Nicholson-type systems, for example. Our results extend and improve some recent conclusions in the literature [2,3,9,16,17,18], which mostly deal with scalar DDEs and/or cooperative n-dimensional models. Special attention will be given to the case where each component f i of f = ( f 1 , , f n ) in (1) depends only on t and on the component i of the solution,
f ( t , ϕ ) = ( f 1 ( t , ϕ 1 ) , , f n ( t , ϕ n ) ) for t 0 , ϕ = ( ϕ 1 , , ϕ n ) C ,
since sharper results will be obtained for such models.
We now introduce some standard notation. In what follows, R + = [ 0 , ) , the matrix I n , or simply I, denotes the n × n identity matrix and 1 = ( 1 , , 1 ) R n . For τ > 0 , the set C + = C ( [ τ , 0 ] ; ( R + ) n ) is the cone of nonnegative functions in C and ≤ the usual partial order generated by C + : ϕ ψ if and only if ψ ϕ C + . A vector v R n is identified in C with the constant function ψ ( θ ) = v for τ θ 0 . In this way, for u , v R n with u v , the interval I = [ u , v ] in C is the set { ϕ C : u ϕ v } . For τ = 0 , we take C = R n , C + = ( R + ) n ; a vector v R n is positive if all its components are positive, and we write v > 0 . We write ϕ < ψ if ψ ( θ ) < ϕ ( θ ) for θ [ τ , 0 ] ; the relations ≥ and > are also defined in the usual way.
For nonlinear DDEs (1) under conditions of existence and uniqueness of solutions, x ( t , σ , ϕ ) denotes the solution of (1) with initial condition x σ = ϕ , for ( σ , ϕ ) R + × C . For models inspired by real world applications, we shall consider:
C 0 + = { ϕ C + : ϕ ( 0 ) > 0 }
as the set of admissible initial conditions. Without loss of generality, we shall restrict the analysis to solutions x ( t , 0 , ϕ ) with ϕ C 0 + , and assume that f is sufficiently regular so that such solutions are defined on R + . If the set C 0 + is (positively) invariant for (1), the notions of (uniform) persistence, permanence, and stability always refer to solutions with initial conditions in C 0 + . In this way, we say that the system is uniformly persistent (in C 0 + ) if there exists a positive uniform lower bound for all solutions with initial conditions in C 0 + , i.e., there is m > 0 such that all solutions x ( t ) = x ( t , 0 , ϕ ) with ϕ C 0 + defined on R + and satisfy x i ( t , 0 , ϕ ) m for t 1 and i = 1 , , n . The system (1) is said to be permanent if there exist positive constants m , M such that all solutions x ( t ) = x ( t , 0 , ϕ ) with ϕ C 0 + are defined on R + and satisfy m x i ( t ) M for t 1 and i = 1 , , n . As usual, the expression t 1 means “for t > 0 sufficiently large". In some contexts, it will be apparent that it is more convenient to consider a proper and invariant subset S of C 0 + as the set of admissible initial conditions for (1) (see e.g., Section 3.3). Here, a DDE x ( t ) = F ( t , x t ) is said to be cooperative if F = ( F 1 , , F n ) satisfies the quasi-monotone condition (Q) in [19] (p. 78): if ϕ , ψ C + and ϕ ψ , then F i ( t , ϕ ) F i ( t , ψ ) for t 0 , whenever ϕ i ( 0 ) = ψ i ( 0 ) for some i.
The remainder of this paper is divided into three sections. In Section 2, we establish sufficient conditions for uniform persistence and permanence for a large family of nonlinear system (1). To illustrate the results, several applications are given in Section 3, such as to generalized Nicholson systems, Mackey–Glass systems, and a competitive chemostat model. Some examples and counter-examples, showing the necessity of some hypotheses, will also be presented in Section 3. The paper ends with a short section of conclusions and open problems.

2. Persistence and Permanence for a Class of Nonautomous DDEs

In this section, we establish explicit and easily verifiable criteria for both the persistence and permanence of systems (1).
Let C : = C ( [ τ , 0 ] ; R n ) with the supremum norm be the phase space. We start with a general nonutonomous linear differential equation in C,
x ( t ) = L ( t ) x t ,
where L : R + L ( C ; R n ) , L ( C ; R n ) is the usual space of bounded linear operators from C to R n equipped with the operator norm, and t L ( t ) ϕ is the Borel measurable for each ϕ , with L ( t ) bounded on R + by a function m ( t ) in L loc 1 ( R + ; R ) .
Assuming the exponential asymptotic stability of (4), the next theorem provides conditions for the dissipativeness and extinction of perturbed nonlinear systems. Its proof is easily deduced from the variation of constant formula [20] and arguments similar to the ones for ODEs (ordinary differential equations), thus it is omitted.
Theorem 1.
Assume that (4) is exponentially asymptotically stable, and consider the perturbed equation:
x ( t ) = L ( t ) x t + f ( t , x t ) , t 0 ,
where f : R + × S R n is continuous and S is a (positively) invariant set for (5).
(i)
If f is bounded, then (5) is dissipative i.e., all solutions of (5) are defined on R + and there exists M > 0 such that any solution x ( t ) of (5) satisfies lim sup t | x ( t ) | M ;
(ii)
If there exists β : R + R + measurable with β ( s ) d s < such that | f ( t , ϕ ) | β ( t ) ϕ , t 1 , then all solutions x ( t ) of (5) satisfy lim sup t x ( t ) = 0 .
For (4), we now suppose that L = ( L 1 , , L n ) is given by:
L i ( t ) ϕ = d i ( t ) ϕ i ( 0 ) + j = 1 n L i j ( t ) ϕ j , t 0 , ϕ = ( ϕ 1 , , ϕ n ) C , i = 1 , , n ,
with d i ( t ) continuous and L i j ( t ) bounded linear functionals. Although it is not relevant for our results, we may assume that L i i ( t ) is non-atomic at zero (see [20] for a definition). For (4), define the n × n matrix-valued functions:
D ( t ) = diag ( d 1 ( t ) , , d n ( t ) ) , A ( t ) = a i j ( t ) ,
where
a i j ( t ) = L i j ( t ) , t 0 , i , j { 1 , , n } .
For (4), the general hypotheses below will be considered:
(H1)
The functions d i : R + R , L i j : R + L ( C ( [ τ , 0 ] ; R ) ; R ) are continuous (for some τ 0 ), i , j = 1 , , n ;
(H2)
There exist a vector v > 0 and a constant δ > 0 such that D ( t ) A ( t ) δ I n v 0 for t 1 .
Instead of (H2), one may assume:
(H2*)
There exist a vector v > 0 and a constant α > 1 such that D ( t ) v α A ( t ) v for t 1 .
Note that the sign of the function d i ( t ) is not prescribed in (H1), however, if (H2) is fulfilled, then d i ( t ) > 0 for t > 0 large. With the notation in (7), e.g., assumption (H2) above translates as: There exist a vector v = ( v 1 , , v n ) > 0 and T 0 , δ > 0 such that d i ( t ) v i j = 1 n a i j ( t ) v j δ v i for all t T , i = 1 , , n .
Next theorem gives some stability results selected from [21].
Theorem 2.
Consider the system (4) under (H1), and assume one of the following sets of conditions:
(i)
(H2) is satisfied and a i j ( t ) are bounded functions on R + for all i , j = 1 , , n ;
(ii)
(H2*) is satisfied and lim inf t d i ( t ) > 0 for i = 1 , , n ;
(iii)
(4) is the ODE system x i ( t ) = d i ( t ) x i ( t ) + j i d i j ( t ) x j ( t ) , 1 i n , and (H2) is satisfied with a i j ( t ) = | d i j ( t ) | . Then, (4) is exponentially asymptotically stable, in other words, there exist k , α > 0 such that:
| x ( t , t 0 , φ ) | k e α ( t t 0 ) φ for all t t 0 0 , φ C .
Proof. 
The result follows from the criteria in [21] (Theorem 3.1). □
Henceforth, we consider delay differential systems written as:
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n L i j ( t ) x j , t + f i ( t , x t ) , t 0 , i = 1 , , n ,
with the linear functionals L i j ( t ) nonnegative (i.e., L i j ( t ) ϕ j 0 for ϕ j 0 ) and f i ( t , ϕ ) continuous and satisfying same requirements formulated below. In general ϕ f i ( t , ϕ ) is not monotone for the order ≤ associated with C + , and therefore (8) is not cooperative. Recall that, by the Riesz representation theorem, the nonnegative bounded functionals L i j ( t ) have a representation:
L i j ( t ) ϕ j = a i j ( t ) τ 0 ϕ j ( s ) d s ν i j ( t , s ) ,
where a i j ( t ) = L i j ( t ) , the functions ν i j ( t , s ) are defined for ( t , s ) R + × [ τ , ) , are continuous in t, left-continuous and nondecreasing in s, and normalized so that τ 0 d s ν i j ( t , s ) = 1 , t 0 . In the case of no delays in (9), then L i j ( t ) x j , t = a i j ( t ) x j ( t ) with a i j ( t ) 0 . Clearly, this framework includes the particular case of DDEs with multiple time discrete delays in both the linear and nonlinear terms.
Systems (8) are sufficiently general to encompass many relevant models from mathematical biology and other fields. In some contexts, they are interpreted as structured models for populations distributed over n different classes or patches, with migration among the patches, where x i ( t ) is the density of the species on class i, a i j ( t ) ( j i ) is the migration coefficient from class j to class i, d i ( t ) the coefficient of instantaneous loss for class i, and f i ( t , ϕ ) is the growth function for class i. DDEs where the delays intervene in the linear terms have deserved the attention of a number of researchers, e.g. as patch structured population or SIS (susceptible-infective-susceptible) multi-strain epidemic models with time delays for the dispersal among patches [9,11,17]. We also refer the reader to [1,2,3,15,19], for real interpretation of the DDEs under consideration and more applications.
In what follows, for ϕ = ( ϕ 1 , , ϕ n ) C + we use the notations:
ϕ i ̲ = min s [ τ , 0 ] ϕ i ( s ) , ϕ ̲ = min 1 i n ϕ i ̲ .
To obtain the uniform persistence of (8), we need to prescribe the behavior of f ( t , ϕ ) when ϕ 0 + . We shall impose that the nonlinearities satisfy the following conditions:
(H3)
The functions f i : R + × C + R + are completely continuous and locally Lipschitzian in the second variable, i { 1 , , n } ;
(H4)
There exist r > 0 , and continuous functions β i : R + ( 0 , ) , h i : R + R + , such that, for i { 1 , , n } :
(i)
h i ( 0 ) = 0 , h i ( x ) > 0 for 0 < x r and there exists ( h i ) ( 0 + ) = 1 ;
(ii)
For any positive solution x ( t ) of (8),
f i ( t , x t ) β i ( t ) h i ( x ̲ t ) if x t ̲ r , for t 1 ;
(H5)
There exists a constant δ > 0 such that [ M ( t ) δ I n ] 1 0 , t 1 , where M ( t ) is the matrix-valued function defined by:
M ( t ) = B ( t ) + A ( t ) D ( t ) , where B ( t ) = diag ( β 1 ( t ) , , β n ( t ) ) , t 0 ,
with D ( t ) , A ( t ) as above and β i ( t ) as in (H4).
Instead of or together with (H5), and with the some notations, we shall often assume:
(H5*)
There exists a constant α > 1 such that B ( t ) 1 α D ( t ) A ( t ) 1 for t 1 .
Some comments about these assumptions are given in the remarks below.
Remark 1.
If the coefficients β i ( t ) are bounded, then (H5) implies (H5*). Indeed, if (H5) holds and there exists M = max 1 i n sup t 0 β i ( t ) (as a matter of fact, it suffices that d i ( t ) j a i j ( t ) M for some M), then (H5*) is satisfied with any constant α ( 0 , 1 + δ / M ] . The converse is also true if β i ( t ) are all bounded from below by a positive constant, since in this case (H5*) implies that (H5) is satisfied with δ ( 0 , ( 1 α 1 ) c ] for c > 0 such that c = min 1 i n inf t 0 β i ( t ) . Similarly, one easily verifies (conf. [21]) that when the coefficients d i ( t ) are bounded from below by a positive constant, then (H2*) implies (H2) and that, if a i j ( t ) are all bounded, (H2) implies (H2*).
In the study of stability for nonautonous DDEs, a condition as (H2*) with v = 1 has been often presented (see e.g., [5,9]) in the equivalent form (for a positive denominator) lim inf t d i ( t ) j = 1 n a i j ( t ) > 1 . Analogously, (H5*) can be written as lim inf t β i ( t ) d i ( t ) j = 1 n a i j ( t ) > 1 , if d i ( t ) j = 1 n a i j ( t ) > 0 for t 1 .
Remark 2.
If (10) holds with a function h i satisfying h i ( 0 ) = 0 , h i ( x ) > 0 on ( 0 , ) and ( h i ) ( 0 + ) = c i > 0 , by replacing h i ( x ) , β i ( t ) by h ¯ i ( x ) = c i 1 h i ( x ) , β ¯ i ( t ) = c i β i ( t ) , respectively, we may always assume that ( h i ) ( 0 + ) = 1 .
The main criterion for the uniform persistence and permanence of (8) is now established. Although its proof follows along the main ideas in [6] (Theorem 3.3), new arguments are used to take into account the more general form of (8): Namely, delays are allowed in the linear part, the nonlinearity f = ( f 1 , , f n ) need not have the form (3), the coefficients d i ( t ) , β i ( t ) are not required to be bounded below or above by positive constants, nor a i j ( t ) if there are no delays in L i j ( t ) , and d i ( t ) may actually change sign.
Moreover, we should mention that there was an incomplete argument in [6], finished here, since the case of a solution x ( t ) with min 1 j n lim inf t x j ( t ) = lim t x i ( t ) for some i with x i ( t ) strictly increasing for large t > 0 , was not addressed in [6]. See Step 4 of the proof below and Corollary 1, for the treatment of this situation.
Theorem 3.
For (8), assume (H1), (H3), and (H4). Furthermore, let the following conditions hold:
(i)
Either L i j ( t ) ϕ j = a i j ( t ) ϕ j ( 0 ) with a i j ( t ) 0 , for all i , j = 1 , , n and t 0 (in other words, there are no delays in (4)), or L i j ( t ) are nonnegative and a i j ( t ) = L i j ( t ) are bounded on R + , i , j = 1 , , n ;
(ii)
(H5) and (H5*) are both satisfied.
Then (8) is persistent (in C 0 + ). If in addition (H2) holds and f ( t , ϕ ) is bounded, the system is permanent.
Proof. 
The proof follows in several steps.
Step 1. Write (8) as x ( t ) = F ( t , x t ) , where F = ( F 1 , , F n ) ,
F i ( t , ϕ ) = d i ( t ) ϕ i ( 0 ) + j = 1 n L i j ( t ) ϕ j + f i ( t , ϕ ) ( 1 i n ) .
It is clear that F is continuous, locally Lipschitzian in the second variable and bounded on bounded sets of R + × C + . Observe that the solutions of (8) satisfy the ordinary differential inequalities x i ( t ) d i ( t ) x i ( t ) ( 1 i n ) , thus the solutions x ( t , 0 , ϕ ) with ϕ C 0 + are positive for t 0 .
From (ii), there are δ 0 > 0 , α 0 > 1 and T 0 > 0 such that:
β i ( t ) d i ( t ) j a i j ( t ) + δ 0 , β i ( t ) α 0 d i ( t ) j a i j ( t ) , t T 0 , i = 1 , , n .
Summing up these inequalities, we obtain:
α 1 β i ( t ) d i ( t ) j a i j ( t ) + δ t T 0 , i = 1 , , n ,
for α = ( 1 + α 0 ) / 2 > 1 , δ 0 = δ / 2 .
Consider a positive solution x ( t ) = x ( t , 0 , ϕ ) of (8). For r > 0 as in (H4) choose m ( 0 , r ] such that h i ( x ) is strictly increasing with α 1 x < h i ( x ) on the interval ( 0 , m ] . In this way, f i ( t , x t ) β i ( t ) h i x t ̲ α 1 β i ( t ) x t ̲ if x t ̲ m . Replacing h i ( x ) by the function min { h i ( x ) , x } , we may also assume that h i ( x ) x for all x 0 .
We now derive the uniform persistence of (8) by showing in the next steps that, for any solution x ( t ) = x ( t , T 0 , ϕ ) , there exists T T 0 such that:
x i ( t ) m for all t T , 1 i n .
Step 2. We first prove that the ordered interval [ m , ) n : = { ϕ C : ϕ i m , i = 1 , , n } C is invariant for (8) for t [ T 0 , ) .
Note that the operators L i j ( t ) are nondecreasing and h i ( x ) > α 1 x on ( 0 , m ] . If ϕ [ m , ) n and ϕ i ( 0 ) = m for some i, from (13) we therefore obtain, for t T 0 ,
F i ( t , ϕ ) m d i ( t ) + j a i j ( t ) + β i ( t ) h i ( m ) > m d i ( t ) + j a i j ( t ) + α 1 β i ( t ) 0 .
From [19] (Remark 5.2.1), it follows that the set [ m , ) n C is positively invariant for (8).
Step 3. For T 0 as before, define:
s 0 = min j min t [ T 0 , T 0 + τ ] x j ( t ) , s 1 = min j min t [ T 0 + τ , T 0 + 2 τ ] x j ( t ) .
Let s 1 = x i ( t 0 ) , for some t 0 [ T 0 + τ , T 0 + 2 τ ] and i { 1 , , n } .
We first show that s 1 < m implies that s 1 > s 0 .
If s 1 s 0 , then:
s 1 = x i ( t 0 ) = min j min t [ T 0 , T 0 + 2 τ ] x j ( t )
and x i ( t 0 ) 0 . Assuming that s 1 s 0 and x i ( t 0 ) < m , since x j , t x i ( t 0 ) ( 1 j n ) for t [ t 0 τ , t 0 ] , we get:
0 x i ( t 0 ) d i ( t 0 ) x i ( t 0 ) + j a i j ( t 0 ) x j , t 0 ̲ + f i ( t 0 , x t 0 ) d i ( t 0 ) + j a i j ( t 0 ) x i ( t 0 ) + β i ( t 0 ) h i ( x i ( t 0 ) ) > d i ( t 0 ) + j a i j ( t 0 ) + α 1 β i ( t 0 ) x i ( t 0 ) 0 ,
which is a contradiction. This shows that s 1 > s 0 whenever s 1 < m .
Step 4. Define the sequence:
s k = min j min t [ T 0 + k τ , T 0 + ( k + 1 ) τ ] x j ( t ) , k N 0 .
For the sake of contradiction, assume that s k < m for all k N 0 . Thus, reasoning as in Step 3, ( s k ) is strictly increasing. Let t k I k : = [ T 0 + k τ , T 0 + ( k + 1 ) τ ] be such that s k = x i k ( t k ) , for some i k { 1 , , n } . By jumping some of the intervals I k and considering a subsequence of ( t k ) , still denoted by ( t k ) , we may consider a unique i { 1 , , n } such that s k = x i ( t k ) . Denote = lim k s k = sup k s k > 0 . Clearly, = lim inf t x i ( t ) = min 1 j n lim inf t x j ( t ) .
If = lim t x i ( t ) , then x i ( t ) is bounded and from the fluctuation lemma (taking a subsequence if necessary) we have x i ( t k ) 0 as k . If x i ( t ) does not converge to , x i ( t k ) are local minima for k large, therefore x i ( t k ) = 0 . In both cases, we have x i ( t k ) 0 as k .
Let α > 1 and δ > 0 be as in (14). We now claim that:
s k + 1 α min j h j ( s k ) , k N 0 .
Otherwise, suppose that there is k such that s k + 1 < α min j h j ( s k ) .
We distinguish two situations: Either there are no delays in (4), or a i j ( t ) are all bounded in R + – in which case we suppose that α is chosen so that it also satisfies 1 < α < M 1 M 1 δ , where max 1 j n sup t 0 ( j a i j ( t ) ) M 1 .
First, we treat the case of no delays in the linear part L of (4). In this situation, L i j ( t ) x j , t = a i j x j ( t ) a i j ( t ) s k for t I k ( k N 0 ) . Estimate (14) leads to:
x i ( t k + 1 ) d i ( t k + 1 ) s k + 1 + j a i j ( t k + 1 ) s k + 1 + β i ( t k + 1 ) h i ( s k ) > d i ( t k + 1 ) + j a i j ( t k + 1 ) + α 1 β i ( t k + 1 ) s k + 1 δ s k + 1 δ > 0 ,
which is not possible. Thus, (16) holds.
When delays are allowed in the linear part, we can write L i j x j , t a i j ( t ) s k for t I k + 1 , and s k h i ( s k ) > α 1 s k + 1 , thus we have:
x i ( t k + 1 ) d i ( t k + 1 ) s k + 1 + j a i j ( t k + 1 ) s k + β i ( t k + 1 ) h i ( s k ) > α 1 β i ( t k + 1 ) + j a i j ( t k + 1 ) d i ( t k + 1 ) s k + 1 + ( s k s k + 1 ) j a i j ( t k + 1 ) δ s k + 1 + s k + 1 ( α 1 1 ) M 1 = [ δ + ( α 1 1 ) M 1 ] s k + 1 [ δ + ( α 1 1 ) M 1 ] > 0 ,
which is a contradiction. Thus, claim (16) is proven.
From (16), we obtain m α min j h j ( ) > , which is not possible. Therefore, s k m for some k, and the result follows by Step 2.
If in addition (H2) is satisfied, from Theorem 2 the linear system (4) is exponentially asymptotically stable (for the case of no delays in the linear functionals L i j ( t ) , recall that the boundedness of a i j ( t ) is not required). With f bounded, Theorem 1 implies now that (8) is dissipative, and therefore permanent. □
When f = ( f 1 , , f n ) has the form in (3), not only may one take x i , t ̲ instead of x t ̲ in (H4)(ii), but a slightly stronger version of the above theorem holds.
Theorem 4.
For
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n L i j ( t ) x j , t + f i ( t , x i , t ) , t 0 , i = 1 , , n ,
the assertions in Theorem 3 are valid with (H4)(ii) replaced by:
f i ( t , x i , t ) β i ( t ) h i x i , t ̲ if x i , t ̲ r , for t 1 .
and (H5), (H5*) replaced by: There exist a positive vector v and:
(h5)
A constant δ > 0 such that [ M ( t ) δ I n ] v 0 , t 1 ;
(h5*)
A constant α > 1 such that B ( t ) v α D ( t ) A ( t ) v for t 1 .
Proof. 
For (18) under the above hypotheses, rescaling the variables by x ^ i ( t ) = v i 1 x i ( t ) ( 1 i n ) , where v = ( v 1 , , v n ) > 0 is a vector satisfying simultaneously (h5) and (h5*), we obtain a new system:
x ^ i ( t ) = d i ( t ) x ^ i ( t ) + j = 1 n L ^ i j ( t ) x ^ j , t + f ^ i ( t , x ^ i , t ) , i = 1 , , n , t 0 ,
where a ^ i j ( t ) : = L ^ i j ( t ) = v i 1 a i j ( t ) v j and f ^ i ( t , ϕ i ) : = v i 1 f i ( t , v i ϕ i ) . Hence, (H4) is satisfied with h i ( x ) in (19) replaced by h ^ i ( x ) = v i 1 h i ( v i x ) . In this way, and after dropping the hats for simplicity, we may consider an original system (8) and take v = 1 in (h5), (h5*). The result follows by the proof of Theorem 3. □
In the case of bounded nonlinearities, Theorem 1 shows that the permanence in Theorem 3 is still obtained if one replaces (H2) by the requirement of (4) being exponentially asymptotically stable.
Remark 3.
When the linearities do not have delays, the proof of Theorem 3 requires the use of assumptions (H5*) and (H5) in (17), but not the boundedness of the coefficients a i j ( t ) . Moreover, as explained below, (H5) is not needed at all, unless the solution has an eventually increasing component.
These observations and Theorem 2 (iii) allow us to conclude the following:
Corollary 1.
For
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n a i j ( t ) x j ( t ) + f i ( t , x t ) , i = 1 , , n ,
with d i : R + ( 0 , ) , a i j : R + R + continuous for all i , j , assume (H3), (H4). If
lim inf t d i ( t ) j = 1 n a i j ( t ) > 0 , lim inf t β i ( t ) d i ( t ) j = 1 n a i j ( t ) > 1 for i = 1 , , n ,
then (21) is persistent. If in addition f = ( f 1 , , f n ) is bounded, then (21) is permanent.
Proof. 
Condition (22) implies that (H5*) is satisfied, as well as (H2) with v = 1 , hence we follow up the proof of Theorem 3.
Let α 0 > 1 be such that B ( t ) 1 α 0 D ( t ) A ( t ) 1 for t > 0 large. Take some α ( 0 , α 0 ) in Step 4 of the aforementioned proof. From [11] (Lemma A6), for x i ( t ) , I k , and ( t k ) as in Step 4 of this proof, either x i ( t k ) 0 or x i ( t ) is eventually increasing. Observe that when there are no delays in the linear terms, assumption (H5) was just used in (17), to rule out the situation of x i ( t ) eventually increasing, since in this case, the minima s k = x ( t k ) are attained at t k = min I k and x i ( t k ) 0 .
In fact, if x i ( t k ) 0 , from (17) we get:
x i ( t k + 1 ) > d i ( t k + 1 ) + j a i j ( t k + 1 ) + α 1 β i ( t k + 1 ) s k + 1 0 ,
thus a contradiction. Now, suppose that x i ( t ) is eventually increasing. We just have to reach a contradiction without using (H5). Recall that, in this situation, = lim t x i ( t ) = sup k x i ( t k ) > 0 and x i ( t k ) 0 + . For k large, as in (17) we have:
x i ( t k ) > d i ( t k ) + j a i j ( t k ) + α 1 β i ( t k ) x i ( t k ) = d i ( t k ) + j a i j ( t k ) + α 0 1 β i ( t k ) x i ( t k ) + ( α 1 α 0 1 ) β i ( t k ) x i ( t k ) ( α 1 α 0 1 ) β i ( t k ) x i ( t k ) > 0 .
If lim inf k β i ( t k ) = : b > 0 , we obtain 0 > b ( α 1 α 0 1 ) , a contradiction. Otherwise, for a subsequence of β i ( t k ) (still denoted by β i ( t k ) ) such that β i ( t k ) 0 , by taking limits in the above inequalities we obtain:
0 = lim k d i ( t k ) + j a i j ( t k ) + α 1 β i ( t k ) = lim k d i ( t k ) + j a i j ( t k ) .
Again, this is not possible, since d i ( t ) j a i j ( t ) δ for t 1 and some δ > 0 .
Therefore, we conclude that (21) is uniformly persistent. If f is bounded, the permanence follows from Theorems 1 and 2 (iii). □
It is apparent that a scaling, as affected in Theorem 4, allows us to consider a general positive vector v in (22) if each component f i of f ( t , x t ) depends only on t and x i , t .
Corollary 2.
Consider the DDE:
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n a i j ( t ) x j ( t ) + f i ( t , x i , t ) , i = 1 , , n ,
where d i : R + ( 0 , ) , a i j : R + R + are continuous for all i , j and f = ( f 1 , , f n ) satisfies (H3). Suppose that (H4) is satisfied with (19) instead of (10) and that there exists a positive vector v = ( v 1 , , v n ) such that:
lim inf t d i ( t ) v i j = 1 n a i j ( t ) v j > 0 , lim inf t β i ( t ) v i d i ( t ) v i j = 1 n a i j ( t ) v j > 1 for i = 1 , , n .
Then (24) is persistent. Moreover, if f is bounded, then (24) is permanent.
If the nonlinear terms f i in (8) are not bounded but are sublinear at infinity, a condition stronger than (H2) still gives the dissipativeness of the system.
Theorem 5.
Consider (8), under (H1), (H3). Suppose that there exist functions β i + , h i + : R + R + and a constant R 0 > 0 such that, for ϕ C 0 + and i = 1 , , n :
(i)
0 f i ( t , ϕ ) β i + ( t ) h i + ( ϕ ) for t 1 and ϕ R 0 ;
(ii)
h i + ( x ) x for x R ;
(iii)
N ( t ) 1 0 for t 1 , where N ( t ) = diag ( d 1 ( t ) β 1 + ( t ) , , d n ( t ) β n + ( t ) ) A ( t ) .Then all positive solutions are bounded. Moreover, if:
(iv)
lim inf t β i + ( t ) > 0 and h i + ( x ) < x for x R 0 , then (8) is dissipative.
Proof. 
From (i) and (iii), take T 0 τ such that 0 f i ( t , ϕ ) β i + ( t ) h i + ( ϕ ) and d i ( t ) j a i j ( t ) β i + ( t ) 0 , for ϕ C 0 + with ϕ R 0 and t T 0 , i = 1 , , n .
Consider any ordered interval [ 0 , R ] n in C + , with R R 0 . For F i defined in (12), if t T 0 and ϕ [ 0 , R ] n with ϕ i ( 0 ) = R , we have
F i ( t , ϕ ) = d i ( t ) R + j L i j ( t ) ϕ j + f i ( t , ϕ ) d i ( t ) j a i j ( t ) β i + ( t ) R 0 .
Hence, from [19] (Remark 5.2.1) the interval [ 0 , R ] n is positively invariant for t T 0 . In particular, all solutions are bounded.
Next, suppose that (iv) holds. Let x ( t ) be a solution with initial condition in C 0 + . We claim that lim sup t | x ( t ) | R 0 . Otherwise, there is i and a sequence t k such that x t k = x i ( t k ) > R 0 and x i ( t k ) 0 . This implies:
x i ( t k ) d i ( t k ) x i ( t k ) + j a i j ( t k ) x j , t k + β i + ( t k ) h i + ( x t k ) < d i ( t k ) x i ( t k ) + j a i j ( t k ) x i ( t k ) + β i + ( t k ) h i + ( x i ( t k ) ) = [ d i ( t k ) j a i j ( t k ) β i + ( t k ) ] x i ( t k ) + β i + ( t k ) h i + ( x i ( t k ) ) x i ( t k ) β i + ( t k ) h i + ( x i ( t k ) ) x i ( t k ) .
As x i ( t k ) 0 , taking a subsequence of ( x i ( t k ) ) if necessary, we obtain 0 lim inf t β i + ( t ) h i + ( M ) M < 0 , where M = lim k x i ( t k ) , a contradiction. □
Remark 4.
Once again, in the case of (18), after a scaling one can replace the unit vector 1 in (iii) of the above theorem by some vector u > 0 . On the other hand, the above proof shows that Theorem 5 is still true with the assumption [ N ( t ) δ I n ] 1 0 for some δ > 0 and t 1 , instead of (iv).
We end this section with two remarks, leading to more precise and general results.
Remark 5.
More explicitly, we could have written the linear DDE (4) as:
x i ( t ) = d i ( t ) x i ( t ) + j i d i j ( t ) x j ( t ) + j = 1 n a i j ( t ) τ 0 x j ( t + s ) d s ν i j ( t , s ) , 1 i n ,
with d i j ( t ) continuous ( j i ) and d i ( t ) , a i j ( t ) , ν i j ( t , s ) as above, with s ν i j ( t , s ) non atomic at zero, and apply more precise criteria for its exponential asymptotic stability, see [21]. Namely, the criteria in Theorem 2 hold with the matrix D ( t ) = diag ( d 1 ( t ) , , d n ( t ) ) replaced by D ˜ ( t ) = [ d ˜ i j ( t ) ] , where d ˜ i ( t ) = d i ( t ) and d ˜ i j ( t ) = | d i j ( t ) | for j i . Naturally, in this case, the condition lim inf t d i ( t ) > 0 in (ii) of Theorem 2 should be replaced by lim inf t ( d i ( t ) v i j | d i j ( t ) | v j ) > 0 , for all i. This means that the criterion for permanence in Theorem 3 remains valid with these changes.
Remark 6.
Consider nonlinearites which also incorporate a strictly sublinear negative feedback term of the form K i ( t , x i ( t ) ) , so that (8) reads as:
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n L i j ( t ) x j , t + f i ( t , x t ) K i ( t , x i ( t ) ) , t 0 , i = 1 , , n ,
where K i ( t , x ) 0 are continuous and K i ( t , x ) κ i ( t ) g i ( x ) for some continuous functions κ i , g i : R + R + with κ i ( t ) bounded, g i ( 0 ) = 0 and with right-hand derivative ( g i ) ( 0 + ) = 0 , and (H1), (H3) hold. With f i bounded functions, solutions of (26) satisfy the inequalities d i ( t ) x i ( t ) K i ( t , x i ( t ) ) x i ( t ) d i ( t ) x i ( t ) + j = 1 n L i j ( t ) x j , t + C ( 1 i n ) , where C > 0 is such that f i ( t , ϕ ) C on R + × C + for all i. By comparing below and above the solutions of (26) with solutions of cooperative systems [19] and from Theorem 1, it follows that (26) is dissipative and that C 0 + is forward invariant for (26). On the other hand, for any fixed ε > 0 small, there is m 0 > 0 such that 0 K i ( t , x ) ε x for x [ 0 , m 0 ] . A careful analysis shows that the arguments in the proof of Theorem 3 carry over to (26) if one chooses ε ( 0 , δ ) , for δ > 0 as in (H5), so that (14) is satisfied with d i ( t ) replaced by d i ( t ) + ε . In this way, one may conclude that the permanence results stated in Theorems 3, 4 and Corollaries 1, 2 are still valid for (26). This more general framework allows in particular to consider structured models with harvesting.

3. Applications

We now apply the above results of persistence and permanence to some selected models.

3.1. Nicholson and Mackey–Glass-Type Systems

We first apply our results to generalized Nicholson and Mackey–Glass systems. The literature is recent, though extensive, and here we mention a few selected references dealing with the persistence and permanence for multidimensional Nicholson equations [3,6,7,10,18] and Mackey–Glass equations [2,6,15], and references therein. Meaningful applications were given in the early work of Berezansky et al. [3], where planar Nicholson systems were proposed as a compartment model for leukemia (with a delay representing the time required for the response of leukemia cells to growth) and a marine model with fishing and protected areas (with a maturation delay for the species).
Consider systems given by:
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n L i j ( t ) x j , t + k = 1 m i b i k ( t ) t τ i k ( t ) t λ i k ( s ) g i k ( s , x i ( s ) ) d s , i = 1 , , n ,
where L i j ( t ) x j , t = a i j ( t ) σ i j ( t ) 0 x j ( t + s ) d s ν i j ( t , s ) are as in (9) with σ i j ( t ) 0 d s ν i j ( t , s ) = 1 and either:
g i k ( t , x ) = x e c i k ( t ) x
or
g i k ( t , x ) = x 1 + c i k ( t ) x α i ( α i 1 ) .
The functions d i ( t ) , a i j ( t ) , b i k ( t ) , σ i j ( t ) , τ i k ( t ) , λ i k ( t ) , c i k ( t ) are assumed to be continuous and nonnegative, with σ i j ( t ) , τ i k ( t ) [ 0 , τ ] (for some τ > 0 ), and d i ( t ) > 0 , c i k ( t ) c 0 > 0 , for all i , j , k , and t 0 . For g i k as in (28) a modified Nicholson-type system is obtained, whereas the choice (29) provides a Mackey–Glass-type system. We suppose that the linear operators L i j ( t ) are nonnegative, thus a i j ( t ) = L i j ( t ) as before, and define:
β i ( t ) : = k = 1 m i b i k ( t ) t τ i k ( t ) t λ i k ( s ) d s > 0 , t 0 i = 1 , , n .
From Theorem 4, we derive sufficient conditions for the permanence of (27).
Theorem 6.
For (27) under the general conditions above, let the matrices D ( t ) , A ( t ) , M ( t ) be as in (7) and (11), and assume that:
(i)
The functions a i j ( t ) , c i k ( t ) , β i ( t ) are bounded on R + for all i , j = 1 , , n , k = 1 , , m i ;
(ii)
There are positive vectors u , v , and δ > 0 such that [ D ( t ) A ( t ) δ I ] u 0 , [ M ( t ) δ I ] v 0 .
Then system (27) is permanent.
Proof. 
Clearly, system (27) has the form (18), with f i ( t , ϕ i ) = k = 1 m i b i k ( t ) τ i k ( t ) 0 λ i k ( t + s ) g i k ( t + s , ϕ i ( s ) ) d s , for which (H1), (H3) hold. Condition (ii) expresses that (H2), (h5) are also satisfied. Since g i k and β i ( t ) are bounded, from Theorem 1 the system is dissipative. Note also that (h5) implies (h5*) (conf. Remark 1).
Choose M > 0 such that 0 < x i ( t ) M , t R + , i = 1 , , n , for any solution x ( t ) with an initial condition in C 0 + . Let 0 < c 0 c i k ( t ) c i ¯ for t 0 and all i , k . Define h i ( x ) : = x e c i ¯ x , respectively h i ( x ) : = x 1 + c i ¯ x α i , for Nicholson, and respectively Machey–Glass systems. We always have h i ( 0 ) = 1 . Moreover, the functions h i are unimodal with h i ( ) = 0 , with the exception of the increasing and bounded function x 1 + c i ¯ x (when α i = 1 in (29)). Now, choose m > 0 such that all the functions h i ( x ) are increasing in [ 0 , m ] and: h i ( m ) h i ( M ) , and define h i ( x ) : = h i ( x ) , 0 x m h i ( m ) , x > m . Thus, (19) is satisfied with β i ( t ) as above and these choices of h i ( x ) . The result follows from Theorem 4. □
Remark 7.
As mentioned previously, Theorems 3, 4, and 6 are still valid if one replaces (H2) by the assumptions (H2*) and lim inf t d i ( t ) > 0 , for all i.
For the situation without delays in the linear part, the boundedness of a i j ( t ) is not required:
Corollary 3.
For g i k as in (28) or (29), consider the system:
x i ( t ) = d i ( t ) x i ( t ) + j = 1 n a i j ( t ) x j ( t ) + k = 1 m i b i k ( t ) t τ i k ( t ) t λ i k ( s ) g i k ( s , x i ( s ) ) d s , i = 1 , , n ,
under the above conditions on the coefficients and delays, and assume that:
(i)
The functions c i k ( t ) , β i ( t ) are bounded on R + for all i = 1 , , n , k = 1 , , m i ;
(ii)
There are positive vectors u , v , and δ > 0 such that [ D ( t ) A ( t ) δ I ] u 0 , [ M ( t ) δ I ] v 0 .
Then (30) is permanent.
We emphasize that this corollary gives a sharper criterion for permanence than the one in [6], and moreover applies to a much larger family of delayed structured models.
When β i ( t ) are not bounded, the results in the previous section allow us to derive criteria for persistence. This is illustrated in the example below.
Example 1.
Consider the Mackey–Glass-type system:
x i ( t ) = d i i t η x i ( t ) + j = 1 , j i n d i j t η x j ( t ) + j = 1 n b i j t η τ i j ( t ) 0 x j ( t + s ) d s + β i ( t ) σ i ( t ) 0 k i ( s ) x i ( t + s ) 1 + c i ( t ) x i ( t + s ) ν i d s , 1 i n ,
where η > 0 , 0 < ν i 1 , b i j , d i j R + with d i : = d i i > 0 for all i, and the delays τ i j ( t ) , σ i ( t ) are continuous with 0 τ i j ( t ) r i j , 0 < σ i ( t ) R i for some constants r i j , R i > 0 , β i , c i : R + ( 0 , ) are continuous with c i ( t ) bounded below and above by positive constants, and k i : [ R i , 0 ] R + are integrable with σ i ( t ) 0 k i ( s ) d s = 1 , t 0 , i , j = 1 , , n . Choose h i ( x ) = x x + c ¯ i x ν i for x in some vicinity [ 0 , m ] of 0 and c ¯ i > 0 such that c i ( t ) c ¯ i ( 1 i n ) .
With the previous notation, we have d i ( t ) = d i t η , a i j ( t ) = ( ( 1 δ i j ) d i j + b i j τ i j ( t ) ) t η [ ( 1 δ i j ) d i j + b i j r i j ] t η , where δ i j = 1 if i = j and δ i j = 0 if i j . Define the n × n matrices D = diag ( d 1 , , d n ) , A = ( 1 δ i j ) d i j + b i j r i j , so that D ( t ) = t η D , A ( t ) = t η A . Assume that N : = D A is a non-singular M-matrix, or, in other words, that there exists a positive vector v such that u : = N v > 0 . For δ > 0 small such that δ A v u , we have D v ( 1 + δ ) A v , hence (H2*) is satisfied. From Theorem 2 (ii) we deduce that the linear system:
x i ( t ) = d i i t η x i ( t ) + j = 1 , j i n d i j t η x j ( t ) + j = 1 n b i j t η τ i j ( t ) 0 x j ( t + s ) d s , i = 1 , , n ,
is exponentially asymptotically stable.
Next, suppose that b i j = 0 for all i , j . If α > 1 such that β i ( t ) v i α t η u i for t 1 , then also β i ( t ) v i t η u i + δ , for some δ > 0 and for t 1 , 1 i n . Under this condition, both (h5) and (h5*) hold. From Theorem 4 we conclude that (31) is uniformly persistent. Observe however that none of the coefficients d i ( t ) , a i j ( t ) , β i ( t ) is bounded.

3.2. A Structured Population Model

In [22], Arino et al. proposed an alternative model to the usual delayed logistic equation, given by:
x ( t ) = γ μ e μ τ x ( t τ ) μ + κ ( 1 e μ τ ) x ( t τ ) μ x ( t ) κ x 2 ( t ) ,
where γ , μ , κ , τ > 0 , and addressed the existence and stability of a positive carrying capacity. In this model, the growth rate at time t depends on the number of individuals over the time interval [ t τ , t ] that have survived until time t and the growth delay τ is consistent with the decay rate. It was assumed in [22] that τ < τ H , where τ H : = 1 μ log γ μ > 0 is a threshold time for survival: If τ < τ H , x ( t ) approaches the carrying capacity as t , otherwise the species dies out. Recently, Lin et al. [23] argued that a distributed delay should be considered to account for different delays for different individuals of the population, and derived a generalization of the modified logistic law (32) as follows:
x ( t ) = γ 0 μ e μ s x ( t s ) k ( s ) μ + κ ( 1 e μ s ) x ( t s ) d s μ x ( t ) κ x 2 ( t ) ,
where γ , μ , κ > 0 , the distribution kernel k ( s ) 0 is normalized so that 0 k ( s ) d s = 1 and satisfies τ : = 0 s k ( s ) d s < , where τ is called the mean delay. See also in [Section 4] of [23], for a two species competition model. Meanwhile, a general nonautonomous n-dimensional version of the model in (32), which can be interpreted as a structured population model, was considered and studied in [5,9]:
x i ( t ) = k = 1 m i α i k ( t ) x i ( t τ i k ( t ) ) 1 + c i k ( t ) x i ( t τ i k ( t ) ) + j = 1 n a i j ( t ) x j ( t σ i j ( t ) ) d i ( t ) x i ( t ) κ i ( t ) x i 2 ( t ) , t 0 , i = 1 , , n
where all the coefficients α i k ( t ) , c i k ( t ) , a i j ( t ) , μ i ( t ) , κ i ( t ) , and the delays τ i k ( t ) , σ i j ( t ) are nonnegative and continuous, 0 < c i k ̲ c i k ( t ) c i k ¯ for t 0 , k = 1 , , m i , i , j = 1 , , n . Contrarily to the setting in [9], the delays τ i k ( t ) , σ i j ( t ) are not assumed to be bounded in [5]. See [9,22,23] and references therein, for a derivation of the models and further biological explanations.
Motivated by the aforementioned works, consider now a system with bounded distributed delays:
x i ( t ) = k = 1 m i α i k 0 ( t ) t τ i k ( t ) t x i ( s ) K i k ( t s ) 1 + c i k ( t s ) x i ( s ) d s + j = 1 n a i j 0 ( t ) t σ i j ( t ) t x j ( s ) G i j ( t s ) d s d i ( t ) x i ( t ) κ i ( t ) x i 2 ( t ) , t 0 , i = 1 , , n ,
where α i k 0 , c i k , a i j 0 : R + R + , τ i k , σ i j : R + [ 0 , τ ] ( τ > 0 ) are continuous, c i k ( t ) are bounded below and above by positive constants, the kernels K i k , G i j : [ 0 , τ ] R + are integrable, with:
α i k ( t ) : = α i k 0 ( t ) 0 τ i k ( t ) K i k ( s ) d s , a i j ( t ) : = a i j 0 ( t ) 0 σ i j ( t ) G i j ( s ) d s
continuous, for all i , j , k . Denote also β i ( t ) : = k = 1 m i α i k ( t ) , i = 1 , , n .
System (35) has the form (26), with L i j ( t ) nonnegative linear operators, a i j ( t ) = L i j ( t ) , K i ( t , x ) = κ i ( t ) x 2 and:
f i ( t , x i , t ) = k = 1 m i α i k 0 ( t ) t τ i k ( t ) t x i ( s ) K i k ( t s ) 1 + c i k ( t s ) x i ( s ) d s .
Proceeding as for (27), we derive that f i ( t , x i , t ) β i ( t ) h i ( x ) , with h i ( x ) = x 1 + c i ¯ x , where c i ¯ = max k sup t 0 c i k ( t ) . In fact, when κ i 0 , (35) is a particular case of (27), with nonlinearities of the form (29) and α i = 1 . Thus (see also Remark 6), when κ i ( t ) are bounded, Theorem 4 applies to (35). Namely, if lim inf t β i ( t ) > 0 and lim inf t β i ( t ) v i d i ( t ) v i j = 1 n a i j ( t ) v j > 1 for some v i > 0 and 1 i n , then both (h5), (h5*) are satisfied, and (35) is persistent. Since (35) is a cooperative system, the theory of monotone systems for DDEs [19] can be used to derive that the system is dissipative if in addition β i ( t ) , a i j ( t ) are all bounded (see [5] for details).
For systems (34), assuming that all the coefficients are bounded (and κ i ( t ) κ 0 > 0 ), its permanence was shown in [5] (Theorem 5.3), provided that there exists a vector v > 0 such that:
M v > 0 , where M = d i a g ( β ̲ 1 d ¯ 1 , , β ̲ n d ¯ n ) + a ̲ i j
and the notations f ¯ = sup t 0 f ( t ) , f ̲ = inf t 0 f ( t ) for f : R + R + bounded were used, a much stronger hypothesis. On the other hand, a criterion for the permanence of (34) was obtained by Győri et al. [9] (Corollary 3.7) under some different assumptions and not assuming the a priori boundedness of all coefficients: Instead of lim inf t d i ( t ) > 0 , it was imposed that d i ( t ) > 0 with 0 d i ( t ) d t = , and that there are uniform bounds 1 < α i β i ( t ) / d i ( t ) C i , a i j / d i ( t ) M i j , κ i ( t ) / d i ( t ) C i for all i , j , and t 0 ; nevertheless, the assumptions in [9] are much more restrictive when lim inf t d i ( t ) > 0 .
In particular, our results allow to conclude that the scalar model:
x i ( t ) = k = 1 m α k 0 ( t ) t τ k ( t ) t x ( s ) K k ( t s ) 1 + c k ( t s ) x ( s ) d s d ( t ) x ( t ) κ ( t ) x 2 ( t ) ,
is permanent if β ( t ) : = k = 1 m α k 0 ( t ) 0 τ k ( t ) K k ( s ) d s β ̲ > 0 and β ( t ) / d ( t ) α > 1 for t large. As expected, in the case of the autonomous Equation (33) with finite delay τ , this requirement translates as γ 0 τ e μ s k ( s ) d s > μ , shown in [23] to be a sharp condition for the existence and global stability of a positive equilibrium. Of course, it would be interesting to extend our results for (35) and (36) to the case of infinite delay.

3.3. A Chemostat Model

We consider a delayed nonautonomous version of the chemostat model with two microorganisms in competition studied by Ellermeyer [24]:
S ( t ) = d ( t ) ( S 0 S ( t ) ) b 1 ( t ) p 1 ( S ( t ) ) x 1 ( t ) b 2 ( t ) p 2 ( S ( t ) ) x 2 ( t ) x i ( t ) = d ( t ) x i ( t ) + e t τ i ( t ) t d ( s ) d s b i ( t ) p i ( S ( t τ i ( t ) ) ) x i ( t τ i ( t ) ) , i = 1 , 2 ,
where S 0 > 0 , the functions d , b i : R + ( 0 , ) , τ i , p i : R + R + are continuous, τ i ( t ) are bounded, and p i ( 0 ) = 0 , p i ( x ) > 0 for x > 0 , i = 1 , 2 . Here, S 0 is the initial concentration of the nutrient in the chemostat vessel free of microorganisms, S ( t ) is the concentration of the nutrient, x i ( t ) is the biomass of the species i in the growth vessel, and τ i ( t ) is the time for the conversation of nutrient into its biomass ( i = 1 , 2 ), at time t. Typically, models with Michaelis–Menton functional responses p i ( x ) = m i x 1 + k i x or other strictly increasing, smooth functions have been considered. See [24,25,26,27] and references therein, for further explanations of the model.
One easily sees that the ordered interval S : = [ 0 , S 0 ] × R + × R + C + is positively invariant for (37). After the change of variable u ( t ) = S 0 S ( t ) , (37) reads as:
u ( t ) = d ( t ) u ( t ) + b 1 ( t ) p 1 ( S 0 u ( t ) ) x 1 ( t ) + b 2 ( t ) p 2 ( S 0 u ( t ) ) x 2 ( t ) x i ( t ) = d ( t ) x i ( t ) + e t τ i ( t ) t d ( s ) d s b i ( t ) p i ( S 0 u ( t τ i ( t ) ) ) x i ( t τ i ( t ) ) , i = 1 , 2 .
For the latter system, the nonlinearites are given by:
f 1 ( t , ϕ ) = b 1 ( t ) p 1 ( S 0 ϕ 1 ( 0 ) ) ϕ 2 ( 0 ) + b 2 ( t ) p 2 ( S 0 ϕ 1 ( 0 ) ) ϕ 3 ( 0 ) f i + 1 ( t , ϕ ) = e t τ i ( t ) t d ( s ) d s b i ( t ) p i ( S 0 ϕ 1 ( τ i ( t ) ) ) ϕ i + 1 ( τ i ( t ) ) , i = 1 , 2 ,
for t 0 , ϕ = ( ϕ 1 , ϕ 2 , ϕ 3 ) S . If lim inf t d ( t ) = : d ̲ > 0 and b 1 ( t ) , b 2 ( t ) are bounded, Theorem 1 implies that (38) is dissipative (in S ). Moreover, for any solution ( S ( t ) , x 1 ( t ) , x 2 ( t ) ) in S of (37), if there is a sequence t k such that S ( t k ) 0 , S ( t k ) 0 , we get:
S ( t k ) = d ( t k ) ( S 0 S ( t k ) ) b 1 ( t k ) p 1 ( S ( t k ) ) x 1 ( t k ) b 2 ( t k ) p 2 ( S ( t k ) ) x 2 ( t k ) 0 S 0 d ̲ > 0 ,
a contradiction. Thus S ( t ) persists, i.e., lim inf t S ( t ) > 0 . Assuming in addition the uniform persistence of S ( t ) , the uniform persistence of the microorganisms is also obtained:
Theorem 7.
For (37) under the above conditions, assume that S ( t ) uniformly persists, with lim inf t S ( t ) s 0 > 0 for any solution ( S ( t ) , x 1 ( t ) , x 2 ( t ) ) with initial condition in S . If b 1 ( t ) , b 2 ( t ) are bounded and:
lim inf t π i e t τ i ( t ) t d ( s ) d s b i ( t ) d ( t ) > 0 , i = 1 , 2 ,
where π i = min x [ s 0 , S 0 ] p i ( x ) , then both species x 1 ( t ) , x 2 ( t ) uniformly persist.
Proof. 
Let δ > 0 be such π i e t τ i ( t ) t d ( s ) d s b i ( t ) d ( t ) > δ for t 1 , i = 1 , 2 . Fix ε > 0 , and define π i ( ε ) = min x [ s 0 ε , S 0 ] p i ( x ) ( i = 1 , 2 ) . For solutions X ( t ) = ( u ( t ) , x 1 ( t ) , x 2 ( t ) ) of (38) with initial condition in S , for t sufficiently large we obtain the estimates f i ( t , X t ) β i ( t ) h i ( X t ̲ ) , where:
β 1 ( t ) = π 1 ( ε ) b 1 ( t ) + π 2 ( ε ) b 2 ( t ) , β i + 1 ( t ) = π i ( ε ) e t τ i ( t ) t d ( s ) d s b i ( t ) , i = 1 , 2 ,
and h i ( x ) = x , i = 1 , 2 , 3 . With the notation in Section 2, D ( t ) = diag ( d ( t ) , d ( t ) , d ( t ) ) , A ( t ) 0 and M ( t ) = diag ( β 1 ( t ) d ( t ) , β 2 ( t ) d ( t ) , β 3 ( t ) d ( t ) ) . For ε > 0 sufficiently small and t large, we obtain β i ( t ) d ( t ) δ , i = 2 , 3 , and thus β 1 ( t ) d ( t ) δ also holds. In this way, (H5) is satisfied, as well (H5*) (conf. Remark 1). From Theorem 3, (38) is uniformly persistent. □
For the autonomous environment model in [24], it was shown that the co-existence of two different competing species is not possible: Either the washout equilibrium ( S 0 , 0 , 0 ) is globally attractive, or only one species persists. For the nonautonomous version (38), the above result leads to the persistence of both microorganisms, provied that π i b i ( t ) / d ( t ) α > 1 and the time-delays are small enough, t 1 .

3.4. A Planar System with Distributed Delays

We consider the planar system studied in [28]:
x ( t ) = d 1 ( t ) t τ 1 ( t ) t f 1 ( y ( s ) ) d s ν 1 ( t , s ) x ( t ) y ( t ) = d 2 ( t ) t τ 2 ( t ) t f 2 ( x ( s ) ) d s ν 2 ( t , s ) y ( t ) ,
where:
(a1)
d i , f i , τ i : R + R + are continuous, f i ( x ) > 0 for x > 0 , τ i ( t ) [ 0 , τ ] for some τ > 0 , ν i ( t , s ) are defined for ( t , s ) R + × [ τ , ) , continuous in t, left-continuous and nondecreasing in s, and such that t τ i ( t ) t d s ν i ( t , s ) = 1 , for i = 1 , 2 , t 0 ;
(a2)
f i ( x ) is nondecreasing on [ 0 , r ] for some r > 0 and lim inf x f i ( x ) > 0 , i = 1 , 2 ;
(a3)
For each i = 1 , 2 , either f i ( 0 ) > 0 or f i ( 0 ) = 0 and there exists f i ( 0 ) > 1 ;
(a4)
lim inf t d i ( t ) > 0 , i = 1 , 2 .
For several real applications of this model, such as neuron networks and Lotka–Volterra cooperative systems, see [28] and references therein.
Theorem 8.
If (a1), (a4) hold, (39) is uniformly persistent. If in addition d i ( t ) , f i ( x ) are bounded, then the system is permanent.
Proof. 
If l i : = lim inf x f i ( x ) f i ( 0 ) , replace f i ( x ) by a continuous function f ^ i ( x ) such that f ^ i ( x ) = α x on [ 0 , ε ) (for a small ε > 0 and α > 1 ), and f ^ i ( x ) f i ( x ) on R + . In this way, we may assume that r > 0 in (a2) is such that f ^ i ( r ) < lim inf x f ^ i ( x ) , i = 1 , 2 . With the previous notations, where L i j ( t ) 0 ( i , j = 1 , 2 ) , (H1), (H3) are satisfied and d i ( t ) t τ i ( t ) t f i ( u ( s ) ) d s ν i ( t , s ) d i ( t ) f ^ i ( u ̲ ) , for a positive function u ( t ) with u ̲ : = inf t u ( t ) r . Choosing a smaller positive r > 0 if necessary, from (a3) there is α > 1 such that (H4) holds with β i ( t ) = α d i ( t ) and h i ( x ) = x . Hence (H5), (H5*) are satisfied, and Theorem 3 gives the result. □
In fact, an interesting method was developed in [28] to study the global behavior of (39), under some other strong assumptions: Instead of (a2) and (a3), it was assumed that
(a2*)
f i ( x ) are strictly increasing on R + , with f ( x ) > 0 for x > 0 , i = 1 , 2 ;
(a3*)
The equation f 1 1 ( x ) = f 2 ( x ) has a unique positive solution K > 0 , with f 2 ( x ) f 1 1 ( x ) ( x K ) < 0 for x > f 1 ( 0 ) , x K ;
(a4*)
0 d i ( t ) = , i = 1 , 2 .
Note that (a2*), (a3*) imply that (39) is cooperative and that there exists a unique positive equilibrium, given by ( K , f 2 ( K ) ) . Under (a1), (a2*)–(a4*), the global attractivity of such an equilibrium was established in [28] [Theorem 3.1]. We stress however that our result on persistence and permanence does not a priori assume the existence of a positive equilibrium, nor that f i are increasing on R + .

3.5. On the Necessity of Some Assumptions

Some illustrative examples showing the necessity of our assumptions are now presented.
Example 2.
Consider the planar system:
x 1 ( t ) = t η x 1 ( t ) + ( t η 1 ) x 2 ( t τ 1 ( t ) ) + β h 1 ( t , x 1 ( t σ 1 ( t ) ) x 2 ( t ) = t η x 2 ( t ) + ( t η 1 ) x 1 ( t τ 2 ( t ) ) + β h 2 ( t , x 2 ( t σ 2 ( t ) ) , t 1 ,
with nonlinearities of either the Mackey–Glass or Nicholson type,
h i ( t , x ) = x 1 + c i ( t ) x ν i with ν i 1 , or h i ( t , x ) = x e c i ( t ) x ,
and η > 0 , β > 1 , where the delays τ i ( t ) , σ i ( t ) are nonnegative, continuous, and bounded, c i ( t ) are positive, continuous, and 0 < c i ( t ) c ¯ i for some constants c ¯ i , i = 1 , 2 . With the previous notation, d i ( t ) = t η , a i i ( t ) = 0 , β i ( t ) β > 1 , i = 1 , 2 , and a 12 ( t ) = a 21 ( t ) = t η 1 , thus:
D ( t ) = diag ( t η , t η ) , A ( t ) = 0 t η 1 t η 1 0 , M ( t ) = β t η t η 1 t η 1 β t η .
As,
[ D ( t ) A ( t ) ] 1 1 = 1 1 , M ( t ) 1 1 = β 1 β 1 ,
(H2), (h5) (and thus also (h5*)) hold with the common vector 1 = ( 1 , 1 ) . Since a 12 ( t ) , a 21 ( t ) are not bounded, we cannot deduce that (40) is permanent (nor that (H2*) is satisfied). Nevertheless, if there are no delays in the linear part, i.e., τ i ( t ) = τ 2 ( t ) 0 , Corollary 2 yields the permanence of (40).
Example 3.
Consider the scalar equation:
x ( t ) = d ( t ) x ( t ) + a ( t ) x ( t τ ( t ) ) + β ( t ) h ( x ( t ) ) , t 0 ,
where τ ( t ) is continuous with 0 τ ( t ) τ , h ( x ) = x 2 if 0 x 1 1 if x 1 and
a ( t ) = μ t + C τ ( t ) τ ( t ) , β ( t ) = μ 1 t + C t + C 1 , d ( t ) = a ( t ) + ( β ( t ) + 1 ) 1 t + C + μ .
for some constants C > max ( τ , 1 ) , μ > 0 , μ 1 > μ + 1 / C . One easily sees that this equation has the solution φ ( t ) = 1 t + C , thus (41) is not persistent. Note that β ( t ) is bounded and d ( t ) a ( t ) > μ , β ( t ) d ( t ) + a ( t ) = μ 1 μ 1 t + C μ 1 μ 1 C > 0 , therefore (H2), (H5), and (H5*) are satisfied. However, hypothesis (H4) is not fulfilled, because h ( 0 ) = 0 (conf. Remark 2).
Example 4.
Consider the system:
x 1 ( t ) = ( a ( t ) + d 1 ( t ) ) x 1 ( t ) + a ( t ) x 2 ( t ) + β ( t ) h ( x 1 ( t τ ) ) , x 2 ( t ) = ( a ( t ) + d 1 ( t ) ) x 2 ( t ) + a ( t ) x 1 ( t ) + β ( t ) h ( x 2 ( t τ ) ) , t 0 ,
where τ > 0 , a ( t ) is continuous, nonnegative and bounded, h ( x ) = x 1 + x and
d 1 ( t ) = μ t + C 1 τ , β ( t ) = t + C + 1 τ t + C d 1 ( t ) 1 t + C ,
for some τ ( 0 , 1 ) and C > τ . It is easy to verify that x ( t ) = 1 t + C , 1 t + C is a solution of (42), thus this system is not permanent. Note that d 1 ( t ) μ C 1 τ , β ( t ) d 1 ( t ) = μ t + C + 1 τ ( t + C ) 2 μ / 2 for t 1 , hence (H2) and (H5) are satisfied. However Theorem 3 does not apply, because (H5*) does not hold.

3.6. A Numerical Example

A concrete numerical simulation is now included. For simplicity, a scalar model is chosen.
In their celebrated paper [29], Mackey and Glass proposed appropriate models to describe the hematopoiesis (the process of production, multiplication, and specialization of blood cells in bone marrow). Here, we consider a 1-periodic hematopoiesis model with monotone production rate:
y ( t ) = 1 + 1 2 cos ( 2 π t ) y ( t ) + η 1 1 + 1 2 cos ( 2 π t ) 1 + y ( t 1 cos ( 2 π t ) ) n + η 2 1 + 1 2 sin ( 2 π t ) 1 + y ( t 2 ) n , t 0 ,
with n > 1 and η 1 , η 2 > 0 . The model possesses two discrete delays, τ 1 ( t ) = 1 + cos ( 2 π t ) and τ 2 ( t ) 2 . Define d ( t ) = 1 + 1 2 cos ( 2 π t ) , b 1 ( t ) = η 1 ( 1 + 1 2 cos ( 2 π t ) ) , b 2 ( t ) = η 2 ( 1 + 1 2 sin ( 2 π t ) ) , b ( t ) = b 1 ( t ) + b 2 ( t ) . Since all the coefficients are bounded, there is some M > 0 such that all solutions of (43) satisfy lim sup t y ( t ) < M . Write (43) as y ( t ) = d ( t ) y ( t ) + f ( t , y t ) . Let η = η 1 + 1 3 η 2 , fix any α > 1 and r ( 0 , M ) such that α r = η 1 + M n . With β ( t ) = α b ( t ) , h ( x ) = x on [ 0 , r ] , h ( x ) = η α ( 1 + M n ) on [ r , ) , one has f ( t , ϕ ) β ( t ) h ( ϕ ̲ ) whenever ϕ C 0 + , and β ( t ) / d ( t ) α for all t. Hence, (H4), (H5) are satisfied. From Theorem 3, (43) is permanent.
Take e.g., η 1 = 0.2 , η 2 = 0.3 , and n = 3 . From [2], (43) has a positive 1-periodic solution. The simulations below suggest that the positive solutions of (43) approach such a periodic solution as t (see Figure 1).

4. Discussion and Open Problems

In this paper, we have proven the persistence and permanence of delayed differential systems (8) which incorporate distributed delays in both the linear and nonlinear parts and are in general noncooperative. Moreover, not all the coefficients are required to be bounded. The main theorem, Theorem 3, extends known results in recent literature [2,3,6,7,9,16,18], as it applies to a broad family of nonautonomous delay differential systems.
Once the permanence of (8) is guaranteed, several open questions arise and should be addressed. First, it would be interesting to have explicit lower and upper uniform bounds for all positive solutions, as investigated in [3,4,5,8,9,16] for cooperative scalar or n-dimensional DDEs and in [7,18] for noncooperative systems. Secondly, the global stability of DDEs is a matter of crucial importance in applications, therefore a relevant task is to propose sufficient conditions forcing x ( t ) y ( t ) 0 as t , for any two positive solutions x ( t ) , y ( t ) of (8). In the case of nonautonomous noncooperative models, it is however clear that the response to these two questions depends on the specific nonlinearities. In a forthcoming paper, these topics will be addressed for generalized Nicholson systems. For periodic n-dimensional DDEs, it has been proven [30] that in some settings the permanence implies the existence of a positive periodic solution. In this context, a stability result will show that such a periodic solution is a global attractor of all positive solutions.
It is worthwhile mentioning that, in the last few years, the stability of nonautonomous linear DDEs has received a great deal of attention, and several methods have been used to obtain explicit sufficient conditions for the asymptotic and exponential asymptotic stability of a general linear system (4), see e.g., [21,31] and references therein. Actually, both delay independent and delay-dependent criteria for the stability of linear DDEs were given in e.g., [21,31], the latter also with possible infinite delays. Since the exponential stability of (4) is a key ingredient to show the permanence of (8), this leads us to two natural lines of future research, explained below.
The first one is to replace assumption (H2) or (H2*) – which forces (4) to possess diagonal terms without delay which dominate the effect of the delayed terms – by a condition depending on the size of delays, in such a way that (4) maintains the exponential asymptotic stability, and further analyze how such a condition interplays with the assumption (H5).
Another open problem is to study the persistence and permanence of systems of the form (8) with unbounded delays. DDEs with infinite delay are surely more challenging: Not only an admissible phase space satisfying some fundamental set of axioms should be chosen, but most techniques for finite delays do not apply for such equations. There has been some recent work on the permanence for scalar nonautonomous DDEs with infinite delay, see e.g., [23]. In the case of multidimensional DDEs with infinite delay, the work in [5] only contemplates situations of cooperative systems, namely of the form x i ( t ) = F i ( t , x t ) x i ( t ) G i ( x t ) ( 1 i n ) with F i , G i cooperative and F i ( t , x ) sublinear in x R + . For the case of possible unbounded coefficients and nonmonotone nonlinearities in (8), it is clear that the technique developed in the proof of Theorem 3 does not apply to systems with infinite delay, since it relies on a step-wise iterative argument on intervals of lenght τ , where τ is the supremum of all delays. Thus, new tools and arguments to tackle the difficulty must be proposed. This open problem is a strong motivation for future investigation.
The treatment of mixed monotonicity models, in what concerns questions of permanence, is another topic deserving attention, since they appear naturally in real-world applications. In fact, there has been an increasing interest in DDEs with mixed monotonicity, where the nonlinear terms involve one or more functions with different delays e.g., of the form f ( t , x ( t τ ( t ) ) , x ( t σ ( t ) ) ) , with f ( t , x , y ) monotone increasing in the variable x and monotone decreasing in y. As illustrated by Berezansky and Braveman [15], though small delays are in general harmless, the presence of two or more delays in the same nonlinear function may change drastically the global properties of the solutions. The permanence and stability of DDEs with nonlinearities of mixed monotonicity have been analyzed in [1,8,15,16,32]. As far as the author knows, only the case of discrete delays has been dealt with. New tools are required to handle the case of systems with mixed monotonicity in the nonlinear terms.

Funding

This research was funded by by National Funding from FCT-Fundação para a Ciência e a Tecnologia (Portugal) under project UIDB/04561/2020.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks the reviewers for their valuable comments and suggestions, which led to the significant improvements of this manuscript.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the writing of the manuscript.

Abbreviations

The following abbreviation is used in this manuscript:
DDEdelay differential equation
ODEordinary differential equation
PDEpartial differential equation
SISsusceptible-infective-susceptible

References

  1. Berezansky, L.; Braverman, E. Boundedness and persistence of delay differential equations with mixed nonlinearity. Appl. Math. Comput. 2016, 279, 154–169. [Google Scholar] [CrossRef] [Green Version]
  2. Berezansky, L.; Braverman, E.; Idels, L. Mackey-Glass model of hematopoiesis with monotone feedback revisited. Appl. Math. Comput. 2013, 219, 4892–4907. [Google Scholar] [CrossRef]
  3. Berezansky, L.; Idels, L.; Troib, L. Global dynamics of Nicholson-type delay systems with applications. Nonlinear Anal. Real Word Appl. 2011, 12, 436–445. [Google Scholar] [CrossRef]
  4. Faria, T. A note on permanence of nonautonomous cooperative scalar population models with delays. Appl. Math. Comput. 2014, 240, 82–90. [Google Scholar] [CrossRef]
  5. Faria, T. Persistence and permanence for a class of functional differential equations with infinite delay. J. Dyn. Diff. Equ. 2016, 28, 1163–1186. [Google Scholar] [CrossRef] [Green Version]
  6. Faria, T.; Obaya, R.; Sanz, A.M. Asymptotic behaviour for a class of nonmonotone delay differential systems with applications. J. Dyn. Diff. Equ. 2018, 30, 911–935. [Google Scholar] [CrossRef] [Green Version]
  7. Faria, T.; Rőst, G. Persistence, permanence and global stability for an n-dimensional Nicholson system. J. Dyn. Diff. Equ. 2014, 26, 723–744. [Google Scholar] [CrossRef] [Green Version]
  8. Győri, I.; Hartung, F.; Mohamady, N.A. On a nonlinear delay population model. Appl. Math. Comput. 2015, 270, 909–925. [Google Scholar] [CrossRef]
  9. Győri, I.; Hartung, F.; Mohamady, N.A. Boundedness of positive solutions of a system of nonlinear delay equations. Discrete Contin. Dyn. Syst. Ser. B 2018, 23, 809–836. [Google Scholar] [CrossRef] [Green Version]
  10. Obaya, R.; Sanz, A.M. Is uniform persistence a robust property in almost periodic models? A well-behaved family: Almost-periodic Nicholson systems. Nonlinearity 2018, 31, 388–413. [Google Scholar] [CrossRef] [Green Version]
  11. Thieme, H.R. Mathematics in Population Dynamics; Princeton University Press: Princeton, NJ, USA, 2003. [Google Scholar]
  12. Volterra, V. LeÇons Sur la Théorie Mathématique de la Lutte Pour la Vie; Gauthier-Villars: Paris, France, 1931. [Google Scholar]
  13. Ionescu, C.M.; Lopes, A.; Copot, D.; Machado, J.A.T.; Bates, J.H.T. The role of fractional calculus in modeling biological phenomena: A review. Commun. Nonlinear Sci. Numer. Simul. 2017, 51, 141–159. [Google Scholar] [CrossRef]
  14. Rihan, F.A.; Al-Mdallal, Q.M.; AlSakaji, H.J.; Hashhish, A. A fractional-order epidemic model with rtime-delay and nonlinera incidence. Chaos Solit. Fract. 2019, 126, 97–105. [Google Scholar] [CrossRef]
  15. Berezansky, L.; Braverman, E. A note on stability of Mackey-Glass equations with two delays. J. Math. Anal. Appl. 2017, 450, 1208–1228. [Google Scholar] [CrossRef]
  16. Győri, I.; Hartung, F.; Mohamady, N.A. Permanence in a class of delay differential equations with mixed monotonicity. Electron. J. Qual. Theory Differ. Equ. 2018, 53, 1–21. [Google Scholar] [CrossRef]
  17. Takeuchi, Y.; Cui, J.; Miyazaki, R.; Saito, Y. Permanence of delayed population model with dispersal loss. Math. Biosci. 2006, 201, 143–156. [Google Scholar] [CrossRef]
  18. Liu, B. Global stability of a class of Nicholson’s blowflies model with patch structure and multiple time-varying delays. Nonlinear Anal. Real World Appl. 2010, 11, 2557–2562. [Google Scholar] [CrossRef]
  19. Smith, H.L. Monotone Dynamical Systems. In An Introduction to the Theory of Competitive and Cooperative Systems. Mathematical Surveys and Monographs; American Mathematical Society: Providence, RI, USA, 1995. [Google Scholar]
  20. Hale, J.K.; Verduyn Lunel, S.M. Introduction to Functional Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  21. Faria, T. Stability for nonautonomous linear differential systems with infinite delay. J. Dyn. Diff. Equ. 2020. [Google Scholar] [CrossRef]
  22. Arino, J.; Wang, L.; Wolkowicz, G.S.K. An alternative formulation for a delayed logistic equation. J. Theor. Biol. 2006, 241, 109–119. [Google Scholar] [CrossRef]
  23. Lin, C.J.; Wang, L.; Wolkowicz, G.S.K. An alternative formulation for a distributed delayed logistic equation. Bull. Math. Biol. 2018, 80, 1713–1735. [Google Scholar] [CrossRef]
  24. Ellermeyer, S.F. Competition in the chemostat: Global asymptotic behavior of a model with delayed response in growth. SIAM J. Appl. Math. 1994, 54, 456–465. [Google Scholar] [CrossRef]
  25. Smith, H.L.; Waltman, P. The Theory of the Chemostat: Dynamics of Microbial Competition; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  26. Fan, G.; Smith, H.L.; Thieme, H.R. Competition in the chemostat with time-dependent differential removal rates. Vietnam J. Math. 2017, 45, 153–178. [Google Scholar] [CrossRef]
  27. Smith, H.L.; Thieme, H.R. Chemostats and epidemics: Competition for nutrients/hosts. Math. Biosci. Eng. 2013, 10, 1635–1650. [Google Scholar] [PubMed]
  28. Berezansky, L.; Braverman, E. On stability of cooperative and hereditary systems with distributed delays. Nonlinearity 2015, 18, 1745–1760. [Google Scholar] [CrossRef] [Green Version]
  29. Mackey, M.C.; Glass, L. Oscillation and chaos in physiological control system. Science 1977, 197, 287–289. [Google Scholar] [CrossRef] [PubMed]
  30. Zhao, X.-Q. Permanence implies the existence of interior periodic solutions for FDEs. Intern. J. Qual. Theory Differ. Equ. Appl. 2008, 2, 125–137. [Google Scholar]
  31. Berezansky, L.; Diblík, J.; Svoboda, Z.; Smarda, Z. Exponential stability of linear delayed differential systems. Appl. Math. Comput. 2018, 320, 474–484. [Google Scholar] [CrossRef]
  32. El-Morshedy, H.A.; Ruiz-Herrera, A. Criteria of global attraction in systems of delay differential equations with mixed monotonicity. J. Differ. Equ. 2020, 268, 5945–5968. [Google Scholar] [CrossRef]
Figure 1. Numerical simulation of three solutions of (43) where η 1 = 0.2 , η 2 = 0.3 and n = 3 , with initial condition φ ( θ ) = 0.3 ( cos ( θ ) + 1 ) , φ ( θ ) = 0.8 e θ , and φ ( θ ) = 0.3 ( sin ( θ ) + 1 ) , for θ [ 2 , 0 ] .
Figure 1. Numerical simulation of three solutions of (43) where η 1 = 0.2 , η 2 = 0.3 and n = 3 , with initial condition φ ( θ ) = 0.3 ( cos ( θ ) + 1 ) , φ ( θ ) = 0.8 e θ , and φ ( θ ) = 0.3 ( sin ( θ ) + 1 ) , for θ [ 2 , 0 ] .
Mathematics 09 00263 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Faria, T. Permanence for Nonautonomous Differential Systems with Delays in the Linear and Nonlinear Terms. Mathematics 2021, 9, 263. https://doi.org/10.3390/math9030263

AMA Style

Faria T. Permanence for Nonautonomous Differential Systems with Delays in the Linear and Nonlinear Terms. Mathematics. 2021; 9(3):263. https://doi.org/10.3390/math9030263

Chicago/Turabian Style

Faria, Teresa. 2021. "Permanence for Nonautonomous Differential Systems with Delays in the Linear and Nonlinear Terms" Mathematics 9, no. 3: 263. https://doi.org/10.3390/math9030263

APA Style

Faria, T. (2021). Permanence for Nonautonomous Differential Systems with Delays in the Linear and Nonlinear Terms. Mathematics, 9(3), 263. https://doi.org/10.3390/math9030263

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop