This lecture is expanded from:
|
|
- John Newton
- 5 years ago
- Views:
Transcription
1 This lecture is expanded from: HIGH VOLUME JOB SHOP SCHEDULING AND MULTICLASS QUEUING NETWORKS WITH INFINITE VIRTUAL BUFFERS INFORMS, MIAMI Nov 2, 2001 Gideon Weiss Haifa University (visiting MS&E, Stanford) Joint work with: Jim Dai, Michal Penn and students: Tammy Boudoukh, Yoni Nazarathy, Anat Kopzon.
2 Job Shop Scheduling Machines: i = 1, L, I Routes: r = 1, L, R Steps: (, ro), o= 1, L, K on machines: σr, 1, L, σr, K Processing times: X,, X r, 1 L r, K Classical Job Shop Problem: ONE job on each route, Schedule to minimize MAKESPAN: T OPT = Earliest time to complete all jobs NP-Hard: 10x10 problem took 10 years to solve. Machine Lower Bound: * * Ti = Xr, o, T = C = max T σ i i ro, = Job lower bound: X = X, X = max X r K * r, o o= 1 r r i MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
3 Step# sum This is the famous 10x10 job-shop problem, First Appeared in 1963 book of Muth and Thompson. Solution found by Carlier and Pinson, 1988, Solution is 930 Lower bound 631. job 1 job 2 job 3 job 4 job job 6 job 7 job 8 job 9 job 10 mach bound Machine lower bound is 631 (Machine 4). Job lower bound is 655. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
4 More About Job-Shops: The problem is NP-hard. In fact, Williamosn et al showed that unless P=NP there is no polynomial time algorithm that can construct a schedule that is within 5/4 of the optimum. Heuristics can solve Job-shops approximately. Pupular is the shifting bottleneck heuristic which can deal with 10x10 in a few seconds (it finds optimum in 20 seconds but without proof of optimality). These heuristics fail to functions for more then a few dozen jobs on 20 machines. Sevastyanov (see also Baranyi) has a polynomial heuristic which achieve a worst case schedule better than * 2 T + ( K 1)( IK + 2K 1) max Xrk, ( j) For the 10x10 problem this gives a gap of ~800,000. Recent heuirstics exploit division of operations to long/short/tiny MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
5 High Volume Job Shops - Identical N IDENTICAL JOBS on each route, processing times X,, X r, 1 L r, K Machine Lower Bound: * * Ti = N Xr, o, T = NC = max T σ i i ro, = Job lower bound: For job j on route r: X = X, X = max X r K * r, o o= 1 r r i MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
6 P i p e l i n i n g : (Tammy Boudoukh, Bertsimas & Gamarnik) Cycle 1: Do operation 1 on job 1 of each route. Cycle 2: Do operation 2 on job 1 and operation 1 on job 2 of each route. Cycle 3: Do operation 3 on job 1, operation 2 on job 2, and operation 1 on job 3 of each route. Cycle K 1: Do operation K 1 on job 1, operation K 2 on job 2,... etc... and operation 1 on job K 1 of each route. Pipeline is full: Do N K+1 Full Cycles. Runout: K 1 Partial Cycles to complete all jobs. * * OPT H * * NC = T T T NC + ( K 1) C H Gap: T * T = O() 1 MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
7 High Volume Job Shops - Similar N SIMILAR JOBS on each route, processing times X (), j L, X () j j =, K, N. r, 1 r, K 1 Average: m = E( X ( j)) or m ro, ro, T = X ( j), T = max T, i N r, o σ = i j= ro, 1 Bottleneck Machine: i = arg max m ro, * * i σ 1 = N X ro j ro,, () j i ro, = i MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
8 Cyclic Scheduling with Safety Stocks To keep Bottleneck machine busy throughout full cycles: Build up SAFETY STOCKS, S in each buffer. Then perform Cycles, initiated by the Bottleneck Machine. Bottleneck idles if maximal queue length exceeds S. Theorem (Dai & W): Xro, () j i.i.d. with means m ro, and (1) i * is a unique, single, bottleneck. (2) Xro, () j possess exponential moments. Let T *, T H be the (random) lower bound and heuristic makespans. There exist C1, C2 such that for all N: With safety stocks S = C 1 log N, H * PT ( T > C log N) < / N 2 1 MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
9 Multi-Class Queueing/Fluid Networks: Nodes (machines): i= 1, K, I Classes, queues, buffers, steps: k = 1, K, K Constituencies: i= σ( k), k C i Jobs j = 1, K, N Arrivals: Aj (), renewal process, rate α Processing: Xk (), j i.i.d. mean m k Routing: k l, Bernoulli, P kl, Qt () = Q( 0) + A() t ( I P ) S( T( t)) 0 ST ( ( t)) FLUID LIMIT UNDER FLUID SCALING OF MCQN Rescale time and space equally by initial state which qt ( ) lim 1 = n Qnt ( ), n = Q( 0 ) n FLUID RELAXATION of JOB-SHOP or MCQN Integrality of Jobs: qk () t need not be integer. 1 Integrality of Machines: Flow rate 0 uk () t. mk Integrality of Work: Part of Job completed at Buffer k can already be processed at Buffer k + 1. qt () = q( 0) + α t ( I P )() usds 0 mu() t 1, ut () 0 k C i k k t 0 MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
10 Minimum Fluid Job Shop Makespan The cyclic policy with safety stocks S serves N S( K 1) full cycles, in which one job is processed in each buffer, i.e each (, ro) This corresponds to u () t N/ T for all buffers. ro, Fluid N * Route R N Route 2 N Route 1 T * Time MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
11 MCFN Minimum Makespan Q ( 0) = q ( 0), k = 1, K, K k k i 1 i i * T = max Ti i + 1 * T { CM( I P ) ( q( 0) + αt)} Q ( 0) = ( I P ) ( Q( 0) + αt ) u () t = Q ( 0)/ T k + k q + t k () t = Q + k ( 0) 1 *, qk() t T Fluid * = t Qk ( 0 ) 1 * T K K-1 1 T * Time MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
12 Fluid Imitation Algorithms Base decision on the gap between fluid and actual queues. Compare the nominal qk + () t to the actual queue Qk + () t Give priority to k for which Q + k () t q + k () t is large. Equivalently: k t t k k 0 0 k Let q () t = q () s ds, Q () t = Q () s ds, Give priority to k for which q () t Q () t is large. k k qk (), t Qk () t are easy to keep track of: They are the cumulative amount that of fluid or of actual jobs that flowed out of buffer k. q Priority measure ρ k () t Q k () t k () t = q () t Cyclic Fluid Algorithm (CFA): Machine i works on buffer: * k = argmax ρ ( t) : k C i k i k { } Greedy Fluid Algorithm (GFA): Machine i works on buffer: * k = argmax ρ ( t) : k C, Q ( t) > 0 { } i k i k MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
13 Comparing the policies: It can be shown that, with safety stocks S, GFA is at least as good as CFA which is at least as good as Dai&W Cyclic algorithm. GFA can function with no Safety Stocks. For identical Jobs: MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
14 For Similar Jobs, GFA with no Safety Stocks does not show logarithmic gap! MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
15 Job Shop Simulation Study (Yoni Nazarathy) MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
16 Empirical Conclusions: (1) The gap does not grow like log( N ), but remains O( 1) (2) The distribution of the processing time has little effect (3) Gap of order O( 1 ) can occur with multiple bottlenecks Experiments with MT10x10 and other examples. N ranging from 1 to Enough replicates for clear evidence. Disttibution: CV=0, CV=0.25, CV=1.0, Weibull 1/2, Pareto 3, Pareto 2. MT10 Re-entrant: N jobs on a single, 100 step route. MT10 Balanced: All machines require 631 time per cycle MT10 Round-Bal: Each route starts on a different machine Differing Nr, r = 1, K, R Policies of: Random buffer, Random job, LBFS, FBFS performed much worse, with gaps up to ON ( ) MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
17 M T 1 0 : MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
18 MT10 (cont1): MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
19 MT10 (cont2) MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
20 MT10 Round-Balanced: MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
21 MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
22 O(1) Gap for Re-entrant Lines: Surprising: O( 1 ) gap says: Bottleneck machine never idles between some random time T starve independent of N and completion of its work When bottleneck machine completes its work, other machines have runout time T runout independent of N. In a re-entrant line, all N jobs follow the same route. In that case, q + () t = N( 1 t/ T * ), q () t = Nt/ T *, k k If Qk() t > 0 Qk () t < Qk 1 () t, k has priority over k 1 Hence, GFA is actually simply LBFS. We know that LBFS is stable, for input rate ρ < 1 MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
23 GFA LBFS Fluid Solution Look at the fluid solution under cyclic policy, and compare to that under LBFS: Cyclic Fluid Solution GFA LBFS Solution a (L) a (1) Algorithm Buffer a () 1 is first buffer of the bottleneck machine, * Slope is N/ T * Flow of N/ T keeps buffers a () 1 + 1, K, K empty. () 1 Fluid out of buffers 1, K, a 1is faster, which defines a secondary bottleneck machine. a ( 2) is first buffer of secondary bottleneck machine. Continue until a ( L) = 1 MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
24 ReEntrant QN with infinite 1st buffer Assume the first buffer has an unlimited supply of jobs. Consider the system under LBFS policy: Buffer 1 is served only if all the other buffers of machine σ 1 are empty. The fluid picture now is: Queue in buffer 1 is ( L 1) ( 1) Buffers a, K, a : ρ > 1 All other queues are stable No starvation after T starve Stable time for T runout N a (L) a (1) MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
25 Multiple Bottlenecks: We found that high volume Job Shops with more than one bottleneck may also be scheduled with a gap of O( 1 ). We also found that high volume job shop problems with random routes, where each route has K processing steps, and the processing times of all steps are i.i.d. (i.e. all machines are bottlenecks), can also be scheduled with a gap of O( 1 ). How is that possible? We now look at MCQN with Infinite Virtual Buffers. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
26 Queueing Systems v. Manufacturing Systems The Single Server Queue: λ As long as λ < µ, ρ = <, µ Server is idle 1 ρ of his time. 1 system will be stable. λ λ µ A Single Machine Manufacturing System: Infinite supply of work. There is a virtual queue. Machine works all the time. Virtual Queue µ µ A Two Machine Manufacturing System: First machine works all the time. Second machine is a single server queue. System λ stable if λ < µ, ρ = < 1, and second machine idles 1 ρ. µ Virtual Queue λ λ µ λ MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
27 Queueing Systems v. Manufacturing Systems Efficient Operation of the Two Machine Manufacturing System: To operate the system efficiently we add work for the second, faster machine: λ λ 1 Virtual 1 Queue λ 1 µ 1 λ 2 Virtual Queue Second machine is a single server queue with vacations: λ Stable, non-idling, production from bottom stream µ λ A Two Machine Manufacturing System with Two Production Routes: Consider Virtual Queue λ 1 µ 1 ν 1 ν 2 µ 2 λ 2 Virtual Queue We can operate this system to be stable and non-idling. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
28 The Push Pull System (Anat Kopzon): Consider the following Two Machine Two Routes System: 1 2 λ 1 µ 2 ν 1 ν 2 µ 1 λ 2 Top Route goes from machine 1 to machine 2. Bottom Route goes from machine 2 to machine 1. Machine i works feeds machine 3 i at rate λ i. Machine i serves its own queue at rate µ i. Flow balance: To keep machines busy all the time and be stable: α i proportion of time machine i is feeding, 1 α i it is serving. ν1 = α1λ1 = ( 1 α 2) µ 2 The flow ratessatisfy: ν2 = α2λ2 = ( 1 α 1) µ 1 Hence: λµ µ λ ν 1 2( 1 2) 1 = µµ λλ This be achieved by a stable policy! ν λ 2µ 1( µ 2 λ = 1) µµ λλ MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
29 The Push Pull System Operating Rules: Assume: λ < µ, λ < µ The policy for Machine i is: produce for 3-i rate λi serve queue i rate µi part leaves system is queue i empty? Empty Not empty send part to queue 3-i keep part produced for 3-i Lemma: After an initial period, exactly one of the queues is non-empty at any time. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
30 Steady State of M/M/ System: Assume: production and service time are momoryless - The system is Markovian: Define the following states: ia, both machines feeding, machine 1 has i in queue ia, machine 1 is serving its queue, and has i in queue ib, both machines feeding, machine 2 has i in queue ib, machine 2 is serving its queue, and has i in queue State transition rates are: λ λ 2 λ 2 λ 1 λ 1 λ 1 2 3,a 2,a 1,a 1,B 2,B 3,B µ 1 λ 1 λ 1 λ 1 µ 1 µ 1 3,A 2,A λ 2 λ 2 λ 2 1,A µ 2 µ 1 1,b µ 2 µ 2 λ 2 λ 2 λ 2 µ 2 2,b 3,b λ 1 λ 1 λ 1 These are solved by: pia, = ΠaP( X = i 1) X ~ Geom( λ1/ λ1+ λ2) pia, = Π AP( X+ Z = i 1) Z ~ Ge om( µ 1 λ2 / λ2) pib, = ΠbP( Y = i 1) Y ~ Geom( λ2 / λ1+ λ2) pib, = ΠBP( Y + W = i 1) W ~ Ge om( µ 2 λ1/ λ1) Here Πa, ΠA, Πb, ΠB are the times proportions of states aabb,,, MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
31 Cycles in the Operations of the System The sample path of the push pull system will cycle through 4 periods: Queue 1 nonempty Queue 2 nonempty x a A b B x exp(λ 1 ) Μ/Μ/1 λ 2, µ 1 exp(λ 2 ) Μ/Μ/1 λ 1, µ 2 x The times proportions of states aabb,,, Π Π Π Π a b A B λ = 2( µ 1 λ 2)( µ 2 λ1) ( λ1+ λ 2)( µ 1µ 2 λ1λ2) λ µ λ µ λ = 1( 1 2)( 2 1) ( λ1+ λ 2)( µ 1µ 2 λ1λ2) λ µ λ = 2( 2 1) ( µµ 1 2 λλ 1 2) λ µ λ = 1( 1 2) ( µµ λλ ) MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
32 Generalization to M/G/ Assume now that the service times are general i.i.d. with rates µ µ production time are still memoryless with rates λ, λ , 2 while the The cycles of operation at each machine can be regarded as an M/G/1 queue with vacation. λ arrivals 2 Server 1 is working Server 1 vacation - no arrivals x a A b B x exp(λ 1 ) Μ/G/1 λ 2, µ 1 exp(λ 2 ) Μ/G/1 λ 1, µ 2 x A small trick is needed: Replace the bb, period, where the queue of server 1 is zero by a single exp( λ 2 ) time which is the wait for first customer in the M/G/1 after end of busy period. Vacation Server 1 is working x exp(λ ) 2 x a A exp(λ ) 1 Μ/G/1 λ, µ 1 2 x We obtain explicit steady state from that. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
33 Generalization J machines Assume there are J machines. We assume that at any time at most one machine i is "serving" i.e. performing a second operation on a job, at rate µ i. All the other machines are feeding machine i, at a total rate Λ i = λ. When the queue of machine i is empty, it will send its previously produced kept from start of busy period job to machine j with probability p ij Machine j then becomes the serving machine. If Λ i < µ i this system is again stable. j i j µ3 Λ 1,C 2,C Λ 3,C Λ µ3 µ3 µ3 λ 3 λ 3 λ 3 1,c 2,c 3,c Λ Λ Λ Λ µ 1 1 Λ1 Λ ,a 2,a 1,a µ 1,B 2,B 3,B 2 µ 2 µ 2 µ 2 λ 1 λ 1 λ 1 µ 1 µ 1 3,A 2,A Λ Λ Λ ,A p p µ 1 p 21 p 13 p 32 λ 2 λ 2 λ 2 1,b p 12 Λ Λ Λ 2,b 3,b Λ Λ Λ MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
34 Conclusion - What have we done? We have analyzed a multiclass queueing system. The policy under which we operated the system broke the sample paths into very simple periods. As a result we were able to obtain explicit expression for the system in steady state. Surprise: This system seems to operate at ρ i = 1 for all the machines, and yet it is stable! MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
35 5 machine 20 buffer re-entrant example I II III IV V Starting from some initial fluid level, and minimizing total holding costs until the system is empty is an SCLP. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
36 S o l u t i o n : In this solution we use piecewise constant flow rates out of the buffers. In each interval several machines are fully utilized, so this has multiple bottlenecks. We hope that this can be operated by a stable fluid imitation heuristic. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
37 M e s s a g e : Manufacturing Systems, with Ample Supply of Work, under Fluid Control can: Work with Full Utilization (ρ=1) Be Stable (no heavy traffic) Optimize Costs. MS&E324, Stanford University, Spring Gideon Weiss manufacturing & control
A Semiconductor Wafer
M O T I V A T I O N Semi Conductor Wafer Fabs A Semiconductor Wafer Clean Oxidation PhotoLithography Photoresist Strip Ion Implantation or metal deosition Fabrication of a single oxide layer Etching MS&E324,
More informationGideon Weiss University of Haifa. Joint work with students: Anat Kopzon Yoni Nazarathy. Stanford University, MSE, February, 2009
Optimal Finite Horizon Control of Manufacturing Systems: Fluid Solution by SCLP (separated continuous LP) and Fluid Tracking using IVQs (infinite virtual queues) Stanford University, MSE, February, 29
More informationControl of multi-class queueing networks with infinite virtual queues
Control of multi-class queueing networks with infinite virtual queues Erjen Lefeber (TU/e) Workshop on Optimization, Scheduling and Queues Honoring Gideon Weiss on his Retirement June 8, 2012 Where innovation
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network. Haifa Statistics Seminar May 5, 2008
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy Gideon Weiss Haifa Statistics Seminar May 5, 2008 1 Outline 1 Preview of Results 2 Introduction Queueing
More informationA fluid approach to job shop scheduling: theory, software and experimentation
A fluid approach to job shop scheduling: theory, software and experimentation Yoni Nazarathy, TU/e Eindhoven, The Department of Mechanical Engineering and EURANDOM, Gideon Weiss, The University of Haifa,
More informationAdvanced Computer Networks Lecture 3. Models of Queuing
Advanced Computer Networks Lecture 3. Models of Queuing Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/13 Terminology of
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy a,1, Gideon Weiss a,1 a Department of Statistics, The University of Haifa, Mount Carmel 31905, Israel.
More informationA Push Pull Network with Infinite Supply of Work
DOI 10.1007/s11134-009-9121-1 A Push Pull Network with Infinite Supply of Work Anat Kopzon Yoni Nazarathy Gideon Weiss Received: 31 January 2007 / Revised: 28 July 2008 Springer Science+Business Media,
More information5 Lecture 5: Fluid Models
5 Lecture 5: Fluid Models Stability of fluid and stochastic processing networks Stability analysis of some fluid models Optimization of fluid networks. Separated continuous linear programming 5.1 Stability
More informationGlossary availability cellular manufacturing closed queueing network coefficient of variation (CV) conditional probability CONWIP
Glossary availability The long-run average fraction of time that the processor is available for processing jobs, denoted by a (p. 113). cellular manufacturing The concept of organizing the factory into
More informationLinear Model Predictive Control for Queueing Networks in Manufacturing and Road Traffic
Linear Model Predictive Control for ueueing Networks in Manufacturing and Road Traffic Yoni Nazarathy Swinburne University of Technology, Melbourne. Joint work with: Erjen Lefeber (manufacturing), Hai
More informationSolutions to Homework Discrete Stochastic Processes MIT, Spring 2011
Exercise 6.5: Solutions to Homework 0 6.262 Discrete Stochastic Processes MIT, Spring 20 Consider the Markov process illustrated below. The transitions are labelled by the rate q ij at which those transitions
More informationComputer Networks More general queuing systems
Computer Networks More general queuing systems Saad Mneimneh Computer Science Hunter College of CUNY New York M/G/ Introduction We now consider a queuing system where the customer service times have a
More informationOn the Resource/Performance Tradeoff in Large Scale Queueing Systems
On the Resource/Performance Tradeoff in Large Scale Queueing Systems David Gamarnik MIT Joint work with Patrick Eschenfeldt, John Tsitsiklis and Martin Zubeldia (MIT) High level comments High level comments
More informationCPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017
CPSC 531: System Modeling and Simulation Carey Williamson Department of Computer Science University of Calgary Fall 2017 Motivating Quote for Queueing Models Good things come to those who wait - poet/writer
More informationCDA5530: Performance Models of Computers and Networks. Chapter 4: Elementary Queuing Theory
CDA5530: Performance Models of Computers and Networks Chapter 4: Elementary Queuing Theory Definition Queuing system: a buffer (waiting room), service facility (one or more servers) a scheduling policy
More informationIntroduction to Markov Chains, Queuing Theory, and Network Performance
Introduction to Markov Chains, Queuing Theory, and Network Performance Marceau Coupechoux Telecom ParisTech, departement Informatique et Réseaux marceau.coupechoux@telecom-paristech.fr IT.2403 Modélisation
More informationPerformance Evaluation of Queuing Systems
Performance Evaluation of Queuing Systems Introduction to Queuing Systems System Performance Measures & Little s Law Equilibrium Solution of Birth-Death Processes Analysis of Single-Station Queuing Systems
More informationCHAPTER 4. Networks of queues. 1. Open networks Suppose that we have a network of queues as given in Figure 4.1. Arrivals
CHAPTER 4 Networks of queues. Open networks Suppose that we have a network of queues as given in Figure 4.. Arrivals Figure 4.. An open network can occur from outside of the network to any subset of nodes.
More informationBIRTH DEATH PROCESSES AND QUEUEING SYSTEMS
BIRTH DEATH PROCESSES AND QUEUEING SYSTEMS Andrea Bobbio Anno Accademico 999-2000 Queueing Systems 2 Notation for Queueing Systems /λ mean time between arrivals S = /µ ρ = λ/µ N mean service time traffic
More informationSlides 9: Queuing Models
Slides 9: Queuing Models Purpose Simulation is often used in the analysis of queuing models. A simple but typical queuing model is: Queuing models provide the analyst with a powerful tool for designing
More informationOn the static assignment to parallel servers
On the static assignment to parallel servers Ger Koole Vrije Universiteit Faculty of Mathematics and Computer Science De Boelelaan 1081a, 1081 HV Amsterdam The Netherlands Email: koole@cs.vu.nl, Url: www.cs.vu.nl/
More informationChapter 6 Queueing Models. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter 6 Queueing Models Banks, Carson, Nelson & Nicol Discrete-Event System Simulation Purpose Simulation is often used in the analysis of queueing models. A simple but typical queueing model: Queueing
More informationIntro Refresher Reversibility Open networks Closed networks Multiclass networks Other networks. Queuing Networks. Florence Perronnin
Queuing Networks Florence Perronnin Polytech Grenoble - UGA March 23, 27 F. Perronnin (UGA) Queuing Networks March 23, 27 / 46 Outline Introduction to Queuing Networks 2 Refresher: M/M/ queue 3 Reversibility
More informationNon Markovian Queues (contd.)
MODULE 7: RENEWAL PROCESSES 29 Lecture 5 Non Markovian Queues (contd) For the case where the service time is constant, V ar(b) = 0, then the P-K formula for M/D/ queue reduces to L s = ρ + ρ 2 2( ρ) where
More informationOPTIMAL CONTROL OF A FLEXIBLE SERVER
Adv. Appl. Prob. 36, 139 170 (2004) Printed in Northern Ireland Applied Probability Trust 2004 OPTIMAL CONTROL OF A FLEXIBLE SERVER HYUN-SOO AHN, University of California, Berkeley IZAK DUENYAS, University
More informationExercises Stochastic Performance Modelling. Hamilton Institute, Summer 2010
Exercises Stochastic Performance Modelling Hamilton Institute, Summer Instruction Exercise Let X be a non-negative random variable with E[X ]
More informationDynamic Control of Parallel-Server Systems
Dynamic Control of Parallel-Server Systems Jim Dai Georgia Institute of Technology Tolga Tezcan University of Illinois at Urbana-Champaign May 13, 2009 Jim Dai (Georgia Tech) Many-Server Asymptotic Optimality
More informationSince D has an exponential distribution, E[D] = 0.09 years. Since {A(t) : t 0} is a Poisson process with rate λ = 10, 000, A(0.
IEOR 46: Introduction to Operations Research: Stochastic Models Chapters 5-6 in Ross, Thursday, April, 4:5-5:35pm SOLUTIONS to Second Midterm Exam, Spring 9, Open Book: but only the Ross textbook, the
More informationSTABILITY OF MULTICLASS QUEUEING NETWORKS UNDER LONGEST-QUEUE AND LONGEST-DOMINATING-QUEUE SCHEDULING
Applied Probability Trust (7 May 2015) STABILITY OF MULTICLASS QUEUEING NETWORKS UNDER LONGEST-QUEUE AND LONGEST-DOMINATING-QUEUE SCHEDULING RAMTIN PEDARSANI and JEAN WALRAND, University of California,
More informationManagement of demand-driven production systems
Management of demand-driven production systems Mike Chen, Richard Dubrawski, and Sean Meyn November 4, 22 Abstract Control-synthesis techniques are developed for demand-driven production systems. The resulting
More informationChapter 2. Poisson Processes. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 2. Poisson Processes Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Outline Introduction to Poisson Processes Definition of arrival process Definition
More informationFluid Models of Parallel Service Systems under FCFS
Fluid Models of Parallel Service Systems under FCFS Hanqin Zhang Business School, National University of Singapore Joint work with Yuval Nov and Gideon Weiss from The University of Haifa, Israel Queueing
More informationQueuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem. Wade Trappe
Queuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem Wade Trappe Lecture Overview Network of Queues Introduction Queues in Tandem roduct Form Solutions Burke s Theorem What
More informationQuiz 1 EE 549 Wednesday, Feb. 27, 2008
UNIVERSITY OF SOUTHERN CALIFORNIA, SPRING 2008 1 Quiz 1 EE 549 Wednesday, Feb. 27, 2008 INSTRUCTIONS This quiz lasts for 85 minutes. This quiz is closed book and closed notes. No Calculators or laptops
More informationStability and Asymptotic Optimality of h-maxweight Policies
Stability and Asymptotic Optimality of h-maxweight Policies - A. Rybko, 2006 Sean Meyn Department of Electrical and Computer Engineering University of Illinois & the Coordinated Science Laboratory NSF
More informationFlow Shop and Job Shop Models
Outline DM87 SCHEDULING, TIMETABLING AND ROUTING Lecture 11 Flow Shop and Job Shop Models 1. Flow Shop 2. Job Shop Marco Chiarandini DM87 Scheduling, Timetabling and Routing 2 Outline Resume Permutation
More informationData analysis and stochastic modeling
Data analysis and stochastic modeling Lecture 7 An introduction to queueing theory Guillaume Gravier guillaume.gravier@irisa.fr with a lot of help from Paul Jensen s course http://www.me.utexas.edu/ jensen/ormm/instruction/powerpoint/or_models_09/14_queuing.ppt
More information5/15/18. Operations Research: An Introduction Hamdy A. Taha. Copyright 2011, 2007 by Pearson Education, Inc. All rights reserved.
The objective of queuing analysis is to offer a reasonably satisfactory service to waiting customers. Unlike the other tools of OR, queuing theory is not an optimization technique. Rather, it determines
More informationTutorial: Optimal Control of Queueing Networks
Department of Mathematics Tutorial: Optimal Control of Queueing Networks Mike Veatch Presented at INFORMS Austin November 7, 2010 1 Overview Network models MDP formulations: features, efficient formulations
More informationOperations Research Letters. Instability of FIFO in a simple queueing system with arbitrarily low loads
Operations Research Letters 37 (2009) 312 316 Contents lists available at ScienceDirect Operations Research Letters journal homepage: www.elsevier.com/locate/orl Instability of FIFO in a simple queueing
More informationSystems Simulation Chapter 6: Queuing Models
Systems Simulation Chapter 6: Queuing Models Fatih Cavdur fatihcavdur@uludag.edu.tr April 2, 2014 Introduction Introduction Simulation is often used in the analysis of queuing models. A simple but typical
More informationMarkov Processes and Queues
MIT 2.853/2.854 Introduction to Manufacturing Systems Markov Processes and Queues Stanley B. Gershwin Laboratory for Manufacturing and Productivity Massachusetts Institute of Technology Markov Processes
More informationBuzen s algorithm. Cyclic network Extension of Jackson networks
Outline Buzen s algorithm Mean value analysis for Jackson networks Cyclic network Extension of Jackson networks BCMP network 1 Marginal Distributions based on Buzen s algorithm With Buzen s algorithm,
More informationQueuing Theory. Using the Math. Management Science
Queuing Theory Using the Math 1 Markov Processes (Chains) A process consisting of a countable sequence of stages, that can be judged at each stage to fall into future states independent of how the process
More informationQueues and Queueing Networks
Queues and Queueing Networks Sanjay K. Bose Dept. of EEE, IITG Copyright 2015, Sanjay K. Bose 1 Introduction to Queueing Models and Queueing Analysis Copyright 2015, Sanjay K. Bose 2 Model of a Queue Arrivals
More informationIntroduction to queuing theory
Introduction to queuing theory Claude Rigault ENST claude.rigault@enst.fr Introduction to Queuing theory 1 Outline The problem The number of clients in a system The client process Delay processes Loss
More information57:022 Principles of Design II Final Exam Solutions - Spring 1997
57:022 Principles of Design II Final Exam Solutions - Spring 1997 Part: I II III IV V VI Total Possible Pts: 52 10 12 16 13 12 115 PART ONE Indicate "+" if True and "o" if False: + a. If a component's
More informationScheduling Multiclass Queueing Networks via Fluid Models
Scheduling Multiclass Queueing Networks via Fluid Models John Hasenbein OR/IE, UT-Austin Ron Billings OR/IE, UT-Austin Leon Lasdon MSIS, UT-Austin Gideon Weiss Statistics, Univ of Haifa www.me.utexas.edu/~has
More informationDynamic Matching Models
Dynamic Matching Models Ana Bušić Inria Paris - Rocquencourt CS Department of École normale supérieure joint work with Varun Gupta, Jean Mairesse and Sean Meyn 3rd Workshop on Cognition and Control January
More informationQueuing Theory. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011
Queuing Theory Richard Lockhart Simon Fraser University STAT 870 Summer 2011 Richard Lockhart (Simon Fraser University) Queuing Theory STAT 870 Summer 2011 1 / 15 Purposes of Today s Lecture Describe general
More information1 Ordinary Load Balancing
Comp 260: Advanced Algorithms Prof. Lenore Cowen Tufts University, Spring 208 Scribe: Emily Davis Lecture 8: Scheduling Ordinary Load Balancing Suppose we have a set of jobs each with their own finite
More informationMaximizing throughput in zero-buffer tandem lines with dedicated and flexible servers
Maximizing throughput in zero-buffer tandem lines with dedicated and flexible servers Mohammad H. Yarmand and Douglas G. Down Department of Computing and Software, McMaster University, Hamilton, ON, L8S
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationQueuing Analysis. Chapter Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall
Queuing Analysis Chapter 13 13-1 Chapter Topics Elements of Waiting Line Analysis The Single-Server Waiting Line System Undefined and Constant Service Times Finite Queue Length Finite Calling Problem The
More informationContinuous-time Markov Chains
Continuous-time Markov Chains Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 23, 2017
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationGI/M/1 and GI/M/m queuing systems
GI/M/1 and GI/M/m queuing systems Dmitri A. Moltchanov moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2716/ OUTLINE: GI/M/1 queuing system; Methods of analysis; Imbedded Markov chain approach; Waiting
More informationChapter 11. Approximation Algorithms. Slides by Kevin Wayne Pearson-Addison Wesley. All rights reserved.
Chapter 11 Approximation Algorithms Slides by Kevin Wayne. Copyright @ 2005 Pearson-Addison Wesley. All rights reserved. 1 Approximation Algorithms Q. Suppose I need to solve an NP-hard problem. What should
More informationLecture 7: Simulation of Markov Processes. Pasi Lassila Department of Communications and Networking
Lecture 7: Simulation of Markov Processes Pasi Lassila Department of Communications and Networking Contents Markov processes theory recap Elementary queuing models for data networks Simulation of Markov
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/6/16. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 24 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/6/16 2 / 24 Outline 1 Introduction 2 Queueing Notation 3 Transient
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/25/17. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 26 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/25/17 2 / 26 Outline 1 Introduction 2 Queueing Notation 3 Transient
More information1.225 Transportation Flow Systems Quiz (December 17, 2001; Duration: 3 hours)
1.225 Transportation Flow Systems Quiz (December 17, 2001; Duration: 3 hours) Student Name: Alias: Instructions: 1. This exam is open-book 2. No cooperation is permitted 3. Please write down your name
More informationAlgorithms. Outline! Approximation Algorithms. The class APX. The intelligence behind the hardware. ! Based on
6117CIT - Adv Topics in Computing Sci at Nathan 1 Algorithms The intelligence behind the hardware Outline! Approximation Algorithms The class APX! Some complexity classes, like PTAS and FPTAS! Illustration
More informationOPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS. w 1. v 2. v 3. Ruth J. Williams University of California, San Diego
OPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS v 2 w3 v1 w 2 w 1 v 3 Ruth J. Williams University of California, San Diego 1 PERSPECTIVE MQN SPN Sufficient conditions
More informationComputer Systems Modelling
Computer Systems Modelling Computer Laboratory Computer Science Tripos, Part II Michaelmas Term 2003 R. J. Gibbens Problem sheet William Gates Building JJ Thomson Avenue Cambridge CB3 0FD http://www.cl.cam.ac.uk/
More informationLoad Balancing in Distributed Service System: A Survey
Load Balancing in Distributed Service System: A Survey Xingyu Zhou The Ohio State University zhou.2055@osu.edu November 21, 2016 Xingyu Zhou (OSU) Load Balancing November 21, 2016 1 / 29 Introduction and
More informationDynamic Scheduling of Multiclass Queueing Networks. Caiwei Li
Dynamic Scheduling of Multiclass Queueing Networks A Thesis Presented to The Academic Faculty by Caiwei Li In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Industrial
More informationControl of Fork-Join Networks in Heavy-Traffic
in Heavy-Traffic Asaf Zviran Based on MSc work under the guidance of Rami Atar (Technion) and Avishai Mandelbaum (Technion) Industrial Engineering and Management Technion June 2010 Introduction Network
More informationSandwich shop : a queuing net work with finite disposable resources queue and infinite resources queue
Sandwich shop : a queuing net work with finite disposable resources queue and infinite resources queue Final project for ISYE 680: Queuing systems and Applications Hongtan Sun May 5, 05 Introduction As
More informationStochastic Models of Manufacturing Systems
Stochastic Models of Manufacturing Systems Ivo Adan Systems 2/49 Continuous systems State changes continuously in time (e.g., in chemical applications) Discrete systems State is observed at fixed regular
More informationContents Preface The Exponential Distribution and the Poisson Process Introduction to Renewal Theory
Contents Preface... v 1 The Exponential Distribution and the Poisson Process... 1 1.1 Introduction... 1 1.2 The Density, the Distribution, the Tail, and the Hazard Functions... 2 1.2.1 The Hazard Function
More informationSingle Machine Problems Polynomial Cases
DM204, 2011 SCHEDULING, TIMETABLING AND ROUTING Lecture 2 Single Machine Problems Polynomial Cases Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline
More informationEE 368. Weeks 3 (Notes)
EE 368 Weeks 3 (Notes) 1 State of a Queuing System State: Set of parameters that describe the condition of the system at a point in time. Why do we need it? Average size of Queue Average waiting time How
More informationQueueing Theory I Summary! Little s Law! Queueing System Notation! Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K "
Queueing Theory I Summary Little s Law Queueing System Notation Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K " Little s Law a(t): the process that counts the number of arrivals
More informationM/G/1 and M/G/1/K systems
M/G/1 and M/G/1/K systems Dmitri A. Moltchanov dmitri.moltchanov@tut.fi http://www.cs.tut.fi/kurssit/elt-53606/ OUTLINE: Description of M/G/1 system; Methods of analysis; Residual life approach; Imbedded
More informationEfficient Nonlinear Optimizations of Queuing Systems
Efficient Nonlinear Optimizations of Queuing Systems Mung Chiang, Arak Sutivong, and Stephen Boyd Electrical Engineering Department, Stanford University, CA 9435 Abstract We present a systematic treatment
More informationState-dependent and Energy-aware Control of Server Farm
State-dependent and Energy-aware Control of Server Farm Esa Hyytiä, Rhonda Righter and Samuli Aalto Aalto University, Finland UC Berkeley, USA First European Conference on Queueing Theory ECQT 2014 First
More informationStochastic Networks and Parameter Uncertainty
Stochastic Networks and Parameter Uncertainty Assaf Zeevi Graduate School of Business Columbia University Stochastic Processing Networks Conference, August 2009 based on joint work with Mike Harrison Achal
More informationRCPSP Single Machine Problems
DM204 Spring 2011 Scheduling, Timetabling and Routing Lecture 3 Single Machine Problems Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. Resource
More informationStability and Heavy Traffic Limits for Queueing Networks
Maury Bramson University of Minnesota Stability and Heavy Traffic Limits for Queueing Networks May 15, 2006 Springer Berlin Heidelberg New York Hong Kong London Milan Paris Tokyo Contents 1 Introduction...............................................
More informationIOE 202: lectures 11 and 12 outline
IOE 202: lectures 11 and 12 outline Announcements Last time... Queueing models intro Performance characteristics of a queueing system Steady state analysis of an M/M/1 queueing system Other queueing systems,
More informationStability, Capacity, and Scheduling of Multiclass Queueing Networks
Stability, Capacity, and Scheduling of Multiclass Queueing Networks A THESIS Presented to The Academic Faculty by John Jay Hasenbein In Partial Fulfillment of the Requirements for the Degree of Doctor
More informationStability of queueing networks. Maury Bramson. University of Minnesota
Probability Surveys Vol. 5 (2008) 169 345 ISSN: 1549-5787 DOI: 10.1214/08-PS137 Stability of queueing networks Received June 2008. Maury Bramson University of Minnesota e-mail: bramson@math.umn.edu Contents
More informationPBW 654 Applied Statistics - I Urban Operations Research
PBW 654 Applied Statistics - I Urban Operations Research Lecture 2.I Queuing Systems An Introduction Operations Research Models Deterministic Models Linear Programming Integer Programming Network Optimization
More informationFirefly algorithm in optimization of queueing systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 60, No. 2, 2012 DOI: 10.2478/v10175-012-0049-y VARIA Firefly algorithm in optimization of queueing systems J. KWIECIEŃ and B. FILIPOWICZ
More informationDesigning load balancing and admission control policies: lessons from NDS regime
Designing load balancing and admission control policies: lessons from NDS regime VARUN GUPTA University of Chicago Based on works with : Neil Walton, Jiheng Zhang ρ K θ is a useful regime to study the
More informationLecture 20: Reversible Processes and Queues
Lecture 20: Reversible Processes and Queues 1 Examples of reversible processes 11 Birth-death processes We define two non-negative sequences birth and death rates denoted by {λ n : n N 0 } and {µ n : n
More informationQueueing Theory and Simulation. Introduction
Queueing Theory and Simulation Based on the slides of Dr. Dharma P. Agrawal, University of Cincinnati and Dr. Hiroyuki Ohsaki Graduate School of Information Science & Technology, Osaka University, Japan
More informationMaximum pressure policies for stochastic processing networks
Maximum pressure policies for stochastic processing networks Jim Dai Joint work with Wuqin Lin at Northwestern Univ. The 2011 Lunteren Conference Jim Dai (Georgia Tech) MPPs January 18, 2011 1 / 55 Outline
More informationChapter 11. Output Analysis for a Single Model Prof. Dr. Mesut Güneş Ch. 11 Output Analysis for a Single Model
Chapter Output Analysis for a Single Model. Contents Types of Simulation Stochastic Nature of Output Data Measures of Performance Output Analysis for Terminating Simulations Output Analysis for Steady-state
More informationQueueing Theory (Part 4)
Queueing Theory (Part 4) Nonexponential Queueing Systems and Economic Analysis Queueing Theory-1 Queueing Models with Nonexponential Distributions M/G/1 Model Poisson input process, general service time
More informationFCFS Infinite Bipartite Matching with Applications
FF Infinite Bipartite Matching with Applications Workshop on ongestion Games, IM, NU ingapore December 2015 Gideon Weiss University of Haifa urvey of work since 2007, with Ed Kaplan, Rene aldentey, Ivo
More information6 Solving Queueing Models
6 Solving Queueing Models 6.1 Introduction In this note we look at the solution of systems of queues, starting with simple isolated queues. The benefits of using predefined, easily classified queues will
More informationSOLUTIONS IEOR 3106: Second Midterm Exam, Chapters 5-6, November 8, 2012
SOLUTIONS IEOR 3106: Second Midterm Exam, Chapters 5-6, November 8, 2012 This exam is closed book. YOU NEED TO SHOW YOUR WORK. Honor Code: Students are expected to behave honorably, following the accepted
More informationAn M/M/1 Retrial Queue with Unreliable Server 1
An M/M/1 Retrial Queue with Unreliable Server 1 Nathan P. Sherman 2 and Jeffrey P. Kharoufeh 3 Department of Operational Sciences Air Force Institute of Technology Abstract We analyze an unreliable M/M/1
More informationSession-Based Queueing Systems
Session-Based Queueing Systems Modelling, Simulation, and Approximation Jeroen Horters Supervisor VU: Sandjai Bhulai Executive Summary Companies often offer services that require multiple steps on the
More informationA cµ Rule for Two-Tiered Parallel Servers
1 A cµ Rule for Two-Tiered Parallel Servers Soroush Saghafian, Michael H. Veatch. Abstract The cµ rule is known to be optimal in many queueing systems with memoryless service and inter-arrival times. We
More informationAsymptotic Coupling of an SPDE, with Applications to Many-Server Queues
Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues Mohammadreza Aghajani joint work with Kavita Ramanan Brown University March 2014 Mohammadreza Aghajanijoint work Asymptotic with
More informationQueueing Theory II. Summary. ! M/M/1 Output process. ! Networks of Queue! Method of Stages. ! General Distributions
Queueing Theory II Summary! M/M/1 Output process! Networks of Queue! Method of Stages " Erlang Distribution " Hyperexponential Distribution! General Distributions " Embedded Markov Chains M/M/1 Output
More information