Transient behaviour in highly dependable Markovian systems: new regimes, many paths.

Size: px
Start display at page:

Download "Transient behaviour in highly dependable Markovian systems: new regimes, many paths."

Transcription

1 Transient behaviour in highly dependable Markovian systems: new regimes, many paths. Daniël Reijsbergen Pieter-Tjerk de Boer Werner Scheinhardt University of Twente RESIM, June 22nd, 2010

2 Outline Problem Setting Path-based Importance Sampling Asymptotic Regimes Many Paths for Fast Repairs What to Remember

3 Problem Setting Outline Problem Setting Path-based Importance Sampling Asymptotic Regimes Many Paths for Fast Repairs What to Remember

4 Problem Setting Model Checking (1) Automated verification of probabilistic systems. System model is specified in high-level language (e.g. Petri net). Model checker then verifies performance properties. E.g. the probability of overflow before time τ is smaller than p.

5 Problem Setting Model Checking (2) Model is translated to continuous time Markov chain with: transition rates q ss, starting state s 0, time bound τ, set of overflow states ( >n components have failed). Our interest: P(s 0, τ) = P(overflow before τ start in s 0 )

6 Problem Setting s 0

7 Problem Setting Model Checking (3) Probability P(s 0, τ) is estimated using simulation. Problem: in reliable systems, P(s 0, τ) is small. Solution: we use importance sampling.

8 Path-based Importance Sampling Outline Problem Setting Path-based Importance Sampling Asymptotic Regimes Many Paths for Fast Repairs What to Remember

9 Path-based Importance Sampling Importance Sampling (1) New measure tries to imitate the zero variance measure. If we are in state s with t time units remaining, then the ideal density of jumping to state s after δ time units is based on ratio P(s, t δ). P(s, t) P(s, t) is not known explicitly, so we use an approximation ˆP(s, t) and normalise.

10 Path-based Importance Sampling Importance Sampling (2) We use an approximation based on probability of the most likely path (or paths) to overflow from a state s. The result is a state-dependent and possibly time-dependent simulation density.

11 Path-based Importance Sampling Most Likely Paths (1) s 0 System with two parallel queues.

12 Path-based Importance Sampling Most Likely Paths (2) s 0 Most likely paths are the straight paths to overflow.

13 Path-based Importance Sampling Most Likely Paths (3) s 0 Most likely paths are the straight paths to overflow.

14 Path-based Importance Sampling Less Likely Paths (1) s 0 Path has an unnecessary arrival in queue 2.

15 Path-based Importance Sampling Less Likely Paths (2) s 0 Path has an unnecessary cycle from the empty state.

16 Path-based Importance Sampling Advantages Advantages of path-based importance sampling: Easy to generalise. Paths are found quickly using high-level model description. Performs well when probability contribution of other paths vanishes quickly.

17 Asymptotic Regimes Outline Problem Setting Path-based Importance Sampling Asymptotic Regimes Many Paths for Fast Repairs What to Remember

18 Asymptotic Regimes Causes of Rarity Probability contribution of non-straight paths may vanish because: 1. the extra jumps occur with vanishing probability 2. the extra sojourns consume an increasingly unlikely amount of time.

19 Asymptotic Regimes Asymptotic Regimes (1) τ 0 0 n τ 0: short mission time. 0: highly dependable system. n : high redundancy. : highly recoverable system.

20 Asymptotic Regimes τ 0 τ 0 0 n Extra sojourns consume too much time. Possible (partial) solution: forcing. Nicola, Nakayama, Heidelberger and Goyal. Fast simulation of highly dependable systems with general failure and repair processes. IEEE Trans on Comp, Nakayama and Shahabuddin. Quick simulation methods for estimating the unreliability of regenerative models of large, highly reliable systems. PEIS, 2004.

21 Asymptotic Regimes τ 0 s 0 Extra sojourns take too much time.

22 Asymptotic Regimes 0 τ 0 0 n Sojourns in empty state consume too much time. Failures before repairs occur with vanishing probability. Possible solution: forcing + (balanced) failure baising. Nicola, Nakayama, Heidelberger and Goyal. Fast simulation of highly dependable systems with general failure and repair processes. IEEE Trans on Comp, de Boer, L Ecuyer, Rubino and Tuffin. Estimating the probability of a rare event over a finite time horizon. Winter Simulation Conference, 2007.

23 Asymptotic Regimes 0 s 0 Extra jumps increasingly unlikely.

24 Asymptotic Regimes n τ 0 0 n Extra sojourns consume too much time. Often considered in combination with and. Possible solutions: based on large deviations.

25 Asymptotic Regimes τ 0 0 n Failures before repairs occur with vanishing probability. So far, only bounds (based on estimates) for large τ. Nakayama and Shahabuddin. Quick simulation methods for estimating the unreliability of regenerative models of large, highly reliable systems. PEIS, 2004.

26 Asymptotic Regimes s 0 If : contribution does not vanish asymptotically!

27 Many Paths for Fast Repairs Outline Problem Setting Path-based Importance Sampling Asymptotic Regimes Many Paths for Fast Repairs What to Remember

28 Many Paths for Fast Repairs Many Paths for Fast Repairs (1) Probability of less likely paths does not vanish when. Solution: use the regenerative behaviour of the model. Probability of the straight paths remains a good guess for π = P( overflow during busy cycle).

29 Many Paths for Fast Repairs Many Paths for Fast Repairs (2) When, duration of busy cycle = sojourn time in s 0. This is exponentially distributed, with rate equal to exit rate of empty state. # busy cycles before overflow geometrically distributed(π). Distribution of time until overflow: exponential.

30 Many Paths for Fast Repairs Many Paths for Fast Repairs (3) ˆP(s 0, t) = 1 e tπ tπ when π small For all states s s 0, we set ˆP(s, t) to be the sum of 1. P(straight paths from s) 2. the probability of the path that falls back to the empty state, multiplied by ˆP(s 0, t) Resulting measure, combined with forcing, leads to measure that works well in regimes 0, τ 0 and.

31 Conclusions Outline Problem Setting Path-based Importance Sampling Asymptotic Regimes Many Paths for Fast Repairs What to Remember

32 Conclusions What to Remember Path-based importance sampling can efficiently be applied in computerised model checking. Challenges due to different forms of rarity. Future work: Quick search for most likely paths. Exploit regenerative behaviour in more general scenarios.

Hypothesis Testing for Rare-Event Simulation: Limitations and Possibilities

Hypothesis Testing for Rare-Event Simulation: Limitations and Possibilities Hypothesis Testing for Rare-Event Simulation: Limitations and Possibilities Daniël Reijsbergen 1, Pieter-Tjerk de Boer 2(B), and Werner Scheinhardt 2 1 University of Edinburgh, Edinburgh, Scotland, UK

More information

On the estimation of the mean time to failure by simulation

On the estimation of the mean time to failure by simulation On the estimation of the mean time to failure by simulation Peter Glynn, Marvin Nakayama, Bruno Tuffin To cite this version: Peter Glynn, Marvin Nakayama, Bruno Tuffin. On the estimation of the mean time

More information

On Derivative Estimation of the Mean Time to Failure in Simulations of Highly Reliable Markovian Systems

On Derivative Estimation of the Mean Time to Failure in Simulations of Highly Reliable Markovian Systems On Derivative Estimation of the Mean Time to Failure in Simulations of Highly Reliable Markovian Systems Marvin K. Nakayama Department of Computer and Information Science New Jersey Institute of Technology

More information

IMPORTANCE SAMPLING SIMULATIONS OF PHASE-TYPE QUEUES

IMPORTANCE SAMPLING SIMULATIONS OF PHASE-TYPE QUEUES Proceedings of the 29 Winter Simulation Conference M. D. Rossetti, R. R. Hill, B. Johansson, A. Dunkin, and R. G. Ingalls, eds. IMPORTANCE SAMPLING SIMULATIONS OF PHASE-TYPE QUEUES Poul E. Heegaard Department

More information

Efficient Simulation of a Tandem Queue with Server Slow-down

Efficient Simulation of a Tandem Queue with Server Slow-down Efficient Simulation of a Tandem Queue with Server Slow-down D.I. Miretskiy, W.R.W. Scheinhardt, M.R.H. Mandjes Abstract Tandem Jackson networks, and more sophisticated variants, have found widespread

More information

THE BALANCED LIKELIHOOD RATIO METHOD FOR ESTIMATING PERFORMANCE MEASURES OF HIGHLY RELIABLE SYSTEMS

THE BALANCED LIKELIHOOD RATIO METHOD FOR ESTIMATING PERFORMANCE MEASURES OF HIGHLY RELIABLE SYSTEMS Proceedings of the 1998 Winter Simulation Conference D.J. Medeiros, E.F. Watson, J.S. Carson and M.S. Manivannan, eds. THE BALANCED LIKELIHOOD RATIO METHOD FOR ESTIMATING PERFORMANCE MEASURES OF HIGHLY

More information

14 Random Variables and Simulation

14 Random Variables and Simulation 14 Random Variables and Simulation In this lecture note we consider the relationship between random variables and simulation models. Random variables play two important roles in simulation models. We assume

More information

IE Advanced Simulation Experiment Design and Analysis. Hong Wan Purdue University

IE Advanced Simulation Experiment Design and Analysis. Hong Wan Purdue University IE Advanced Simulation Experiment Design and Analysis Hong Wan Purdue University Spring, 2007, based on Prof. Barry L. Nelson s notes for IEMS465, Northwestern University 1 SIMULATION OVERVIEW One view

More information

Juan A. Carrasco a a Departament d'enginyeria Electrònica, Universitat Politècnica de. Catalunya, Barcelona, Spain Published online: 15 Feb 2007.

Juan A. Carrasco a a Departament d'enginyeria Electrònica, Universitat Politècnica de. Catalunya, Barcelona, Spain Published online: 15 Feb 2007. This article was downloaded by: [Mary Ann Muller] On: 21 January 2014, At: 07:06 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer

More information

Rare Events, Splitting, and Quasi-Monte Carlo

Rare Events, Splitting, and Quasi-Monte Carlo Rare Events, Splitting, and Quasi-Monte Carlo PIERRE L ECUYER and VALÉRIE DEMERS Université de Montréal, Canada and BRUNO TUFFIN IRISA INRIA, Rennes, France In the context of rare-event simulation, splitting

More information

MAT SYS 5120 (Winter 2012) Assignment 5 (not to be submitted) There are 4 questions.

MAT SYS 5120 (Winter 2012) Assignment 5 (not to be submitted) There are 4 questions. MAT 4371 - SYS 5120 (Winter 2012) Assignment 5 (not to be submitted) There are 4 questions. Question 1: Consider the following generator for a continuous time Markov chain. 4 1 3 Q = 2 5 3 5 2 7 (a) Give

More information

Contents LIST OF TABLES... LIST OF FIGURES... xvii. LIST OF LISTINGS... xxi PREFACE. ...xxiii

Contents LIST OF TABLES... LIST OF FIGURES... xvii. LIST OF LISTINGS... xxi PREFACE. ...xxiii LIST OF TABLES... xv LIST OF FIGURES... xvii LIST OF LISTINGS... xxi PREFACE...xxiii CHAPTER 1. PERFORMANCE EVALUATION... 1 1.1. Performance evaluation... 1 1.2. Performance versus resources provisioning...

More information

Bounded Normal Approximation in Highly Reliable Markovian Systems

Bounded Normal Approximation in Highly Reliable Markovian Systems Bounded Normal Approximation in Highly Reliable Markovian Systems Bruno Tuffin To cite this version: Bruno Tuffin. Bounded Normal Approximation in Highly Reliable Markovian Systems. [Research Report] RR-3020,

More information

Markov Model. Model representing the different resident states of a system, and the transitions between the different states

Markov Model. Model representing the different resident states of a system, and the transitions between the different states Markov Model Model representing the different resident states of a system, and the transitions between the different states (applicable to repairable, as well as non-repairable systems) System behavior

More information

This paper discusses the application of the likelihood ratio gradient estimator to simulations of highly dependable systems. This paper makes the foll

This paper discusses the application of the likelihood ratio gradient estimator to simulations of highly dependable systems. This paper makes the foll Likelihood Ratio Sensitivity Analysis for Markovian Models of Highly Dependable Systems Marvin K. Nakayama IBM Thomas J. Watson Research Center P.O. Box 704 Yorktown Heights, New York 10598 Ambuj Goyal

More information

ASYMPTOTIC ROBUSTNESS OF ESTIMATORS IN RARE-EVENT SIMULATION. Bruno Tuffin. IRISA/INRIA, Campus de Beaulieu Rennes Cedex, France

ASYMPTOTIC ROBUSTNESS OF ESTIMATORS IN RARE-EVENT SIMULATION. Bruno Tuffin. IRISA/INRIA, Campus de Beaulieu Rennes Cedex, France Submitted to: Proceedings of the 2007 INFORMS Simulation Society Research Workshop ASYMPTOTIC ROBUSTNESS OF ESTIMATORS IN RARE-EVENT SIMULATION Jose H. Blanchet Department of Statistics Harvard University,

More information

AN EFFICIENT IMPORTANCE SAMPLING METHOD FOR RARE EVENT SIMULATION IN LARGE SCALE TANDEM NETWORKS

AN EFFICIENT IMPORTANCE SAMPLING METHOD FOR RARE EVENT SIMULATION IN LARGE SCALE TANDEM NETWORKS Proceedings of the Winter Simulation Conference E. Yücesan, C.-H. Chen, J. L. Snowdon, and J. M. Charnes, eds. AN EFFICIENT IMPORTANCE SAMPLING METHOD FOR RARE EVENT SIMULATION IN LARGE SCALE TANDEM NETWORKS

More information

57:022 Principles of Design II Final Exam Solutions - Spring 1997

57: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 information

Stochastic models in product form: the (E)RCAT methodology

Stochastic models in product form: the (E)RCAT methodology Stochastic models in product form: the (E)RCAT methodology 1 Maria Grazia Vigliotti 2 1 Dipartimento di Informatica Università Ca Foscari di Venezia 2 Department of Computing Imperial College London Second

More information

Sensitivity of System Reliability to Usage Profile Changes

Sensitivity of System Reliability to Usage Profile Changes Sensitivity of System Reliability to Usage Profile Changes Kim Weyns Department of Communication Systems, Lund University PO Box 118 SE-211 00 LUND, Sweden kimweyns@telecomlthse Per Runeson Department

More information

Availability. M(t) = 1 - e -mt

Availability. M(t) = 1 - e -mt Availability Availability - A(t) the probability that the system is operating correctly and is available to perform its functions at the instant of time t More general concept than reliability: failure

More information

IMPORTANCE SAMPLING USING THE SEMI-REGENERATIVE METHOD. Marvin K. Nakayama

IMPORTANCE SAMPLING USING THE SEMI-REGENERATIVE METHOD. Marvin K. Nakayama Proceedings of the 2001 Winter Simulation Conference B.A.Peters,J.S.Smith,D.J.Medeiros,andM.W.Rohrer,eds. IMPORTANCE SAMPLING USING THE SEMI-REGENERATIVE METHOD James M. Calvin Department of Computer Science

More information

Stochastic Modelling of a Computer System with Software Redundancy and Priority to Hardware Repair

Stochastic Modelling of a Computer System with Software Redundancy and Priority to Hardware Repair Stochastic Modelling of a Computer System with Software Redundancy and Priority to Hardware Repair Abstract V.J. Munday* Department of Statistics, Ramjas College, University of Delhi, Delhi 110007 (India)

More information

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1 Irreducibility Irreducible every state can be reached from every other state For any i,j, exist an m 0, such that i,j are communicate, if the above condition is valid Irreducible: all states are communicate

More information

Markov Reliability and Availability Analysis. Markov Processes

Markov Reliability and Availability Analysis. Markov Processes Markov Reliability and Availability Analysis Firma convenzione Politecnico Part II: Continuous di Milano e Time Veneranda Discrete Fabbrica State del Duomo di Milano Markov Processes Aula Magna Rettorato

More information

MODELLING DYNAMIC RELIABILITY VIA FLUID PETRI NETS

MODELLING DYNAMIC RELIABILITY VIA FLUID PETRI NETS MODELLING DYNAMIC RELIABILITY VIA FLUID PETRI NETS Daniele Codetta-Raiteri, Dipartimento di Informatica, Università di Torino, Italy Andrea Bobbio, Dipartimento di Informatica, Università del Piemonte

More information

THE STATIC network reliability problem, which consists

THE STATIC network reliability problem, which consists 1 Approximate Zero-Variance Importance Sampling for Static Network Reliability Estimation Pierre L Ecuyer, Gerardo Rubino, Samira Saggadi and Bruno Tuffin Abstract We propose a new Monte Carlo method,

More information

An Introduction to Stochastic Modeling

An Introduction to Stochastic Modeling F An Introduction to Stochastic Modeling Fourth Edition Mark A. Pinsky Department of Mathematics Northwestern University Evanston, Illinois Samuel Karlin Department of Mathematics Stanford University Stanford,

More information

Description Framework for Proxel-Based Simulation of a General Class of Stochastic Models

Description Framework for Proxel-Based Simulation of a General Class of Stochastic Models Description Framework for Proxel-Based Simulation of a General Class of Stochastic Models Sanja Lazarova-Molnar, Graham Horton Institute for Simulation and Graphics University of Magdeburg Universitaetsplatz

More information

Representation and Analysis of Transfer Lines. with Machines That Have Different Processing Rates

Representation and Analysis of Transfer Lines. with Machines That Have Different Processing Rates March, 1985 Revised November, 1985 LIDS-P-1446 Representation and Analysis of Transfer Lines with Machines That Have Different Processing Rates by Stanley B. Gershwin 35-433 Massachusetts Institute of

More information

An M/M/1 Queue in Random Environment with Disasters

An M/M/1 Queue in Random Environment with Disasters An M/M/1 Queue in Random Environment with Disasters Noam Paz 1 and Uri Yechiali 1,2 1 Department of Statistics and Operations Research School of Mathematical Sciences Tel Aviv University, Tel Aviv 69978,

More information

MIT Manufacturing Systems Analysis Lectures 6 9: Flow Lines

MIT Manufacturing Systems Analysis Lectures 6 9: Flow Lines 2.852 Manufacturing Systems Analysis 1/165 Copyright 2010 c Stanley B. Gershwin. MIT 2.852 Manufacturing Systems Analysis Lectures 6 9: Flow Lines Models That Can Be Analyzed Exactly Stanley B. Gershwin

More information

Quantitative evaluation of Dependability

Quantitative evaluation of Dependability Quantitative evaluation of Dependability 1 Quantitative evaluation of Dependability Faults are the cause of errors and failures. Does the arrival time of faults fit a probability distribution? If so, what

More information

APPROXIMATE ZERO-VARIANCE SIMULATION

APPROXIMATE ZERO-VARIANCE SIMULATION Proceedings of the 2008 Winter Simulation Conference S. J. Mason, R. R. Hill, L. Moench, O. Rose, eds. APPROXIMATE ZERO-VARIANCE SIMULATION Pierre L Ecuyer DIRO, Université de Montreal C.P. 6128, Succ.

More information

QUICK SIMULATION METHODS FOR ESTIMATING THE UNRELIABILITY OF REGENERATIVE MODELS OF LARGE, HIGHLY RELIABLE SYSTEMS

QUICK SIMULATION METHODS FOR ESTIMATING THE UNRELIABILITY OF REGENERATIVE MODELS OF LARGE, HIGHLY RELIABLE SYSTEMS Probability in the Engineering and Informational Sciences, 18, 2004, 339 368+ Printed in the U+S+A+ QUICK SIMULATION METHODS FOR ESTIMATING THE UNRELIABILITY OF REGENERATIVE MODELS OF LARGE, HIGHLY RELIABLE

More information

EE 445 / 850: Final Examination

EE 445 / 850: Final Examination EE 445 / 850: Final Examination Date and Time: 3 Dec 0, PM Room: HLTH B6 Exam Duration: 3 hours One formula sheet permitted. - Covers chapters - 5 problems each carrying 0 marks - Must show all calculations

More information

CPSC 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 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 information

STEADY-STATE BEHAVIOR OF AN M/M/1 QUEUE IN RANDOM ENVIRONMENT SUBJECT TO SYSTEM FAILURES AND REPAIRS. S. Sophia 1, B. Praba 2

STEADY-STATE BEHAVIOR OF AN M/M/1 QUEUE IN RANDOM ENVIRONMENT SUBJECT TO SYSTEM FAILURES AND REPAIRS. S. Sophia 1, B. Praba 2 International Journal of Pure and Applied Mathematics Volume 101 No. 2 2015, 267-279 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v101i2.11

More information

Chapter 3 Balance equations, birth-death processes, continuous Markov Chains

Chapter 3 Balance equations, birth-death processes, continuous Markov Chains Chapter 3 Balance equations, birth-death processes, continuous Markov Chains Ioannis Glaropoulos November 4, 2012 1 Exercise 3.2 Consider a birth-death process with 3 states, where the transition rate

More information

Dependable Computer Systems

Dependable Computer Systems Dependable Computer Systems Part 3: Fault-Tolerance and Modelling Contents Reliability: Basic Mathematical Model Example Failure Rate Functions Probabilistic Structural-Based Modeling: Part 1 Maintenance

More information

Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk

Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk ANSAPW University of Queensland 8-11 July, 2013 1 Outline (I) Fluid

More information

Link Models for Circuit Switching

Link Models for Circuit Switching Link Models for Circuit Switching The basis of traffic engineering for telecommunication networks is the Erlang loss function. It basically allows us to determine the amount of telephone traffic that can

More information

Single-part-type, multiple stage systems. Lecturer: Stanley B. Gershwin

Single-part-type, multiple stage systems. Lecturer: Stanley B. Gershwin Single-part-type, multiple stage systems Lecturer: Stanley B. Gershwin Flow Line... also known as a Production or Transfer Line. M 1 B 1 M 2 B 2 M 3 B 3 M 4 B 4 M 5 B 5 M 6 Machine Buffer Machines are

More information

6 Solving Queueing Models

6 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 information

Quantitative evaluation of Dependability

Quantitative evaluation of Dependability Quantitative evaluation of Dependability 1 Quantitative evaluation of Dependability Faults are the cause of errors and failures. Does the arrival time of faults fit a probability distribution? If so, what

More information

Robustesse des techniques de Monte Carlo dans l analyse d événements rares

Robustesse des techniques de Monte Carlo dans l analyse d événements rares Institut National de Recherche en Informatique et Automatique Institut de Recherche en Informatique et Systèmes Aléatoires Robustesse des techniques de Monte Carlo dans l analyse d événements rares H.

More information

Introduction to rare event simulation

Introduction to rare event simulation 1/35 rare event simulation (some works with colleagues H. Cancela, V. Demers, P. L Ecuyer, G. Rubino) INRIA - Centre Bretagne Atlantique, Rennes Aussois, June 2008, AEP9 Outline 2/35 1 2 3 4 5 Introduction:

More information

Introduction to Queuing Networks Solutions to Problem Sheet 3

Introduction to Queuing Networks Solutions to Problem Sheet 3 Introduction to Queuing Networks Solutions to Problem Sheet 3 1. (a) The state space is the whole numbers {, 1, 2,...}. The transition rates are q i,i+1 λ for all i and q i, for all i 1 since, when a bus

More information

CSE 4201, Ch. 6. Storage Systems. Hennessy and Patterson

CSE 4201, Ch. 6. Storage Systems. Hennessy and Patterson CSE 4201, Ch. 6 Storage Systems Hennessy and Patterson Challenge to the Disk The graveyard is full of suitors Ever heard of Bubble Memory? There are some technologies that refuse to die (silicon, copper...).

More information

N.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a

N.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure

More information

Stochastic Petri Nets. Jonatan Lindén. Modelling SPN GSPN. Performance measures. Almost none of the theory. December 8, 2010

Stochastic Petri Nets. Jonatan Lindén. Modelling SPN GSPN. Performance measures. Almost none of the theory. December 8, 2010 Stochastic Almost none of the theory December 8, 2010 Outline 1 2 Introduction A Petri net (PN) is something like a generalized automata. A Stochastic Petri Net () a stochastic extension to Petri nets,

More information

Supervisory Control: Advanced Theory and Applications

Supervisory Control: Advanced Theory and Applications Supervisory Control: Advanced Theory and Applications Dr Rong Su S1-B1b-59, School of EEE Nanyang Technological University Tel: +65 6790-6042, Email: rsu@ntu.edu.sg EE6226, Discrete Event Systems 1 Introduction

More information

Queues and Queueing Networks

Queues 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 information

Proxel-Based Simulation of Stochastic Petri Nets Containing Immediate Transitions

Proxel-Based Simulation of Stochastic Petri Nets Containing Immediate Transitions Electronic Notes in Theoretical Computer Science Vol. 85 No. 4 (2003) URL: http://www.elsevier.nl/locate/entsc/volume85.html Proxel-Based Simulation of Stochastic Petri Nets Containing Immediate Transitions

More information

Adam Caromicoli. Alan S. Willsky 1. Stanley B. Gershwin 2. Abstract

Adam Caromicoli. Alan S. Willsky 1. Stanley B. Gershwin 2. Abstract December 1987 LIDS-P-1727 MULTIPLE TIME SCALE ANALYSIS OF MANUFACTURING SYSTEMS Adam Caromicoli Alan S. Willsky 1 Stanley B. Gershwin 2 Abstract In this paper we use results on the aggregation of singularly

More information

Analysis of the Durability of Replicated Distributed Storage Systems

Analysis of the Durability of Replicated Distributed Storage Systems Analysis of the Durability of Replicated Distributed Storage Systems Sriram Ramabhadran Joseph Pasquale ABSRAC In this paper, we investigate the roles of replication vs. repair to achieve durability in

More information

Outline. Finite source queue M/M/c//K Queues with impatience (balking, reneging, jockeying, retrial) Transient behavior Advanced Queue.

Outline. Finite source queue M/M/c//K Queues with impatience (balking, reneging, jockeying, retrial) Transient behavior Advanced Queue. Outline Finite source queue M/M/c//K Queues with impatience (balking, reneging, jockeying, retrial) Transient behavior Advanced Queue Batch queue Bulk input queue M [X] /M/1 Bulk service queue M/M [Y]

More information

of an algorithm for automated cause-consequence diagram construction.

of an algorithm for automated cause-consequence diagram construction. Loughborough University Institutional Repository Development of an algorithm for automated cause-consequence diagram construction. This item was submitted to Loughborough University's Institutional Repository

More information

Session-Based Queueing Systems

Session-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 information

Reliable Computing I

Reliable Computing I Instructor: Mehdi Tahoori Reliable Computing I Lecture 5: Reliability Evaluation INSTITUTE OF COMPUTER ENGINEERING (ITEC) CHAIR FOR DEPENDABLE NANO COMPUTING (CDNC) National Research Center of the Helmholtz

More information

Glossary availability cellular manufacturing closed queueing network coefficient of variation (CV) conditional probability CONWIP

Glossary 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 information

Stochastic process. X, a series of random variables indexed by t

Stochastic process. X, a series of random variables indexed by t Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,

More information

Varieties of Stochastic Calculi

Varieties of Stochastic Calculi Research is what I'm doing when I don't know what I'm doing. Wernher Von Braun. Artificial Biochemistry Varieties of Stochastic Calculi Microsoft Research Trento, 26-5-22..26 www.luca.demon.co.uk/artificialbiochemistry.htm

More information

COMBINATION OF CONDITIONAL MONTE CARLO AND APPROXIMATE ZERO-VARIANCE IMPORTANCE SAMPLING FOR NETWORK RELIABILITY ESTIMATION

COMBINATION OF CONDITIONAL MONTE CARLO AND APPROXIMATE ZERO-VARIANCE IMPORTANCE SAMPLING FOR NETWORK RELIABILITY ESTIMATION Proceedings of the 2010 Winter Simulation Conference B. Johansson, S. Jain, J. Montoya-Torres, J. Hugan, and E. Yücesan, eds. COMBINATION OF CONDITIONAL MONTE CARLO AND APPROXIMATE ZERO-VARIANCE IMPORTANCE

More information

SPLITTING FOR RARE-EVENT SIMULATION. ˆγ n has relative error. RE[ ˆγ n ] = (Var[ ˆγ n]) 1/2

SPLITTING FOR RARE-EVENT SIMULATION. ˆγ n has relative error. RE[ ˆγ n ] = (Var[ ˆγ n]) 1/2 Proceedings of the 2006 Winter Simulation Conference L. F. Perrone, F. P. Wieland, J. Liu, B. G. Lawson, D. M. Nicol, and R. M. Fujimoto, eds. SPLITTING FOR RARE-EVENT SIMULATION Pierre L Ecuyer Valérie

More information

Fast Estimation of the Statistics of Excessive Backlogs in Tandem Networks of Queues.

Fast Estimation of the Statistics of Excessive Backlogs in Tandem Networks of Queues. Fast Estimation of the Statistics of Excessive Backlogs in Tandem Networks of Queues. M.R. FRATER and B.D.O. ANDERSON Department of Systems Engineering Australian National University The estimation of

More information

Stochastic Models 3 (1998)

Stochastic Models 3 (1998) CALCULATING TRANSIENT CHARACTERISTICS OF THE ERLANG LOSS MODEL BY NUMERICAL TRANSFORM INVERSION J. ABATE W. WHITT 900 Hammond Road AT&T Labs Research Ridgewood, NJ 07450-2908 Room 2C-178 Murray Hill, NJ

More information

Advanced Queueing Theory

Advanced Queueing Theory Advanced Queueing Theory 1 Networks of queues (reversibility, output theorem, tandem networks, partial balance, product-form distribution, blocking, insensitivity, BCMP networks, mean-value analysis, Norton's

More information

57:022 Principles of Design II Midterm Exam #2 Solutions

57:022 Principles of Design II Midterm Exam #2 Solutions 57:022 Principles of Design II Midterm Exam #2 Solutions Part: I II III IV V Total Possible Pts: 20 15 12 16 12 75 PART ONE Indicate "+" if True and "O" if False: _+_a. If a component's lifetime has exponential

More information

CASPA - A Tool for Symbolic Performance Evaluation and Stochastic Model Checking

CASPA - A Tool for Symbolic Performance Evaluation and Stochastic Model Checking CASPA - A Tool for Symbolic Performance Evaluation and Stochastic Model Checking Boudewijn R. Haverkort 1, Matthias Kuntz 1, Martin Riedl 2, Johann Schuster 2, Markus Siegle 2 1 : Universiteit Twente 2

More information

A sequential hypothesis test based on a generalized Azuma inequality 1

A sequential hypothesis test based on a generalized Azuma inequality 1 A sequential hypothesis test based on a generalized Azuma inequality 1 Daniël Reijsbergen a,2, Werner Scheinhardt b, Pieter-Tjerk de Boer b a Laboratory for Foundations of Computer Science, University

More information

J. MEDHI STOCHASTIC MODELS IN QUEUEING THEORY

J. MEDHI STOCHASTIC MODELS IN QUEUEING THEORY J. MEDHI STOCHASTIC MODELS IN QUEUEING THEORY SECOND EDITION ACADEMIC PRESS An imprint of Elsevier Science Amsterdam Boston London New York Oxford Paris San Diego San Francisco Singapore Sydney Tokyo Contents

More information

Performance Evaluation of Queuing Systems

Performance 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 information

NON-STATIONARY QUEUE SIMULATION ANALYSIS USING TIME SERIES

NON-STATIONARY QUEUE SIMULATION ANALYSIS USING TIME SERIES Proceedings of the 2003 Winter Simulation Conference S. Chick, P. J. Sánchez, D. Ferrin, and D. J. Morrice, eds. NON-STATIONARY QUEUE SIMULATION ANALYSIS USING TIME SERIES Rita Marques Brandão Departamento

More information

Reliability Analysis of Two-Unit Warm Standby System Subject to Hardware and Human Error Failures

Reliability Analysis of Two-Unit Warm Standby System Subject to Hardware and Human Error Failures Reliability Analysis of Two-Unit Warm Standby System Subject to Hardware and Human Error Failures Pravindra Singh 1, Pankaj Kumar 2 & Anil Kumar 3 1 C. C. S. University, Meerut,U.P 2 J.N.V.U Jodhpur &

More information

Proceedings of the 2017 Winter Simulation Conference W. K. V. Chan, A. D Ambrogio, G. Zacharewicz, N. Mustafee, G. Wainer, and E. Page, eds.

Proceedings of the 2017 Winter Simulation Conference W. K. V. Chan, A. D Ambrogio, G. Zacharewicz, N. Mustafee, G. Wainer, and E. Page, eds. Proceedings of the 207 Winter Simulation Conference W. K. V. Chan, A. D Ambrogio, G. Zacharewicz, N. Mustafee, G. Wainer, and E. Page, eds. ON THE ESTIMATION OF THE MEAN TIME TO FAILURE BY SIMULATION Peter

More information

A STUDY OF ASYMPTOTIC AVAILABILITY MODELING FOR A FAILURE AND A REPAIR RATES FOLLOWING A WEIBULL DISTRIBUTION

A STUDY OF ASYMPTOTIC AVAILABILITY MODELING FOR A FAILURE AND A REPAIR RATES FOLLOWING A WEIBULL DISTRIBUTION A STUDY OF ASYMPTOTIC AVAILABILITY MODELING FOR A FAILURE AND A REPAIR RATES FOLLOWING A WEIBULL DISTRIBUTION Salem Bahri a, Fethi Ghribi b, Habib Ben Bacha a,c a Electro Mechanical systems laboratory

More information

A FAST MATRIX-ANALYTIC APPROXIMATION FOR THE TWO CLASS GI/G/1 NON-PREEMPTIVE PRIORITY QUEUE

A FAST MATRIX-ANALYTIC APPROXIMATION FOR THE TWO CLASS GI/G/1 NON-PREEMPTIVE PRIORITY QUEUE A FAST MATRIX-ANAYTIC APPROXIMATION FOR TE TWO CASS GI/G/ NON-PREEMPTIVE PRIORITY QUEUE Gábor orváth Department of Telecommunication Budapest University of Technology and Economics. Budapest Pf. 9., ungary

More information

Poisson Processes. Stochastic Processes. Feb UC3M

Poisson Processes. Stochastic Processes. Feb UC3M Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written

More information

THIS paper introduces a simple practical correction to

THIS paper introduces a simple practical correction to VOLOVOI: NON-MARKOVIAN ASYMPTOTIC EFFECTS CORRECTION 1 Correcting for Non-Markovian Asymptotic Effects Using Markovian Representation Vitali Volovoi arxiv:175.17v2 [cs.pf] 16 May 217 Abstract Asymptotic

More information

An M/M/1 Retrial Queue with Unreliable Server 1

An 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 information

A Study on M x /G/1 Queuing System with Essential, Optional Service, Modified Vacation and Setup time

A Study on M x /G/1 Queuing System with Essential, Optional Service, Modified Vacation and Setup time A Study on M x /G/1 Queuing System with Essential, Optional Service, Modified Vacation and Setup time E. Ramesh Kumar 1, L. Poornima 2 1 Associate Professor, Department of Mathematics, CMS College of Science

More information

Markov Models for Reliability Modeling

Markov Models for Reliability Modeling Markov Models for Reliability Modeling Prof. Naga Kandasamy ECE Department, Drexel University, Philadelphia, PA 904 Many complex systems cannot be easily modeled in a combinatorial fashion. The corresponding

More information

A Hostile model for network reliability analysis

A Hostile model for network reliability analysis Croatian Operational Research Review 489 CRORR 8(2017), 489 498 A Hostile model for network reliability analysis Daniel Lena 1, Franco Robledo 1 and Pablo Romero 1, 1 Facultad de Ingenería, Universidad

More information

MARKOV PROCESSES. Valerio Di Valerio

MARKOV PROCESSES. Valerio Di Valerio MARKOV PROCESSES Valerio Di Valerio Stochastic Process Definition: a stochastic process is a collection of random variables {X(t)} indexed by time t T Each X(t) X is a random variable that satisfy some

More information

Modeling Common Cause Failures in Diverse Components with Fault Tree Applications

Modeling Common Cause Failures in Diverse Components with Fault Tree Applications Modeling Common Cause Failures in Diverse s with Fault Tree Applications Joseph R. Belland, Isograph Inc. Key Words: Common Cause Failures, Fault Trees, Beta Factor SUMMARY & CONCLUSIONS A common cause

More information

Fault Tolerance. Dealing with Faults

Fault Tolerance. Dealing with Faults Fault Tolerance Real-time computing systems must be fault-tolerant: they must be able to continue operating despite the failure of a limited subset of their hardware or software. They must also allow graceful

More information

Readings: Finish Section 5.2

Readings: Finish Section 5.2 LECTURE 19 Readings: Finish Section 5.2 Lecture outline Markov Processes I Checkout counter example. Markov process: definition. -step transition probabilities. Classification of states. Example: Checkout

More information

The square root rule for adaptive importance sampling

The square root rule for adaptive importance sampling The square root rule for adaptive importance sampling Art B. Owen Stanford University Yi Zhou January 2019 Abstract In adaptive importance sampling, and other contexts, we have unbiased and uncorrelated

More information

Uniform random numbers generators

Uniform random numbers generators Uniform random numbers generators Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2707/ OUTLINE: The need for random numbers; Basic steps in generation; Uniformly

More information

Queuing Theory. Using the Math. Management Science

Queuing 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 information

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University

More information

System Reliability Analysis. CS6323 Networks and Systems

System Reliability Analysis. CS6323 Networks and Systems System Reliability Analysis CS6323 Networks and Systems Topics Combinatorial Models for reliability Topology-based (structured) methods for Series Systems Parallel Systems Reliability analysis for arbitrary

More information

Proceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds.

Proceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds. Proceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds. ON THE ROBUSTNESS OF FISHMAN S BOUND-BASED METHOD FOR THE NETWORK

More information

Reliability Analysis of a Fuel Supply System in Automobile Engine

Reliability Analysis of a Fuel Supply System in Automobile Engine ISBN 978-93-84468-19-4 Proceedings of International Conference on Transportation and Civil Engineering (ICTCE'15) London, March 21-22, 2015, pp. 1-11 Reliability Analysis of a Fuel Supply System in Automobile

More information

Reliability of Technical Systems

Reliability of Technical Systems Reliability of Technical Systems Main Topics 1. Short Introduction, Reliability Parameters: Failure Rate, Failure Probability, etc. 2. Some Important Reliability Distributions 3. Component Reliability

More information

Answers to selected exercises

Answers to selected exercises Answers to selected exercises A First Course in Stochastic Models, Henk C. Tijms 1.1 ( ) 1.2 (a) Let waiting time if passengers already arrived,. Then,, (b) { (c) Long-run fraction for is (d) Let waiting

More information

The range of tree-indexed random walk

The range of tree-indexed random walk The range of tree-indexed random walk Jean-François Le Gall, Shen Lin Institut universitaire de France et Université Paris-Sud Orsay Erdös Centennial Conference July 2013 Jean-François Le Gall (Université

More information

Markov Chain Profit Modelling and Evaluation between Two Dissimilar Systems under Two Types of Failures

Markov Chain Profit Modelling and Evaluation between Two Dissimilar Systems under Two Types of Failures Available at http://pvamu.edu/aam Appl. Appl. Math. SSN: 9-966 Vol., ssue (December 06), pp. 89-50 Applications and Applied Mathematics: An nternational Journal (AAM) Markov Chain Profit Modelling and

More information

STOCHASTIC MODELLING OF A COMPUTER SYSTEM WITH HARDWARE REDUNDANCY SUBJECT TO MAXIMUM REPAIR TIME

STOCHASTIC MODELLING OF A COMPUTER SYSTEM WITH HARDWARE REDUNDANCY SUBJECT TO MAXIMUM REPAIR TIME STOCHASTIC MODELLING OF A COMPUTER SYSTEM WITH HARDWARE REDUNDANCY SUBJECT TO MAXIMUM REPAIR TIME V.J. Munday* Department of Statistics, M.D. University, Rohtak-124001 (India) Email: vjmunday@rediffmail.com

More information