## A Polling Model with Reneging at Polling Instants ∗

## M.A.A. Boon †

marko@win.tue.nl

April, 2010

## Abstract

In this paper we consider a single-server, cyclic polling system with switch-over times and Poisson arrivals. The service disciplines that are discussed, are exhaustive and gated service. The novel contribution of the present paper is that we consider the reneging of customers at polling instants. In more detail, whenever the server starts or ends a visit to a queue, some of the customers waiting in each queue leave the system before having received service. The probability that a certain customer leaves the queue, depends on the queue in which the customer is waiting, and on the location of the server. We show that this system can be analysed by introducing customer subtypes, depending on their arrival periods, and keeping track of the moment when they abandon the system. In order to determine waiting time distributions, we regard the system as a polling model with varying arrival rates, and apply a generalised version of the distributional form of Little's law. The marginal queue length distribution can be found by conditioning on the state of the system (position of the server, and whether it is serving or switching).

Keywords: Polling, reneging, varying arrival rates, queue lengths, waiting times

## 1 Introduction

A polling system is a queueing system that consists of multiple queues being served by one server, generally in a fixed, cyclic, order. There is a vast literature on polling systems, motivated by many real-life applications. These applications are frequently found in production environments, where one machine produces different part types. Typically, after the production of several parts of the same type, the machine is reconfigured and starts producing parts of the next type. The performance measures of interest are, e.g., the mean throughput of the machine, the number of different product orders that are waiting to be processed, and the mean order lead time (i.e., the time between the placement of the order and the completion of the last item in the order). Other typical application areas of polling systems are telecommunications, where several protocols use a round-robin principle for the communication of data packets between multiple devices, and transportation. We recommend surveys of, e.g.,

∗ 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

Takagi [24], Levy and Sidi [18] and Vishnevskii and Semenova [25], for a better overview of applications of polling systems, and techniques to analyse them.

In queueing systems, waiting is an inevitable nuisance. When waiting times become too large, the impatience grows and customers might decide to leave the queue and possibly return another time. This phenomenon, which occurs in many real-life situations, is referred to in queueing literature as reneging . Alternatively, the terms abandonment or impatience are used. The first paper on the subject of abandonment due to impatience, has been written by Palm [21] who studies annoyance of customers in telecommunications. In [21], as in most of the literature on reneging, impatience is modelled as a timer that starts running at the moment that a customer joins the queue. When this timer reaches a certain (usually random) value while the customer is still waiting in the queue, he abandons the system immediately. The system studied in [21] is an M/M/n queue with exponentially distributed customer patience, which is used frequently to model call centers, and is referred to as the Erlang-A ( M/M/n + M ) queueing system. This model, and generalisations to other patience distributions, are studied in more detail in by, e.g., Riordan [23], Haugen and Skogan [15], Baccelli and Hebuterne [4], and Boxma and De Waal [8]. The vast majority of papers on reneging focusses on the application to call centers, studying the loss probability and comparing different staffing rules. See [17, 20, 28] for some recent developments in reneging in queueing models for call centers. In the application area of computer systems, processor sharing is an important discipline and reneging has been studied in this context as well. Assaf and Haviv [3] consider a model where customers may decide to abandon the system, depending on the number of customers that are in service simultaneously. Gromoll et al. [14] consider an overloaded processor sharing queue with impatient customers, and find a scaling procedure that makes the model analytically tractable.

One common aspect of all models considered in the aforementioned literature, is that whenever customers are available in the system, (at least) one of them is in service. In polling models, customer impatience might be increased by absence of the server at the queue of arrival. This kind of behaviour has been studied in single-queue systems by Altman and Yechiali [2], who discuss M/G/ 1 and M/M/c queues with server vacations, and customers growing impatient while the server is away. Zhang et al. [29] study a similar system, but with an M/M/ /N 1 queue. Their model also includes balking, which means that customers may decide not to enter the system at all, depending on the number of customers present in the queue. Madan [19] studies a system where server vacations may start at arbitrary moments, even when customers are being served or when the system is idle. Whenever a vacation starts, a random number of customers abandons the system immediately. The fact that more than one customer can leave the system at the same time, and only at specified moments (in this case the beginning of a vacation) makes this model quite different from most of the other papers dealing with customer impatience. We refer to this kind of abandonment as synchronised reneging , a term that is introduced by Adan et al. [1], who consider a model where each customer has the same probability of abandoning the system at synchronised reneging epochs. They consider a queueing system with server vacations that start as soon as the queue becomes empty, distinguishing between two cases. In the first case, which is called the Unique Abandonment Epoch (UAE) model, customers leave the queue at visit beginnings only. In the second case, referred to as the Multiple Abandonment Epochs (MAE) model, impatient customers also abandon at (randomly) specified, synchronised moments during the server vacation.

In the present paper we study synchronised reneging in polling systems , with abandonments taking place at visit beginnings and endings. Basically, this means that we extend the UAE model of [1] to systems with multiple queues, thereby increasing the number of synchronised reneging epochs. The higher frequency of abandonment epochs significantly increases the complexity of the analysis. Although the reneging policy considered in the current paper is based on [1], the analysis is different. It is based on new techniques, developed in a recent paper on polling models with so-called smart customers , cf. [5], to find waiting time and queue length distributions in polling systems with varying arrival rates. In the present paper we use and extend these techniques, so they can be applied to a polling model with synchronised reneging. This makes it possible to extend the results of [1] in three new directions. Firstly, this allows us to study systems with multiple queues. Secondly, we can consider other service disciplines than exhaustive service (i.e., serve customers in a queue until it is empty). In the present paper we discuss not only exhaustive, but also gated service (i.e., serve all customers present at the server's arrival at the queue). The third new contribution of the present paper is that we can compute other relevant performance measures as well, like the distributions of the cycle times and the waiting times. These extensions provide the main motivation to study this model, which is one of the first attempts to introduce reneging in polling systems. The only related work is by Vishnevsky and Semenova [26] who study a two-queue polling system with exponential service times and exhaustive service, with the more conventional timer to model the patience of the customer. They illustrate how the Power-series algorithm can be used to find the equilibrium state probabilities, but no explicit performance measures are computed.

Since the analysis in the present paper relies heavily on [5], we briefly summarise their results. In a polling system with smart customers, the arrival rates of the different customer types depend on the state of the server, where state is defined as a combination of its location, and whether it is working or switching. In this situation, it is no longer possible to apply the distributional form of Little's law in its standard form (see, for example, [16]), but it requires a generalization, developed in [5]. In the present paper we apply this model with smart customers to determine the Laplace-Stieltjes Transform (LST) of the waiting time distribution of each customer type. This generalised version of the distributional form of Little's law is applied to the joint queue length distribution at departure epochs of customers that have not abandoned the system prematurely, which means that this waiting time is determined only for customers that are actually served. To determine the Probability Generating Function (PGF) of the marginal queue length distribution of each customer type, we do need to take into account the impatient customers that abandon the system before being served. This requires a different approach, as will be shown later in this paper.

The structure of the present paper is as follows. In the next section we describe the model and the notation in more detail. In Section 3 we introduce an alternative model with smart customers that is used to determine the cycle time and waiting time of served customers in the original model. The stability condition is also presented in this section. The marginal queue length distributions are studied in Section 4, because this requires a different approach. Section 5 discusses a special case of the model under consideration, a polling system with only one queue and exhaustive service. This queueing system with server vacations, has been studied in [1]. We show that the queue length PGF obtained in Section 4 agrees with their result, and we mention some further results that have not been discussed in [1], like the cycle time and sojourn times of all customers, which reduce to elegant expressions when the system

consists of only one queue. The last section discusses a numerical example to illustrate typical features of a polling model with synchronised reneging at polling instants.

## 2 Model and notation

The polling model under consideration contains N queues, Q , . . . , Q 1 N . These queues are served by one server in a fixed, cyclic order. The time that is required to switch from Q i to Q i +1 is denoted by S i , which is called a switch-over time, with LST σ i ( ). · Throughout the paper, all indices are modulo N , so Q N +1 refers to Q 1 and so on. The arrival process of customers in Q i , denoted by type i customers, is a Poisson process with parameter λ i . The service time of a type i customer is denoted by B i , with LST β i ( ). · The switch-over times, interarrival times and service times are all assumed to be independent of each other. The service discipline of each queue determines when the server switches to the next queue. The following property, which is defined in [22] and [12], plays a key role in the analysis of polling systems.

Property 2.1 (Branching Property) If the server arrives at Q i to find k i customers there, then during the course of the server's visit, each of these k i customers will effectively be replaced in an i.i.d. manner by a random population having probability generating function h i ( z , . . . , z 1 N ), which can be any N -dimensional probability generating function.

glyph[negationslash]

Performance measures, like queue length distributions and waiting times, can be determined for polling systems with all queues satisfying Property 2.1, whereas only very few, exceptional, polling systems can be analysed if the service disciplines do not satisfy this property. In the present paper we discuss the two most common service disciplines satisfying this branching property , exhaustive and gated service. A queue with exhaustive service is served until it is completely empty. In a queue with gated service only those customers are served, that are present at the beginning of a visit to this queue. The PGF h i ( z , . . . , z 1 N ) in Property 2.1 is β i ( ∑ N j =1 λ j (1 -z j ) ) for gated service, and π i ( ∑ j = i λ j (1 -z j ) ) for exhaustive service, where π i ( ) is the LST of a busy period distribution in an · M/G/ 1 system with only type i customers, so it is the root in (0 , 1] of the equation π i ( ω ) = β i ( ω + λ i (1 -π i ( ω ))), ω ≥ 0 (cf. [10], p. 250). Define θ i ( ) as the LST of the time that the server spends at · Q i due to the presence of one customer there. For gated service θ i ( ) = · β i ( ), and for exhaustive service · θ i ( ) = · π i ( ). ·

A cycle consists of the visit times of all queues, denoted by V , . . . , V 1 N , and the switch-over times S , . . . , S 1 N . The distribution of the length of one cycle depends on the starting point of this cycle. We use the notation C i for the time between two successive visit beginnings to Q i , with LST γ i ( ), and · C ∗ i for the time between two successive visit completions to Q i , with LST γ ∗ i ( ). · When studying a queue with gated service, it turns out that C i plays an important role, whereas C ∗ i is used in the analysis of queues receiving exhaustive service. The intervisit time I i is the time between a visit completion of Q i and the next visit beginning at Q i .

The model discussed in the present paper is different from models in existing literature because customers grow impatient and may decide to leave the waiting line before actually being served. This is called reneging . The moments at which customers are allowed to leave the system, are only those moments when a new visit or switch-over time starts. For this reason, we refer to this model as a polling model with synchronised reneging at polling instants .

Now we give a more formal description. As stated before, a cycle consists of the periods V , S 1 1 , . . . , V N , S N . Now let P ∈ { V , S 1 1 , . . . , V N , S N } . At the moment that P starts, each customer waiting in Q i immediately leaves the system with probability p ( P ) i , i = 1 , . . . , N . We denote the probability that a customer stays by q ( P ) i = 1 -p ( P ) i . The difficulty in the analysis of this model, is that customers in a certain queue may leave the system at more than just one occasion. We use different ways to circumvent this problem in order to find the performance measures of interest. An important part in the analysis, is the fact that we artificially split each visit time and switch-over time into two parts, a and b . Visit time V i is split into V ia and V ib . We consider V ia as the subperiod in which all the customers abandon the system, right before the start of the actual service of the type i customers, which takes place during V ib . So during V ia , first each type 1 customer abandons the system with probability p ( V i ) 1 , followed by the type 2 customers, and so on, until the reneging of the type N customers. This requires no time, so E [ V ia ] = 0. During V ib the service of the type i customers that remain in the system takes place, so E [ V ib ] = E [ V i ]. Similarly, the switch-over time S i is also split into S ia and S ib , with E [ S ia ] = 0, and E [ S ib ] = E [ S i ]. During S ia the reneging of the customers that abandon the system before the beginning of S i takes place. We need this way of looking at the system, because the queue lengths at the beginning of V ia are different from the queue lengths at the beginning of V ib (and the same holds for the switch-over times). In fact, one can regard V ia , for i = 1 , . . . , N , as separate visit periods during which subsequently type 1 , . . . , N customers are served with probability p ( V i ) 1 , . . . , p ( V i ) N (or probability p ( S i ) 1 , . . . , p ( S i ) N for S ia ), with all service times equal to 0.

## 3 Cycle time, visit times and waiting time distributions

In the present section we study the LSTs of the cycle time distribution, visit time distributions, and of the waiting time distribution of each customer type. The section ends with a note on the stability condition of the model. The waiting time of a customer is the time between the moment of arrival, and the moment that the customer is taken into service. The waiting time is only determined for customers that have not prematurely abandoned the system. The time that all , including reneging, customers spend in the system, requires a different approach and is not discussed in detail in the present paper. We only show in an example in Section 5 how this can be done.

In this section we introduce a different way of looking at the system. Obviously, the length of a visit V i is solely determined by those type i customers that have not abandoned the system at any of the moments where this had been possible. A logical consequence is that only customers that are eventually served, contribute to the cycle time and determine whether the system is stable or not. This observation forms the basis of the analysis in this section. If we remove reneging customers from the system and focus on the remaining customers only, we can show that the system can be viewed as a polling system where the arrival rates of the N customer types depend on the state of the server, i.e., its location and whether it is serving or switching. This type of model is called a polling system with smart customers, introduced in [7], and analysed in more detail in [5]. The present section uses results from these papers and applies them to a polling model with reneging at polling instants.

Westart by introducing the joint queue length PGF at the beginning of all subperiods, denoted

by ˜ LB ( P ) ( z , . . . , z 1 N ), where subperiod P ∈ { V 1 a , V 1 b , S 1 a , S 1 b , . . . , V Na , V Nb , S Na , S Nb } . The PGFs of the joint queue length distributions at the beginnings of the various subperiods in the cycle can be related to each other in the following way:

$$\widetilde { L B } ^ { ( V _ { i b } ) } ( \mathbf z ) = \widetilde { L B } ^ { ( V _ { i a } ) } ( q _ { 1 } ^ { ( V _ { i } ) } z _ { 1 } + p _ { 1 } ^ { ( V _ { i } ) }, \dots, q _ { N } ^ { ( V _ { i } ) } z _ { N } + p _ { N } ^ { ( V _ { i } ) } ),$$

$$\widetilde { L } B ^ { ( S _ { i a } ) } ( \mathbf z ) = \widetilde { L } B ^ { ( V _ { i b } ) } ( z _ { 1 }, \dots, z _ { i - 1 }, h _ { i } ( \mathbf z ), z _ { i + 1 }, \dots, z _ { N } ),$$

$$\widetilde { L } B ^ { ( S _ { i b } ) } ( { \mathbf z } ) = \widetilde { L } B ^ { ( S _ { i a } ) } ( q _ { 1 } ^ { ( S _ { i } ) } z _ { 1 } + p _ { 1 } ^ { ( S _ { i } ) }, \dots, q _ { N } ^ { ( S _ { i } ) } z _ { N } + p _ { N } ^ { ( S _ { i } ) } ),$$

$$\widetilde { L B } ^ { ( V _ { ( i + 1 ) a ) } ) } ( \mathbf z ) = \widetilde { L B } ^ { ( S _ { i b ) } ) } ( \mathbf z ) \, \sigma _ { i } ( \sum _ { j = 1 } ^ { N } \lambda _ { j } ( 1 - z _ { j } ) ),$$

where we use the shorthand notation z for the vector ( z , . . . , z 1 n ). Successive substitution leads to a recursive expression for the joint queue length PGF at an arbitrary polling epoch. In, e.g., [22] it is discussed how this recursive expression leads to the PGF of the joint queue length distribution at polling epochs, written as an infinite product. The recursive equation itself can be used to compute the moments of this joint queue length distribution explicitly. For now, we are more interested in a nice property of the Poisson arrival processes. During the actual visit period V ib , the type i customers that have not abandoned the system are served. For the moment, assuming exhaustive service, we focus on the end of V i , when there are no type i customers present in the system. Then the PGF of the number of type i customers present at the end of S i (which coincides with the beginning of V i +1 ) is

$$\widetilde { L B } ^ { ( V _ { ( i + 1 ) a } ) } ( 1, \dots, 1, z _ { i }, 1, \dots, 1 ) = \sigma _ { i } \left ( \lambda _ { i } ( 1 - z _ { i } ) \right ).$$

Each of these customers abandons the system before the start of V i +1 with probability p ( V i +1 ) i , so the PGF of the number of type i customers that are still in the system at the beginning of V i +1 is:

$$\widetilde { L B } ^ { ( V _ { ( i + 1 ) b ) } ) } ( 1, \dots, 1, z _ { i }, 1, \dots, 1 ) & = \widetilde { L B } ^ { ( V _ { ( i + 1 ) a ) } } ( 1, \dots, 1, q _ { i } ^ { ( V _ { i + 1 } ) } z _ { i } + p _ { i } ^ { ( V _ { i + 1 } ) }, 1, \dots, 1 ) \\ & = \sigma _ { i } \left ( \lambda _ { i } ( 1 - q _ { i } ^ { ( V _ { i + 1 } ) } z _ { i } - p _ { i } ^ { ( V _ { i + 1 } ) } ) \right ) \\ & = \sigma _ { i } \left ( \lambda _ { i } q _ { i } ^ { ( V _ { i + 1 } ) } ( 1 - z _ { i } ) \right ).$$

This short example illustrates that, as far as the joint queue lengths at polling epochs is concerned, and only focussing on the customers that did not abandon the system prematurely, we can view the system as a polling model with Poisson arrivals, but with varying arrival rates (in the example, we have that the new arrival rate is λ q i ( V i +1 ) i during S i ). Just before the start of S i +1 each type i customer in the system reneges with probability p ( S i +1 ) i . This implies that a customer that arrived during S i is still in the system at the beginning of S i +1 with probability q ( V i +1 ) i q ( S i +1 ) i . Hence, the number of type i customers present at the beginning of S i +1 is the same as in a polling system without reneging, but with arrival rates q ( V i +1 ) i q ( S i +1 ) i λ i during S i , and q ( S i +1 ) i λ i during V i +1 . This observation makes it possible to analyse the polling system with reneging by regarding a dual system, without reneging but with varying arrival rates. For the remainder of this section, we consider this dual system. A system with arrival rates that depend on the location of the server, is studied in [5], where it is referred to as a

polling system with smart customers. We apply their results and adopt their notation. Let λ ( P ) i denote the arrival intensity of type i customers during period P ∈ { V , S 1 1 , . . . , V N , S N } . In order to create a similar system as the original polling system with reneging, we define these arrival intensities in the following way:

$$\text{these arrival intensities in the following way:} \\ \lambda _ { i } ^ { ( S _ { i - 1 } ) } & = \lambda _ { i } q _ { i } ^ { ( V _ { i } ) }, \\ \lambda _ { i } ^ { ( V _ { i - 1 } ) } & = \lambda _ { i } q _ { i } ^ { ( S _ { i - 1 } ) } q _ { i } ^ { ( V _ { i } ) }, \\ \lambda _ { i } ^ { ( S _ { i - 2 } ) } & = \lambda _ { i } q _ { i } ^ { ( V _ { i - 1 } ) } q _ { i } ^ { ( S _ { i - 1 } ) } q _ { i } ^ { ( V _ { i } ) }, \\ & \quad \vdots & & & & ( 3. 5 ) \\ \lambda _ { i } ^ { ( S _ { i - N } ) } & = \lambda _ { i } q _ { i } ^ { ( V _ { i } ) } \prod _ { j = 1 } ^ { N - 1 } q _ { i } ^ { ( V _ { i - j } ) } q _ { i } ^ { ( S _ { i - j } ) }, \\ \lambda _ { i } ^ { ( V _ { i - N } ) } & = \begin{cases} \lambda _ { i } \prod _ { j = 1 } ^ { N } q _ { i } ^ { ( V _ { i - j } ) } q _ { i } ^ { ( S _ { i - j } ) }, & \text{ if } Q _ { i } \text{ receives gated service,} \\ \lambda _ { i }, & & & & \text{ if } Q _ { i } \text{ receives exhaustive service.} \end{cases} \\ \text{The only difference for gated service, compared to exhaustive service, is that type } i \text{ customers}$$

The only difference for gated service, compared to exhaustive service, is that type i customers arriving during V i have to wait until the next visit period of type i customers before they are served. Thus, type i customers that arrive during V i , abandon the system before the start of S i with probability p ( S i ) i , whereas for exhaustive service all of these customers would be served during the visit period in which they arrive.

## Cycle time

The cycle time distribution of this dual system is the same as in the original system with reneging. Theorem 5 1 from [5], applied to the dual system with arrival intensities as defined . in (3.5), gives the LSTs of distributions of the cycle time C 1 and the intervisit time I 1 :

$$\mathbb { E } \left [ e ^ { - \omega C _ { 1 } } \right ] & = \widetilde { L B } ^ { ( V _ { 1 b } ) } \left ( \theta _ { 1 } ( \psi ^ { ( V _ { 1 } ) } ( \omega ) ), \dots, \theta _ { N } ( \psi ^ { ( V _ { N } ) } ( \omega ) ) \right ) \prod _ { i = 1 } ^ { N } \sigma _ { i } \left ( \psi ^ { ( S _ { i } ) } ( \omega ) \right ), \\ \mathbb { E } \left [ e ^ { - \omega I _ { 1 } } \right ] & = \widetilde { L B } ^ { ( S _ { 1 b } ) } \left ( 1, \theta _ { 2 } ( \psi ^ { ( V _ { 2 } ) } ( \omega ) ), \dots, \theta _ { N } ( \psi ^ { ( V _ { N } ) } ( \omega ) ) \right ) \prod _ { i = 1 } ^ { N } \sigma _ { i } \left ( \psi ^ { ( S _ { i } ) } ( \omega ) \right ), \\ & \quad \text{--}.$$

where the functions ψ ( P ) ( ω ) are defined in the following, recursive way:

$$\text{where} \, \text{$\psi^{\prime}\nu$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$ } \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \  \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \ text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{ $\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\ varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi`} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$}$$

$$k = i + 1$$

Remark 3.1 Specifically for exhaustive and gated service, more compact expressions for the LSTs of the cycle time and intervisit time distributions are found in Theorem 5 2 in [5]. . These expressions follow from an analysis based on subtypes of customers, where arrivals during the different periods within a cycle mark the customer subtypes. E.g., a type i ( V j ) customer is a type i customer that arrives during V j . Using the analysis based on customer subtypes, we can express the LSTs of the cycle time and intervisit time in terms of the PGF of the joint queue length distribution at polling instants of all customer subtypes, which we denote here as ˜ LB ( V i ) ( z ( V 1 ) 1 , . . . , z ( S N ) 1 , . . . , z ( V 1 ) N , . . . , z ( S N ) N ) . We do not repeat the complete analysis on how to obtain this PGF, but instead refer to Section 5 of [5].

For exhaustive service, the LSTs of the distributions of C ∗ i and I i are:

$$\mathbb { E } [ \mathrm e ^ { - \omega C _ { i } ^ { \ast } } ] = \widetilde { L B } ^ { ( V _ { i } ) } ( 1, \dots, 1, \pi _ { i } ( \omega ) - \frac { \omega } { \lambda _ { i } ^ { ( V _ { 1 } ) } }, \dots, \pi _ { i } ( \omega ) - \frac { \omega } { \lambda _ { i } ^ { ( S _ { N } ) } }, 1, \dots, 1 ), \quad \ ( 3. 6 )$$

$$\mathbb { E } [ e ^ { - \omega I _ { i } } ] = \widetilde { L } B ^ { ( V _ { i } ) } ( 1, \dots, 1, 1 - \frac { \omega } { \lambda _ { i } ^ { ( V _ { 1 } ) } }, \dots, 1 - \frac { \omega } { \lambda _ { i } ^ { ( S _ { N } ) } }, 1, \dots, 1 ).$$

If Q i receives gated service , the LST of the cycle time distribution C i , and the LST of the intervisit time distribution I i , are given by:

$$\Phi \left [ e ^ { - \omega C _ { i } } \right ] & = \widetilde { L } B ^ { ( V _ { i } ) } ( 1, \dots, 1, 1 - \frac { \omega } { \lambda _ { i } ^ { ( V _ { 1 } ) } }, \dots, 1 - \frac { \omega } { \lambda _ { i } ^ { ( S _ { N } ) } }, 1, \dots, 1 ), \\ \mathbb { E } \left [ e ^ { - \omega I _ { i } } \right ] & = \widetilde { L } B ^ { ( V _ { i } ) } ( 1, \dots, 1, 1, 1 - \frac { \omega } { \lambda _ { i } ^ { ( S _ { 1 } ) } }, \dots, 1 - \frac { \omega } { \lambda _ { i } ^ { ( S _ { N } ) } }, 1, \dots, 1 ).$$

$$\cdot$$

## Visit time

The LSTs of the distributions of the visit times V , i i = 1 , . . . , N , can directly be determined for any branching-type service discipline using the function θ i ( ), · and the joint queue length distribution (without subtypes) at the beginning of subperiod V ib :

$$\mathbb { E } [ \mathrm e ^ { - \omega V _ { i } } ] = \widetilde { L B } ^ { ( V _ { i b } ) } ( 1, \dots, 1, \theta _ { i } ( \omega ), 1, \dots, 1 ).$$

The mean cycle time E [ C ] and mean visit times E [ V i ], which are needed later in this paper, can be obtained by differentiating the corresponding LSTs. A numerically more efficient way to compute them, is using MVA for polling systems with smart customers, which is described in more detail in [5].

## Waiting time

The waiting time of customers in the dual system also has the same distribution as the time that customers in the original system have to wait before being taken into service (not taking the impatient customers into account). Note that the marginal queue length distribution at departure epochs of customers that did not renege, is not the same as the marginal queue length distribution at arbitrary epochs , because the arrival intensities change during the cycle. This implies that we cannot use PASTA, and the standard distributional form of Little's law, as discussed by, e.g., Keilson and Servi [16], cannot be used to obtain the waiting time

distribution from the queue length distribution. What we can use though, is a slightly generalised version of the distributional form of Little's law, that can be applied to the joint queue length distribution at departure epochs, as discussed in the proof of Theorem 4 3 in [5]. . We mention the result here, but refer to [5], Section 4, for details on how to obtain the PGF of the joint queue length distributions of all type i customer subtypes at a departure epoch from

$$\begin{array} { l l$$

The LST of the distribution of the waiting time W i of a type i customer, i = 1 , . . . , N , is:

$$\mathbb { E } \left [ e ^ { - \omega W _ { i } } \$$

In the next section we study the queue lengths of all customers that enter the system, including those that renege before the start of their service.

## Stability condition

The stability condition of a polling model with reneging at polling instants, is the same as in the model with smart customers, discussed in the present section. It is clear that customers leaving the system without being served, do not contribute to the workload of the server. The stability condition of a model with smart customers is discussed in [5, 9]. In [9] it is shown that a necessary and sufficient condition for ergodicity is that the Perron-Frobenius eigenvalue of the matrix R -I N should be less than 0, where I N is the N × N identity matrix, and R is an N × N matrix containing elements ρ ij := λ ( V j ) i E [ B i ].

## Proportion of customers served

In queueing models with reneging, the expected proportion of customers that do not abandon the system prematurely, denoted by r , is an important quality measure of the system. In some systems, this might be difficult to compute. In the model considered in the present paper, this quantity is relatively easily obtained. The probability that an arbitrary customer is of type i , is obviously λ / i Λ, where Λ = ∑ N j =1 λ j . The fraction of customers arriving during period P ∈ { V , S 1 1 , . . . , V N , S N } is E [ P / ] E [ C ]. Conditioning on the type of an arbitrary arriving customer, and the location of the server upon his arrival, one can determine the probability that he will not abandon the system prematurely similarly to determining the arrival intensities (3.5). Denote by r i the probability that an arbitrary type i customer is eventually served. Then it is easily seen that

$$\Phi _ { i } & = \sum _ { j = 1 } ^ { N } \left$$

## 4 Queue length distributions

In the previous section, we divided each visit period and switch-over period into two subperiods, part a where impatient customers decide to abandon the system, and part b , where the server is actually serving (or switching in the case of a switch-over period). The PGFs of the joint queue length distributions at the beginnings of all these subperiods are given implicitly by (3.1)-(3.4). In the present section we show how the marginal queue length distribution at an arbitrary epoch can be expressed in terms of these PGFs. We denote the marginal queue length of a type i customer by L i , i = 1 , . . . , N . The PGF of the distribution of L i is determined by conditioning on the subperiod during which the queue is observed. The number of type i customers at an arbitrary moment in subperiod P ∈ { V 1 a , V 1 b , S 1 a , S 1 b , . . . , V Na , V Nb , S Na , S Nb } is denoted by L ( P ) i . By conditioning on P , we have

$$\Re [ z ^ { L _ { i } } ] = \sum _ { j = 1 } ^$$

where we used that E [ V ja ] = E [ S ja ] = 0, E [ V jb ] = E [ V j ], and E [ S jb ] = E [ S j glyph[negationslash]

]. i i

Since S , j j = 1 , . . . , N and V , j j = , are i non-serving intervals for customers of type i , we use a standard result (see, e.g., [6]) to find the PGFs of L ( V jb ) and L ( S jb ) respectively:

glyph[negationslash]

$$\mathbb { E } \left [ z ^ { L _ { i } ^ { ( V _ {$$

$$\mathbb { E } \left [ z ^ { L _ { i } ^ { ( S _ { j b } ) } } \right ] = \frac { \overset { \cdot } { \mathbb { E } } \left [ z ^ { L B _ { i } ^ { ( S _ { j b } ) } } \right ] - \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( V _ { ( j + 1 ) a } ) } } \right ] ^ { \cdot } } { ( 1 - z ) \left ( \mathbb { E } [ L B _ { i } ^ { ( V _ { ( j + 1 ) a } ) } ] - \mathbb { E } [ L B _ { i } ^ { ( S _ { j b } ) } ] \right ) }, \quad \ \ i, j = 1, \dots, N,$$

$$\mathbb { E } \left [ z ^ { L _ { i } ^ { ( V _ { j b } ) } } \right ] & = \frac { \mathbb { E } \left [ z ^ { \cdot \cdot \cdot } \right ] ^ { z } \right ] ^ { z } \right ] ^ { z } \right ] } { ( 1 - z ) \left ( \mathbb { E } [ L B _ { i } ^ { ( S _ { j a } ) } ] - \mathbb { E } [ L B _ { i } ^ { ( V _ { j b } ) } ] \right ) }, & i = 1, \dots, N ; j \neq i, \quad ( 4. 2 ) \\ \mathbb { E } \left [ z ^ { L _ { i } ^ { ( S _ { j b } ) } } \right ] & = \frac { \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( S _ { j b } ) } } \right ] - \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( V _ { j + 1 } ) a } } \right ] } { ( 1 - z ) \left ( \mathbb { E } [ L B _ { i } ^ { ( V _ { ( j + 1 ) a } ) } ] - \mathbb { E } [ L B _ { i } ^ { ( S _ { j b } ) } ] \right ) }, & i, j = 1, \dots, N,$$

where LB ( V jb ) i and LB ( S ja ) i are the number of type i customers at respectively a visit beginning and completion at Q j . Their PGFs are given by ˜ LB ( V jb ) (1 , . . . , 1 , z, 1 , . . . , 1) and ˜ LB ( S ja ) (1 , . . . , 1 , z, 1 , . . . , 1), where z is the element at position i . Differentiation of these PGFs and substituting z = 1 gives the mean values. Similarly, LB ( S jb ) i and LB ( V ( j +1) a ) i are the number of type i customers at respectively the beginning and ending of S j .

It remains to compute E [ z L ( V i ) i ] , i = 1 , . . . , N , i.e. the PGF of the number of type i customers at an arbitrary epoch within V i . In order to do this, we temporarily look at a polling system without reneging, and focus on type i customers. As far as the marginal queue length of type i customers is concerned, the system can be viewed as a vacation queue where the intervisit time I i corresponds to the server vacation. In this 'ordinary' polling model, we can use the Fuhrmann-Cooper decomposition [13], which states that

$$\mathbb { E } [ z ^ { L _ { i } } ] = \frac { ( 1 - \lambda _ { i } \mathbb { E } [ B _ { i } ] ) ( 1 - z ) \beta _ { i } \left ( \lambda _ { i } ( 1 - z ) \right ) } { \beta _ { i } \left ( \lambda _ { i } ( 1 - z ) \right ) - z } \times \frac { \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( S _ { i } ) } } \right ] - \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( V _ { i } ) } } \right ] } { ( 1 - z ) \left ( \mathbb { E } [ L B _ { i } ^ { ( V _ { i } ) } ] - \mathbb { E } [ L B _ { i } ^ { ( S _ { i } ) } ] \right ) }, \quad ( 4. 4 )$$

where LB ( V i ) i and LB ( S i ) i denote the number of type i customers at respectively the beginning and completion of V i . The first part in this decomposition is the PGF of the marginal queue

length of an M/G/ 1 queue with type i customers only. The second part, which is independent of the first, is the PGF of the number of type i customers at an arbitrary epoch during the intervisit time I i , which we denote by E [ z L ( I i ) i ]. Now we focus on the visit and intervisit time of Q i separately, using the relation E [ z L i ] = E [ V i ] E [ C ] E [ z L ( V i ) i ] + E [ I i ] E [ C ] E [ z L ( I i ) i ]. Plugging this relation into (4.4), leads to:

$$\mathbb { E } [ z ^ { L _ { i } ^ { ( V _ { i } ) } } ] = \frac { 1 - \lambda _ { i } \mathbb { E } [ B _ { i } ] } { \lambda _ { i } \mathbb { E } [ B _ { i } ] } \frac { z \left ( 1 - \beta _ { i } ( \lambda _ { i } ( 1 - z ) ) \right ) } { \beta _ { i } ( \lambda _ { i } ( 1 - z ) ) - z } \times \frac { \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( S _ { i } ) } } \right ] - \mathbb { E } \left [ z ^ { L B _ { i } ^ { ( V _ { i } ) } } \right ] } { ( 1 - z ) \left ( \mathbb { E } [ L B _ { i } ^ { ( V _ { i } ) } ] - \mathbb { E } [ L B _ { i } ^ { ( S _ { i } ) } ] \right ) }.$$

The second part of this decomposition is, again, the PGF of the number of customers at an arbitrary moment during the intervisit time I i . The first part can be recognised as the PGF of the queue length of an M/G/ 1 queue with type i customers only, at an arbitrary epoch during a busy period.

Now we return to the model with synchronised reneging. The key observation is that during a visit period , this system behaves exactly as a polling system without reneging. Equation (4.5) no longer depends on anything that happens during the intervisit time, because this is all captured in LB ( V i ) i , the number of type i customers at the beginning of a visit to Q i . This implies that (4.5) also holds for the system considered in the present paper. The only difference is that the interpretation of (4.5) is different. Obviously, the first part in (4.5) still is the PGF of the queue length of an M/G/ 1 queue at an arbitrary epoch during a busy period. However, the last term can no longer be interpreted as the PGF of the distribution of the number of type i customers at an arbitrary moment during the intervisit time I i . For this reason, the Fuhrmann-Cooper decomposition does not hold in a polling model with synchronised reneging. The condition that has to be satisfied for the Fuhrmann-Cooper decomposition, is:

glyph[negationslash]

$$\sum _ { j \neq i } \left ( \frac { \mathbb { E } [ V _ { j } ] } { \mathbb { E } [ C ] } \mathbb { E } \left [ z _ { i } ^ { ( V _ { j } b ) } \right ] \right ) + \sum _ { j = 1 } ^ { N } \left ( \frac { \mathbb { E } [ S _ { j } ] } { \mathbb { E } [ C ] } \mathbb { E } \left [ z _ { i } ^ { ( S _ { j } b ) } \right ] \right ) = \frac { \mathbb { E } [ I _ { i } ] } { \mathbb { E } [ C ] } \frac { \mathbb { E } \left [ z _ { i } ^ { L B _ { i } ^ { ( S _ { i } ) } } \right ] - \mathbb { E } \left [ z _ { i } ^ { L B _ { i } ^ { ( V _ { i } ) } } \right ] } { ( 1 - z ) \left ( \mathbb { E } [ L B _ { i } ^ { ( V _ { i } ) } ] - \mathbb { E } [ L B _ { i } ^ { ( S _ { i } ) } ] \right ) }, \\. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

the left-hand side of (4.6) being E [ I i ] E [ C ] E [ z L ( I i ) i ]. However, substitution of (4.2) and (4.3) in (4.6), and using (3.5), shows that (4.6) is only true if q ( P ) i = 1 for all P ∈ { V , S 1 1 , . . . , V N , S N } , because only in that case all terms in the numerator, except for E [ z LB ( S i ) i ] and E [ z LB ( V i ) i ] cancel out.

Substitution of (4.2), (4.3), and (4.5) in (4.1) gives the desired expression for the PGF of the marginal queue length in Q i .

## Additional remarks

The marginal queue length distribution at an arbitrary epoch is given by (4.1). It is noteworthy that, unlike in Section 3, we can use PASTA now because we focus on all customers - including the impatient ones. This implies that the marginal queue length distribution at arrival and departure epochs of type i customers is also given by (4.1). It should be noted that, when studying the queue length at departure epochs, we assume that reneging customers leave the

system in order of arrival, even though several of them might leave at the same reneging epoch. When looking at it this way, we do not really have group departures, but consecutive departures that might take place during an interval of zero length. Determining the LST of the sojourn time distribution of type i customers, including those that abandon the system before being served, remains difficult. One cannot use the distributional form of Little's law, because there are multiple occasions within a cycle during which a type i customer might leave the system. In the present paper we do not show exactly how to compute the sojourn time distribution of an arbitrary customer, because it requires a lot of bookkeeping. For each customer type, one needs to keep track of when a customer entered the system, and at which point he is going to leave the system. There are N customer types, entering the system during 2 N subperiods, and leaving the system at 2 N + 1 occasions (2 N reneging moments plus one visit period). This gives a maximum of 2(2 N +1) N 2 customer subtypes, although it is determined by the service disciplines in the different queues, which of these subtypes are actually needed. For all customer subtypes, the joint queue lengths at departure moments have to be determined in order to find the sojourn time distributions. We only show how this is done for a vacation model, in Section 5. Although not applicable in its distributional form, Little's law can still be used to determine the mean sojourn time of a type i customer, T i , for i = 1 , . . . , N :

$$\mathbb { E } [ T _ { i } ] = \mathbb { E } [ L _ { i } ] / \lambda _ { i }.$$

## 5 Vacation system with exhaustive service

This section discusses the special case N = 1 and exhaustive service in more detail. The results obtained in the previous sections reduce to nice, compact expressions. When N = 1, there is only one queue being served, and the switch-over time between successive visits is called a server vacation . This model is studied in Adan et al. [1], where it is referred to as the Unique Abandonment Epoch (UAE) model. In [1] only the queue length PGF is determined. The methods used in the present paper also make it possible to find the LSTs of the cycle time and the waiting time distribution. Although it is not discussed explicitly in this section, of course it is also possible to analyse the vacation system with gated service.

We use the same notation as in the rest of the paper, which is slightly different from common notation in vacation models. Since there is only one queue, the indices i = 1 , . . . , N , are dropped. A vacation is the switch-over period S with LST σ ( ), · whereas V denotes the visit period. The analysis in the present section is a slightly more extended version of the one in the previous sections, because we aim at finding the sojourn time of an arbitrary customer, including those that renege, as well. This requires distinguishing not only between moments at which customers abandon the system, but also between their arrival (sub)periods. Therefore, the cycle is divided into four subperiods: V , V a ( S ) b , V ( V ) b and S . During V a the impatient customers abandon the system. All of these customers have arrived during S . The remaining customers, that have also arrived during S , are served during V ( S ) b . During V ( V ) b , all customers that have entered the system during V ( S ) b , and newly arriving customers, are served until the system is empty. Since the system is empty at the beginning of S , we do not split S into subperiods, and we use the notation p for the probability that a customer abandons the system before being served, and q = 1 -p . The astute reader has noticed that we distinguish between three customer types. Type a customers are those that abandon the system during V a , type

b ( S ) customers are those that enter during the vacation S without abandoning the system, type b ( V ) customers enter during the visit period and are always served. The advantage of considering both the arrival epoch and the departure epoch of each customer type, is that it enables us to combine the techniques from Sections 3 and 4. The PGFs of the joint queue length distributions, at the start of the four subperiods, follow from (3.1)-(3.4), and the LST of the cycle time C ∗ follows from (3.6). In this vacation model, it reduces to:

$$\gamma ^ { * } ( \omega ) = \widetilde { L B } ^ { ( V _ { b } ^ { ( S ) } ) } ( 1, \pi ( \omega ) - \frac { \omega } { q \lambda }, 1 ) = \sigma ( \omega + q \lambda ( 1 - \pi ( \omega ) ) ).$$

The mean cycle time and the mean visit time are:

$$\mathbb { E } [ C ] & = \frac { 1 - p \rho } { 1 - \rho } \mathbb { E } [ S ], & \mathbb { E } [ V ] & = \frac { q \rho } { 1 - \rho } \mathbb { E } [ S ],$$

where ρ = λ E [ B ]. In [1], the PGF of the marginal queue length distribution is obtained using similar techniques as in Section 4. In this section we use a different approach, based on the joint queue length distribution at departure epochs. This approach is similar to the one used in Section 3, but now including the customers that renege from the system. We follow the steps taken by Borst [6], who extends an idea of Eisenberg [11], to find the PGF of the joint distribution of the queue lengths and state of the server at departure epochs, M ( P ) ( z , z a ( S ) b , z ( V ) b ). The state of the server is identified by the subperiod P , which can be any of the three periods during which customers depart from the system, V , V a ( S ) b and V ( V ) b .

$$\text{any of the three periods during which customers depart from the system, V} _ { a }, V _ { b } ^ { \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ M ^ { ( V _ { a } ) } ( z _ { a }, z _ { b } ^ { ( S ) }, z _ { b } ^ { ( V ) } ) = & \frac { 1 } { \Lambda \mathbb { E } [ C ] } \frac { 1 } { z _ { a } - 1 } \left ( \widetilde { L B } ^ { ( V _ { a } ) } ( z _ { a }, z _ { b } ^ { ( S ) }, z _ { b } ^ { ( V ) } ) - \widetilde { L B } ^ { ( V _ { b } ^ { ( S ) } ) } ( z _ { a }, z _ { b } ^ { ( S ) }, z _ { b } ^ { ( V ) } ) \right ), \\ M ^ { ( V _ { b } ^ { ( S ) } ) } ( z _ { a }, z _ { b } ^ { ( S ) }, z _ { b } ^ { ( V ) } ) = & \frac { 1 } { \Lambda \mathbb { E } [ C ] } \frac { \beta \left ( \lambda ( 1 - z _ { b } ^ { ( V ) } ) \right ) } { z _ { b } ^ { ( V ) } - \beta \left ( \lambda ( 1 - z _ { b } ^ { ( V ) } ) \right ) } \left ( \widetilde { L B } ^ { ( V _ { b } ^ { ( V ) } ) } ( z _ { a }, z _ { b } ^ { ( S ) }, z _ { b } ^ { ( V ) } ) \right ), \\ \text{The PGF of the joint queue length distribution at an arbitrary departure epoch is simply}$$

The PGF of the joint queue length distribution at an arbitrary departure epoch is simply the sum of these three PGFs. Using PASTA and an up-and-down crossing argument, we find that the marginal queue length distribution at an arbitrary moment is:

$$\mathbb { E } [ z ^ { L } ] = M ^ { ( V _ { a } ) } ( z, z, z ) + M ^ { ( V _ { b } ^ { ( S ) } ) } ( z, z, z ) + M ^ { ( V _ { b } ^ { ( V ) } ) } ( z, z, z ).$$

The sojourn time of an arbitrary customer in a polling model with synchronised reneging has not been discussed in the present paper because of the effort it takes to keep track of all customers and their arrival and departure moments. For this vacation model, it is not too complicated though. It requires applying the generalised distributional form of Little's law, as discussed in [5], to the joint queue length distribution at departure epochs. This leads to

the following LST of the distribution of the sojourn time T of an arbitrary customer:

$$\intertext { n o w o n v e n t o n o f t o n e s o j o u r m e i o f o n a n a r o u r y c u s t o m e r $ } \mathbb { E } [ e ^ { - \omega T } ] & = M ^ { ( \nu _ { \omega } ) } ( 1 - \frac { \omega } { p \lambda }, 1, 1 ) + M ^ { ( \nu _ { b } ^ { ( S ) } ) } ( 1, 1 - \frac { \omega } { q \lambda }, 1 - \frac { \omega } { \lambda } ) + M ^ { ( \nu _ { b } ^ { ( V ) } ) } ( 1, 1 - \frac { \omega } { q \lambda }, 1 - \frac { \omega } { \lambda } ) \\ & = \frac { p ( 1 - \rho ) } { 1 - p \rho } \frac { 1 - \sigma ( \omega ) } { \omega \mathbb { E } [ S ] } + \frac { q ( 1 - \rho ) } { 1 - p \rho } \frac { \sigma ( \omega ) - \sigma ( q \lambda ( 1 - \beta ( \omega ) ) ) } { ( q \lambda ( 1 - \beta ( \omega ) ) - \omega ) \mathbb { E } [ S ] } \beta ( \omega ) \\ & \quad + \frac { 1 - \rho } { 1 - p \rho } \frac { \sigma ( q \lambda ( 1 - \beta ( \omega ) ) ) - 1 } { ( \lambda ( 1 - \beta ( \omega ) ) - \omega ) \mathbb { E } [ S ] } \beta ( \omega ).$$

The waiting time of customers that did not renege the system, is obtained in the same way, but without taking into account the type a customers:

$$\mathbb { E } [ e ^ { - \omega W _ { b } } ] = \frac { \lambda \mathbb { E } [ C ] } { q \lambda \mathbb { E } [ S ] + \lambda \mathbb { E } [ V ] } \left ( M ^ { ( V _ { b } ^ { ( s ) } ) } ( 1, 1 - \frac { \omega } { q \lambda }, 1 - \frac { \omega } { \lambda } ) + M ^ { ( V _ { b } ^ { ( V ) } ) } ( 1, 1 - \frac { \omega } { q \lambda }, 1 - \frac { \omega } { \lambda } ) \right ) \frac { 1 } { \beta ( \omega ) }.$$

## 6 Numerical example

In this section we study the impact of the reneging probabilities on the mean queue lengths in a two-queue polling system. The service times of all customers, in both queues, are exponentially distributed with mean 1. The arrival processes are Poisson with rates 1 10 for Q 1 and 7 10 for Q 2 . We deliberately choose an imbalanced system to study differences between a heavily and a lightly loaded queue. The switch-over times are also exponentially distributed. We compare a system with small switch-over times, E [ S i ] = 1, with a system having larger switch-over times, E [ S i ] = 10. Finally, two different combinations of reneging probabilities are taken:

$$\text{Case } 1 \colon \quad \ p _ { i } ^ { ( V _ { j } ) } = p _ { i } ^ { ( S _ { j } ) } = p _ { i },$$

$$\text{Case 2} \colon \quad \ p _ { i } ^ { ( S _ { i } ) } = \frac { 8 } { 1 0 } p _ { i }, p _ { i } ^ { ( V _ { i + 1 } ) } = \frac { 6 } { 1 0 } p _ { i }, p _ { i } ^ { ( S _ { i + 1 } ) } = \frac { 4 } { 1 0 } p _ { i }, p _ { i } ^ { ( V _ { i } ) } = \frac { 2 } { 1 0 } p _ { i },$$

for i, j = 1 2. , In Case 1, the reneging probabilities per customer type are the same for all reneging moments. In Case 2, which might be considered as more realistic, the probabilities of reneging decrease as the moment of being served comes nearer. The parameters p i are varied, independently for i = 1 2, , between 0 and 1. Furthermore, we study all possible combinations of gated and exhaustive service for each queue. The results are depicted in Figures 1 - 4, where the mean queue lengths E [ L 1 ] and E [ L 2 ] are plotted against p 1 and p 2 . As expected, Q 2 dominates the behaviour of the system, because of its heavy load compared to Q 1 . For this reason, results are omitted for Q 1 receiving exhaustive service, because they hardly deviate from the results where Q 1 receives gated service. A conclusion that can be drawn from a comparison between Figures 2 and 4, is that the lengths of the switch-over times hardly influence the impact of p 1 and p 2 on the mean queue lengths if the reneging probabilities are decreasing as in Case 2. For constant reneging probabilities, as in Case 1, the behaviour of the mean queue lengths changes when the mean switch-over times become larger. The non-monotonic behaviour that was noted in one of the examples studied in [1], is also visible in Figure 3. If switch-over times are relatively large, higher values of p 2 may result in an increase in the mean number of customers in Q 2 , but also in Q 1 if Q 2 receives

gated service. Furthermore, it is interesting to observe in Figures 2 and 4 that in Case 2, both p 1 and p 2 have a high impact on E [ L 1 ] and E [ L 2 ]. In contrast, for Case 1, the influence of p 1 and p 2 varies per queue. E.g., Figures 1 and 3 illustrate that E [ L 2 ] is mainly influenced by p 2 , whereas E [ L 1 ] is influenced by both parameters.

Figure 1: Mean queue lengths in the polling system in Case 1 of Example 2, versus p 1 and p 2 . The reneging probabilities are constant, E [ S i ] = 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0131v1-figure-0af543b4-1384-4670-94d3-05e481e8654e.jpg)

## 7 Conclusions and topics for future research

In the present section we summarise our findings and conclusions, and discuss possible future extensions of the model under consideration. We have extended some results on a vacation model with synchronised reneging, as presented in [1], to a polling system consisting of multiple queues. Using techniques from a polling model with varying arrival rates, depending on the server location (cf. [5]) we have been able to find the LST of the cycle time distribution, visit and intervisit time distributions, and of the waiting time distribution of customers that get served eventually. An adaptation of these techniques leads to the PGF of the marginal queue length distribution of all customers in each queue. It also leads to the sojourn time distribution, but this requires lengthy, cumbersome computations that have only been discussed for a system consisting of one queue.

When comparing the model of the present paper with existing literature, one of the most striking differences is the non-monotonic behaviour of the mean queue lengths that might occur in polling (or vacation) models. As illustrated in the numerical example, the presence of large switch-over times might cause the mean queue lengths to increase if reneging probabilities increase. Another notable difference is that it is relatively easy to compute the

Figure 2: Mean queue lengths in the polling system in Case 2 of Example 2, versus p 1 and p 2 . The reneging probabilities are decreasing, E [ S i ] = 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0131v1-figure-836c7e95-4a64-44c8-b50c-8d6d0422b0b4.jpg)

Figure 3: Mean queue lengths in the polling system in Case 1 of Example 2, versus p 1 and p 2 . The reneging probabilities are constant, E [ S i ] = 10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0131v1-figure-3078771e-e240-4dd3-b1ef-e28af0db7778.jpg)

Figure 4: Mean queue lengths in the polling system in Case 2 of Example 2, versus p 1 and p 2 . The reneging probabilities are decreasing, E [ S i ] = 10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0131v1-figure-de996cd8-4b6c-4882-a383-8644b2501bea.jpg)

proportion of customers served. This metric is a relevant quantity in reneging literature. In our model, given the customer type and server location upon his arrival, one knows exactly how many possibilities this customer will have to abandon the system, and the corresponding reneging probabilities.

Several research topics, beyond the scope of the present paper, are worth studying. An interesting extension is to allow (synchronised) reneging at various epochs during switch-over and visit periods. Possibly, the analysis of the Multiple Abandonment Epochs (MAE) model in [1] (synchronised reneging), and the analysis of Altman and Yechiali [2] (customers growing impatient during vacations) might be extended to polling models. Another extension, relevant from a practical point of view, is to develop numerically more efficient algorithms to compute performance measures of interest. In the present paper, we use the buffer occupancy method, but it would be more efficient to extend the Mean Value Analysis (MVA) framework for polling systems (cf. [27]) to a polling model with reneging at polling instants. In [1], as well as in [5], MVA has been used to find the mean queue lengths. These two implementations should give a good indication of how to implement MVA for the model discussed in the present paper. Note that difficulties in finding the sojourn time distribution do not occur when studying the mean sojourn time, which can simply be found using Little's law.

## Acknowledgements

The author wishes to thank Ivo Adan and Onno Boxma for their valuable comments on earlier drafts of the present paper.

## References

- [1] I. J. B. F. Adan, A. Economou, and S. Kapodistria. Synchronized reneging in queueing systems with vacations. Queueing Systems , 62:1-33, 2009.
- [2] E. Altman and U. Yechiali. Analysis of customers' impatience in queues with server vacations. Queueing Systems , 52(4):261-279, 2006.
- [3] D. Assaf and M. Haviv. Reneging from processor sharing systems and random queues. Mathematics of Operations Research , 15(1):129-138, 1990.
- [4] F. Baccelli and G. Hebuterne. On queues with impatient customers. In F. Kylstra, editor, Performance '81 , pages 159-179. North-Holland, Amsterdam, 1981.
- [5] M. A. A. Boon, A. C. C. van Wijk, I. J. B. F. Adan, and O. J. Boxma. A polling model with smart customers. Eurandom report 2009-038, Eurandom , 2009.
- [6] S. C. Borst. Polling Systems , volume 115 of CWI Tracts . 1996.
- [7] O. J. Boxma. Polling systems. In From universal morphisms to megabytes: A Baayen space odyssey. Liber amicorum for P.C. Baayen , pages 215-230. CWI, Amsterdam, 1994.
- [8] O. J. Boxma and P. R. de Waal. Multiserver queues with impatient customers. In J. Labetoulle and J. W. Roberts, editors, The Fundamental Role of Teletraffic in the Evolution of Telecommunications Networks , Proceedings 14th International Teletraffic Congress (ITC 14) Antibes Juan-les-Pins, France, June 6-10, 1994, pages 743-756. Elsevier, Amsterdam, 1994.
- [9] O. J. Boxma, J. Ivanovs, K. Kosi´ski, and M. Mandjes. L´vy-driven polling systems and n e continuous-state branching processes. Eurandom report 2009-026, Eurandom , 2009.
- [10] J. W. Cohen. The Single Server Queue . North-Holland, Amsterdam, revised edition, 1982.
- [11] M. Eisenberg. Queues with periodic service and changeover time. Operations Research , 20(2):440-451, 1972.
- [12] S. W. Fuhrmann. Performance analysis of a class of cyclic schedules. Technical memorandum 81-59531-1, Bell Laboratories, March 1981.
- [13] 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.
- [14] H. C. Gromoll, P. Robert, A. P. Zwart, and R. F. Bakker. The impact of reneging in processor sharing queues. ACM SIGMETRICS Performance Evaluation Review , 34(1): 87-96, 2006.
- [15] S. E. Haugen, R. B. Queueing systems with stochastic time out. IEEE Transactions on Communications , 28(12):1984-1989, 1980.
- [16] J. Keilson and L. D. Servi. The distributional form of Little's Law and the FuhrmannCooper decomposition. Operations Research Letters , 9(4):239-247, 1990.

- [17] G. Koole and A. Mandelbaum. Queueing models of call centers: An introduction. Annals of Operations Research , 113:41-59, 2002.
- [18] H. Levy and M. Sidi. Polling systems: applications, modeling, and optimization. IEEE Transactions on Communications , 38:1750-1760, 1990.
- [19] K. C. Madan. On a vacation queue with instantaneous customer loss. Microelectronics and Reliability , 36(2):179-188, 1996.
- [20] A. Mandelbaum and S. Zeltyn. The impact of customers patience on delay and abandonment: some empirically-driven experiments with the M/M/n + G queue. OR Spectrum , 26:377-411, 2004.
- [21] C. Palm. Methods of judging the annoyance caused by congestion. Tele. , 4:189-208, 1953.
- [22] J. A. C. Resing. Polling systems and multitype branching processes. Queueing Systems , 13:409 - 426, 1993.
- [23] J. Riordan. Stochastic service systems . Wiley, New York, 1962.
- [24] H. Takagi. Queuing analysis of polling models. ACM Computing Surveys (CSUR) , 20: 5-28, 1988.
- [25] V. M. Vishnevskii and O. V. Semenova. Mathematical methods to study the polling systems. Automation and Remote Control , 67(2):173-220, 2006.
- [26] V. M. Vishnevsky and O. V. Semenova. The power-series algorithm for two-queue polling system with impatient customers. In International Conference on Telecommunications 2008 (ICT 2008) , pages 1-3, 2008.
- [27] E. M. M. Winands, I. J. B. F. Adan, and G.-J. van Houtum. Mean value analysis for polling systems. Queueing Systems , 54:35-44, 2006.
- [28] S. Zeltyn and A. Mandelbaum. Call centers with impatient customers: Many-server asymptotics of the M/M/n + G queue. Queueing Systems , 51:361-402, 2005.
- [29] Y. Zhang, D. Yue, and W. Yue. Analysis of an M/M/ /N 1 queue with balking, reneging and server vacations. Proc. of the 5th International Symposium on Operations Research and its Applications (ISORA'05) , pages 37-47, 2005.
Marko Boon
A Polling Model with Reneging at Polling Instants
In this paper we consider a single-server, cyclic polling system with
switch-over times and Poisson arrivals. The service disciplines that are
discussed, are exhaustive and gated service. The novel contribution of the
present paper is that we consider the reneging of customers at polling
instants. In more detail, whenever the server starts or ends a visit to a
queue, some of the customers waiting in each queue leave the system before
having received service. The probability that a certain customer leaves the
queue, depends on the queue in which the customer is waiting, and on the
location of the server. We show that this system can be analysed by introducing
customer subtypes, depending on their arrival periods, and keeping track of the
moment when they abandon the system. In order to determine waiting time
distributions, we regard the system as a polling model with varying arrival
rates, and apply a generalised version of the distributional form of Little's
law. The marginal queue length distribution can be found by conditioning on the
state of the system (position of the server, and whether it is serving or
switching).