Tutorials › Stochastic Processes › Worked Problems in Regenerative Maintenance and Queues

Renewal-Reward and Regenerative Processes · Tutorial 621 of 994

Worked Problems in Regenerative Maintenance and Queues

Learn to turn regenerative cycles into exact long-run performance measures for maintenance systems and a capacity-one queue.

Advanced 11 min read

What You'll Learn

  • Identify cycle boundaries that make successive cycles independent and identically distributed.
  • Choose cycle rewards for availability, maintenance cost, and time spent in a queue state.
  • Apply the renewal-reward ratio with conditions controlling the unfinished cycle.
  • Compute exact availability and maintenance cost rates from a replacement policy.
  • Calculate queue occupancy and rejected-arrival rates from idle-busy cycles.
  • Distinguish a ratio of expected cycle quantities from an average of cycle-by-cycle ratios.

From Regeneration to a Performance Rate

In an operating system, a useful long-run quantity is often a fraction of time or a cost per unit time. Regeneration makes these quantities calculable: choose epochs at which the system starts afresh, measure the time and reward in one cycle, and compare their expectations. The choice of cycle boundary matters. It must make the cycles identically distributed and independent, and the reward must represent the quantity being measured.

For a continuous-time process \((X_t)_{t\ge0}\) on a state space \(S\), let \(0=T_0<T_1<T_2<\cdots\) be regeneration epochs, and write \(C_k=T_k-T_{k-1}\) for cycle lengths. Assume the cycles have the same distribution and are independent, with \(C_k>0\) almost surely and \(0<\mathbb{E}C_1<\infty\). Let \(R_k\) be the net reward in cycle \(k\). The cycle reward may be signed, as for profit minus cost. For a reward accrued over time, also track its total variation \(V_k\): the total amount of positive and negative reward accumulated in absolute value during cycle \(k\).

Cycle-rate result: Suppose \((C_k,R_k,V_k)\) are independent and identically distributed, \(\mathbb{E}C_1\in(0,\infty)\), \(\mathbb{E}|R_1|<\infty\), and \(\mathbb{E}V_1<\infty\). Assume the absolute reward accumulated in any part of cycle \(k\) is at most \(V_k\). Then the reward accumulated by time \(t\), divided by \(t\), converges almost surely to \(\mathbb{E}R_1/\mathbb{E}C_1\).

Proof. Let \(S_n=\sum_{k=1}^n C_k\), and let \(N(t)=\max\{n:S_n\le t\}\), the number of completed cycles by time \(t\). The standard iid strong law gives \(S_n/n\to\mathbb{E}C_1\) and \(\sum_{k=1}^n R_k/n\to\mathbb{E}R_1\) almost surely. Since the cycle lengths are positive and have finite positive mean, \(N(t)\to\infty\). The inequalities \(S_{N(t)}\le t<S_{N(t)+1}\), together with \(S_n/n\to\mathbb{E}C_1\), imply \(N(t)/t\to1/\mathbb{E}C_1\). Therefore the reward from completed cycles, divided by \(t\), converges to \(\mathbb{E}R_1/\mathbb{E}C_1\).

It remains to control the part of the current cycle that has not been completed. The strong law applied to \(V_k\) gives \(\sum_{k=1}^n V_k/n\to\mathbb{E}V_1\). Subtracting the corresponding limit for \(n-1\) from that for \(n\) shows \(V_n/n\to0\). Thus \(V_{N(t)+1}/t\to0\): indeed, \(V_{N(t)+1}/(N(t)+1)\to0\) while \((N(t)+1)/t\to1/\mathbb{E}C_1\). The unfinished-cycle reward is bounded in absolute value by \(V_{N(t)+1}\), so its contribution divided by \(t\) vanishes. Adding it to the completed-cycle reward proves the claim. \(\square\)

The total-variation condition is important for signed rewards: a small net reward for a cycle does not necessarily mean little reward was exchanged within it. In the examples below, rewards are nonnegative quantities such as operating time, cost, or time in a state. Their within-cycle accumulated amount is bounded by their cycle total, so the condition is automatic when that total has finite expectation.

Worked Example: Age-Based Maintenance

Worked Example: Availability and Maintenance Cost

A machine begins operating at the end of each maintenance action. Its lifetime \(L\) is 2 or 5 time units, each with probability \(\tfrac12\). If it fails at time 2, corrective maintenance takes 1 time unit and costs 5 units. If it survives to the planned replacement age 3, preventive maintenance takes \(\tfrac12\) time unit and costs 2 units. A cycle ends when maintenance is complete. Find the long-run availability and maintenance cost per unit time.

Take the cycle start as the regeneration epoch. Cycles are iid because the lifetime and maintenance outcome are newly sampled in the same way each time. The operating reward is the time operating in a cycle; the cost reward is the maintenance cost paid in that cycle. The two possible cycle outcomes are:

  • Failure: probability \(\tfrac12\), operating time \(2\), maintenance time \(1\), cycle length \(2+1=3\), cost \(5\).
  • Survival to age 3: probability \(\tfrac12\), operating time \(3\), maintenance time \(\tfrac12\), cycle length \(3+\tfrac12=\tfrac72\), cost \(2\).

The expected cycle length is \(\tfrac12(3)+\tfrac12(\tfrac72)=\tfrac{13}{4}\). Expected operating time is \(\tfrac12(2)+\tfrac12(3)=\tfrac52\). Therefore availability, the long-run fraction of time operating, is

$$ \frac{\mathbb{E}[\text{operating time in a cycle}]}{\mathbb{E}[\text{cycle length}]} =\frac{\tfrac52}{\tfrac{13}{4}}=\frac{10}{13}. $$

Expected maintenance cost per cycle is \(\tfrac12(5)+\tfrac12(2)=\tfrac72\). The long-run cost rate is consequently

$$ \frac{\mathbb{E}[\text{cost per cycle}]}{\mathbb{E}[\text{cycle length}]} =\frac{\tfrac72}{\tfrac{13}{4}}=\frac{14}{13}. $$

Checks: Expected maintenance time is \(\tfrac12(1)+\tfrac12(\tfrac12)=\tfrac34\), so its fraction of time is \((\tfrac34)/(\tfrac{13}{4})=\tfrac3{13}\). This complements the operating fraction, since \(\tfrac{10}{13}+\tfrac3{13}=1\). Both outcomes have nonnegative bounded cycle rewards, and the expected cycle length is positive and finite, as required.

State Occupancy as a Cycle Reward

Occupancy is calculated by using time in a state as the cycle reward. For a finite state space, let \(A_i\) be the time spent in state \(i\) during one cycle. Each \(A_i\) is nonnegative and at most the cycle length. The same cycle-rate result applies to every state.

Occupancy formula: If the process regenerates in iid cycles of positive finite mean length and \(S\) is finite, then the long-run fraction of time in state \(i\) is \(\mathbb{E}A_i/\mathbb{E}C_1\). These fractions sum to 1.

Proof. Apply the cycle-rate result to the reward that accrues at rate 1 while the process is in state \(i\), and at rate 0 otherwise. Its cycle reward is \(A_i\), and its total variation is also \(A_i\le C_1\), so its time average converges to \(\mathbb{E}A_i/\mathbb{E}C_1\). In every cycle, the state occupation times partition the cycle, so \(\sum_{i\in S}A_i=C_1\). Taking expectations and dividing by the positive finite \(\mathbb{E}C_1\) gives \(\sum_{i\in S}\mathbb{E}A_i/\mathbb{E}C_1=1\). \(\square\)

Worked Example: Occupancy in a Capacity-One Queue

Consider a continuous-time queue with state space \(S=\{0,1\}\), where 0 means empty and 1 means one customer present. Arrivals occur according to a Poisson process of rate \(\lambda=\tfrac23\). An arrival when the system is empty starts service; arrivals during service are rejected. Service times are independent exponential random variables with rate \(\mu=\tfrac43\), independent of arrivals. Find the long-run fraction of time occupied and the rate of rejected arrivals.

Regenerate whenever a service completion leaves the system empty. From that epoch, the idle time \(I\) until the next arrival is exponential with rate \(\lambda\), and the ensuing busy time \(B\) is exponential with rate \(\mu\). Independent increments and the model's fresh service times make successive idle-busy cycles iid. Their expected durations are \(\mathbb{E}I=1/\lambda=\tfrac32\) and \(\mathbb{E}B=1/\mu=\tfrac34\), so the expected cycle length is \(\tfrac32+\tfrac34=\tfrac94\). The cycle reward for occupancy is \(B\). Thus the fraction of time occupied is

$$ \frac{\mathbb{E}B}{\mathbb{E}(I+B)} =\frac{\tfrac34}{\tfrac94} =\frac13. $$

To count rejected arrivals, use the number of arrivals during the busy portion as the cycle reward. Conditional on \(B=b\), the expected number is \(\lambda b\), because arrivals are Poisson with rate \(\lambda\) and independent of service. Hence its expectation is \(\lambda\mathbb{E}B=(\tfrac23)(\tfrac34)=\tfrac12\). Dividing by the expected cycle length gives the rejected-arrival rate \((\tfrac12)/(\tfrac94)=\tfrac29\) per unit time.

Checks: The empty fraction is \((\tfrac32)/(\tfrac94)=\tfrac23\), which complements the occupied fraction \(\tfrac13\). Rejected arrivals occur only while busy, and their rate \(\tfrac29\) is less than the total arrival rate \(\tfrac23\). The expected busy time and expected rejected count are both finite, so the cycle rewards meet the integrability conditions.

Why the Ratio Is Not an Average of Cycle Rates

The cycle-rate result uses the ratio of expected reward to expected duration. In general, it does not use the expected value of the reward-to-duration ratio within a single cycle. Short cycles and long cycles occupy different amounts of calendar time, so averaging their individual rates gives them the wrong weights.

Worked Example: A Variable-Cycle Cost Rate

A maintenance operation has iid cycles. With probability \(\tfrac12\), a cycle lasts 1 time unit and costs 2 units; with probability \(\tfrac12\), it lasts 3 time units and costs 3 units. Find the long-run cost rate, and compare it with the average of the two cycle-by-cycle cost rates.

The expected cost is \(\tfrac12(2)+\tfrac12(3)=\tfrac52\). The expected cycle length is \(\tfrac12(1)+\tfrac12(3)=2\). The long-run cost rate is therefore \((\tfrac52)/2=\tfrac54\) cost units per time unit. By contrast, the average of the rates within the two cycle types is \(\tfrac12(2/1)+\tfrac12(3/3)=\tfrac32\), which is not the long-run rate.

Check: Over two representative cycles, one of each type, total cost is \(2+3=5\) and total time is \(1+3=4\), giving \(5/4\). This agrees with the renewal-reward calculation. The cycle-by-cycle average \(3/2\) gives the short and long cycles equal weight instead of weighting them by the time they contribute.

Common Pitfall: Net Reward Can Hide Large Within-Cycle Fluctuations

A completed-cycle formula for signed rewards needs care if one wants the reward accumulated at every time \(t\), not just at cycle boundaries. For instance, a cycle might accrue a large positive reward early and an equal negative reward later. Its net cycle reward would be zero, but the reward partway through the cycle could be large. Thus finite expected absolute net reward alone does not ensure that the unfinished-cycle contribution divided by time vanishes.

A concrete incorrect calculation would set the expected net cycle reward to zero, conclude that the long-run accumulated reward per unit time is zero, and ignore the current unfinished cycle. That conclusion does not follow from the net reward alone. The correction is to use a condition that controls within-cycle accumulation, such as finite expected total variation, as in the cycle-rate result above. For nonnegative rewards bounded by the cycle total, the condition holds directly. If one only asks for the sum of rewards over completed cycles at the regeneration epochs, the standard iid strong law gives the corresponding cycle-based ratio without needing to describe within-cycle fluctuations.

In applications, identify the cycle boundary first, specify exactly what reward is counted, and check its integrability and the expected cycle length before taking a ratio. For availability, count operating time rather than the number of operating periods. For cost, include each cost at the cycle in which it is incurred. For occupancy, count time in the state, not merely the number of visits.

Check Your Understanding

For each question, identify the cycle length and reward before calculating a rate.

  1. A regenerative system has expected cycle length 6 and expected operating time per cycle 4. What is its long-run availability?
  2. Why does a cycle reward equal to time spent in a state have total variation at most the cycle length?
  3. In the capacity-one queue example, what expected cycle reward counts rejected arrivals, and why is its expectation \(\lambda\mathbb{E}B\)?
  4. For the two-type cost example, calculate the long-run rate using expected cost divided by expected duration. Why is the average of the two cycle-specific rates different?
  5. What additional control is needed to infer a time-indexed long-run rate for signed rewards from net cycle rewards?