In this paper we describe a perfect simulation algorithm for the stable M/G/c queue. Sigman (2011) showed how to build a dominated coupling-from-the-past algorithm for perfect simulation of the super-stable M/G/c queue operating under first-come-first-served discipline. Sigman's method used a dominating process provided by the corresponding M/G/1 queue (using Wolff's sample path monotonicity, which applies when service durations are coupled in order of initiation of service). The method exploited the fact that the workload process for the M/G/1 queue remains the same under different queueing disciplines, in particular under the processor sharing discipline, for which a dynamic reversibility property holds. We generalise Sigman's construction to the stable case by comparing the M/G/c queue to a copy run under random assignment. This allows us to produce a naïve perfect simulation algorithm based on running the dominating process back to the time it first empties. We also construct a more efficient algorithm that uses sandwiching by lower and upper processes constructed as coupled M/G/c queues started respectively from the empty state and the state of the M/G/c queue under random assignment. A careful analysis shows that appropriate ordering relationships can still be maintained, so long as service durations continue to be coupled in order of initiation of service. We summarise statistical checks of simulation output, and demonstrate that the mean run-time is finite so long as the second moment of the service duration distribution is finite.

A new Markov process is introduced, describing growth or spread in two dimensions, via the aggregation of particles or the filling of cells. States of the process are configurations of part of the boundary of the growing aggregate, and transitions are captures or escapes of single particles. For suitably chosen transition rates, the process is dynamically reversible, leading to an explicit stationary distribution and a statistical description of the boundary. The growth rate is calculated and growth behaviour described. Different asymptotic relations between transition rates lead to different growth patterns or regimes. Besides the regimes familiar in polymer crystal growth, several new ones are described. The aggregate can have a porous structure resembling thin solid films deposited from vapour. Two measures of porosity, one for the boundary and one for the bulk, are calculated. The process is relevant to growing colonies of bacteria or the like, to the spread of epidemics and grass or forest fires, and to voter models.