## Closed-Form Waiting Time Approximations for Polling Systems ∗

M.A.A. Boon †

marko@win.tue.nl

E.M.M. Winands emm.winands@few.vu.nl

‡

I.J.B.F. Adan †

iadan@win.tue.nl

A.C.C. van Wijk § a.c.c.v.wijk@tue.nl

December, 2010

## Abstract

A typical polling system consists of a number of queues, attended by a single server in a fixed order. The vast majority of papers on polling systems focusses on Poisson arrivals, whereas very few results are available for general arrivals. The current study is the first one presenting simple closed-form approximations for the mean waiting times in polling systems with renewal arrival processes, performing well for all workloads. The approximations are constructed using heavy traffic limits and newly developed light traffic limits. The closed-form approximations may prove to be extremely useful for system design and optimisation in application areas as diverse as telecommunication, maintenance, manufacturing and transportation.

Keywords: Polling, waiting times, queue lengths, approximation

## 1 Introduction

Polling systems are queueing systems consisting of multiple queues, visited by a single server -typically in a fixed, cyclic order. They find their origin in many real-life applications, e.g. (computer) communication, production and manufacturing environments, traffic and transportation. For a good literature overview of polling systems and their applications, we refer to surveys of, e.g., Takagi [21], Levy and Sidi [14], and Vishnevskii and Semenova [24]. When studying literature on polling systems, it rapidly becomes apparent that the computation of the distributions and moments of the waiting times and marginal queue lengths is very cumbersome. Closed form expressions do not exist, and even when one specifies the number of queues and solves the set of equations that leads to the mean waiting times, the obtained

∗ The research was done in the framework of the BSIK/BRICKS project, and of the European Network of Excellence Euro-NF.

† Eurandom and Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands

‡ Department of Mathematics, Section Stochastics, VU University, De Boelelaan 1081a, 1081HV Amsterdam, The Netherlands

§ Department of Industrial Engineering &amp; Innovation Sciences and Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands

expressions are still too lengthy and complicated to interpret directly. Numerical procedures, both approximate and exact, have been developed in the past to compute these performance measures. However, these methods have several drawbacks. Firstly, they are not transparent and act as a kind of black box. It is, for instance, rather difficult to study the impact of parameters like the occupation rate and the service level. Secondly, these procedures are computationally complex and hard, if not impossible, to implement in a standard spreadsheet program commonly used on the work floor. Finally, the vast majority of standard methods focusses on Poisson arrival processes, which may not be very realistic in many application areas. In the present paper we study polling systems in which the arrival streams are not (necessarily) Poisson, i.e., the interarrival times follow a general distribution. The goal is to derive closed-form approximate solutions for the mean waiting times and mean marginal queue lengths, which can be computed by simple spreadsheet calculations.

Our approach in developing an approximation for the mean waiting times uses novel developments in polling literature. Recently, a heavy traffic (HT) limit has been developed for the mean waiting times as the system becomes saturated [6, 17, 23]. In the present paper we derive an approximation for the light traffic (LT) limit, i.e. as the load decreases to zero, which is exact for Poisson arrivals. The main idea is to create an interpolation between the LT limit and the HT limit. This interpolation yields good results, and has several nice properties, like satisfying the Pseudo Conservation Law (PCL), and being exact for symmetric systems with Poisson arrivals and in many limiting cases. These properties are described in more detail in the present paper. In polling literature, several alternative approximations have been developed before, most of which assume Poisson arrivals. For polling systems with Poisson arrivals and gated or exhaustive service, the best results, by far, are obtained by an approximation based on the PCL (see, e.g., [3, 7, 11]). Fischer et al. [8] study an approximation for the mean waiting times in polling systems, which is also based on an interpolation between (approximate) LT and HT limits. Their approach, however, is applied to a system with Poisson arrivals and time-limited service. Hardly any closed-form approximations exist for non-Poisson arrivals. The few that exist, perform well in specific limiting cases, e.g., under HT conditions [17, 23], or if switch-over times become very large [27, 28], but performance deteriorates rapidly if these limiting conditions are abandoned, in contrast to the approximation developed in the present paper. We show in an extensive numerical study that the quality of our approximation can be compared to the PCL approximation for systems with Poisson arrivals, but provides good results as well for systems with renewal arrivals.

Because of its simple form, the approximation function is very suitable for optimisation purposes and implementation in a spreadsheet. Although only the mean waiting times of systems with exhaustive or gated service are studied, the results can be extended to higher moments and general branching-type service disciplines. Polling systems with polling tables and/or batch service can also be analysed in a similar manner.

The structure of the present paper is as follows: the next section introduces the model and the required notation, and states the main result. Section 3 illustrates how this main result is obtained, while Section 4 provides results on the accuracy of the approximation for a large set of combinations of input parameter values. The last section discusses further research topics and possible extensions of the model.

## 2 Model description and main result

The model under consideration is a polling system consisting of N queues, Q , . . . , Q 1 N , with renewal arrival processes. Indices throughout the present paper are understood to be modulo N : Q N +1 actually refers to Q 1 . Whenever a server switches from Q i to Q i +1 , a random switch-over time S i is incurred. The generic service requirement of a customer arriving in Q i , also referred to as a type i customer, is denoted by the random variable B i . We make the usual independence assumptions for polling systems; the interarrival times, service times and switch-over times are all independent. The moment at which the server switches from one queue to the next queue, is determined by the service discipline of the queue that is being served. In the present paper we focus on polling systems in which each queue is either served according to the gated service discipline, which states that during the course of a visit of the server to Q i , only those type i customers are served that were present at the beginning of that visit, or according to the exhaustive service discipline, which means that the server keeps on serving type i customers until Q i is empty, before switching to Q i +1 .

We regard several variables as a function of the load ρ in the system. Scaling is done by keeping the service time distributions fixed, and varying the interarrival times. For each variable x that is a function of the load in the system, ρ , its value evaluated at ρ = 1 is denoted by ˆ. x For ρ = 1, the generic interarrival time of the stream in Q i is denoted by ˆ . A i Reducing the load ρ is done by scaling the interarrival times, i.e., taking the random variable A i := ˆ A /ρ i as generic interarrival time at Q i . After scaling, the load at Q i becomes ρ i = ρ E [ B i ] E [ ˆ ] A i . The (scaled) rate of the arrival stream at Q i is defined as λ i = 1 / E [ A i ]. Similarly, we define arrival rates ˆ λ i = 1 / E [ ˆ ], A i and proportional load at Q i , ˆ ρ i = ρ i ρ ('proportional' because ∑ N i =1 ˆ ρ i = 1). The system is assumed to be stable, so ρ is varied between 0 and 1. We use B to denote the generic service requirement of an arbitrary customer entering the system, with E [ B k ] = ∑ N i =1 ˆ λ i E [ B k i ] ∑ N j =1 ˆ λ j for any integer k &gt; 0, and S = ∑ N i =1 S i denotes the total switch-over time in a cycle. Finally, the (equilibrium) residual length of a random variable X is denoted by X res , with E [ X res ] = 1 E [ X / 2 ] E [ X ].

2

We now present the main result of this paper, which is a closed-form approximation formula for the mean waiting time E [ W i ] of a type i customer as a function of ρ :

$$\mathbb { E } [ W _ { i, a p p } ] = \frac { K _ { 0, i } + K _ { 1, i } \rho + K _ { 2, i } \rho ^ { 2 } } { 1 - \rho }, \quad i = 1, \dots, N.$$

The constants K 0 ,i , K 1 ,i , and K 2 ,i depend on the input parameters and the service discipline. If all queues receive exhaustive service, the constants become:

$$K _ { 0, i } = \mathbb { E } [ S ^ { \text{res} } ],$$

$$\nu, \quad \Gamma ^ { \nu } _ { 1, i } = & \hat { \rho } _ { i } ( \mathbb { E } [ \hat { A } _ { i } ] \hat { \rho } _ { i } ( 0 ) - 1 ) \mathbb { E } [ B _ { i } ^ { r e s } ] + \mathbb { E } [ B ^ { r e s } ] + \hat { \rho } _ { i } ( \mathbb { E } [ S ^ { r e s } ] - \mathbb { E } [ S ] ) - \frac { 1 } { \mathbb { E } [ S ] } \sum _ { j = 0 } ^ { N - 1 } \sum _ { k = 0 } ^ { j } \hat { \rho } _ { i + k } \mathbb { V } \tilde { \Gamma } [ S _ { i + j } ], \\ & \quad \,. \,. \,. \,. \,. \,.$$

$$K _ { 2, i } = & \frac { 1 - \hat { \rho } _ { i } } { 2 } \left ( \frac { \sum _ { j = 1 } ^ { N } \hat { \lambda } _ { j } \left ( \mathbb { V } { \text{ar} } [ B _ { j } ] + \hat { \rho } _ { j } ^ { 2 } \mathbb { V } { \text{ar} } [ \hat { A } _ { j } ] \right ) } { \sum _ { j = 1 } ^ { N } \hat { \rho } _ { j } ( 1 - \hat { \rho } _ { j } ) } + \mathbb { E } [ S ] \right ) - K _ { 0, i } - K _ { 1, i }.$$

If all queues receive gated service, we get:

$$K _ { 0, i } = \mathbb { E } [ S ^ { \text{res} } ],$$

$$\Phi _ { 1, i } = & \hat { \rho } _ { i } ( \mathbb { E } [ \hat { A } _ { i } ] \hat { \rho } _ { i } ( 0 ) - 1 ) \mathbb { E } [ B _ { i } ^ { r e s } ] + \mathbb { E } [ B ^ { r e s } ] + \hat { \rho } _ { i } \mathbb { E } [ S ^ { r e s } ] - \frac { 1 } { \mathbb { E } [ S ] } \sum _ { j = 0 } ^ { N - 1 } \sum _ { k = 0 } ^ { j } \hat { \rho } _ { i + k } \text{Var} [ S _ { i + j } ], \\ & \quad \Phi _ { 1, i } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phis _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j } \Phi _ { j + k } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psi _ { j } \Psigma _ { j } \Psi _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j + k } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psi _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j + k } \Psigma _ { j } \Psigma _ _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j} \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma _ { j } \Psigma }$$

$$K _ { 2, i } = & \frac { 1 + \hat { \rho } _ { i } } { 2 } \left ( \frac { \sum _ { j = 1 } ^ { N } \hat { \lambda } _ { j } \left ( \mathbb { V } { \text{ar} [ B _ { j } ] } + \hat { \rho } _ { j } ^ { 2 } \mathbb { V } { \text{ar} [ \hat { A } _ { j } ] } \right ) } { \sum _ { j = 1 } ^ { N } \hat { \rho } _ { j } ( 1 + \hat { \rho } _ { j } ) } + \mathbb { E } [ S ] \right ) - K _ { 0, i } - K _ { 1, i }.$$

The term ˆ ( ) is the density of g i t ˆ , the interarrival times at A i ρ = 1. This term is discussed in more detail in the next section, but for practical purposes it is useful to know that E [ ˆ ]ˆ (0) A g i i can be very well approximated by

$$\mathbb { E } [ \hat { A } _ { i } ] \hat { g } _ { i } ( 0 ) \approx \begin{cases} 2 \frac { c v _ { A _ { i } } ^ { 2 } } { c v _ { A _ { i } } ^ { 2 } + 1 } & \quad \text{ if } c v _ { A _ { i } } ^ { 2 } > 1, \\ \left ( c v _ { A _ { i } } ^ { 2 } \right ) ^ { 4 } & \quad \text{ if } c v _ { A _ { i } } ^ { 2 } \leq 1, \end{cases}$$

where cv 2 A i is the squared coefficient of variation (SCV) of A i (and, hence, also of ˆ ). A i Note that this simplification results in an approximation that requires only the first two moments of each input variable (i.e., service times, switch-over times, and interarrival times).

Remark 2.1 In case of Poisson arrivals, the constants K 1 ,i and K 2 ,i simplify considerably. E.g., for exhaustive service they simplify to:

$$K _ { 1, i } ^ { P o i s s o n } = & \mathbb { E } [ B ^ { r e s } ] + \hat { \rho } _ { i } \left ( \mathbb { E } [ S ^ { r e s } ] - \mathbb { E } [ S ] \right ) - \frac { 1 } { \mathbb { E } [ S ] } \sum _ { j = 0 } ^ { N - 1 } \sum _ { k = 0 } ^ { j } \hat { \rho } _ { i + k } \mathbb { V } \text{ar} [ S _ { i + j } ],$$

$$K _ { 2, i } ^ { P o i s s o n } = & ( 1 - \hat { \rho } _ { i } ) \left ( \frac { \mathbb { E } [ B ^ { r e s } ] } { \sum _ { j = 1 } ^ { N } \hat { \rho } _ { j } ( 1 - \hat { \rho } _ { j } ) } + \frac { \mathbb { E } [ S ] } { 2 } \right ) ^ { \cdot } - K _ { 0, i } - K _ { 1, i } ^ { P o i s s o n }.$$

$$\AA _ { i + j }, \\ \AA _ { i + j },$$

The derivation of this approximative formula for the mean waiting time is the topic of the next section. An approximation for the mean queue length at Q i , E [ L i ] is obtained by application of Little's Law to the sojourn time of type i customers, i.e. the waiting time plus the service time. As a function of ρ , we have

$$\mathbb { E } [ L _ { i, a p p } ] = \rho \frac { \mathbb { E } [ W _ { i, a p p } ] + \mathbb { E } [ B _ { i } ] } { \mathbb { E } [ \hat { A } _ { i } ] }.$$

## 3 Derivation of the approximation

Approximation (2.1) is an interpolation approximation based on LT and HT limits. In the next subsection we first provide a motivation for this approach.

## 3.1 Generic interpolation function

In its generic form, the interpolation approximation proceeds as follows; see [18, 20]. Consider an open queueing system with load ρ . Let f ( ρ ), 0 ≤ ρ &lt; 1, be some function of the queueing system (such as the mean waiting time), which is assumed to be analytic on [0 1), i.e., by , Taylor's Theorem f ( ρ ) can be expressed as

$$f ( \rho ) = \sum _ { n = 0 } ^ { \infty } \frac { f ^ { ( n ) } ( 0 ) } { n! } \rho ^ { n }, \ \ 0 \leq \rho < 1,$$

where f ( n ) ( ρ ) denotes the n th derivative f ( ρ ). Usually, f ( ρ ) is intractable, but it may be possible to derive partial information about f ( ρ ), such as the light traffic limits f ( n ) (0) for n = 0 1 , , . . . , k and the 'canonical' heavy traffic limit h = lim ρ → 1 (1 -ρ f ) ( ρ ). For examples, see [9, 18, 25], where based on the partial information, an approximation for f ( ρ ) is constructed of the form

$$\tilde { f } ( \rho ) = \frac { q ( \rho ) } { 1 - \rho },$$

where q ρ ( ) is the ( k +1)st degree polynomial, uniquely determined by the requirement that ˜ ( f ρ ) has to match everything that is known about f ( ρ ), i.e., ˜ f ( n ) (0) = f ( n ) (0) for n = 0 1 , , . . . , k and lim ρ → 1 (1 -ρ f ) ˜ ( ρ ) = h . The heavy traffic limit implies (see Prop. 1 in [20]),

$$\lim _ { n \to \infty } \frac { f ^ { ( n ) } ( 0 ) } { n! } = h.$$

This suggests that, in the Taylor series, f ( n ) (0) for n &gt; k can be approximated by n h ! . Thus, combined with knowledge of the light traffic limits f ( n ) (0) for n = 0 1 , , . . . , k , the following new approximation can be produced,

$$\bar { f } ( \rho ) = \sum _ { n = 0 } ^ { k } \frac { f ^ { ( n ) } ( 0 ) } { n! } \rho ^ { n } + h \frac { \rho ^ { k + 1 } } { 1 - \rho }.$$

Interestingly, this seemingly different approximation ¯ ( f ρ ) is identical to the interpolation approximation ˜ ( f ρ ), confirming the notion that they are the 'natural' approximation for f ( ρ ), given the partial information. In [18] the interpolation approximation (3.1) is shown to work extremely well for several examples. The present paper applies approximation (3.1) to the new setting of polling systems with general renewal arrivals, for which no analytic expressions are known for the mean waiting times. Choosing f ( ρ ) as the mean waiting time of a type i customer, we derive new (approximations for the) light traffic limits f (0) and the first derivative f ′ (0), which together with the heavy traffic limit, yield an interpolation with a quadratic polynomial q ρ ( ); see (2.1).

In Sections 3.2 and 3.3 we derive the LT and HT limits, respectively. The interpolation approximation, matching these limits, is presented in Section 3.4. Finally, in Section 3.5, we show that the interpolation approximation also matches known exact results for mean waiting times in polling systems. In particular, we prove that the form (3.1) is crucial for satisfying the pseudo-conservation law for all loads ρ . An extensive numerical validation is the topic of Section 4, showing that the interpolation approximation works extremely well for polling systems.

## 3.2 Light traffic

The mean waiting times in the polling model under consideration in light-traffic, have been studied in Blanc and Van der Mei [2], under the assumption of Poisson arrivals. They obtain expressions for the mean waiting times in light traffic that are exact up to (and including) firstorder terms in ρ . These expressions have been found by carefully inspecting numerical results obtained with the Power-Series Algorithm, but no proof is provided. In the present section we shall not only prove the correctness of the light-traffic results in a system with Poisson arrivals, but also use them as base for an approximation for the mean waiting times in polling systems with renewal interarrival times. The key ingredient to the LT analysis of a polling system, is the well-known Fuhrmann-Cooper decomposition [10]. It states that in a vacation system with Poisson arrivals the queue length of a customer is the sum of two independent random variables: the number of customers in an isolated M/G/ 1 queue, and the number of customers during an arbitrary moment in the vacation period. The distributional form of Little's Law [12] can be used to translate this result to waiting times. Since no independence is required between the length of a vacation and the length of the preceding visit period, this decomposition also holds for polling systems with Poisson arrivals. We introduce V i to denote the length of a visit period to Q i , and I i to denote the length of the intervisit period, i.e. the time that the server is away between two successive visits to Q i . Using C i to denote the cycle time, starting at a visit beginning to Q i , we have E [ V i ] = ρ i E [ C i ] and E [ I i ] = (1 -ρ i ) E [ C i ]. It is well-known that the mean cycle time in polling systems, unlike higher moments, does not depend on the starting point: E [ C i ] = E [ C ] = E [ S ] 1 -ρ .

The Fuhrmann-Cooper decomposition, applied to the mean waiting time, results in:

$$\text{exhaustive} \colon & & \mathbb { E } [ W _ { i } ] = \mathbb { E } [ W _ { i, M / G / 1 } ] + \mathbb { E } [ I _ { i } ^ { r e s } ], & & ( 3. 2 )$$

$$\mathbb { E } [ W _ { i } ] = \mathbb { E } [ W _ { i, M / G / 1 } ] + \mathbb { E } [ I _ { i } ^ { r e s } ],$$

gated:

$$\mathbb { E } [ W _ { i } ] = \mathbb { E } [ W _ { i, M / G / 1 } ] + \mathbb { E } [ I _ { i } ^ { r e s } ] + \frac { \mathbb { E } [ V _ { i } I _ { i } ] } { \mathbb { E } [ I _ { i } ] }. \quad \quad ( 3. 3 )$$

For our approximation, we assume that this decomposition also holds for renewal arrival processes in light traffic. Determining the LT limit of the mean waiting time, E [ W LT i ], in a polling system with exhaustive or gated service is based on the following two-step approach. The first step is to find the LT limit of E [ W i,GI/G/ 1 ], the mean waiting time of a GI/G/ 1 queue with only type i customers in isolation, i = 1 , . . . , N . The second step is determining E [ I res i ], the mean residual intervisit time of Q i , and E [ V I i i ] E [ I i ] , the mean visit time of Q i given that it is being observed at a random epoch during the following intervisit time.

Remark 3.1 Bertsimas and Mourtzinou [1] state that the decompositions (3.2) and (3.3) also hold for polling systems with Mixed Generalised Erlang arrivals. However, simulation and exact analysis of some simple cases indicate that the decomposition result is not valid for the mean waiting times.

For the LT limit of the mean waiting time in a GI/G/ 1 queue, we use Whitt's result (Equation (16) in [25]), which gives:

$$\lim _ { \rho _ { i } \downarrow 0 } \frac { \mathbb { E } [ W _ { i, G I / G / 1 } ] } { \rho _ { i } } = \frac { 1 + c \upsilon _ { B _ { i } } ^ { 2 } } { 2 } \mathbb { E } [ \hat { A } _ { i } ] \hat { g } _ { i } ( 0 ) \mathbb { E } [ B _ { i } ],$$

where cv 2 B i is the SCV of the service times, and ˆ ( ) is the density of the interarrival times g i t ˆ . A i For practical purposes, it may be more convenient to express ˆ (0) in terms of the density of g i

A i , the generic interarrival time of Q i in the scaled situation.The relation between the density of the scaled interarrival times A i (= ˆ A /ρ i ), denoted by g i ( ), and the density of t ˆ , A i ˆ ( ), is g i t simply: g i ( ) = t ρg ˆ ( i ρt ). This means that the term E [ ˆ ]ˆ (0) can be rewritten as A g i i

$$\mathbb { E } [ \hat { A } _ { i } ] \hat { g } _ { i } ( 0 ) = \mathbb { E } [ A _ { i } ] g _ { i } ( 0 ).$$

Because of this equality, in the remainder of the paper we might use either notation. Since determining E [ ˆ ]ˆ (0) is a required step in the computation of our approximation for A g i i E [ W i ], we give some practical examples.

Example 1 If the scaled interarrival times A i are exponentially distributed with parameter λ i := 1 / E [ A i ], we have g i ( ) = t λ i e -λ t i . This implies that E [ A g i ] i (0) = 1.

Example 2 In this example we assume that A i follows a H 2 distribution with balanced means. The SCV of A i is denoted by cv 2 A i . The density of this hyper-exponential distribution is (see, e.g., [22])

$$g _ { i } ( t ) = p \mu _ { 1 } \text{e} ^ { - \mu _ { 1 } t } + ( 1 - p ) \mu _ { 2 } \text{e} ^ { - \mu _ { 2 } t },$$

with

$$p \ & = \ \frac { 1 } { 2 } \left ( 1 + \sqrt { \frac { c v _ { A _ { i } } ^ { 2 } - 1 } { c v _ { A _ { i } } ^ { 2 } + 1 } } \right ), \\ \mu _ { 1 } \ & = \ \frac { 1 } { \mathbb { E } \left [ A _ { i } \right ] } \left ( 1 + \sqrt { \frac { c v _ { A _ { i } } ^ { 2 } - 1 } { c v _ { A _ { i } } ^ { 2 } + 1 } } \right ), \\ \mu _ { 2 } \ & = \ \frac { 1 } { \mathbb { E } \left [ A _ { i } \right ] } \left ( 1 - \sqrt { \frac { c v _ { A _ { i } } ^ { 2 } - 1 } { c v _ { A _ { i } } ^ { 2 } + 1 } } \right ). \\ \dots ^ { 2 } \ \, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

This leads to E [ A g i ] i (0) = 1 + cv 2 A -1 cv 2 A +1 = 2 cv 2 A cv 2 A +1 .

Example 3 Now we assume that the interarrival times follow a mixed Erlang distribution. The density of the scaled interarrival times is:

$$g _ { i } ( t ) = p \frac { \mu ^ { k - 1 } t ^ { k - 2 } } { ( k - 2 )! } \text{e} ^ { - \mu t } + ( 1 - p ) \frac { \mu ^ { k } t ^ { k - 1 } } { ( k - 1 )! } \text{e} ^ { - \mu t },$$

i.e., a mixture of an Erlang( k -1) and an Erlang( k ) distribution with

$$\Pi \text{Elam} ( \varkappa - 1 ) \text{ and } \Pi \text{Elam} ( \varkappa ) \text{ using the output} \, \text{ with} \, \\ k \quad = \quad \left \lceil \frac { 1 } { c v _ { A _ { i } } ^ { 2 } } \right \rceil, \\ p \quad = \quad \frac { k \, c v _ { A _ { i } } ^ { 2 } - \sqrt { k ( 1 + c v _ { A _ { i } } ^ { 2 } ) - k ^ { 2 } \, c v _ { A _ { i } } ^ { 2 } } } { 1 + c v _ { A _ { i } } ^ { 2 } }, \\ \mu \quad = \quad \frac { k - p } { \mathbb { E } [ A _ { i } ] }. \\ \text{to $\mathbb{F} [ A ] } a ( 0 ) = 0$$

If k &gt; 2, this leads to E [ A i ] g i (0) = 0.

The distributions in Examples 1 -3 are typical distributions to be used in a two-moment fit if the SCV of the interarrival times is respectively 1, greater than 1, and less than 1 (cf. [22]). The examples illustrate how E [ A g i ] i (0) can be computed if the density of the (scaled) interarrival times is known. If no information is available about the complete density, but the first two moments of A i are known, Whitt suggests to use the following approximation for E [ A g i ] i (0):

$$\mathbb { E } [ A _ { i } ] g _ { i } ( 0 ) = \begin{cases} 2 \frac { c v _ { A _ { i } } ^ { 2 } } { c v _ { A _ { i } } ^ { 2 } + 1 } & \quad \text{ if } c v _ { A _ { i } } ^ { 2 } > 1, \\ \left ( c v _ { A _ { i } } ^ { 2 } \right ) ^ { 4 } & \quad \text{ if } c v _ { A _ { i } } ^ { 2 } \leq 1, \end{cases}$$

where cv 2 A i is the squared coefficient of variation of the interarrival times of Q i . This approximation is exact for cv 2 A i &gt; 1, if the interarrival time distribution is a hyper-exponential distribution as discussed in Example 2. For cv 2 A i ≤ 1, the approximation is rather arbitrary, but Example 3 shows that E [ A g i ] i (0) becomes small (or even zero) very rapidly as cv 2 A i gets smaller.

Summarising, the LT limit of a GI/G/ 1 queue (ignoring O ( ρ 2 i ) terms and higher) is:

$$\mathbb { E } [ W _ { i, G I / G / 1 } ^ { L T } ] = \rho _ { i } \mathbb { E } [ A _ { i } ] g _ { i } ( 0 ) \mathbb { E } [ B _ { i } ^ { r e s } ].$$

For Poisson arrivals ( E [ A g i ] i (0) = 1), it is known that E [ W i,M/G/ 1 ] = ρ i 1 -ρ i E [ B res i ] = ρ i E [ B res i ]+ O ( ρ 2 i ), which is consistent with our approximation.

The second step in determining the LT limit of the mean waiting time of a type i customer in a polling system, is finding the LT limits of E [ I res i ], the mean residual intervisit time of Q i , and (for gated service only) E [ V I i i ] E [ I i ] , the mean visit time V i given that it is observed from the following intervisit time I i . In this LT analysis we need to focus on first order terms only. Noting the fact that I i = S i + V i +1 + S i +1 + · · · + V i + N -1 + S i + N -1 , we condition on the moment at which I i is observed. We distinguish between two cases. The moment of observation either takes place during a visit time, or during a switch-over time:

$$\Phi _ { \Phi } \Phi _ { i } ^ { \Phi } & = \sum _ { j = 1 } ^ { N - 1 } \frac { \mathbb { E } [ V _ { i + j } ] } { \mathbb { E } [ I _ { i } ] } \mathbb { E } [ I _ { i } ^ { L T, r e s } | \text{observed during } V _ { i + j } ] \\ & \quad + \sum _ { j = 0 } ^ { N - 1 } \frac { \mathbb { E } [ S _ { i + j } ] } { \mathbb { E } [ I _ { i } ] } \mathbb { E } [ I _ { i } ^ { L T, r e s } | \text{observed during } S _ { i + j } ].$$

glyph[negationslash]

Observation during visit time. The probability that a random observation epoch takes place during a visit time, say V j , is E [ V j ] E [ I i ] , for any j = . However, we are only interested in i order ρ terms, so this probability simplifies to

$$\frac { \mathbb { E } [ V _ { j } ] } { \mathbb { E } [ I _ { i } ] } = \frac { \rho _ { j } \mathbb { E } [ C ] } { ( 1 - \rho _ { i } ) \mathbb { E } [ C ] } = \rho _ { j } + \mathcal { O } ( \rho ^ { 2 } ).$$

The fact that this probability is O ( ρ ), implies that all further O ( ρ ) terms can be ignored in E [ I LT,res i | observed during V j ], because in LT we focus on first order terms only.

The length of the residual intervisit time is the length of the residual visit period of type j customers, V res j , plus all switch-over times S j + · · · + S i -1 , plus all visit times V j +1 + · · · + V i -1 .

The first term simplifies to E [ V res j ] = E [ B res j ] + O ( ρ ). The terms E [ V k | observed from V j ] , k = j + 1 , . . . , i -1, in light traffic, are all O ( ρ ). Summarising, the mean residual intervisit period when observed during V j is simply a mean residual service time E [ B res j ], plus all mean switch-over times E [ S j + · · · + S i -1 ], plus O ( ρ ) terms:

$$\mathbb { E } [ I _ { i } ^ { L T, r e s } | \text{observed during } V _ { j } ] = \mathbb { E } [ B _ { j } ^ { r e s } ] + \sum _ { k = j } ^ { i - 1 } \mathbb { E } [ S _ { k } ] + \mathcal { O } ( \rho ).$$

The intuition behind this equation is that, in light traffic, the probability of having another service during the residual cycle is negligible, i.e., O ( ρ ). Hence, the length of the residual intervisit time is solely determined by the residual service time and the remaining switch-over times in the cycle.

Observation during switch-over time. We continue by determining the mean residual intervisit period, conditioned on a random observation epoch during a switch-over time, say S j , j = 1 , . . . , N . The probability that such an epoch takes place during S j , is

$$\frac { \mathbb { E } [ S _ { j } ] } { \mathbb { E } [ I _ { i } ] } = \frac { \mathbb { E } [ S _ { j } ] } { ( 1 - \rho _ { i } ) \mathbb { E } [ C ] } = \frac { \mathbb { E } [ S _ { j } ] } { \mathbb { E } [ S ] } \frac { 1 - \rho } { 1 - \rho _ { i } } = \frac { \mathbb { E } [ S _ { j } ] } { \mathbb { E } [ S ] } ( 1 - \rho + \rho _ { i } ) + \mathcal { O } ( \rho ^ { 2 } ).$$

It becomes apparent from this expression that things get slightly more complicated now, because order ρ terms in the conditional residual intervisit time may no longer be neglected. The residual intervisit time now consists of the residual switch-over time S res j , plus the switchover times S j +1 + · · · + S i -1 , plus all visit periods V j +1 + · · · + V i -1 . The length of a visit period V k , for k &gt; j , is the sum of the busy periods of all type k customers that have arrived during S k -N , . . . , S j -1 , S past j , S res j , and S j +1 , . . . , S k -1 . By S past j we denote the elapsed switch-over time during which the intervisit period is observed, which has the same distribution as the residual switch-over time S res j . Compared to an observation during a visit time, it is more difficult to determine the conditional mean length of a busy period E [ V k | observed during S j ] under LT. We use a heuristic approach, which is exact if the arrival process of type k customers is Poisson, and approximate it by:

glyph[negationslash]

$$\mathbb { E } [ V _ { k } | \text{observed during } S _ { j } ] \approx \rho _ { k } \left ( \sum _ { l \neq j } \mathbb { E } [ S _ { l } ] + \mathbb { E } [ S _ { j } ^ { p a s t } ] + \mathbb { E } [ S _ { j } ^ { r e s } ] \right ) + \mathcal { O } ( \rho ^ { 2 } ), \quad k = j + 1, \dots, i - 1.$$

If A k is exponentially distributed, the above expression is exact. Nevertheless, numerical experiments have shown that this approximative assumption has no or at least negligible impact on the accuracy of the approximated mean waiting times. Summarising:

$$\ell [ I _ { i } ^ { L T, r e s } | \text{observed during } S _ { j } ] & \approx \sum _ { k = i } ^ { j - 1 } \mathbb { E } [ S _ { k } ] ( \sum _ { l = j + 1 } ^ { i + N - 1 } \rho _ { l } ) + \mathbb { E } \left ( S _ { j } ^ { p a s t } \right ) ( \sum _ { k = j + 1 } ^ { i + N - 1 } \rho _ { k } ) \\ & \quad + \mathbb { E } \left ( S _ { j } ^ { r e s } \right ) ( 1 + \sum _ { k = j + 1 } ^ { i + N - 1 } \rho _ { k } ) + \sum _ { k = j + 1 } ^ { i + N - 1 } \mathbb { E } [ S _ { k } ] ( 1 + \sum _ { l = j + 1 } ^ { i + N - 1 } \rho _ { l } ) + \mathcal { O } ( \rho ^ { 2 } ).$$

$$( 3. 8 )$$

The expression for I res i under light traffic conditions now follows from substituting (3.7) and (3.8) in (3.6). The result can be rewritten to:

$$& \text{The expression for } I _ { i } ^ { r e s } \text{ under light traffic conditions now follows from substituting (3.7) and } ( 3. 8 ) \text{ in (3.6).  The result can be rewritten to: } \\ & \text{[T](I_T, r e s ] } \approx \sum _ { j = i + 1 } ^ { i + N - 1 } \rho _ { j } \text{E[B]} ^ { r e s ] } + \sum _ { j = i + 1 } ^ { i + N - 1 } \rho _ { j } \sum _ { k = j } ^ { \text{E[S,k]} } \\ & \quad \quad \quad + \sum _ { j = i } ^ { i + N - 1 } \frac { 1 } { 2 \text{E[S] } } \left [ \begin{array} { c } \text{$E(S_{j})$ (1 - \rho + \rho_{i}+2 \sum _ { k=j+1})$ } \rho _ { k}) \end{array} \right ] \\ & \quad \quad + \sum _ { j = i + 1 } ^ { i + N - 1 } \left [ \begin{array} { c } j - 1 \\ j - 1 \\ \sum _ { j = i } ^ { \text{E[S,j]}$\text{E[S,k]}$( \sum _ { j = i + 1 } ^ { i + N - 1 } \frac { i + N - 1 } { 2 \text{E[S,j]}$\text{E[S,k]}$(1 - \rho + \rho_{i}+2 \sum _ { l=j+1})$ } \rho _ { l=j+1 } \end{array} \right ] \\ & \quad \quad + \overline { \text{E[S] } } \sum _ { j = i + 1 } ^ { i + N - 1 } \left [ \begin{array} { c } j - 1 \\ j - 1 \\ \sum _ { k = j + 1 } ^ { \text{E[S,j]}$\text{E[S,k]}$( \sum _ { j = i + 1 } ^ { i + N - 1 } \frac { i + N - 1 } { 2 \text{E[S,j]}$\text{E[S,k]}$(1 - \rho + \rho_{i}+2 \sum _ { l=j+1})$ } \rho _ { l=j+1 } \end{array} \right ] \\ & \quad \quad + \overline { \text{O(r}^{2)}$} \\ & \quad \quad \quad$$

for i = 1 , . . . , N . The last step in (3.9) follows after some straightforward (but tedious) rewriting.

The Fuhrmann-Cooper decomposition of the mean waiting time for customers in a polling system with gated service (3.3), also requires the computation of E [ V I i i ] E [ I i ] under LT conditions.

Here, again, we have to resort to using a heuristic and use E [ V I i i ] E [ I i ] = ρ i E [ S ] + O ( ρ 2 ), because this value is exact in the case of Poisson arrivals. This term is the mean length of the visit time V i given that it is observed during the following intervisit time I i . The term appears because, contrary to exhaustive service, type i customers arriving during V i are not served until the next cycle. However, it is easier to consider E [ V I i i ] / E [ V i ] instead, and to use the relation

$$\frac { \mathbb { E } [ V _ { i } I _ { i } ] } { \mathbb { E } [ I _ { i } ] } = \frac { \mathbb { E } [ V _ { i } I _ { i } ] } { \mathbb { E } [ V _ { i } ] } \times \frac { \mathbb { E } [ V _ { i } ] } { \mathbb { E } [ I _ { i } ] }.$$

The term E [ V I i i ] / E [ I i ] is the mean length of the intervisit time I i following V i , given that it is observed during this visit time V i . Firstly, we note that

$$\frac { \mathbb { E } [ V _ { i } ] } { \mathbb { E } [ I _ { i } ] } = \frac { \rho _ { i } } { 1 - \rho _ { i } } = \rho _ { i } + \mathcal { O } ( \rho ^ { 2 } ).$$

This implies that we can ignore all O ( ρ ) terms in E [ V I i i ] / E [ V i ], which means that only the switch-over times play a role,

$$\frac { \mathbb { E } [ V _ { i } I _ { i } ] } { \mathbb { E } [ V _ { i } ] } = \mathbb { E } [ S ] + \mathcal { O } ( \rho ).$$

Concluding, we have

$$\frac { \mathbb { E } [ V _ { i } I _ { i } ] } { \mathbb { E } [ I _ { i } ] } = \rho _ { i } \mathbb { E } [ S ] + \mathcal { O } ( \rho ^ { 2 } ),$$

in the case of Poisson arrivals. If the arrival process is not Poisson, this is not exact, but we use it as an approximation.

Having made all required preparations, we are ready to formulate the main result of the present subsection. Under light traffic, an approximation for the mean waiting time of a type i customer in a polling model with general arrivals and respectively exhaustive and gated service in Q i , is:

$$\mathbb { E } [ W _ { i } ^ { L T, e x h } ] & \approx \mathbb { E } [ S ^ { r e s } ] + \rho _ { i } ( \mathbb { E } [ \hat { A } _ { i } ] \hat { g } _ { i } ( 0 ) - 1 ) \mathbb { E } [ B _ { i } ^ { r e s } ] + \rho \mathbb { E } [ B ^ { r e s } ] + ( \rho - \rho _ { i } ) \left ( \mathbb { E } [ S ] - \mathbb { E } [ S ^ { r e s } ] \right ) \\ & + \frac { 1 } { \mathbb { E } [ S ] } \sum _ { k = i + 1 } ^ { i + N - 1 } \rho _ { k } \sum _ { j = i } ^ { k - 1 } \mathbb { V } \text{ar} [ S _ { j } ] + \mathcal { O } ( \rho ^ { 2 } ), \quad i = 1, \dots, N, \\ \mathbb { w } \colon \colon L T. \text{atched} \quad \mathbb { w } \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \col. \end{matrix}$$

$$\mathbb { E } [ W _ { i } ^ { L T, g a t e d _ { i } } ] \approx \mathbb { E } [ W _ { i } ^ { L T, e x h _ { i } } ] + \rho _ { i } \mathbb { E } [ S ],$$

where ˆ ( ) is the density of the interarrival times of type g i t i customers at ρ = 1. Equation (3.11) follows from substitution of (3.5) and (3.9) in

$$\mathbb { E } [ W _ { i } ] \approx \mathbb { E } [ W _ { i, G I / G / 1 } ] + \mathbb { E } [ I _ { i } ^ { r e s } ], \quad i = 1, \dots, N.$$

For Poisson arrivals, (3.11) and (3.12) are exact. The LT limit for polling systems with Bernoulli service (and Poisson arrivals) has been experimentally found in [2] and, indeed, it can be shown that their result for exhaustive service, which is a special case of Bernoulli service, agrees with our result after substituting E [ ˆ ]ˆ (0) = 1 in (3.11). A g i i

## 3.3 Heavy traffic

Heavy traffic limits in polling systems have been studied by Coffman et al. [5, 6], and by Olsen and Van der Mei [16, 17]. In these papers, the HT limits of the waiting time distributions are found under the assumption of Poisson arrivals. For general renewal arrivals, a proof is given for the special case N = 2 (cf. [5, 6]), and a strong conjecture for larger values of N (cf. [17]). In [23], the following result for the mean waiting time is proven rigorously for polling systems with renewal arrivals:

$$\mathbb { E } [ W _ { i } ^ { H T } ] = \frac { \omega _ { i } } { 1 - \rho } + o ( ( 1 - \rho ) ^ { - 1 } ), \quad \rho \uparrow 1.$$

Obviously, in HT, all queues become unstable and, thus, E [ W i ] tends to infinity for all i . The rate at which E [ W i ] tends to infinity as ρ ↑ 1 is indicated by ω i , which is referred to as the mean asymptotic scaled delay at queue i , and depends on the service discipline. For exhaustive service,

$$\omega _ { i } = \frac { 1 - \hat { \rho } _ { i } } { 2 } \left ( \frac { \sigma ^ { 2 } } { \sum _ { j = 1 } ^ { N } \hat { \rho } _ { j } ( 1 - \hat { \rho } _ { j } ) } + \mathbb { E } [ S ] \right ), \quad i = 1, \dots, N,$$

with

$$\sigma ^ { 2 } \coloneqq \sum _ { i = 1 } ^ { N } \hat { \lambda } _ { i } \left ( \mathbb { V } \arho [ B _ { i } ] + \hat { \rho } _ { i } ^ { 2 } \mathbb { V } \arho [ \hat { A } _ { i } ] \right ).$$

Here, the limits are taken such that the arrival rates are increased, while keeping the servicetime distributions fixed, and keeping the distributions of the interarrival times A i ( i = 1 , . . . , N ) fixed up to a common scaling constant ρ . Notice that in the case of Poisson arrivals we have σ 2 = E [ B 2 ] / E [ B ].

For gated service, we have

$$\omega _ { i } = \frac { 1 + \hat { \rho } _ { i } } { 2 } \left ( \frac { \sigma ^ { 2 } } { \sum _ { j = 1 } ^ { N } \hat { \rho } _ { j } ( 1 + \hat { \rho } _ { j } ) } + \mathbb { E } [ S ] \right ).$$

## 3.4 Interpolation

Now that we have the expressions for the mean delay in both LT and HT, we can determine the constants K 0 ,i , K 1 ,i , and K 2 ,i in approximation formula (2.1). We simply impose the requirements that approximation (2.1) results in the same mean waiting time for ρ = 0 as the LT limit, and for ρ ↑ 1 as the HT limit. Since (3.11) (and (3.12) for gated service) has been determined up to the first order of ρ terms, we also add the requirement that the derivative with respect to ρ , taken at ρ = 0, of our approximation is equal to the derivative of the LT limit. A more formal definition of these requirements is presented below:

$$\mathbb { E } [ W _ { i, a p p } ] \Big | _ { \rho = 0 } & = \mathbb { E } [ W _ { i } ] \Big | _ { \rho = 0 }, \\ \frac { d } { d \rho } \mathbb { E } [ W _ { i, a p p } ] \Big | _ { \rho = 0 } & = \frac { d } { d \rho } \mathbb { E } [ W _ { i } ] \Big | _ { \rho = 0 }, \\ ( 1 - \rho ) \mathbb { E } [ W _ { i, a p p } ] \Big | _ { \rho = 1 } & = ( 1 - \rho ) \mathbb { E } [ W _ { i } ] \Big | _ { \rho = 1 }.$$

This leads to (2.1) as approximation for E [ W i ] in a polling system with general arrivals. Constants K 0 ,i , K 1 ,i , and K 2 ,i are defined in (2.2)-(2.4) for systems with exhaustive service, or (2.5)-(2.7) for gated service.

## 3.5 Matching properties

A desirable property of an approximation is that it matches known exact results. In the present section we discuss several cases where the interpolation approximation yields exact results. Most cases require Poisson arrivals, but it is shown that also in two limiting cases where exact results are available for general arrivals, i.e., heavy traffic and large switch-over times, the approximation is exact. It further turns out that, in case of Poisson arrivals, the approximated mean waiting times satisfy the pseudo-conservation law , implying that the weighted sum ∑ N i =1 ρ i E [ W i, app ] is exact for each load 0 ≤ ρ &lt; 1. In fact, this appears to be true for any interpolation approximation for the mean waiting times of the form (3.1) matching the HT limit and the LT limits of order 0 up to order k , provided k &gt; 0. These properties of the interpolation approximation indicate that it is the 'natural' approximation, given the HT and LT limits.

Light and heavy traffic. The light traffic limit of E [ W i ], given by (3.11) for exhaustive service and by (3.12) for gated service, is exact for Poisson arrivals. The heavy traffic limit (3.14) of E [ W i ] is even exact for renewal arrivals. An appropriate choice of constants K 0 ,i , K 1 ,i

and K 2 ,i can reduce (2.1) to either (3.11), (3.12), or (3.14). Since the LT and HT limits have been used in the set of equations that determine the coefficients of the approximation, it goes without saying that E [ W i, app ] is equal to (3.11) (or (3.12) for gated service) and (3.14), for ρ ↓ 0 and ρ ↑ 1 respectively. This implies that the LT limit of our approximation is exact for Poisson arrivals, and the HT limit is exact for general arrivals.

Symmetric system. If ˆ ρ i = 1 N for all i = 1 , . . . , N , all B i have the same distribution, and the variances V ar[ S i ] of all switch-over times are equal, then our approximation is exact if all interarrival distributions are exponential. For exhaustive service, we obtain

$$K _ { 1, i } = & \mathbb { E } [ B ^ { r e s } ] + \frac { N - 1 } { N } \mathbb { E } [ S ] - \left ( 2 - \frac { 1 } { N } \right ) \mathbb { E } [ S ^ { r e s } ] + \frac { 1 } { \mathbb { E } [ S ] } \sum _ { k = i + 1 } ^ { i + N - 1 } \hat { \rho } _ { k } \sum _ { j = i } ^ { k - 1 } \mathbb { V } _ { \ } [ S _ { j } ] \\ = & \mathbb { E } [ B ^ { r e s } ] + \frac { N - 1 } { N } \mathbb { E } [ S ] - \left ( 2 - \frac { 1 } { N } \right ) \mathbb { E } [ S ^ { r e s } ] + \frac { N - 1 } { N } \frac { \text{Var} [ S ] } { 2 \mathbb { E } [ S ] } \\ = & \mathbb { E } [ B ^ { r e s } ] + \left ( 1 - \frac { 1 } { N } \right ) \frac { \mathbb { E } [ S ] } { 2 } - \mathbb { E } [ S ^ { r e s } ], \\ \dots \pi \nu \quad \partial \, \text{Tr.\colon \dots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \ddots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \vdots \pi \cup \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \, \delta \}$$

and K 2 ,i = 0. This implies that E [ W i, app ] = E [ W i, symm ], since

$$\mathbb { E } [ W _ { i, s y m m } ] = \frac { \rho } { 1 - \rho } \mathbb { E } [ B ^ { r e s } ] + \mathbb { E } [ S ^ { r e s } ] + \frac { \rho ( 1 - \frac { 1 } { N } ) } { 1 - \rho } \frac { \mathbb { E } [ S ] } { 2 }.$$

Note that E [ W i, symm ] is of the form (3.1) with q ρ ( ) being a linear polynomial. In fact, this immediately implies that the interpolation is exact, given the (exact) LT and HT limits. Further, we do not require that the mean switch-over times E [ S i ] are equal. One can verify that the same holds for gated service.

Single queue (vacation model). An immediate consequence of the fact that our approximation is exact in symmetric polling systems with Poisson arrivals, is that it also gives exact results for the mean waiting time of customers in a single-queue polling system with Poisson arrivals. A polling system consisting of only one queue, but with a switch-over time between successive visits to this queue, is generally referred to as a queueing system with multiple server vacations.

Large switch-over times. For S deterministic, S →∞ , and, again, under the assumption of Poisson arrivals, it is proven in [27, 28] that E [ W i ] S → 1 -ρ i 2(1 -ρ ) for exhaustive service. It can easily be verified that our approximation has the same limiting behaviour:

$$\lim _ { S \to \infty } \frac { \mathbb { E } [ W _ { i, a p p } ] } { S } = \frac { 1 - \rho _ { i } } { 2 ( 1 - \rho ) }.$$

For gated service, E [ W i, app ] S → 1+ ρ i 2(1 -ρ ) , which is also the exact limit (see, e.g., [26]).

Olsen [15] studies the effects of large switch-over times in polling systems operating under HT conditions. She discovers that, under these conditions, the limiting behaviour (3.15) is also exhibited in polling systems with general renewal arrivals .

Miscellaneous other exact results. The approximation is also exact in several other cases, all with Poisson arrivals, when the parameter values are carefully chosen. The relations between the input parameters that yield exact approximation results become very complicated, especially in polling systems with more than two queues. We only mention one interesting example here: our approximation gives exact results for a two-queue polling system with exhaustive service and

$$\mathbb { E } [ B _ { 1 } ] = \mathbb { E } [ B _ { 2 } ], \mathbb { E } [ S _ { 1 } ] = \mathbb { E } [ S _ { 2 } ], c \nu _ { A _ { 1 } } ^ { 2 } = c \nu _ { A _ { 2 } } ^ { 2 }, c \nu _ { B _ { 1 } } ^ { 2 } = c \nu _ { B _ { 2 } } ^ { 2 }, c \nu _ { S _ { 1 } } ^ { 2 } = c \nu _ { S _ { 2 } } ^ { 2 },$$

if the following constraint is satisfied:

$$\rho = \frac { 1 + I _ { A _ { i } } ^ { 2 } } { 2 I _ { A _ { i } } } - \frac { c v _ { S _ { i } } ^ { 2 } } { 1 + c v _ { B _ { i } } ^ { 2 } } \cdot \frac { \mathbb { E } [ S _ { i } ] } { \mathbb { E } [ B _ { i } ] },$$

where I A i = ˆ ρ 1 ˆ ρ 2 is the ratio of the loads of the two queues. Obviously, if I A i = 1, the system is symmetric and our approximation gives exact results regardless of the other parameter settings.

Pseudo-conservation law. Awell-known result in polling literature, is the so-called pseudoconservation law , derived by Boxma and Groenendijk [4] using the concept of work decomposition. This law gives the following exact expression for the weighted sum of the mean waiting times in a polling system with Poisson arrivals:

$$\sum _ { i = 1 } ^ { N } \rho _ { i } \mathbb { E } [ W _ { i } ] = \frac { \rho ^ { 2 } } { 1 - \rho } \mathbb { E } [ B ^ { r e s } ] + \rho \mathbb { E } [ S ^ { r e s } ] + \frac { \mathbb { E } [ S } { 2 } \frac { \rho ^ { 2 } - \sum _ { j = 1 } ^ { N } \rho _ { j } ^ { 2 } } { 1 - \rho } + \sum _ { j = 1 } ^ { N } \mathbb { E } [ Z _ { j j } ], \quad \ ( 3. 1 8 )$$

where E [ Z jj ] denotes the mean amount of work left behind in Q j at the completion of a visit of the server to Q j . It is shown in [4] that E [ Z jj ] is the only term that depends on the service discipline. For exhaustive service E [ Z jj ] = 0, for gated service E [ Z jj ] = ρ 2 j E [ S ] 1 -ρ . It can be shown that the interpolation approximation (2.1) satisfies the pseudo-conservation law in the case of Poisson arrivals: if E [ ˆ ]ˆ (0) = 1 for A g i i i = 1 , . . . , N , then ∑ N i =1 ρ i E [ W i, app ] equals the right-hand side of (3.18). In fact, in case of Poisson arrivals, any interpolation approximation for the mean waiting times of the form (3.1), matching the HT limit and the LT limits of order 0 up to order k , satisfies (3.18), provided k &gt; 0. To establish this result, we first rewrite (3.18) in the form Q ρ ( ) = 0, by moving the terms at right-hand side of (3.18) to the left. Now let ˜ ( Q ρ ) denote the version of Q ρ ( ) where the mean waiting times have been replaced by their interpolation approximation. Clearly, the interpolation approximation and the right-hand side of (3.18) are both of the form (3.1), and thus ˜ ( Q ρ ) is also of the form

$$\tilde { Q } ( \rho ) = \frac { \tilde { q } ( \rho ) } { 1 - \rho },$$

where ˜( q ρ ) is a polynomial of degree k +2. Since the interpolation approximation matches the LT limits of order 0 up to order k , it follows that ˜ Q ( n ) (0) = Q ( n ) (0) = 0, and hence, ˜ q ( n ) (0) = 0 for all n = 0 1 , , . . . , k +1. Further, by the HT limit, ˜(1) = lim q ρ → 1 (1 -ρ Q ρ ) ˜ ( ) = lim ρ → 1 (1 -ρ Q ρ ) ( ) = 0. This implies that ˜( q ρ ) = 0, and thus ˜ ( Q ρ ) = 0 for all ρ .

## 4 Numerical study

## 4.1 Initial glance at the approximation

Before we study the accuracy of the approximation to a large test bed of polling systems, we just pick a rather arbitrary, simple system to compare the approximation with exact results in order to get some initial insights. Consider a three-queue polling system with loads of Q 1 , Q 2 , and Q 3 divided as follows: ˆ ρ 1 = 0 1 . , ρ ˆ 2 = 0 3, and . ˆ ρ 3 = 0 6. . All service times and switch-over times are exponentially distributed, with mean 1. The interarrival times have SCV cv 2 A i = 3 for i = 1 2 3. In Figure 1 we plot the approximated mean waiting time of , , Q 2 , E [ W 2 , app ], versus the load of the system ρ . Since this system cannot be analysed analytically, we compare the approximated values with simulated values. Both in the approximation and in the simulation we fit a H 2 distribution as described in Example 2.

The errors are largest for Q 2 , which is the reason why we chose this queue in particular in Figure 1. The most important information that this figure reveals, is that even though the accuracy of the approximation is worst for this queue (a relative error of -4 47% for . ρ = 0 7), . the shape of the approximation function is very close to the shape of the exact function, which makes it very suitable for optimisation purposes. The maximum relative errors of Q 1 and Q 3 are 3 10% and 2 90% respectively. . .

In order to get more insight in the numerical accuracy of the approximation for a wide variety of different parameter settings, we create a large test bed in the next subsection and compare the approximation with exact or simulated results. It turns out that the maximum relative errors for most of the polling systems are smaller than the one selected in the above example.

Figure 1: Approximated and simulated mean waiting time E [ W 2 ] of Q 2 of the example in subsection 4.1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0134v1-figure-2327208e-0847-45e8-9f99-eb6287b05be4.jpg)

## 4.2 Accuracy of the approximation

In the present section we study the accuracy of our approximation. We compare the approximated mean waiting times of customers in various polling systems to the exact values. The complete test bed of polling systems that are analysed, contains 2304 different combinations of parameter values, all listed in Table 1. We show detailed results for exhaustive service first, and discuss polling systems with gated service at the end of this section. We have varied

Table 1: Test bed used to compare the approximation to exact results.

| Parameter                           | Notation   | Values                                         |
|-------------------------------------|------------|------------------------------------------------|
| Number of queues                    | N          | 2 , 3 , 4 , 5                                  |
| Load                                | ρ          | 0 . 1 , 0 . 3 , 0 . 5 , 0 . 7 , 0 . 9 , 0 . 99 |
| SCV interarrival times              | cv 2 A i   | 0 . 25 , 1 , 2                                 |
| SCV service times                   | cv 2 B i   | 0 . 25 , 1                                     |
| SCV switch-over times               | cv 2 S i   | 0 . 25 , 1                                     |
| Imbalance interarrival times        | I A i      | 1 , 5                                          |
| Imbalance service times             | I B i      | 1 , 5                                          |
| Ratio service and switch-over times | I S i /B i | 1 , 5                                          |

glyph[negationslash]

the load between 0 1 and 0 9 with steps of 0 2, and included . . . ρ = 0 99 to analyse the limiting . behaviour of our approximation when the load tends to 1. The SCV of the interarrival times, cv 2 A i , is varied between 0 25 and 2. . In case of non-Poisson arrivals, i.e. cv 2 A i = 1, the exact values have been established through extensive simulation because they cannot be obtained in an analytic way. In these simulations we fit a phase-type distribution to the first two moments of the interarrival times, as described in Examples 2 and 3. For service times and switch-over times, only SCVs of 0 25 and 1 are considered. . SCVs greater than 1 are less common in practice and are discussed separately from the test bed later in this section. The imbalance in interarrival times and service times, I A i and I B i , is the ratio between the largest and the smallest mean interarrival/service time. The interarrival times are determined in such a way, that the overall mean is always 1, λ 1 is the largest and λ N the smallest, and the steps between the λ i are linear. E.g., for N = 5 and I A i = 5 we get λ i = 2 -i/ 3 , i = 1 , . . . , 5. The mean service times E [ B i ] increase linearly in i = 1 , . . . , N , with E [ B N ] = I B i E [ B 1 ] (so E [ B 1 ] is the smallest mean service time). They follow from the relation ∑ N i =1 λ i E [ B i ] = ρ . E.g., for N = 5, and I A i = I B i = 5 we get E [ B /ρ i ] = 3 i/ 35. The last parameter that is varied in the test bed, is the ratio between the mean switch-over times and the mean service times, I S /B i i = E [ S i ] E [ B i ] . The total number of systems analysed is 4 × 6 × 3 × 2 5 = 2304. A system consisting of N queues results N mean waiting times, E [ W ,..., 1 ] E [ W N ], so in total these 2304 systems yield 8064 mean waiting times. The absolute relative errors, defined as | o - | e /e , where o stands for observed (approximated) value, and e stands for expected (exact) value, are computed for all these 8064 queues. Table 2 shows these relative errors (times 100%) categorised in bins of 5%. In this table, and in all other tables, results for systems with a different number of queues are displayed in separate rows. The reader should keep in mind that the statistics in each row are based on 1 4 × 2304 × N absolute relative errors, where N is the number of queues used in the specified row. Table 2 shows that, e.g., 98 84% of the approximated mean . waiting times in polling systems consisting of 3 queues deviate less than 5% from their true

values. From Table 2 it can be concluded that the approximation accuracy increases with the number of queues in a polling system. More specifically, for systems with more than 2 queues, no approximation errors are greater than 10%, and the vast majority is less than 5%. The mean relative errors for N = 2 , . . . , 5 are respectively 2 18% 0 93% 0 70% and 0 57%. . , . , . . It is also noteworthy, that 193 out of the 2304 systems yield exact results. All of these 193 systems have Poisson input, and all of them - except for one - are symmetric. The only asymmetric case for which our approximation yields an exact result, happens to satisfy constraints (3.16) and (3.17).

In Table 3 the mean relative error percentages are shown for a combination of input parameter settings. The number of queues is always varied per row, while per column another input parameter is varied. This way we can find in more detail which (combinations of) parameter settings result in large approximation errors. In Table 3(a) the load ρ is varied, and it can be seen that for a load of ρ = 0 7 the approximation is least accurate. . E.g., the mean relative error of all approximated waiting times in polling systems consisting of 3 queues with a load of ρ = 0 7 is 1 69%. . . Table 3(b) shows the impact of the SCV of the interarrival times on the accuracy. Especially for systems with more than 2 queues the accuracy is very satisfactory, in particular for the case cv 2 A i = 1. In Table 3(c) the impact of imbalance in a polling system on the accuracy is depicted, and, as could be expected, it can be concluded that a high imbalance in either service or interarrival times has a considerable, negative, impact on the approximation accuracy. Polling systems with more than 2 queues are much less bothered by this imbalance than polling systems with only 2 queues.

Table 2: Errors of the approximation applied to the 2304 test cases with exhaustive service, as described in Section 4, categorised in bins of 5%.

|   N | 0 - 5%   | 5 - 10%   | 10 - 15%   | 15 - 20%   |
|-----|----------|-----------|------------|------------|
|   2 | 86 . 46  | 10 . 24   | 2 . 78     | 0 . 52     |
|   3 | 98 . 84  | 1 . 16    | 0 . 00     | 0 . 00     |
|   4 | 99 . 78  | 0 . 22    | 0 . 00     | 0 . 00     |
|   5 | 99 . 93  | 0 . 07    | 0 . 00     | 0 . 00     |

## 4.3 Miscellaneous other cases

In this subsection we discuss several cases that are left out of the test bed because they might not give any new insights, or because the combination of parameter values might be rarely found in practice.

More queues. Firstly, we discuss polling systems with more than 5 queues briefly. Without listing the actual results, we mention here that the approximations become more and more accurate when letting N grow larger, and still varying the other parameters in the same way as is described in Table 1. For N = 10 already, all relative errors are less than 5%, with an average of less than 0 5% and it only gets smaller as . N grows further.

More variation in service times and switch-over times. In the test bed we only use SCVs 0 25 and 1 for the service times and switch-over times, because these seem more relevant .

| N        | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   |
|----------|--------------|--------------|--------------|--------------|--------------|
| 0 . 10   | 0 . 30       | 0 . 50       | 0 . 70       | 0 . 90       | 0 . 99       |
| 2 0 . 31 | 1 . 81       | 3 . 41       | 4 . 17       | 2 . 70       | 0 . 67       |
| 3 0 . 16 | 0 . 84       | 1 . 44       | 1 . 69       | 1 . 07       | 0 . 39       |
| 4 0 . 13 | 0 . 68       | 1 . 14       | 1 . 28       | 0 . 73       | 0 . 25       |
| 5 0 . 11 | 0 . 57       | 0 . 94       | 1 . 03       | 0 . 57       | 0 . 22       |

(a)

| N   | SCV interarrival 0 . 25   | 1      | times ( cv A 2   |
|-----|---------------------------|--------|------------------|
| 2   | 2 . 27                    | 1 . 76 | 2 . 50           |
| 3   | 1 . 36                    | 0 . 52 | 0 . 92           |
| 4   | 1 . 13                    | 0 . 29 | 0 . 69           |
| 5   | 0 . 97                    | 0 . 19 | 0 . 56           |
|     |                           | (b)    |                  |

| N   | Imbalance interarrival and service times   | Imbalance interarrival and service times   | Imbalance interarrival and service times   | Imbalance interarrival and service times   |
|-----|--------------------------------------------|--------------------------------------------|--------------------------------------------|--------------------------------------------|
|     | I A i = 1 , I B i = 1                      | I A i = 1 , I B i = 5                      | I A i = 5 , I B i = 1                      | I A i = 5 , I B i = 5                      |
| 2   | 0 . 69                                     | 2 . 92                                     | 2 . 80                                     | 2 . 30                                     |
| 3   | 0 . 65                                     | 1 . 27                                     | 0 . 75                                     | 1 . 06                                     |
| 4   | 0 . 56                                     | 0 . 89                                     | 0 . 62                                     | 0 . 73                                     |
| 5   | 0 . 49                                     | 0 . 69                                     | 0 . 53                                     | 0 . 59                                     |

(c)

Table 3: Mean relative approximation error, categorised by number of queues ( N ) and total load of the system (a), SCV interarrival times (b), and imbalance of the interarrival and service times (c).

from a practical point of view. As the coefficient of variation grows larger, our approximation will become less accurate. E.g., for Poisson arrivals we took cv 2 B i ∈ { 2 5 , } , cv 2 S i ∈ { 2 5 , } , and varied the other parameters as in our test bed (see Table 1). This way we reproduced Table 2. The result is shown in Table 4 and indicates that the quality of our approximation deteriorates in these extreme cases. The mean relative errors for N = 2 , . . . , 5 are respectively 3 58%, 1 78%, 1 07%, and 0 77%, which is still very good for systems with such high variation . . . . in service times and switch-over times. For non-Poisson input, no investigations were carried out because the results are expected to show the same kind of behaviour.

Table 4: Errors of the approximation applied to the 768 test cases with Poisson arrival processes and high SCVs of the service times and switch-over times, categorised in bins of 5%.

|   N | 0 - 5%   | 5 - 10%   | 10 - 15%   | 15 - 20%   |
|-----|----------|-----------|------------|------------|
|   2 | 74 . 22  | 14 . 84   | 6 . 51     | 2 . 08     |
|   3 | 89 . 76  | 7 . 29    | 2 . 08     | 0 . 69     |
|   4 | 94 . 53  | 4 . 56    | 0 . 91     | 0 . 00     |
|   5 | 97 . 71  | 2 . 19    | 0 . 10     | 0 . 00     |

Small switch-over times. Systems with small switch-over times, in particular smaller than the mean service times, also show a deterioration of approximation accuracy - especially in systems with 2 queues. In Figure 2 we show an extreme case with N = 2, service times and switch-over times are exponentially distributed with E [ B i ] = 9 40 and E [ S i ] = 9 200 for i = 1 2, which makes the mean switch-over times 5 times , smaller than the mean service times. Furthermore, the interarrival times are exponentially distributed with λ 1 = 5 λ 2 . In Figure 2 the mean waiting times of customers in both queues are plotted versus the load of the system. Both the approximation and the exact values are plotted. For customers in Q 1 the mean waiting time approximations underestimate the true values, which leads to a maximum relative error of -11 2% for . ρ = 0 7 ( . E [ W 1 , app ] = 0 43, whereas . E [ W 1 ] = 0 49). . The approximated mean waiting time for customers in Q 2 is systematically overestimating the true value. The maximum relative error is attained at ρ = 0 5 and is 28 8% ( . . E [ W 1 , app ] = 0 41, . whereas E [ W 1 ] = 0 52). . Although the relative errors are high in this situation, the absolute errors are still rather small compared to the mean service time of an individual customer. This implies that the mean sojourn time is already much better approximated. Nevertheless, this example illustrates one of the situations where our approximation gives unsatisfactory results.

Figure 2: Approximated and exact mean waiting times for a two-queue polling system with small switch-over times.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0134v1-figure-adb81657-6294-461f-93f5-5f9cb2f2a8a6.jpg)

## 4.4 Comparison with existing approximations

The only alternative approximations that exist for polling systems with non-exponential interarrival times, only perform well under extreme, limiting conditions. In [17, 23] it is suggested to use the HT limit (3.14) as an approximation, but the accuracy is only found to be acceptable for ρ &gt; 0 8. . Another approximation for the mean waiting time in polling systems with non-exponential interarrival times uses the limit for S →∞ [27, 28]. This approximation is usable if either the total switch-over time in the system is large and the switch-over times have low variance, or if the total switch-over time in the system is large and the system is in heavy traffic. For completeness, we mention that it is also possible to construct an approximation that is purely based on the LT limit, developed in the present paper:

$$\mathbb { E } [ W _ { i } ^ { L T } ] \approx \frac { K _ { 0, i } + ( K _ { 1, i } - K _ { 0, i } ) \rho } { 1 - \rho }.$$

We do not wish to go into further details on this topic, because the accuracy of (4.1) turns out to be worse for high loads. This makes our approximation, which is exact in all these

limiting cases, the only one which can be applied under all circumstances. We support this statement by reproducing Table 3(a) for the three alternative approximations, based on HT limit, large switch-over times, and LT limit. A comparison of Table 5, which displays these results, to Table 3(a), clearly indicates the drawbacks of using approximations based on one limiting case only.

| N         | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   |
|-----------|--------------|--------------|--------------|--------------|--------------|
| 0 . 10    | 0 . 30       | 0 . 50       | 0 . 70       | 0 . 90       | 0 . 99       |
| 2 58 . 97 | 51 . 03      | 40 . 34      | 27 . 47      | 11 . 34      | 1 . 46       |
| 3 30 . 62 | 25 . 66      | 19 . 63      | 12 . 61      | 4 . 56       | 0 . 62       |
| 4 22 . 49 | 18 . 61      | 14 . 02      | 8 . 80       | 3 . 08       | 0 . 42       |
| 5 18 . 06 | 14 . 83      | 11 . 10      | 6 . 88       | 2 . 36       | 0 . 37       |

(a)

(b)

| N   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   |
|-----|--------------|--------------|--------------|--------------|--------------|--------------|
|     | 0 . 10       | 0 . 30       | 0 . 50       | 0 . 70       | 0 . 90       | 0 . 99       |
| 2   | 27 . 13      | 31 . 08      | 35 . 95      | 41 . 32      | 46 . 78      | 49 . 23      |
| 3   | 19 . 28      | 21 . 68      | 24 . 51      | 27 . 52      | 30 . 33      | 31 . 44      |
| 4   | 15 . 05      | 16 . 75      | 18 . 76      | 20 . 92      | 22 . 93      | 23 . 67      |
| 5   | 12 . 37      | 13 . 67      | 15 . 24      | 16 . 96      | 18 . 56      | 19 . 17      |

Table 5: Mean relative approximation error, categorised by number of queues ( N ) and total load of the system, for three alternative approximations: based on the HT limit (a), based on large switch-over times (b), and based on the LT limit (c).

| N Load ( ρ )   | N Load ( ρ )   | N Load ( ρ )   | N Load ( ρ )   | N Load ( ρ )   | N Load ( ρ )   |
|----------------|----------------|----------------|----------------|----------------|----------------|
|                | 0 . 10         | 0 . 30         | 0 . 50 0 . 70  | 0 . 90         | 0 . 99         |
| 2              | 0 . 34 2 . 65  | 7 . 00         | 13 . 42        | 21 . 90        | 26 . 82        |
| 3              | 0 . 24 1 . 68  | 3 . 90         | 6 . 70         | 9 . 62         | 10 . 82        |
| 4              | 0 . 19 1 . 31  | 2 . 97         | 4 . 87         | 6 . 75         | 7 . 46         |
| 5              | 0 . 16 1 . 09  | 2 . 45         | 4 . 02         | 5 . 54         | 6 . 12         |

For polling systems with Poisson arrivals, several alternative approximations have been developed in existing literature. The best one among them (see, e.g., [3, 7, 11]) uses the relation E [ W i ] = (1 ± ρ i ) E [ C res i ], where C i is the cycle time, starting at a visit completion to Q i when service is exhaustive, and starting at a visit beginning for gated service. By ± we mean -for exhaustive service, and + for gated service. The mean residual cycle time, E [ C res i ], is assumed to be equal for all queues, i.e. E [ C res i ] ≈ E [ C res ], and can be found by substituting E [ W i ] ≈ (1 ± ρ i ) E [ C res ] in the pseudo-conservation law (3.18). We have used this PCL-based approximation to estimate the mean waiting times of all queues in the test bed described in Table 1, but taking only the 768 cases where cv 2 A i = 1. Table 6 shows the mean relative

errors for our approximation (a) and the PCL approximation (b), categorised in bins of 5% as was done before in Table 2. From these tables (and from other performed experiments that are not mentioned for the sake of brevity) it can be concluded that for N &gt; 2 both approximations have almost the same accuracy, our approximation being slightly better for small values of ρ , and the PCL approximation being slightly better for high values of ρ (both methods are asymptotically exact as ρ ↑ 1). However, for N = 2 our method suffers greatly from imbalance in the system, whereas the PCL approximation proves to be more robust.

Table 6: Errors of the approximation applied to the 768 test cases with Poisson input, categorised in bins of 5%. In (a) the percentages of mean relative errors in each bin are shown for our approximation, in (b) results are shown for the PCL approximation.

| N   | 0 - 5%   | 5 - 10%   | 10 - 15%   |
|-----|----------|-----------|------------|
| 2   | 89 . 32  | 9 . 11    | 1 . 56     |
| 3   | 100 . 00 | 0 . 00    | 0 . 00     |
| 4   | 100 . 00 | 0 . 00    | 0 . 00     |
| 5   | 100 . 00 | 0 . 00    | 0 . 00     |
|     |          | (a)       |            |

| N   | 0 - 5%   | 5 - 10%   | 10 - 15%   |
|-----|----------|-----------|------------|
| 2   | 96 . 09  | 2 . 86    | 1 . 04     |
| 3   | 99 . 31  | 0 . 69    | 0 . 00     |
| 4   | 100 . 00 | 0 . 00    | 0 . 00     |
| 5   | 100 . 00 | 0 . 00    | 0 . 00     |
|     |          | (b)       |            |

## 4.5 Gated service

glyph[negationslash]

Until now we have only shown and discussed approximation results for polling systems with exhaustive service. The complete test bed described in Table 1 has also been analysed for polling systems where each queue receives gated service. As can be seen in Table 7, the overall quality of the approximation is good, but worse than for polling systems with exhaustive service. More details on the reason for these inaccuracies can be found in Table 8, which is the equivalent of Table 3 for gated service. Table 8(b) illustrates that there is now a huge difference between systems with Poisson arrivals, and systems with non-Poisson arrivals. For the cases with cv 2 A i = 1, the approximation is extremely accurate, even for two-queue polling systems. The accuracy in cases with cv 2 A i = 1 is worse, which is caused by the assumptions that are made to approximate the LT limit (3.12). Firstly, the decomposition (3.3) does not hold for non-Poisson arrivals, and secondly, the terms E [ I res i ] and E [ V I i i ] E [ I i ] in this decomposition have only been approximated. For exhaustive service, these assumptions do not have much negative impact on the accuracy, but apparently, for gated service, they do. The mean relative errors for N = 2 , . . . , 5 queues are respectively 2 70%, 2 25%, 1 90%, and 1 63%. The . . . . imbalance of the mean interarrival and service times hardly influences the accuracy of the approximation, as can be concluded from Table 8(c).

If we consider the 768 cases with Poisson arrivals only, the mean relative errors of our approximation for N = 2 , . . . , 5 are respectively 0 34% . , 0 17% 0 10%, and 0 08%. This accuracy . , . . is even better than the one achieved by the PCL approximation.

|   N | 0 - 5%   | 5 - 10%   | 10 - 15%   | 15 - 20%   |
|-----|----------|-----------|------------|------------|
|   2 | 82 . 55  | 12 . 33   | 2 . 95     | 1 . 56     |
|   3 | 85 . 42  | 10 . 53   | 3 . 13     | 0 . 81     |
|   4 | 88 . 85  | 8 . 46    | 2 . 43     | 0 . 26     |
|   5 | 92 . 22  | 6 . 60    | 1 . 15     | 0 . 03     |

Table 7: Errors of the approximation applied to the 2304 test cases with gated service, as described in Subsection 4.5, categorised in bins of 5%.

| N        | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   | Load ( ρ )   |
|----------|--------------|--------------|--------------|--------------|--------------|
| 0        | . 10 0 . 30  | 0 . 50       | 0 . 70       | 0 . 90       | 0 . 99       |
| 2 2 .    | 64 4 . 55    | 4 . 31       | 3 . 10       | 1 . 25       | 0 . 37       |
| 3 2 . 03 | 3 . 78       | 3 . 68       | 2 . 68       | 1 . 04       | 0 . 30       |
| 4 1 . 62 | 3 . 14       | 3 . 13       | 2 . 32       | 0 . 92       | 0 . 28       |
| 5 1 .    | 35 2 . 67    | 2 . 71       | 2 . 03       | 0 . 81       | 0 . 21       |

(a)

Table 8: For gated service: mean relative approximation error, categorised by number of queues ( N ) and total load of the system (a), SCV interarrival times (b), and imbalance of the interarrival and service times (c).

| N   | SCV interarrival times ( cv 2 A i )   | SCV interarrival times ( cv 2 A i )   | SCV interarrival times ( cv 2 A i )   |
|-----|---------------------------------------|---------------------------------------|---------------------------------------|
|     | 0 . 25                                | 1                                     | 2                                     |
| 2   | 4 . 72                                | 0 . 34                                | 3 . 05                                |
| 3   | 4 . 06                                | 0 . 17                                | 2 . 53                                |
| 4   | 3 . 45                                | 0 . 10                                | 2 . 16                                |
| 5   | 2 . 98                                | 0 . 08                                | 1 . 84                                |
|     |                                       | (b)                                   |                                       |

| N   | Imbalance interarrival and service times   | Imbalance interarrival and service times   | Imbalance interarrival and service times   | Imbalance interarrival and service times   |
|-----|--------------------------------------------|--------------------------------------------|--------------------------------------------|--------------------------------------------|
|     | I A i = 1 , I B i =                        | 1 I A i = 1 , I B i = 5                    | I A i = 5 , I B i = 1                      | I A i = 5 , I B i = 5                      |
| 2   | 2 . 76                                     | 2 . 64                                     | 2 . 81                                     | 2 . 59                                     |
| 3   | 2 . 28                                     | 2 . 25                                     | 2 . 27                                     | 2 . 21                                     |
| 4   | 1 . 93                                     | 1 . 91                                     | 1 . 90                                     | 1 . 87                                     |
| 5   | 1 . 64                                     | 1 . 66                                     | 1 . 64                                     | 1 . 58                                     |
| (c) | (c)                                        | (c)                                        | (c)                                        | (c)                                        |

## 5 Further research topics

The research that is done in the present paper can be extended in many different directions. In this section we discuss some possibilities that we find most relevant.

Higher moments. Firstly, a logical follow-up step would be to use the same approach to find approximations for higher moments of the waiting time distribution as well. This might prove to be a hard exercise, since the LT limit of E [ W 2 i ] is unknown and, although its derivation might follow the same lines as in Section 3, it probably requires substantially more effort. In [13], explicit expressions for the second moments of the waiting time distributions are given, but only for symmetric systems with N = 2 3 , , and 4, and under the assumption of Poisson arrivals. Also, the HT limit of E [ W 2 i ] is unknown, although some research in this area

has already been done and in [17] a strong conjecture is given for the limiting distribution of W i as ρ ↑ 1.

Another question that remains to be investigated, is the required form of the interpolation, as (2.1) is surely not adequate to approximate higher moments of E [ W i ].

Other service disciplines. In the present paper, only exhaustive and gated service are discussed. In order to obtain results for polling systems with some queues receiving exhaustive service, and others receiving gated service, only minor modifications should be made, but we leave this to the reader. It would be more challenging to generalise the approximation to a wider variety of service disciplines. In particular, it would be nice to have one expression for the mean waiting time of customers in a queue with an arbitrary branching-type service discipline (cf. [19]). The exhaustiveness of a branching-type service discipline (cf. [26]) might appear in this expression. Gated and exhaustive are both branching type service disciplines, but are discussed separately in the present paper. The HT limit can most likely be established for arbitrary branching type service disciplines (see conjectures in [17]), so the question that remains is whether the LT limit can be found in a similar way.

Optimisation. One of the main reasons to choose (2.1) as form of the interpolation, besides its asymptotic correctness, is its simplicity. Having this exact and simple expression for the approximate mean waiting times, makes it very useful for optimisation purposes. In production environments, one can, for example, determine what the optimal strategy is to combine orders of different types (i.e., determine what queue customers should join). Because general arrivals are supported, one can determine optimal sizes of batches in which items are grouped and sent to a specific machine. The simplicity of (2.1) makes it possible for a manager to create a handy Excel sheet that can be used by operators to compute all kind of optimal parameter settings. No difficult computations are required at all, so a large variety of users can use the approximation.

In the present paper the accuracy of the approximation has been investigated and has been found to be very good in most situations. Another advantage of our approximation regarding optimisation purposes, is that the general shape of the approximated curve follows the exact curve very closely. Even in cases where the relative errors are rather large, like in Figure 1, the shape of the actual curves is still very well approximated. This means that plugging our approximation, instead of an exact expression if it had been available, in an optimisation function yields an optimum that should be close to the true optimum.

Polling Table. The interpolation based approximation can also be extended to polling systems where the visiting order of the queues is not cyclic. Waiting times in polling systems with so-called polling tables can be obtained in the same way as shown in the present paper. Both the LT and HT limits are not difficult to determine in this situation, and the interpolation follows directly from these limits.

Model. The form of the interpolation might be changed to improve the accuracy of approximations for cases that give less satisfactory results in the present form. E.g., one could try other functions than a second-order polynomial as numerator of (2.1). Alternatively, one

could try to find a correction term which could be added to (2.1) to obtain better results for, e.g., two-queue polling systems. But most of all, if an exact LT limit of the mean waiting time in a polling system with non-Poisson arrivals could be found, the accuracy of the approximation in the case of gated service might be improved.

## Acknowledgements

The authors wish to thank Onno Boxma for valuable discussions and for useful comments on earlier drafts of the present paper.

## References

- [1] D. Bertsimas and G. Mourtzinou. Decomposition results for general polling systems and their applications. Queueing Systems , 31:295-316, 1999.
- [2] J. P. C. Blanc and R. D. van der Mei. Optimization of polling systems with Bernoulli schedules. Performance Evaluation , 22:139-158, 1995.
- [3] O. J. Boxma. Workloads and waiting times in single-server systems with multiple customer classes. Queueing Systems , 5:185-214, 1989.
- [4] O. J. Boxma and W. P. Groenendijk. Pseudo-conservation laws in cyclic-service systems. Journal of Applied Probability , 24(4):949-964, 1987.
- [5] E. G. Coffman, Jr., A. A. Puhalskii, and M. I. Reiman. Polling systems with zero switchover times: A heavy-traffic averaging principle. The Annals of Applied Probability , 5(3):681-719, 1995.
- [6] E. G. Coffman, Jr., A. A. Puhalskii, and M. I. Reiman. Polling systems in heavy-traffic: A Bessel process limit. Mathematics of Operations Research , 23:257-304, 1998.
- [7] D. Everitt. Simple approximations for token rings. IEEE Transactions on Communications , COM-34(7):719-721, 1986.
- [8] M. J. Fischer, C. M. Harris, and J. Xie. An interpolation approximation for expected wait in a time-limited polling system. Computers &amp; Operations Research , 27:353-366, 2000.
- [9] P. J. Fleming and B. Simon. Interpolation approximations of sojourn time distributions. Operations Research , 39(2):251-260, 1991.
- [10] S. W. Fuhrmann and R. B. Cooper. Stochastic decompositions in the M/G/ 1 queue with generalized vacations. Operations Research , 33(5):1117-1129, 1985.
- [11] W. P. Groenendijk. Waiting-time approximations for cyclic-service systems with mixed service strategies. In Proc. 12th ITC , pages 1434-1441. North-Holland Publ. Co., Amsterdam, 1989.

| [12]   | J. Keilson and L. D. Servi. The distributional form of Little's Law and the Fuhrmann- Cooper decomposition. Operations Research Letters , 9(4):239-247, 1990.                    |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [13]   | S. Kudoh, H. Takagi, and O. Hashida. Second moments of the waiting time in symmetric polling systems. Journal of the Operations Research Society of Japan , 43(2):306-316, 1996. |
| [14]   | H. Levy and M. Sidi. Polling systems: applications, modeling, and optimization. IEEE Transactions on Communications , 38:1750-1760, 1990.                                        |
| [15]   | T. L. Olsen. Limit theorems for polling models with increasing setups. Probability in the Engineering and Informational Sciences , 15(1):35-55, 2001.                            |
| [16]   | T. L. Olsen and R. D. van der Mei. Polling systems with periodic server routeing in heavy traffic: distribution of the delay. Journal of Applied Probability , 40:305-326, 2003. |
| [17]   | T. L. Olsen and R. D. van der Mei. Polling systems with periodic server routing in heavy traffic: renewal arrivals. Operations Research Letters , 33:17-25, 2005.                |
| [18]   | M. I. Reiman and B. Simon. An interpolation approximation for queueing systems with Poisson input. Operations Research , 36(3):454-469, 1988.                                    |
| [19]   | J. A. C. Resing. Polling systems and multitype branching processes. Queueing Systems , 13:409 - 426, 1993.                                                                       |
| [20]   | B. Simon. A simple relationship between light and heavy traffic limits. Operations Research , 40(Supplement 2):S342-S345, 1992.                                                  |
| [21]   | H. Takagi. Queuing analysis of polling models. ACM Computing Surveys (CSUR) , 20: 5-28, 1988.                                                                                    |
| [22]   | H. C. Tijms. Stochastic models: an algorithmic approach . Wiley, Chichester, 1994.                                                                                               |
| [23]   | R. D. van der Mei and E. M. M. Winands. A note on polling models with renewal arrivals and nonzero switch-over times. Operations Research Letters , 36:500-505, 2008.            |
| [24]   | V. M. Vishnevskii and O. V. Semenova. Mathematical methods to study the polling systems. Automation and Remote Control , 67(2):173-220, 2006.                                    |
| [25]   | W. Whitt. An interpolation approximation for the mean workload in a GI/G/ 1 queue. Operations Research , 37(6):936-952, 1989.                                                    |
| [26]   | E. M. M. Winands. Polling, Production & Priorities . PhD thesis, Eindhoven University of Technology, 2007.                                                                       |
| [27]   | E. M. M. Winands. On polling systems with large setups. Operations Research Letters , 35:584-590, 2007.                                                                          |
| [28]   | E. M. M. Winands. Branching-type polling systems with large setups. To appear in OR Spectrum , 2009.                                                                             |
Marko Boon Erik Winands Ivo Adan Sandra van Wijk
Closed-Form Waiting Time Approximations for Polling Systems
A typical polling system consists of a number of queues, attended by a single
server in a fixed order. The vast majority of papers on polling systems
focusses on Poisson arrivals, whereas very few results are available for
general arrivals. The current study is the first one presenting simple
closed-form approximations for the mean waiting times in polling systems with
renewal arrival processes, performing well for ALL workloads. The
approximations are constructed using heavy traffic limits and newly developed
light traffic limits. The closed-form approximations may prove to be extremely
useful for system design and optimisation in application areas as diverse as
telecommunication, maintenance, manufacturing and transportation.