Tutorials › Stochastic Processes › Why Performance Changes Sharply Near Saturation

The M/M/1 Queue · Tutorial 642 of 994

Why Performance Changes Sharply Near Saturation

Quantify the sharp increase in mean queue length and delay near saturation, and compare stable M/M/1 queues using exact calculations.

Intermediate 9 min read

What You'll Learn

  • Identify the stability condition needed for stationary M/M/1 performance formulas.
  • Use utilization and service rate to calculate mean queue length and delay.
  • Explain the divergence of these means as utilization approaches one.
  • Compare performance at several exact, stable traffic intensities.
  • Detect the mistaken assumption that performance grows linearly with utilization.

Why a Small Change in Utilization Can Matter

In a stable queue, increasing utilization means that the server is busy for a larger fraction of time. But the effect on congestion is not proportional to that increase: when utilization is already high, a small further increase can produce a large change in the mean queue length and delay. The reason is that random variation in arrivals and service has less spare capacity to be absorbed.

Consider an M/M/1 queue on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), with continuous time \(t\ge0\). Let \(X_t\in S=\{0,1,2,\ldots\}\) be the number of customers in the system at time \(t\). Arrivals occur at rate \(\lambda>0\); service times are exponential with rate \(\mu>0\), and there is one server. Write \(\rho=\lambda/\mu\) for utilization. In this tutorial, performance means stationary averages, so we require the stability condition \(0<\rho<1\).

Write \(L=\mathbb{E}[X_t]\) for the stationary mean number in the system, \(L_q\) for the stationary mean number waiting in queue, \(W\) for a customer's mean time in the system, and \(W_q\) for the mean time waiting before service. The standard M/M/1 formulas established in the earlier material are

$$ L=\frac{\rho}{1-\rho},\qquad L_q=\frac{\rho^2}{1-\rho},\qquad W=\frac{1}{\mu-\lambda}=\frac{1}{\mu(1-\rho)},\qquad W_q=\frac{\rho}{\mu-\lambda}=\frac{\rho}{\mu(1-\rho)}. $$

These are stationary means, not guarantees for individual customers. In particular, a finite mean at any fixed \(\rho<1\) does not mean that every customer experiences a short delay. The formulas quantify how the averages behave as the stability boundary is approached.

The Divergence Near Saturation

Theorem: For an M/M/1 queue with fixed \(\mu>0\) and \(0<\rho<1\), the stationary means satisfy \(\lim_{\rho\uparrow1}(1-\rho)L=1\), \(\lim_{\rho\uparrow1}(1-\rho)L_q=1\), \(\lim_{\rho\uparrow1}(1-\rho)W=1/\mu\), and \(\lim_{\rho\uparrow1}(1-\rho)W_q=1/\mu\). Thus all four means diverge, with leading growth proportional to \(1/(1-\rho)\).

Proof. Substitute the formulas above and simplify for each \(0<\rho<1\): \((1-\rho)L=(1-\rho)\rho/(1-\rho)=\rho\), which tends to \(1\) as \(\rho\uparrow1\). Likewise, \((1-\rho)L_q=\rho^2\), which tends to \(1\); \((1-\rho)W=1/\mu\), which is constant in \(\rho\); and \((1-\rho)W_q=\rho/\mu\), which tends to \(1/\mu\). Since \(1-\rho\) tends to zero through positive values and each product has a strictly positive limit, each corresponding mean tends to infinity. For example, if \((1-\rho)L\to1\), then for \(\rho\) sufficiently close to one, \(L>(1/2)/(1-\rho)\), which grows without bound. The same argument applies to the other three means. \(\square\)

The factor \(1/(1-\rho)\) is the key scale. It is not a finite value at \(\rho=1\): the theorem describes a limit taken through stable queues, each with \(\rho<1\). At fixed \(\mu\), the mean sojourn time \(W\) grows exactly as the reciprocal of the unused service fraction \(1-\rho\). The other three means have the same leading growth, although their exact numerators differ.

Result (exact sensitivity): Holding \(\mu>0\) fixed, for \(0<\rho<1\), \(\frac{dL}{d\rho}=\frac{1}{(1-\rho)^2}\) and \(\frac{dW}{d\rho}=\frac{1}{\mu(1-\rho)^2}\). For \(0<\rho_1<\rho_2<1\), the amplification ratios are \(\frac{L(\rho_2)}{L(\rho_1)}=\frac{\rho_2(1-\rho_1)}{\rho_1(1-\rho_2)}\) and \(\frac{W(\rho_2)}{W(\rho_1)}=\frac{1-\rho_1}{1-\rho_2}\).

Proof. With \(\mu\) fixed, \(L(\rho)=\rho/(1-\rho)\). The quotient rule gives \(\frac{dL}{d\rho}=[(1-\rho)+\rho]/(1-\rho)^2=1/(1-\rho)^2\). Also \(W(\rho)=1/[\mu(1-\rho)]\), so differentiation gives \(\frac{dW}{d\rho}=1/[\mu(1-\rho)^2]\). For the ratios, divide the positive finite means at \(\rho_2\) by those at \(\rho_1\). For \(L\), this yields \([\rho_2/(1-\rho_2)]/[\rho_1/(1-\rho_1)] =\rho_2(1-\rho_1)/[\rho_1(1-\rho_2)]\). For \(W\), the common factor \(1/\mu\) cancels, leaving \([1/(1-\rho_2)]/[1/(1-\rho_1)]=(1-\rho_1)/(1-\rho_2)\). All denominators are positive under the stated assumptions. \(\square\)

The derivative of each mean becomes large as \(\rho\) approaches one. This is a precise sense in which performance changes sharply near saturation: an equal-sized increase in utilization has a larger effect when starting from a higher utilization. For \(L\), the exact ratio also includes the numerator \(\rho\); for delay \(W\), the ratio depends only on the remaining capacity fractions.

Comparing Stable Traffic Intensities

Worked Example: Four Utilizations at One Service Rate

Suppose a server processes customers at rate \(\mu=12\) per hour. Compare stationary mean queue length and time measures at \(\rho=\tfrac12,\tfrac45,\tfrac9{10},\tfrac{19}{20}\). Each utilization is strictly between zero and one, so each queue is stable. The corresponding arrival rates, calculated as \(\lambda=\rho\mu\), are \(6,\tfrac{48}{5},\tfrac{54}{5},\tfrac{57}{5}\) customers per hour, respectively.

At \(\rho=\tfrac12\), \(1-\rho=\tfrac12\), so \(L=(\tfrac12)/(\tfrac12)=1\) and \(L_q=(\tfrac12)^2/(\tfrac12)=\tfrac12\). The mean time in the system is \(W=1/[12(\tfrac12)]=\tfrac16\) hour, or \(10\) minutes; the mean queue wait is \(W_q=(\tfrac12)/[12(\tfrac12)]=\tfrac1{12}\) hour, or \(5\) minutes.

At \(\rho=\tfrac45\), \(1-\rho=\tfrac15\), so \(L=(\tfrac45)/(\tfrac15)=4\) and \(L_q=(\tfrac{16}{25})/(\tfrac15)=\tfrac{16}{5}\). Here \(W=1/[12(\tfrac15)]=\tfrac5{12}\) hour, or \(25\) minutes, and \(W_q=(\tfrac45)/[12(\tfrac15)]=\tfrac13\) hour, or \(20\) minutes.

At \(\rho=\tfrac9{10}\), \(1-\rho=\tfrac1{10}\), so \(L=(\tfrac9{10})/(\tfrac1{10})=9\) and \(L_q=(\tfrac{81}{100})/(\tfrac1{10})=\tfrac{81}{10}\). Also, \(W=1/[12(\tfrac1{10})]=\tfrac56\) hour, or \(50\) minutes, and \(W_q=(\tfrac9{10})/[12(\tfrac1{10})]=\tfrac34\) hour, or \(45\) minutes.

At \(\rho=\tfrac{19}{20}\), \(1-\rho=\tfrac1{20}\), so \(L=(\tfrac{19}{20})/(\tfrac1{20})=19\) and \(L_q=(\tfrac{361}{400})/(\tfrac1{20})=\tfrac{361}{20}\). Further, \(W=1/[12(\tfrac1{20})]=\tfrac53\) hours, or \(100\) minutes, and \(W_q=(\tfrac{19}{20})/[12(\tfrac1{20})]=\tfrac{19}{12}\) hours, or \(95\) minutes.

Check: For each utilization, the difference \(L-L_q\) is the mean number in service: the four differences are \(1-\tfrac12=\tfrac12\), \(4-\tfrac{16}{5}=\tfrac45\), \(9-\tfrac{81}{10}=\tfrac9{10}\), and \(19-\tfrac{361}{20}=\tfrac{19}{20}\). Each equals \(\rho\). Also, the four values of \(W\) in minutes are \(10,25,50,100\); each is \(60/[12(1-\rho)]\), as required by the formula.

The comparison shows why the final increments are consequential. Increasing utilization from \(\tfrac9{10}\) to \(\tfrac{19}{20}\) is an increase of only \(\tfrac1{20}\), but mean time in the system rises from \(50\) to \(100\) minutes. The service rate stays fixed; the change comes from halving the remaining fraction \(1-\rho\).

Worked Example: The Same Utilization Increase at Different Starting Points

Take \(\mu=10\) customers per hour and compare a change from \(\rho_1=\tfrac45\) to \(\rho_2=\tfrac9{10}\). The arrival rates are \(\lambda_1=10(\tfrac45)=8\) and \(\lambda_2=10(\tfrac9{10})=9\) per hour. Both queues are stable.

At \(\rho_1=\tfrac45\), \(L_1=(\tfrac45)/(\tfrac15)=4\) and \(W_1=1/[10(\tfrac15)]=\tfrac12\) hour. At \(\rho_2=\tfrac9{10}\), \(L_2=(\tfrac9{10})/(\tfrac1{10})=9\) and \(W_2=1/[10(\tfrac1{10})]=1\) hour. Thus this increase of \(\tfrac1{10}\) in utilization raises the mean system size by \(9-4=5\), and raises mean time in the system by \(1-\tfrac12=\tfrac12\) hour.

Check by the exact ratios: The mean-size ratio is \((\tfrac9{10})(\tfrac15)/[(\tfrac45)(\tfrac1{10})]=(\tfrac9{50})/(\tfrac2{25})=\tfrac94\), which equals \(9/4\). The delay ratio is \((\tfrac15)/(\tfrac1{10})=2\), which equals \(1/(\tfrac12)\). Both ratios agree with direct calculation.

The denominator explains the difference between “a small utilization increase” and “a small performance change.” Utilization itself changes by the same amount, but its effect depends on the starting slack. The exact derivative \(dW/d\rho=1/[\mu(1-\rho)^2]\) is four times as large at \(\rho=\tfrac9{10}\) as at \(\rho=\tfrac45\), since \((\tfrac15)^2/(\tfrac1{10})^2=4\).

Worked Example: Halving the Spare-Capacity Fraction

Let \(\mu=20\) customers per hour. Compare \(\rho=\tfrac9{10}\), with spare-capacity fraction \(1-\rho=\tfrac1{10}\), and \(\rho=\tfrac{19}{20}\), with spare-capacity fraction \(1-\rho=\tfrac1{20}\). Both are stable, and the second queue has half the spare-capacity fraction of the first.

At \(\rho=\tfrac9{10}\), \(L=(\tfrac9{10})/(\tfrac1{10})=9\), and \(W=1/[20(\tfrac1{10})]=\tfrac12\) hour. At \(\rho=\tfrac{19}{20}\), \(L=(\tfrac{19}{20})/(\tfrac1{20})=19\), and \(W=1/[20(\tfrac1{20})]=1\) hour. The mean system size increases by \(19-9=10\), while mean time in the system doubles from \(\tfrac12\) hour to \(1\) hour.

Check: The exact delay ratio is \((\tfrac1{10})/(\tfrac1{20})=2\), matching \(1/(\tfrac12)=2\). For mean system size, the exact ratio is \(19/9\), not \(2\); the formula \(L=1/(1-\rho)-1\) shows why. Halving the spare fraction doubles \(1/(1-\rho)\), while subtracting one leaves a ratio slightly below two.

Stability Boundary and Common Mistake

The condition \(\rho<1\) is essential: the displayed formulas describe a stationary queue only on the stable side of the boundary. For example, if \(\lambda=\mu\), then \(\rho=1\). The birth-and-death balance equations for a stationary M/M/1 queue require \(\pi_n=\pi_0\rho^n\) for \(n\ge0\). At \(\rho=1\), these weights are all \(\pi_0\); if \(\pi_0>0\), their sum is infinite, and if \(\pi_0=0\), their sum is zero. Neither choice gives a probability distribution. Thus there is no stationary distribution to which the stationary mean formulas could apply at \(\rho=1\). The limit theorem concerns stable queues with \(\rho<1\); it does not assign a stationary mean at the boundary.

Common mistake: treating congestion as linear in utilization. A mistaken approximation is \(L=\rho\), which would give \(L=\tfrac9{10}\) at \(\rho=\tfrac9{10}\). The correct substitution is \(L=\rho/(1-\rho)=(\tfrac9{10})/(\tfrac1{10})=9\). The incorrect calculation drops the denominator that records spare capacity; near saturation, that denominator is precisely what drives the sharp increase.

When comparing queues, first confirm that each utilization is below one and then use the same performance measure and units. In particular, distinguish the number waiting, \(L_q\), from the number in the system, \(L\), and the queue wait, \(W_q\), from the total time in the system, \(W\). The formulas make each comparison exact, but only if the quantities being compared are defined consistently.

Check Your Understanding

Use the stable M/M/1 formulas and verify each numerical answer by substitution.

  1. With fixed \(\mu\), what happens to \((1-\rho)W\) as \(\rho\uparrow1\), and why does this imply that \(W\) diverges?
  2. For \(\mu=8\) per hour and \(\rho=\tfrac34\), calculate \(L\) and \(W\). State the time unit for \(W\).
  3. At fixed \(\mu\), compare \(W\) at \(\rho=\tfrac45\) and \(\rho=\tfrac9{10}\) using the exact amplification ratio.
  4. Why can the stationary mean formulas not be applied at \(\rho=1\)?
  5. Explain why halving \(1-\rho\) doubles \(W\), but does not in general exactly double \(L\).