CONTINUITY AND MONOTONICITY OF SOLUTIONS TO A GREEDY MAXIMIZATION PROBLEM
Łukasz Kruk
Maria Curie-Skłodowska
University, Lublin, Poland
Abstract
Motivated by an application to resource sharing network modelling, we consider a problem of greedy maximization (i.e., maximization of the consecutive minima) of a vector in Rn, with the admissible set indexed by the time parameter. The structure of the constraints depends on the underlying network topology. We investigate continuity and monotonicity of the resulting maximizers with respect to time. Our results have important consequences for fluid models of the corresponding networks which are optimal, in the appropriate sense, with respect to handling real-time transmission requests.
Keywords partial order, greedy maximization,
continuity, monotonicity,
resource sharing, fluid model
Mathematics subject classification 06A06, 26A15, 26A48,
68M20, 90B10, 90C35.
1 Introduction
Starting from seminal papers of Rybko and Stolyar [14], Dai [5], Bramson [2, 3, 4] and other authors, fluid models have become a standard tool in investigating long-time behaviour of complicated
queueing systems.
Such models are useful for establishing stability and obtaining hydrodynamic or diffusion limits for multiclass queueing networks and resource sharing networks with various service protocols. Using a similar
methodology in the case of real-time Earliest Deadline First (EDF) networks with resource sharing is hindered by the
lack of suitable fluid model equations. To overcome this difficulty, Kruk [13] suggested the definition of fluid models for these systems by means of an optimality property, called local edge minimality, which is known to
characterize the EDF discipline in stochastic
resource sharing networks. It turns out that the success of this approach depends on establishing suitable local properties of a vector-valued mapping F:[0,∞)→Rn,
resulting from a greedy maximization (i.e., maximization of the consecutive minima) of a vector in Rn over the admissible set At, depending on the underlying network topology and indexed by the time parameter.
A convenient way to describe the value F(t) for a given t≥0 is to define it as the maximal element of At with respect to a suitable “min-sensitive” partial ordering.
Roughtly speaking, the function F, when well behaved, determines the so-called frontiers (i.e., the left endpoints of the supports) of the states in the corresponding locally edge minimal fluid model. The idea of using frontiers for the asymptotics of EDF systems dates back to the paper of
Doytchinov et al. [7] on a G/G/1 queue, and it has been used several times since then. However, both the application of this idea to resource sharing networks and our approach to determining
the frontiers by finding
maximal elements of partially ordered sets
appear to be new.
In this paper, for any t≥0, we
construct the value of
F(t) as a solution of a nested sequence of max-min problems in a t-dependent admissible set. While none of these max-min problems is hard to solve, their number and forms vary with t in a complex, discontinuous way, making the analysis of the resulting mapping F on [0,∞) rather involved.
We investigate key properties of F,
namely, its continuity and monotonicity.
Our main results are described in more detail in Section 2.4, to follow, after the introduction of indispensable notation.
These results are used in
our forthcoming paper to establish fundamental properties of locally edge minimal fluid models, like their existence, uniqueness and stability.
We hope that the
theory developed in this paper
will be useful not only in the asymptotic analysis of EDF-like
disciplines, but also in the case of other “greedy“
scheduling policies for resource sharing networks, for example
Longest Queue First [6] or Shortest Remaining
Processing Time [15].
1.1 Notation
For sets A, B, we write A⊊B if A is a proper subset of B.
For a finite set A, let ∣A∣ denote the cardinality of A.
Let N denote the
set of positive integers and let R denote the set of real
numbers. For a,b∈R, we write a∨b (a∧b) for the
maximum (minimum) of a
and b, a+ for a∨0. Vector inequalities are to be interpreted componentwise, i.e., for a,b∈Rn, a=(a1,..,an), b=(b1,...,bn), a≤b if and only if ai≤bi for all i=1,...,n.
For a=(a1,..,an)∈Rn, we write mina for mini=1,...,nai and
Argmina for {i∈{1,...,n}:ai=mina}.
For a=(a1,..,an)∈Rn and a set A⊊{1,...,n} with {1,...,n}∖A={i1,...,ik}, where k=n−∣A∣ and i1<i2<...<ik, we identify (ai)i∈/A with (ai1,...aik)∈Rk.
By
convention,
a sum of the form ∑i=nm (⋃i=nm) with n>m or, more generally,
a sum of numbers (resp., sets)
over the empty set of indices equals zero (resp., \mathchar31).
For a set A⊆R, let A denote the closure of A.
2 The mapping F
2.1 “Min-sensitive” partial ordering on Rn
We define the relation “⪕” inductively on Rn as follows.
Definition 2.1
For a,b∈R, we write a⪕b iff a≤b.
If n≥2, for a,b∈Rn, a=(a1,..an), b=(b1,...,bn), we write a⪕b if one
of the following four cases holds.
(i) a=b,
(ii) mina<minb,
(iii) mina=minb and Argminb⊊Argmina,
(iv) mina=minb, Argmina=Argminb⊊{1,...,n} and (ai)i∈/Argmina⪕(bi)i∈/Argminb in R∣Argmina∣.
Remark 2.2
Clearly, if a⪕b, then mina≤minb. The converse is, in general, false, unless n=1. For example, the relation (0,1)⪕(1,0) does not hold.
**
The proof of the following lemmma is elementary and it is left to the reader.
Lemma 2.3
The relation “⪕” is a partial ordering on Rn.
Remark 2.4
For a,b∈Rn, the inequality a≤b implies a⪕b. In dimensions greater than one the converse is in general, false, for example, (0,2)⪕(1,1).
**
2.2 The mapping definition
Let I,J∈N and let I={1,...,I}, J={1,...,J}. For i∈I, let hi:R→R be a
continuous, nonnegative, nondecreasing function with
limx→∞hi(x)=∞ and
[TABLE]
In other words, each hi is the cumulative distribution function
of an atomless, σ-finite measure in R, with finite,
nonpositive infimum of its support. Let Gj, j∈J, be a
family of distinct, nonempty subsets
of I (not necessarily pairwise disjoint) such that
I=⋃j∈JGj.
Definition 2.5
For t≥0, we denote by F(t)=(Fi(t))i∈I the maximal element of the set
[TABLE]
with respect to the relation “⪕”.
Somewhat informally, F(t) may be thought of as the result of “greedy” maximization of a vector a∈RI, subject to the constraints defining the set At, in the following sense. We first maximize mina over a∈At, then we maximize the “next minimum” min{ai,i∈/Argmina} over the set of maximizers of the previous problem, and we continue in this way until all the entries of the maximizer a∗=F(t) are determined.
In Section 3.1 we formalize and describe in detail this nested max-min procedure
which implies,
in particular, the existence and uniqueness of the maximizer F(t).
Remark 2.6
A seemingly more general version of Definition 2.5, in which
for a fixed (possibly positive) t0∈R and t≥t0,
the set At is replaced by
[TABLE]
and the half-line (−∞,0] in (2.1) is replaced by (−∞,t0],
may be easily reduced to the case considered in Definition 2.5 by the change of variables y=x−t0, s=t−t0 and by using h~i(y)=hi(x)=hi(y+t0) instead of hi.
**
2.3 Motivation: fluid models for resource sharing networks
The need for investigating the properties of the mapping F defined above arises in the theory of fluid models for real-time networks with resource sharing.
Below, we briefly (and somewhat informally) describe this
connection. The reader may consult [13] for more details and references.
Consider a network with a finite number of resources (nodes), labelled by j=1,...,J, and a finite set of routes, labelled by i=1,...,I.
Let I={1,...,I}, J={1,...,J}. For j∈J, let Gj⊆I be the set of routes using the
resource j. For convenience, we assume that all the resources have
a unit service rate. By a flow on route i we mean a continuous
transmission of a file through the resources used by this route. We
assume that a flow takes simultaneous possession of all the
resources on its route during the transmission. Each flow in the
network has a deadline for transmission completion. Networks of this
type may be used to model, e.g., voice and video transmission,
manufacturing systems with order due dates or emergency health care
services. In what follows, by the lead time of a flow we mean the
difference between its deadline and the current time.
As in [12, 13], the time evolution of such a system may be described by the process
X(t,s)=(Z(t,s),D(t,s),T(t,s),Y(t,s)), t≥0,s∈R, where the component processes Z, D, T, Y are defined as follows.
For t≥0 and s∈R, Z(t,s)=(Zi(t,s))i∈I, where
Zi(t,s) is the number of flows on route i with lead times at time t less than or equal to s which are still present in the system at that time.
Similarly, the vectors
D(t,s)=(Di(t,s))i∈I,
T(t,s)=(Ti(t,s))i∈I denote the number of departures (i.e., transmission completions) and the cumulative transmission time
by time t corresponding to each route i of flows with lead times at time t
less than or equal to s. Let Yi(t,s)=t−Ti(t,s), i∈I, denote the cumulative idleness by time t with regard to
transmission of flows on route i with lead times at time t less than or equal to s and let Y(t,s)=(Yi(t,s))i∈I.
The process X satisfies the following network equations, valid for t~≥t≥0, s∈R:
[TABLE]
where E(t,s)=(Ei(t,s))i∈I is the corresponding external arrival process and
Si(t′,t,s) denotes the number of
transmission completions of flows on route i∈I having lead times at time t
less than or equal to s, by the time the system has spent t′ units of time transmitting these flows.
Fluid models are deterministic, continuous analogs of resource sharing networks, in which individual flows are replaced by a divisible commodity (fluid), moving along I routes with J resources (nodes).
They usually arise from formal functional law of large numbers approximations of the corresponding stochastic flow level models.
The analogs of the network equations (2.2)-(2.4) are the following fluid model equations, valid for t~≥t≥0, s∈R:
[TABLE]
where mi is the mean transmission time of a flow on route i and α=(αi)i∈I is the vector of flow arrival rates.
A system
[TABLE]
with continuous, nonnegative components,
satisfying the equations (2.5)-(2.7), together with some natural monotonicity assumptions, is called a fluid model for the resource sharing network under consideration.
To proceed further, we will introduce a class of fluid models which is, in some sense, optimal with respect to handling real-time transmission requests.
To this end,
we define a partial ordering “≪” on the space of real functions on R, which is extremely sensitive to the behaviour of the functions under comparison for small arguments.
Definition 2.7** ([13], Definition 5)**
Let f,g:R→R be such that for some a∈R we have
f≡g on (−∞,a] and let
c=sup{a∈R:f(x)=g(x)∀x≤a}.
We write f≪g if either c=∞ (i.e., f≡g on R), or
c<∞ and there exists ϵ>0 such that f≤g on
[c,c+ϵ].
Definition 2.8** (see [13], Definition 11)**
A fluid model X of the form (2.8) for a resource sharing network with ∑i∈IZi(0,⋅)≡0 on (−∞,c] for some c∈R is called locally edge minimal at a time
t0≥0 if there exists h>0 such that for any fluid model
X′ with the same
α, mi, Gi and the same initial state
Z′(0,⋅)=Z(0,⋅),
satisfying
X′(t0,⋅)=X(t0,⋅), we have ∑i∈IYi(t,⋅)≪∑i∈IYi′(t,⋅)
(equivalently,
∑i∈IZi(t,⋅)≪∑i∈IZi′(t,⋅))
for every t∈(t0,t0+h).
The fluid model X is called locally edge
minimal, if it is locally edge minimal at every t0≥0.
The intuition behind these notions is that a locally edge minimal fluid model tries to transmit as much “customer mass” corresponding to the earliest deadlines as possible, and hence its idleness with respect to such mass is as small as possible. Accordingly, such a model may be thought of as a “macroscopic” counterpart of a resource sharing network working under the Earliest-Deadline-First (EDF) protocol. Indeed, the EDF service discipline in such a network may be characterized by
an analogous notion of local edge minimality, see [13], Definition 8 and Theorems 5-7.
In a forthcoming paper, we show that existence and local uniqueness of a fluid model for given data α, mi, Gi, and initial state
Z(0,⋅),
which is locally edge minimal at a time t0≥0,
is closely related to local monotonicity of the mapping F introduced in Section 2.2, with
[TABLE]
and suitably defined, not necessarily nonpositive, xi∗, depending on t0 (see Remark 2.6).
It turns out that for large t0, the points xi∗ are also large, so under a natural assumption that the supports of
Zi(0,⋅) are bounded above, the functions hi in (2.9) are linear on [xi∗,∞). Consequently, the
linear case, investigated in Section 5.1 of this paper, is of considerable importance, because it determines the long-time
behaviour of the corresponding locally edge minimal fluid model, for
example, its stability or the form of its invariant manifold. In
particular, one of remarkable implications of the formulae developed in Section 5.1 is stability of locally edge
minimal fluid models, regardless of the underlying resource sharing
network topology. Results along this line may be found in
our forthcoming paper
2.4 Overwiev of the main results
In this paper, we investigate key properties of the mapping F:[0,∞)→RI, in particular those
which are relevant to the theory of locally edge minimal fluid models.
In Section 3, we present a detailed construction of F(t) and we provide some illustrating examples.
The main result of Section 4 is Theorem 4.2, stating that if each function hi is strictly increasing in [xi∗,∞), then the corresponding map F is continuous. This basic regularity result, together with the method of partions, introduced in Section 3.3, is useful in proving various refinements, like Lipschitz continuity of F for Lipschitz hi (Section 5.2) or an upgrade of the local monotonicity result from linear to C1 functions hi (Section 5.3). The main contribution of Section 5 is the explicit evaluation of the mapping F near zero in the linear case, implying, in particular, its local monotonicity in a neighbourhood of zero for piecewise linear hi.
As we have already mentioned, the latter fact is then generalized to (piecewise) C1 functions.
Finally, in Section 5.4, we show that, somewhat surprisingly, the mapping F may fail to be globally monotone on [0,∞), even if the corresponding functions hi are linear in [xi∗,∞).
3 The mapping construction algorithm
3.1 Construction
Fix t≥0. We define the vector F(t)=(Fi(t))i∈I as follows.
Let f(1)=f(1)(t) be the supremum of x≤t satisfying the
constraints
[TABLE]
If f(1)=t, we take Fi(t)=t for each i∈I, I(1)=I, J(1)=J, N(1)=\mathchar31 and kmax=kmax(t)=1. In what follows, we
assume that f(1)<t. By continuity of hi, x=f(1)
satisfies (3.1) and the set
{\bf J}^{(1)}=\big{\{}j\in{\bf J}:\sum_{i\in G_{j}}h_{i}(f^{(1)})=t\big{\}}
of active constraints is nonempty. (Indeed, if J(1)=\mathchar31, then
f=f(1)+ϵ also satisfies (3.1) for ϵ>0
small enough, which contradicts the definition of f(1).) Let I(1)=⋃j∈J(1)Gj and
[TABLE]
For i∈I(1), put
Fi(t)=sup{x≤t:hi(x)=hi(f(1))}. If I(1)=I, this completely determines the vector F(t). In this case, let
kmax=1. Otherwise, let K(1)=J∖(J(1)∪N(1)) and let f(2)=f(2)(t) be the
supremum of x≤t satisfying the constraints
[TABLE]
Note that Gj∖I(1)=\mathchar31 for j∈K(1), so the second sum in (3.3) is taken over a
nonempty set of indices. We also have f(1)<f(2) by
definition. If f(2)=t, we take I(2)=I∖I(1), J(2)=K(1), N(2)=\mathchar31, Fi(t)=t for each i∈I(2) and
kmax=2. If f(2)<t,
then x=f(2) satisfies (3.3) and the set J(2) of active constraints in (3.3) (i.e., those j∈K(1), for which equality holds in (3.3) with
x=f(2)) is nonempty. In this case, let I(2)=⋃j∈J(2)Gj∖I(1), N(2)={j∈J:Gj⊆I(1)∪I(2)}∖(J(1)∪J(2)∪N(1)) and put
Fi(t)=sup{x≤t:hi(x)=hi(f(2))}
for i∈I(2). If I(1)∪I(2)=I, the definition of the vector F(t) is complete and we take
kmax=2, otherwise we let K(2)=J∖(J(1)∪J(2)∪N(1)∪N(2)),
and we continue our construction as follows.
Suppose that for some k≥2, we have defined numbers f(1)<f(2)<...<f(k)<t,
nonempty, disjoint subsets J(1),…,J(k) of J, disjoint (not necessarily nonempty) subsets N(1),…,N(k) of J with
[TABLE]
and nonempty, disjoint subsets I(1), …,I(k) of I with ⋃l=1kI(l)=I such that for l=1,..,k,
[TABLE]
where
[TABLE]
and K(l)=J∖⋃p=1l(J(p)∪N(p)). Note that K(l)=\mathchar31 by (3.4) and Gj∖⋃p=1lI(p)=\mathchar31 for j∈K(l), l=1,...,k,
by (3.5), so that the second sum in (3.6) is taken
over a nonempty set of indices. (Such numbers and sets were defined in the last paragraph
for k=2.) Let
f(k+1)=f(k+1)(t) be the supremum of x≤t satisfying the constraints
[TABLE]
The inequality (3.6) implies that f(k+1)>f(k).
If f(k+1)=t, we take
I(k+1)=I∖⋃l=1kI(l), J(k+1)=K(k), N(k+1)=\mathchar31 and
Fi(t)=t for each
i∈I(k+1),
so the definition of the vector
F(t) is complete. In this case, we put kmax=k+1. If f(k+1)<t,
then x=f(k+1) satisfies (3.8)
and the set J(k+1) of active constraints in (3.8) (i.e., these j∈K(k), for which equality holds in (3.8) with x=f(k+1)) is nonempty. In this case, let I(k+1)=⋃j∈J(k+1)Gj∖⋃l=1kI(l),
define N(k+1) by (3.5) with l=k+1
and put Fi(t)=sup{x≤t:hi(x)=hi(f(k+1))} for i∈I(k+1).
This ends the k+1-th step of our construction.
If ⋃l=1k+1I(l)=I, the definition of the vector
F(t) is complete. In this case, put kmax=k+1. Otherwise, we make another (i.e., the k+2-th) step of our algorithm, taking k+1 instead of k and proceeding as above.
When the construction terminates after kmax steps, we have defined the vector F(t).
Remark 3.1
The index kmax
and the sets I(k), J(k), N(k), K(k) defined above depend on the time t. In what follows, when we want to stress this
dependence, we write kmax(t), I(k)(t), J(k)(t), N(k)(t), K(k)(t),
respectively.**
Remark 3.2
Let hˉ(x)=maxj∈J∑i∈Gjhi(x), x∈R. If hˉ is strictly increasing in [mini∈Ixi∗,∞), then, by definition,
[TABLE]
Remark 3.3
If i∈I(k)(t) for some k∈{1,...,kmax(t)} and if the function hi is strictly increasing in [xi∗,∞), then (compare (3.7)),
[TABLE]
Remark 3.4
In general, some of the sets N(k), k=1,...,kmax, may be nonempty, see Example 3.5, to follow. If j∈N(k) for some k, then
[TABLE]
This strict inequality may be interpreted as an indication of “unavoidable bottleneck idleness” in the corresponding locally edge minimal fluid network - transferring higher priority fluids by other resources does not allow j to use its full capacity on the time interval [0,t]. This phenomenon is well known in the theory of resource sharing networks and it was discussed in detail, e.g., by Gurvich and Van Mieghem [8, 9].
A mild sufficient condition for all the sets N(k) to be empty is that
for each j∈J, Gj∖⋃j′=jGj′=\mathchar31.
This corresponds to the so-called local traffic condition for the underlying network topology, under which every resource has at least one route using only that resource, see [10, 11]. The latter requirement is satisfied, for example, by linear networks, for which I=J+1 and Gj={1,j+1}, j=1,...,J.
There are, however, some important systems that do not satisfy the local traffic assumption, for example ring networks, used as counterexamples for stability of the LQF protocol [1, 6].
**
3.2 Examples
In this subsection, we provide two examples illustrating the construction of the mapping F:[0,∞)→RI defined in Subsection 3.1. The first one has relatively
simple structure, yielding time-independent kmax, J(k), I(k), N(k) in (0,∞) and linear function F. The second one, in which kmax, J(k), I(k) vary in time and F is nonlinear, indicates some of the difficulties encountered in
more general
situations.
Example 3.5
Let hi(x)=x+ for x∈R and i∈I. Then xi∗=0 for each i∈I. For t≥0 and j∈J, (3.1) takes the form
∣Gj∣x+≤t, so its maximal solution is
f(1)(t)=t/max{∣Gj∣:j∈J}. Hence, for t=0 we have f(1)(0)=0,
J(1)(0)=J, I(1)(0)=I, N(1)(0)=\mathchar31, kmax(0)=1 and Fi(0)=0 for all i. In what follows, we assume that t>0. Then
[TABLE]
and Fi(t)=max{∣Gj∣:j∈J}t for i∈I(1)(t)=⋃j∈J(1)(t)Gj.
In the remainder of this example, the time argument in J(k)(t), I(k)(t), N(k)(t), K(k)(t), will be skipped.
If ∣Gj∣=∣Gj′∣ for every j,j′∈J,
then J(1)=J, N(1)=\mathchar31 and
the definition of F is complete. Otherwise, let N(1) be given by (3.2). For j belonging to the set
[TABLE]
the constraint
(3.3) takes the form
max{∣Gj′∣:j′∈J}∣Gj∩I(1)∣t+∣Gj∖I(1)∣x+≤t,
yielding the maximal solution
[TABLE]
Moreover, J(2) is the set of j∈K(1) attaining the minimum in (3.11) and Fi(t)=f(2)(t) for i∈I(2)=⋃j∈J(2)Gj∖I(1). If I(1)∪I(2)=I (in particular, if the minimum in (3.11) is attained at every j∈K(1)), then the construction of F is complete. Otherwise we proceed similarly, until we get f(3)(t),…,f(kmax)(t), and hence all Fi(t), i∈I,
in the form of linear functions of t, with slopes depending (in an increasingly complicated way) on the sets Gj, describing the topology of the corresponding network.
Observe that in this example some sets N(k) may, in general, be nonempty. The simplest such case occurs for I=J=2, G1={1} and G2={1,2}=I. Then J(1)={2}, I(1)=I, and hence N(1)={1} and kmax=1.
**
In general, the main difficulty in analyzing the properties of the mapping F, already indicated in Remark 3.1, is the time-dependence of the index kmax and the sets I(k), J(k), k=1,...,kmax. The following example, corresponding to a simple linear network topology, illustrates this point. Note that, by Remark 3.4, in this case the sets N(k) are empty.
Example 3.6
Let I=3, J=2, G1={1,2} and G2={1,3}. Furthermore, for x∈R, let
h1(x)=(x+2)+, h2(x)=(x+1)+ and h3(x)=5x+, so that x1∗=−2, x2∗=−1, x3∗=0.
Let 0≤t≤1. Then the maximal solution of (3.1) is
f(1)(t)=t−2
and we have J(1)(t)=J. Thus, I(1)(t)=I, kmax(t)=1, and (3.10) implies that
[TABLE]
For 1<t<7/2, the maximal solution of (3.1) is
f(1)(t)=t/2−3/2, J(1)(t)={1}, so I(1)(t)=G1={1,2}.
Furthermore, K(1)(t)={2} and the maximal solution of (3.3) is f(2)(t)=(t−1)/10, so J(2)(t)={2}, I(2)(t)={3}, kmax(t)=2 and (3.10) yields
[TABLE]
For t=7/2 we still have
f(1)(t)=t/2−3/2, but J(1)(t)=J, I(1)(t)=I, kmax(t)=1. However, (3.13) still holds. Note that F1(7/2)=F2(7/2)=F3(7/2)=1/4.
Finally, for t>7/2, the maximal solution of (3.1) is
f(1)(t)=t/6−1/3, J(1)(t)={2} and I(1)(t)=G2={1,3}. Thus, K(1)(t)={1} and the maximal solution of (3.3) is f(2)(t)=5t/6−8/3, so J(2)(t)={1}, I(2)(t)={2}, kmax(t)=2 and (3.10) yields
[TABLE]
It is easy to see that the mapping F given by (3.12)-(3.14) is continuous and nondecreasing on [0,∞).
This example also indicates that there is, in general, no hope for obtaining any global
results about the auxiliary functions f(p), p>1. Here, f(2) is linear and strictly increasing in (1,7/2) and (7/2,∞), but it does not even exist in [0,1]∪{7/2}.
**
3.3 Partitions and inverses
Let J be the set of “ordered partitions” of the set J, i.e., finite sequences of subsets of J in the form (J1,...,Jk), where Ji=\mathchar31, i=1,...,k, Ji∩Jl=\mathchar31, 1≤i<l≤k, and ⋃i=1k(Ji∪Ni)=J,
(⋃i=1kJi)∩(⋃i=1kNi)=\mathchar31,
where the sets N1,...,Nk
are defined recursively as
[TABLE]
(compare the first relation in (3.4) and (3.5)).
Fix T>0 and let
[TABLE]
We label the nonempty sets of the form TD, D∈J, as T1,...Td. Clearly, ⋃i=1dTi=[0,T].
For a function g:R→[0,∞) with limx→∞g(x)=∞, let g−1 denote its (generalized) right-continuous inverse (see, e.g., [16], Section 13.6),
[TABLE]
If the function g is nondecreasing and g(x0)=0 for some x0∈R,
then
[TABLE]
Fix D=(J1,...,Jk)∈J. For t∈TD, we can write
f(p)(t), p=1,...,kmax(t), (and hence F(t)) in closed form,
using the inverse functions introduced above. To see this, let t∈TD and note that k=kmax(t). Choose j1,...,jk∈J so that jp∈Jp=J(p)(t) for each j=1,...,k. Then
[TABLE]
and for p=2,...,k, we have the recursive formulae
[TABLE]
where
[TABLE]
(Actually, the minimization with t in (3.16)-(3.17) is necessary only in the case of f(k)(t).)
4 Continuity
In general,
the function F may have jumps, as the following one-dimensional
example indicates.
Example 4.1
Let I=J=1 and let h1(x)=(x+∧1)+(x−2)+, x∈R. Then F1(t)=f(1)(t)=t/2 for 0≤t<2 and F1(t)=f(1)(t)=t/2+1 for t≥2.
**
Clearly, the discontinuity of F1=f(1) at t=2 in this example is caused by a “flat spot”, i.e., the interval [1,2] on which h1 takes the constant value 2, resulting in the jump of h1−1 (see (3.16)). In our queueing application, this corresponds to the lack of “customer mass” with deadlines in the interval [1,2] in the system, causing the frontier to
jump over this empty interval. The following theorem assures that in the absence of such “flat spots”, the function F is actually continuous.
Theorem 4.2
Assume that for every i∈I, the function hi is strictly increasing in [xi∗,∞). Then the mapping F is continuous in [0,∞).
In the proof of this result, we will use the following elementary lemma.
For the sake of completeness, we provide its justification.
Lemma 4.3
Let (X,d) be a metric space. Let y∈X and let {xn} be a sequence of elements of X such that every subsequence {xnk} of {xn} contains a further subsequence {xnkl} converging to y. Then limn→∞xn=y.
Proof. Suppose that the sequence {xn} does not converge to y. This means that there exist ϵ>0 and a subsequence {xnk} of {xn} such that d(xnk,y)≥ϵ for all k. However, we have assumed the existence of a subsequence {xnkl} of {xnk} such that liml→∞d(xnkl,y)=0, so we have a contradiction.
\hfill□
The proof of Theorem 4.2 is long and somewhat involved, so we only sketch it here,
moving most of the technical details to the appendix.
It is convenient to introduce additional
notation
[TABLE]
Clearly, the sets D(p)=D(p)(t) depend on the time t, see Remark 3.1.
Sketch of the proof of Theorem 4.2.
Let tn>0 be such that tn→t0 as n→∞.
We may assume that tn≤t0+1 for all n.
Our aim is to show that F(tn)→F(t0), i.e., for every i∈I, we have
[TABLE]
Without
loss of generality (passing to a subsequence if necessary) we may
assume that for every m,n≥1, we have
kmax(tm)=kmax(tn) and J(p)(tm)=J(p)(tn) (hence I(p)(tm)=I(p)(tn), N(p)(tm)=N(p)(tn)) for
p=1,...,kmax(tm). Consequently, in what follows, we will simply write kmax, instead of kmax(tn), n≥1, J(p) instead of J(p)(tn), n≥1,
e.t.c..
By definition, for n≥1 and p=1,...,kmax,
mini∈Ixi∗≤f(p)(tn)≤tn≤t0+1. Hence, by Lemma 4.3, without loss of generality (passing to a subsequence if necessary) we may assume that the sequences {f(p)(tn)},
p=1,...,kmax, converge. Let
[TABLE]
Remark 3.2 implies that f(1) is continuous, hence
[TABLE]
First suppose that t0=0. We have
[TABLE]
J(1)(0)=J, I(1)(0)=I, kmax(0)=1 and Fi(0)=xi∗ for all i∈I. Let i∈I. Then i∈I(k) for some k∈{1,...,kmax}. By the definition of f(k), 0≤hi(f(k)(tn))≤tn, and hence, as tn↓0, by (2.1),
[TABLE]
so Fi is continuous at [math].
Hence, in what follows, we will assume that t0>0.
We will first consider the case in which f(1)(tn)=tn for
all n≥1. Then f(1)(t0)=t0 and thus, for each i∈I,
Fi(tn)=tn→t0=Fi(t0). Similarly, if f(1)(t0)=t0 then f(1)(tn)→t0, so for each i∈I, the inclusion Fi(tn)∈[f(1)(tn),tn] implies that Fi(tn)→t0=Fi(t0). Therefore, we may additionally assume
[TABLE]
We prove (4.2) inductively for i∈I(l)(t0), l=1,...,kmax(t0). In the l-th inductive step, we define
[TABLE]
It is easy to check that
[TABLE]
where Bp(l)=I(p)∩I(l)(t0).
We show that there exists pˉl∈{bl,...,pl} such that
[TABLE]
By (3.10) and (4.10), for i∈Bp(l), p=bl,...,pˉl, as n→∞, we have
[TABLE]
We also argue that if pˉl<pl, then for i∈Bp(l), p=pˉl+1,...,pl, we have
f(l)(t0)<f∞(p)≤xi∗,
so by (3.10), as n→∞,
[TABLE]
Finally, the equations (4.9), (4.11)-(4.12) imply
(4.2) for every i∈I(l)(t0).
The details of the above inductive argument may be found in the Appendix.
5 Monotonicity
In this section, we investigate monotonicity of the mapping F.
It turns out that, in general, F fails to be globally
nondecreasing, even if the functions hi are piecewise linear,
see Section 5.4. However, under suitable assumptions on
hi, monotonicity of F
in some neighbourhood of
[math]
may be established. More precisely, our goal is to find T>0 such that
for every 0≤t<t~<T and i∈I, we have
[TABLE]
This is done in Section 5.1 for piecewise linear functions hi, i∈I.
5.1 Local monotonicity in the linear case
In this subsection we assume that for all i∈I,
[TABLE]
where ρi, i∈I, are given positive constants.
Without loss of generality we also assume that
[TABLE]
Let m∗∈{1,...I} and 0=n0<1≤n1<n2<...<nm∗=I be such that
[TABLE]
Let yk∗=xnk∗, k=1,...,m∗,
and let ym∗+1∗=∞. By (4.5) and
(5.4), we have f(1)(0)=x1∗=y1∗. It follows from
Remark 3.2 that f(1)(⋅) is continuous and strictly
increasing in [0,∞).
Let t1∗=(f(1))−1(y2∗) if m∗≥2 and t1∗=∞ otherwise.
Let
[TABLE]
We have
[TABLE]
(see Remark 3.2). Using (5.7),
we get t1∗=a(1)(y2∗−y1∗).
Note that J(1) (and hence I(1), N(1)) is constant in (0,t1∗).
If x1∗=0 and a(1)≤1, then m∗=1 and t1∗=∞. In this case, (5.7) implies that f(1)(t)=t for t≥0, so Fi(t)=t for all t≥0, i∈I. In the remainder of this section we assume that either x1∗<0, or a(1)>1. In particular, we have
[TABLE]
where tˉ1=∞ if a(1)≥1 and tˉ1=a(1)x1∗/(a(1)−1) otherwise.
In what follows, we consider only
t∈[0,t1∗∧tˉ1).
If I(1)=I, then kmax=1, and hence the numbers f(k)(t), k=1,..,kmax, have already been defined.
Assume that I(1)=I. Then
we continue our construction by induction as follows.
Assume that for some k≥1, there exist strictly positive
numbers tl∗, tˉl, l=1,...,k, such that for all 0<t<tk:=min1≤l≤k(tl∗∧tˉl) the
sets I(l)=I(l)(t), J(l)=J(l)(t), N(l)=N(l)(t), l=1,...,k, do not
depend on t and ⋃l=1kI(l)=I,
i.e., kmax(t)>k.
For l=1,...,k, let i(l)=min(I∖⋃p=1l−1I(p)) and let
m(l)∈{1,...,m∗} be such that xi(l)∗=ym(l)∗.
By definition, 1=i(1)<i(2)<...<i(k)≤I and
1=m(1)≤m(2)≤...≤m(k)≤m∗. Furthermore, we assume that for each l=1,...,k, there exists a constant a(l)>0 such that
[TABLE]
(For notational convenience, for l>1, we have defined f(l)(0) in (5.9) by continuity, i.e., as f(l)(0+)=xi(l)∗, although J(1)(0)=J and hence f(l)(0) has not been defined by the algorithm from Section 3.1.)
Note that all the above assumptions have already been verified for k=1.
By (3.10), (5.2) and (5.9),
for t∈(0,tk), j∈K(k),
the equation (3.6) implies
[TABLE]
while the equation (3.8) takes the form
[TABLE]
which, in turn, is equivalent to
[TABLE]
where
[TABLE]
By (5.10),
bj(k+1)>0
for each j∈K(k).
Let i(k+1)=min(I∖⋃p=1kI(p)), let
m(k+1)∈{1,...,m∗} be such that xi(k+1)∗=ym(k+1)∗
and let
[TABLE]
Recall that f(k+1)(t) is the supremum of x≤t satisfying the constraints (5.12).
By an argument similar to the one used for f(1) above, one may check that f(k+1)(⋅) is continuous and strictly increasing in (0,tk).
Let tk+1∗=∞ if either m(k+1)=m∗,
or f(k+1)(tk−)≤ym(k+1)+1∗, and
tk+1∗=(f(k+1))−1(ym(k+1)+1∗)
otherwise.
By definition, f(k+1)(0):=f(k+1)(0+)=xi(k+1)∗=ym(k+1)∗
(see the notational remark following (5.9)),
and, more generally,
for t≤tk∧tk+1∗,
we have
[TABLE]
Using (5.13),
we get tk+1∗=a(k+1)(ym(k+1)+1∗−ym(k+1)∗)
unless tk+1∗=∞. Note that J(k+1) (and hence
I(k+1), N(k+1)) is constant in
(0,tk∧tk+1∗).
If ym(k+1)∗=0 and a(k+1)<1, then m(k+1)=m∗ and tk+1∗=∞.
In this case,
(5.13) implies that f(k+1)(t)=t for t∈[0,tk), so kmax(t)=k+1, Fi(t)=t for such t and i∈I∖⋃p=1kI(p).
In what follows, we assume that either ym(k+1)∗<0, or a(k+1)≥1. In particular,
[TABLE]
where tˉk+1=∞ if a(k+1)≥1 and tˉk+1=a(k+1)xi(k+1)∗/(a(k+1)−1) otherwise.
This ends the k+1-th step of our construction.
If ⋃l=1k+1I(l)=I, then for 0<t<tk+1=min1≤l≤k+1(tl∗∧tˉl), we have
kmax(t)=k+1 and
the definition of
f(p)(t), p=1,...,kmax(t), is complete. Otherwise, we make another (i.e., the k+2-th) step of our algorithm, taking k+1 instead of k and proceeding as above.
When the construction terminates, in the time interval (0,tkmax), we have the sets
I(l)=I(l)(t), l=1,...,kmax, constant in t,
and f(l)(t), l=1,...,kmax, in the form of strictly increasing linear functions. Let i∈I and let k∈{1,...,kmax} be such that i∈I(k). Let 0≤t<t~<tkmax. Then, by (3.10),
Fi(t)=f(k)(t)∨xi∗≤f(k)(t~)∨xi∗=Fi(t~),
and (5.1) follows.
Remark 5.1
Since the above argument is local in time, it actually requires only that for each i∈I,
(5.2) holds in some neighborhood of xi∗, i.e., for all x<Xi∗, where Xi∗>xi∗ are given constants.
(Without loss of generality we may further assume that Xi′∗=Xi′′∗ if i′,i′′∈I and xi′∗=xi′′∗.)
We only have to
restrict t in each step of our construction to the interval [0,tk′) with
tk′=min1≤l≤k(tl∗∧tˉl∧tˉˉl), where
tˉˉl=∞ if f(l)(tl−)≤Xi(l)∗ and tˉˉl=(f(l))−1(Xi(l)∗) otherwise.
In this case, (5.1) holds for
all i∈I and
0≤t<t~<tkmax′. **
5.2 Reduction lemmas and Lipschitz continuity
In this subsection we consider the case of general hi, assuming only that
each hi is strictly increasing in [xi∗,∞). Fix T>0 and
recall the sets T1,...Td from Subsection 3.3. The following lemma reduces the problem of establishing
monotonicity of F to showing its monotonicity under the additional assumption that kmax and the sets J(p), I(p), p=1,...,kmax, are constant in t.
Lemma 5.2
Assume that for k=1,...,d, the mapping F is nondecreasing on Tk. Then F is nondecreasing on [0,T].
Proof.
First note that F is nondecreasing on Tk for each k=1,...,d. Indeed, let t, t~ be such that t<t~ and t,t~∈Tk for some k. Then for each n∈N, there exist tn,t~n∈Tk such that tn<t~n and tn→t, t~n→t~ as n→∞. By assumption, F(tn)≤F(t~n) for each n.
Letting n→∞,
by Theorem 4.2 we get F(t)≤F(t~), so F is indeed monotone
on Tk.
Let 0≤t<t~≤T and let k0 be such that t∈Tk0. Let t1=sup{s≤t~:s∈Tk0}. Then t1∈Tk0, so F(t)≤F(t1). If t1=t~, we have (5.1) and the proof is complete. Assume that t1<t~. In this case, (t1,t~]∩Tk0=\mathchar31 by the definition of t1. However, there exist k1=k0 and a sequence sn∈Tk1 such that sn↓t1. Let t2=sup{s≤t~:s∈Tk1}. Then
t1,t2∈Tk1, so F(t1)≤F(t2), and hence
F(t)≤F(t2). If t2=t~, the proof is complete, otherwise
(t2,t~]∩(Tk0∪Tk1)=\mathchar31. Let k2∈/{k0,k1} and a sequence sn∈Tk2 be such that sn↓t2. Put t3=sup{s≤t~:s∈Tk2}. As above, we have F(t2)≤F(t3), and hence
F(t)≤F(t3). After a finite number l≤d of such steps we get F(t)≤F(tl) and tl=t~, so (5.1)
holds.
\hfill□
A slight modification of the above argument yields
Lemma 5.3
Assume that for k=1,...,d, the mapping F is Lipschitz continuous on Tk. Then F is Lipschitz on [0,T].
This lemma, together with (3.16)-(3.18), may be used to provide a simple proof of the following result.
Theorem 5.4
Assume that there exist 0<c<C<∞ such that
[TABLE]
Then F is Lipschitz continuous in [0,∞).
Proof.
Fix T>0. We will argue that F is Lipschitz in [0,T] (with the Lipschitz constant independent on T). According to Lemma 5.3, it suffices to show
this in each Tl, l=1,...,d.
Fix D=(J1,...,Jk)∈J.
By (5.15), for 0≤x≤y, we have
(y−x)/(CI)≤g(y)−g(x)≤(y−x)/c,
where g is any of the inverses appearing in (3.16)-(3.17). Hence, the function f(1) defined by (3.16) is Lipschitz in TD. Using this fact, together with (3.17), (5.15), and proceeding by induction, we show Lipschitz continuity of f(p), p=2,...,k, in TD. From this, we get Lipschitz continuity of F in TD by (3.10).
\hfill□
Example 4.1 indicates that the lower bound in (5.15) is necessary for Theorem 5.4 to hold.
5.3 Local monotonicity in the C1 case
In this subsection, we assume that for each i∈I, hi∈C1([xi∗,∞)) and
[TABLE]
where hi′(xi∗) denotes the right derivative of hi at xi∗.
Again, with no loss of generality we may assume
(5.3). Define m∗, n0,...,nm∗,
y1∗,...,ym∗∗ as in Section 5.1. Because our
concern here is local monotonicity of F, without loss of
generality we may assume that hi′(x)>0 for all i∈I and
x≥xi∗ (compare Remark 5.1).
The main idea of our analysis for this case
is to consider it as a small perturbation of the linear problem considered in Section 5.1, with ρi given by (5.16). In particular, let tkmax be as in in Section 5.1 and
let kmaxL be the constant kmax(t), 0<t<tkmax,
defined there. Furthermore for k=1,...,kmaxL, let a(k)
be as in
Section 5.1 and let JL(k), IL(k), NL(k), KL(k), fL(k), denote the sets J(k)(t), I(k)(t), N(k)(t), K(k)(t), t∈(0,tkmax), and the function f(k) for the
above-mentioned linear problem, respectively.
Fix 0<T≤tkmax/2 and the partition
D=(J1,...,Jk)∈J such that there exists
a sequence tn↓0 such that tn∈TD for all n. Since
f(1)(t)<...<f(k−1)(t)<t for t∈TD,
using (3.16)-(3.17) and proceeding by
induction, it is easy to see that f(1),…, f(k−1) are the truncations of C1 functions
to TD, and f(k) is the truncation to
TD of a function in the form g(t)∧t,
where g is C1. Therefore, in order to prove that the functions f(p),
p=1,...,k=kmax, (and hence, by (3.10), Fi, i∈I) are nondecreasing in an intersection of
TD and a neighborhood of zero, it suffices to
check that
[TABLE]
where f(p)(0):=f(p)(0+)=limn→∞f(p)(tn). Consequently, by Lemma 5.2, if we verify
(5.17), then the proof of monotonicity of the mapping
F in a neighborhood of zero is complete. In what follows, we
consider only t∈TD.
We will first consider the case of x1∗=y1∗=0 (hence m∗=1, n1=I) and a(1)≤1. Then for every j∈J and t>0 small enough,
[TABLE]
with equality for j∈JL(1).
If a(1)<1, then (5.18) implies that f(1)(t)=t, and hence Fi(t)=t, i∈I, for t small enough.
If a(1)=1, then, by (5.18), f(1)(t)=t+o(t). Since f(1)(t)<f(p)(t)≤t for p=2,...,k, this implies (f(1))′(0)=(f(2))′(0)=...=(f(k))′(0)=1. Therefore, in the remainder of the proof we may assume that either x1∗<0, or a(1)>1.
By (4.5), we have f(1)(0)=x1∗=y1∗.
We claim that
[TABLE]
(compare (5.8)).
Suppose that (5.19) is false and let j∈J1∖JL(1),
j′∈JL(1).
Then, by (5.5)-(5.6), for small
t>0 we have
[TABLE]
We have obtained a contradiction, proving (5.19). Using (5.6) and (5.19), for j∈J1, for small t>0, we get
[TABLE]
yielding (5.20). If k=1, then (5.20) implies
(5.17) and our proof
is complete, otherwise we proceed by induction as follows.
Assume that for some 1≤l<k there are indices r∈{1,..,l∧kmaxL},
0=p0<1≤p1<p2<...<pr=l such that
[TABLE]
where the sets Ip are as in (3.18).
We also assume that for p=ps−1+1,...ps, s=1,...,r, we have
[TABLE]
By (5.19)-(5.20), the above assumptions hold for l=1 (hence r=1, p1=1). Note that in this case (5.22)
is vacuously true.
For p=0,...,k−1, let Kp={j′∈J:Gj∖⋃q=1pIq=\mathchar31}. Clearly, K0=J. Moreover, (3.18) implies that Kp=K(p)(t) for p=1,...k−1. We claim that
[TABLE]
Indeed,
the inclusion in (5.24) is obvious, since pr−1<l. By (5.22),
⋃p=1pr−1Ip=⋃s=1r−1IL(s), hence ⋃p=1pr−1D(p)(t)=⋃s=1r−1DL(s), yielding the equality in (5.24).
By (3.18) and (5.21) with s=r, ⋃p=pr−1+1lIp⊆⋃j∈JL(r)Gj⊆⋃s=1rIL(s).
This, together with (5.22) and the fact that the sets Ip
are disjoint,
yields
[TABLE]
We first assume that the inclusion in (5.25) is strict. Then
[TABLE]
Indeed, if (5.26) is false, then JL(r)⊆⋃p=1lD(p)(t), and hence IL(r)⊆⋃p=1lIp. The latter inclusion, together with (5.22) and the fact that the sets Ip are disjoint, yields IL(r)⊆⋃p=pr−1+1lIp, contrary to the case assumption, so (5.26) follows.
By
(2.1), (3.7)-(3.8), (5.16),
(5.22)-(5.23) and (5.25), for small t>0, f(l+1)(t) is the supremum of x≤t satisfying the
constraints
[TABLE]
[TABLE]
for j∈Kl.
Recall that in the corresponding linear problem, for small t>0,
[TABLE]
for j∈KL(r−1), with equality for j∈JL(r) (compare (5.9), (5.1)). Note that the second equality in (5.28) follows from (5.25). Comparing (5.27) to (5.28) and using (5.24), (5.26), we get
the inclusion Jl+1⊆JL(r) and
(5.23) for p=l+1, s=r.
This ends the inductive step in the case of strict inclusion in (5.25).
It remains to analyze the case in which
[TABLE]
Then (5.22) holds for s=1,...,r, so
⋃p=1lIp=⋃s=1rIL(p) and
Kl=KL(r) (this is the equality in (5.24), with r in the place of r−1). Also, (5.22) for s≤r, together with the inequality l<k, implies that r<kmaxL.
The counterpart of (5.27) in this case is
[TABLE]
for j∈Kl. If r+1=kmaxL, ym(r+1)=ym∗=0 and a(r+1)<1, then fL(r+1(t)=t for small t>0 and the counterpart of (5.28) is
[TABLE]
for j∈KL(r)=Kl. Comparing (5.30) to (5.31) we see that f(l+1)(t)=t and l+1=kmax(t)=k, so our inductive proof of (5.17) is complete.
If ym(r+1)<0 or a(r+1)≥1, then the counterpart of (5.28) is
[TABLE]
for j∈KL(r)=Kl, with equality for j∈JL(r). Comparing (5.30) to (5.32), we
get (5.23) for p=l+1, s=r+1, and Jl+1⊆JL(r+1), yielding (5.21) for s=r+1, and the proof of
the inductive step is complete.
5.4 Lack of global monotonicity
It is not hard to prove that, without any additional assumptions on the functions hi, (5.1) holds for every 0≤t<t~ and i∈I(1)(t). In spite of this, the mapping F is, in general, not monotone on [0,∞), even if hi, i∈I, are given by (5.2), as the following example shows.
Example 5.5
Let I=7, J=4, G1={1,3,6}, G2={2,4,7}, G3={3,4,5,6,7} and G4={6,7}. Next, let
[TABLE]
so that x1∗=x2∗=x3∗=x4∗=−11, x5∗=0, x6∗=x7∗=−10.
One may easily check that for 0≤t≤8, J(1)(t)={1,2}, J(2)(t)={3},
I(1)(t)={1,2,3,4,6,7}, I(2)(t)={5}, N(1)(t)={4}
and N(2)(t)=\mathchar31. Moreover,
[TABLE]
hence
[TABLE]
and the mapping F fails to be nondecreasing.
The “network topology” in the above example may be somewhat
simplified, at the price of making some of the functions hi(x), x≥xi∗, nonlinear. Namely,
let I=5, J=3, G1={1,3}, G2={2,4} and G3={3,4,5}. Assume (5.33)-(5.34) and let
h3(x)=h4(x)=(x+11)++2(x+10)+, so that x1∗=x2∗=x3∗=x4∗=−11, x5∗=0.
(Somewhat informally, this network structure has been obtained from the previous one by removing the fourth server and merging the routes 3,6 (resp., 4,7) into a single route 3 (resp., 4).
It is easy to verify that in this case for 0≤t≤8, J(1)(t)={1,2}, J(2)(t)={3},
I(1)(t)={1,2,3,4}, I(2)(t)={5} and (5.35)-(5.37) still hold.
Note that this network satisfies the local traffic condition and hence
N(1)(t)=N(2)(t)=\mathchar31 for all t.
**
6 Further research directions
The results obtained in this paper may be regarded as
introductory in nature and there are several important issues regarding our mapping F that remain to be addressed.
First, one would like to relax the assumptions on regularity of hi necessary for local monotonicity of F. Example 5.5 shows that
“kinks” of hi may create problems in this regard, so it is not immediately clear that the monotonicity result of Section 5.3 may be carried over even to the Lipschitz case. A remedy for this problem may be creating a “differential” version of the algorithm from Section 3.1, determining the derivatives of f(p), rather than their values, for Lipschitz hi, in a way similar to our analysis for the linear case. This, however, in the absence of C1 regularity of hi, yields an ODE system with discontinuous right-hand side, so even establishing existence of solutions to such a system may be challenging.
Another direction that appears to be important for applications is to skip the assumption of strict monotonicity of hi. As Example 4.1 indicates, this results in jumps of the corresponding process F. However, from the point of view of the queueing application described in Section 2.3, with the functions hi given by (2.9), this is not necessarily a problem, because some Fi may “jump over the flat spots”, containing no mass of the corresponding initial distributions Zi(0,⋅), without causing discontinuity of the resulting locally edge minimal fluid model.
Similarly, it may be useful to investigate functions hi with upward jumps, corrresponding to distributions with atoms. This would open an avenue to using techniques similar to those developed in our forthcoming paper, but for pre-limit stochastic networks, rather than
for the corresponding fluid limits.
Example 5.5 shows that there is no hope for global monotonicity of F in the general case. However, it is plausible that for some simple network topologies (e.g., linear or tree networks), the mapping F is monotone on [0,∞). This would greatly simplify the analysis of the corresponding fluid limits, and aid the investigation of the pre-limit stochastic networks.
Finally, it may be interesting to replace the relation “⪕” and/or the set At in the definition of F(t) by a different partial ordering and/or admissible set, and to investigate properties and possible applications of the resulting mappings.
7 Appendix: Inductive proof of Theorem 4.2
We continue the argument starting in Section 4. As we have already explained, we may assume that t0>0 and that (4.3), (4.6) hold.
7.1 The base case: i∈I(1)(t0)
For j∈J(1), we have
∑i∈Gjhi(f(1)(tn))=tn, so ∑i∈Gjhi(f(1)(t0))=t0.
Consequently, J(1)⊆J(1)(t0), and thus I(1)⊆I(1)(t0).
Let i∈I(1). Then i∈I(1)(t0), so, by (3.10),
Fi(tn)=f(1)(tn)∨xi∗→f(1)(t0)∨xi∗=Fi(t0)
as n→∞.
We have shown that for i∈I(1), we have (4.2).
If I(1)=I(1)(t0) (in particular, if I(1)=I), we have (4.2) for every i∈I(1)(t0). Assume that I(1)=I(1)(t0). Then
[TABLE]
because ⋃j∈D(1)Gj=I(1) (see
(4.1)), while ⋃j∈J(1)(t0)Gj=I(1)(t0). Let
[TABLE]
By (7.1), we have p1≥2. For p=1,...,p1, let
[TABLE]
Clearly,
[TABLE]
and ⋃p=1p1Bp(1)⊆I(1)(t0). We will check
that
[TABLE]
Indeed, let i∈I(1)(t0). Then i∈Gj for some
j∈J(1)(t0), so, by (7.4), j∈Ap(1)⊆D(p) for some p∈{1,...,p1}.
Consequently, for some q∈{1,...,p}, we have i∈I(q) and hence i∈Bq(1). Thus, I(1)(t0)⊆⋃p=1p1Bp(1) and
(7.5) follows. Moreover, the above argument justifies
the inclusion
[TABLE]
where the equality follows from the inclusion Ap(1)⊆D(p).
Let p∈{2,...,p1} be such that Ap(1)=\mathchar31 and
let j∈Ap(1). Since j∈J(1)(t0) and (4.6) holds,
we have
[TABLE]
If f(p)(tn)=tn for infinitely many n, then,
by the definition of f(p) (compare
(3.7)-(3.8)) and the inequality f(1)<f(p), for these n we have
[TABLE]
Letting n→∞ in (7.1), we get
[TABLE]
The relations (7.7) and (7.9), together with
the inequalities f(1)(t0)≤t0, t0>0≥maxi∈Ixi∗, imply that f(1)(t0)=t0, which contradicts
(4.6). Consequently, since Ap1(1)=\mathchar31 by
(7.2)-(7.3),
in the remainder of our proof we may assume that for p≤p1,
[TABLE]
Since j∈Ap(1)⊆D(p)=J(p)∪N(p)⊆K(p−1), for every n we have
[TABLE]
(compare (3.6)-(3.7)), and
[TABLE]
with equality for j∈Cp(1):=J(p)∩J(1)(t0)⊆Ap(1). Note that (7.11) and the
above-mentioned equality in (7.12) imply that for n≥1,
j∈Cp(1),
[TABLE]
By (4.3) and (7.13),
[TABLE]
Letting n→∞ in (7.12) for j∈Ap(1) and using (4.3), we get
[TABLE]
This, together with (7.7), (4.4) and monotonicity
of hi, yields the equality in (7.15) for each j∈Ap(1) and the equations
[TABLE]
By (7.6) and (7.16) with p=p1, we have
[TABLE]
The relations (4.4), (7.13)-(7.14), (7.16) (with k=p) and the properties of the functions hi, in turn, imply that either
[TABLE]
or
[TABLE]
provided that Cp(1)=\mathchar31.
In the latter case, fj,p∗ does not depend on j∈Cp(1) and hence it will be denoted by fp∗.
The above argument implies that if
f(1)(t0)≥minj∈Cp(1)fj,p∗ (or if fj,p∗=fj′,p∗ for some j,j′∈Cp(1)),
then
(7.19) cannot hold, and hence (7.18) holds.
Similarly, if Cp(1)=\mathchar31 and
[TABLE]
then, by (7.14) and
(7.19),
[TABLE]
Let
pˉ1=max{p∈{1,...,p1}:f∞(p)=f(1)(t0)}. We have shown that
[TABLE]
and that if pˉ1<p1, then for every p>pˉ1 such that Cp(1)=\mathchar31, the relations (7.20)-(7.21) hold.
It may be the case that
pˉ1<p1 and
for some p~∈{pˉ1+1,...,p1}, we have Cp~(1)=\mathchar31. In this case, if i∈Bp~(1), then
the relations (4.4) and
(7.17)
imply that
[TABLE]
Fix i∈I(1)(t0).
If i∈Bp(1) for some p≤pˉ1, then, by (3.10) and (7.22),
[TABLE]
so (4.2) holds.
If i∈Bp(1) for some p>pˉ1, then there exists k∈{p,...,p1} such that i∈Gj∩I(p) for some j∈Ak(1) (see (7.6)). If Cp(1)=\mathchar31, then, by (7.21), f∞(p)=fp∗ and, by the same argument as in (7.23),
[TABLE]
Thus, since i∈I(1)(t0),
by (2.1), (3.10), (7.21) and (7.24),
as n→∞,
[TABLE]
and again (4.2) holds.
Finally, if Cp(1)=\mathchar31, then reasoning as in (7.25), but using (7.23) instead of (7.24), we get
limn→∞Fi(tn)=Fi(t0).
Since, by (7.5), every i∈I(1)(t0) belongs to some Bp(1), p≤p1, the above argument shows that (4.2) is true for all i∈I(1)(t0).
If D(1)(t0)=J, then I(1)(t0)=I and (4.2) holds for each i∈I. In what follows, we assume that
K(1)(t0)=\mathchar31.
For future reference note that
[TABLE]
Indeed, let j∈J(p) for some p≤pˉ1. By (7.10), ∑i∈Gjhi(Fi(tn))=∑k=1p∑i∈Gj∩I(k)hi(f(k)(tn))=tn, which, together with (7.22), implies that
∑i∈Gjhi(f(1)(t0))=∑k=1p∑i∈Gj∩I(k)hi(f(1)(t0))=t0,
and hence j∈J(1)(t0). We also have
⋃p=1pˉ1D(p)⊆D(1)(t0). Indeed, (7.26) implies ⋃p=1pˉ1I(p)=⋃p=1pˉ1⋃j∈J(p)Gj⊆⋃j∈J(1)(t0)Gj=I(1)(t0).
Thus, for j∈D(p), p≤pˉ1, we have Gj⊆⋃k=1pI(k)⊆I(1)(t0), and hence j∈D(1)(t0).
7.2 Inductive assumption
Fix m∈{1,...,kmax(t0)−1}. For l=1,...,m, let bl, pl be given by (4.7)-(4.8).
By definition, b1=1 and bl≤pl for each l=1,...,m. For l=1,...,m−1,
let ql=bl+1−1, so that by (4.7) we have
[TABLE]
For p=bl,...,pl, l=1,...,m, let
[TABLE]
The definitions (4.7)-(4.8) imply that
[TABLE]
(see also (7.27)).
Using (7.27) and (7.29), one may easily check
(4.9).
Suppose that for l=1,...,m, the following assertions hold.
There exist indices pˉl∈{bl,...,pl} such that (4.10) holds and
[TABLE]
Note that (4.7) and (7.30) imply
[TABLE]
If pˉl<pl, then for every p∈{pˉl+1,...,pl}
such that Cp(l)=\mathchar31, we
have
[TABLE]
where fp∗ is the common value of fj,p∗ (defined in
(7.13)) for j∈Cp(l) and, moreover,
fp∗≤xi∗ for i∈Bp(l).
If
Cp~(l)=\mathchar31 for some p~∈{pˉl+1,...,pl}, then
[TABLE]
For j∈Ap(l), p=bl,...,pl and r=1,..,l,
[TABLE]
For r=1,...,m and k=br,...,pr, Bk(r)⊆I(r)(t0)⊆⋃j∈J(r)(t0)Gj=⋃p=brpr⋃j∈Ap(r)Gj,
where the equality follows from (7.29). Hence, (7.34) implies
[TABLE]
Because f(1)(t0)<...<f(l)(t0)≤f∞(k) for k=bl,...,pl,
(7.34) implies also that for p=bl,...,pl,
[TABLE]
Finally, we assume that (4.2) holds for all i∈⋃k=1mI(k)(t0).
Note that all the above assumptions have already been verified in the case of m=1.
7.3 Inductive step
Define bm+1 by (4.7) with l=m+1 and let qm=bm+1−1.
By definition, we have (7.27) for l=m,
so K(m)(t0)⊆K(qm).
The inclusion
(7.30) implies that
[TABLE]
Also, for j∈K(qm) (and hence for j∈K(m)(t0)), we have
[TABLE]
If
[TABLE]
for some n,
then kmax=bm+1,
J(bm+1)=K(qm), N(bm+1)=\mathchar31 and Fi(tn)=tn for all
i∈I(bm+1)=I∖⋃k=1qmI(k). If (7.39) holds for infinitely many n, then
letting n→∞ in (7.38) along this subsequence, by (7.27) for l=m, monotonicity of hi, the inequality Fi(t0)≤t0 for all i∈I, and the inductive assumption (4.2) for i∈⋃k=1mI(k)(t0), we get
[TABLE]
for all j∈K(m)(t0). This, in turn, implies that Fi(t0)=t0 for all i∈I(m+1)(t0)=I∖⋃p=1mI(p)(t0)⊆I∖⋃k=1qmI(k), where the inclusion follows from the second inclusion in (7.27) for l=m. Hence, for all i∈I(m+1)(t0) (and thus for all i∈I), we have (4.2) along a subsequence {tn} satisfying (7.39).
Consequently, in the remainder of the proof we may assume that
f(bm+1)(tn)<tn for each
n≥1, and thus equality holds in (7.38)
for j∈J(bm+1).
We claim that
[TABLE]
Letting n→∞ in (7.38), we get
[TABLE]
for j∈K(qm), with equality for j∈J(bm+1) (compare (7.40)). For k=1,...,m, by
(4.10) and
(7.37), we have
[TABLE]
so for i∈I(k)(t0),
hi(Fi(t0))=hi(f(k)(t0))≤hi(f∞(bm+1)). This, together with (7.27) for l=m and (7.43) implies
[TABLE]
for j∈K(qm) (hence for j∈K(m)(t0)), and thus
[TABLE]
In order to prove the opposite inequality, we will first check that
[TABLE]
Suppose, to the converse, that
[TABLE]
Then
I(bm+1)⊆⋃j∈J(bm+1)Gj⊆⋃p=1m⋃j∈J(p)(t0)Gj=⋃p=1mI(p)(t0),
and thus, by the second inclusion in (7.27) for l=m,
⋃p=1bm+1I(p)⊆⋃p=1mI(p)(t0). Consequently, for j∈N(bm+1), Gj⊆⋃p=1mI(p)(t0), and hence N(bm+1)⊆⋃p=1mD(p)(t0). This, together with (7.48), implies that D(bm+1)⊆⋃p=1mD(p)(t0), which contradicts the definition of bm+1. We have proved (7.47).
Fix j∈J(bm+1)∖⋃p=1mJ(p)(t0). Then equality holds in (7.43). We will check that equality holds in (7.45) as well.
Suppose that i∈Gj∩⋃p=1mI(p)(t0)∖⋃p=1qmI(p). Then i∈I(k)(t0)∩I(r) for some k≤m and r≥bm+1.
Letting n→∞ in the equality
hi(Fi(tn))=hi(f(r)(tn)), and using the inductive
assumption (4.2) for i∈⋃p=1mI(p)(t0), we get hi(f(k)(t0))=hi(Fi(t0))=hi(f∞(r)), while f(k)(t0)≤f(m)(t0)≤f∞(bm+1)≤f∞(r) by (7.44).
Hence hi(Fi(t0))=hi(f∞(bm+1)), so (7.27) for l=m and the equality in
(7.43) imply that equality holds in
(7.45) as well. Arguing as in (7.13), we get f(m+1)(t0)>min{xi∗:i∈Gj∖⋃p=1mI(p)(t0)}, so (7.46)
and the equality in (7.45) imply
that f∞(bm+1)=f(m+1)(t0).
We have proved
(7.41). Replacing f∞(bm+1) by f(m+1)(t0) in the equality in (7.45) we get j∈J(m+1)(t0), so (7.42) holds as well.
Let i∈I(bm+1)∖⋃p=1mI(p)(t0).
We will show that (4.2) holds. Since i∈I(bm+1), we have i∈Gj for some j∈J(bm+1). If j∈⋃p=1mD(p)(t0), then i∈⋃p=1mI(p)(t0), contrary to the choice of i, so j∈J(bm+1)∖⋃p=1mD(p)(t0)⊆J(bm+1)∖⋃p=1mJ(p)(t0) and hence, by (7.42), j∈J(m+1)(t0).
Then Gj⊆⋃p=1m+1I(p)(t0), so i∈I(m+1)(t0). Consequently, by (3.10) and (7.41),
[TABLE]
as n→∞, so (4.2) holds.
If I(m+1)(t0)⊆I(bm+1), we have (4.2) for each i∈I(m+1)(t0). Assume that
[TABLE]
and
let pm+1 be given by (4.8) with l=m+1. By (7.27) for l=m and (7.50), bm+1<pm+1. For p=bm+1,...,pm+1, define the sets Ap(m+1), Bp(m+1), Cp(m+1) by (7.28) with l=m+1. By the definitions of bm+1, pm+1 and (7.27) for l=m, we have (7.29) for l=m+1.
Clearly, ⋃p=bm+1pm+1Bp(m+1)⊆I(m+1)(t0). We will check
that (4.9) holds
for l=m+1.
Indeed, let i∈I(m+1)(t0). Then i∈Gj for some
j∈J(m+1)(t0), so, by (7.29) with l=m+1, j∈Ap(m+1)⊆D(p) for some p∈{bm+1,...,pm+1}.
Consequently, we have i∈I(q) for some q≤p. However, q>qm by the second inclusion in (7.27) with l=m, and hence i∈Bq(m+1) for some q∈{bm+1,...,pm+1}. Thus, I(m+1)(t0)⊆⋃p=bm+1pm+1Bp(m+1) and (4.9)
with l=m+1 follows.
We will first consider the case in which f(m+1)(t0)=t0. Then I(m+1)(t0)=I∖⋃p=1mI(p)(t0) and Fi(t0)=t0 for each i∈I(m+1)(t0). Moreover, by (7.41),
[TABLE]
Let i∈I(m+1)(t0). Then i∈I(p) for some p≥bm+1 by the second inclusion in (7.27) with l=m, and thus
f(bm+1)(tn)∨xi∗≤Fi(tn)≤tn,
which,
by (7.51), implies that Fi(tn)→t0=Fi(t0),
so (4.2) holds for all i.
It remains to consider the case of
[TABLE]
Let p∈{bm+1,...,pm+1} be such that Ap(m+1)=\mathchar31 and
let j∈Ap(m+1). Since j∈J(m+1)(t0) by (7.28) with l=m+1, the relation (7.52)
implies
[TABLE]
Suppose that f(p)(tn)=tn for infinitely many n, and thus p=kmax.
Then N(p)=\mathchar31, so J(p)=D(p) and Ap(l)=Cp(l) for each l≤m+1.
If pl=kmax
for some l≤m, then pl=p≥bm+1, and Cpl(l)=Apl(l)=\mathchar31 by (4.8), (7.28).
Hence,
by (7.31), (7.37), (7.32), as n→∞, we have
0<t0←tn=f(pl)(tn)→fpl∗≤0, contradiction. Thus, p>max{pl,l≤m}, so by (7.28) and (4.9),
I(p)∩⋃l=1mI(l)(t0)=\mathchar31.
Consequently, for j∈Ap(m+1)⊆J(m+1)(t0), Gj∩I(p)⊆I(m+1)(t0) and the definition of f(p) (compare
(3.8)), together with (4.9) for l=1,...,m+1, imply
[TABLE]
Letting n→∞ and using (4.10), (7.41), we get
[TABLE]
which, together with (7.53), implies that f(m+1)(t0)=t0, contrary to (7.52).
Consequently, since Apm+1(m+1)=\mathchar31 by (4.8) and (7.28) with l=m+1, in the remainder of the proof we may assume that (7.10) holds for p=bm+1,...,pm+1.
This, in turn, implies (7.11), (7.13) and the equality in (7.12) for j∈Cp(m+1), p∈{bm+1,...,pm+1}, by the same argument as in the case of j∈Cp(1), p≤p1.
The relation
(7.13) implies
[TABLE]
For j∈Ap(m+1), p=bm+1,...,pm+1,
letting n→∞ in (7.12),
we get (7.15).
This, in turn, together with (7.53), (4.4), (4.7)-(4.8), (4.10), (7.41) and monotonicity of hi,
yields (7.34) for j∈Ap(m+1),
p=bm+1,...,pm+1 and r=1,..,m+1. Consequently,
(7.36) holds with l=m+1 and p=bm+1,...,pm+1.
Using (7.54) and (7.36) (with l, p as above) instead of (7.14) and (7.16), respectively, in analogy to the development of
(7.18)-(7.23),
(7.24), we can derive the following facts.
Let
[TABLE]
Then (4.10) holds with l=m+1.
If pˉm+1<pm+1, then for every p∈{pˉm+1+1,...,pm+1} such that Cp(m+1)=\mathchar31, we have (7.32) with l=m+1,
where fp∗ is the common value of fj,p∗ for j∈Cp(m+1). Moreover, for such p,
[TABLE]
If
Cp~(m+1)=\mathchar31 for some p~∈{pˉm+1+1,...,pm+1}, then (7.33) holds with l=m+1.
Let i∈I(m+1)(t0). We will check that (4.2) holds. By (4.9)
with l=m+1, i∈Bp(m+1) for some p∈{bm+1,...,pm+1}. If p≤pˉm+1,
then (4.2) follows from (4.10) with l=m+1 by an argument similar to (7.49). If pˉm+1<p≤pm+1, then
the argument is analogous to the proof of the corresponding case
pˉ1<p≤p1 for i∈I(1)(t0) (see
(7.25) and the paragraph surrounding it); we just use
(7.32) with l=m+1, (7.55),
(7.33) with l=m+1, instead of (7.21),
(7.24), (7.23), respectively.
If m+1=kmax(t0), then (4.2) holds for all i∈I and we are done. Assume that K(m+1)(t0)=\mathchar31. For the sake of the next inductive step, we will show
(7.30) for l=m+1.
Let j∈J(p)∖⋃k=1mJ(k)(t0) for some bm+1≤p≤pˉm+1. Since (7.10) holds for this p, we have
tn=∑k=1p∑i∈Gj∩I(k)hi(f(k)(tn)). Letting n→∞ and using (7.28), (4.9), (4.7) with l=m+1, (7.35) and (4.10) with l=m+1, we get
[TABLE]
so j∈J(m+1)(t0). We have shown that
⋃p=bm+1pˉm+1J(p)⊆⋃p=1m+1J(p)(t0), so
[TABLE]
Let j∈D(p) for some bm+1≤p≤pˉm+1. By (7.27) for l=m and (7.56),
[TABLE]
Thus, j∈⋃p=1m+1D(p)(t0)
and (7.30) with l=m+1 follows.
In summary, we have shown that either (4.2) holds for all i∈I, or all the assertions of the inductive assumption, including (4.2) for i∈I(l), are true for l=1,...,m+1. This ends the inductive proof of (4.2) for each i∈I.