UvA-DARE (Digital Academic Repository)
Extremes and fluid queues
Dieker, A.B.
Publication date 2006 Document Version Final published version
Link to publication
Citation for published version (APA): Dieker, A. B. (2006). Extremes and fluid queues. Ponsen & Looijen.
General rights It is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), other than for strictly personal, individual use, unless the work is under an open content license (like Creative Commons).
Disclaimer/Complaints regulations If you believe that digital publication of certain material infringes any of your rights or (privacy) interests, please let the Library know, stating your reasons. In case of a legitimate complaint, the Library will make the material inaccessible and/or remove it from the website. Please Ask the Library: https://uba.uva.nl/en/contact, or a letter to: Library of the University of Amsterdam, Secretariat, Singel 425, 1012 WP Amsterdam, The Netherlands. You will be contacted as soon as possible.
UvA-DARE is a service provided by the library of the University of Amsterdam (https://dare.uva.nl)
Download date:30 Sep 2021 Extremes and uid queues Extremes and uid queues / Antonius Bernardus Dieker, 2006 Proefschrift Universiteit van Amsterdam Met lit. opg. – Met samenvatting in het Nederlands Omslagontwerp door Tobias Baanders Gedrukt door Ponsen & Looijen BV ISBN 90-5776-151-3
Dit onderzoek kwam tot stand met steun van NWO Extremes and uid queues
ACADEMISCH PROEFSCHRIFT
ter verkrijging van de graad van doctor aan de Universiteit van Amsterdam op gezag van de Rector Magni cus prof. mr. P. F. van der Heijden ten overstaan van een door het College voor Promoties ingestelde commissie, in het openbaar te verdedigen in de Aula der Universiteit op donderdag 9 maart 2006, te 10.00 uur
door
Antonius Bernardus Dieker
geboren te Amsterdam Promotiecommissie:
Promotor: prof. dr. M. Mandjes
Overige leden: dr. A. A. Balkema prof. dr. ir. O. J. Boxma prof. dr. C. A. J. Klaassen prof. dr. I. Norros prof. dr. T. Rolski dr. ir. W. R. W. Scheinhardt dr. P. J. C. Spreij
Faculteit der Natuurwetenschappen, Wiskunde en Informatica Dankwoord
Dit boekje is het resultaat van vier jaar zweten. Toen ik aan dit project begon kon ik me nauwelijks voorstellen wat het schrijven ervan allemaal zou behelzen, maar het uiteindelijke resultaat vervult mij met trots. Hoewel alleen mijn naam op de omslag staat, had ik het niet kunnen schrijven zonder de hulp van velen. Graag richt ik het woord tot enkele personen in het bijzonder, en haal enkele speciale herinneringen naar boven. Michel, jij bent zonder twijfel de persoon die de meeste invloed heeft gehad op mijn werk; de hoofdstukken 8, 9, 10 en 14 hebben wij samen geschreven. Niet alleen deze samenwerking heb ik als zeer prettig ervaren, maar ook het vertrouwen dat jij altijd in mij gesteld hebt. De deur stond en staat altijd open voor een uitgebreide discussie over uiteenlopende onderwerpen. Veel dank hiervoor. Dit proefschrift heeft ook erg gepro teerd van vele discussies met collega’s uit binnen- en buitenland, waarvan ik hier slechts enkelen kan noemen. Hoofdstuk 13 is gezamenlijk werk met Krzy s en Tomasz; jullie ben ik in het bijzonder erkentelijk. Bert, jij hebt mij altijd van goede adviezen voorzien. Misja, bedankt voor je ongezouten kritiek op een deel van het manuscript. Verder gaat mijn dank uit naar Thyra en Richard voor alle Twentse inspanningen en het instituut EURANDOM voor de Brabantse gastvrijheid. Ik kan niet aan de afgelopen jaren denken zonder dat mooie herinneringen naar boven komen aan de plaatsen waar dit werk mij heeft gebracht. Ik denk daarbij niet alleen aan de grote conferenties in Brazili e en Canada, maar ook aan de werkbezoeken aan Polen en Ierland (waarvoor dank, Krzy s, Tomasz en Neil). Een langere tijd heb ik in Parijs doorgebracht, op deze mogelijkheid geattendeerd door Peter. In Parijs heb ik de basis gelegd voor een artikel met Marc [109]. Ren e, van alle kamergenoten neem jij een speciale plaats in. Vier jaar lang wist jij het uit te houden tussen mijn ongeloo ijke troep, waarvoor respect op z’n plaats is. Na ruim twee jaar werd onze kamer aangevuld met Regina; enigszins verrast door de saaie werkomgeving van twee mannen leukte jij de kamer al snel op met een paar planten. Lo c, jou wil ik bedanken voor jouw geduld als ik Frans probeerde te praten. Mijn familie heeft mij van heel dichtbij zien worstelen aan dingen die ik hen nooit precies heb kunnen uitleggen. Pap, mam, Jos, Miranda en Marjon, heel erg bedankt voor alle steun. Jasper, jij weet altijd een lach op m’n gezicht te toveren. Jullie weten dat de afgelopen vier jaar niet altijd zo makkelijk was als dit gelikte boekje doet vermoeden. De wiskundige stellingen en talloze formules uit dit proefschrift verbleken bij mijn grootste ontdekking van de afgelopen vier jaar: Cyndi. Lieverd, ik weet niet hoeveel punten ik je moet geven voor alles wat je voor me hebt gedaan! Ik ben blij dat we samen van Ierland kunnen gaan genieten.
Amsterdam, januari 2006 Ton
Contents
1 Introduction 1 1.1 The uid queue ...... 3 1.2 Fluid networks ...... 6 1.3 Approximations for the input process ...... 8 1.4 Outline of the thesis ...... 10
2 Techniques 13 2.1 Regular variation ...... 13 2.2 Weak convergence ...... 15 2.3 Large deviations ...... 17 2.4 Tail asymptotics ...... 20
A Gaussian queues 23
3 Background on Gaussian processes 25 3.1 Gaussian measures and large-deviation principles ...... 26 3.2 Two fundamental inequalities ...... 28 3.3 The double sum method ...... 29 3.4 Outline of Part A ...... 32
4 Conditional limit theorems 33 4.1 Introduction ...... 33 4.2 Queueing results ...... 36 4.3 Preliminaries ...... 38 4.4 Weak convergence results ...... 41 4.5 Proofs for Section 4.2 ...... 47 4.6 Proof of Proposition 4.5 ...... 48
5 Extremes of Gaussian processes 51 5.1 Introduction ...... 51 5.2 Description of the results and assumptions ...... 53 5.3 Special cases: stationary increments and self-similar processes ...... 57 5.4 A variant of Pickands’ lemma ...... 64 5.5 Four cases ...... 67 5.6 Upper bounds ...... 73 5.7 Lower bounds ...... 77 iv Contents
6 Reduced-load equivalence 85 6.1 Introduction ...... 85 6.2 Description of the result ...... 87 6.3 Examples ...... 88 6.4 Proof of Theorem 6.2 ...... 89
B Simulation 93
7 Background on simulation 95 7.1 Importance sampling and asymptotic e ciency ...... 96 7.2 Examples: random walks ...... 98 7.3 Outline of Part B ...... 101
8 Conditions for asymptotic e ciency 103 8.1 Introduction ...... 103 8.2 The Varadhan conditions for e ciency of exponential twisting ...... 104 8.3 Comparison with Sadowsky’s conditions ...... 107 8.4 Uniqueness of the exponential twist ...... 110 8.5 An example ...... 111
9 Random walks exceeding a nonlinear boundary 113 9.1 Introduction ...... 113 9.2 Sample-path large deviations and exceedance probabilities ...... 114 9.3 Path-level twisting ...... 115 9.4 Step-level twisting ...... 120 9.5 A comparison ...... 124
10 Simulation of a Gaussian queue 127 10.1 Introduction ...... 127 10.2 The bu er-content probability ...... 129 10.3 Simulation methods ...... 133 10.4 Evaluation ...... 139 10.5 Concluding remarks ...... 144 10.A Appendix: proofs ...... 145
C L evy-driven uid queues and networks 149
11 Background on L evy processes 151 11.1 Random walks ...... 151 11.2 L evy processes ...... 155 11.3 Outline of Part C ...... 157
12 Factorization embeddings 159 12.1 Introduction ...... 159 12.2 On factorization identities ...... 161 12.3 L evy processes with phase-type upward jumps ...... 164 12.4 Perturbed risk models ...... 166 12.5 Asymptotics of the maximum ...... 173 Contents v
13 Quasi-product forms 179 13.1 Introduction ...... 179 13.2 Splitting times ...... 181 13.3 The P-distribution of (X, F ) ...... 184 P 13.4 The k↓-distribution of H ...... 186 13.5 Multidimensional Skorokhod problems ...... 189 13.6 Tandem networks and priority systems ...... 194 13.A Appendix: a compound Poisson process with negative drift ...... 200
14 Extremes of Markov-additive processes 205 14.1 Introduction ...... 205 14.2 A discrete-time process and its extremes ...... 207 14.3 Markov-additive processes and their extremes ...... 223 14.4 The uid queue: theory ...... 233 14.5 The single queue: examples ...... 237 14.6 Tandem networks with Markov-additive input ...... 242 14.7 Concluding remarks ...... 244
Bibliography 247
Author Index 261
Samenvatting (Summary) 265
About the author 269
CHAPTER 1
Introduction
In queueing theory, the main entity is a station where service is provided; examples include counters, call centers, elevators, and tra c lights. Customers seeking this service arrive at the system, where they wait if the service facility cannot immediately allocate the required amount of service. They depart after being served. In most of these applications, it is intrinsically uncertain at what time customers arrive and how long they need to be served. This explains why probability theory plays an important role in the analysis of queueing systems. In order to design these systems optimally, it is desirable to know how this randomness in- uences the performance of the system. Therefore, it is crucial to study system characteristics that re ect this performance. For instance, one may wish to analyze the probability distribu- tions of the queue length or the waiting time of a customer, as a function of the random arrival process and service requirements. The investigation of such a performance measure starts with a reformulation of the problem into mathematical terms. On an appropriate level of abstraction, the analysis no longer involves queues; rather a problem of (applied) probability theory needs to be solved. The system characteristics that we encounter in this thesis are all related to so-called extremes as a result of this translation. An inherent advantage of investigating queues through examining extremes is the wide applicability of the resulting theory. For instance, many of our results are also relevant for risk theory and nancial mathematics; even though it may be unclear upfront how these elds relate to customers waiting in a line, extremes also play a pivotal role in these theories. Since our results are often illustrated with queueing examples, we rst brie y discuss how queues are related to extremes.
Extremes
To see how the connection between queues and extremes arises, we set V0 := 0, and consider Lindley’s recursion V := max(0, V + B r), (1.1) n+1 n n for n 0, where r > 0 is given and the Bn are nonnegative. This recursion plays a role in a variet y of queueing situations, for instance when analyzing waiting times and the remaining amount of work in a queueing system. In those situations, the sequence B := B : n 0 is { n } random, and V := Vn : n 0 inherits this property; we then say that B and V are stochastic processes. { } In the absence of the maximum in (1.1), the solution of this recursion can immediately be 2 Chapter 1 Introduction given. Indeed, if V˜ := V˜ + B r for n 0 while V˜ = 0, then we have for n 1, n+1 n n 0 n 1 V˜ = B rn. (1.2) n ` X`=0
We call V˜ := V˜n : n 0 the free process. Before solving the recursion with the maximum, it is insightful to{interpret this} maximum in terms of the free process. As long as the constrained process V := Vn : n 0 is strictly positive, it behaves exactly in the same way as the free process. How{ever, if the }free process becomes negative, the constrained process is ‘pushed’ back to zero. The process Vn : n 0 is called the re ected process; occasionally, the term regulated process is encountered{ in the literature.} In Section 1.1.2, we apply the recursion in (1.1) n times to obtain the identity
Vn = V˜n min V˜k. (1.3) 0 k n
For many underlying stochastic processes Bn , V˜n min0 k n V˜k has the same distribution as { } max0 k n V˜k. This constitutes the remarkable (and perhaps even counterintuitive) fact that the re ected process is equal in distribution to the running maximum of the free process. In other words, the probability that Vn exceeds x equals the probability that V˜ reaches x no later than time n. As a result, there are essentially two ways to analyze the solution of Lindley’s recursion. A ‘direct’ approach, which is based on the queueing interpretation, studies the evolution of the re ected process V . Alternatively, one can analyze the running maximum of V˜ ; this approach is based on the extremes of the free process V˜ . This shows that extremes and queues are closely related. Throughout, we use the term ‘extreme’ to refer to the maximum (or minimum) of a random function. We are also interested in other quantities, such as the epoch at which the maximum is attained. These objects, as well as the closely related hitting and rst-passage times, have been investigated ever since the foundations of modern probability theory. In this broader context, the term uctuation theory is sometimes used.
Fluid queues In several applications of queueing theory, individual customers are so small that they can hardly be distinguished. Instead of customers, it is then easier to imagine a continuous stream of work that ows into the system. The resulting queueing models are called uid queues. Fluid queues are closely related to dams. A dam can be modeled as a reservoir, in which water builds up due to rainfall, is temporarily stored, and then released according to some release rule. Consequently, a uid queue can be viewed as a dam in which work is bu ered until enough capacity becomes available. Alternatively, the water in the dam can be interpreted as the amount of goods stored, and the name storage model is therefore sometimes used. Modern communication networks provide further motivation for studying uid queues. In such networks, small data packets are sent to routers (or switches), where they are queued up, subsequently inspected, and then forwarded to (other routers closer to) their destination. For network design purposes it is desirable to gain insight into the amount of work that builds up in router bu ers, rather than the individual waiting times of the packets. Instead of interpreting packets as customers, the aggregate packet ow can thus be viewed as uid. Technically, shifting from standard queues to uid queues amounts to working with a continuous-time version of Lindley’s recursion. Even though a considerable body of theory is needed to specify this ‘continuous recursion’, the free and re ected processes are simply continuous-time analogues of (1.2) and (1.3) respectively. This can be intuitively understood 1.1 The uid queue 3 by imagining increasingly smaller time units. Most importantly, a distributional duality again often relates the re ected process to the running maximum of its free counterpart. This chapter is organized as follows. First, in Section 1.1, we describe uid queues and their relationship to extremes of stochastic processes in more detail. Then, a brief introduction to uid networks is given in Section 1.2. Finally, Section 1.3 motivates the investigation of several special uid queues that are encountered in this thesis.
1.1 The uid queue
1.1.1 Model description
Let A := At : t 0 be a continuous-time stochastic process such that for any t 0, At is the amoun{t of work o ered} to the system in the interval [0, t]. Throughout this thesis, we use At and A(t) interchangeably; this notational convention is employed for all continuous-time stochastic processes. We suppose that A has right-continuous sample paths with left-hand limits, which is often summarized by saying that A has c adl ag sample paths. The bu er can be interpreted as a uid reservoir, to which input is o ered according to the input process A. The bu er is drained at a constant rate r, i.e., a tap at the bottom of the uid reservoir releases uid at rate r as long as the bu er is nonempty. After the uid is processed, it immediately leaves the system. Throughout, we suppose that the bu er capacity is unlimited. Although its name suggests that sample paths of the input process are always nondecreasing, we do not impose this condition. Indeed, if the input (over a certain time interval) is negative while the bu er stays nonempty, this input is interpreted as additional out ow. In this way, nonconstant (e.g., random) release rules can be incorporated into the model. On the other hand, since it is impossible to drain uid from an empty reservoir, nothing happens if the bu er is empty and the system dynamics dictate that uid be released. We write Wt for the amount of work in the bu er at epoch t, and call this the bu er content. The bu er-content process is also known as a (stochastic) storage process. For simplicity, we suppose that there is initially no work in the system, i.e., that W0 = 0. The uid queue naturally arises from the standard queueing model with individual cus- tomers, since it can be used to represent the un nished work in such a model. Indeed, the input corresponds to service requirements at arrival epochs of customers, and the un nished work declines at unit rate when service is provided. In this context, the bu er content is sometimes referred to as workload or virtual waiting time.
1.1.2 Lindley’s recursion To gain intuition for the general uid queue, we rst consider the special case where work only arrives at the epochs 0, 1, 2, . . .. Let Bn be the (nonnegative) amount of work that arrives at epoch n. Therefore, A is a nondecreasing (random) step function with jump Bn at epoch n. It is our aim to study the bu er content just before a new batch of work arrives: we analyze Vn := Wn . It follows from the description of the model that Vn satis es Lindley’s recursion (1.1). Since we start with an empty system, we set V0 := 0. Lindley’s recursion shows that Vn can be expressed in terms of the Bn:
V = max(0, B r), 1 0 V = max(0, V + B r) = max(0, B r, B + B 2r), 2 1 1 1 1 0 . .
Vn = max(0, Bn 1 r, Bn 1 + Bn 2 2r, . . . , Bn 1 + . . . + B0 rn). 4 Chapter 1 Introduction
The latter equation can be rewritten in terms of the An as
Vn = An 1 rn min(0, A0 r, A1 2r, . . . , An 1 rn), which is Equation (1.3). For a wide class of random step functions A, the distribution of A0, . . . , An 1 equals the { } distribution of An 1 An 2, . . . , An 1 A0, An 1 . In that case, A is called time-reversible, { } and Vn then has the same distribution as
max(0, A0 r, A1 2r, . . . , An 1 rn). (1.4) In queueing theory, one is often interested in the behavior of a system after initial e ects have disappeared, i.e., when the system is in steady-state. This amounts to letting n grow large, but further assumptions are required to ensure that Vn does not blow up. In other words, we need assumptions for stability of the system. In view of (1.4), this means that lim supn An 1 rn < . →∞ ∞ Under this stability assumption, as indicated by (1.4), the random quantity Vn converges in distribution (see De nition 2.9) to
V := sup An 1 rn, (1.5) n 0 where A 1 should be interpreted as zero. If A fails to be time-reversible, a related formula can be giv en by looking backward in time; this concept is often called Loynes’ construction. Equation (1.5) is the main reason why extremes play an important role in the context of queues. There is a vast body of literature on queueing theory in general, and Lindley’s recur- sion in particular. Some textbooks (in alphabetical order) are Asmussen [19], Bene s [38], Borovkov [55], Cohen [79], Kleinrock [191], Prabhu [263], Robert [272], Tak acs [293], and Tijms [295]; a list that is by no means exhaustive. Baccelli and Br emaud [33] is a basic reference on stochastic recursions in queueing theory; for stability issues, we refer to the survey paper by Foss and Konstantopoulos [134]. We also mention Aldous and Bandyopadhyay [9], who give an overview of max-type stochastic recursions in a more general context.
1.1.3 The bu er content Motivated by the above representation of V as the maximum of the free process in case the input process A is a step function, we now focus on the distribution of the bu er content Wt as t for a general input process A. It is our aim to give a continuous-time analogue of (1.5).→ ∞ Since the analysis is much more technical than in the discrete-time case, we do not give a derivation of the resulting formula, but refer to Section 13.5 for more details. We suppose that A is time-reversible; in the continuous-time case, this means that the distributions of At limu (t s) Au : 0 s t and As : 0 s t are equal. The steady-state bu er {conten t then↑ has the same distribution } {as }
sup At rt, (1.6) t 0 a representation that is often attributed to Reich [268]. It is interesting to see the connection between Reich’s formula and the bu er-content process W := Wt : t R on the whole real line. To give the de nition of W , suppose that time is indexed{ by R ∈instead} of R , and for t 0, let A be the amount of work fed into the uid + t c c c 1.1 The uid queue 5 reservoir during the time interval [t, 0]. The bu er content is then also de ned at negative epochs, and Wt can be thought of as the amount of work in the system at time t if the bu er was empty at time . Since W satis es the continuous-time analogue of Lindley’s recursion, ∞ it can be shocwn that Wct = sup At As r(t s). (1.7)