Next Article in Journal
Mathematical Model Describing the Hardening and Failure Behaviour of Aluminium Alloys: Application in Metal Shear Cutting Process
Next Article in Special Issue
Renewable k-Out-of-n System with the Component-Wise Strategy of Preventive System Maintenance
Previous Article in Journal
Relation-Theoretic Weak Contractions and Applications
Previous Article in Special Issue
On the Control over the Distribution of Ticks Based on the Extensions of the KISS Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ergodicity and Related Bounds for One Particular Class of Markovian Time—Varying Queues with Heterogeneous Servers and Customer’s Impatience

1
Department of Applied Mathematics, Vologda State University, 15 Lenina Str., 160000 Vologda, Russia
2
Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences, 44-2 Vavilov Str., 119133 Moscow, Russia
3
Department of Applied Probability and Informatics, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Str., 117198 Moscow, Russia
4
Vologda Research Center of the Russian Academy of Sciences, 556A Gorky Str., 160014 Vologda, Russia
5
Moscow Center for Fundamental and Applied Mathematics, Moscow State University, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(9), 1979; https://doi.org/10.3390/math11091979
Submission received: 29 March 2023 / Revised: 18 April 2023 / Accepted: 20 April 2023 / Published: 22 April 2023

Abstract

:
We consider a non-standard class of Markovian time: varying infinite capacity queues with possibly heterogeneous servers and impatience. We assume that during service time, a customer may switch to the faster server (with no delay), when such a server becomes available and no other customers are waiting. As a result, customers in the queue may become impatient and leave it. Under this setting and with certain restrictions on the intensity functions, the quantity of interest, the total number of customers in the system, is the level-dependent birth-and-death process (BPD). In this paper, for the first time in the literature, explicit upper bounds for the distance between two probability distributions of this BDP are obtained. Using the obtained ergodicity bounds in combination with the sensitivity bounds, we assess the stability of BDP under perturbations. Truncation bounds are also given, which allow for numerical solutions with guaranteed truncation errors. Finally, we provide numerical results to support the findings.

1. Introduction

In this paper, consideration is given to certain aspects of the performance evaluation of one particular class of Markovian time: varying queues with heterogeneous servers and customer impatience. Both single- and multiple-class multiserver queues, either with or without impatience, have been the subjects of extensive research for many decades and remain so. Therefore, substantial literature on the topic exists. For the current state of art from various points of view, both theoretical and practical, one can refer, for example, to the introduction sections of [1,2,3,4,5,6,7,8,9,10,11,12,13].
The class of queues, investigated in this paper, consists solely of M t / M t / S / queues with customer impatience. All service times are exponential, and the service intensity of server i is μ i ( t ) . Therefore, servers are allowed to be heterogeneous. Moreover, they are assumed to be numbered as 1 , 2 , , S without repetition by decreasing service speed. To make the latter possible we assume that μ i ( t ) is of the form μ i ( t ) = f ( t ) μ i , where f ( t ) is a locally integrable (non necessarily continuous) bounded positive function and 0 < μ i < . Customers arrive one by one with the intensity λ ( t ) and either occupy one server, if there is one idle, or one place in the queue. The basic quantities under consideration will be X ( t ) : the total number of customer in the system at the instant t, and its distributions. Thus, the queuing discipline is irrelevant as long as it remains work-conserving and non-preemptive (for example, FIFO, LIFO, or RANDOM).
The state space of X ( t ) is { 0 , 1 , 2 } and, from the description given above, it is clear that X ( t ) is the level-dependent BDP. Let us denote its time-dependent intensity matrix (generator) by A ( t ) = ( a i j ( t ) ) i , j = 0 (its form is given in Section 2). The distribution of X ( t ) can be represented as a probability vector p ( t ) = p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , T .
There exist various techniques to calculate performance measures for time-varying queues (see, for example, Introduction in [14,15,16,17]). Probably the most widely applicable (or at least most popular nowadays, due to increasing computer power) set of techniques consists of numerical methods for systems of ordinary differential equations (ODEs). Whenever the queues’ capacities are infinite, the following two questions usually arise, when one attempts to compute stationary (or limiting) performance characteristics of queues using an ODE technique: (i) how to solve the Kolmogorov system d d t p ( t ) = A ( t ) p ( t ) , p ( t ) T 1 = 1 , having infinitely many equations? (ii) how to understand that the limiting regime is reached? Traditionally, (i) is circumvented using truncation. Yet, there is no “rule of thumb” for choosing the truncation threshold, except for “the larger the better”. With respect to (ii), stationary (or limiting) characteristics are usually considered to be identical to the solution of the ODEs at some distant time interval. However, usually it cannot be specified in advance how to choose the position and the length of the “distant time interval”, on which the solution of the system must be found. It can happen that the steady-state is detected prematurely (see [18]).
There are methods in the literature that can solve (i) and (ii) for certain continuous-time Markov chains (see [19,20,21,22,23]). In those cases where the stationary distribution of the Markov chain can be efficiently calculated in advance, solution techniques equipped with the steady-state detection exist (see, for example, [24]). However, in any case, an ODE technique, used to compute the distribution of the queue size or its moments, benefits from the answer to question (ii). This answer follows from the convergence (ergodic) properties of Markov chains. In addition, even though those properties have been the subject of many research papers, it remains difficult to obtain practically useful general ergodicity bounds. In this paper, we show that one well-known method, which is based on the notion of the logarithmic norm of a linear operator and utilizes the properties of linear systems of differential equations, may yield explicit ergodicity bounds for certain time-varying multi-server queues with impatience. Using these bounds, both questions (i) and (ii) can be answered. It must be noticed that the ergodicity results available in the literature for level-(in)dependent nonhomogeneous BDP (and X ( t ) is such a Markov chain) can do the job as well (see, for example, [25,26]). However, being usually quite general, and requiring guessing or estimation of certain quantities, it can prove unclear how to apply them in specific use cases. The novel explicit ergodicity bounds provided in this paper constitute its main contribution. As an additional utility of these bounds, we show that, in combination with the sensitivity bounds, they allow assessment of the stability of X ( t ) under perturbations.
In the description of the class of queues under investigation given above, two ingredients are missed: a server selection rule for the newly arriving customers and the impatience mechanism. When the servers are not identical and an arriving customer sees more than one idle server, they need a rule to select a server. Various selection rules can be found in the literature. Probably the most popular (since [27]) is random selection, according to which arrivals choose randomly any of the idle servers. When arriving customers choose the first idle server, this is the ordered entry rule (see [28,29,30,31]). Assuming the servers are numbered, if a customer occupies the idle server with the lowest number, then this is an ordered hunt rule (see [32]). General selection rules, where arrivals select one of the idle servers with a certain probability, have also been considered (see [33,34]). In this paper, we assume that the ordered hunt rule is adopted. Moreover, in order to avoid the state space collapse, being the main obstacle in the analysis of multi-server heterogeneous queues, we additionally assume that the following service mechanism (known in the literature; see, for example, [10]) is implemented in the system: during the service time, a customer switches to the fastest server (without any delay) when such a server becomes available and there are no other customers waiting. The choice of the impatience mechanism is dictated by the limitations of the method used to obtain the ergodicity bounds. As will be seen, the method heavily relies on the notion of the logarithmic norm of a linear operator, which in the considered case is equal to sup i a i i ( t ) + j i | a j i ( t ) | . For it to be finite, the generator A ( t ) has to be bounded (element-wise) for any t. This condition conflicts with the classic impatience mechanism, according to which the impatience intensity grows as the queue-size increases. Therefore, in this paper we assume that the probability that exactly one customer leaves the queue during the interval ( t , t + Δ ) , given that there are i customers in the system, is ζ i ( t ) Δ + o ( Δ ) , i > S , and ζ S + 1 ( t ) , ζ S + 2 ( t ) , with a single finite upper bound.
In what follows, by · we denote the l 1 -norm, i.e., if x is an ( l + 1 ) -dimensional column vector then x = k = 0 l | x k | . If x is a probability vector, and then x = 1 . The choice of operator norms will be the one induced by the l 1 -norm in column vectors, i.e., A = sup 0 j l i = 0 l | a i j | for a linear operator A.

2. Model Description

Let X ( t ) denote the total number of customers in the model at time t 0 . Then X ( t ) is the continuous-time BDP with the state space { 0 , 1 , 2 } . Denote the birth intensity by λ ( t ) and the death intensity by ν k ( t ) if k customers are in the model at the instant t. We assume that ν k ( t ) can be represented in the form
ν k ( t ) = f ( t ) i = 1 k μ k , if   1 k S , f ( t ) i = 1 S μ k + ζ ( t ) θ k , if   k > S ,
where μ 1 μ S > 0 and { θ k , k S + 1 } is a bounded monotonically non-decreasing sequence of positive numbers i.e., θ k θ k + 1 and there exists some M ( 0 , ) such that θ k M for all k. Another restriction imposed on the intensities λ ( t ) and { ν k ( t ) , k 1 } is that they are locally integrable on [ 0 , ) (not necessarily continuous) functions and bounded in the following sense: there exists some L ( 0 , ) such that sup k 1 | λ ( t ) + ν k ( t ) | L for almost all t 0 .
Since X ( t ) is the BDP, then its transposed generator, further denoted by A ( t ) = ( a i j ( t ) ) i , j = 0 , has the form:
A t = λ ( t ) ν 1 ( t ) 0 0 λ ( t ) λ ( t ) + ν 1 ( t ) ν 2 ( t ) 0 0 λ ( t ) λ ( t ) + ν 2 ( t ) ν 3 ( t ) 0 0 λ ( t ) λ ( t ) + ν 3 ( t ) .
Given an appropriate initial condition, the Kolmogorov forward equations for the distribution p ( t ) = p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , T of X ( t ) can be written as
d d t p ( t ) = A ( t ) p ( t ) .
Since A ( t ) = 2 sup k ( λ ( t ) + ν k ( t ) ) 2 L < for almost all t 0 , the linear operator A ( t ) is bounded and locally integrable for t [ 0 , ) . Therefore (2) is the system of differential equations in the space l 1 with the bounded linear operator and thus it has a unique solution for arbitrary initial conditions.

3. Ergodicity Bounds

All the ergodicity bounds, as shown in this section, are obtained through the application of the same method, known as the “logarithmic norm” method. It based on the notion of the logarithmic norm (see, for example, [14,35]) and utilizes the properties of linear systems of differential equations and thus is directly applicable to (2).

3.1. Null Ergodicity

Fix any ε ( 0 , 1 ) and define the decreasing sequence of positive numbers { d k , k 0 } by d k = ( 1 ε ) k . Let D = d i a g ( d 0 , d 1 , d 2 , ) . Since the inverse D 1 always exists, the system of Equation (2) can be rewritten in the following equivalent form:
d d t D p ( t ) = D A ( t ) D 1 D p ( t ) ,
where
D A ( t ) D 1 = λ ( t ) ν 1 ( t ) 1 ε 0 0 ( 1 ε ) λ ( t ) λ ( t ) + ν 1 ( t ) ν 2 ( t ) 1 ε 0 0 ( 1 ε ) λ ( t ) λ ( t ) + ν 2 ( t ) ν 3 ( t ) 1 ε 0 0 ( 1 ε ) λ ( t ) λ ( t ) + ν 3 ( t ) .
Denote by α k ( t ) the sum of all elements in the kth column of this matrix. By direct inspection, it can be seen that α 0 ( t ) = λ ( t ) ε and α k ( t ) = ε λ ( t ) ν k ( t ) 1 ε for k 1 . Since for a fixed t 0 the sequence of transition rates ν k ( t ) is non-decreasing in k and bounded, then ν ( t ) = lim k ν k ( t ) exists. Therefore, α k ( t ) α k + 1 ( t ) and
α * ( t ) = inf k 0 α k ( t ) = ε λ ( t ) ν ( t ) 1 ε , t 0 .
Theorem 1.
If 0 α * ( t ) d t = + for some ε ( 0 , 1 ) , then the process { X ( t ) , t 0 } is null ergodic, and for any non-negative integer N it holds that
P ( X ( t ) N ) min 1 , e 0 t α * ( u ) d u k = 0 ( 1 ε ) k N p k ( 0 ) , t 0 .
Proof. 
Recall, that D = d i a g ( d 1 , d 2 , ) , where d k = ( 1 ε ) k . From the logarithmic norm method, it follows that D p ( t ) e 0 t α * ( u ) d u D p ( 0 ) for any initial condition p ( 0 ) . Fix any non-negative integer N. Since d k + 1 < d k < 1 for k 1 , then
d N k = 0 N p k ( t ) k = 0 N d k p k ( t ) D p ( t ) .
In the homogeneous case, i.e., when λ ( t ) = λ and ν k ( t ) = ν k for all k 1 , the null ergodicity condition of the Theorem 1 is reduced to λ / ν > 1 , i.e., to the natural condition where the offered load exceeds 1. The right part of (4) takes a simpler form, once the initial probability distribution { p k ( 0 ) , k 0 } has finite support. Under the conditions of the Theorem 1 it is clearly seen that the probability that the process X ( t ) is concentrated in a finite set of states [ a , b ] ( 0 , ) is equal to
P ( a X ( t ) b ) = e 0 t α * ( u ) d u ( 1 ε ) a 1 ( 1 ε ) b ( 1 ε ) a + b 1 k = 0 ( 1 ε ) k p k ( 0 ) ,
and gradually diffuses, i.e., approaches 0 as the time increases indefinitely. Concentration inequalities can be used to construct bounds on the moments of X ( t ) ; for example, the Markov inequality gives the following lower bound for the expected number E ( X ( t ) ) of customers in the model at a time t 0 :
E ( X ( t ) ) N max 0 , 1 e 0 t α * ( u ) d u k = 0 ( 1 ε ) k N + 1 p k ( 0 ) .

3.2. Weak Ergodicity

Put p 0 ( t ) = 1 i 1 p i ( t ) . Consider now the reduced forward Kolmogorov system (2) in the form
d z ( t ) d t = B t z ( t ) + f t , t 0 ,
where f t = λ ( t ) , 0 , 0 , T , z t = p 1 ( t ) , p 2 ( t ) , T and
B t = 2 λ ( t ) + ν 1 ( t ) ν 2 ( t ) λ ( t ) λ ( t ) λ ( t ) λ ( t ) λ ( t ) + ν 2 ( t ) ν 3 ( t ) 0 0 λ ( t ) λ ( t ) + ν 3 ( t ) ν 4 ( t ) .
If instead of the sequence { d k = ( 1 ε ) k , k 0 } used in the previous section, one considers an arbitrary sequence of positive numbers { d k , k 0 } and the triangular (instead of the diagonal) matrix
D = d 0 d 0 d 0 0 d 1 d 1 0 0 d 2 ,
then for d 0 = 1 (as shown in ([25], Section 3)) the sum α k ( t ) of all elements in the kth column of the matrix D A ( t ) D 1 is equal to
α k ( t ) = λ ( t ) + ν k + 1 ( t ) d k + 1 d k λ ( t ) d k 1 d k ν k ( t ) , k 0 , d 1 = 1 .
According to the definition of weak ergodicity (see, for example, [14]), the process { X ( t ) , t 0 } will be weakly ergodic, if p * ( t ) p * * ( t ) 0 as t for any initial conditions p * ( 0 ) and p * * ( 0 ) , where p * ( t ) and p * * ( t ) are the corresponding solutions of (2). As shown in ([25] Theorem 3.1), the quantity p * ( t ) p * * ( t ) for any regular birth and death process (and thus for X ( t ) as well) can be upper bounded by applying the logarithmic norm method; specifically, the relation (3.11) gives the explicit form of the bound:
p * ( t ) p * * ( t ) e 0 t α * ( u ) d u i = 1 4 k = 0 i 1 d k inf k 1 d k p i * ( 0 ) p i * * ( 0 ) ,
where α * ( t ) = inf k 0 α k ( t ) . The next theorem gives the sequence { d k , k 0 } and the sufficient condition, under which the right part of the inequality (7) tends to 0 as time increases indefinitely, and therefore X ( t ) is weakly ergodic.
Theorem 2.
If 0 ν S ( t ) d λ ( t ) d t = + for a d ( 1 , 1 + 1 + μ S / i = 1 S 1 μ i ] , then the process { X ( t ) , t 0 } is weakly ergodic and for any initial conditions p * ( 0 ) and p * * ( 0 ) it holds that
p * ( t ) p * * ( t ) 4 e d 1 d 0 t ν S ( u ) d λ ( u ) d u i = 1 1 + d i 1 1 d 1 p i * ( 0 ) p i * * ( 0 ) , t 0 .
Proof. 
Assume that 0 ν S ( t ) d λ ( t ) d t = + and d ( 1 , 1 + 1 + μ S / i = 1 S 1 μ i ] . Define the increasing sequence of positive numbers { d k , k 0 } by d 0 = 1 , d k = d k 1 for k 1 . Then the expression for α k ( t ) given by (6) is reduced to α k ( t ) = ν k + 1 ( t ) d 1 ν k ( t ) ( d 1 ) λ ( t ) . Now, depending on whether k < S or k S , the values of α k ( t ) can be lower bounded. Indeed, from the exact form of ν k ( t ) , given by (1), and the assumption θ k θ k + 1 , it follows that
α k ( t ) 1 d 1 i = 1 k μ i + μ k + 1 f ( t ) d 1 λ ( t ) , if   0 k S 1 , 1 d 1 i = 1 S μ i f ( t ) d 1 λ ( t ) , if   k S .
Consider the following two functions on ( 1 , ) :
ζ 1 ( x ) = min μ 1 , μ 2 + ( 1 x 1 ) i = 1 1 μ i , , μ S + ( 1 x 1 ) i = 1 S 1 μ i ,
ζ 2 ( x ) = 1 x 1 i = 1 S μ i .
Both are increasing and bounded, that is μ S < ζ 1 ( x ) < μ 1 and 0 < ζ 2 ( x ) < i = 1 S μ i . Moreover ζ 2 ( x ) < μ S as long as x ( 1 , 1 + μ S / i = 1 S 1 μ i ] . Therefore, as long as x varies in this region, we have ζ 1 ( x ) > ζ 2 ( x ) and, therefore, inf k 0 α k ( t ) α S ( t ) . Now the statement of the theorem follows from (7). □
Corollary 1.
Under the conditions of Theorem 2 the process X ( t ) has the limiting mean and if inf i 1 d i i > 0 then
| i = 0 i p i ( t ) m * ( t ) | 1 inf i 1 d i 1 i D ( p t p * t ) .
where p * ( t ) = p 0 * ( t ) , p 1 * ( t ) , p 2 * ( t ) , T is the limiting distribution of X ( t ) .
Proof. 
Since in Theorem 2 d i = d i 1 , i 1 and thus D p * t < for any t 0 we can see that the limiting mean
i = 1 i p i * ( t ) = i = 1 d i 1 d i i p i * ( t ) 1 inf i 1 d i i D p * t
exists. Now for any fixed X ( 0 ) we find that
| i = 0 i ( p i ( t ) m ( t ) ) | 1 inf i 1 d i 1 i i = 1 d i | p i ( t ) p i * ( t ) |
which is identical to the right part of (9). □
Note that if all the intensity functions are periodic and the smallest common multiple of the periods is T, then the sufficient weak ergodicity condition of Theorem 2 is reduced to 0 T ν S ( t ) d λ ( t ) d t > 0 , and from the Corollary 1 it follows that the limiting mean will have a period equal to T.
Corollary 2.
Assume λ ( t ) = λ and ν k ( t ) = ν k for all k 1 . If λ < ν , then X ( t ) is exponentially ergodic. If π = ( π 0 , π 1 , ) T is its stationary distribution, then
p ( t ) π 4 e α t i = 1 k = 0 i 1 d k p i ( 0 ) π i , t 0 ,
where the sequence { d k , k 0 } and the parameter α > 0 can be given explicitly.
Proof. 
Put d 0 = 1 and d k = i = 1 k ( 1 + ε i 1 ) , k 1 , for some positive sequence of real numbers ε 0 , ε 1 , . Then from (6) it follows that
α k ν k ε k 1 1 + ε k 1 λ ε k .
Let λ < ν . Then, there must exist a positive ε such that ( 1 + ε ) λ < ν or, equivalently, ν ε 1 + ε λ ε > 0 . Since ν k is a non-decreasing sequence, then one can find the smallest k * for which simultaneously
ν k * + 1 ε 1 + ε λ ε > 0   and   ν k * ε 1 + ε λ ε 0 .
Put ε k = ε for k k * . Then
inf k k * + 1 α k ν k * + 1 ε 1 + ε λ ε = α .
It can be shown that if ε and ε k * 1 , , ε 1 are chosen consecutively from the following intervals:
ε 0 , i = 1 k * ν i i = 1 k * λ i j = 1 k * i ν i , ε k λ ε k + 1 ν k + 1 λ ε k + 1 , i = 1 k ν i i = 1 k λ i j = 1 k i ν i ,
then it is guaranteed that α k α for all 1 k k * 1 . Hence, inf k 1 α k α , and thus (10) follows from (7). □

4. Perturbation and Truncation Bounds

Since X ( t ) is the BDP, the results regarding perturbation and truncation bounds, obtained for quite a general setting in [36,37,38]), are applicable.
We start with the perturbation bounds. Consider a Markov chain, say X ¯ ( t ) , with a state space identical to X ( t ) and with such a generator, say A ¯ ( t ) = a ¯ i j ( t ) i , j = 0 , that A ¯ ( t ) = A ( t ) + A ^ ( t ) holds to all t 0 . Here A ^ ( t ) = a ^ i j ( t ) i , j = 0 is such a matrix that is small or A ^ ( t ) 0 as t increases. Denote the state probabilities of the “perturbed” chain by p ¯ i ( t ) . Using the normalization condition p ¯ 0 ( t ) = 1 i 1 p ¯ i ( t ) , the system corresponding X ¯ ( t ) can be written as follows:
d d t z ¯ ( t ) = B ¯ ( t ) z ¯ ( t ) + f ¯ ( t ) .
Here the matrix B ¯ ( t ) with the elements b ¯ i j ( t ) = a ¯ i j ( t ) a ¯ i 0 ( t ) . Let the perturbed intensities λ ¯ ( t ) and ν ¯ k ( t ) be such that | λ ( t ) λ ¯ ( t ) | ϵ ^ , | ν k ( t ) ν ¯ k ( t ) | ϵ ^ for t 0 . Then for all t 0 we have the upper bounds:
B ( t ) 1 D ( 1 + d ) | λ ( t ) | + 1 + d d | ν k ( t ) | ( 1 + d ) L ,
f ( t ) 1 D = | λ ( t ) | Λ ,
B ( t ) B ¯ ( t ) 1 D ( 1 + d ) | λ ( t ) λ ¯ ( t ) | + 1 + d d | ν k ( t ) ν ¯ k ( t ) | ( 1 + d ) 2 d ϵ ^ ,
f ( t ) f ¯ ( t ) 1 D | λ ( t ) λ ¯ ( t ) | ϵ ^ .
In the article in [37], the following theorem was proven.
Theorem 3.
If a Markov chain X ( t ) is 1 D -exponentially weakly ergodic
e d 1 d 0 t ν S ( u ) d λ ( u ) d u M e a ( t s )
for all 0 s t , then X ¯ ( t ) is also 1 D -exponentially weakly ergodic and the following perturbation estimate in the 1 D -norm holds:
lim sup t p ( t ) p ¯ ( t ) 1 D M ϵ ^ M ( 1 + d ) 2 Λ + a d a a d M ( 1 + d ) 2 ϵ ^ .
If W = inf i 1 d i i > 0 , then both chains X ( t ) and X ¯ ( t ) have limiting means and
lim sup t | ϕ ( t ) ϕ ¯ ( t ) | M ϵ ^ M ( 1 + d ) 2 Λ + a d W a a d M ( 1 + d ) 2 ϵ ^ .
Assume that the conditions of the theorem are fulfilled, i.e.,
B ( t ) B ¯ ( t ) 1 D χ ( t ) , f ( t ) f ¯ ( t ) 1 D χ ( t ) ,
where χ ( t ) 0 as t . Then the bound from [36]
z ^ ( t ) M e a t z ^ ( 0 ) + M e a ( t t * ) ε 0 a + ε a 1 + M z ¯ ( 0 ) + M ( F + ε 0 ) a ε 0 ,
where
f ( t ) 1 D F , χ ( 0 ) = ε 0 , ε ( 0 , ε 0 ) , χ ( t * ) = ε ,
yields the perturbation bound for the considered X ( t ) .
Let us consider the truncation bound. Consider the copy of the process X ( t ) but with the finite state space { 0 , 1 , 2 , N } (denote it by X * ( t ) ). Let its birth intensity λ * ( t ) = λ ( t ) and death intensity by ν k * ( t ) = ν k ( t ) , when there are, in total, k customers in the model. Denote its transposed generator by A * ( t ) and its probability vector by p * ( t ) . For the Markov chain X * ( t ) , instead of (2), we have
d d t p * ( t ) = A * t p * t .
Assume that for t 0
λ ( t ) Λ L < , ν n ( t ) Δ n L < .
In order to obtain a truncation bound, one needs to find proper sequences { d k } and { d k * } . Assume that for the sequence { d k } there exist positive numbers M and α , such that
e s t α ( τ ) d τ M e α ( t s ) ,
for all 0 s t , where α t = inf k 1 α k t and
α k t = λ t + ν 1 t d 2 d 1 λ t t , k = 1 , λ t + ν k + 1 t d k + 2 d k + 1 λ t d k d k + 1 ν k t , k > 1 .
Assume that for the sequence { d k * } there exist positive numbers M * and α * , such that
e s t α * ( τ ) d τ M * e α * ( t s )
for all 0 s t , where α * t = min α k * t and
α k * t = λ t + ν 1 t d 2 * d 1 * λ t t , k = 1 , λ t + ν k + 1 t d k + 2 * d k + 1 * λ t d k * d k + 1 * ν k t , k > 1 .
Denote
G k = j = 1 k d j , G k * = j = 1 k d j * , d = d 1 , W = inf k > 0 G k k .
The following theorem holds.
Theorem 4.
Assume the Markov chains X ( t ) and X * ( t ) satisfy (20) and (21) and X ( 0 ) = X * ( 0 ) = 0 . Then the following bounds of truncations hold:
p t p * t 4 M M * Λ d i + 1 * d α α * G N + 1 Λ G N * ,
p t p * t 1 E 4 M M * Λ d i + 1 * W α α * G N + 1 Λ G N * .

5. Examples

In this section, we will consider three different examples, demonstrating the applicability of the obtained results. The values of the parameters used do not follow any real data but were chosen for illustration purposes only. As Example 1, let us consider the homogeneous system with three servers (i.e., S = 3 ), periodic service intensity μ 1 ( t ) = μ 2 ( t ) = μ 3 ( t ) = μ ( t ) = 2 + cos 2 π t , and periodic arrival intensity λ ( t ) = 3 + sin 2 π t . Assume that the impatience intensity is time-independent and has the form ζ k ( t ) = 0.4 ( 1 k 1 ) , k > S . Thus, ζ = 0.1 and θ k = 4 ( 1 k 1 ) . Then from (6) it follows that α k ( t ) = μ ( t ) + 1 d 1 k μ ( t ) d 1 λ ( t ) . The conditions of Theorem 2 are satisfied if d = 3 / 2 and the function α * ( t ) in (7) is equal to
α * ( t ) = α 0 ( t ) = α 3 ( t ) = μ ( t ) 0.5 λ ( t ) = 0.5 + cos 2 π t 0.5 sin 2 π t .
In Figure 1, it is shown how the probabilities p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , p 3 ( t ) , and p 4 ( t ) behave as t increases, if initially the system was empty. In this case, the functions quite quickly approach their limiting values. Convergence to the limiting value of the empty system probability p 0 ( t ) is shown in Figure 2. It can be seen that starting from approximately t = 40 = t * the system “forgets” its initial state and the distribution of X ( t ) for t > t * can be regarded as limiting. Moreover, since the limiting distribution of X ( t ) is periodic, it is sufficient to solve (numerically) the system of ODEs only in the interval [ 0 , t * + T ] , where T is the smallest common multiple of the periods of λ ( t ) and μ ( t ) i.e., T = 1 . The probability distribution of X ( t ) in the interval [ t * , t * + T ] is the estimate of the limiting probability distribution of X ( t ) . Convergence of the conditional mean number of customers in the system is shown in Figure 3. Starting from the instant t = t * , the system “forgets” its initial state and the value of i = 0 i p i ( t ) can be regarded as the limiting value of the mean number of customers.
If one increases (some of) the service intensities, then the rate of convergence is not expected to slow. This is illustrated by the Example 2 (see Figure 4, Figure 5 and Figure 6), in which the service intensities are changed to
μ 1 ( t ) = 1.5 2 + cos 2 π t , μ 2 ( t ) = 2 + cos 2 π t , μ 3 ( t ) = 0.5 2 + cos 2 π t .
The conditions of the Theorem 2 are still satisfied if d = 3 / 2 , but the function α * ( t ) in (7) is now equal to
α * ( t ) = 0.5 + cos 2 π t 0.5 sin 2 π t .
As the last example, let us consider a heterogeneous system with three servers, periodic arrival intensity λ ( t ) = 3 + sin 2 π t , and periodic service intensities (24). Let k m a x be a positive integer, η > 0 and q ( 0 , 1 ] . If we put ζ = η q k m a x and
θ k = k S k m a x , k = S + 1 , S + 2 , . . . , S + k m a x , 1 , k > S + k m a x
then, as long as k < S + k m a x , the time-independent impatience intensity ζ k ( t ) = ζ θ k is equal to ( k η ) q . In other words, the adopted impatience mechanism mimics the classic one (with impatience intensity k η and probability of successful impatience q), as long as the queue size stays below S + k m a x . Figure 6 and Figure 7 show the convergence and the limiting value of the empty system probability p 0 ( t ) for q = 0.2 and various values of k m a x = 2 , 3 , 10 , 200 . The limiting values of the conditional mean are shown in the Figure 8 and Figure 9. For k m a x > 10 the results are almost indistinguishable from the case k m a x = 10 . The rate of convergence here is identical to in Example 2, and for d = 1.9 we have
α * ( t ) = 9 19 0.3 + 3 cos 2 π t 1.9 sin 2 π t .
By fixing α = 0.142 and M = 3 , we obtain the following bound on the rate of convergence:
p * ( t ) p * * ( t ) 12 e 0.142 t i = 1 1 + 1 . 9 i 1 1 0.9 p i * ( 0 ) p i * * ( 0 ) , t 0 .
The perturbation bound
lim sup t p ( t ) p ¯ ( t ) 1 D 1126 ϵ ^ 0.142 13.3 ϵ ^ .
and thus, if ϵ ^ = 0.000001 , we have that lim sup t p ( t ) p ¯ ( t ) 1 D 0.02 .
In order to obtain a truncation bound, let us consider the sequence { d * } with d = 100 99 . Then
α * ( t ) = 49 1650 + 3 100 cos ( 2 π t ) 1 99 sin ( 2 π t )
By putting α * = 49 1650 , M * = exp ( 3 / 100 + 1 / 99 π ) , Λ = 4 , we have that G 200 * 6 × 10 55 and G 200 639 . Therefore
p t p * t 4 × 3 × 2 × 4 0.14 ( 49 / 1650 ) 639 × 4 6 × 10 55 10 50 .

6. Conclusions

It is well known that numerical solution techniques for the computation of time-dependent probialbity distributions of Markovian systems can benefit from methods providing ergodicity bounds, because the latter can indicate how to choose the position and the length of the “distant time interval” (in the periodic case) on which the solution has to be computed. They are also helpful whenever state space truncation is required, since they allow for guaranteed truncation error. In this paper, explicit ergodicity, perturbation, and truncation bounds for a particular class of Markovian time-varying queue with heterogeneous servers and customer impatience were obtained. Further research in this area could concentrate on the relaxation of the assumptions made about the transition intensities. For example, for the presented solution, the monotonicity of θ k is crucial (it guarantees the existence of the lim k ν k ( t ) ). If, in addition, the transition intensities are time-independent, the condition λ lim k ν k gives exponential/weak/null ergodicity. If the monotonicity assumption is dropped, other justifications need to be found.

Author Contributions

Conceptualization, supervision, A.Z.; methodology, A.Z. and Y.S.; investigation, software, validation, visualization, Y.S. and I.K.; writing, R.R. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Ministry of Science and Higher Education of the Russian Federation, project No. 075-15-2020-799.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Haight, F.A. Queueing with reneging. Metrika 1959, 2, 186–197. [Google Scholar] [CrossRef]
  2. Hasenbein, J.; Perry, D. Introduction: Queueing systems special issue on queueing systems with abandonments. Queueing Syst. 2013, 75, 111–113. [Google Scholar] [CrossRef]
  3. Van Ommeren, J.K.; Baer, N.; Mishra, N.; Roy, D. Batch service systems with heterogeneous servers. Queueing Syst. 2020, 95, 251–269. [Google Scholar] [CrossRef]
  4. Klimenok, V.; Dudin, A.; Vishnevsky, V. Priority Multi-Server Queueing System with Heterogeneous Customers. Mathematics 2020, 8, 1501. [Google Scholar] [CrossRef]
  5. Harchol-Balter, M.; Osogami, T.; Scheller-Wolf, A.; Wierman, A. Multi-Server Queueing Systems with Multiple Priority Classes. Queueing Syst. 2005, 51, 331–360. [Google Scholar] [CrossRef]
  6. Jouini, O.; Roubos, A. On multiple priority multi-server queues with impatience. J. Oper. Res. Soc. 2014, 65, 616–632. [Google Scholar] [CrossRef]
  7. Puha, A.L.; Ward, A.R. Scheduling an Overloaded Multiclass Many-Server Queue with Impatient Customers. In Operations Research & Management Science in the Age of Analytics; INFORMS TutORials: New York, NY, USA, 2019; pp. 189–217. [Google Scholar] [CrossRef]
  8. Dong, J.; Ibrahim, R. SRPT Scheduling Discipline in Many-Server Queues with Impatient Customers. Manag. Sci. 2021, 67, 7708–7718. [Google Scholar] [CrossRef]
  9. Efrosinin, D.; Stepanova, N.; Sztrik, J. Algorithmic Analysis of Finite-Source Multi-Server Heterogeneous Queueing Systems. Mathematics 2021, 9, 2624. [Google Scholar] [CrossRef]
  10. Kumar, R.; Sharma, S. Transient analysis of a Markovian queuing model with multiple-heterogeneous servers, and customers’ impatience. Opsearch 2021, 58, 540–556. [Google Scholar] [CrossRef]
  11. Melikov, A.Z.; Ponomarenko, L.A.; Mekhbaliyeva, E.V. Analyzing the Models of Systems with Heterogeneous Servers. Cybern. Syst. Anal. 2020, 56, 89–99. [Google Scholar] [CrossRef]
  12. Dudin, A.; Dudina, O.; Dudin, S.; Samouylov, K. Analysis of Multi-Server Queue with Self-Sustained Servers. Mathematics 2021, 9, 2134. [Google Scholar] [CrossRef]
  13. Bhati, A.; Pillai, S.R.B.; Vaze, R. On the Age of Information of a Queuing System with Heterogeneous Servers. In Proceedings of the National Conference on Communications (NCC), Kanpur, India, 27–30 July 2021; pp. 1–6. [Google Scholar] [CrossRef]
  14. Zeifman, A.; Satin, Y.; Kovalev, I.; Razumchik, R.; Korolev, V. Facilitating Numerical Solutions of Inhomogeneous Continuous Time Markov Chains Using Ergodicity Bounds Obtained with Logarithmic Norm Method. Mathematics 2021, 9, 42. [Google Scholar] [CrossRef]
  15. Arns, M.; Buchholz, P.; Panchenko, A. On the numerical analysis of inhomogeneous continuous-time Markov chains. Informs J. Comput. 2010, 22, 416–432. [Google Scholar] [CrossRef]
  16. Schwarz, J.A.; Selinka, G.; Stolletz, R. Performance analysis of time-dependent queueing systems: Survey and classification. Omega 2016, 63, 170–189. [Google Scholar] [CrossRef]
  17. Sidje, R.B.; Burrage, K.; Macnamara, S. Inexact uniformization method for computing transient distributions of Markov chains. Siam J. Sci. Comput. 2007, 29, 2562–2580. [Google Scholar] [CrossRef]
  18. Zapreev, I.S.; Katoen, J.-P. Safe on-the-fly steady-state detection for time-bounded reachability. In Proceedings of the 3rd International Conference on the Quantitative Evaluation of Systems, Riverside, CA, USA, 11–14 September 2006. [Google Scholar] [CrossRef]
  19. Down, D.; Meyn, S.P.; Tweedie, R.L. Exponential and Uniform Ergodicity of Markov Processes. Ann. Probab. 1995, 23, 1671–1691. [Google Scholar] [CrossRef]
  20. Meyn, S.P.; Tweedie, R.L. Markov Chains and Stochastic Stability; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  21. Zeifman, A.; Satin, Y.; Kiseleva, K.; Korolev, V. On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process. Mathematics 2019, 7, 477. [Google Scholar] [CrossRef]
  22. Masuyama, H. Error bounds for augmented truncations of discrete-time block-monotone Markov chains under geometric drift conditions. Adv. Appl. Probab. 2015, 47, 83–105. [Google Scholar] [CrossRef]
  23. Tweedie, R.L. Truncation approximations of invariant measures for Markov chains. J. Appl. Probab. 1998, 35, 517–536. [Google Scholar] [CrossRef]
  24. Burak, M.R.; Korytkowski, P. Inhomogeneous CTMC Birth-and-Death Models Solved by Uniformization with Steady-State Detection. Acm Trans. Model. Comput. Simul. (Tomacs) 2020, 30, 1–18. [Google Scholar] [CrossRef]
  25. Zeifman, A.I. Upper and lower bounds on the rate of convergence for nonhomogeneous birth and death processes. Stoch. Process. Their Appl. 1995, 59, 157–173. [Google Scholar] [CrossRef]
  26. Granovsky, B.L.; Zeifman, A. Nonstationary Queues: Estimation of the Rate of Convergence. Queueing Syst. 2004, 46, 363–388. [Google Scholar] [CrossRef]
  27. Gumbel, H. Waiting Lines with Heterogeneous Servers. Oper. Res. 1960, 8, 504–511. [Google Scholar] [CrossRef]
  28. Cooper, R.B. Queues with ordered servers that work at different rates. Opsearch 1976, 13, 69–78. [Google Scholar]
  29. Nawijn, W.M. A note on many-server queueing systems with ordered entry with applications to conveyor theore. J. Appl. Probab. 1983, 20, 144–152. [Google Scholar] [CrossRef]
  30. Razumchik, R.; Zaryadov, I. Stationary Blocking Probability in Multi-server Finite Queuing System with Ordered Entry and Poisson Arrivals. In Distributed Computer and Communication Networks, Proceedings of the DCCN 2015—Communications in Computer and Information Science, Moscow, Russia, 19–22 October 2016; Vishnevsky, V., Kozyrev, D., Eds.; Springer: Cham, Switzerland, 2016; Volume 601, pp. 344–357. [Google Scholar] [CrossRef]
  31. Meykhanadzhyan, L.; Matyushenko, S.; Pyatkina, D.; Razumchik, R. Revisiting joint stationary distribution in two finite capacity queues operating in parallel. Inform. Primen. 2017, 11, 106–112. [Google Scholar] [CrossRef]
  32. Baxley, R.V.N. The Multiple-Server Queue with Heterogeneous Service Times. Ph.D. Thesis, Georgia Institute of Technology, Atlanta, GA, USA, 1973. [Google Scholar]
  33. Yu, O.S.; Neuts, M.F. The steady state solution of a heterogeneous-server queue with Erlang service times. Tims Stud. Manag. Sci. 1977, 7, 199–213. [Google Scholar]
  34. Grassmann, W.K.; Zhao, Y.Q. Heterogeneous Multiserver Queues with General Input. Infor Inf. Syst. Oper. Res. 1997, 35, 208–224. [Google Scholar] [CrossRef]
  35. Zeifman, A.; Razumchik, R.; Satin, Y.; Kovalev, I. Ergodicity Bounds for the Markovian Queue with Time-Varying Transition Intensities, Batch Arrivals and One Queue Skipping Policy. Appl. Math. Comput. 2021, 395, 125846. [Google Scholar] [CrossRef]
  36. Zeifman, A.I.; Korolev, V.Y.; Razumchik, R.V.; Satin, Y.A.; Kovalev, I.A. Limiting Characteristics of Queueing Systems with Vanishing Perturbations. Dokl. Math. 2022, 106, 375–379. [Google Scholar] [CrossRef]
  37. Zeifman, A.; Korolev, V.; Satin, Y. Two approaches to the construction of perturbation bounds for continuous-time Markov chains. Mathematics 2020, 8, 253. [Google Scholar] [CrossRef]
  38. Satin, Y.A.; Razumchik, R.V.; Zeifman, A.I.; Kovalev, I.A. Upper bound on the rate of convergence and truncation bound for non-homogeneous birth and death processes on Z. Appl. Math. Comput. 2022, 423, 127009. [Google Scholar] [CrossRef]
Figure 1. Example 1. Behavior of the probabilities p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , p 3 ( t ) , and p 4 ( t ) given that initially the system was empty (i.e., p 0 ( 0 ) = 1 ).
Figure 1. Example 1. Behavior of the probabilities p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , p 3 ( t ) , and p 4 ( t ) given that initially the system was empty (i.e., p 0 ( 0 ) = 1 ).
Mathematics 11 01979 g001
Figure 2. Example 1. Convergence and the limiting value of the empty system probability p 0 ( t ) given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Figure 2. Example 1. Convergence and the limiting value of the empty system probability p 0 ( t ) given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Mathematics 11 01979 g002
Figure 3. Example 1. Convergence of the conditional mean i = 0 i p i ( t ) number of customers in the system given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Figure 3. Example 1. Convergence of the conditional mean i = 0 i p i ( t ) number of customers in the system given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Mathematics 11 01979 g003
Figure 4. Example 2. Behavior of the probabilities p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , p 3 ( t ) and p 4 ( t ) , given that initially the system was empty (i.e., p 0 ( 0 ) = 1 ).
Figure 4. Example 2. Behavior of the probabilities p 0 ( t ) , p 1 ( t ) , p 2 ( t ) , p 3 ( t ) and p 4 ( t ) , given that initially the system was empty (i.e., p 0 ( 0 ) = 1 ).
Mathematics 11 01979 g004
Figure 5. Example 2. Convergence and the limiting value of the empty system probability p 0 ( t ) , given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Figure 5. Example 2. Convergence and the limiting value of the empty system probability p 0 ( t ) , given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Mathematics 11 01979 g005
Figure 6. Example 2. Convergence of the conditional mean i = 0 i p i ( t ) number of customers in the system, given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Figure 6. Example 2. Convergence of the conditional mean i = 0 i p i ( t ) number of customers in the system, given two different initial conditions: p 0 ( 0 ) = 1 (blue line), p 99 ( 0 ) = 1 (orange line).
Mathematics 11 01979 g006
Figure 7. Example 3. The perturbation bounds (green lines) for the limiting values of p 0 ( t ) for k = 200 .
Figure 7. Example 3. The perturbation bounds (green lines) for the limiting values of p 0 ( t ) for k = 200 .
Mathematics 11 01979 g007
Figure 8. Example 3. Behavior of the limiting values of p 0 ( t ) , given that initially the system was empty (i.e., p 0 ( 0 ) = 1 ) for various values of k m a x .
Figure 8. Example 3. Behavior of the limiting values of p 0 ( t ) , given that initially the system was empty (i.e., p 0 ( 0 ) = 1 ) for various values of k m a x .
Mathematics 11 01979 g008
Figure 9. Example 3. Behavior of the limiting conditional mean i = 0 i p i ( t ) number of customers in the system for various values of k m a x .
Figure 9. Example 3. Behavior of the limiting conditional mean i = 0 i p i ( t ) number of customers in the system for various values of k m a x .
Mathematics 11 01979 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Satin, Y.; Razumchik, R.; Kovalev, I.; Zeifman, A. Ergodicity and Related Bounds for One Particular Class of Markovian Time—Varying Queues with Heterogeneous Servers and Customer’s Impatience. Mathematics 2023, 11, 1979. https://doi.org/10.3390/math11091979

AMA Style

Satin Y, Razumchik R, Kovalev I, Zeifman A. Ergodicity and Related Bounds for One Particular Class of Markovian Time—Varying Queues with Heterogeneous Servers and Customer’s Impatience. Mathematics. 2023; 11(9):1979. https://doi.org/10.3390/math11091979

Chicago/Turabian Style

Satin, Yacov, Rostislav Razumchik, Ivan Kovalev, and Alexander Zeifman. 2023. "Ergodicity and Related Bounds for One Particular Class of Markovian Time—Varying Queues with Heterogeneous Servers and Customer’s Impatience" Mathematics 11, no. 9: 1979. https://doi.org/10.3390/math11091979

APA Style

Satin, Y., Razumchik, R., Kovalev, I., & Zeifman, A. (2023). Ergodicity and Related Bounds for One Particular Class of Markovian Time—Varying Queues with Heterogeneous Servers and Customer’s Impatience. Mathematics, 11(9), 1979. https://doi.org/10.3390/math11091979

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop