1. Introduction
It has been more than a century since Erlang published his pioneering work [
1] in 1909. Since then, many authors have studied the M/D/c queue with Poisson arrivals, deterministic (constant) service time and c multiple servers. The M/D/c queue can be applied in wireless networks like the 5G (Fifth Generation) cellular network with fixed-size machine-to-machine packets. Crommelin studied the delay probability of the M/D/c system in [
2], while Franx worked out a simple solution for the waiting time distribution of this system. In addition, Jansen and Van Leeuwaarden have derived explicit solutions for the roots of the equation for the M/D/c systems [
3]. Tijms presented the mean number of customers for an M/D/c system in terms of the unsaturated probabilities of the number of customers in the system [
4]. Furthermore, Knessl has closed-form expressions for the unsaturated probabilities of number of customers in the M/D/2 system [
5]. On the other hand, Shortle et al. presented a method of solving the M/D/c system numerically [
6]. Seelen, Tijms and Van Hoorn have also published exact numerical results on the mean number of customers in the queue for the M/D/c system [
7].
The contributions of this paper are as follows. An alternate form of explicit root(s) of an equation for the M/D/c queue is/are presented in terms of the Lambert-W function. In addition, the unsaturated probabilities of the number of customers in the system are expressed in terms of the root(s). Finally, the mean number of customers in the system and the mean queueing delay in the system are derived in terms of the root(s). The mean queueing delay formula of the M/D/c system is not only exact and of closed-form, it is also a product of the form of the M/D/1 queueing formula and another term to account for the multiple servers. These are the main contributions of this paper.
2. The M/D/c Queueing Delay Formula
We start to find the solution(s) of the roots of the following equation for the roots of the M/D/c system with c servers and the system utilization, denoted by
ρ, from [
6].
and
λ is the arrival rate to the system and
b is the deterministic (constant) service time. Let
W(
x) be the Lambert-W function and it is given by
The solution(s) of the roots of equation (1),
zk,
k = 1, 2, …,
c – 1, are as follows. If the value of c is even,
If the value of c is odd,
The roots for
c = 1, 2, 3, 4, 5, 6 are shown in
Table 1. One of the advantages of using the Lambert-W function is that it is readily available in the Matlab software for numerical computation.
Let
pk be the unsaturated probabilities of the number of customers in the system, where
k = 0, 1, 2, …,
c – 1. From [
6],
the probabilities of
pk ’s are solved using the following equation.
where
Solving equation (8), we have the unsaturated probabilities of the number of customers in the M/D/c system for c = 1, 2, 3, 4, 5, 6. They are tabulated in
Table 2,
Table 3,
Table 4,
Table 5 and
Table 6.
Let us define the followings.
From [
5], the mean number of customers in the M/D/c queue,
Nq, is given by
Expressing the probabilities of
pk’s in terms of the
gk’s in equations (11) and (12), we have
Furthermore, substituting them into equation (14), we have
Applying Little’s Law, the mean queueing delay of the M/D/c system, denoted by
W, is given by
Simplifying and rearranging some terms, we have
Note that the square of the deterministic (constant) service time is equal to the second moment of its service time. Substituting equation (16) into equation (18), we have
Equation (19) is the exact closed-form M/D/c queueing delay formula. They are tabulated in
Table 7 for the number of servers,
c = 2, 3, 4, 5, 6. Note that
c = 2, 3, …, in general.
On further simplification of (19), we have
Note that the service time, b, needs not be set to 1 and there is no need for any rescaling after that.
3. Results
We present the results of the M/D/c system with the number of servers,
c = 2, 3, 4, 5, 6 and the constant service time,
b = 1, to be consistent with [
7] in this results section. The mean number of customers in the system,
Nq, from analysis in
Section 2 is compared with known exact numerical results from [
7], while the exact mean queueing delay,
W, from the analysis in
Section 2 is also presented. Note that the constant service time needs not be set to 1 in (20) in general.
Table 8 presents the results for the mean number of customers with our results as compared to available exact numerical results from [
7]. Both results match excellently. Only available numerical results from [
7] are presented for each system utilization value.
Table 9 presents our results for the mean queueing delay with our exact closed-form M/D/c queueing delay formula. The results agree with Little’s Law and can also be derived from
Table 8.
4. Discussion
The result in
Section 2 provides an exact closed-form M/D/c queueing delay formula. It is in a form that is extending the form of M/D/1 queue expression as well. In addition, this result is dependent on the roots which are expressed in terms of the Lambert-W function. Furthermore, this function is readily available in software like Matlab for ease of numerical computation. The M/D/c queueing system also has applications in wireless networks and other fields.
5. Conclusions
An exact closed-form M/D/c queueing delay formula has been derived. The final form of the result does not require the solution of a set of linear equations anymore. This formula depends on the roots of the solution to an equation for the M/D/c system. Furthermore, the roots are presented in an alternate explicit form in terms of the Lambert-W function for ease of numerical computation. Numerical results show that these results are in excellent agreement with published exact numerical results.
The significance of this queueing delay formula presented in this paper is that it moves us a step closer to an exact closed-form M/G/c queueing system as the exact closed form M/M/c queueing delay formula is also well-known [
8].
Author Contributions
The author has sole contribution to this paper.
Data Availability Statement
The numerical data is available in [
7].
Acknowledgments
The author would like to thank Prof. Jon W Mark, Prof. Kin Mun Lye and Prof. Kee Chaing Chua for the queueing theory/communication networks courses that they have taught the author.
Conflicts of Interest
The author declares no conflicts of interest.
References
- Erlang, A.K. The Theory of Probabilities and telephone conversations. Nyt. Tidsskr. Mat. Ser. 1909, 20, 33–39. [Google Scholar]
- Crommelin, C.D. Delay Probability Formulae When the Holding Times are Constant. POEE 1932, 25, 41–50. [Google Scholar]
- Janssen, A.J.E.M.; Van Leeuwaarden, J.S.H. Back to the Roots of the M/D/s Queue and the Works of Erlang, Crommelin and Pollaczek. Stat. Neer 2008, 62, 299–313. [Google Scholar] [CrossRef]
- Tijms, H.C. Stochastic Models: An Algorithmic Approach, 1st, ed.; Wiley, 1994; p. 290. [Google Scholar]
- Knessl, C.; Matkowsky, B.J.; Schuss, Z.; Tier, C. An Integral Equation Approach to the M/G/2 Queue. OR 1990, 38, 506–518. [Google Scholar] [CrossRef]
- Shortle, J.F.; Thompson, J.M.; Gross, D.; Harris, C.M. Fundamentals of Queueing Theory; John Wiley & Sons, 2018; pp. 330–332. [Google Scholar]
- Seelen, L.P.; Tijms, H.C.; Van Hoorn, M.H. Tables for Multi-Sever Queues. North–Holland, 1985. [Google Scholar]
- Wong, D.T.C.; Kong, P.-Y.; Liang, Y.C.; Chua, K.C.; Mark, J.W. Wireless Broadband Networks; John Wiley & Sons, 2009; pp. 493–495. [Google Scholar]
Table 1.
Roots of the M/D/c system with the number of servers, c = 2, 3, 4, 5, 6.
Table 1.
Roots of the M/D/c system with the number of servers, c = 2, 3, 4, 5, 6.
c |
z1
|
z2
|
z2
|
z4
|
z5
|
2 |
|
|
|
|
|
3 |
|
|
|
|
|
4 |
|
|
|
|
|
5 |
|
|
|
|
|
6 |
|
|
|
|
|
Table 2.
The unsaturated probability of the number of customers in the M/D/c system with the number of servers, c = 2.
Table 2.
The unsaturated probability of the number of customers in the M/D/c system with the number of servers, c = 2.
k |
pk |
1 |
|
Table 3.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 3.
Table 3.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 3.
k |
pk |
1 |
|
2 |
|
Table 4.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 4.
Table 4.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 4.
k |
pk |
1 |
|
2 |
|
3 |
|
Table 5.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 5.
Table 5.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 5.
k |
pk |
1 |
|
2 |
|
3 |
|
4 |
|
Table 6.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 6.
Table 6.
The unsaturated probabilities of the number of customers in the M/D/c system with the number of servers, c = 6.
k |
pk |
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
Table 7.
The M/D/c Queueing Delay Formulae with the number of servers, c = 2, 3, 4, 5, 6.
Table 7.
The M/D/c Queueing Delay Formulae with the number of servers, c = 2, 3, 4, 5, 6.
c |
M/D/c Queueing Delay Formulae |
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
Table 8.
The Mean Number of Customers for the M/D/c system with the number of servers, c = 2, 3, 4, 5, 6.
Table 8.
The Mean Number of Customers for the M/D/c system with the number of servers, c = 2, 3, 4, 5, 6.
ρ\c |
2 |
3 |
4 |
5 |
6 |
|
Ours |
[7] |
Ours |
[7] |
Ours |
[7] |
Ours |
[7] |
Ours |
[7] |
0.1 |
0.0012 |
|
0.0003 |
|
0.0001 |
|
0.0000 |
|
0.0000 |
|
0.2 |
0.0097 |
0.0097 |
0.0039 |
0.0039 |
0.0016 |
|
0.0007 |
|
0.0003 |
|
0.3 |
0.0332 |
0.0332 |
0.0181 |
0.0181 |
0.0101 |
0.0101 |
0.0058 |
0.0058 |
0.0033 |
|
0.4 |
0.0826 |
0.0826 |
0.0540 |
0.0540 |
0.0363 |
0.0363 |
0.0249 |
0.0249 |
0.0172 |
0.0172 |
0.5 |
0.1767 |
0.1767 |
0.1308 |
0.1308 |
0.0993 |
0.0993 |
0.0766 |
0.0766 |
0.0598 |
0.0598 |
0.6 |
0.3516 |
0.3516 |
0.2851 |
0.2851 |
0.2361 |
0.2361 |
0.1982 |
0.1982 |
0.1680 |
0.1680 |
0.7 |
0.6911 |
0.6911 |
0.6011 |
0.6011 |
0.5312 |
0.5312 |
0.4743 |
0.4743 |
0.4267 |
0.4267 |
0.8 |
1.4453 |
1.445 |
1.3294 |
1.329 |
1.2355 |
1.236 |
1.1562 |
1.156 |
1.0875 |
1.088 |
0.9 |
3.8645 |
3.864 |
3.7204 |
3.720 |
3.5999 |
3.600 |
3.4951 |
3.495 |
3.4018 |
3.402 |
Table 9.
The Mean Queueing Delay for the M/D/c system with the number of servers, c = 2, 3, 4, 5, 6.
Table 9.
The Mean Queueing Delay for the M/D/c system with the number of servers, c = 2, 3, 4, 5, 6.
ρ\c |
2 |
3 |
4 |
5 |
6 |
|
Ours |
Ours |
Ours |
Ours |
Ours |
|
0.1 |
0.0062 |
0.0009 |
0.0002 |
0.0000 |
0.0000 |
0.2 |
0.0242 |
0.0066 |
0.0021 |
0.0007 |
0.0002 |
0.3 |
0.0553 |
0.0201 |
0.0085 |
0.0039 |
0.0019 |
0.4 |
0.1033 |
0.0450 |
0.0227 |
0.0124 |
0.0072 |
0.5 |
0.1767 |
0.0872 |
0.0497 |
0.0307 |
0.0199 |
0.6 |
0.2930 |
0.1584 |
0.0984 |
0.0661 |
0.0467 |
0.7 |
0.4936 |
0.2862 |
0.1897 |
0.1355 |
0.1016 |
0.8 |
0.9033 |
0.5539 |
0.3861 |
0.2891 |
0.2266 |
0.9 |
2.1469 |
1.3779 |
1.0000 |
0.7767 |
0.6300 |
|
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. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).