A model for a network of conveyor belts with discontinuous speed and capacity
|
|
- Ophelia Shepherd
- 5 years ago
- Views:
Transcription
1 A model for a network of conveyor belts with discontinuous speed and capacity Adriano FESTA Seminario di Modellistica differenziale Numerica work in collaboration with M. Pfirsching, S. Goettlich
2 Introduction
3 Motivations Conveyor Belts - Simple scenario - Transport equation Micro-meso-macro models - Simple interaction rules Networks and discontinuities - Complex behaviors Basic issue for network control, output tracking, etc. Figure: Inside a brewery 1 a pressureless-accumulation device. 1 Image courtesy of Sidel Blowing & Services SAS.
4 Basic model A production system is described by a conservation law with discontinuous flux function and constant speed a > 0. More precisely for a closed set Ω R and calling ρ : Ω [0, T ] [0, ] the product density, the evolution of the system is { t ρ(x, t) + x (ah( ρ(x, t))ρ(x, t)) = 0, ρ(x, 0) = ρ 0 (x), where H is the Heaviside function and the initial data ρ 0 is a BV function s.t. ρ 0 (x).
5 On a single arc: On a single arc the equation simply reduces to an advection equation and therefore the solution is simply ρ(x, t) = ρ 0 (x at). The presence of discontinuities in the flux function, anyway, poses some problems in the case of ρ 0 (x) = for some x Ω.
6 On a single arc: On a single arc the equation simply reduces to an advection equation and therefore the solution is simply ρ(x, t) = ρ 0 (x at). The presence of discontinuities in the flux function, anyway, poses some problems in the case of ρ 0 (x) = for some x Ω.
7 Weak solutions: Obtained using some Friedrichs mollifiers (classic). We denote by f the following multivalued function f (ρ) = aρ if ρ, f ( ) = [0, a ]. Definition A function ρ L (R [0, T ]) is called weak solution to the Cauchy problem (4) if there exists a function v L (R [0, T ]) such that v(x, t) f (ρ) a.e. and T 0 Ω ρ φ T t dxdt + v φ 0 Ω x dxdt + ρ 0 (x)φ(x, 0)dx = 0 Ω for each φ C 1 c (R [0, T ]) (where φ C 1 c means φ C 1 with compact support).
8 Entropy Weak Solutions: Denote by H the following multivalued function H(ρ) = H(ρ) if ρ 0, H(0) = [0, 1]. Definition A weak solution ρ of the Cauchy problem (4) is called an entropy weak solution if, for each entropy η C 1 (R), η convex, there exists a function w L (R [0, T ]) such that w(x, t) H(ρ(x, t)) a.e. and t η(ρ) + x F (ρ) η ( ) w x 0, where F (ρ) = a ρ 0 η (s)h( s)ds.
9 One arc case We remark that that the solution ρ(x, t) = ρ 0 (x at). is a weak entropy solution. This can be shown by choosing v(x, t) = aρ f (ρ), for ρ [0, ] and (Kruskov Entropy functions) η(ρ) = ρ k, for ρ [0, ], w(x, y) 1 for any constant k R. We observe that Entropy Definition reduces to the advection equation with constant speed a andthe function above provides the solution to the problem for every initial condition ρ 0 (x).
10 Analytic discussion of the Network case
11 Extension to networks: We now extend the model to a network Γ = (V, E) with V the set of vertices and E V V the set of arcs. For a fixed node v V, the sets δv, δ v + denote the ingoing and outgoing arcs, respectively. We consider the following network problem: t ρ(x, t) + x (f i(ρ(x, t))) = 0, x Ω i, t [0, T ] ρ(x, 0) = ρ 0 i (x), x Ω i f i (ρ(v, t)) = f i (ρ(v, t)), v V, where i δ v i δ + v f i (ρ(x, t)) = a i H( i ρ(x, t))ρ(x, t), x Ω i, t [0, T ] is the flux function on arc i E and a i R + the transport velocity.
12 One-to-one junction We consider the problem on Ω = Ω 1 {0} Ω 2 = (, 0) {0} (0, ), where the intersection is located at x = 0. The model equations are given by (10), where the junction condition is specified as f 1 (ρ(0, t)) = f 2 (ρ(0, t)) with flux function (10) on arc i = 1, 2.
13 The solution can be easily computed if no congestion occurs during the transportation, i.e. a 1 ρ 0 1(x) a 2 2 x Ω 1. For x Ω 1, the characteristics of the problem are simply the straight lines y(t) = x a 1 t, which lead to the solution Analogously, we find ρ(x, t) = ρ 0 1(x a 1 t) for (x, t) : x < 0. ρ(x, t) = ρ 0 2(x a 2 t) for { (x, t) Ω 2 (0, T ] x > a 2 t } which corresponds to the density initially placed on Ω 2. The solution for (x, t) in the case 0 < x < a 2 t is found by f 1 (ρ(x +, t)) = f 2 (ρ(x, t)) where x ± denotes lim h 0 x ± h. This leads to ρ(x, t) = a 1 a 2 ρ(x +, t) at the interface x = 0.
14 t x = a 2 t + c x = a 1 t + c 0 Figure: Characteristics in the non-congested case The solution is then described by ρ 0 1 (x a 1t), )) (x, t) Ω 1 (0, T ] a ρ(x, t) = 1 a 2 ρ 0 1 ( a 1 (t x, (x, t) Ω a2 2 (0, T ] : 0 < x a 2 t ρ 0 2 (x a 2t), (x, t) Ω 2 (0, T ] : x > a 2 t. x
15 If a 1 ρ 0 1 (x) a 2 2 is not satisfied, congestion arises and the problem is more complex. Let t 0 denote the first time of congestion: { t 0 = inf t 0 such that ρ 0 1( a 1 t) > a } 2 2. a 1 We track the interface describing the congested area (Λ). The interface is a time-dependent function g(t). We can track g(t) by integrating the difference between the fluxes entering and exiting the region Λ as well as the current density. t t 0 0 ( a1 ρ 0 1( a 1 s) a 2 ) 2 ds = g(t) ( 1 ρ 0 1(y a 1 t))dy.
16 Rearranging the terms, we can describe the congested region Λ as { } Λ := (x, t) Ω 1 [t 0, t0 E ], such that g(t) x 0 with interface g defined as g : [0, ) R 0, where t is mapped to the solution of x + (t t 0 ) ρmax 2 1 a 2 1 if this is negative, and to zero otherwise. The final time of congestion is t E = min {t t 0 such that g(t) = 0}. 1 a 1 t 0 x a 1 t ρ 0 1(s)ds = 0,
17 t t E 0 x = a 2 t + c x = a 1 t + c g(t) Λ 0 Trajectories of the problem with congestion (a 1 > a 2 ). Inside Λ, the transport velocity ā is t 0 2 ā = a 2 1. x
18 Ze derive a general solution discussing the Riemann problem on [0, T ) R with the initial data { ρl, x < 0, ρ(x, 0) = ρ r, x 0. 1 ρ 2 ρ l ρ r Shock, s l < 0 x = 0 Shock, s r > 0 Shock, s = 0 x
19 We allow waves with negative velocity and we distinguish: a) f 1 f2 max : no congestion arises standard advection. b) f 1 > f2 max : a congestion arises the solution consists of three shock waves starting at x = 0: i) At x = 0, a shock with velocity s = 0 arises, where the solution ρ jumps from 1 to 2. For the special case 1 = 2, there is no jump in the solution at x = 0. ii) Left-going shock wave, where the density jumps from ρ l to 1 with negative velocity s l < 0 (Rankine-Hugoniot c.): s l = āρmax 1 f 1 (ρ l ) = f max 2 f 1 (ρ l ) 1 ρ l < 0. 1 ρ l This shock wave describes the left boundary of the congested area Λ, meaning that it follows the function g(t). iii) Right-going shock wave, where the density jumps from 2 to ρ r with positive velocity s r = a 2 > 0 (Rankine-Hugoniot c.): s r = f 2(ρ r ) f 2 ( 2 ) ρ r 2 = a2ρ r a 2 ρmax 2 ρ r 2 = a 2 > 0.
20 In the general case with non-constant initial data, we obtain the following solution: ρ 0 1 (x a 1t), (x, t) Ω 1 \ Λ (0, T ] 1, (x, t) Λ (0, T ] 2, (x, t) Ω 2 (0, T ] ρ(x, t) = )) x a 2 t, g(t x a 2 ) 0 a 1 a 2 ρ 0 1 ( a 1 (t x, (x, t) Ω a2 2 (0, T ] x a 2 t, g(t x a 2 ) = 0 ρ 0 2 (x a 2t), (x, t) Ω 2 (0, T ] : x > a 2 t. If g(t) = 0 for all t [0, T ], we have Λ = and recover the previously described solution. Contrary to traffic flow models, rarefaction waves do not appear in the conveyor belt problem.
21 One-to-two junction We denote by i = 1 the incoming and by i = 2, 3 the outgoing arcs. We consider two kind of intersections: a) if the intersection is passive, we impose a fixed rate between the outgoing fluxes. b) if the intersection is active the ratio changes in order to reduce congestion. a) f 1 b) f 1 D f 2 f 3 f 2 f 3
22 We consider the problem on Ω = Ω 1 e 1 {0} Ω 2 e 1 Ω 3 e 2 = (, 0) e 1 {0} (0, ) e 1 (0, ) e 2, where (e 1, e 2 ) is the standard base of R 2, and we identify the element x Ω 1 with the vector (x, 0) T, and analogously for arcs i = 2, 3. The equations that we consider are as before, where the junction condition is specified as f 1 (ρ(0, t)) = f 2 (ρ(0, t)) + f 3 (ρ(0, t)), equipped with the discontinuous flux function on arc i = 1, 2, 3.
23 Case 1: passive junction. An intersection device D keeps the ratio of the two outgoing fluxes constant, i.e., for a fixed distribution parameter µ [0, 1], it holds f 2 = µf 1, f 3 = (1 µ)f 1, even if only one outgoing arc is congested. No congestion arises iff { 1 ρ 0 a 2 1(x) min µ a 1 2, } 1 a 3 3, x Ω 1. 1 µ a 1 We have the actual density on arc 2 and 3 as { } { } µ a 3 1 µ ρ 2 = min 3, a 2 2, ρ 3 = min 2, 3. 1 µ a 2 µ a 3
24 The first time of congestion t 0 can be determined as { { 1 t 0 = inf t 0 such that ρ 0 a 2 1( a 1 t) > min µ a 1 We can introduce the interface g(t). The interface is g : [0, ) R 0 mapped to the solution of 2, }} 1 a µ a 1 g(t) + (t t 0 ) ρmax 2 1 a a 1 t 0 g(t) a 1 t ρ 0 1(s)ds = 0, if this is negative, and to zero otherwise. The region of congestion Λ on Ω 1 is given by { } Λ := (x, t) Ω 1 [t 0, t0 E ], such that g(t) x 0, analogously to the 1-to-1 case.
25 The general solution ρ on Ω (passive junction) is given by ρ 0 1 (x a 1t), (x, t) Ω 1 \ Λ (0, T ] 1, (x, t) Λ (0, T ] ρ i, (x, t) Ω i (0, T ] : x a i t, g(t x a ρ(x, t) = i ) 0, i = 2, 3 α i a 1 a i ρ 1 0 ( a 1 (t x ai )), (x, t) Ω i (0, T ] : x a i t, g(t x a i ) = 0, i = 2, 3 ρ 0 i (x a it), (x, t) Ω i (0, T ] : x > a i t, i = 2, 3, where α 2 = µ and α 3 = 1 µ for a compact notation.
26 Case 2: active junction. If congestion arises, the device adapts the fluxes to ensure the maximal total outgoing flux, i.e., f 2 = f2 max and f 3 = f3 max. We fix a parameter µ [0, 1]. We have (no congestion) µf 1 = f 2 and (1 µ)f 1 = f 3. To get a unique solution also in the congested case, we define β i with the following properties: a) The flux conservation f 1 = f 2 + f 3 at the coupling is satisfied. b) The parameters β i are equal to α i, i = 2, 3 if no congestion. c) If arc i is congested, i.e., ρ 0 1 (x a 1t) > a i α i a 1 i parameter β i changes to β i = a i a 1 d) If the value β i = a i i a 2 2 +a 3 3 i ρ 0 1 ( a 1t). at time t the is reached the congestion starts.
27 No congestion arises as long as the following inequality holds true: ρ 0 1(x) a 2 a a 3 a 1 3, x Ω 1. If not, the time t 0 (first time of congestion) is { t 0 = inf t 0 such that ρ 0 1( a 1 t) > a a } 3 3. a 1 a 1 The interface g is the solution of ( ρ max 2 g(t) + (t t 0 ) 1 a 2 + ρmax 3 1 a 3 ) 1 1 a 1 t 0 g(t) a 1 t ρ 0 1(s)ds = 0, if this solution is negative, and to zero otherwise. The region of congestion Λ on Ω 1 is { } Λ := (x, t) Ω 1 [t 0, t0 E ], such that g(t) x 0.
28 We define for i = 2, 3 and α 2 = µ, α 3 = 1 µ, { { β i (t) := min max α i, a i } i a i } i a 1 ρ 0 1 (x a, 1t) a a 3 3 The general solution ρ on Ω (active junction) is ρ 0 1 (x a 1t), (x, t) Ω \ Λ (0, T ] 1, (x, t) Λ (0, T ] i, (x, t) Ω i (0, T ] x a i t, g(t x a ρ(x, t) = i ) 0 )) β i (t x a 1 ) a 1 a i ρ 1 0 ( a 1 (t x, (x, t) Ω ai i (0, T ] x a i t, g(t x a i ) = 0 with i = 2, 3. ρ 0 i (x a it), (x, t) Ω i (0, T ] x > a i t
29 Two-to-one junction We know from traffic flow that in the free flow regime no additional information is needed. Conversely, in the congested case, we need a priority rule between the two incoming arcs, i.e., how to use released capacities of the outgoing arc. We denote by i = 1, 2 the incoming and by i = 3 the outgoing arcs. We consider the problem on Ω = Ω 1 e 1 Ω 2 e 2 {0} Ω 3 e 1 = (, 0) e 1 (, 0) e 2 {0} (0, ) e 1 with the same interpretation as before, where the system is given by conservation laws as before, the coupling condition reads as f 1 (ρ(0, t)) + f 2 (ρ(0, t)) = f 3 (ρ(0, t))..
30 The solution can be directly computed if it holds a 1 ρ 0 1( a 1 t) + a 2 ρ 0 2( a 2 t) a 3 3, i.e. no congestion arises. Then, the solution is given by ρ 0 i (x a it), (x, t) Ω i (0, T ], i = 1, 2 2 )) a ρ(x, t) = i a 3 ρ 0 i ( a i (t x, (x, t) Ω a3 3 (0, T ] : x a 3 t i=1 ρ 0 3 (x a 3t), (x, t) Ω 3 (0, T ] : x > a 3 t and the incoming mass can be totally absorbed by the outgoing arc. This is independent of the capacity of the incoming arcs, and no further priority rule is needed to obtain a unique solution.
31 In the congested case, to obtain a unique solution, we impose a priority rule. We want to use the whole capacity of the outgoing arc i = 3, which implies f 3 = f3 max = a 3 3. We set the merging parameter q [0, 1] such that This leads to f 1 = qf max 3 and f 2 = (1 q)f max 3. f 2 = 1 q q f 1, describing a ratio of the actual fluxes on the corresponding arcs.
32 a) f max 2 f 2 f 1 + f 2 = f max 3 Θ P f max 1 f 2 = q 1 q f 1 f2 max f 1 b) f 2 f 1 + f 2 = f max 3 Figure: Choice of the merging parameter q Θ f 2 = q Q f max 1 P 1 q f 1 f 1 The admissible region for the fluxes is Θ = {(f 1, f 2 ) : 0 f 1 f max 1, 0 f 2 f max 2, 0 f 1 + f 2 f max 3 } and shaded gray in the Figure above.
33 Merging parameters: a) The intersection point P between the maximal outgoing flux (the line f 1 + f 2 = f3 max ) and the priority ratio is inside the admissible set Θ. In this case, we keep q (0, 1) fix. b) The intersection point P is outside Θ. We choose the closest point Q inside Θ on the line f 1 + f 2 = f3 max, which guarantees maximal throughput, i.e., f 1 = f1 max. The merging parameter changes to q = f1 max /f3 max. Then, the resulting fluxes are f 1 = f1 max and f 2 = f3 max f1 max.
34 Λ i is the congested region on arc i = 1, 2, g i (t) its interface ) solution of g i (t) + (t t 0 ) a 3 1 a i t 0 ρ 0 i (s)ds = 0 ( qi 3 i i g i (t) a i t and q i {q, 1 q} the corresponding merging parameters. The solution is then (for i = 1, 2) ρ 0 i (x a it), (x, t) Ω i \ Λ i (0, T ], i, (x, t) Λ i (0, T ], 3, (x, t) { Ω 3 (0, T ] ( ) } x a 3 t, max g i t x ρ(x, t) = i=1,2 a 3 0 )) a i a 3 ρ 0 i ( a i (t x, (x, t) Ω a3 3 (0, T ] i=1,2 { ( ) } x a 3 t, max g i t x i=1,2 a 3 = 0 ρ 0 3 (x a 3t), (x, t) Ω 3 (0, T ] : x > a 3 t.
35 Numerical Tools
36 Numerical approximation We discretize Ω i = (ā i, b i ) X i = (x i,0 = ā i, x i,1,..., x i,mi = b i ) with constant discretization x = x i,j x i,j 1, j = 1,..., m i. The spatial grid cells are defined as C i,j = (x i,j 1/2, x i,j+1/2 ) R. We discretize [0, T ] with {t n = n t, n = 0, 1,..., T / t}, where t [0, T ]. We define the piecewise constant approximation of the solution ρ as ρ n i,j ρ(x i,j, t n ) with ρ n constant in each grid cell C i,j. The scheme to update the approximation in each time step is for all i E, j 2,..., m i 1 given by { ( ρ n+1 i,j = ρ n i,j x t h ( ρ n i,j, i,j+1) ρn h ( ρ n i,j 1, ρ n i,j) ), ρ 0 i,j = ρ0 (x i,j ). For j = 1 and all i E.
37 The inflow and outflow fluxes must be adapted at the junction a) For an one-to-one junction, i.e. δ v = δ + v = 1, we set h n,out i = h n,in î = h ( ρ n i,m i, ρ ṋ i,1) for i δ v, î δ + v. b) For an one-to-two junction, i.e. 1 = δ v δ + v = 2, we choose in the non-congested case for i δ v, î δ + v h n,out i = f ( ρ n i,m i ) for i δ v, and h n,in î b.i) In the congested case: on a passive junction { h n,in 2 = min f2 max, h n,out i = î δ + v µ 1 µ f 3 max h n,in î for i δ v, }, h n,in 3 = min = αîf i ( ρ n i,m i ). { f3 max, 1 µ } µ f 2 max.
38 b.ii) On an active junction h n,out i = î δ + v f max î for i δ v, and h n,in î = f max î for î δ + v. c) For a two-to-one junction, i.e. 2 = δv δ v + = 1, in the non-congested case = h n,out i = f i ( ρ n i,m i ) for i δ v, and h n,in î i δ v f i ( ρ n i,m i ) for î δ + v. In the congested case, we apply the merging parameter q i and set h n,out i = q i f max î for i δ v, î δ + v, and h n,in î = f max î for î δ + v.
39 To define h( ρ n i,j, ρn i,j+1 ) we look for a function h satisfying m (ũ) ũ h(0, 0) = h( i, i ) = 0, h(ũ, u) 0 u h(ũ, u) m +(u), with continuous function m : R R and m = min(m, 0), m + = max(m, 0).
40 Numerical Mollifiers Since f i is discontinuous in ρ, we need a suitable regularization. We define a Friedrichs mollifier ϕ C0 (R) ϕ( y) = ϕ(y), ϕ(y)dy = 1. In our case, we use the mollifier ϕ(y) := max(0, 1 y ) and define ϕ ξ (y) := 2 ξ ϕ( 2y ξ ) for a small parameter ξ > 0, which implies that ϕ ξ (y) has compact support in [ ξ 2, ξ 2 ]. For each arc i E, we introduce the following smooth regularization of the flux function f ξ,i (ρ) := a i ρ ( 1 ρ i R ( ϕ ξ y i ξ ) ) dy 2,
41 f ξ,i (ρ) ā ρ a i i 0 i ρ i ρ + ξ The function coincides with the original flux function in x [0, i ] but there is a continuously differentiable connection to the value f ξ,i ( i + ξ) = 0.
42 We choose the numerical flux function h as Godunov flux min f ξ,i(z), if ρ n h( ρ n i,j, ρ n z [ ρ i,j+1) = n i,j, ρn i,j+1 ] i,j ρn i,j+1 max f ξ,i(z), if ρ n z [ ρ n i,j+1, ρn i,j ] i,j ρn i,j+1. The scheme is stable, if the following CFL condition is fulfilled. x t max v V δ v m L (0, +ξ) In particular, if max δ v 2 for a fixed v V, the CFL condition reduces to t x ξ 4.
43 Tests
44 We compute the numerical solution based on scheme described before. The space step size is x = , the time step size t = 10 5 and the smoothing parameter is fixed to ξ = We choose = 1. One-to-one junction The linear network is given by Ω 1 = ( π, 0) and Ω 2 = (0, π), i.e. the coupling is at x = 0. We fix the initial solution ρ 0 on Ω 1 Ω 2 as ( ρ 0 (x) = exp 3 (x + 3 ) ) 2 5 π.
45 t = 0 t = 1 t = 1.4 t = 1.8 t = 2.2 t = 2.6 Figure: Test 1: non-congested case with a 1 = 1 and a 2 = 2
46 One-to-one congested case t = 0 t = 0.5 t = 0.8 t = 1.1 t = 1.4 t = 1.7 Figure: Test 2: congested case with a 1 = 2 and a 2 = 1
47 Numerical solution x t error Analytical solution ξ error Table: space-time diagram for the congested case (above) and decreasing step sizes (left), decreasing smoothing parameter ξ (right)
48 One-to-two passive junction with µ = 0.5 t = 0 : t = 0.5 : t = 1.0 :
49 One-to-two active junction with µ = 0.5 t = 0 : t = 0.5 : t = 1.0 :
50 Two-to-one junction with q = 0.3 t = 0 : t = 2 : t = 4 :
51 Further steps: Control problems related to this model 1 Tracking of a production output (K. Lux SCICOM) 2 Optimization of the network Adding noise and stochastic disturbance 1 Appearance of viscosity terms 2 Evolution prediction in complex networks
52 D. Armbruster, S. Göttlich, and M. Herty, A scalar conservation law with discontinuous flux for supply chains with finite buffers, SIAM J. Appl. Math., 71 (2011), pp A. Festa, M. Pfirsching, and S. Göttlich, A model for a network of conveyor belts with discontinuous speed and capacity, ArXiv (2018). M. Garavello, K. Han, and B. Piccoli, Models for vehicular traffic on networks, vol. 9, American Institute of Mathematical Sciences (AIMS), Springfield, MO, S. Göttlich, A. Klar, and P. Schindler, Discontinuous conservation laws for production networks with finite buffers, SIAM J. Appl. Math., 73 (2013), pp M. Herty, C. Joerres, and B. Piccoli, Existence of solution to supply chain models based on pde with discontinuous flux function, J. Math. Anal. Appl., 401 (2013), pp Thank you!
M. 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 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 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 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 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 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 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 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 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 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 informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
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 informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationMULTI CLASS TRAFFIC MODELS ON ROAD NETWORKS
COMM. MATH. SCI. Vol. 4, No. 3, pp. 59 68 c 26 International Press MULTI CLASS TRAFFIC MODELS ON ROAD NETWORKS M. HERTY, C. KIRCHNER, AND S. MOUTARI Abstract. We consider a multi class (resp. source destination)
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationNumerical discretization of tangent vectors of hyperbolic conservation laws.
Numerical discretization of tangent vectors of hyperbolic conservation laws. Michael Herty IGPM, RWTH Aachen www.sites.google.com/michaelherty joint work with Benedetto Piccoli MNCFF 2014, Bejing, 22.5.2014
More informationThe flux f(ρ) is given by v ρ where v is the average speed of packets among nodes, derived considering the amount of packets that may be lost.
A FLUID DYNAMIC MODEL FOR TELECOMMUNICATION NETWORKS WITH SOURCES AND DESTINATIONS CIRO D APICE, ROSANNA MANZO AND BENEDETTO PICCOLI Abstract. This paper proposes a macroscopic fluid dynamic model dealing
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 informationSupply chains models
Supply chains models Maria Ivanova Technische Universität München, Department of Mathematics, Haupt Seminar 2016 Garching, Germany July 3, 2016 Overview Introduction What is Supply chain? Traffic flow
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
More informationA MACROSCOPIC TRAFFIC MODEL FOR ROAD NETWORKS WITH A REPRESENTATION OF THE CAPACITY DROP PHENOMENON AT THE JUNCTIONS
A MACROSCOPIC TRAFFIC MODEL FOR ROAD NETWORKS WITH A REPRESENTATION OF THE CAPACITY DROP PHENOMENON AT THE JUNCTIONS B. Haut G. Bastin Y. Chitour Aspirant FNRS, haut@auto.ucl.ac.be, Centre for Systems
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
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 informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationarxiv: v2 [math.na] 11 Jul 2014
Runge-Kutta Discontinuous Galerkin Method for Traffic Flow Model on Networks Suncica Canic 1, Benedetto Piccoli, Jing-Mei Qiu 3, Tan Ren 4 arxiv:1403.3750v [math.na] 11 Jul 014 Abstract. We propose a bound-preserving
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 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 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 informationRegularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws
1 Regularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws Sébastien Blandin, Xavier Litrico, Maria Laura Delle Monache, Benedetto Piccoli and Alexandre Bayen Abstract
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University. Email: wxs27@psu.edu November 5, 2015 Abstract We study two systems of conservation
More informationNumerical Approximations of a Traffic Flow Model on Networks
Numerical Approximations of a Traffic Flow Model on Networks Gabriella Bretti, Roberto Natalini, Benedetto Piccoli Abstract We consider a mathematical model for fluid-dynamic flows on networks which is
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 informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
More informationThe LWR model on a network
Mathematical Models of Traffic Flow (October 28 November 1, 2007) Mauro Garavello Benedetto Piccoli DiSTA I.A.C. Università del Piemonte Orientale C.N.R. Via Bellini, 25/G Viale del Policlinico, 137 15100
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationStabilité des systèmes hyperboliques avec des commutations
Stabilité des systèmes hyperboliques avec des commutations Christophe PRIEUR CNRS, Gipsa-lab, Grenoble, France Séminaire EDP, Université de Versailles 9 octobre 2014 1/34 C. Prieur UVSQ, Octobre 2014 2/34
More informationOn finite time BV blow-up for the p-system
On finite time BV blow-up for the p-system Alberto Bressan ( ), Geng Chen ( ), and Qingtian Zhang ( ) (*) Department of Mathematics, Penn State University, (**) Department of Mathematics, University of
More informationSOURCE-DESTINATION FLOW ON A ROAD NETWORK
COMM. MATH. SCI. Vol. 3, No. 3, pp. 261 283 c 2005 International Press SOURCE-DESTINATION FLOW ON A ROAD NETWORK MAURO GARAVELLO AND BENEDETTO PICCOLI Abstract. We construct a model of traffic flow with
More informationarxiv: v1 [physics.soc-ph] 1 Aug 2017
OFF-RAMP COUPLING CONDITIONS DEVOID OF SPURIOUS BLOCKING AND RE-ROUTING NAJMEH SALEHI, JULIA SOMERS, AND BENJAMIN SEIBOLD arxiv:1708.00529v1 [physics.soc-ph] 1 Aug 2017 Abstract. When modeling vehicular
More informationLARGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMERICS FOR 1-D SCALAR CONSERVATION LAWS
LAGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMEICS FO 1-D SCALA CONSEVATION LAWS L. I. IGNAT, A. POZO, E. ZUAZUA Abstract. In this paper we analyze the large time asymptotic behavior of the discrete
More informationIndex. Chromatography, 3 Condition of C 0 compatibility, 8, 9 Condition of C 1 compatibility, 8, 9, 33, 36, 38, 39, 41, 43, 50, 53, 56, 58, 59, 62
References 1. Bressan, A., Čanić, S., Garavello, M., Herty, M., Piccoli, B.: Flows on networks: recent results and perspectives. EMS Surv. Math. Sci. 1, 47 111 (2014) 2. Coron, J.-M., Wang, Z.: Controllability
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationMULTICOMMODITY FLOWS ON ROAD NETWORKS
MULTICOMMODITY FLOWS ON ROAD NETWORKS M. HERTY, C. KIRCHNER, S. MOUTARI, AND M. RASCLE Abstract. In this paper, we discuss the multicommodity flow for vehicular traffic on road networks. To model the traffic,
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University Abstract We study two systems of conservation laws for polymer flooding in secondary
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
More informationMathematical Modelling of Traffic Flow at Bottlenecks
Mathematical Modelling of Traffic Flow at Bottlenecks Cathleen Perlman Centre for Mathematical Sciences Lund Institute of Technology Advisor: Stefan Diehl June 10, 2008 Abstract This master s thesis gives
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationSTABILITY ESTIMATES FOR SCALAR CONSERVATION LAWS WITH MOVING FLUX CONSTRAINTS. Maria Laura Delle Monache. Paola Goatin
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX STABILITY ESTIMATES FO SCALA CONSEVATION LAWS WITH MOVING FLUX CONSTAINTS Maria Laura Delle Monache Department
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 informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationADER schemes and high order coupling on networks of hyperbolic conservation laws
ADER schemes and high order coupling on networks of hyperbolic conservation laws Raul Borsche a,, Jochen Kall a a Erwin Schrödinger Straße, TU Kaiserslautern, Building 48, 67663 Kaiserslautern, Germany
More informationOn a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws
On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws Stefano Bianchini and Alberto Bressan S.I.S.S.A., Via Beirut 4, Trieste 34014 Italy. E-mail addresses: bianchin@mis.mpg.de,
More informationWEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 94, pp. 1 14. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu WEAK ASYMPTOTIC
More informationStabilité des systèmes hyperboliques avec des commutations
Stabilité des systèmes hyperboliques avec des commutations Christophe PRIEUR CNRS, Gipsa-lab, Grenoble, France Contrôles, Problèmes inverses et Applications Clermont-Ferrand Septembre 2014 1/34 C. Prieur
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 253 Monotonicity of Implicit Finite Difference Methods for Hyperbolic Conservation Laws Michael
More informationThe 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy
The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,
More informationSimulation of diffusion. processes with discontinuous coefficients. Antoine Lejay Projet TOSCA, INRIA Nancy Grand-Est, Institut Élie Cartan
Simulation of diffusion. processes with discontinuous coefficients Antoine Lejay Projet TOSCA, INRIA Nancy Grand-Est, Institut Élie Cartan From collaborations with Pierre Étoré and Miguel Martinez . Divergence
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationA GENERAL PHASE TRANSITION MODEL FOR VEHICULAR TRAFFIC
A GENEAL PHASE TANSITION MODEL FO VEHICULA TAFFIC S. BLANDIN, D. WOK, P. GOATIN, B. PICCOLI, AND A. BAYEN Abstract. An extension of the Colombo phase transition model is proposed. The congestion phase
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationHysteresis rarefaction in the Riemann problem
Hysteresis rarefaction in the Riemann problem Pavel Krejčí 1 Institute of Mathematics, Czech Academy of Sciences, Žitná 25, 11567 Praha 1, Czech Republic E-mail: krejci@math.cas.cz Abstract. We consider
More informationModeling and Numerical Approximation of Traffic Flow Problems
Modeling and Numerical Approximation of Traffic Flow Problems Prof. Dr. Ansgar Jüngel Universität Mainz Lecture Notes (preliminary version) Winter Contents Introduction Mathematical theory for scalar conservation
More informationRadon measure solutions for scalar. conservation laws. A. Terracina. A.Terracina La Sapienza, Università di Roma 06/09/2017
Radon measure A.Terracina La Sapienza, Università di Roma 06/09/2017 Collaboration Michiel Bertsch Flavia Smarrazzo Alberto Tesei Introduction Consider the following Cauchy problem { ut + ϕ(u) x = 0 in
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
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 informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationSOME PRELIMINARY NOTES ON HYPERBOLIC CONSERVATION LAWS
SOME PRELIMINARY NOTES ON HYPERBOLIC CONSERVATION LAWS KAYYUNNAPARA T. JOSEPH CENTRE FOR APPLICABLE MATHEMATICS TATA INSTITUTE OF FUNDAMENTAL RESEARCH BANGALORE 5665 1. Introduction Many of the physical
More informationMarte Godvik. On a Traffic Flow Model. Thesis for the degree philosophiae doctor
Marte Godvik On a Traffic Flow Model Thesis for the degree philosophiae doctor Norwegian University of Science and Technology Faculty for Information Technology, Mathematics and Electrical Engineering
More informationNumerical simulation of some macroscopic mathematical models of traffic flow. Comparative study
Numerical simulation of some macroscopic mathematical models of traffic flow. Comparative study A. Lakhouili 1, El. Essoufi 1, H. Medromi 2, M. Mansouri 1 1 Hassan 1 st University, FST Settat, Morocco
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationNumerical Simulation of Traffic Flow via Fluid Dynamics Approach
International Journal of Computing and Optimization Vol. 3, 2016, no. 1, 93-104 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijco.2016.6613 Numerical Simulation of Traffic Flow via Fluid Dynamics
More informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More informationResearch Article Splitting of Traffic Flows to Control Congestion in Special Events
International Journal of Mathematics and Mathematical Sciences Volume 0, Article ID 567, 8 pages doi:0.550567 Research Article Splitting of Traffic Flows to Control Congestion in Special Events Ciro D
More informationDynamics of Propagation and Interaction of Delta-Shock Waves in Conservation Law Systems
Dynamics of Propagation and Interaction of Delta-Shock Waves in Conservation Law Systems V. G. Danilov and V. M. Shelkovich Abstract. We introduce a new definition of a δ-shock wave type solution for a
More informationIntroduction to hyperbolic PDEs and Shallow Water Equations
Introduction to hyperbolic PDEs and Shallow Water Equations Technische Universität München Fundamentals of Wave Simulation - Solving Hyperbolic Systems of PDEs Supervisor Leonhard Rannabauer 8/12/2017
More informationNUMERICAL METHOD FOR THE COMPUTATION OF TANGENT VECTORS TO 2 2 HYPERBOLIC SYSTEMS OF CONSERVATION LAWS
NUMERICAL METHOD FOR THE COMPUTATION OF TANGENT VECTORS TO HYPERBOLIC SYSTEMS OF CONSERVATION LAWS MICHAEL HERTY AND BENEDETTO PICCOLI Abstract. We are interested in the development of a numerical method
More informationSimultaneous vs. non simultaneous blow-up
Simultaneous vs. non simultaneous blow-up Juan Pablo Pinasco and Julio D. Rossi Departamento de Matemática, F..E y N., UBA (428) Buenos Aires, Argentina. Abstract In this paper we study the possibility
More informationNon-Oscillatory Central Schemes for a Traffic Flow Model with Arrhenius Look-Ahead Dynamics
Non-Oscillatory Central Schemes for a Traffic Flow Model with Arrhenius Look-Ahead Dynamics Alexander Kurganov and Anthony Polizzi Abstract We develop non-oscillatory central schemes for a traffic flow
More informationVarious lecture notes for
Various lecture notes for 18306. R. R. Rosales (MIT, Math. Dept., 2-337) April 14, 2014 Abstract Notes for MIT s 18.306 Advanced PDE with applications. Contents 1 First lecture. 3 Tomography. Use filtering
More informationVISCOSITY SOLUTIONS OF HAMILTON JACOBI EQUATIONS WITH DISCONTINUOUS COEFFICIENTS
Dept. of Math. Univ. of Oslo Pure Mathematics No. 2 ISSN 86 2439 January 25 VISCOSITY SOLUTIONS OF HAMILTON JACOBI EQUATIONS WITH DISCONTINUOUS COEFFICIENTS GIUSEPPE MARIA COCLITE AND NILS HENRIK RISEBRO
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
More informationRegularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws
Regularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws Sébastien Blandin, Xavier Litrico, Maria Laura Delle Monache, Benedetto Piccoli and Alexandre Bayen Abstract
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 informationGlobal Existence of Large BV Solutions in a Model of Granular Flow
This article was downloaded by: [Pennsylvania State University] On: 08 February 2012, At: 09:55 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
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 informationIntroduction to nonlinear wave models
Introduction to nonlinear wave models Marko Nedeljkov Department of Mathematics and Informatics, University of Novi Sad Trg Dositeja Obradovića 4, 21000 Novi Sad, Serbia marko.nedeljkov@dmi.uns.ac.rs http://www.dmi.uns.ac.rs
More information