Modeling and Numerical Approximation of Traffic Flow Problems
|
|
- Rachel Long
- 5 years ago
- Views:
Transcription
1 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 laws 4 3 Traffic flow models 3 4 Numerical approximation of scalar conservation laws 6 5 The Gudonov method 7 eferences 3 Introduction The goal of these lecture notes is to introduce to the modeling of car traffic flow using conservation laws and to the numerical discretization of the resulting hyperbolic differential equations. Consider the traffic flow of cars on a highway with only one lane (i.e., overtaking is impossible). Instead of modeling the cars individually, we use the density ρ(x, t) of cars (in vehicles per kilometer, say) in x at time t. The number of cars which are in the interval (x, x ) at time t is x x ρ(x, t) dx.
2 Let v(x, t) denote the velocity of the cars in x at time t. The number of cars which pass through x at time t (in unit length) is ρ(x, t)v(x, t). We want to derive an equation for the evolution of the car density. The number of cars in the interval (x, x ) changes according to the number of cars which enter or leave this interval (see Figure.): d dt x x ρ(x, t) dx = ρ(x, t)v(x, t) ρ(x, t)v(x, t). Integration of this equation with respect to time and assuming that ρ and v are regular functions yields t x t x t ρ(x, t) dx dt = = t (ρ(x, t)v(x, t) ρ(x, t)v(x, t)) dx dt t t x Since x, x, t, t > are arbitrary, we conclude t x x (ρ(x, t)v(x, t)) dx dt. ρ t + (ρv) x =, x, t >. (.) This is a partial differential equation. It has to be supplemented by the initial condition ρ(x, ) = ρ (x), x. (.) ρ (x,t) v(x,t) cars ρ (x,t) v(x,t) Figure.: Derivation of the conservation law. We also need an equation for the velocity v. We assume that v only depends on ρ (see Section 3 for other choices). If the highway is empty (ρ = ), we will drive with maximal velocity v = v max ; in heavy traffic we will slow down and will stop (v = ) in a tailback where the cars are bumper to bumper (ρ = ρ max ). The simplest model is the linear relation ( v(ρ) = v max ρ ), ρ ρ max. ρ max Equation (.) then becomes [ ( ρ t + v max ρ ρ )] =, ρ max x x, t >. (.3)
3 This equation is a so-called conservation law since it expresses the conservation of the number of cars. Indeed, integrating (.3) formally over x gives [ ( )] d ρ(x, t) ρ(x, t) dx = v max ρ(x, t) dx =, dt x ρ max and thus the number of cars in is constant for all t. Equation (.3) belongs to the class of hyperbolic equations. We call a system of equations t u + x f(u) =, x, t >, (.4) u(x, ) = u (x), x, (.5) with f : m m hyperbolic if and only if for any u m, f (u) can be diagonalized and has only real eigenvalues. A function u : [, ) m is called a classical solution of (.4) (.5) if u C ( (, )) C ( [, )) and u solves (.4) (.5) pointwise. Equation (.3) can be simplified by bringing it in a dimensionless form. Let L and τ be a typical length and time, respectively, such that L/τ = v max. Introducing we obtain x s = x L, t s = t τ, u = ρ ρ max, t ρ = [ τ ρmax ] t s ( u) = ρ max τ t s u, ] ( x [v max ρ ρ )] = L [ ρ ρ max xs v max max ( u) ( + u) and hence we can write (.) (.3) as (with (x, t) instead of (x s, t s )) = ρ max τ x s ( u ( ) u u t + x =, x, t >, (.6) u(x, ) = u (x), x, (.7) with u = ρ /ρ max. If the highway is empty (ρ = ), we have u = ; in a tailback (ρ = ρ max ), u = holds. Equation (.6) is also called the inviscid Burgers equation. Let In order to see how the solutions of (.6) may look like, consider the following example. : x < u (x) = x : x < : x, ), (.8) i.e., there are initially no cars in x < and a moderate traffic in x. What happens for t >? It is easy to check that u(x, t) = : x < t t x t : t x < t : x t 3
4 is a solution of (.6) (.7) for x, t < / which satisfies (.4) pointwise except at x = and x = t (see Figure.). The solution becomes discontinuous at t =. t = t = t = 4 x Figure.: Solution of (.6) (.8). The cars which were initially in (, ) drive from the left to the right until no cars were remaining in that interval. As u = in x <, no other cars are coming. This gives rise to the following questions: Does there exist any solution for t /? If yes, how does the solution for t / look like? How can we solve (.6) (.7) for more complicated intial data numerically? It is clear that we need some theory before we can turn to the numerical treatment of hyperbolic conservation laws. Mathematical theory for scalar conservation laws In this section we study the problem u t + f(u) x =, x, t >, (.) u(x, ) = u (x), x, (.) for some function f :. This problem can be solved using the method of characteristics. 4
5 Definition. Let u : [, T ) be a (classical) solution of (.) (.). The solutions χ of the initial-value problem are called the characteristics of (.). χ (t) = f (u(χ(t), t)), t >, χ() = x The main property of the characteristics is that u is constant along them: d dt u(χ(t), t) = u t(χ(t), t) + u x (χ(t), t)χ (t) =, t >, and hence u(χ(t), t) = const. for t >. Example. Let f(u) = u /. The characteristics of (.) are given as the solutions of χ (t) = u(χ(t), t), t >, χ() = x. As u(χ(t), t) = const. for t >, χ is a straight line in the (x, t)-plane through x with slope /u (χ(t)) = /u (x ). Characteristics allow to illustrate the solutions of (.6) (.7) in a compact form. An example with u as in (.8) is given in Figure.. The characteristics should not be confused with the vehicle trajectories; they rather illustrate the propagation of the density values. x Figure.: Characteristics of (.6) (.8). Example. shows again that solutions of (.) (.) may develop discontinuities after finite time. Therefore we need a solution concept including discontinuous functions. For 5
6 this, let u be a classical solution of (.) (.); multiplying (.) by φ C ( ) = {φ C ( ) : φ has compact support} and integrating over yields = (u t + f(u) x ) φ dx dt = (uφ t + f(u)φ x ) dx dt u(x, )φ(x, ) dx. In order to define the last two integrals we only need integrable functions u! This motivates the following definition. Definition.3 The function u : (, T ) is called a weak solution of (.) (.) if for all φ C ( ) (uφ t + f(u)φ x ) dx dt = u (x)φ(x, ) dx. Implicitely, this definition includes some regularity requirements on u (for instance, integrability of u and f(u)), but we do not specify them. It is not difficult to check that each classical solution is a weak solution, but the inverse does not need to be true. Another weak formulation of (.) can be obtained as follows. Let u be a classical solution of (.) (.). Integrating (.) over (a, b) (s, t) for any a, b and s, t > gives b a u(x, t) dx b a u(x, s) dx = t s f(u(b, τ)) dτ + t s f(u(a, τ)) dτ. (.3) Then we can say that u is a weak solution of (.) (.) if u(, ) = u and (.3) is satisfied for all a, b, s, t >. It is possible to prove that any weak solution in the sense of Definition.3 also satisfies (.3), but we do not present the (technical) proof. We now consider conservation laws with special discontinuous initial data. Definition.4 The problem (.) (.) with initial datum { u l : x < u (x) = u r : x (.4) and u l, u r is called a iemann problem. We want to solve the iemann problem (.) (.), (.4). For this, we observe that with u(x, t) also u(αx, αt) is a solution of (.) (.), (.4) for any α >. Therefore, u only depends on ξ = x/t, i.e. u = u(ξ). This implies = u t + f(u) x = x t u (ξ) + f (u(ξ))u (ξ) t = t u (ξ)(f (u(ξ)) ξ). Now several possibilities can occur: 6
7 u (ξ) = = u(ξ) = const. f (u(ξ)) = ξ = u(ξ) = (f ) (ξ) (if the inverse of f exists; a sufficient condition is f > in ). u is discontinuous along ξ = x/t, i.e., u (ξ) does not exist. This observation motivates to consider three cases: Case : u l = u r. This gives the solution u(x, t) = u r = u l for all x, t >. Case : u l > u r. For f(u) = u /, a traffic flow interpretation is that the vehicle density in x > is larger than in x <. As the flow direction is from the left to the right, we expect a shock line, i.e. a discontinuity curve x = ψ(t). discountinuous function u(x, t) = { u l : x < st u r : x st We claim that the (.5) is a weak solution of (.) (.), (.4). Then the discontinuity line is given by x = ψ(t) = st and s = ψ (t) is the shock speed which has to be determined. In order to prove our claim let φ C ( ). Then, since u = const. except on x = st, uφ t dx dt = = and, by integration by parts, f(u)φ x dx dt = ( st = uφ t dx + ( t st + t st st ) uφ t dx dt uφ dx s u(st, t)φ(st, t) ) uφ dx + s u(st +, t)φ(st, t) dt u(x, )φ(x, ) dx s (u l u r ) ( st st = (f(u l ) f(u r )) We conclude (uφ t + f(u)φ x ) dx dt = φ(st, t) dt f(u) x φ dx + f(u(st, t))φ(st, t) ) f(u) x φ dx f(u(st +, t))φ(st, t) dt u (x)φ(x, ) dx φ(st, t) dt. [s (u l u r ) (f(u l ) f(u r ))] φ(st, t) dt. 7
8 Thus, choosing we have proved the claim. s = f(u l) f(u r ) u l u r, (.6) The choice (.6) is called the ankine-hugoniot condition. For a iemann problem, s is always a constant, i.e., the discontinuity curve is always a straight line. This may not be true for general initial data. In this situation, the ankine-hugoniot condition generalizes to [4, Sec. 6] s(t) = ψ (t) = f(u l(t)) f(u r (t)), (.7) u l (t) u r (t) where u l (t) = lim u(x, t), u r(t) = lim u(x, t). (.8) x ψ(t) x ψ(t) It can be shown that (.5) is the unique weak solution of (.) (.) (see Theorem.9). Example.5 Let f(u) = u / and u l =, u r =. Then the shock speed is s = (u l + u r ) =, and the solution of (.) (.), (.4) is illustrated in Figure.. The traffic flow interpretation is the following. In x < st = t/, cars with a moderate density are coming from the left and stop at the tailback at x = t/. In x > t/, the density is maximal and the cars are bumper to bumper and are not moving. x = t/ denotes the beginning of the tailback. The shock line t x Figure.: Characteristics for (.6) (.7) with u l =, u r =. 8
9 Case 3: u l < u r. For this case we assume that f > in. One solution is given by (.5), as the above proof does not need a sign on u l u r : { u l : x < st u (x, t) = u r : x st. It can be shown that also u l : x < f (u l )t u (x, t) = (f ) (x/t) : f (u l )t x f (u r )t u r : x > f (u r )t is a weak solution (motivated by the computation before Case ; see Figure.3). In fact, it is possible to show that the problem (.) (.), (.4) possesses infinitely many weak solutions! (See, for instance, [5, Sec. 3.5].) What is the physically relevant solution? We present two approaches. t t x x Figure.3: Characteristics for (.) (.), (.4) with f(u) = u / and u l =, u r =, corresponding to u (left) and u (right). In the traffic flow interpretation, the condition u l < u r means that there are more cars (per kilometer) in {x < } than in {x > }. The solution u would mean that all cars in {x < st} drive with the same velocity v(u l ) = v max ( + u l )/, whereas all drivers in {x > st} move with the velocity v(u r ) > v(u l ), and there is a shock at x = st. It would be more realistic if the drivers just before the shock line (x < ) tried to drive as fast as the drivers behind the shock (x > ). After some time there are drivers with velocities v(u l ) and v(u r ) far away from x = and some drivers with velocities between v(u l ) and v(u r ) near x =. Thus, u seems to be the preferable physical solution. The solution u is called rarefaction wave. Is it possible to formulate a general principle which allows to select the physically correct solution? The answer is yes and this leads to the notion of entropy condition. 9
10 Definition.6 A weak solution u : (, T ) of (.) (.), (.4) satisfies the entropy condition of Oleinik if and only if along each discontinuity curve x = ψ(t), f(u l (t)) f(v) u l (t) v ψ (t) f(u r(t)) f(v) u r (t) v (.9) for all t (, T ) and for all v between u l (t) and u r (t), where u l (t), u r (t) are defined in (.8). Does u satisfy the entropy condition (.9)? Since ψ (t) = s(t) = f(u r) f(u l ) u r u l and f is assumed to be strictly convex, we obtain for any u l < v < u r f(u l ) f(v) u l v < f(u l) f(u r ) u l u r = s < f(u r) f(v), u r v which contradicts (.9). Thus, u does not satisfy the entropy condition (.9). The function u of Case, defined in (.5), however, satisfies (.9) (if f is convex). As the function u is continuous, we do not need to check (.9) for this function. The second approach uses the notion of entropy. We call a function η C () an entropy and ψ C () an entropy flux if and only if η is strictly convex and if for any classical solution u of (.) (.) it holds η(u) t + ψ(u) x =, x, t >. (.) The idea of the second approach is to consider the conservation law as an idealization of a diffusive problem given by the equation u t + f(u) x = εu xx, x, t >, (.) where ε > is the diffusion coefficient. This equation, together with the initial condition (.) has a unique smooth solution u ε (by the theory of parabolic differential equations), and we assume that u ε u pointwise in (, T ) for ε, η (u ε )u ε,x L ( (,T )) c, (.) where the constant c > is independent of ε. The limit ε is called the vanishing viscosity limit. It is possible to prove that u is a solution of (.) (.), and we say that u is the physically relevant solution. Then something happens with the entropy equation (.). We multiply (.) by η (u ε ) and choose ψ = f η : η(u ε ) t + ψ(u ε ) x = εη (u ε )u ε,xx = ε(η (u ε )u ε,x ) x εη (u ε )u ε,x.
11 Multiplying this equation by φ C( ), φ, and integrating over (, ) gives: (η(u ε ) t + ψ(u ε ) x ) φ dx dt = ε η (u ε )u ε,x φ x dx dt ε η (u ε )u ε,x φ dx dt ε η (u ε )u ε,x L ( (, )) φ x L ( (, )) (as ε ), using η (u ε ) > and (.). As φ is arbitrary, we deduce the entropy inequality η(u) t + ψ(u) x. Clearly, this only holds for smooth solutions. From the definition of weak solutions follows that the entropy inequality for weak solutions writes (η(u) φ t + ψ(u) φ x )dx dt η(u (x))φ(x, ) dx φ C( ). (.3) We define: Definition.7 Let u : (, T ) be a weak solution of (.) (.). Then u is called an entropy solution if and only if for all convex entropies η and corresponding entropy fluxes ψ, the inequality (.3) holds. The function u satisfies the entropy equation (.) almost everywhere since u is continuous and we can define the derivatives in a weak sense (it is also possible to prove (.) in the weak form similarly to (.3)). Thus, u is an entropy solution. Does this also hold for u? We give the (negative) answer in the following example. Example.8 Let f(u) = u /, η(u) = u (hence, ψ(u) = (/3)u 3 ), and let φ C( ), φ. Then, since s = (u l + u r )/, (u φ t + 3 u3 φ x )dx dt = ( t st u φ dx su lφ(st, t) + t u φ dx + su rφ(st, t) + 3 u3 lφ(st, t) ) 3 u3 rφ(st, t) dt = u (x) φ(x, ) dx (u l + u r )(u l u r) st + 3 (u3 l u 3 r) φ(st, t) dt = u (x) φ(x, ) dx + 6 (u l u r ) 3 φ(st, t) dt η(u (x))φ(x, ) dx φ(st, t) dt
12 if and only if u l u r. Thus, u is not an entropy solution. Example.8 shows that the two equations ( ) u u t + = and (u ) t + 3 (u3 ) x = are equivalent only for classical solutions. x The above considerations motivate that only rarefaction waves u are the physically relevant solutions of the iemann problem (.) (.), (.4) if u l < u r. For u l > u r we have to expect discontinuous solutions with shocks. We summarize the above results in the following theorem. Theorem.9 Let f C () with f > in. () Let u l > u r and set s = (f(u l ) f(u r ))/(u l u r ). Then { u l : x < st u(x, t) = u r : x > st is a weak solution of (.) (.), (.4) satisfying the entropy condition (.9) of Oleinik. () Let u l < u r. Then u l : x < f (u l )t u(x, t) = (f ) (x/t) : f (u l )t x f (u r )t u r : x > f (u r )t is a weak entropy solution of (.) (.), (.4). Above we have written that the problem (.) (.), (.4) with u l < u r possesses infinitely many solutions and that the solution u does neither satisfy the entropy condition of Oleinik nor is an entropy solution. However, is u the only solution satisfying (.9) and (.3)? The answer is affirmative but not easy to see, so we only cite the result (see [6]): Theorem. Let f C () and u L (). Then there exists at most one entropy solution of (.) (.) satisfying the entropy condition (.9). 3 Traffic flow models We present some simple traffic flow models. The first model has been already presented in Section :
13 () Lighthill-Whitham-ichards model: ( ρ t + (ρv(ρ)) x =, v(ρ) = v max ρ ), ρ ρ max. ρ max In Section we have shown that this equation can be simplified by means of the transformation u = ρ/ρ max to u t + (u ) x =. (3.) () Greenberg model: In this model it is assumed that the velocity of the vehicles can be very large for low densities: This implies ρ t + (ρv(ρ)) x =, v(ρ) = v max ln ρ max ρ, < ρ ρ max. ρ t v max (ρ ln ρ) x =. (3) Payne-Whitham model: ρ t + (ρv) x =, (ρv) t + (ρv + p(ρ)) x =. This model mimics the flow of gas particles. In fact, the above equations are known as the Euler equations of gas dynamics with pressure p(ρ) = aρ γ, a >, γ. The disadvantage of this model is that there may be solutions for which the velocity v is negative [3]. (4) Aw-ascle model: ρ t + (ρv) x =, (ρv + ρp(ρ)) t + (ρv + ρvp(ρ)) x =. This model has been proposed as an improvement of the Payne-Whitham model []. It has been derived from microscopic models []. In the following, we will consider the Lighthill-Whitham-ichards model in some detail. More precisely, we want to study the influence of a temporal disturbance due to a traffic interruption (for instance, a traffic light). Suppose that the highway is initially filled with cars with uniform density ρ in the interval (, ), that there is a red light at x =, and that the highway is empty in (, ). We work with the transformed equation (3.) and the uniform density u = ρ/ρ max. Clearly, there will be a tailback in front of the traffic light. At time t = ω, the traffic light changes from red to green, and 3
14 the cars move from the left to the right. We want to know what happens with the tailback at time t > ω. First step: red phase ( t < ω). We solve the Burgers equation (3.) for x < with the initial function u (x) = u, x <, and boundary condition u(, t) = (which models the red traffic light). The solution is (see Figure 3.) { u : x < st u(x, t) = : x > st, where x <, < t < ω, and the shock speed s equals s = u l + u r = u. The solution for x > is given by assumption by u(x, t) =, x >, < t < ω. t t green phase u= u = u u= w red phase u= st u = u sw u= u= x Figure 3.: Characteristics for the traffic flow problem. Second step: green phase (t ω). We solve the Burgers equation (3.) in, with initial datum u (x) = u(x, ω). Since u l = u > = u r we get a shock ψ(t) = st, t > ω, with speed s = (u )/ <. u l = < = u r. The solution is given by u : x < st, : st x < ω t, u(x, t) = x : ω t x t ω, t ω : x > t ω, A rarefaction wave develops at x =, since 4
15 x, t > ω. This solution makes sense as long as st < ω t or, equivalently, t < t := ω/(s + ) = ω/(u + ). Third step: green phase (t > t ). What happens with the shock for t > t? The shock speed is given by the generalized anking-hugoniot condition (.7): s(t) = ψ (t) = (u(ψ(t) +, t) + u(ψ(t), t)) = ( ) ψ(t) t ω + u, t > t. This is a linear ordinary differential equation with inital condition ψ(t ) = st = ω u u +. The solution is ψ(t) = u(t ω) t ω ω( u ), t t. Now there are two cases. If u, ψ(t) as t, i.e., there is a shock for all time and the shock moves into the negative x-direction with shock speed ψ (t) u (t ). The velocity of the shock line is thus reduced from s = (u )/ to u s. Clearly, the discontinuity step at ψ(t) goes to zero as t. However, the drivers observe the shock even a long time after the traffic disturbance. This agrees with practical experience. ψ(t) t ψ(t) t t slope /s slope / u t t x x Figure 3.: Shock curve ψ(t) for u (left) and u > (right). In the other case, u >, ψ(t) as t, i.e., the shock line moves into the positive x-direction. Suppose that after the time t = ω the traffic light changes from green to red. How long should the green phase be to eliminate the shock, i.e., we need t ω where ψ(t ) =. The equation ψ(t ) = has the (unique) solution t = ω/u. 5
16 Thus, t ω is satisfied if and only if u /. The tailback in front of the traffic light will vanish during the green phase if and only if u / or, in original density, ρ ρ := ρ max ( ).46ρ max, independently of the length of the green phase! Even for moderate vehicle density ρ > ρ, the tailback in front of the traffic light will become larger and larger as time increases. 4 Numerical approximation of scalar conservation laws Before studying numerical methods for nonlinear scalar equations we start with a very simple linear equation: u t + au x =, x, t >, (4.) u(x, ) = u (x), x, (4.) where a >. This problems has the explicit solution u(x, t) = u (x at) which is a weak solution (at least if u is smooth enough). We discretize the (x, t)-plane by the mesh (x i, t n ) with x i = ih (i Z), t n = nk (n N ) and h, k >. For simplicity of presentation we take a uniform mesh with h and k constant, but the discussed methods can be easily extended to non-uniform meshes. We are looking for finite difference approximations u n i to the solution u(x i, t n ) at the discrete grid points. The idea is to replace the partial derivatives in (4.) by difference quotients. For example, (4.) can be written by Taylor expansion as u(x i, t n+ ) u(x i, t n ) k which motivates the first numerical scheme: + O(k) = a u(x i+) u(x i ) h Central scheme: Consider the approximation of (4.3) i u n i k We can write this scheme as i = a un i+ u n i, n, i Z. h + O(h ), (4.3) = u n i ak h (un i+ u n i ). (4.4) 6
17 As we can compute i from the data u n i explicitely, this is an explicit scheme. Another idea would be to use the scheme or, equivalently, i ak u n i k h un+ i+ + un+ i = a un+ i+ un+ i h ak h un+ i = un i. This is an implicit scheme. In each time step, a linear system has to be solved. As for time-dependent hyerbolic equations, implicit schemes are rarely used, we consider in the following only explicit schemes. In practice we must compute on a finite spatial domain, say x Nh, and we require appropriate boundary conditions at x =, x = N h, respectively. For the above equation we choose periodic boundary conditions or, for the approximations u n i, u(, t) = u(nh, t), t >, u n = u n N, n. We can determine u n and u n N using (4.4). In fact, setting i = or i = N in (4.4) would require to determine u n or u n N+, We consider these points as artificial points with u n = u n N and un N = un, by the periodicity. Therefore, our central scheme reads as follows u i = u (x i ), i =,..., N, i = u n i ak h (un i+ u n i ), i =,..., N, = u n ak h (un u n N ), N = u n N ak h (un u n N ). Thus, if u n = u n N then un+ = N, and the periodicity is preserved. This scheme has the advantage that it is mass-preserving, i.e. N i= u n i = N i= u (x i ), n. 7
18 Indeed, N i= i = = = = N i= N i= N i= [ u n i ak u n i ak h ] h (un i+ u n i ) + u n ak h (un u n N ) ( N ) N u n i+ i= u n i ak h (un N u n ) N u n i. i= i= u n i Figure 4. (first row, left) shows the numerical solution using the above central scheme for a =, h =., N =, k =., with initial data u (x) = sin(πx), x, (4.5) at time t = (broken line) and t =.5 (solid line). The result is an approximation of sin(π(x.5)). We see that the solution is oscillating. This can be improved by using an arithemetic average in the approximation of the time derivative and leads to the following scheme. Lax-Friedrichs scheme. In this scheme the time derivative is approximated by (u(x, t + k) ) k (u(x + h, t) + u(x h, t)), i.e. i = (un i+ + u n i ) ak h (un i+ u n i ), i =,..., N. (4.6) The initial and boundary data are approximated by the same idea as for the central scheme. It is not difficult to check that the Lax-Friedrichs scheme is mass-preserving. Its numerical solution in Figure 4. (first row, right) with the same parameters as above shows that there are no oscillations, but now the solution is too smeared out. happened? error What We use the exact solution in the scheme (4.6) to compute the so-called truncation T E(x, t) = k (u(x, t + k) ) (u(x + h, t) + u(x h, t)) + a (u(x + h, t) u(x h, t)). h 8
19 Central Lax Friedrichs Downwind Upwind Upwind Upwind Figure 4.: Various numerical schemes for (4.) with a =, h =., k =. and smooth initial data (4.5). The numerical solution for t = (broken line) and t =.5 (solid line) is shown. First row: central (left), Lax-Friedrichs (right); second row: downwind (left), upwind (right); last row: upwind but with h =. (left) and k =.5 (right). A Taylor expansion yields at (x, t) T E = [( u + u t k + ) k u ttk + O(k 3 ) ( u + uxx h + O(h 4 ) )] + a h (u xh + O(h 3 )) = u t + au x + ) h (u ( h tt k u xx + O(k 4 ) ) + O + O(h ). k k Since u tt = au xt = a u xx and assuming that k/h = const., we conclude T E(x, t) = k ( ( h ) ) a u xx (x, t) + O(h ) = O(k). (4.7) k 9
20 This shows first that the Lax-Friedrichs scheme is a first-order method (as the truncation error satisfies T E(x, t) Ck for all (x, t)) and second that we approximate up to an error of the form ku xx. Now, spatial second derivatives are modeling diffusion phenomena, and we expect the discrete solutions to be smeared out justified by the numerical experiments. The numerical solution generated by the Lax-Friedrichs scheme can serve as an approximation of the advection-diffusion equation u t + au x = k ( ( h ) ) a u xx, x, t >. k For h and k, the solutions of this modified equation converge (at least formally) to the solution of u t + au x =. This is related to the vanishing viscosity limit discussed in Section 3. The Lax-Friedrichs approximation becomes better and better for smaller k >. Moreover, the artificial diffusion (also called artificial viscosity) avoids oscillations. Downwind scheme. The Lax-Friedrichs scheme gives accurate approximations only if k (or, equivalently, h if k/h = const.) is sufficiently small. One-sided finite difference approximations avoid too much artificial diffusion. Therefore, we choose the approximation i = u n i ak h (un i+ u n i ), i =,..., N, with initial and periodic boundary conditions analogously as above. Also this scheme is mass-preserving. The numerical solution shows that the scheme is unstable (Figure 4. (second row, left); same parameters as above). Why? The solution describes a wave from the left to the right. The spatial derivative at x i, however, uses the information at x i+ where the wave will go in the next time step. This does not make sense. In fact, the above scheme is useless. It would be more reasonable to use the information at x i where the wave comes from. This is done in the following scheme. Upwind scheme. This scheme reads as follows i = u n i ak h (un i u n i ), i =,..., N. (4.8) Again, the scheme is mass-preserving. Its numerical solution with the same parameter values as in the previous schemes in Figure 4. (second row, right) shows the correct solution, no oscillations, but with artificial diffusion. In fact, the truncation error T E(x, t) = k (u(x, t + k) u(x, t)) + a (u(x, t) u(x h, t)) h can be written as T E(x, t) = ak ( a h ) u xx + O(h ) + O(k ). (4.9) k
21 The corresponding equation with artificial diffusion reads here: u t + au x = ah ( a h ) u xx. (4.) k For the values of a, k and h used here, the diffusion term is.45u xx, whereas the corresponding term in the Lax-Friedrich scheme is.495u xx. Thus, we expect upwind to be less diffusive than Lax-Friedrichs, as confirmed by the numerical experiments. Choosing h > smaller, the diffusion becomes less apparent and for h =. the numerical solution is close to the exact solution (Figure 4. (last row, left)). For small k >, we need a lot of time steps to compute the solution at a certain fixed time. Can we accelerate the computations by choosing k larger? Figure 4. (last row, right) shows the result for h =. and k =.5 (a and u are as above). Thus, the answer is no. Why? The modified equation (4.) is well-posed only if the diffusion coefficient is non-negative: ah ( a h ) ak k h. (4.) In the converse case we do not have diffusion but concentration. In the above situation we have ak/h =.5 >, and we expect the solution to break down after some time. The condition (4.) is known as the Courant-Friedrich-Levy (CFL) condition. Strictly speaking, the CFL condition is not defined by (4.) but equivalent to it, and we will identify the CFL condition and (4.). This condition imposes a severe restriction on the choice of the time step; the more accurate the solution is computed (we mean: h > small), the smaller k > has to be chosen and the longer the computations take. In Figure 4. we illustrate the described behavior using discontinuous initial data { : x < / u (x) = (4.) : / x and the parameters a =, t =.5, k =.. We restrict the computational domain to [, ]. We do not choose periodic boundary conditions here, but inflow- and outflow-type conditions. In the traffic flow interpretation, the traffic is heavy in [, /) and light in [/, ]. So, at x =, cars are entering the domain and are leaving at x =. Hence, in the discrete formulation, = u n and N = un N, n. (4.3) We observe again that the central scheme is oscillatory with damped oscillations for smaller k (Figure 4. (first row, left); h =.); the Lax-Friedrichs scheme is quite diffusive but not oscillatory (Figure 4. (first row, right)); and the upwind scheme is less diffusive than the Lax-Friedrichs scheme (Figure 4. (last row, left)). Choosing the mesh
22 .5 Central.5 Lax Friedrichs.5.5 k =. k =.5 h =. h = Upwind Lax Friedrichs/Upwind h =. h =. Lax Friedrichs Upwind Figure 4.: Various numerical schemes for (4.) with a = and discontinous data (4.). First row, left: central; first row, right: Lax-Friedrichs; last row, left: upwind; last row, right: Lax-Friedrichs (broken line) and upwind (solid line) for h =., k =.. size h = k =., both schemes produce a solution which is very close to the exact solution (Figure 4. (last row, right)). This is clear since for this choice, a h/k = holds and thus, the artificial diffusion in (4.7) and (4.9) vanishes. All these schemes are able to compute the correct shock speed. Summarizing the above results we observe that the central scheme tends to produce oscillations whereas the other diffusive schemes have usually too much artificial viscosity except for special mesh sizes. The traffic flow models are nonlinear, so we want to study what can happen when we discretize nonlinear equations. We choose the scaled Lighthill-Whitham-ichards model (which is the same as the inviscid Burgers equation) u t + uu x =, x, t >, (4.4) u(x, ) = u (x), x, (4.5)
23 with initial data (4.). A natural generalization of the linear upwind scheme could be i = u n i k h un i (u n i u n i ), t Z, n. Since u i = if i N/ and u i u i = for i < N/, u i = u i for all i. This happens in every time step and so u n i = u i for all i. As the grid is refined, the numerical solution converges to the function u(x, t) = u (x). This is not a weak solution of (4.4) (4.5)! One might think that the upwind scheme i = u n i k h un i (u n i u n i ), i Z, n, (4.6) gives better results since here we use the information from the left where the wave is coming from. The numerical solution of this scheme (with boundary conditions (4.3) and k =., t =.5) is depicted in Figure 4.3 (left). For smaller mesh size, the numerical solution converges nicely to a function of type u (x at). The exact solution of (4.4) (4.5), (4.) is u(x, t) = { : x < st : x > st with s = (u l + u r ) =. At t =, the discontinuity should be at x = 3. Thus, the 4 numerical solution propagates with the wrong speed!.5 Upwind.5 Lax Friedrichs.5.5 exact h =. h =. exact h =. h = Figure 4.3: Exact and numerical solutions for the inviscid Burger equation using the upwind scheme (4.6) (left) and the Lax-Friedrichs scheme (4.7) (right). A better behavior is given by the Lax-Friedrichs scheme i = (un i+ + u n i ) k ( (u n 4h i+ ) (u n i ) ), i Z, n. (4.7) 3
24 From Figure 4.3 (right) we see that the scheme is non-oscillatory and the shock speed is correct. What is the reason for the different shock speeds? The upwind scheme (4.6) is a discretization of the quasilinear equation u t + uu x =, whereas the scheme (4.7) is an approximation of the equation in conservation form ( ) u u t + =. For smooth solutions both equations are the same, but we know from Section (see, for instance, Example.8) that this may not be true for weak solutions. In the following we consider only numerical methods being in conservation form, meaning that the scheme is of the form i x = u n i k h [F (un i p,..., u n i+q) F (u n i p,..., u n i +q)] for some function F of p + q + arguments. We call F the numerical flux function. The simplest case is for p = and q =, where i = u n i k h [F (un i, u n i+) F (u n i, u n i )]. (4.8) This expression can be interpreted as a cell average. More precisely, we know that the weak solution of satisfies the integral form h xi+/ u t + f(u) x =, x, t >, (4.9) xi+/ u(x, t n+ ) dx = u(x, t n ) dx (4.) x i / h x i / k [ tn+ f(u(x i+/, t)) dt tn+ ] f(u(x i /, t)) dt, h k k t n where x i±/ = (i ± /)h denotes the cell middle points (see Figure 4.4). Interpretating u n i as approximations of the cell averages, u n i h xi+/ x i / u(x, t n ) dx, and F (u n i, u n i+) as approximations of the average flux through x i+/ over the time interval (t n, t n+ ), F (u n i, u n i+) k tn+ we obtain the approximation (4.8) from (4.). t n f(u(x i+/, t)) dt, We give now some examples of numerical methods in conservation form: t n 4
25 x. x x x ẋ x x i-/ i+/ i- i i+ N Figure 4.4: The spatial mesh. Example 4. () Lax-Friedrichs scheme: i = (un i+ + u n i ) k h (f(un i+) f(u n i )). This method can be written in the conservation form (4.8) by taking F (u n i, u n i+) = h k (u i u i+ ) + (f(u i) + f(u i+ )). It can be proved that this scheme is a first-order method (i.e., the truncation error T E satisfies T E = O(k) as k ). () Lax-Wendroff scheme: i = u n i k h (f(un i+) f(u n i )) + k h [f (u n i+/)(f(u n i+) f(u n i )) f (u n i /)(f(u n i ) f(u n i ))], where u n i±/ = (un i + u n i±/ ). This is a second-order method (i.e. T E = O(k ) as k ). It has the disadvantage that it requires evaluating f in two points at each time step which is expensive when we deal with systems of conservation laws (and f is the Jacobi matrix). This is avoided in the following schemes. (3) ichtmyer two-step Lax-Wendroff scheme: u i = (un i + u n i+) k h [f(un i+) f(u n i )], i = u n i k h [f(u i ) f(u i )]. This is a second-order method. Notice that the derivation of the local truncation error T E(x, t) needs smooth solutions, so all these second-order methods are of second order for smooth solutions. around shocks. The methods are usually only of first order 5
26 (4) Mac-Cornack scheme: u i = u n i k h [f(un i+) f(u n i )], i = (un i + u i ) k h [f(u i ) f(u i )]. This scheme uses first forward differencing and then backward differencing. It is also a second-order method. Notice that the schemes () (4) reduce to the Lax-Wendroff method if f is linear. All the described schemes in Example 4. yield numerical solutions which converge to a weak solution of (4.9) if the mesh discretization parameters h and k tend to zero. In order to formulate this result more generally we need the following definitions. From now on, we assume that f : n n is a vector-valued function, i.e., we consider systems of conservation laws (4.9). Definition 4. () A difference scheme of the form i = u n i k h [F (un i p,..., u n i+q) F (u n i p,..., u n i +q)] for some function F : ( n ) p+q+ n is called conservative. () A conservative scheme is called consistent if F is locally Lipschitz continuous and F (u,..., u) = f(u) u. (3) The total variation T V (v) of a function v : n is defined by T V (v) = sup N v(ξ i ) v(ξ i ), i= where the supremum is taken over all subdivisions = ξ < ξ <... < ξ n = of the real line. Notice that for the total variation to be finite v must approach constant values as x ±. For differentiable v the definition reduces to T V (v) = v (x) dx. Theorem 4.3 (Lax-Wendroff) Let u j (x, t) be a numerical solution computed with a consistent and conservative method on a mesh with mesh size h j and time step k j, with h j, k j as j. (The function u j can be, for instance, the constant extension of u n i in the cells.) Assume that there exists a function u(x, t) such that 6
27 () for all a, b, T >, T b a u j (x, t) u(x, t) dxdt as j, () for all T > there is a number K > such that T V (u j (, t)) K t T, j N. Then u(x, t) is a weak solution of (4.9). For the proof we refer to [5, Sec. ]. Theorem 4.3 does not guarantee that weak solutions obtained in this way satisfy the entropy condition. This is true under the following additional condition. Theorem 4.4 Let (η, ψ) be an entropy-entropy flux pair such that η (u) > and η (u) div f(u) = div ψ(u) for all u n. Furthermore, let Ψ be a numerical flux function consistent with ψ in the sense of Definition 4. (). Let the assumptions of Theorem 4.3 hold and, additionally, η( i ) η(u n i ) k h [Ψ(un i p,..., u n i+q) Ψ(u n i p,..., u n i +q)] for all i, n. Then u(x, t) (obtained in Theorem 4.3) satisfies the entropy inequality (.3). Again, we refer to [5, Sec. ] for a proof. 5 The Gudonov method In Section 4 we have seen that a nonlinear conservation law (or a system of conservation laws) can be numerically approximated by the Lax-Friedrichs scheme. For the linear advection equation (4.) we have also seen that the Lax-Friedrichs scheme is generally more dissipative than the upwind method and gives less accurate solutions. In the scalar case, a natural generalization of the upwind scheme (4.8) is i = u n i k h [F (un i, u n i+) F (u n i, u n i )] (5.) with F (v, w) = { f(v) : (f(v) f(w))/(v w) f(w) : (f(v) f(w))/(v w) <. 7
28 For linear equations, F (v, w) = f(v) and (5.) reduces to (4.8). However, the above scheme may not give the correct approximation. For instance, let f(u) = u / and choose the initial data u i = { : i + : i >. Since F (u i, u i+) and F (u i, u i ) are either equal to f( ) or f() and since f( ) = f(), we obtain u i = u i for all i and hence u n i = u i for all i. This is not the correct discrete solution. In this section we derive a conservative and consistent generalization of the upwind scheme which avoids the above problem, the so-called Gudonov scheme. Let f be a convex C function. The idea of the method is to approximate the solution u(x, t n ) of the conservation law (scalar or system) u t + f(u) x =, x, t >, by a piecewise constant function ũ n (x, t n ) and to determine the approximate solution ũ n (x, t) by solving the iemann problem in the interval t [t n, t n+ ]. After obtaining this solution, we define the approximate solution i at time t n+ by averaging this exact solution at time t n+ : i = h xi+/ x i / ũ n (x, t n+ ) dx, (5.) where x i±/ = (i±/)h. These values are then used to define the new piecewise constant data ũ n+ (x, t n+ ) by and the process repeats. ũ n+ (x, t n+ ) = i if x i / x < x i+/, In practice, this algorithm is considerably simplified since the above integral can be computed explicitely. Since ũ n is assumed to be the exact weak solution, it satisfies the integral formulation (.3) divided by h: h xi+/ x i / xi+/ = h k h ũ n (x, t n+ ) dx x i / tn+ [ k ũ n (x, t n ) dx t n f(ũ n (x i+/, t)) dt k tn+ t n ] f(ũ n (x i /, t)) dt. From (5.) follows i = u n i k h [ F (u n i, u n i+) F (u n i, u n i ) ], (5.3) 8
29 where the numerical flux function F is given by F (u n i, u n i+) = k tn+ t n f(ũ n (x i+/, t)) dt. Thus, the Gudonov scheme is conservative (see Definition 4.). The function ũ n is constant on the line x = x i+/, t n t t n+ (see Figure 5.). We denote this value by u (u n i, u n i+). t n+ ψ (t) ~ u n (x,t) n u n i- n u u i+ i w(t) t n x i-/ x i+/ x Figure 5.: Illustration for the Gudonov scheme. There is a shock through x i / and a rarefaction wave starting at x i+/. Then the flux reduces to F (u n i, u n i+) = f(u (u n i, u n i+)) and the Gudonov scheme becomes i = u n i k h [f(u (u n i, u n i+)) f(u (u n i, u n i ))]. The scheme is consistent (see Definition 4.) since F (u n i, u n i ) = f(u n i ) and f is assumed to be smooth. There are two questions remaining: For large t t n, the solution may not remain constant at x i+/ because of the effects of waves arising from neighboring iemann problems. How large can k = t n+ t n be chosen? How can we determine u from u n i, u n i+? We answer the first question. Assume the situation of Figure 5., i.e., let ψ(t) be the shock line through x i / and let w(t) be the left end of the rarefaction wave starting at x i+/. The time t n+ is determined by the requirement ψ(t n+ ) w(t n+ ). Since ψ(t) = x i / + s(t t n ) with s = f(un i+) f(u n i ) u n i+ un i and w(t) = x i+/ + f (u n i )(t t n ), 9
30 this yields h = x i+/ x x / (s f (u n i ))(t n+ t n ) = (s f (u n i ))k. (5.4) As f is assumed to be convex, s lies between f (u n i ) and f (u n i+), i.e. Thus, if we infer s max{ f (u n i ), f (u n i+) }. ν := sup i,n f (u n i )k h s f (u n i+) s + f (u n i ) sup f (u n i ) h i,n k, (5.5) and (5.4) is satisfied. The condition (5.5) ensures that the shock and the rarefaction wave do not interact in the mesh cell [x i /, x i+/ ] [t n, t n+ ). We obtain the same condition (5.5) if there is a rarefaction wave at x i / and a shock at x i+/ or if there are two shocks at x i±/ since the wave speeds are always bounded by sup f (u n i ). In fact, we can allow the waves to interact during the time step, provided the interaction is entirely contained within a mesh cell. This leads to the condition sup f (u n i ) h/k or ν = sup i,n f (u n i ) k. (5.6) h This condition can be interpreted as a generalization of the CFL condition (4.) for linear conservation laws. In the case of systems of conservation laws we replace (5.6) by sup λ p (u n i ) k, h i,n,p where λ p (u n i ) are the eigenvalues of the Jacobian f (u n i ). This answers the first question. Concerning the second question, we need to solve the iemann problem to determine u (u n i, u n i+). For systems of conservation laws, this may be expensive. Notice however, that most of the structure of the iemann solver is not used in the Gudonov scheme and therefore, approximate iemann solvers have been devised to improve the efficiency of the Gudonov method. We do not describe these approximate solvers but refer to [5, Sec. 4.]. For scalar conservation laws, we can determine u easily from the sign of f (u n i ) and f (u n i+). Indeed, we consider the following four cases: () f (u n i ) and f (u n i+) : In this case there is a rarefaction wave starting at x i+/, and from Figure 5. (a) we see that u = u n i. () f (u n i ) < and f (u n i+) < : Again, there is a rarefaction wave starting at x i+/ but now u = u n i+ (see Figure 5. (b)). 3
31 (3) f (u n i ) and f (u n i+) < : There is a shock through x i+/ and (Figure 5. (c)) { u u n i : s = u n i+ : s <. (4) f (u n i ) < and f (u n i+) : There is a rarefaction wave starting at x i+/ and u is the unique solution of f (u ) = (since f is convex). As the wave speed passes through zero within the wave it is called transonic rarefaction wave (Figure 5. (d)). The value u is called sonic point. (a) (b) (c) (d) s> n u i n u i u n i+ n u n i+ u i u n i+ n u i u n i+ x i+/ x i+/ x i+/ x i+/ Figure 5.: Computation of u. The resulting flux function can be written in the simplified form min f(u) : u n F (u n i, u n u i+) = n i u n i+ i u un i+ max f(u) : u n u n i > u n i+. i+ u un i Surprisingly, this expression also holds for general conservation laws, even non-convex ones, and gives the correct Gudonov flux corresponding to the weak solution satisfying the entropy condition (.9) of Oleinik. Example 5. For f(u) = au (a > ) the flux function becomes F (u n i, u n i+) = au n i, and hence, from (5.3), the Gudonov scheme is equal to the upwind scheme for linear equations (see Section 4), i = u n i ak h (un i i ). In this sense, the Gudonov method is a generalization of the upwind scheme to nonlinear equations. We summarize the Gudonov scheme for scalar conservation laws: 3
32 Initialize u i = h xi+/ x i / u (x) dx. For all i, n do: if f (u n i ) and f (u n i+) then u i = u n i ; if f (u n i ) < and f (u n i+) < then u i = u n i+; if f (u n i ) and f (u n i+) < then u i = u n i (if s ) or u i = u n i+ (if s < ); if f (u n i ) < and f (u n i+) then u i is the unique solution of f (u i ) =. Set i = u n i k h [f(u i ) f(u i )]. eferences [] A. Aw and M. ascle: esurrection of second order models of traffic flow. SIAM J. Appl. Math. 6 (), [] A. Aw, A. Klar, T. Materne, and M. ascle: Derivation of continuum traffic flow models from microscopic follow-the-leader models. Submitted for publication,. [3] C. Daganzo: equiem for second-order fluid approximation to traffic flow. Transport. es. B 9 (995), [4] M. Hanke-Bourgeois: Grundlagen der Numerischen Mathematik und des Wissenschaftlichen echnens. Teubner, Stuttgart,. [5]. LeVeque: Numerical Methods for Conservation Laws. Birkhuser, Basel, 99. [6] G. Warnecke: Analytische Methoden in der Theorie der Erhaltungsgleichungen. Teubner, Stuttgart
Notes: 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 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 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 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 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 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 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 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 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 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 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 informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
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 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 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 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 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 informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
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 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 informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
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 informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
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 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 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 informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
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 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 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 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 informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
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 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 informationThe Aw-Rascle traffic flow model with phase transitions
The traffic flow model with phase transitions Paola Goatin Laboratoire d Analyse Non linéaire Appliquée et Modélisation I.S.I.T.V., Université du Sud Toulon - Var B.P. 56, 83162 La Valette du Var Cedex,
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 informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
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 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 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 informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
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 informationHyperbolic Conservation Laws Past and Future
Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
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 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 numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationSolving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation
Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation W.L. Jin and H.M. Zhang August 3 Abstract: In this paper we study the Payne-Whitham (PW) model as
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationDefinition and Construction of Entropy Satisfying Multiresolution Analysis (MRA)
Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 6 Definition and Construction of Entropy Satisfying Multiresolution Analysis (MRA Ju Y. Yi Utah State University
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationhal , version 1-18 Nov 2010
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
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 informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationHyperbolic Partial Differential Equations
Chapter 1 Hyperbolic Partial Differential Equations We begin our study of finite difference methods for partial differential equations by considering the important class of partial differential equations
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 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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
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 informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationON A DIFFUSIVELY CORRECTED KINEMATIC-WAVE TRAFFIC MODEL WITH CHANGING ROAD SURFACE CONDITIONS
Dept. of Math. University of Oslo Pure Mathematics ISBN 8 55 85 No. 9 ISSN 86 49 June ON A DIFFUSIVELY CORRECTED KINEMATIC-WAVE TRAFFIC MODEL WITH CHANGING ROAD SURFACE CONDITIONS R. BÜRGERA AND K.H. KARLSEN
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationNumerical methods for conservation laws with a stochastically driven flux
Numerical methods for conservation laws with a stochastically driven flux Håkon Hoel, Kenneth Karlsen, Nils Henrik Risebro, Erlend Briseid Storrøsten Department of Mathematics, University of Oslo, Norway
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
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 informationAffordable, entropy-consistent, Euler flux functions
Affordable, entropy-consistent, Euler flux functions (with application to the carbuncle phenomenon) Phil Roe Aerospace Engineering University 0f Michigan Ann Arbor Presented at HYP 2006 1 Entropy pairs
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
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 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 informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
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 informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
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 informationNotes on Computational Mathematics
Notes on Computational Mathematics James Brunner August 19, 2013 This document is a summary of the notes from UW-Madison s Math 714, fall 2012 (taught by Prof. Shi Jin) and Math 715, spring 2013 (taught
More information