The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

Size: px
Start display at page:

Download "The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria"

Transcription

1 ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Delta and Singular Delta Locus for One Dimensional Systems of Conservation Laws M. Nedeljkov Vienna, Preprint ESI 837 (000) February 9, 000 Supported by Federal Ministry of Science and Transport, Austria Available via

2 DELTA AND SINGULAR DELTA LOCUS FOR ONE DIMENSIONAL SYSTEMS OF CONSERVATION LAWS Marko Nedeljkov Abstract. A condition for existence of singular and delta shock waves for systems of conservation laws is given in the paper. The systems considered here have uxes which are linear in one of the dependent variables. The condition obtained here is analogous to the one for the standard Hugoniot locus. Three dierent solution concept are used in the paper: associated solution in Colombeau sense, limits of nets of smooth functions together with Rankin-Hugoniot conditions and a kind of a measure valued solutions. 1. Introduction The aim of this paper is to give some criterion for which certain classes of function nets, called delta and singular shock waves, satisfy the following system of conservation laws (1) () u t + (f 1 (u)v + f (u)) x = 0 v t + (g 1 (u)v + g (u)) x = 0; in an approximated sense. Here, f i, g i, i = 1; are smooth functions, polynomially bounded together with all their derivatives, and u; v : R R +! R satisfy (3) u(x; 0) = u 0 ; x < 0; u(x; 0) = u 1 ; x > 0 v(x; 0) = v 0 ; x < 0; v(x; 0) = v 1 ; x > 0: A pair of nets (u " ; v " ) D 0, " (0; 1), is a solution in an approximated sense if hu "t + (f 1 (u " )v " + f (u " )) x ; i! 0 hv "t + (g 1 (u " )v " + g (u " )) x ; i! 0 as "! 0, for every D(R R + ). In order to have a well dened composition of functions we use families of smooth function which can be viewed as representatives of Colombeau functions ([1], [9]) or families of piecewise constant functions which are called box approximations. For the rst case we will use the space of generalized functions, G g, and the concept of a approximate solution dened in [10]. In the second case we use the standard Rankin-Hugoniot conditions for systems of conservation laws like in [6]. Roughly speaking, delta and singular shock wave are pairs of nets (u " ; v " ) which converge to sums of step and delta functions in D 0. They will be precisely dened bellow. Let us remark that singular shock waves have similarities to innite narrow 1 Typeset by AMS-TEX

3 M. NEDELJKOV solitons which are introduced by Maslov and Omel'anov [8]. Further results in this direction for systems of conservation laws can be found in the paper [3] from Danilov, Maslov and Shelkovich. For a given point (u 0 ; v 0 ) R, a delta (singular delta) locus is a subset of R consists of points (u 1 ; v 1 ) for which there exists a delta (singular) shock wave joining (u 0 ; v 0 ) and (u 1 ; v 1 ). This denition is an analogue of the denition of the classical Hugoniot locus in the case of piecewise constant solutions (shock waves) for Riemann problem. Admissibility condition for such shocks could be the overcompressiveness like in [6], [11], [4] and [5], i.e. (u 0 ; v 0 ) > 1 (u 0 ; v 0 ) c (u 1 ; v 1 ) > 1 (u 1 ; v 1 ); where c is a speed of the delta or singular shock wave, 1 ; are the eigenvectors for the system, (u 0 ; v 0 ) and (u 1 ; v 1 ) denotes left and right hand side initial value, respectively. Let us remark that the hyperbolicity condition is not used in the construction of approximate solutions. Like in the classical theory of Hugoniot locus, the (singular) delta locus can be used for construction of solutions containing also a rarefaction wave on the left or on the right hand side of a singular or delta shock. Our investigations are motivated by Keytz and Kranzer ([6]) who found singular shock wave solutions to the system (4) u t + (u? v) x = 0 1 v t + 3 u3? u = 0; r t in the form u " (x; t) = G " (x? ct) + a a t x x? ct " ", v " (x; t) = H " (x? ct) + x? ct ", where G " and H " converge to appropriate step functions dened " by the Riemann initial data (3), " := "?1 (=") converges to the delta distribution and i ", i = 1; 3 converge to zero in D 0 as "! 0. This system is included in our general system (1), (). Our paper shows that their approach is suciently general for solving (1), (). An example of a system with a delta shock wave solution is (5) u t + (u ) x = 0 v t + (uv) x = 0 examined in [11]. The solution is of the form u " (x; t) = G " (x? ct), v " (x; t) = H " (x? ct) + st " (x? ct), where G " and H " are the same as before and " is a delta net. By adding a viscosity term Tan, Zhang and Zheng ([11]) proved that the limit of vanishing viscosity solutions is a solution of this form. Let us notice that there is no additional "singular" term in u ". Here, the delta shock wave will be used to denote that there are no other terms in a solution except approximations of step functions and the delta distribution. Systems which are to follow have the similar form of solutions as (5), that is they do not contain term added to G " and a delta net is a part of v ".

4 DELTA AND SINGULAR DELTA LOCUS 3 Oberguggenberger ([9]) proved that for (6) u t + (u =) x = 0 v t + (uv) x = 0 the viscosity limit is a delta shock wave. By using Le Flock and Vol'pert denition of the product, Hayes and Le Flock ([5]) found a delta shock wave solution to (7) u t + (u ) x = 0 v t + ((u? 1)v) x = 0: They also proved that the vanishing viscosity limit exists and it is equal to the solution of (7). Finally, Ercole ([4]) proved that the system (8) u t + f(u) x = 0 v t + (g(u)v) x = 0; with some minor assumptions on f and g, has a delta shock solution as a limit of vanishing viscosity solutions. The dierences between system (4) and systems (5)-(8) are obvious. Only system (4) has second dependent variable in the rst ux function and it can not be solved with delta shock waves. In this paper we will not use the viscosity approximation. The existence and convergence of vanishing viscosity solutions is not a trivial task even for some specic choice of a function in the ux. An example of this one can nd in Dafermos- DiPerna paper []. All the systems considered here have a ux which is a linear function in one variable and this will be a general form of a ux for which we nd a delta locus. The interpretation of a limit of " where " is a delta net brings diculties. Colombeau proved in [1] that " denes an element in the space of the generalized functions which is not associated with any classical distribution. In all the cited papers a delta locus is a subset of R with the non-zero Lebesgue measure. In general, this set is just a curve, but for systems (5)-(8) our procedure gives the same locus as it was already obtained. Singular delta locus is an set with non-zero Lebesgue measure in general. In the sequel we shall consider only solutions which equals to initial data for every " (also after regularization when we use smooth nets) at the initial time. This means that the \value" of the Heaviside function in zero does not have any inuence, for example. The paper is organised as follows. In Section the algebra of generalized functions, delta and singular shock wave solution are precisely dened. In Section 3 we give a form of approximate solutions which will be used in nding delta locus's for smooth functions f i, g i, i = 1;. Roughly speaking, the main property of approximations is a distribution of a \mass" of singular part on the left and right-hand side from zero. Properly placed right and left hand sided masses determines the delta locus. For the use of a delta locus in constructions of admissible solutions one has to remove the points for which the overcompressiveness do not hold.

5 4 M. NEDELJKOV A box approximation solution is used in Section 4 and it gives the same delta locus as in the case of Colombeau generalized functions. Theorem 1 applied to systems (5)-(8) gives the results previously obtained in the quoted papers. System (4) can not be solved by this procedure. A concept and construction of a solution in a certain measure space is given in Section 5. Approximated solutions constructed in Section 3 converge to this solution. In Section 6 systems with polynomial uxes are examined. There exists a broad class of polynomials for which we can nd singular delta locus of the non-zero Lebesgue measure modifying u " by adding a net which converge to zero but some of its powers do not. In fact we use an idea of [6] adopted for arbitrary polynomials. In contrast to the singular shock wave solution obtained in [6], we have to assume that singular parts of approximations u " and v " have disjoint supports. This can be omitted if f 1 and g 1 are constants, or linear functions in u. There are many open problems concerning the system (1), (). For example one can try to obtain a limit of viscosity self-similar solutions to the system like it was done in Dafermos and DiPerna's paper []. Also, one can try to describe singular shock wave solution formally obtained as a net of approximate solutions by using a weighted measure space as it was done by Keytz and Kranzer in [6].. Definitions We shall briey recall some of the denitions of Colombeau algebra given in [10]. Let R + = R (0; 1), R + = R [0; 1) and Cb 1 () be the algebra of smooth functions on bounded together to all its derivatives. Let C 1 b (R + ) be a set of all (R (0; T )) for every T > 0. Let functions u C 1 (R +) such that ujr(0;t ) Cb 1 us remark that every element of Cb 1(R +) has a smooth extension up to the line ft = 0g, i.e. Cb 1(R +) = Cb 1(R +). This is also true for C 1 b (R +). E M;g (R +) is the set of all maps G : (0; 1) R +! R, ("; x; t) 7! G " (x; t), smooth for every " (0; 1), which satisfy: For every (; ) N 0 and T > 0, there exists N N such that sup j@ t G " (x; t)j = O("?N ); as "! 0: (x;t)r(0;t ) N g (R +) is the set of all G " E M;g (R +) which satisfy: For every (; ) N 0, a R and T > 0 sup j@ t G " (x; t)j = O(" a ); as "! 0: (x;t)r(0;t ) The dierential algebra G g (R +) of generalized functions is the factor algebra G g (R + ) = E M;g(R + )=N g(r + ). By G will be denoted a class in G g which has a representative G " E M;g. Since C 1 b (R +) = C 1 b (R +) one can dene a restriction of a generalized function to ft = 0g as follows. For given G G g (R + ), its restriction Gj t=0 G g (R) is the class which contain a function G " (x; 0). Similarly, G(x? ct) G g (R) is dened with G " (x? ct) E M;g (R). Let us note that if f is a smooth function polynomially bounded together with all its derivatives, then the composition f(g), G G g is well-dened.

6 DELTA AND SINGULAR DELTA LOCUS 5 Let u D 0 L1(R). Fix an element D(R), (x) 0, x R, R (x)dx = 1, supp [?1; 1]. Let " (x) = "?1 (x="), x R. Then : u! class of u " denes an embedding of D 0 L1(R) into G g (R). Obviously, commutes with derivative. We shall call a generalized function G a step function if it equals (u), where u is a step function in the usual sense. The delta generalized function is given by a representative ". In the sequel we shall omit the subscript. An element G G g () is associated with u D 0 () if for some (and hence every) representative G " of G, G "! u in D 0 () as "! 0. It will be denoted by G u. Two generalized functions G and H are associated, G H if G? H 0. A generalized function G is of the bounded type if sup jg " (x; t)j = O(1) as "! 0; (x;t)r(0;t) for every T > 0. The initial data for the Riemann problem are regularized in the following way. Let u0 ; x < 0 v0 ; x < 0 1 (x) = and (x) = u 1 ; x > 0; v 1 ; x > 0: The initial data are now (G; H) where G " = 1 " and H " = ". This means that a solution to system (1), () satises u " (x; 0) = G " (x); G " (x) = u 0 ; x <?"; G " (x) = u 1 ; x > "; v " (x; 0) = H " (x); H " (x) = v 0 ; x <?"; H " (x) = v 1 ; x > "; where G " and H " are of the bounded type, i.e. they are uniformly bounded with respect of ". Generalized functions of this type are called the generalized step functions. One will see that the results obtained here will be valid for all representatives of step functions which satisfy these conditions. Let us remark that if u " is a function satisfying the above conditions, this will be true for every its power u m, m N, also. Denition 1. (Generalized functions) (U; V ) G(R + ) is a singular shock wave solution to (1), () and (3) if a) (9) (10) (11) U t + (f 1 (U)V + f (U)) x 0 V t + (g 1 (U)V + g (U)) x 0: Uj t=0 = G; V j t=0 = H: b) U(x; t) 1 (x? ct), where c is a constant (speed of the shock) and V (x; t) (x? ct) + s(t)(x? ct), where c is the same constant, s C 1 ([0; 1)), s(0) = 0 and is Dirac delta function. If in addition U is of the bounded type, the solution is called delta shock wave.

7 6 M. NEDELJKOV Denition. (Box approximations) A pair of nets (u " ; v " ) consisting of piecewise smooth functions is called an approximated singular shock wave solution to (1), () and (3) if a) limhu " ; t i + hf 1 (u " )v " + f (u " ); x i = 0 "!0 lim hv "; t i + hg 1 (u " )v " + g (u " ); x i = 0: "!0 b) (u " ; v " ) satisfy (3) for " small enough. c) u " converge in the distributional sense to 1 (x?ct), where c is a constant (speed of the shock) and 1 a step function. v " converge in the distributional sense to (x? ct) + s(t)(x? ct), where c is the same constant, s C 1 ([0; 1)), s(0) = 0, ) is a step function and is Dirac delta function. If in addition u " is bounded independently on ", the solution is called delta shock wave. 3. Generalized functions and delta locus Theorem 1. a) Let f 1 6 const. Then a delta shock wave solution to (9), (10) and (11) exists if u 0 6= u 1, f 1 (u 0 ) 6= f 1 (u 1 ) and (1) c = f 1(u 1 )v 1 + f (u 1 )? f 1 (u 0 )v 0? f (u 0 ) u 1? u 0 = g 1(u 0 )f 1 (u 1 )? g 1 (u 1 )f 1 (u 0 ) ; f 1 (u 1 )? f 1 (u 0 ) where c is the velocity of the delta shock. The set of all pairs of right hand states (u 1 ; v 1 ) such that (1) holds is called the delta locus for (1), () and (3). b) If f 1 (u 0 ) = f 1 (u 1 ) = 0 (specially, if f 1 0) and g 1 6 const, then the delta locus is the set of all points (u 1 ; v 1 ) such that g 1 (u 0 ) 6= g 1 (u 1 ). c) If f 1 0 and g 1 b R, then the delta locus is the set of all points (u 1 ; v 1 ) such that b(u 1? u 0 ) = f (u 1 )? f (u 0 ). Proof. Let C 1 0 (R), supp [?1; 1], R (x)dx = 1, 0 and let D " = D? " A?1 x + A+1 " A?1! + D+ " A?1 x? A+1 " A?1 where D? + D + = 1 Let s 1 (t); s (t) be smooth functions for t 0, such that s 1 (0) = s (0) = 0 and A > 0. Since in the general we do not have any information on the behaviour of f 1, f, g 1 and g we shall try to nd delta shock wave solution in the form (13) u " (x; t) = G " (x? ct); v " (x; t) = H " (x? ct) + s(t)d " (x? ct); where D " is dened as above and s is from Denition 1. Let us note that D " converge to delta function as "! 0. After a substitution of (13) into (1) we have (14) G "t + (f 1 (G " )H " + f (G " )) x + s(t)(f 1 (G " )D " ) x 0:! ;

8 DELTA AND SINGULAR DELTA LOCUS 7 A missing argument of a function means that it equals x? ct. [?X; X] [?T; T ], for some X; T > 0. Then ZZ Let supp (@ t G " (x? ct) x (f 1 (G " (x? ct))h " (x? ct) + f (G " (x? ct))) (x; t)dxdt ZZ = (cg " (x?ct)? f 1 (G " (x?ct))h " (x?ct)? f (G " (x?ct)))@ x (x; t)dxdt = = = Z T Z ct?"?t?x Z T Z ct+"?t ct?" Z T Z X?T ct+" =I 1 + I + I 3 : At rst let us notice that (cg " (x?ct)? f 1 (G " (x?ct))h " (x?ct)? f (G " (x?ct)))@ x (x; t)dxdt (cg " (x?ct)? f 1 (G " (x?ct))h " (x?ct)? f (G " (x?ct)))@ x (x; t)dxdt (cg " (x?ct)? f 1 (G " (x?ct))h " (x?ct)? f (G " (x?ct)))@ x (x; t)dxdt ji j Z T?T "CC dt! 0; as "! 0; since all of G ", H " and are bounded independently on ". Next, I 1 + I 3 = + Z T Z ct?"?t?x Z T Z ct?"?t?x (cu 0? f 1 (u 0 )v 0? f (u 0 ))@ x (x; t)dxdt (cu 1? f 1 (u 1 )v 0? f (u 1 ))@ x (x; t)dxdt = (cu 0? f 1 (u 0 )v 0? f (u 0 ))? (cu 1? f 1 (u 1 )v 1? f (u 1 )) Z T?T Z T?T (ct? "; t)dt (ct + "; t)dt!?(c[u]? [f 1 (u)v + f (u)])j x=ct ; as "! 0: This means that the sum of rst two members of (14) is associated to (x? ct) multiplied by a constant, i.e. G "t + (f 1 (G " )H " + f (G " )) x?(c[u]? [f 1 (u)v + f (u)]): In the sequel we shall omit this type of calculations. The function (f 1 (G " )D " ) x is a derivative of a product of approximations for step and delta function. This product is associated with the delta function multiplied by a constant, and the constant depend of a representative of the delta function. This constant as well as (c[u]? [f 1 (u)v + f (u)]) has to be zero. (c[u]? [f 1 (u)v + f (u)]) = 0 gives the Rankine-Hugoniot conditions for G " and H " of the rst equation, i.e. c = f 1(u 1 )v 1 + f (u 1 )? f 1 (u 0 )v 0? f (u 0 ) u 1? u 0 :

9 8 M. NEDELJKOV Now we shall nd a representative of the delta function such that f 1 (G " )D " 0. Let us note that supp D [?A;?1] [ [1; A], A is a real number greater than 1, R D(y)dy = 1 and D" (y) = "?n D(y="). The condition on supports ensures that the product of a step function and will not depend on the regularization of the initial data. It will depend only on values of a step function away of zero, R 0?1 D(y)dy R?1 and D(y)dy. Thus, we have the following system 0 Z?1?A f 1 (u 0 ) The solution is given by D? := Z?1?a D(y)dy = D(y)dy + Z?1?A Z A 1 D(y)dy = 1 D(y)dy + f 1 (u 1 ) Z A f 1 (u 1 ) f 1 (u 1 )? f 1 (u 0 ) ; D+ := 1 Z b 1 D(y)dy = 0: D(y)dy =?f 1 (u 0 ) f 1 (u 1 )? f 1 (u 0 ) ; if f 1 (u 1 ) 6= f 1 (u 0 ). In the case when f 1 const 6= 0, we shall nd singular shock solutions of system (1), () but for some special forms of f, g 1 and g. There are no delta shock solutions of this form (only classical shocks, i.e. s = 0). If f 1 0, then D? and D + can be choosen such that cd 0 " (g 1 (G " )D " ) x : The only neccesary condition is that g 1 (u 0 ) 6= g 1 (u 1 ). If g 1 is also a constant then c equals to this constant. The set of all u 1 for which this is true determines the delta locus. When we substitute u " and v " into (), H "t + s 0 (t)d "? cs(t)d "t + (g 1 (G " )H " + g (G " ) + g 1 (G " )s(t)d " ) x =H "t + (g 1 (G " )H " + g (G " )) x + s 0 (t)d "? cs(t)d "t + s(t)(g 1 (G " )D " ) x : Since G " and H " are representatives of classical step functions, it follows that H "t + (g 1 (G " )H " + g (G " )) x?; where is Rankine-Hugoniot decit ([6]) given by = c(v 1? v 0 )? (g 1 (u 1 )v 1 + g (u 1 )? g 1 (u 0 )v 0? g (u 0 )): By putting = 0 one can see that (u 1 ; v 1 ) has to belong to the standard Rankine- Hugoniot locus (then s = 0, since only the term s 0 (t)d " is associated with the delta function multiplied by a constant and all others are associated with the derivative of the delta function multiplied by some constant). In the case 6= 0 the condition s 0 (t) = should be satised, which implies that s(t) = t and This implies that cd 0 " (g 1 (G " )D " ) x : (15) c = g 0(u 1 )f 1 (u 1 )? g 1 (u 1 )f 1 (u 0 ) f 1 (u 1 )? f 1 (u 0 ) (since the representative of delta function D " (i.e. D? and D + ) is already determinated). That means that the delta locus is dened as the set of all (u 1 ; v 1 ) such that (15) holds. Remark. In the case when f 1 const 6= 0 the delta locus is the empty set.

10 DELTA AND SINGULAR DELTA LOCUS 9 4. Box approximations and delta locus Theorem 1 is also true for the box approximations. In this case D " (y) consists of two boxes, D " (y) = 8 >< >: D? " D + ; y [?"; 0] ; y [0; "] " 0; otherwise: After the substitution of u and v into (1), by using the standard Rankine- Hugoniot conditions for step functions, we have (?c[g] + [f 1 (G)H + f (G)]) 0 + s(t)([f 1 (G)D " ]?"?" +[f 1 (G)D " ] [f 1 (G)D " ] " " ) = 0; where h a ; i = (a) and [] denotes a jump in a point. Again the same condition for the speed holds, c = [f 1(G)H + f (G)] ; [G] and h[f 1 (G)D " ]?"?" ; i = f 1 (u 0 )D? "?1 (?") = f 1 (u 0 )D? "?1 ( (0)? " 0 (y " )) = f 1 (u 0 )D? "?1 (0)? f 1 (u 0 )D? 0 (y " ); where y " is some number in [?"; 0]. Similarly, ~y " [0; "]. Finally, Then h[f 1 (G)D " ] " " ; i =?f 1 (u 1 )D + "?1 (0)? f 1 (u 1 )D + 0 (~y " ); h[f 1 (G)D " ] 0 0 ; i = (?f 1 (u 0 )D? + f 1 (u 1 )D + )"?1 (0): h[f 1 (G)D " ]?"?" + [f 1 (G)D " ] [f 1 (G)D " ] " " ; i =? f 1 (u 0 )D? 0 (y " ) + f 1 (u 1 )D + 0 (~y " )! (?f 1 (u 0 )D?? f 1 (u 1 )D + ) 0 (0) =(?f 1 (u 0 )D?? f 1 (u 1 )D + )h 0 ; i; as "! 0. This gives the following two conditions D? + D + = 1 f 1 (u 0 )D? + f 1 (u 1 )D + = 0: These conditions are the same as in the previous case. In the similar way as above one can obtain the third condition which denes the delta locus g 1 (u 0 )D? + g 1 (u 1 )D + = c:

11 10 M. NEDELJKOV The analysis of possible cases is the same as for smooth aproximations. In both cases it seems that one splits delta function in two parts, left delta,?, and right one, +, such that its product with piecewise continuous and bounded functions is dened by (y)f(y) = lim f(y); where is? or + : y!0 In the case when system (1), () is hyperbolic, admissible part of the delta locus is the set of (u 1 ; v 1 ) such that (u 1 ; v 1 ) c 1 (u 0 ; v 0 ) where 1 < are the eigenvalues of (1), (). Remark. Eigenvalues of (1), () equals 1; (u; v) = f0 1(u)v + f(u) 0 + g 1 (u) p (f 0 1 (u)v + f0(u)? g 1(u)) + 4f 1 (u)(g 0 (u)v + 1 g0 (u)) If one uses the box approximation then one can see that the hyperbolicity of the system means 1 6= R, for u = u i, v = v i, i = 1; and for v equals some number depending on " which tends to innity as "! 0. The last will be true if f 0 1(u i ) 6= 0, i = 1;. If one uses smooth approximation, then the question of hyperbolicity is much more dicult and highly depends on functions in the uxes. Overcompressivity condition can be written in the following explicit form f1(u 0 0 )v 0 + f(u 0 0 ) + g 1 (u 0 ) p (f 0? 1 (u 0)v 0 + f 0 (u 0)? g 1 (u 0 )) + 4f 1 (u 0 )(g 0 1 (u 0)v 0 + g 0 (u 0)) g 1(u 0 )f 1 (u 1 )? g 1 (u 1 )f 1 (u 0 ) f 1 (u 1 )? f 1 (u 0 ) f0 1(u 1 )v 1 + f(u 0 1 ) + g 1 (u 1 ) p (f (u 1)v 1 + f 0 (u 1)? g 1 (u 1 )) + 4f 1 (u 1 )(g 0 1 (u 1)v 1 + g 0 (u 1)) 5. Measure valued solutions In order to explain more simply the concept of solution we use the self similar solutions to (1) and (), i.e. u = u(x=t), v = v(x=t). In that case the system is (16) (17)?yu 0 (y) + (f 1 (u(y))v(y) + f (u(y))) 0 = 0?yv 0 (y) + (g 1 (u(y))v(y) + g (u(y))) 0 = 0; (18) u(?1) = u 0 ; u(1) = u 1 ; v(?1) = v 0 ; v(1) = v 1 : Our aim is to nd a solution (u; v), where u is an element of the space of piecewise continuous functions having the right and left hand limits in all points in R, P L (R), while v = (v 0 : : : :; v m+1 ) belonges to the space y0;:::;ymm(r) = M((?1; y 0 ])

12 DELTA AND SINGULAR DELTA LOCUS 11 M([y 0 ; y 1 ]) M([y m ; 1)), where y i, i = 0; : : :; m are some real numbers. In this case the products f 1 (u)v and g 1 (u)v exist and belong to the space y0;:::;ymm(r). Denote by M F (R) the union of all y0;:::;ymm(r), where fy 1 ; : : :; y m g is a nite set of real numbers dened above. Let w = (w 0 ; : : :; w m+1 ) y0;:::;ym M(R), w i = w i + i?1 + (? y i?1 ) + i? (? y i ); i = 1; : : :; m; w 0 = w 0 + 0? (? y 0 ); w m+1 = w m+1 + m + (? y m ); where w i are the measures with the zero mass in the points fy 0 ; : : :; y m g. The mapping : w 7! (w) = W D 0 (R) is dened in the following way. (w) := W + ( )(? y 0 ) + + ( m + m )(? y m ); where W is the image of w = w 0 + +w m+1 M(R) by the usual injection of the space of measures into the space of distributions. Thus, (w) is a sum of images of two measures. The rst one is a measure which has zero mass at points y 0 ; : : :; y m and equals w w m+1 when restricted to Rn fy 0 ; : : :; y m g. The second one is a linear combination of delta measures in the points fy 0 ; : : :; y m g. We say that (u; v) P L (R) M F (R) is a solution to (16), (17) if?y(u) 0 + (f 1 (u)v + f (u)) 0 = 0?y(v) 0 + ((g 1 (u)v + g (u))) 0 = 0 in D 0. With such interpretation, for given (u 0 ; v 0 ), the delta locus will be the set of all pairs (u 1 ; v 1 ) R such that there exists a solution in the above sense which is constant except in one point. Now we will nd the delta locus (with the above denition) for system (16), (17) and (18). Suppose that y 0 = c is the point of discontinuity and let Then u(y) = u0 ; y < c u 1 ; y > c ; v = w +? + + ; where w(y) = f 1 (u)v + f (u) = f 1 (u)w + f (u) + f 1 (u)(? + + ) v0 ; y < c v 1 ; y > c : = f 1 (u)w + f (u) + f 1 (u 0 )? + f 1 (u 1 ) + : After the injection of the above term in D 0, equation (16) becomes?yu 0 + (f 1 (u)w + f (u)) 0 + (f 1 (u 0 ) + f 1 (u 1 )) 0 = 0: The term?yu 0 + (f 1 (u)w + f (u)) 0 = const determines the shock wave speed c as usual, i.e. const = 0. So, the term (f 1 (u 0 ) + f 1 (u 1 )) 0 has to be zero, i.e. we have (19) c = f 1(u 1 )v 1 + f (u 1 )? f 1 (u 0 )v 0? f (u 0 ) u 1? u 0

13 1 M. NEDELJKOV and (0) f 1 (u 0 ) + f 1 (u 1 ) = 0: Equation (19) denes the speed of the shock and (0) gives the rst condition on and. Equation (17) is transformed in the following way. The rst term, yv = yw + y? +y +, is mapped into yw+y(+) D 0. The second term, g 1 (u)v+g (u) = g 1 (u)w+g (u)+g 1 (u 0 )? +g 1 (u 1 ) +, is mapped into g 1 (u)w+g (u)+(g 1 (u 0 ) + g 1 (u 1 )) D 0. After the substitution of these elements of D 0 into (17) we obtain?yw 0 + (g 1 (u)w + g ) 0? c( + ) 0 (y? c) + ( + )(y? c) +(g 1 (u 0 ) + g 1 (u 1 )) 0 (y? c) = 0: Since?yw 0 + (g 1 (u)w + g ) 0 =?(y? c), where is Rankine-Hugoniot decit, = c(v 1? v 0 )? (g 1 (u 1 )v 1 + g (u 1 )? g 1 (u 0 )v 0? g (u 0 )), and satisfy (1) () + = g 1 (u 0 ) + g 1 (u 1 ) = c: The solution to system (0), (1) and () in respect of and exists if f 1 (u 0 ) f 1 (u 1 ) 0 1 1? g 1 (u 0 ) g 1 (u 1 )?c 6= 0: The result is the same as in the case when the approximations are used. One can see that s = +, D? = =s, and D + = =s. 6. Singular shocks In the case when the uxes f 1, f, g 1 and g are polynomials it is often possible to nd a singular shock solution to system (1)-(3) for a larger set of the initial data. This will be done by adding some terms in u " beside G ". When all functions in the ux are polynomials we can control behaviour of these additional terms in u ". Suppose that the maximal degree of all polynomials in the uxes equals m. Let f 1 (y) = mx i=0 a i y i ; f (y) = mx i=0 a i y i ; g 1 (y) = We will try to nd solutions in the form mx i=0 b i y i ; g (y) = mx i=0 b i y i : (3) u " (x; t) = G " (x? ct) + s 1 (t)d " (x? ct) v " (x; t) = H " (x? ct) + s (t)d " (x? ct); where s i, u = 1; and D " are of the same form as before. A net d " is a representative of the generalized function associated with zero, and it will be described in the proof of the following theorem.

14 DELTA AND SINGULAR DELTA LOCUS 13 Theorem. A constant c = f 1(u 1 )v 1 + f (u 1 )? f 1 (u 0 )v 0? f (u 0 ) u 1? u 0 = [f 1(u)v + f (u)] [u] will be called a speed of the singular shock, and = c(v 1? v 0 )? (g 1 (u 1 )v 1 + g (u 1 )? g 1 (u 0 )v 0? g (u 0 )) = c[v]? [g 1 (u)v + g ] will be called Rankin-Hugoniot decit. A singular shock wave solution to (9)-(11) in the form (3) exists in the following cases. a) Let m be an even number. Then there must be possible to nd > 0 such that there exists a solution (D? ; d? ) R R + to the system (3) [f 1 (u)]d? + [v]a m d? = a m + f 1 (u 1 ) + v 1 a m [g 1 (u)]d? + [v]b m d? = b m + g 1 (u 1 ) + v 1 b m? c: b) Let m be an odd number. One of the two following conditions has to hold. (i) There exists Rn f0g such that there exists a solution (D? ; d? ) R to system (3). (ii) There exists > 0 such that there exists a solution (D? ; d? ) R R + to the system (4) [f 1 (u)]d? + (a m?1 [u] + ma m [uv] + ma m [u])d? =(a m v 1 + ma m u 1 v 1 + ma m u 1 ) + f 1 (u 1 ) + a m?1 [g 1 (u)]d? + (b m?1 [u] + mb m [uv] + mb m [u])d? =(b m v 1 + mb m u 1 v 1 + mb m u 1 ) + g 1 (u 1 ) + b m?1? c: The set of all points (u 1 ; v 1 ) R for which there exists a singular shock wave solution is called the singular delta locus. d " = d? " b?a Proof. a) Let C 1 0 (R), supp [?1; 1], R (x)dx = 1, 0 and let x + a+b " b?a!! 1=m + d + " b?a x? a+b " b?a!! 1=m ; where (d? ) m + (d + ) m = 1. Let s 1 (t); s (t) be smooth functions for t 0, s 1 (0) = s (0) = 0, and let D " be of the same form as in Theorem 1. If G " is a representative of a generalized function of the bounded type (some power of a step function, for example), then Z G " (x? ct)d j "(x? ct) (x; t)dxdt C GC Z " 1?j=m 1=m (y)dydt The last term converge to zero as "! 0, for j < m. Since the multiplication with a smooth function preserves association, there holds (G " + s(t)d " ) j G j "; j < m and (G " + s(t)d " ) m G m " + s m (t)d m " :

15 14 M. NEDELJKOV Let G " (y) be a generalized step function which equals u 0 for y <?" and equals u 1 for y > " One can see that G " d m " (u 0 d? + u 1 d + )(x? ct) from the construction of d ". Also, the P construction of D " i d " implies that D " d " 0. m If?(y) = i=0 iy i, then?(g " + s(t)d " )?(G " ) + s m (t) m d m " : After substitution of U " = G " +s 1 (t)d " and V " = H " +s (t)d " in system (9)-(11), one can see that from the rst equation that c = [f1(u)v+f(u)] (Rankin-Hugoniot [u] condition) and from the second that s (t) = t, where = c[v]? [g 1 (u)v + g ] is a Rankin-Hugoniot decit. Further on, after association procedure like in the proof of Theorem 1 and grouping coecients in respect of and 0 one nds that s 1 (t) = ~st 1=m and (recall that D? + D + = (d? ) m + (d + ) m = 1) [f 1 (u)]d? + [v]a m ~s m d? = a m ~s m + f 1 (u 1 ) + v 1 a m ~s m [g 1 (u)]d? + [v]b m ~s m d? = b m ~s m + g 1 (u 1 ) + v 1 b m ~s m? c: This proves the rst part of the theorem. b) (i) The proof of the assertion is the same as the one of a). The only dierence is that ~s m has not to be positive. (ii) The construction of d " in the second case is a slightly dierent. Let be as above and d? 1=(m?1) x + (a + 3b)=4 d? 1=(m?1) x + (3a + b)=4 d " = "(b? a)? "(b? a)=4 "(b? a) "(b? a)=4 d + x? 1=(m?1) (3a + b)=4 d + x? 1=(m?1) (a + 3b)=4 + "(b? a) + "(b? a)=4 "(b? a) ; "(b? a)=4 where (d? ) m?1 + (d + ) m?1 = 1. Like in the previous case, one can see that (G " + s(t)d " ) j G j "; j < m? 1 and (G " + s(t)d " ) m?1 G m?1 " + s m?1 (t)d m?1 " for every generalized function G " of the bounded type. Because of the specic denition of d " one can see that G " d m " 0 if G " is a generalized step function. The above means that?(g " + s 1 (t)d " )?(G " ) + m?1 s m?1 1 (t)d m?1 " + m m s m?1 1 (t)g " d m?1 " : Substituting U " = G " + s 1 (t)d " and V " = H " + s (t)d " in system (9)-(11) one can see that c = [f1(u)v+f(u)], s [u] (t) = t, where = c[v]? [g 1 (u)v + g ] like in a). Also, one nds that s 1 (t) = ~st 1=(m?1) and [f 1 (u)]d? + (a m?1 [u] + ma m [uv] + ma m [u])~s m?1 d? =~s m?1 (a m v 1 + ma m u 1 v 1 + ma m u 1 ) + f 1 (u 1 ) + a m?1 ~s m?1 [g 1 (u)]d? + (b m?1 [u] + mb m [uv] + mb m [u])~s m?1 d? =~s m?1 (b m v 1 + mb m u 1 v 1 + mb m u 1 ) + g 1 (u 1 ) + b m?1 ~s m?1? c: This proves the theorem. In the following corollaries we shall describe singular delta locus in the case f 1 a where a is a nonzero constant. This is the case when the delta locus is the empty set (see Theorem 1).

16 DELTA AND SINGULAR DELTA LOCUS 15 Corollary 1. Let f 1 a, a 6= 0 is a constant, and let g 1 be non constant function. a) If m is an even number, then (u 1 ; v 1 ) is in a delta singular locus if a m 6= 0 (i.e. deg f = m), there exists > 0 such that sgn([v]a m ) = sgn(a + v 1 a m ) and [g 1 (u)] 6= 0. b) If m is an odd number and a m 6= 0, then the delta singular locus is the set R n fu 1 R: [g 1 (u)] = 0g. c) Let m be an odd number, a m = 0 and a m?1 6= 0 (i.e. deg f = m? 1). Then (u 1 ; v 1 ) is in a delta singular locus if a=a m?1 < 0 and [g 1 (u)] 6= 0. Let us note that in the case deg f < m? 1 the delta singular locus is the empty set ( = 0 is the condition for Rankin-Hugoniot locus). Proof. In all of these cases, the condition [g 1 (u)] 6= 0 ensures that there always exists D? if appropriate d? exists. a) The proof easily follows from (4) since D? is an arbitrary real number. b) The proof is an immediate consequence of b) (i) in the previous theorem. c) One can see from the rst equation in (4) that =?a=a m?1 > 0. Corollary. Let f 1 a and g 1 b, where a 6= 0 and b are some constants. a) Let a m and b m be dierent from zero. Then = ((b? c)a m? ab m ) a m?1 b m? a m b m?1 has to be greater than zero, and (u 1 ; v 1 ) is in delta singular locus if d? = ma mu 1 + a + a m?1 ma m [u] [0; 1] for dened above. b) Let a m 6= 0 and b m = 0. Then = c? b b m?1 has to be positive, and (u 1 ; v 1 ) is in a delta locus if d? = ma mu 1 + a + a m?1 ma m [u] [0; 1] for dened above. c) Let a m = 0 and b m 6= 0. Then =?a a m?1 has to be positive, and (u 1 ; v 1 ) is in a delta locus if d? = mb mu 1 + b + b m?1? c ma m [u] [0; 1] for dened above.

17 16 M. NEDELJKOV d) Let a m = b m = 0. Then a delta singular locus exist, and it is R, only if c? b b m?1 =?a a m?1 > 0: Proof. In all of these cases we are using the second denition of d " and the assertions follow by solving the system of equations (4) with the unknowns and d?. Remark 1. As in the rst section all the calculation can be made for box approximations. Remark. In [6] the singular parts of functions u and v have nonempty intersection (even d " = D " ). This is possible because in the system (4) there are no multiplication of d " and D ", i.e. f 1 =?1, g 1 = 0 are the constants. It is also possible to choose d " and D " with non-disjoint supports if ux functions are linear in u. In general case, problems with terms D " d j " can appear. The solution of Keytz and Kranzer is the same as in the case c) of Corrolary for m = 3, a =?1, b = 0, a m?1 = 1, b m = 1=3 and b m?1 = 0. When all points for which the overcompressibility condition does not hold are extracted, one obtain the same area as them (it is denoted by Q 7 in their paper). Acknowledgement The author is grateful to Barbara Keytz, Michael Oberguggenberger and Stevan Pilipovic for their great help in the preparation of this paper. He would like to thank to Austrian Academic Exchange Service ( OAD) for nancing his visit to Austria during which the paper was written. Also he is thankful to Institut f. Geometrie u. Mathematik, Univeristat Innsbruck and E. Schrodinger Institute, Wien, for their great hospitality. References 1. J. F. Colombeau, Elementary Introduction in New Generalized Functions, North Holland, C. M. Dafermos, R. J. DiPerna, The Riemann problem for certain classes of hyperbolic systems of conservation laws, J. Dierential Equations 0 (1976), V. G. Danilov, V. P. Maslov, V. M. Shelkovich, Algebras of singularities of singular solutions to rst-order quasilinear strictly hyperbolic systems, Preprint (1998). 4. G. Ercole, Delta shock waves as self-similar viscosity limits, Preprint (1997). 5. B. T. Hayes, P. G. Le Flock, Measure solutions to a strictly hyperbolic system of conservation laws, Nonlinearity 9 (1996), B. L. Keytz, H. C. Kranzer, Spaces of weighted measures for conservation laws with singular shock solutions, J. Di. Eq. 118, (1995), P. D. Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM, Philadelphia, V. P. Maslov, G. A. Omel'anov, Asymptotic solitonlike solutions of equations with small dispersion, Russian Math. Surveys 36 (1981), M. Oberguggenberger, Multiplication of Distributions and Applications to Partial Dierential Equations, Pitman Res. Not. Math. 59, Longman Sci. Techn., Essex, M. Oberguggenberger, Y.-G. Wang, Generalized solutions to conservation laws, Zeitschr. Anal. Anw. 13 (1994), D. Tan, T. Zhang, Y. Zheng, Delta-shock waves as limits of vanishing viscosity for hyperbolic systems of conservation laws, J. Di. Eq. 11 (1994), Y. Zheng, Systems of conservation laws with incomplete sets of eigenvectors everywhere, in A Volume Dedicated to Professor Xiaqi Ding's 70th Birthay (to appear).

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-19 Wien, Austria The Negative Discrete Spectrum of a Class of Two{Dimentional Schrodinger Operators with Magnetic

More information

Pointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang

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

WEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II

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

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On the Point Spectrum of Dirence Schrodinger Operators Vladimir Buslaev Alexander Fedotov

More information

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

Green s Functions and Distributions

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

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On an Extension of the 3{Dimensional Toda Lattice Mircea Puta Vienna, Preprint ESI 165 (1994)

More information

Blowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari

Blowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari Blowup for Hyperbolic Equations Helge Kristian Jenssen and Carlo Sinestrari Abstract. We consider dierent situations of blowup in sup-norm for hyperbolic equations. For scalar conservation laws with a

More information

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................

More information

DELTA SHOCK WAVE FORMATION IN THE CASE OF TRIANGULAR HYPERBOLIC SYSTEM OF CONSERVATION LAWS

DELTA SHOCK WAVE FORMATION IN THE CASE OF TRIANGULAR HYPERBOLIC SYSTEM OF CONSERVATION LAWS DELTA SHOCK WAVE FORMATION IN THE CASE OF TRIANGULAR HYPERBOLIC SYSTEM OF CONSERVATION LAWS V.G. DANILOV AND D. MITROVIC Abstract. We describe δ-shock wave generation from continuous initial data in the

More information

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-090 Wien, Austria On the Rigid Body with Two Linear Controls Mircea Craioveanu Mircea Puta Vienna, Preprint

More information

ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL

ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,

More information

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

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

8 Singular Integral Operators and L p -Regularity Theory

8 Singular Integral Operators and L p -Regularity Theory 8 Singular Integral Operators and L p -Regularity Theory 8. Motivation See hand-written notes! 8.2 Mikhlin Multiplier Theorem Recall that the Fourier transformation F and the inverse Fourier transformation

More information

STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION

STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION Electronic Journal of Differential Equations, Vol. 216 (216, No. 126, pp. 1 16. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN

More information

Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D

Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, U.K. R.J. Douglas, Isaac Newton

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Common Homoclinic Points of Commuting Toral Automorphisms Anthony Manning Klaus Schmidt

More information

Applications of the compensated compactness method on hyperbolic conservation systems

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

growth rates of perturbed time-varying linear systems, [14]. For this setup it is also necessary to study discrete-time systems with a transition map

growth rates of perturbed time-varying linear systems, [14]. For this setup it is also necessary to study discrete-time systems with a transition map Remarks on universal nonsingular controls for discrete-time systems Eduardo D. Sontag a and Fabian R. Wirth b;1 a Department of Mathematics, Rutgers University, New Brunswick, NJ 08903, b sontag@hilbert.rutgers.edu

More information

Contents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces

Contents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and

More information

Math Partial Differential Equations 1

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

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria Spherically Symmetric Solutions of Nonlinear Schrodinger Equations Roman Cherniha Vienna,

More information

r( = f 2 L 2 (1.2) iku)! 0 as r!1: (1.3) It was shown in book [7] that if f is assumed to be the restriction of a function in C

r(  = f 2 L 2 (1.2) iku)! 0 as r!1: (1.3) It was shown in book [7] that if f is assumed to be the restriction of a function in C Inverse Obstacle Problem: Local Uniqueness for Rougher Obstacles and the Identication of A Ball Changmei Liu Department of Mathematics University of North Carolina Chapel Hill, NC 27599 December, 1995

More information

290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f

290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica

More information

Congurations of periodic orbits for equations with delayed positive feedback

Congurations of periodic orbits for equations with delayed positive feedback Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Group Analysis of Dierential Equations and Generalized Functions M. Kunzinger M. Oberguggenberger

More information

Math 145. Codimension

Math 145. Codimension Math 145. Codimension 1. Main result and some interesting examples In class we have seen that the dimension theory of an affine variety (irreducible!) is linked to the structure of the function field in

More information

MATH 220: Problem Set 3 Solutions

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

Delta-shock Wave Type Solution of Hyperbolic Systems of Conservation Laws

Delta-shock Wave Type Solution of Hyperbolic Systems of Conservation Laws Delta-shock Wave Type Solution of Hyperbolic Systems of Conservation Laws V. G. Danilov and V. M. Shelkovich Abstract. For some classes of hyperbolic systems of conservation laws we introduce a new definition

More information

2 Introduction T(,t) = = L Figure..: Geometry of the prismatical rod of Eample... T = u(,t) = L T Figure..2: Geometry of the taut string of Eample..2.

2 Introduction T(,t) = = L Figure..: Geometry of the prismatical rod of Eample... T = u(,t) = L T Figure..2: Geometry of the taut string of Eample..2. Chapter Introduction. Problems Leading to Partial Dierential Equations Many problems in science and engineering are modeled and analyzed using systems of partial dierential equations. Realistic models

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Hypersymmetries of Extremal Equations D.P. Zhelobenko Vienna, Preprint ESI 391 (1996) October

More information

Abstract. A front tracking method is used to construct weak solutions to

Abstract. A front tracking method is used to construct weak solutions to A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria A Nonlinear Theory of Tensor Distributions J.A. Vickers J.P. Wilson Vienna, Preprint ESI

More information

Part V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory

Part V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory Part V 7 Introduction: What are measures and why measurable sets Lebesgue Integration Theory Definition 7. (Preliminary). A measure on a set is a function :2 [ ] such that. () = 2. If { } = is a finite

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Noncommutative Contact Algebras Hideki Omori Yoshiaki Maeda Naoya Miyazaki Akira Yoshioka

More information

Remarks on the embedding of spaces of distributions into spaces of Colombeau generalized functions

Remarks on the embedding of spaces of distributions into spaces of Colombeau generalized functions Remarks on the embedding of spaces of distributions into spaces of Colombeau generalized functions Antoine Delcroix To cite this version: Antoine Delcroix. Remarks on the embedding of spaces of distributions

More information

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition) Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria On the Foundations of Nonlinear Generalized Functions I E. Farkas M. Grosser M. Kunzinger

More information

MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ

MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ GENERALIZED UNIFORMLY CONTINUOUS SEMIGROUPS AND SEMILINEAR HYPERBOLIC SYSTEMS WITH REGULARIZED DERIVATIVES MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ Abstract. We aopt the theory of uniformly continuous operator

More information

Spurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics

Spurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics

More information

Downloaded 03/01/17 to Redistribution subject to SIAM license or copyright; see

Downloaded 03/01/17 to Redistribution subject to SIAM license or copyright; see SIAM J. DISCRETE MATH. Vol. 31, No. 1, pp. 335 382 c 2017 Society for Industrial and Applied Mathematics PARTITION CONSTRAINED COVERING OF A SYMMETRIC CROSSING SUPERMODULAR FUNCTION BY A GRAPH ATTILA BERNÁTH,

More information

Analysis on Graphs. Alexander Grigoryan Lecture Notes. University of Bielefeld, WS 2011/12

Analysis on Graphs. Alexander Grigoryan Lecture Notes. University of Bielefeld, WS 2011/12 Analysis on Graphs Alexander Grigoryan Lecture Notes University of Bielefeld, WS 0/ Contents The Laplace operator on graphs 5. The notion of a graph............................. 5. Cayley graphs..................................

More information

Simple Lie subalgebras of locally nite associative algebras

Simple Lie subalgebras of locally nite associative algebras Simple Lie subalgebras of locally nite associative algebras Y.A. Bahturin Department of Mathematics and Statistics Memorial University of Newfoundland St. John's, NL, A1C5S7, Canada A.A. Baranov Department

More information

Discrete Halanay-type inequalities and applications

Discrete Halanay-type inequalities and applications Nonlinear Analysis 55 (2003) 669 678 www.elsevier.com/locate/na Discrete Halanay-type inequalities and applications Eduardo Liz a;, Anatoli Ivanov b, Juan Bosco Ferreiro c a Departamento de Matematica

More information

The Riemann problem. The Riemann problem Rarefaction waves and shock waves

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

Bulletin T.CXXXIII de l Académie serbe des sciences et des arts 2006 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 31

Bulletin T.CXXXIII de l Académie serbe des sciences et des arts 2006 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 31 Bulletin T.CXXXIII de l Académie serbe des sciences et des arts 2006 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 31 GENERALIZED SOLUTIONS TO SINGULAR INITIAL-BOUNDARY HYPERBOLIC

More information

A NOTE ON MATRIX REFINEMENT EQUATIONS. Abstract. Renement equations involving matrix masks are receiving a lot of attention these days.

A NOTE ON MATRIX REFINEMENT EQUATIONS. Abstract. Renement equations involving matrix masks are receiving a lot of attention these days. A NOTE ON MATRI REFINEMENT EQUATIONS THOMAS A. HOGAN y Abstract. Renement equations involving matrix masks are receiving a lot of attention these days. They can play a central role in the study of renable

More information

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

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

More information

STOCHASTIC DIFFERENTIAL EQUATIONS DRIVEN BY GENERALIZED POSITIVE NOISE. Michael Oberguggenberger and Danijela Rajter-Ćirić

STOCHASTIC DIFFERENTIAL EQUATIONS DRIVEN BY GENERALIZED POSITIVE NOISE. Michael Oberguggenberger and Danijela Rajter-Ćirić PUBLICATIONS DE L INSTITUT MATHÉMATIQUE Nouvelle série, tome 7791 25, 7 19 STOCHASTIC DIFFERENTIAL EQUATIONS DRIVEN BY GENERALIZED POSITIVE NOISE Michael Oberguggenberger and Danijela Rajter-Ćirić Communicated

More information

Introduction to Arithmetic Geometry Fall 2013 Lecture #23 11/26/2013

Introduction to Arithmetic Geometry Fall 2013 Lecture #23 11/26/2013 18.782 Introduction to Arithmetic Geometry Fall 2013 Lecture #23 11/26/2013 As usual, a curve is a smooth projective (geometrically irreducible) variety of dimension one and k is a perfect field. 23.1

More information

A Note On Minimum Path Bases

A Note On Minimum Path Bases A Note On Minimum Path Bases Petra M. Gleiss a, Josef Leydold b, Peter F. Stadler a,c a Institute for Theoretical Chemistry and Structural Biology, University of Vienna, Währingerstrasse 17, A-1090 Vienna,

More information

Entrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems

Entrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems September, 27 Solve exactly 6 out of the 8 problems. Prove by denition (in ɛ δ language) that f(x) = + x 2 is uniformly continuous in (, ). Is f(x) uniformly continuous in (, )? Prove your conclusion.

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical hysics A-090 Wien, Austria A Note on the enrose Junction Conditions Michael Kunzinger Roland Steinbauer Vienna, reprint

More information

1 Holomorphic functions

1 Holomorphic functions Robert Oeckl CA NOTES 1 15/09/2009 1 1 Holomorphic functions 11 The complex derivative The basic objects of complex analysis are the holomorphic functions These are functions that posses a complex derivative

More information

DELTA-SHOCK WAVE TYPE SOLUTION OF HYPERBOLIC SYSTEMS OF CONSERVATION LAWS

DELTA-SHOCK WAVE TYPE SOLUTION OF HYPERBOLIC SYSTEMS OF CONSERVATION LAWS QUARTERLY OF APPLIED MATHEMATICS VOLUME LXIII, NUMBER 3 SEPTEMBER 25, PAGES 41 427 S 33-569X5961-8 Article electronically published on August 17, 25 DELTA-SHOCK WAVE TYPE SOLUTION OF HYPERBOLIC SYSTEMS

More information

LEBESGUE INTEGRATION. Introduction

LEBESGUE INTEGRATION. Introduction LEBESGUE INTEGATION EYE SJAMAA Supplementary notes Math 414, Spring 25 Introduction The following heuristic argument is at the basis of the denition of the Lebesgue integral. This argument will be imprecise,

More information

Stochastic dominance with imprecise information

Stochastic dominance with imprecise information Stochastic dominance with imprecise information Ignacio Montes, Enrique Miranda, Susana Montes University of Oviedo, Dep. of Statistics and Operations Research. Abstract Stochastic dominance, which is

More information

STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER. Department of Mathematics. University of Wisconsin.

STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER. Department of Mathematics. University of Wisconsin. STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER Department of Mathematics University of Wisconsin Madison WI 5376 keisler@math.wisc.edu 1. Introduction The Loeb measure construction

More information

DELTA SHOCK WAVES AS A LIMIT OF SHOCK WAVES

DELTA SHOCK WAVES AS A LIMIT OF SHOCK WAVES st Reading October 0, 200 :2 WSPC/JHDE 002 Journal of Hyperbolic Differential Equations Vol. 4, No. 4 (200) 2 c World Scientific Publishing Company DELTA SHOCK WAVES AS A LIMIT OF SHOCK WAVES DARKO MITROVIĆ

More information

Lifting to non-integral idempotents

Lifting to non-integral idempotents Journal of Pure and Applied Algebra 162 (2001) 359 366 www.elsevier.com/locate/jpaa Lifting to non-integral idempotents Georey R. Robinson School of Mathematics and Statistics, University of Birmingham,

More information

Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr

Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL

More information

Markov processes Course note 2. Martingale problems, recurrence properties of discrete time chains.

Markov processes Course note 2. Martingale problems, recurrence properties of discrete time chains. Institute for Applied Mathematics WS17/18 Massimiliano Gubinelli Markov processes Course note 2. Martingale problems, recurrence properties of discrete time chains. [version 1, 2017.11.1] We introduce

More information

Theorem 1. Suppose that is a Fuchsian group of the second kind. Then S( ) n J 6= and J( ) n T 6= : Corollary. When is of the second kind, the Bers con

Theorem 1. Suppose that is a Fuchsian group of the second kind. Then S( ) n J 6= and J( ) n T 6= : Corollary. When is of the second kind, the Bers con ON THE BERS CONJECTURE FOR FUCHSIAN GROUPS OF THE SECOND KIND DEDICATED TO PROFESSOR TATSUO FUJI'I'E ON HIS SIXTIETH BIRTHDAY Toshiyuki Sugawa x1. Introduction. Supposebthat D is a simply connected domain

More information

Part III. 10 Topological Space Basics. Topological Spaces

Part III. 10 Topological Space Basics. Topological Spaces Part III 10 Topological Space Basics Topological Spaces Using the metric space results above as motivation we will axiomatize the notion of being an open set to more general settings. Definition 10.1.

More information

J. Huisman. Abstract. let bxc be the fundamental Z=2Z-homology class of X. We show that

J. Huisman. Abstract. let bxc be the fundamental Z=2Z-homology class of X. We show that On the dual of a real analytic hypersurface J. Huisman Abstract Let f be an immersion of a compact connected smooth real analytic variety X of dimension n into real projective space P n+1 (R). We say that

More information

Notes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow

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 information

2. Function spaces and approximation

2. Function spaces and approximation 2.1 2. Function spaces and approximation 2.1. The space of test functions. Notation and prerequisites are collected in Appendix A. Let Ω be an open subset of R n. The space C0 (Ω), consisting of the C

More information

MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione

MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,

More information

Quantum logics with given centres and variable state spaces Mirko Navara 1, Pavel Ptak 2 Abstract We ask which logics with a given centre allow for en

Quantum logics with given centres and variable state spaces Mirko Navara 1, Pavel Ptak 2 Abstract We ask which logics with a given centre allow for en Quantum logics with given centres and variable state spaces Mirko Navara 1, Pavel Ptak 2 Abstract We ask which logics with a given centre allow for enlargements with an arbitrary state space. We show in

More information

D l. S * S r. S l. Figure 2: Domain of integration. Figure 3: Entropy structure across inviscid shock wave.

D l. S * S r. S l. Figure 2: Domain of integration. Figure 3: Entropy structure across inviscid shock wave. t D l Dr x Figure 2: Domain of integration. S S * S r S l x Figure 3: Entropy structure across inviscid shock wave. 14 10 2 21 18 15 S 12 9 6 3 0 5 3 1 x λ o 1 3 5 Figure 1: Entropy proles through viscous

More information

Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always

Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.

More information

vestigated. It is well known that the initial transition to convection is accompanied by the appearance of (almost) regular patterns. s the temperatur

vestigated. It is well known that the initial transition to convection is accompanied by the appearance of (almost) regular patterns. s the temperatur Symmetry of ttractors and the Karhunen-Loeve Decomposition Michael Dellnitz Department of Mathematics University of Hamburg D-2046 Hamburg Martin Golubitsky Department of Mathematics University of Houston

More information

Lecture 2: Review of Prerequisites. Table of contents

Lecture 2: Review of Prerequisites. Table of contents Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace s Equation. Chapter Mean Value Formulas Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic

More information

SOME MEASURABILITY AND CONTINUITY PROPERTIES OF ARBITRARY REAL FUNCTIONS

SOME MEASURABILITY AND CONTINUITY PROPERTIES OF ARBITRARY REAL FUNCTIONS LE MATEMATICHE Vol. LVII (2002) Fasc. I, pp. 6382 SOME MEASURABILITY AND CONTINUITY PROPERTIES OF ARBITRARY REAL FUNCTIONS VITTORINO PATA - ALFONSO VILLANI Given an arbitrary real function f, the set D

More information

arxiv: v2 [math.ap] 1 Jul 2011

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

2 SCHECTER is a step function, then in the sense of distributions, f(u) x is the measure (f(u + ) f(u ))(x): (1.6) This leads to the Rankine-Hugoniot

2 SCHECTER is a step function, then in the sense of distributions, f(u) x is the measure (f(u + ) f(u ))(x): (1.6) This leads to the Rankine-Hugoniot TRAVELING-WAVE SOLUTIONS OF CONVECTION-DIFFUSION SYSTEMS BY CENTER MANIFOLD REDUCTION STEPHEN SCHECTER Abstract. Traveling waves u(x st) for systems of conservation laws u t + Df(u)u x = (B(u)u x ) x were

More information

LECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1

LECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2  1 LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,

More information

Mathematical Institute, University of Utrecht. The problem of estimating the mean of an observed Gaussian innite-dimensional vector

Mathematical Institute, University of Utrecht. The problem of estimating the mean of an observed Gaussian innite-dimensional vector On Minimax Filtering over Ellipsoids Eduard N. Belitser and Boris Y. Levit Mathematical Institute, University of Utrecht Budapestlaan 6, 3584 CD Utrecht, The Netherlands The problem of estimating the mean

More information

Mathematical Methods for Physics and Engineering

Mathematical Methods for Physics and Engineering Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory

More information

Group Theory. 1. Show that Φ maps a conjugacy class of G into a conjugacy class of G.

Group Theory. 1. Show that Φ maps a conjugacy class of G into a conjugacy class of G. Group Theory Jan 2012 #6 Prove that if G is a nonabelian group, then G/Z(G) is not cyclic. Aug 2011 #9 (Jan 2010 #5) Prove that any group of order p 2 is an abelian group. Jan 2012 #7 G is nonabelian nite

More information

Partial Differential Equations

Partial Differential Equations Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,

More information

Existence Theory for Hyperbolic Systems of Conservation Laws with General Flux-Functions

Existence Theory for Hyperbolic Systems of Conservation Laws with General Flux-Functions Existence Theory for Hyperbolic Systems of Conservation Laws with General Flux-Functions Tatsuo Iguchi & Philippe G. LeFloch Abstract For the Cauchy problem associated with a nonlinear, strictly hyperbolic

More information

PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS

PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix

More information

Classification of homogeneous quadratic conservation laws with viscous terms

Classification of homogeneous quadratic conservation laws with viscous terms Volume 26, N. 2, pp. 251 283, 2007 Copyright 2007 SBMAC ISSN 0101-8205 www.scielo.br/cam Classification of homogeneous quadratic conservation laws with viscous terms JANE HURLEY WENSTROM 1 and BRADLEY

More information

On Coarse Geometry and Coarse Embeddability

On Coarse Geometry and Coarse Embeddability On Coarse Geometry and Coarse Embeddability Ilmari Kangasniemi August 10, 2016 Master's Thesis University of Helsinki Faculty of Science Department of Mathematics and Statistics Supervised by Erik Elfving

More information

Notes on Complex Analysis

Notes on Complex Analysis Michael Papadimitrakis Notes on Complex Analysis Department of Mathematics University of Crete Contents The complex plane.. The complex plane...................................2 Argument and polar representation.........................

More information

The Rationality of Certain Moduli Spaces of Curves of Genus 3

The Rationality of Certain Moduli Spaces of Curves of Genus 3 The Rationality of Certain Moduli Spaces of Curves of Genus 3 Ingrid Bauer and Fabrizio Catanese Mathematisches Institut Universität Bayreuth, NW II D-95440 Bayreuth, Germany Ingrid.Bauer@uni-bayreuth.de,

More information

October 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable

October 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable International Journal of Wavelets, Multiresolution and Information Processing c World Scientic Publishing Company Polynomial functions are renable Henning Thielemann Institut für Informatik Martin-Luther-Universität

More information

INTRODUCTION TO NETS. limits to coincide, since it can be deduced: i.e. x

INTRODUCTION TO NETS. limits to coincide, since it can be deduced: i.e. x INTRODUCTION TO NETS TOMMASO RUSSO 1. Sequences do not describe the topology The goal of this rst part is to justify via some examples the fact that sequences are not sucient to describe a topological

More information

(dim Z j dim Z j 1 ) 1 j i

(dim Z j dim Z j 1 ) 1 j i Math 210B. Codimension 1. Main result and some interesting examples Let k be a field, and A a domain finitely generated k-algebra. In class we have seen that the dimension theory of A is linked to the

More information

On the piecewise smoothness of entropy solutions to scalar conservation laws for a large class of initial data

On the piecewise smoothness of entropy solutions to scalar conservation laws for a large class of initial data On the piecewise smoothness of entropy solutions to scalar conservation laws for a large class of initial data Tao Tang Jinghua Wang Yinchuan Zhao Abstract We prove that if the initial data do not belong

More information

van Rooij, Schikhof: A Second Course on Real Functions

van Rooij, Schikhof: A Second Course on Real Functions vanrooijschikhof.tex April 25, 2018 van Rooij, Schikhof: A Second Course on Real Functions Notes from [vrs]. Introduction A monotone function is Riemann integrable. A continuous function is Riemann integrable.

More information

On the Front-Tracking Algorithm

On the Front-Tracking Algorithm JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 7, 395404 998 ARTICLE NO. AY97575 On the Front-Tracking Algorithm Paolo Baiti S.I.S.S.A., Via Beirut 4, Trieste 3404, Italy and Helge Kristian Jenssen

More information

Non-linear Scalar Equations

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

Denition and some Properties of Generalized Elementary Functions of a Real Variable

Denition and some Properties of Generalized Elementary Functions of a Real Variable Denition and some Properties of Generalized Elementary Functions of a Real Variable I. Introduction The term elementary function is very often mentioned in many math classes and in books, e.g. Calculus

More information

Simple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables

Simple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables s and characteristic decompositions of quasilinear hyperbolic systems in two independent variables Wancheng Sheng Department of Mathematics, Shanghai University (Joint with Yanbo Hu) Joint Workshop on

More information

Solutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness

Solutions 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

A Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract

A Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin

More information

2 Garrett: `A Good Spectral Theorem' 1. von Neumann algebras, density theorem The commutant of a subring S of a ring R is S 0 = fr 2 R : rs = sr; 8s 2

2 Garrett: `A Good Spectral Theorem' 1. von Neumann algebras, density theorem The commutant of a subring S of a ring R is S 0 = fr 2 R : rs = sr; 8s 2 1 A Good Spectral Theorem c1996, Paul Garrett, garrett@math.umn.edu version February 12, 1996 1 Measurable Hilbert bundles Measurable Banach bundles Direct integrals of Hilbert spaces Trivializing Hilbert

More information