TITLE:  Gaussian Skewness Approximation for Dynamic Rate Multi-Server Queues

SPEAKER:  Jamol Pender, Ph.D. Candidate in ORFE at Princeton

ABSTRACT:

Large scale systems such as customer contact centers, like telephone call centers, as well as healthcare centers, like hospitals, have customer inflow-outflow dynamics with many common features. The customer arrival patterns may have time of day or seasonal effects. Moreover, customer population sizes tend to be large where the individual actions are intrinsic and independent of other customer actions and there are multiple service agents, so many customers have access to services in parallel. Finally, arriving customers engaging in service may be delayed if all the available agents are busy. Moreover, these waiting customers may decide to leave the systems if they feel that their delay in receiving service is excessively long.

A fundamental Markov process, with dynamic rates, queueing model for large service systems is a multiserver queue with non-homogeneous Poisson arrivals as well as service and customer abandonment times that are exponentially distributed. An asymptotic scaling for these queues leads to both functional strong law of large numbers and central limit theorems. The first yields a dynamical system that we call a fluid model. The second scaled limit yields a diffusion process model that is Gaussian under conditions involving the fluid model not lingering too closely to the number of servers. Finally, the fluid model coupled with the mean and covariance of the Gaussian model is also a dynamical system. These results are a special case of a general asymptotic theory for Markovian service networks as derived in the seminal work of Mandelbaum, Massey, and Reiman [1998].

In practice, these results yield useful Gaussian approximations to the transient queue length distribution. However, these estimates tend not to work as well when the lingering effect is significant. We can improve these methods by introducing a new technique that is called the Gaussian-skewness approximation (GSA). It is the special case of a general Hermite polynomial expansion for a Gaussian random variable. We then obtain dynamical systems that approximate the mean, variance and higher cumulant moments of the queueing process for a more accurate, non-Gaussian estimation.