Supply chains models
|
|
- Clifton McGee
- 5 years ago
- Views:
Transcription
1 Supply chains models Maria Ivanova Technische Universität München, Department of Mathematics, Haupt Seminar 2016 Garching, Germany July 3, 2016
2 Overview Introduction What is Supply chain? Traffic flow vs. production systems Basic model Notations Basic model Redefining the flux Constitutive relation Asymptotic validity of the conservation law Scaling Interpolation and weak formulation An exact solution for a single bottleneck Concept of virtual processors Numerical experiments Experiment 1: Verifying Theorem 2 Experiment 2: long supply chain with unstructured throughput times and capacities
3 What is Supply chain? 1. A network of suppliers that produce goods, both, for one another and for generic customers. Goods travel from origin-suppliers to destination-customers, possibly visiting intermediate suppliers and being altered or recombined in the process. Conservation rules at each supplier define its outputs as a function of its inputs. [C.Daganzo, A Theory of Supply Chains]
4 What is Supply chain? 2. A finite, connected directed, simple graph consisting of arcs J = {1,..., N} and vertices V = {1,..., N 1}. Each supplier j is modelled by an arc j, which is again parameterized by an interval [a j, b j ]. We use a 1 = and b N = + for the first respectively the last supplier in the supply chain. [S.Goettlich, M.Herty, A.Klar, Network Models for Supply Chains] We consider a chain of M suppliers or processors S 0,..., S M 1. Each supplier takes a certain good (measured in units or parts), processes it, and hands it over to the next supplier in the chain.
5 Traffic flow vs. production systems Conservation Laws: t ρ + x [φ(ρ)] = 0 ρ: Work in progress (WIP) density vs. density of vehicles; x: stage of the production process (degree of completion - DOC). x = 0 : unfinished product, x = X: finished product vs. traffic flow: road; φ: flux function.
6 Similarities between Traffic and Production Models Complexity and Topology: Complex re - entrant production systems vs. Networks of roads. Control: Policies for production systems vs. Traffic control mechanisms. Random behavior. First Principle Models: Discrete Event Simulation (DES), Multi - Agent Models (incorporate stochastic behavior) Kinetic equations for densities (mean field theories, large time asymptotics) Fluid dynamics Rate equations (fluid models). Simulation optimization and control.
7 Differences between Traffic and Production Modeling Capacity Limits: Limited capacity of processors (machines) vs. limited space (capacity of the road) backward wave propagation in traffic flow φ(ρ) ρ > 0 and φ (ρ) < 0. Parameters and Control Variables: Traffic: randomly given influx and individual behavior, distribute road capacities. Production: randomly given demand, choose start rate and (to some extent) topology.
8 Notations n - part index; m - processor index; τ(m, n) - the time at which part number n passes from supplier number m 1 to supplier number m; U(t) - N-curve given by the number of parts which have passed from processor S m 1 to processor S m at time t, i.e., by U(m, t) = H(t τ(m, n)) n=0 F (m, t) - the flux from processor S m 1 into processor S m F (m, t) = d dt U(m, t) = δ(t τ(m, n)), m = 0,..., M; n=0
9 Notations W (m, t) - the work in progress (WIP) of processor S m at time t W (m, t) = U(m, t) U(m + 1, t) + K(m), m = 0,..., M 1. Combining the last two equations yields the conservation law d W (m, t) = F (m, t) F (m + 1, t) dt T (m) - minimal processing time.
10 Assumptions It holds: τ(m + 1, n) τ(m, n) + T (m) WIP W (m, t) can never become negative. Conservation law: t ρ + x min {µ, W T } = 0 x - artificial continuous variable which indexes the suppliers; µ - service rate; ρ(x, t) - product density over x, i.e., W = ρdx.
11 Basic model Transition times recursion: τ(m + 1, n) = max {τ(m, n) + T (m), τ(m + 1, n 1) + 1 µ(m,n 1) }, m = 0,..., M 1, n 1 We assume finite capacities of the processors, so that µ(m, n) C(m), m = 0,..., M 1, n 0 Initial and boundary conditions for the recursion: τ(0, n) = τ A (n), n 0, τ(m, 0) = τ I (m), m = 0,..., M, where τ A (n) - the arrival time of part n in the first processor; τ I (m) - the time the first part has arrived at supplier S m.
12 Redefining the flux Observe for an arbitrary test function ψ(t) τ I (m) ψ(t)f (x m, t) = n=0 ψ(τ(m, n)) Rewrite into a Riemann sum for an integral τ I (m) where ψ(t)f (x m, t) = ψ(τ(m, n)) n τ(m, n)f(x m, τ(m, n)), n=0 (a) n τ(m, n) := τ(m, n + 1) τ(m, n), (b) f(x m, τ(m, n)) := 1 nτ(m,n) On a time scale where n is small, it becomes τ I (m) ψ(t)f (x m, t) τ I (m) ψ(t)f(x m, t)
13 Redefining the flux Consider the case when the arrival times τ would be distributed continuously, i.e., if they were given as a function τ(x, y). Then we have (a) d dy U(x, τ(x, y)) = tu(x, τ) y τ = 1, (b) d dx U(x, τ(x, y)) = xu(x, τ) + t U(x, τ) x τ. So, for continuum τ(x, y) of arrival times, we set f(x, τ) = 1 yτ and ρ(x, τ) = xτ yτ ρ and f satisfy a conservation law t ρ + x f = 0. Approximate density ρ and flux f defined from the arrival times (a) f(x m, τ(m, n)) = 1, m = 0,..., M, n 0, nτ(m,n) mτ(m,n+1) h m nτ(m+1,n) (b) ρ(x m, τ(m + 1, n)) =, m = 0,..., M 1, n 0, (c) m τ(m, n) = τ(m + 1, n) τ(m, n), h m := x m+1 x m.
14 Constitutive relation THEOREM 1. Let the arrival times τ(m, n) satisfy the transition times recursion. Let the approximate density ρ and flux f be defined as above. Then the approximate flux can be written in terms of the approximate density via a constitutive relation of the form f(x m, τ(m, n)) = φ mn (ρ(x m 1, τ(m, n))), m = 1,..., M, n 0, with the flux function φ m given by φ mn (ρ) = min{µ(m 1, n), h m 1 ρ T (m 1) }.
15 Scaling We define by T 0 the average processing time, i.e., T 0 := 1 M 1 M m=0 T (m) Then the following holds for scaled variables τ(m, n) = MT 0 τ s (m, n), T (m) = T 0 T s (x m ), µ(m, n) = µs(xm,τs(m+1,n)) T 0. Consider a regime where M >> 1 holds and set ε = 1 M << 1. With this scaling the transition times recursion becomes (a) τ s (m + 1, n + 1) = max{τ s (m, n + 1) + εt s (x m ), τ s (m + 1, n) ε + µ s(x m,τ s(m+1,n)) }, m = 0,..., M 1, n 0, (b) τ s (0, n) = τs A (n), n 0, τ s (m, 0) = τs I (m), m = 0,..., M.
16 Scaling We set giving n τ(m, n) = τ(m, n + 1) τ(m, n) = T 0 ns τ s (m, n), m τ(m, n) = τ(m + 1, n) τ(m, n) = T 0 ms τ s (m, n), τ s (m + 1, n) = τ s (m, n) + ε ms τ s (m, n), τ s (m, n + 1) = τ s (m, n) + ε ns τ s (m, n). We scale the density ρ and the flux f by f(x, t) = 1 T 0 f s (x, t MT 0 ), ρ(x, t) = M X ρ s(x, t MT 0 ), where X is the length of the DOC interval. This gives us scaled flux and density (a) f s (x m, τ(m, n)) = 1 (b) ρ s (x m, τ s (m + 1, n)) = nsτ s(m,n) ε msτ(m,n+1) h m nsτ s(m+1,n), m = 0,..., M, n 0,, m = 0,..., M 1, n 0
17 Scaling Scaled version of the constitutive relation then reads (a) f s (x m, τ s (m, n)) = φ s (x m 1, τ s (m, n), ρ s (x m 1, τ s (m, n))), h (b) φ s (x m 1, t, ρ s ) = min{µ s (x m 1, t), m 1 ρ s εxt s(x m 1 ) }. (b) suggests a natural choice for the grid in x-direction, namely h m = εxt s (x m ) = XT (m) M 1 m = 0,..., M 1. m=0 T (m),
18 Interpolation and weak formulation Goal: initial boundary value problem for a conservation law t ρ + x f = 0, f = min{µ(x, t), ρ}, f(0, t) = f A (t) Complications: a domain bounded by t > τ I ; we cannot assume any kind of smooth relation between two consecutive processors; the influx function f can become discontinuous. Solution: Set ρ(x, t) = x u(x, t). This gives t u = min{µ(x, t), x u}, lim u(x, t) = x 0 ga (t), d dt ga (t) = f A (t)
19 Interpolation Interpolation in time direction: (a) f 1 (x m, t) = f(x m, τ(m, n)), τ(m, n) t < τ(m, n + 1); (b) ρ 1 (x m, t) = ρ(x m, τ(m + 1, n 1)), τ(m + 1, n 1) t < τ(m + 1, n); (c) f A (t) = 1 nτ A (n), τ A (n) t < τ A (n + 1). Define: u 1 (x m+1, t) = u 1 (x m, t) hm X ρ 1(x m, t), m = 0,..., M 1, u 1 (x 0, t) = t τ(0,0) f A (s)ds Functions of continuous space and time: (a) f 2 (x, t) = f 1 (x m+1, t), x m x < x m+1, m = 0,..., M 1, (b) τ I 2 (x) = τ I (m + 1), x m x < x m+1, m = 0,..., M 1, (c) u 2 (x, t) = u 1 (x m+1, t), x m x < x m+1, m = 0,..., M 1.
20 Weak formulation THEOREM 2. Given the scaled density and flux at the discrete points x m, τ(m, n), let the piecewise constant interpolant u 2 and f 2 be defined as above. Let the scaled throughput times T (x m ) stay uniformly bounded in m. Assume finitely many bottlenecks for a finite amount of time, i.e., let m τ(m, n) be bounded for ε 0 except for a certain number of nodes m and a finite number of parts n, which stays bounded as ε 0. Then, for ε 0 and max h m 0 the interpolated N-function and flux u 2, f 2 satisfy the initial boundary value problem (a) t u 2 = f 2, t > τ2 I (x), 0 < x < X, (b) u 2 (x, τ I (x)) = 0, lim u 2(x, t) = t x 0 τ 2 (0,0) f A (s)ds, in the limit ε 0, weakly in x and t.
21 An exact solution for a single bottleneck Suppose t ρ + x f = 0, f = min{µ(x, t), ρ}, f(0, t) = f A (t) is posed on x [0, 1], with a bottleneck at x = 1 2, f A - prescribed arrival rate { µ1 for 0 < x < 1 µ(x) = 2, 1 µ 2 for 2 < 0 < 1, µ 2 < f A (t) < µ 1.
22 An exact solution for a single bottleneck The solution is given by ρ(x, t) = ρ c (x, t) + q(t)δ(x 1 2 ), where ρ c satisfies t ρ c + x ρ c = 0, x (0, 1 2 ) ( 1 2, 1), ρ c (0, t) = f A (t), ρ c ( 1 2 +, t) = µ 2, whose solution is given by { f ρ c (x, t) = A (t x) for 0 < x < 1 2, 1 µ 2 for 2 < 0 < 1. In order the whole solution ρ(x, t) to be a spatially weak solution of the conservation law, we have to satisfy 1 0 φ(x) tρ(x, t) min{µ(x), ρ(x, t)} x φ(x)dx = = φ(0)f A (t) φ(1) min{µ 2, ρ(1, t)} for any arbitrary smooth test function φ(x).
23 An exact solution for a single bottleneck Integrating by parts separately on the intervals (0, 1 2 ) and ( 1 2, 1) gives 1 0 φ(x) tρ c (x, t)dx + q (t)φ( 1 2 ) + φ(x) x min{µ(x), ρ c (x, t)}dx min{µ 1, ρ c ( 1 2, t)}φ( 1 2 ) + min{µ 1, ρ c (0, t)}φ(0) + min{µ 1, ρ c (0, t)}φ(0) + min{µ 2, ρ c ( 1 2 +, t)}φ( 1 2 ) = φ(0)f A (t) Since ρ c (x, t) < µ(x) everywhere and ρ c (0, t) = f A (t), this reduces to q (t) = f A (t 1 2 ) µ 2 = 0
24 Concept of virtual processors Basic idea: One processor with a processing time T and a service rate µ can be replaced by K virtual processors with the same service rate µ and processing times T K. THEOREM 3. Let the first processor S 0 in the chain be governed by the transition times recursion. If we replace the single processor by K virtual processors with the same processing rates and the same total throughput time, i.e., by ˆτ(m + 1, n + 1) = max{ˆτ(m, n + 1) + T (0) 1 K, ˆτ(m + 1, n) + µ(m,n) }, m = 0,..., K 1, ˆτ(0, n) = τ(0, n), ˆτ(m, 0) = τ(1, 0) (1 m K )T (0), m = 1,..., K, then we obtain the same outflux, i.e., holds. ˆτ(K, n) = τ(1, n), n 0,
25 Experiment 1: Verifying Theorem 2 Algorithm: Solve the transition times recursion; Compute the WIP W, the N-curve and the flux F ; Compare it to ρ, u and f computed from the solution of the hyperbolic equation. The hyperbolic problem for the approximate N-curve u is solved via a standard finite difference scheme of the form (a) u(x m, t n+1 ) = u(x m, t n ) + tf(x m, t n ), m = 1,..., M, t = t n+1 t n, X min{µ(x m 1, t), MT 0 x m 1 [u(x m, t) u(x m 1, t)]}, (b) f(x m, t) = m = 1,..., M, f A (t n ), m = 0
26 Experiment 1: Verifying Theorem 2 Chain of 3 suppliers with respective nodes S 0, S 1, S 2 ; Throughput times: T (0) = 1, T (1) = 3, T (2) = 1; Capacities C(0) = 15, C(1) = 10, C(2) = 15 parts per time unit; Create the regime of Theorem 2 by introducing virtual processors: split the nodes S 0 and S 2 into 10 virtual nodes and node S 1 into 30 virtual nodes. Set the service rates µ equal to the capacities, giving 15 for 0 < x < 0.2 µ(x, t) = 10 for 0.2 < x < for 0.8 < x < 1
27 Experiment 1: Verifying Theorem 2 We first solve the recursion for the transition times τ, starting with all empty queues, i.e., holds, and set τ I (M) = 0. τ I (m + 1) τ I (m) = T (m) = 0.1 We compute the arrival times randomly according to τ A (n + 1) τ A (n) = 1 f A (τ A (n)), τ A (0) = τ I (0) where f A (t) is a prescribed influx rate.
28 Experiment 1: Verifying Theorem 2
29 Experiment 1: Verifying Theorem 2
30 Experiment 1: Verifying Theorem 2
31 Experiment 1: Verifying Theorem 2
32 Experiment 1: Verifying Theorem 2
33 Experiment 1: Verifying Theorem 2
34 Experiment 2: long supply chain with unstructured throughput times and capacities
35 Experiment 2: long supply chain with unstructured throughput times and capacities
36 Experiment 2: long supply chain with unstructured throughput times and capacities
37 Experiment 2: long supply chain with unstructured throughput times and capacities
38 Experiment 2: long supply chain with unstructured throughput times and capacities
39 Experiment 2: long supply chain with unstructured throughput times and capacities
40 References [D. Armbruster, P. Degond and C. Ringhofer]: A model for the dynamics of large queuing networks and supply chains. SIAM J. Applied Mathematics, 66 (2006), pp [S. Goettlich, M. Herty, A. Klar]: Network models for supply chains. Communication in Mathematical Sciences, 3 (2005), pp [C. F. Daganzo]: A theory of supply chains. Lecture Notes in Economics and Mathematical Systems, 526. Berlin: Springer. viii, 123 p, [C. Ringhofer]: Kinetic and fluid dynamic models for complex supply network. chris/cime09/cime09.htm.
c 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 3, pp. 896 92 c 26 Society for Industrial and Applied Mathematics A MODEL FOR THE DYNAMICS OF LARGE QUEUING NETWORKS AND SUPPLY CHAINS D. ARMBRUSTER, P. DEGOND, AND C.
More informationA Model for the Dynamics of large Queuing Networks and Supply Chains
A Model for the Dynamics of large Queuing Networks and Supply Chains D. Armbruster, P. Degond, C. Ringhofer September 26, 25 Abstract We consider a supply chain consisting of a sequence of buffer queues
More informationHyperbolic Models for Large Supply Chains. Christian Ringhofer (Arizona State University) Hyperbolic Models for Large Supply Chains p.
Hyperbolic Models for Large Supply Chains Christian Ringhofer (Arizona State University) Hyperbolic Models for Large Supply Chains p. /4 Introduction Topic: Overview of conservation law (traffic - like)
More informationA scalar conservation law with discontinuous flux for supply chains with finite buffers.
A scalar conservation law with discontinuous flux for supply chains with finite buffers. Dieter Armbruster School of Mathematical and Statistical Sciences, Arizona State University & Department of Mechanical
More informationTime-dependent order and distribution policies in supply networks
Time-dependent order and distribution policies in supply networks S. Göttlich 1, M. Herty 2, and Ch. Ringhofer 3 1 Department of Mathematics, TU Kaiserslautern, Postfach 349, 67653 Kaiserslautern, Germany
More informationM. HERTY, CH. JÖRRES, AND B. PICCOLI
EXISTENCE OF SOLUTION TO SUPPLY CHAIN MODELS BASED ON PARTIAL DIFFERENTIAL EQUATION WITH DISCONTINUOUS FLUX FUNCTION M. HERTY, CH. JÖRRES, AND B. PICCOLI Abstract. We consider a recently [2] proposed model
More informationSome Aspects of First Principle Models for Production Systems and Supply Chains. Christian Ringhofer (Arizona State University)
Some Aspects of First Principle Models for Production Systems and Supply Chains Christian Ringhofer (Arizona State University) Some Aspects of First Principle Models for Production Systems and Supply Chains
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
More informationA HYPERBOLIC RELAXATION MODEL FOR PRODUCT FLOW IN COMPLEX PRODUCTION NETWORKS. Ali Unver, Christian Ringhofer and Dieter Armbruster
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS Supplement 2009 pp. 790 799 A HYPERBOLIC RELAXATION MODEL FOR PRODUCT FLOW IN COMPLEX PRODUCTION NETWORKS Ali Unver, Christian Ringhofer
More informationA model for a network of conveyor belts with discontinuous speed and capacity
A model for a network of conveyor belts with discontinuous speed and capacity Adriano FESTA Seminario di Modellistica differenziale Numerica - 6.03.2018 work in collaboration with M. Pfirsching, S. Goettlich
More informationThermalized kinetic and fluid models for re-entrant supply chains
Thermalized kinetic and fluid models for re-entrant supply chains D. Armbruster, C. Ringhofer May 19, 2004 Abstract Standard stochastic models for supply chains predict the throughput time (TPT) of a part
More informationOptimization For Supply Chain Models With Policies
Optimization For Supply Chain Models With Policies M. Herty and C. Ringhofer December 7, 26 Abstract We develop a methodology to investigate optimal dynamic policies for a large class of supply networks.
More informationModel hierarchies and optimization for dynamic flows on networks
Model hierarchies and optimization for dynamic flows on networks S. Göttlich and A. Klar Department of mathematics, TU Kaiserslautern Fraunhofer ITWM, Kaiserslautern Collaborators: P. Degond (Toulouse)
More informationScalar conservation laws with moving density constraints arising in traffic flow modeling
Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International
More informationKinetic and fluid models for supply chains supporting policy attibutes
Kinetic and fluid models for supply chains supporting policy attibutes D. Armbruster (1), P. Degond (2), C. Ringhofer (3) (1) Department of Mathematics, Arizona State University, Tempe, USA, AZ 85287-1804,
More informationTraffic flow on networks: conservation laws models. Benedetto Piccoli I.A.C. C.N.R. Rome
Traffic flow on networks: conservation laws models Benedetto Piccoli I.A.C. C.N.R. Rome Conservation laws on networks u t + f(u) x=0 Dynamics at nodes? 1. The only conservation at nodes does not determine
More informationUNIVERSITY OF YORK. MSc Examinations 2004 MATHEMATICS Networks. Time Allowed: 3 hours.
UNIVERSITY OF YORK MSc Examinations 2004 MATHEMATICS Networks Time Allowed: 3 hours. Answer 4 questions. Standard calculators will be provided but should be unnecessary. 1 Turn over 2 continued on next
More informationConservation laws and some applications to traffic flows
Conservation laws and some applications to traffic flows Khai T. Nguyen Department of Mathematics, Penn State University ktn2@psu.edu 46th Annual John H. Barrett Memorial Lectures May 16 18, 2016 Khai
More informationIntersection Models and Nash Equilibria for Traffic Flow on Networks
Intersection Models and Nash Equilibria for Traffic Flow on Networks Alberto Bressan Department of Mathematics, Penn State University bressan@math.psu.edu (Los Angeles, November 2015) Alberto Bressan (Penn
More informationKinetic and fluid models for supply chains supporting policy attributes
Kinetic and fluid models for supply chains supporting policy attributes D. Armbruster (1), P. Degond (2), C. Ringhofer (3) (1) Department of Mathematics, Arizona State University, Tempe, USA, AZ 85287-1804,
More informationModeling of Material Flow Problems
Modeling of Material Flow Problems Simone Göttlich Department of Mathematics University of Mannheim Workshop on Math for the Digital Factory, WIAS Berlin May 7-9, 2014 Prof. Dr. Simone Göttlich Material
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
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 informationMinimizing Total Delay in Fixed-Time Controlled Traffic Networks
Minimizing Total Delay in Fixed-Time Controlled Traffic Networks Ekkehard Köhler, Rolf H. Möhring, and Gregor Wünsch Technische Universität Berlin, Institut für Mathematik, MA 6-1, Straße des 17. Juni
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 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 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 informationABHELSINKI UNIVERSITY OF TECHNOLOGY
Modelling the population dynamics and file availability of a P2P file sharing system Riikka Susitaival, Samuli Aalto and Jorma Virtamo Helsinki University of Technology PANNET Seminar, 16th March Slide
More informationStability 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 informationOptimization of order policies in supply networks
Optimization of order policies in supply networs S. Göttlich M. Herty C. Ringhofer August 18, 2008 Abstract The purpose of this paper is to develop a model which allows for the study and optimization of
More informationAnswers to Problem Set Number 02 for MIT (Spring 2008)
Answers to Problem Set Number 02 for 18.311 MIT (Spring 2008) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139). March 10, 2008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics,
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Monday March 12, 2018. Turn it in (by 3PM) at the Math.
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 informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE IV: HOMOGENIZATION OF KINETIC FLOWS IN NETWORKS. Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE IV: HOMOGENIZATION OF KINETIC FLOWS IN NETWORKS. Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More informationTraffic models on a network of roads
Traic models on a network o roads Alberto Bressan Department o Mathematics, Penn State University bressan@math.psu.edu Center or Interdisciplinary Mathematics Alberto Bressan (Penn State) Traic low on
More informationCumulative Count Curve and Queueing Analysis
Introduction Traffic flow theory (TFT) Zhengbing He, Ph.D., http://zhengbing.weebly.com School of traffic and transportation, Beijing Jiaotong University September 27, 2015 Introduction Outline 1 Introduction
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Espen R. Jakobsen Phone: 73 59 35 12 Examination date: December 16, 2017 Examination
More informationTime Reversibility and Burke s Theorem
Queuing Analysis: Time Reversibility and Burke s Theorem Hongwei Zhang http://www.cs.wayne.edu/~hzhang Acknowledgement: this lecture is partially based on the slides of Dr. Yannis A. Korilis. Outline Time-Reversal
More informationMathematical modeling of complex systems Part 1. Overview
1 Mathematical modeling of complex systems Part 1. Overview P. Degond Institut de Mathématiques de Toulouse CNRS and Université Paul Sabatier pierre.degond@math.univ-toulouse.fr (see http://sites.google.com/site/degond/)
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 22 12/09/2013. Skorokhod Mapping Theorem. Reflected Brownian Motion
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 22 12/9/213 Skorokhod Mapping Theorem. Reflected Brownian Motion Content. 1. G/G/1 queueing system 2. One dimensional reflection mapping
More informationCONTINUUM MODELS FOR FLOWS ON COMPLEX NETWORKS
CONTINUUM MODELS FOR FLOWS ON COMPLEX NETWORKS C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OUTLINE 1 Introduction: Network Examples: Some examples of real world networks + theoretical structure.
More informationModeling, Validation and Control of Manufacturing Systems
Modeling, Validation and Control of Manufacturing Systems E. Lefeber, R.A. van den Berg and J.E. Rooda Abstract In this paper we elaborate on the problem of supply chain control in semiconductor manufacturing.
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 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 informationFlow-level performance of wireless data networks
Flow-level performance of wireless data networks Aleksi Penttinen Department of Communications and Networking, TKK Helsinki University of Technology CLOWN seminar 28.8.08 1/31 Outline 1. Flow-level model
More informationMacroscopic Simulation of Open Systems and Micro-Macro Link
Macroscopic Simulation of Open Systems and Micro-Macro Link Ansgar Hennecke 1 and Martin Treiber 1 and Dirk Helbing 1 II Institute for Theoretical Physics, University Stuttgart, Pfaffenwaldring 57, D-7756
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 informationModeling and Analysis of Manufacturing Systems
34 Modeling and Analysis of Manufacturing Systems E. Lefeber Eindhoven University of Technology J.E. Rooda Eindhoven University of Technology 34.1 Introduction... 34-1 34. Preliminaries... 34-34.3 Analytical
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 informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
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 informationExam in TMA4195 Mathematical Modeling Solutions
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Exam in TMA495 Mathematical Modeling 6..07 Solutions Problem a Here x, y are two populations varying with time
More informationMathematical Modeling. Preliminary Lecture Notes
Mathematical Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date April 8, 2013 2 Contents 1 Preface 5 2 Modeling Process 7 2.1 Constructing Models......................... 7 2.1.1 Conservation
More informationNICTA Short Course. Network Analysis. Vijay Sivaraman. Day 1 Queueing Systems and Markov Chains. Network Analysis, 2008s2 1-1
NICTA Short Course Network Analysis Vijay Sivaraman Day 1 Queueing Systems and Markov Chains Network Analysis, 2008s2 1-1 Outline Why a short course on mathematical analysis? Limited current course offering
More informationTraffic Flow Problems
Traffic Flow Problems Nicodemus Banagaaya Supervisor : Dr. J.H.M. ten Thije Boonkkamp October 15, 2009 Outline Introduction Mathematical model derivation Godunov Scheme for the Greenberg Traffic model.
More informationL 1 stability of conservation laws for a traffic flow model
Electronic Journal of Differential Equations, Vol. 2001(2001), No. 14, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu ftp ejde.math.unt.edu (login:
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationOptima and Equilibria for Traffic Flow on Networks with Backward Propagating Queues
Optima and Equilibria for Traffic Flow on Networks with Backward Propagating Queues Alberto Bressan and Khai T Nguyen Department of Mathematics, Penn State University University Park, PA 16802, USA e-mails:
More information2D Traffic Flow Modeling via Kinetic Models
Modeling via Kinetic Models Benjamin Seibold (Temple University) September 22 nd, 2017 Benjamin Seibold (Temple University) 2D Traffic Modeling via Kinetic Models 09/22/2017, ERC Scale-FreeBack 1 / 18
More informationContinuum Modelling of Traffic Flow
Continuum Modelling of Traffic Flow Christopher Lustri June 16, 2010 1 Introduction We wish to consider the problem of modelling flow of vehicles within a traffic network. In the past, stochastic traffic
More informationOn the distribution schemes for determining flows through a merge
On the distribution schemes for determining flows through a merge W. L. Jin and H. M. Zhang April 11, 2002 Abstract: In this paper, we study various distribution schemes for determining flows through a
More informationMetadata of the chapter that will be visualized in SpringerLink
Metadata of the chapter that will be visualized in SpringerLink Book Title Series Title Chapter Title Copyright Year 212 Copyright HolderName Decision Policies f Production Netwks The Production Planning
More informationExact Simulation of the Stationary Distribution of M/G/c Queues
1/36 Exact Simulation of the Stationary Distribution of M/G/c Queues Professor Karl Sigman Columbia University New York City USA Conference in Honor of Søren Asmussen Monday, August 1, 2011 Sandbjerg Estate
More informationA.Piunovskiy. University of Liverpool Fluid Approximation to Controlled Markov. Chains with Local Transitions. A.Piunovskiy.
University of Liverpool piunov@liv.ac.uk The Markov Decision Process under consideration is defined by the following elements X = {0, 1, 2,...} is the state space; A is the action space (Borel); p(z x,
More informationy + p(t)y = g(t) for each t I, and that also satisfies the initial condition y(t 0 ) = y 0 where y 0 is an arbitrary prescribed initial value.
p1 Differences Between Linear and Nonlinear Equation Theorem 1: If the function p and g are continuous on an open interval I : α < t < β containing the point t = t 0, then there exists a unique function
More informationQuality of Real-Time Streaming in Wireless Cellular Networks : Stochastic Modeling and Analysis
Quality of Real-Time Streaming in Wireless Cellular Networs : Stochastic Modeling and Analysis B. Blaszczyszyn, M. Jovanovic and M. K. Karray Based on paper [1] WiOpt/WiVid Mai 16th, 2014 Outline Introduction
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
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 informationRandom Walk on a Graph
IOR 67: Stochastic Models I Second Midterm xam, hapters 3 & 4, November 2, 200 SOLUTIONS Justify your answers; show your work.. Random Walk on a raph (25 points) Random Walk on a raph 2 5 F B 3 3 2 Figure
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationA random perturbation approach to some stochastic approximation algorithms in optimization.
A random perturbation approach to some stochastic approximation algorithms in optimization. Wenqing Hu. 1 (Presentation based on joint works with Chris Junchi Li 2, Weijie Su 3, Haoyi Xiong 4.) 1. Department
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationOn a diophantine inequality involving prime numbers
ACTA ARITHMETICA LXI.3 (992 On a diophantine inequality involving prime numbers by D. I. Tolev (Plovdiv In 952 Piatetski-Shapiro [4] considered the following analogue of the Goldbach Waring problem. Assume
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 informationMathematical Modeling. Preliminary Lecture Notes
Mathematical Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date March 23, 2012 2 Contents 1 Preface 5 2 Modeling Process 7 2.1 Constructing Models......................... 7 2.1.1 Conservation
More informationDelay of Incidents Consequences of Stochastic Incident Duration
Delay of Incidents Consequences of Stochastic Incident Duration Victor L. Knoop 1, Serge P. Hoogendoorn 1, and Henk J. van Zuylen 1 Delft University of Technology & TRIL research School, Stevinweg 1, 68
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 informationThe Riemann Solver for Traffic Flow at an Intersection with Buffer of Vanishing Size
The Riemann Solver for Traffic Flow at an Intersection with Buffer of Vanishing Size Alberto Bressan ( ) and Anders Nordli ( ) (*) Department of Mathematics, Penn State University, University Par, Pa 16802,
More informationMIT Manufacturing Systems Analysis Lectures 18 19
MIT 2.852 Manufacturing Systems Analysis Lectures 18 19 Loops Stanley B. Gershwin Spring, 2007 Copyright c 2007 Stanley B. Gershwin. Problem Statement B 1 M 2 B 2 M 3 B 3 M 1 M 4 B 6 M 6 B 5 M 5 B 4 Finite
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationMath 180C, Spring Supplement on the Renewal Equation
Math 18C Spring 218 Supplement on the Renewal Equation. These remarks supplement our text and set down some of the material discussed in my lectures. Unexplained notation is as in the text or in lecture.
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationModelling, Simulation & Computing Laboratory (msclab) Faculty of Engineering, Universiti Malaysia Sabah, Malaysia
1.0 Introduction Intelligent Transportation Systems (ITS) Long term congestion solutions Advanced technologies Facilitate complex transportation systems Dynamic Modelling of transportation (on-road traffic):
More informationPartial Differential Equations Summary
Partial Differential Equations Summary 1. The heat equation Many physical processes are governed by partial differential equations. temperature of a rod. In this chapter, we will examine exactly that.
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationANALYSIS OF THE LAW OF THE ITERATED LOGARITHM FOR THE IDLE TIME OF A CUSTOMER IN MULTIPHASE QUEUES
International Journal of Pure and Applied Mathematics Volume 66 No. 2 2011, 183-190 ANALYSIS OF THE LAW OF THE ITERATED LOGARITHM FOR THE IDLE TIME OF A CUSTOMER IN MULTIPHASE QUEUES Saulius Minkevičius
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationOptimal control of a linearized continuum model for re-entrant manufacturing production systems
Optimal control of a linearized continuum model for re-entrant manufacturing production systems Xiaodong Xu 1 and Stevan Dubljevic 1 1 Department of Chemical and Materials Engineering, University of Alberta,
More informationEstimating the Loss of Efficiency due to Competition in Mobility on Demand Markets
Estimating the Loss of Efficiency due to Competition in Mobility on Demand Markets Thibault Séjourné 1 Samitha Samaranayake 2 Siddhartha Banerjee 2 1 Ecole Polytechnique 2 Cornell University November 14,
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More informationBRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 15: Crump-Mode-Jagers processes and queueing systems with processor sharing
BRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 5: Crump-Mode-Jagers processes and queueing systems with processor sharing June 7, 5 Crump-Mode-Jagers process counted by random characteristics We give
More informationOne dimensional Maps
Chapter 4 One dimensional Maps The ordinary differential equation studied in chapters 1-3 provide a close link to actual physical systems it is easy to believe these equations provide at least an approximate
More informationControl of Continuum Models of Production Systems
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 55, NO. 11, NOVEMBER 2010 2511 Control of Continuum Models of Production Systems Michael La Marca, Dieter Armbruster, Member, IEEE, Michael Herty, and Christian
More information