An approximate numerical solution for multiclass preemptive priority queues with general service time distributions

In this paper an approximate numerical solution for a multiclass preemptive priority single server queue is developed. The arrival process of each class follows a Poisson distribution. The service time distribution must have a rational Laplace transform, but is otherwise arbitrary and may be different for different classes. The work reported here was motivated by a desire to compute the equilibrium probability distribution of networks containing preemptive priority servers. Such networks are frequently encountered when modeling computer systems, medical care delivery systems and communication networks. We wish to use an iterative technique which constructs a series of two station networks consisting of one station from the original network and one “complementary” station whose behavior with respect to the original station mimics that of the rest of the network. At each iteration, it is necessary to compute the equilibrium probability distribution of one or more preemptive priority queues. Although such queues have been studied for some time, the resulting solutions have most often been developed utilizing transforms or probability generating functions, e.g. Jaiswal [1968]; in many cases of interest, inversion has not been attempted. Miller [1981] presented explicit solutions for two class priority queues but Miller's work, which is based on that of Neuts, is limited to exponential service times and two classes. The approach presented here is applicable to many classes and to more general service time distributions than have previously been considered. The algorithm utilizes a bootstrap approach, a concept borrowed from dynamic programming. The solution for class 1 is trivial. Once we have solved the system with k different classes, we have available all the necessary information to solve the system with k+l classes. We shall assume that each class has a distinct service time distribution G k , with mean g k and variance s 2 k . Let class k have preemptive priority over class l if k Successive steps of the algorithm are based upon the machine breakdown and repair model, previously used by Keilson [1962]; White and Christie [1958]; Gaver [1962]; and Avi-Itzhak and Noar [1963] to model preemptive priority queues. When we are solving for the equilibrium probability distribution of class k jobs, k > 1, we consider a model with one machine whose service time distribution is G k . The breakdown rate of the machine is the sum of the arrival rates of all higher priority jobs. The downtime or repairtime of the machine has mean γ k-l and variance s 2 k-l ; these parameters are the mean and variance of the busy period in a preemptive priority system with k-l classes. The first step in the solution of all the machine breakdown and repair models considered herein is to construct the infinitesimal generator