Abstract

This paper presents a discrete-time Geo/G/1 retrial queue with two waiting buffers to model an ATM network, in which the server begins a single vacation in cases where the system is empty at the instant of a service completion. New arriving customer who finds the server being on vacation can decide to either enter the retrial buffer with some probability or leave the system with complementary probability . But the new arriving customer can begin its service immediately if he finds the server idle and join the original buffer if he finds the server busy. We first carry out an extensive analysis of the model by using the supplementary variable method and the generating function approach, and give some performance measures, such as server’s state probabilities and mean queue lengths in the original buffer, retrial buffer, and in the system. Secondly, we give the generating function of the sojourn time of a customer in the system and prove that Little’s law still holds in our model. Sensitivity analysis and cost optimization are finally given for illustrative purposes.

1. Introduction

1.1. Motivation

In this work, we consider a discrete-time retrial queue based on ATM networks. In ATM networks, the basic time unit is a binary code, all information units are transferred through the network in fixed size units, called cells (also called as “customers” in queueing terminology), and each cell contains 48 bytes of user information and 5 bytes of header. Therefore, discrete-time queueing models are more realistic models which include the discrete-time nature and other basic properties of the ATM cell, and readers are referred to Bruneel and Kim [1], Takagi [2], Ramaswami and Wang [3], Walraevens et al. [4], and references therein.

In daily life, two phenomena in ATM networks often occur: one is the retrial phenomenon, which often occurs because of the repeated requests for establishing a communication channel (called as “server”) between nodes and the repeated sending of unsuccessfully received cells; the other is the vacation phenomenon, because sometimes the communication channel in ATM networks needs to be closed for some time (called as “vacation”) due to lack of work, which aims to utilize the idle times for different purposes. For example, computer maintenance and testing and economizing operating cost are being assigned to do a secondary job and so forth. However, the arriving cells are often affected by the server vacation in ATM networks. When the server is on vacation, some arriving cells may temporarily leave the service area and are stored in a retrial buffer (called as orbit), where they can repeat their attempt to connect with the server some time later according to a specific retrial policy, and some nonpersistent arriving cells may leave the system forever and seek services elsewhere. Therefore, in such situation, the retrial phenomenon is incurred by the server’s vacation. Based on the above characteristics of the ATM networks, we mainly focus on investigating the performance measures of such a discrete-time ATM network with vacation and two waiting buffers: original buffer and retrial buffer, in which the arriving cells can begin their services if they find the server idle or are willing to join a waiting queue in the service area (called as an original buffer) if they find the server busy. However, when the system becomes empty, the server goes on a vacation, and some arriving cells either leave the system due to the urgent messages to be transmitted or are likely to enter the orbit to seek service later.

It should be noted that one of the motivations of this paper stems from the interesting mathematical property and its distinction from other retrial queues and vacation queues. The other motivations are its applications in the ATM networks. The advantages and contributions of this work are as follows:(i)Model. A novel ATM network retrial queue is presented. A discrete-time Geo/G/1 retrial queue with server vacation and two waiting buffers is applied in the ATM network. This model is more suitable for reflecting the characteristics of information transmission and performance analysis for ATM networks.(ii)Methodology and Results. We adopt the supplementary variable method to have a Markovian system and use the generating function approach to obtain the performance measures of the system and the sojourn time distribution. Little’s law is proved to be true in our model.(iii)Numerical Illustrations. Sensitivity analysis and cost optimization are included in this work for illustrative purposes.

1.2. Related Works

Since the pioneering works published in 1950s, retrial queues have been extensively studied and used to stochastically model many problems in telecommunication, computer networks, and telephone systems. A good survey of results and fundamental methods in the area of retrial queues can be found in the books, for example, Falin and Templeton [5] and Artalejo and Gomez-Corral [6]. The papers of Artalejo [7, 8] presented a classified bibliography of the work on retrial queues, respectively, in the decades 1990–1999 and 2000–2009. Similarly, Gomez-Corral [9] provided a bibliographical guide to the use of matrix-analytic techniques in retrial queues. Kim and Kim [10] dealt with a survey on various continuous-time retrial queueing models, in which they presented analytic results for queue length distributions, waiting time distributions, and tail asymptotics for the queue length and waiting time distributions and considered the stability analysis of retrial queueing models. Recently, motivated by a contention problem in the downlink direction of wireless base stations in cognitive radio networks, Gao et al. [11] presented a repairable M/G/1 retrial queue with the Bernoulli schedule and a general retrial policy. By using the solution of a Riemann boundary value problem, Dimitriou [12] derived the generating function of the stationary distribution of the number of orbiting customers at service completion epochs, and he also gave the explicit expressions for the expected delay in an orbit without solving a boundary value problem. In the study of the continuous-time retrial queueing models, the Markov process is the main mathematical tool. About more works on the Markov process, readers are referred to Dai et al. [13], Shen et al. [14], and Shen et al. [15].

The study of the retrial queues was mainly focused on the continuous-time. However, it is more realistic to use the discrete-time queues to model telecommunication networks and computer systems.

Over the past decades, Geo/G/1 queues and Geo/G/1 retrial queues with a variety of vacation policies have been well developed, such as Takagi [2], Zhang and Tian [16], Gao and Wang [17], Wang [18], Wang et al. [19], Yue and Zhang [20], Luo et al. [21], and Wang [22, 23] (not exhaustive list). However, most of these vacation queues deal with the case of single waiting queues; i.e., there is only a waiting queue in the service area in the Geo/G/1 queue or an orbit waiting queue in the Geo/G/1 retrial queue. However, it is not the case in daily life. We often see that the waiting queue in the service area and the waiting queue in the orbit may exist simultaneously, for example, in ATM networks. The works on continuous-time retrial queues with two waiting buffers with different characteristics (such as finite or infinite waiting buffer, unreliable server, and constant retrial or general retrial policy) can be found in Krishna Kumar et al. [24], Do [25], Falin [26, 27], Sherman and Kharoufeh [28] and Gao et al. [11], and references therein. However, the discrete-time retrial queues with two waiting buffers have not yet received much attention. It should be remarked that an excellent and complete study on discrete-time queueing system with two waiting buffers was presented by Li and Yang [29], in which an arriving customer either joins the waiting space with probability α in the service area or enters the orbit with probability if he finds the server busy upon arrival. Enlightened by the previous works and the applications of discrete-time retrial queue in the ATM networks, this work is devoted to present the performance analysis of a discrete-time retrial queue with server vacation and two waiting buffers. Remarkably, the models are different in Li and Yang [29] and Gao et al.’ study [11]; on the one hand, we consider the case that the customers arriving during the server busy period always join the waiting buffer in the service area; i.e., these customers are willing to wait in the service area for their services but not to choose to retry after some time later. On the other hand, the customer who arrives during the server vacation period either joins a retrial buffer with probability or leaves the system for ever with probability . To the best of our knowledge, so far, no work has been done on studying the discrete-time Geo/G/1 retrial queue with server vacation, two waiting buffers, and general retrial times by using the supplementary variable method and generating the function approach.

The outline is as follows: In Section 2, the mathematical model is introduced. In Section 3, a detailed analysis of the queueing model is presented. In Section 4, the sojourn time analysis is worked out. In Section 5, sensitivity analysis and cost optimization are given. In Section 6, conclusions are provided.

2. Model Description

We consider a discrete-time Geo/G/1 queue with EAS (early arrival system, see Takagi’s study [2] for details), in which there are two waiting buffers in the system: one is the original buffer in the service area and the other is the retrial buffer (called as orbit), and the server can take a single vacation.

Specifically, let the time slots be numbered by . A potential external customer arrives within the interval just after the beginning of slot t, and a potential service completion occurs within the interval just before the beginning of slot t. We assume that the single vacation policy is taken by the server; that is, if the system becomes empty after a service completion, the server leaves for a vacation with random length . When the vacation ends, the server immediately becomes available until the system becomes empty after a service completion once again. The vacation begins at the slot after a service completion with no customer staying in the system and ends at the slot t. If the arriving external customer finds the server busy or idle, he immediately joins the original buffer under first-come, first-served (FCFS) discipline or begins his service. Otherwise, if the arriving external customer finds the server being on vacation, he enters the orbit with probability or leaves the system forever with probability . The customers in the orbit can access the server according to the FCFS schedule when the original buffer is empty and the server is idle. Various time epochs at which events occur are shown in Figure 1.

Thereinafter, we denote for . The detailed model description is as follows:(1)Independent and identically distributed (i.i.d.) interarrival times with generic random variable A follow the geometric distribution ; i.e., λ is the probability that there is an arrival in each slot and is the probability that no customer arrives in each slot.(2)The service times are i.i.d. random variables with generic random variable B, probability mass function (p.m.f.) , probability generating function (p.g.f.) , and j-th factorial moments . The server can handle only one customer at a time according to the FCFS discipline.(3)The vacation time follows general distribution with p.m.f. , p.g.f. , and j-th factorial moments .(4)Only the retrial customer at the head of the orbit can access to the server if he is idle and the original queue is empty. Successive interretrial times are i.i.d. random variables with generic random variable R; p.m.f. , and p.g.f. .(5)Interarrival times, service times, vacation times, and retrial times are mutually independent.

Remark 1 (special case). When the vacation time takes fixed value 1, the model is changed to the standard discrete-time Geo/G/1 queue.

3. Analysis of the System

For the discrete-time retrial queue considered in this paper, we are chiefly concerned with some performance measures: server’s state probabilities, numbers of customers in the original buffer and retrial buffer, busy cycle, and sojourn time. However, our discrete-time retrial queueing model is not a Markovian queue. Thus, in this section, we introduce some supplementary variables (including the remaining retrial time, the remaining service time, and the remaining vacation time) to make the queueing model Markovian and set up a set of Kolmogorov equations satisfied by the stationary distribution of the system. Based on this, we obtain these important performance measures.

At time , the state of the system can be described by the Markovian chain , , where denotes the state of the server, where denotes the number of customers in the original queue (excluding the one in service); represents the number of customers in the orbit; is a supplementary variable, which denotes the remaining retrial time when and , or the remaining service time when or the remaining vacation time when .

The state space of the Markov chain is , where the state denotes that the server is idle and the number of customers in the orbit is zero; the state indicates that the server is idle, customers are in the orbit, and the remaining retrial time is u; the state represents that the server is busy, there are customers in the original buffer and customers in the retrial buffer, and the remaining service time is u; the state denotes that the server is on vacation, customers are in the orbit, and the remaining vacation time is u.

Our first goal is to find the stationary distribution of the Markov chain under the stability condition :

For the stationary distribution of the system, the Kolmogorov equations are as follows:where if or 1 if .

To obtain the performance measures of the system and analyze the sojourn time of an arbitrary customer, we define the following transforms:

Based on (3)–(7), we will derive the expressions of from simple to complex.

Firstly, multiplying (4) by and and summing over u and , we obtain

Taking yields that

Combining (9) and (10), we have

Secondly, it follows from (7) that

Letting in (12) leads to

Inserting (13) into (12), we have that

Thirdly, from (5) and (6), we have

Putting in (15) leads to

Following from (15) and (16), we obtain that

Now, we derive the expressions of and , respectively, from (11) and (17). Letting in (16), we can get

Combining (16), (17), and (18) yields that

Taking in (11) leads to

Moreover, combining (3), (10), (13), (18), and (20) leads to

Then, from (10), (11), and (21), one can have that

Letting in (13) leads to

Substituting (13) and (23) into (21), we can obtain

Inserting (24) into (19) yields that

From (14), (22), and (25), we should find the probability to obtain the expressions of , and . From (3) and (23), we know that

Using the normalization equation , we obtain thatwhere .

Summing the above results, we have the following Theorem 1.

Theorem 1. If , the joint distribution of the number of customers in the original buffer, the number of customers in orbit, the remaining retrial time, the remaining service time, and the remaining vacation time and the server state, is determined by its transform as follows:

Letting in Theorem 1, we have the following corollary.

Corollary 1. The expressions for the partial generating functions when the server is idle with nonempty orbit, busy, and on vacation, respectively, are as follows:

In the following subsections, we give some performance measures of our queueing system.

3.1. Steady-State Probabilities of the Server

The steady-state probabilities that the server is idle, busy, and on vacation, respectively, denoted by , and , are given by

3.2. Mean Number Customers in Original Buffer and in Orbit

Let be the mean number customers in the original buffer when the server is in state 1 and be the mean number customers in orbit when the server is in state ; then, from Corollary 1, we have that

Let be the mean original buffer size. Since no customers exist in the original buffer when the server is idle and on vacation,

Furthermore, let be the mean number of customers in the orbit:

The mean total customers in the system including the one in service if any, denoted as , can be obtained by using

3.3. Busy Cycle Analysis

Let be the time between two consecutive vacations, which is called a busy cycle, , , and , respectively, be the time periods during which the server is on vacation, idle, and busy, then . Using the alternating renewal theory, we have thatwhich leads to

Furthermore, we obtain that

4. Sojourn Time Distribution

In this section, we consider the distribution of the sojourn time in the system of a test customer (TC). Let T be the TC’s sojourn time with p.g.f. , which begins from the TC’s arrival epoch to the TC’s departure epoch.

To derive the expression of , let D be the busy period in the standard Geo/G/1 queue, which begins at the epoch when one customer arrives and finds the server idle to begin his service immediately and ends at the epoch when one service is completed and the server becomes idle again. From Tian et al. [30], we know that the p.g.f. of D satisfies and .

Theorem 2. If , the expression of the p.g.f. can be given bywhere

Proof. Assume that the TC arrives in , according to the state of the system at time , we have thatwhere (with p.g.f ) denotes the sojourn time of the TC who is in the -th position in the orbit, which begins at the epoch when the server completes his vacation and ends at the epoch when the TC leaves the system.

In the following, we finds the expression of . Assume that the vacation ends at the epoch n, considering whether an external customer arrives or not, we have the following.

Case 1. If there is an external arrival in the time interval (with probability λ), the server will begin a busy period with length D, and then the customer at the head of the orbit can immediately begin retrying to access to the server. Let (with p.g.f. ) be the sojourn time of the TC who is in the -th position calculating from the epoch when the customer at the head of the orbit begins retrying to access to the server to the epoch when the TC leaves the system. In this case, we have .

Case 2. If no external customer arrives in the time interval (with probability ), then the customer at the head of the orbit immediately begins retrying to access to the server. In this case, we have .
Then,However, for , in this case, the TC is in the first position of the orbit, there exists two cases: one is that if the next external customer arrives later than the TC does, then equals the sum of the retrial time and the service time of the TC, i.e., ; the other is that if the next external customer arrives earlier than the TC does, then after a busy period with length D, the TC begins accessing to the server again; thus, , where has the same distribution as and , and are independent random variables. Thus, we obtain thatwhich leads toSimilarly, for conditioning on whether the next external customer arrives earlier or later than the customer at the head of the orbit, we havewhich yields thatIt follows from (41)–(45), and we can obtainSubstituting (46) into (40) leads toBy equation (13) and Theorem 1, the above equation can be changed aswhich completes the proof of Theorem 2.
Using Theorem 2 and , we can obtain the mean sojourn time in the system of an arbitrary TC as follows:After some laborious calculations, we haveInserting the above equations into (49) leads toHere, we can verify thatTherefore, we obtain thatwhich shows that Little’s law holds for our discrete-time Geo/G/1 retrial queue with sever vacation and two waiting buffers.

5. Sensitivity Analysis and Cost Optimization

In this section, we present some numerical examples to show the effect of some parameters on performance measures of the system (i.e., sensitivity analysis) and the optimization problem based on a given cost structure. We assume that the arrival rate , service time follows a discrete phase-type distribution , , , and . The retrial time R follows a discrete phase-type distribution , .

In the following subsection, we first give the sensitivity analysis based on the above results.

5.1. Sensitivity Analysis

Considering different distributions of vacation time V, respectively, with deterministic distribution with p.g.f and negative binomial distribution with p.g.f , we present the effect of D and on , , and for different values of , respectively, in Figures 24 and in Figures 57. Some findings are given as follows.(1)If , takes the fixed value no matter what distribution the vacation time follows, which is because when , all customers who arrive at the system during the vacation period enter the orbit and no customer is lost, and the busy probability is the same as that in the standard Geo/G/1 queue, i.e., .(2)Figure 2 shows that when , the value of has no effect on , , and , i.e., , , and , respectively, take fixed values, which are corresponding to the busy probability, mean sojourn time, and mean busy cycle in the standard Geo/G/1 queue.(3), , and are increasing as increases given that vacation time (i.e., D or ) is fixed and is larger than 1. The reason is that as increases, more customers who arrive during the vacation period enter the orbit, which leads to the increase in busy period and thus makes the server’s busy probability, mean sojourn time of the TC, and mean busy cycle increase.(4)Under given , is decreasing as mean vacation time increases, but and are increasing as mean vacation time increases, which agrees with our expectation.

5.2. Cost Model

In the sequel, we will discuss the optimal vacation length under a given cost assumption because determining an optimal policy to achieve the minimum cost is very important in queueing theory. For the discrete-time Geo/G/1 retrial queue with two waiting buffers and a fixed vacation time, we consider the cost elements as follows:Ch1 ≡ unit time cost of every customer in the original bufferCh2 ≡ unit time cost of every customer in the orbitC0 ≡ unit time cost for keeping the server idleC1 ≡ unit time cost for keeping the server busyC2 ≡ unit time cost for the server being on vacationCs ≡ setup cost per cycle

Firstly, assuming that the vacation time follows deterministic distribution with fixed length rather than the variable, then . Taking D as the decision variable, our objective is to find the optimal value D, say to minimize the cost function , where denotes the long-run average cost per unit time and can be expressed as

The cost minimization problem can be mathematically illustrated as follows:

Owing to the highly nonlinear and complex nature of the optimization problem, it is an arduous task to find the optimal solution analytically. In the following, we present some numerical examples for illustrative purposes and find the optimum solution by Matlab. Taking , we depict the effect of D on the cost function for different values of in Figure 8.

Similarly, in the second numerical example, we consider that the vacation time follows negative binomial distribution with p.g.f and take as the decision variable, and then the cost minimization problem can be mathematically illustrated as follows:where

We present the effect of on the cost function for different values of in Figure 9, where .

From Figures 8 and 9, we can see that the optimal values and minimal values exist and are given in Table 1, which can be obtained by applying the computer software Matlab. Table 1 exhibits that (i) and increase with increasing ; (ii) as increases, the optimal vacation time length decreases; i.e., decreases and increases.

6. Conclusions

In this paper, based on the characteristics of the ATM network, we presented a discrete-time Geo/G/1 retrial queue with two waiting buffers and single vacation, where retrial customers are incurred by the server’s single vacation; that is, arriving customers during the server vacation period either choose to join the retrial buffer or leave the system without service. We studied some important system characteristics by using the supplementary variable method and generating the function approach. Distribution of the sojourn time of any arbitrary customer was studied by using the generating function method. Finally, we presented the sensitivity and cost optimization analysis. As further potential future study, we can generalize this retrial queueing model with two waiting buffers to some different cases, as follows:(i)Retrial customers are incurred by the unavailability of the server, in which an arriving customer who find the server busy or on vacation either joins the waiting buffer with probability or joins the retrial buffer with probability or leaves the system for ever with probability where .(ii)Retrial queues with working vacation policy, where retrial customers are incurred by the lower service rate of the server during the working vacations.(iii)In view of the feedback phenomenon in the retrial queue, one can consider the case that retrial customers are incurred by those feedback customers.(iv)In reality, the capacity of the buffer in service systems is often finite. Hence, one can consider that an arriving customer who finds no vacant space either enters a retrial buffer or leaves the system for ever.

Data Availability

The data used to support the findings of this study are included within the article. The reader can directly access the data in the work.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the Natural Science Foundation of Anhui Higher Education Institutions of China (Grant nos. KJ2017A340 and KJ2019A0535), Fuyang Municipal Government-Fuyang Normal University Horizontal Cooperation projects in 2017 (no. XDHXTD201709), and the National Natural Science Foundation of China (no. 61672006).