ARTICLE IN PRESS - Columbia Universityww2040/Liu_Whitt_EJOR_online.pdf · ARTICLE IN PRESS JID: EOR [m5G;August 20, 2016; ... Available online xxx Keywords: ... addition we consider
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
European Journal of Operational Research 0 0 0 (2016) 1–14
Contents lists available at ScienceDirect
European Journal of Operational Research
journal homepage: www.elsevier.com/locate/ejor
Theory and Methodology Paper
Stabilizing performance in a service system with time-varying arrivals
and customer feedback
Yunan Liu
a , ∗, Ward Whitt b
a Department of Industrial Engineering, North Carolina State University, Raleigh, NC 27695-7906, United States b Department of Industrial Engineering and Operations Research, Columbia University, New York, NY 10027-6699, United States
a r t i c l e i n f o
Article history:
Received 1 June 2015
Accepted 12 July 2016
Available online xxx
Keywords:
Queueing
Staffing algorithms for service systems
Time-varying arrival rates
Queues with feedback
Stabilizing performance
a b s t r a c t
Analytical offered-load and modified-offered-load (MOL) approximations are developed to determine
staffing levels that stabilize performance at designated targets in a non-Markovian many-server queueing
model with time-varying arrival rates, customer abandonment from queue and random feedback with
additional feedback delay in an infinite-server or finite-server queue. To provide a flexible model that can
be readily fit to system data, the model has Bernoulli routing, where the feedback probabilities, service-
time, patience-time and feedback-delay distributions all are general and may depend on the visit number.
Simulation experiments confirm that the new MOL approximations are effective. A many-server heavy-
traffic FWLLN shows that the performance targets are achieved asymptotically as the scale increases.
2 Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
Fig. 1. The (M t / { GI , GI } /s t + { GI , GI } ) + (GI/ ∞ ) model with delayed customer feedback and its Delayed Infinite-Server (DIS) approximation. The approximating offered load
is m (t) = m 1 (t) + m 2 (t) ≡ E[ B 1 (t)] + E[ B 2 (t)] .
p
s
c
m
d
D
d
a
g
m
a
d
g
(
t
w
p
t
c
h
t
t
c
W
b
c
t
t
T
t
w
S
w
c
d
s
s
a
{
focus on call center customers that may return later because
the initial service was unsatisfactory. Second, Yom-Tov and
Mandelbaum (2014) focus on the treatment of patients by a doc-
tor in a hospital that may naturally occur in stages, starting with
an initial screening and continuing later after tests have been or-
dered and completed. Our paper is closely related to Yom-Tov and
Mandelbaum (2014) , where a modified-offered-load (MOL) approx-
imation was proposed to help set staffing levels at a queue with
time-varying arrival rates and Markovian feedback after a delay in
an infinite-server (IS) queue. They showed that the MOL approxi-
mation has great potential for improved performance analysis in
healthcare, where the service times tend to be relatively long, so
that PSA does not apply.
Motivated by these applications, we consider a feedback model
that has appealing flexibility. In particular, instead of the Marko-
vian routing with fixed feedback probability p and one fixed
service-time distribution considered in Yom-Tov and Mandelbaum
(2014) , we consider history-dependent Bernoulli routing, where
there may be any number of visits and the feedback probability
p and the service-time distribution and the subsequent delay dis-
tribution (before returning for a new service) all may vary with
the visit number. We focus on the common important case of at
most one feedback, which seems to be a more realistic model than
Markovian routing, which produces a geometric random number
of feedbacks. It is significant that the approach here also extends
directly to any finite number of feedbacks; we demonstrate by also
considering examples with two feedback opportunities. Our meth-
ods also extend directly to time-dependent feedback probabilities,
but we do not examine that here. (The justification is that a time-
dependent independent thinning of an NHPP is again an NHPP; see
Sections 2.3 and 2.4 of Ross (1996) .)
We also allow customer abandonment, which often tends to be
more realistic for many service systems, as observed by Garnett,
Mandelbaum, and Reiman (2002) . The patience-time distributions
are also allowed to be non-exponential and depend on the visit
number. Just as in Yom-Tov and Mandelbaum (2014) , we use the
general offered-load (OL) method with the MOL refinement, as
reviewed in Jennings et al. (1996) , Green et al. (2007) , Liu and
Whitt (2012c) and Whitt (2013) . There are difference between the
MOL methods designed to stabilize the delay probability and the
abandonment probability, as discussed in Liu and Whitt (2012c) ,
but the main contribution here beyond Yom-Tov and Mandelbaum
(2014) is the new method for computing the time-varying of-
fered load. Because the offered load is the primary determinant of
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
erformance, the performance impact from more faithfully repre-
enting the service and feedback process in a time-varying setting
an be great.
To analyze this new feedback model with customer abandon-
ent, we draw on Liu and Whitt (2012c) in which we developed a
elayed-infinite-server (DIS) offered-load approximation and a new
IS-MOL ( DIS-modified-offered-load ) algorithm to determine time-
ependent staffing levels in order to stabilize expected delays and
bandonment probabilities at specified quality of service (QoS) tar-
ets in a many-server queue with time-varying arrival rates. The
odel in Liu and Whitt (2012c) was M t /GI /s t + GI model, having
rrivals according to an NHPP with arrival rate function λ( t ), in-
ependent and identically distributed (i.i.d.) service times with a
eneral distribution (the first GI ), a time-varying number of servers
the s t , to be determined), i.i.d. patience times with a general dis-
ribution (times to abandon from queue, the final + GI), unlimited
aiting space and the first-come first-served (FCFS) service disci-
line. We included non-exponential service and patience distribu-
ions as well as time-varying arrivals because they commonly oc-
ur; e.g. see Armony et al. (2015) and Brown et al. (2005) .
We refer to the base model with a single feedback considered
ere as (M t / { GI , GI } /s t + { GI , GI } ) + (GI/ ∞ ) . The main queue has
he two service-time cdf’s G i and patience cdf’s F i , depending on
he visit number, while the orbit queue has a single service-time
df H , with all waiting customers entering service in a FCFS order.
e develop approximations for the number of customers waiting
efore service and in service upon each visit and the number of
ustomers in orbit. When we refer to the number of customers in
he system or the waiting time, we do not include the orbit queue.
We also consider the associated (M t / { GI , GI } /s t + { GI , GI } ) +(GI /s t + GI ) model in which the orbit queue has finite capacity; in
hat case, it also has a staffing function and a patience distribution.
he goal is to stabilize expected potential waiting times (the vir-
ual waiting time before starting service on any visit of an arrival
ith infinite patience) at a fixed value w for all time and i = 1 , 2 .
ince these models are special kinds of two-class queueing models,
e also consider the more elementary ∑ 2
i =1 (M t /GI + GI) /s t two-
lass queue, in which the two classes arrive according to two in-
ependent NHPPs with arrival rate functions λ( i ) ( t ) and their own
ervice-time cdf’s G i and patience cdf’s F i , i = 1 , 2 , but there is a
ingle service facility with a time-varying number of servers s ( t ),
gain to be determined.
The approximating DIS model for the (M t / { GI , GI } /s t + GI , GI } ) + (GI/ ∞ ) feedback queue has five IS queues in
service system with time-varying arrivals and customer feedback,
Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14 7
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
0 2 4 6 8 10 12 14 16 18 20
90
100
110A
rriv
al r
ate
0 2 4 6 8 10 12 14 16 18 200
20
40
Exp
ecte
dqu
eue
leng
th
0 2 4 6 8 10 12 14 16 18 2010
15
20
Exp
ecte
dor
bit n
umbe
r
0 2 4 6 8 10 12 14 16 18 200
0.1
0.2
Aba
ndon
men
tpr
obab
ility
(ne
w)
0 2 4 6 8 10 12 14 16 18 200
0.2
0.4
Aba
ndon
men
tpr
obab
ility
(ol
d)
0 2 4 6 8 10 12 14 16 18 20
0.8
0.9
1
Del
aypr
obab
ility
0 2 4 6 8 10 12 14 16 18 200
0.1
0.2
0.3
0.4
0.5
Exp
ecte
dde
lay
0 2 4 6 8 10 12 14 16 18 200
50
100
150
Sta
ffing
Time
Fig. 2. Performance functions in the (M t (0 . 2) / { H 2 (1 , 4) , H 2 (5 , 4) } /s t + { M(2) , M(1) } ) + (0 . 2 , H 2 (1 , 4) / ∞ ) model with the sinusoidal arrival rate in (19) for λ̄ = 100 and r =
0 . 2 , Bernoulli feedback with probability p = 0 . 2 and an IS orbit queue: four cases of high waiting-time (low QoS) targets ( w = 0 . 10 , 0.20, 0.30 and 0.40) and simple DIS
staffing.
c
M
M
t
b
l
a
i
v
f
t
5
i
i
t
ustomers be exponential, but with different means, denoted by
( m ). In particular, we consider the (M t (r) /H 2 (1 , 4) , H 2 (5 , 4) /s t + (2) , M (1)) + (p, H 2 (1 , 4) / ∞ ) model with r = p = 0 . 2 . All service-
ime distributions are H 2 , while all patience distributions are M ,
ut the means vary, so that the complex refined DIS-MOL formu-
as in Section 4 associated with the aggregate model are needed,
nd are tested in these experiments. We also consider correspond-
ng models with non-exponential patience cdf’s in the and larger
alues of r in the appendix. The same stable performance is seen
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
or r = 0 . 5 , but some degradation in performance is seen where
he staffing decreases for r = 0 . 8 .
.2. Results from the simulation experiment
We simulated the model above starting empty over the time
nterval [0, 20]. We estimated the performance functions by tak-
ng averages from 20 0 0 independent replications. (Additional de-
ails are given in the online appendix.)
service system with time-varying arrivals and customer feedback,
8 Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
0 2 4 6 8 10 12 14 16 18 20
90
100
110A
rriv
al r
ate
0 2 4 6 8 10 12 14 16 18 200
2
4
6
Exp
ecte
dqu
eue
leng
th
0 2 4 6 8 10 12 14 16 18 20
14
16
18
20
22
Exp
ecte
dor
bit n
umbe
r
0 2 4 6 8 10 12 14 16 18 200
0.01
0.02
0.03
Aba
ndon
men
tpr
obab
ility
(ne
w)
0 2 4 6 8 10 12 14 16 18 200
0.02
0.04
Aba
ndon
men
tpr
obab
ility
(ol
d)
0 2 4 6 8 10 12 14 16 18 200
0.5
1
Del
aypr
obab
ility
0 2 4 6 8 10 12 14 16 18 200
0.010.020.030.040.05
Exp
ecte
dde
lay
0 2 4 6 8 10 12 14 16 18 200
100
200
Sta
ffing
Time
Fig. 3. Performance functions in the (M t (0 . 2) / { H 2 (1 , 4) , H 2 (5 , 4) } /s t + { M(2) , M(1) } ) + (0 . 2 , H 2 (1 , 4) / ∞ ) model with the sinusoidal arrival rate in (19) for λ̄ = 100 and r =
0 . 2 , Bernoulli feedback with probability p = 0 . 2 and an IS orbit queue: four cases of low waiting-time (high QoS) targets ( w = 0 . 01 , 0.02, 0.03 and 0.04) and DIS-MOL
staffing.
a
t
a
S
a
1
w
w
(
Figs. 2 and 3 show the results of the simulation experiment
for high and low waiting-time targets. In Fig. 2 the waiting-
time targets are w = 0 . 10 , 0 . 20 , 0 . 30 , 0 . 40 , so that the simple DIS
staffing is used, while in Fig. 3 the waiting-time targets are w =0 . 01 , 0 . 02 , 0 . 03 , 0 . 04 , ten times smaller, so that the refined DIS-
MOL staffing is used. The performance functions are averages
based on 20 0 0 independent replications.
Consistent with Liu and Whitt (2012c) and the FWLLN in
Section 3 , with the higher waiting-time targets in Fig. 2 we see
very smooth and accurate plots of the expected waiting times and
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
bandonment probabilities, which are the performance functions
o be stabilized, but strongly fluctuating expected queue lengths
nd delay probabilities, which agree closely with the formulas in
ection 2 . With the higher waiting-time targets, there is higher
bandonment probability, so that the maximum staffing is about
60 instead of about 100 + 100 = 200 in Fig. 3 with the lower
aiting-time targets. There is greater variability with the lower
aiting-time targets.
Fig. 3 shows that, consistent with experience in Feldman et al.
2008) and Liu and Whitt (2012c) , all performance functions tend
service system with time-varying arrivals and customer feedback,
Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14 9
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
2 4 6 8 10 12 14 16 18 200
0.5
1
1.5
2
2.5
Em
piric
al Q
oS β
(t)
Time
w = 0.0025
w = 0.005
w = 0.01
w = 0.02
w = 0.04w = 0.03
w = 0.06
Fig. 4. The empirical Quality of Service (QoS) provided by the DIS-MOL staffing in the (M t (0 . 2) / { H 2 (1 , 4) , H 2 (5 , 4) } /s t + { M(2) , M(1) } ) + (0 . 2 , H 2 (1 , 4) / ∞ ) example of
Fig. 3 as a function of the waiting-time target w .
0 2 4 6 8 10 12 14 16 18 20
Arr
ival
rat
e
90
100
110
Class 1
0 2 4 6 8 10 12 14 16 18 20
Arr
ival
rat
e
50
55
60
65
70
Class 2
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d
queu
e le
ngth
0
10
20
30
40
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d
queu
e le
ngth
0
5
10
15
20
25
0 2 4 6 8 10 12 14 16 18 20
Aba
ndon
men
tpr
obab
ility
0
0.05
0.1
0.15
0.2
0 2 4 6 8 10 12 14 16 18 20
Aba
ndon
men
tpr
obab
ility
0
0.1
0.2
0.3
0.4
0 2 4 6 8 10 12 14 16 18 20
Del
ay p
roba
bilit
y
0.75
0.8
0.85
0.9
0.95
1
0 2 4 6 8 10 12 14 16 18 20
Del
ay p
roba
bilit
y
0.75
0.8
0.85
0.9
0.95
1
Time0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
dde
lay
0
0.1
0.2
0.3
0.4
0.5
Time0 2 4 6 8 10 12 14 16 18 20
Sta
ffing
0
50
100
150
200
Fig. 5. Performance functions in the ∑ 2
i =1 (M t /H 2 (m i , 4) + M(m i ) /s t two-class model with the two sinusoidal arrival-rate functions in (22) , service-time means m 1 = 1 . 0 and
m 2 = 0 . 6 and patience means m 1 = 2 . 0 and m 2 = 1 . 0 : four cases of identical high waiting-time (low QoS) targets ( w = 0 . 10 , 0.20, 0.30 and 0.40) and simple DIS staffing at
both queues.
t
g
l
d
a
1
a
5
S
t
s
c
o be stabilized simultaneously with the lower waiting-time tar-
ets, after an initial startup effect due to starting empty. The de-
ay probability starts at 1 because the stabilizing staffing algorithm
oes not start staffing until time w > 0 . That feature ensures that
ll arrivals wait exactly w in the limiting fluid model (see Section
0 of Liu & Whitt (2012a) ), but it would probably not be used in
pplications.
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
.3. Square root staffing
We emphasize that the DIS OL m ( t ) given explicitly in
ection 2 is the key quantity being computed. The DIS OL quan-
ifies the essential demand, combining the impact of the random
ervice times with the time-varying arrival rate, both of which are
omplicated by the feedback. The relatively complicated DIS-MOL
service system with time-varying arrivals and customer feedback,
10 Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
0 2 4 6 8 10 12 14 16 18 20
Arr
ival
rat
e
90
100
110
Class 1
0 2 4 6 8 10 12 14 16 18 20
Arr
ival
rat
e
50
55
60
65
70
Class 2
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d
queu
e le
ngth
0
2
4
6
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d
queu
e le
ngth
0
1
2
3
0 2 4 6 8 10 12 14 16 18 20
Aba
ndon
men
tpr
obab
ility
0
0.01
0.02
0.03
0 2 4 6 8 10 12 14 16 18 20A
band
onm
ent
prob
abili
ty
0
0.01
0.02
0.03
0.04
0.05
0 2 4 6 8 10 12 14 16 18 20
Del
ay p
roba
bilit
y
0
0.2
0.4
0.6
0 2 4 6 8 10 12 14 16 18 20
Del
ay p
roba
bilit
y
0
0.2
0.4
0.6
Time0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
dde
lay
0
0.01
0.02
0.03
0.04
0.05
Time0 2 4 6 8 10 12 14 16 18 20
Sta
ffing
0
50
100
150
200
Fig. 6. Performance functions in the ∑ 2
i =1 (M t /H 2 (m i , 4) + M(m i ) /s t two-class model with the two sinusoidal arrival-rate functions in (22) , service-time means m 1 = 1 . 0 and
m 2 = 0 . 6 and patience means m 1 = 2 . 0 and m 2 = 1 . 0 : four cases of identical low waiting-time (high QoS) targets ( w = 0 . 01 , 0.02, 0.03 and 0.04) and DIS-MOL staffing at
both queues.
o
f
t
b
D
t
6
i
c
d
p
t
6
t
staffing, which requires an algorithm for computing an approxima-
tion for the steady-state performance in the stationary M/GI/s + GI
model, is of course also important in identifying the exact staffing
level required to stabilize the expected potential waiting times at
the target w . However, except for the specific QoS parameter β , the
same goal could be achieved by applying the simple square root
staffing (SRS) formula
s (t) ≡ m (t) + β√
m (t) , (20)
with this DIS OL m ( t ). Without the DIS-MOL step, we could just
search for the appropriate constant β to use in the SRS formula.
The DIS OL already succeeds in eliminating the dependence on
time.
As in Feldman et al. (2008) , we demonstrate the importance of
the DIS OL in the present context by plotting the implied empirical
QoS,
βDIS−MOL (t) =
s DISMOL (t) − m (t) √
m ( t) (21)
for the example considered in Fig. 3 . Fig. 4 shows that the DIS-MOL
staffing is indeed equivalent to SRS staffing for an appropriate QoS
parameter β , which is given on the y axis on the left, as a function
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
f the target w on the right. We present similar empirical QoS plots
or other examples in the online appendix.
The DIS OL is appropriate for smaller models as well, but then
he actual staffing and the resulting performance are complicated
ecause the discretization becomes very important. However, the
IS OL remains an important first step to identify the effective
ime-dependent demand.
. Other models
In this section we discuss the other two models mentioned
n the introduction. We first discuss the ∑ 2
i =1 (M t /GI + GI) /s t two-
lass queue, in which the two classes arrive according to two in-
ependent NHPPs. We then discuss the (M t / { GI , GI } /s t + { GI , GI } ) +(GI /s t + GI ) feedback model in which the orbit queue has finite ca-
acity. Afterwards, we discuss the model with two feedback oppor-
unities. More examples are discussed in the online appendix.
.1. Two-class queue
In this section we consider the associated
∑ 2 i =1 (M t /GI + GI) /s t
wo-class queue, in particular, the ∑ 2
i =1 (M t /H 2 (m i , 4) + M(m i ) /s t
service system with time-varying arrivals and customer feedback,
Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14 11
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
0 2 4 6 8 10 12 14 16 18 20
Arr
ival
rat
e
90
100
110
Main Queue
0 2 4 6 8 10 12 14 16 18 20
Orb
it
ar
rival
rat
e
0
50
Orbit Queue
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d
queu
e le
ngth
0
20
40
60
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d or
bit
queu
e le
ngth
0
10
20
0 2 4 6 8 10 12 14 16 18 20Aba
ndon
men
t
pr
obab
ility
(ne
w)
0
0.1
0.2
Time0 2 4 6 8 10 12 14 16 18 20A
band
onm
ent
prob
abili
ty (
old)
0
0.2
0.4Time
0 2 4 6 8 10 12 14 16 18 20Aba
ndon
men
t
pr
obab
ility
(or
bit)
0
0.2
0.4
0 2 4 6 8 10 12 14 16 18 20
Del
ay
pr
obab
ility
0.8
0.9
1
0 2 4 6 8 10 12 14 16 18 20
Orb
it D
elay
prob
abili
ty
0.8
0.9
1
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
dde
lay
0
0.2
0.4
0 2 4 6 8 10 12 14 16 18 20
Orb
it
Exp
ecte
d de
lay
0
0.2
0.4
0 2 4 6 8 10 12 14 16 18 20
Sta
ffing
0
50
100
150
0 2 4 6 8 10 12 14 16 18 20
Orb
it S
taffi
ng
0
20
40
60
Fig. 7. Performance functions in the (M t (0 . 2) / { H 2 (1 , 4) , H 2 (10 / 6 , 4) } /s t + { M(2) , M(1) } ) + (0 . 6 , H 2 (1 , 4) /s t + M(1)) model with the sinusoidal arrival rate in (19) for λ̄ = 100
and r = 0 . 2 , Bernoulli feedback with probability p = 0 . 6 and a finite-capacity orbit queue: four cases of identical high waiting-time (low QoS) targets ( w = 0 . 10 , 0.20, 0.30
and 0.40) and simple DIS staffing at both queues.
m
1
m
d
p
λ
λ
T
o
p
O
t
t
o
b
6
6
{
r
s
b
)
m
t
b
l
F
h
2
6
i
h
o
c
m
t
w
f
odel with H 2 ( m , 4) service-time cdf’s for both classes with m 1 = . 0 and m 2 = 0 . 6 and M ( m ) patience cdf’s for both classes with
1 = 2 . 0 and m 2 = 1 . 0 . We let the arrival processes be indepen-
ent NHPPs, but with different sinusoidal arrival-rate functions, in
articular,
1 (t) = 100(1 + 0 . 2 sin (t)) , and
2 (t) = 60(1 + 0 . 2 sin (0 . 8 t + 2)) . (22)
he analysis of this model is more elementary. First, there is no
rbit queue. We get the DIS OL by simply applying the DIS ap-
roximation to the two classes separately. That yields the per-class
L’s m i (t) = E[ B i (t)] for i = 1 , 2 and then we add to get the to-
al OL: m (t) = m 1 (t) + m 2 (t) . Given this overall DIS OL, we apply
he same refined DIS-MOL approximation in Section 4 . The results
f simulation experiments for high and low waiting-time targets,
ased on 20 0 0 independent replications, are shown in Figs. 5 and
. The results are good, just as in Section 5 .
.2. A finite-capacity orbit queue
In this section we consider the associated (M t / { GI , GI } /s t + GI , GI } ) + (GI/s t + GI) model with Bernoulli feedback after a
andom delay in a finite-capacity orbit queue. We use the
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
ame waiting-time targets to set the staffing levels in the or-
it queue and the main queue. In particular, we consider the
(M t (r) /H 2 (1 , 4) , H 2 (10 / 6 , 4) /s t + M(2) , M(1)) + (p, H 2 (1 , 4) /s t + M(1)
odel with r = 0 . 2 and p = 0 . 6 . Just as in Section 5 , all service-
ime distributions are H 2 , while all patience distributions are M ,
ut the means vary, so that the complex refined DIS-MOL formu-
as in Section 4 associated with the aggregate model are needed.
igs. 7 and 8 show the results of the simulation experiment for
igh and low waiting-time targets, respectively, again based on
0 0 0 independent replications, each starting empty.
.3. Two feedback opportunities
In this section we consider a modification of the base model
n which there are two feedback opportunities. Each customer that
as been fed back once returns again with probability p 2 after an-
ther delay in an IS orbit queue with cdf H 2 . Upon return, these
ustomers have service cdf G 3 and patience cdf F 3 . The new DIS
odel has eight IS queues in series, as depicted in Fig. 9 .
Since there are now three customer classes, characterized by
heir class-dependent service-time and patience-time distributions,
e easily generalize results in Theorem 1 to include the formulas
or class 3. We have
service system with time-varying arrivals and customer feedback,
12 Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
0 2 4 6 8 10 12 14 16 18 20
Arr
ival
rat
e
90
100
110
Main Queue
0 2 4 6 8 10 12 14 16 18 20
Orb
it
ar
rival
rat
e
0
20
40
60
80Orbit Queue
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d
queu
e le
ngth
0
2
4
6
8
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
d or
bit
queu
e le
ngth
0
1
2
3
0 2 4 6 8 10 12 14 16 18 20Aba
ndon
men
t
pr
obab
ility
(ne
w)
0
0.01
0.02
0.03
Time0 2 4 6 8 10 12 14 16 18 20A
band
onm
ent
prob
abili
ty (
old)
0
0.02
0.04Time
0 2 4 6 8 10 12 14 16 18 20Aba
ndon
men
t
pr
obab
ility
(or
bit)
0
0.01
0.02
0.03
0 2 4 6 8 10 12 14 16 18 20
Del
ay
pr
obab
ility
0
0.2
0.4
0.6
0 2 4 6 8 10 12 14 16 18 20
Orb
it D
elay
prob
abili
ty
0
0.2
0.4
0.6
0 2 4 6 8 10 12 14 16 18 20
Exp
ecte
dde
lay
0
0.02
0.04
0 2 4 6 8 10 12 14 16 18 20
Orb
it
Exp
ecte
d de
lay
0
0.02
0.04
0.06
0 2 4 6 8 10 12 14 16 18 20
Sta
ffing
0
100
200
0 2 4 6 8 10 12 14 16 18 20
Orb
it
Sta
ffing
0
20
40
60
Fig. 8. Performance functions in the (M t (0 . 2) / { H 2 (1 , 4) , H 2 (10 / 6 , 4) } /s t + { M(2) , M(1) } ) + (0 . 6 , H 2 (1 , 4) /s t + M(1)) model with the sinusoidal arrival rate in (19) for λ̄ = 100
and r = 0 . 2 , Bernoulli feedback with probability p = 0 . 6 and an IS orbit queue: four cases of low waiting-time (high QoS) targets ( w = 0 . 01 , 0.02, 0.03 and 0.04) and DIS-MOL
staffing.
Fig. 9. The DIS approximation for the (M t / { GI , GI , GI } /s t + { GI , GI , GI } ) + (GI/ ∞ ) + (GI/ ∞ ) model with two delayed customer feedback opportunities. Here there are two IS
orbit queues. The approximating offered load is m (t) = m 1 (t) + m 2 (t) + m 3 (t) ≡ E[ B 1 (t)] + E[ B 2 (t)] + E[ B 3 (t)] .
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a service system with time-varying arrivals and customer feedback,
European Journal of Operational Research (2016), http://dx.doi.org/10.1016/j.ejor.2016.07.018
DIS-MOL) approximation to set staffing levels with low waiting-
ime (high QoS) targets. We showed that we can use either the ag-
regate abandonment probability target or the waiting-time target,
ut the waiting-time target tends to produce a faster algorithm,
n part because the abandonment probability target F mol (w ; t) is a
ime-dependent function. We have presented results of simulation
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
xperiments in Sections 5 and 6 showing that the new DIS and
IS-MOL staffing algorithms are effective across a wide range of
oS targets.
The queue with Bernoulli feedback after an additional delay in a
nite-capacity orbit queue is a special case of a network of many-
erver queues with feedback. Our excellent results in this case in-
icate that the methods should apply to more general networks of
ueues, including multiple queues and customer classes, with var-
ous forms of routing, including models with retrials from blocked
rrivals as in the large literature reported in Artalejo (2010) , but
uch more general models remain to be examined carefully.
cknowledgments
We thank the two referees for their valuable comments. We
hank the National Science Foundation for support: NSF grants
MMI 1362310 (first author) and CMMI 10 6 6372 and 1265070 (sec-
nd author). This research began as part of the first author’s doc-
oral dissertation at Columbia University.
upplementary material
Supplementary material associated with this article can be
ound, in the online version, at 10.1016/j.ejor.2016.07.018 .
eferences
rmony, M. , Israelit, S. , Mandelbaum, A. , Marmor, Y. , Tseytlin, Y. , & Yom–
Tov, G. (2015). Patient flow in hospitals: A data-based queueing-science per-spective. Stochastic Systems, 5 , 146–194 .
rtalejo, J. R. (2010). Accessible bibliography on retrial queues: Progress in20 0 0–20 09. Mathematics of Computer Modeling, 51 , 1071–1081 .
rown, L. , Gans, N. , Mandelbaum, A. , Sakov, A. , Shen, H. , Zeltyn, S. , et al. (2005).
Statistical analysis of a telephone call center: A queueing-science perspective.Journal of American Statistical Association, 100 , 36–50 .
efraeye, M. , & van Nieuwenhuyse, I. (2013). Controlling excessive waiting timesin small service systems with time-varying demand: An extension of the ISA
algorithm. Decision Support Systems, 54 , 1558–1567 . ick, S. G. , Massey, W. A. , & Whitt, W. (1993). The physics of the M t / G / ∞ queue.
Operations Research, 41 , 731–742 .
eldman, Z. , Mandelbaum, A. , Massey, W. A. , & Whitt, W. (2008). Staffing of time–varying queues to achieve time-stable performance. Management Science, 54 ,
324–338 . arnett, O. , Mandelbaum, A. , & Reiman, M. I. (2002). Designing a call center
with impatient customers. Manufacturing and Service Operations Management, 4 ,208–227 .
reen, L. V. , Kolesar, P. J. , & Whitt, W. (2007). Coping with time-varying demand
when setting staffing requirements for a service system. Production and Opera-tions Management, 16 , 13–29 .
e, B., Liu, Y., & Whitt, W. (2016). Staffing a service system with non-Poisson non-stationary arrivals. Probability in the Engineering and Informational Sciences,
doi: 10.1017/S02699648160 0 019X . ennings, O. B. , Mandelbaum, A. , Massey, W. A. , & Whitt, W. (1996). Server staffing
to meet time-varying demand. Management Science, 42 , 1383–1394 .
iu, Y. , & Whitt, W. (2012a). The G t /GI /s t + GI many-server fluid queue. QueueingSystems, 71 , 405–4 4 4 .
iu, Y. , & Whitt, W. (2012b). A many-server fluid limit for the G t /GI /s t + GI queue-ing model experiencing periods of overloading. Operations Research Letters, 40 ,
307–312 . iu, Y. , & Whitt, W. (2012c). Stabilizing customer abandonment in many-server
queues with time-varying arrivals. Operations Research, 60 , 1551–1564 .
iu, Y. , & Whitt, W. (2014). Stabilizing performance in networks of queues withtime-varying arrival rates. Probability in the Engineering and Informational Sci-
ences, 28 , 419–449 . iu, Y., & Whitt, W. (2015). Online appendix: Stabilizing performance in many-server
queues with time-varying arrivals and customer feedback. Columbia University,http://www.columbia.edu/ ∼ww2040/StabFeedbackapp.pdf .
assey, W. A. , & Whitt, W. (1993). Networks of infinite-server queues with nonsta-tionary Poisson input. Queueing Systems, 13 , 183–250 .
oss, S. M. (1996). Stochastic processes (2nd). New York: Wiley .
tolletz, R. (2008). Approximation of the nonstationary M ( t )/ M ( t )/ c ( t ) queue usingstationary models: the stationary backlog-carryover approach. European Journal
of Operations Research, 190 , 478–493 . e Vericourt, F. , & Zhou, Y.-P. (2005). Managing response time in a call-routing prob-
lem with service failure. Operations Research, 53 , 968–981 .
service system with time-varying arrivals and customer feedback,
14 Y. Liu, W. Whitt / European Journal of Operational Research 0 0 0 (2016) 1–14
ARTICLE IN PRESS
JID: EOR [m5G; August 20, 2016;17:21 ]
W
Y
Whitt, W. (1982). Approximating a point process by a renewal process: Two basicmethods. Operations Research, 30 , 125–147 .
Whitt, W. (2002). Stochastic-process limits . New York: Springer . Whitt, W. (2005). Engineering solution of a basic call-center model. Management
Science, 51 , 221–235 .
Please cite this article as: Y. Liu, W. Whitt, Stabilizing performance in a
European Journal of Operational Research (2016), http://dx.doi.org/10.10
hitt, W. (2013). Offered load analysis for staffing. Manufacturing and Service Oper-ations Management, 15 , 166–169 .
om-Tov, G. , & Mandelbaum, A. (2014). Erlang R: A time-varying queue with reen-trant customers, in support of healthcare staffing. Manufacturing and Service Op-
erations Management, 16 , 283–299 .
service system with time-varying arrivals and customer feedback,