Balance laws with integrable unbounded sources
|
|
- Malcolm Gardner
- 5 years ago
- Views:
Transcription
1 Università di Milano Bicocca Quaderni di Matematica Balance laws with integrable unbounded sources Graziano Guerra, Francesca Marcellini and Veronika Schleper Quaderno n. 8/28 arxiv: v
2 Stampato nel mese di ottobre 28 presso il Dipartimento di Matematica e Applicazioni, Università degli Studi di Milano-Bicocca, via R. Cozzi 53, 225 Milano, ITALIA. Disponibile in formato elettronico sul sito Segreteria di redazione: Ada Osmetti - Giuseppina Cogliandro tel.: fax: Esemplare fuori commercio per il deposito legale agli effetti della Legge 5 aprile 24 n.6.
3 Balance laws with integrable unbounded sources arxiv: v [math.ap] 6 Sep 28 Graziano Guerra Francesca Marcellini Veronika Schleper September 6, 28 Abstract We consider the Cauchy problem for a n n strictly hyperbolic system of balance laws 8 >< >: u t + fu x = gx,u, x R, t > u,. = u o L BV R; R n, λ iu c > for all i {,..., n}, gx, C 2 Mx L, each characteristic field being genuinely nonlinear or linearly degenerate. Assuming that the L norm of gx, C and u o BVR are small enough, we prove the existence and uniqueness of global entropy solutions of bounded total variation extending the result in [] to unbounded in L sources. Furthermore, we apply this result to the fluid flow in a pipe with discontinuous cross sectional area, showing existence and uniqueness of the underlying semigroup. 2 Mathematics Subject Classification: 35L65, 35L45, 35L6. Keywords: Hyperbolic Balance Laws, Unbounded Sources, Pipes with Discontinuous Cross Sections. Introduction The recent literature offers several results on the properties of gas flows on networks. For instance, in [4, 5, 6, 8] the well posedness is established for the gas flow at a junction of n pipes with constant diameters. The equations governing the gas flow in a pipe with a smooth varying cross section ax are given by see for instance []: This work has been supported by the SPP253 priority program of the DFG and by the DAAD program D/6/9582. Università di Milano Bicocca, Italy, graziano.guerra@unimib.it Università di Milano Bicocca, Italy, f.marcellini@campus.unimib.it TU Kaiserslautern, Germany, sachers@mathematik.uni-kl.de
4 ρ t + q x = a x ax q q t + q2 ρ +p x = a x q 2 ax ρ e t + q e+p ρ x = a x ax q ρ e + p. The well posedness of this system is covered in [] where an attractive unified approach to the existence and uniqueness theory for quasilinear strictly hyperbolic systems of balance laws is proposed. The case of discontinuous cross sections is considered in the literature inserting a junction with suitable coupling conditions at the junction, see for example [4, 5, 9]. One way to obtain coupling conditions at the point of discontinuity of the cross section a is to take the limit of a sequence of Lipschitz continuous cross sections a ε converging to a in L for a different approach see for instance [7]. Unfortunately the results in [] require L bounds on the source term and well posedness is proved on a domain depending on this L bound. Since in the previous equations the source term contains the derivative of the cross sectional area one cannot hope to take the limit a ε a. Indeed when a is discontinuous, the L norm of a ε goes to infinity. Therefore the purpose of this paper is to establish the result in [] without requiring the L bound. More precisely, we consider the Cauchy problem for the following n n system of equations u t + fu x = gx, u, x R, t >, endowed with a suitably small initial data u, x = u o x, x R. 2 belonging to L BV R; R n, the space of integrable functions with bounded total variation Tot.Var. in the sense of [2]. Here ut, x R n is the vector of unknowns, f : Ω R n denotes the fluxes, i.e. a smooth function defined on Ω which is an open neighborhood of the origin in R n. The system is supposed to be strictly hyperbolic, with each characteristic field either genuinely nonlinear or linearly degenerate in the sense of Lax []. Concerning the source term g, we assume that it satisfies the following Caratheodory type conditions: P g : R Ω R n is measurable with respect to w.r.t. x, for any u Ω, and is C 2 w.r.t. u, for any x R; P 2 there exists a L function Mx such that gx, C 2 Mx; P 3 there exists a function ω LR such that gx, C ωx. Remark. Note that the L norm of Mx does not have to be small but only bounded differently from ωx whose norm has to be small see Theorem below. Furthermore condition P 2 replaces the L bound of the C 2 norm 2
5 of g in []. Finally observe that we do not require any L bound on ω. On the other hand we will need the following observation: if we define ε h = sup x R by absolute continuity one has ε h as h. ωx + s ds, 3 Moreover, we assume that a non-resonance condition holds, that is the characteristic speeds of the system are bounded away from zero: λ i u c >, u Ω, i {,..., n}. 4 The following theorem states the well posedness of in the above defined setting. Theorem. Assume P P 3 and 4. If the norm of ω in LR is sufficiently small, there exist a constant L >, a closed domain D of integrable functions with small total variation and a unique semigroup P : [, + D D satisfying i P u = u, P t+s u = P t P s u for all u, v D and t, s ; ii P s u P t v LR L s t + u v LR for all u, v D and t, s ; iii for all u o D the function ut, = P t u o is a weak entropy solution of the Cauchy problem 2 and satisfies the integral estimates 44, 45. Conversely let u : [, T] D be Lipschitz continuous as a map with values in LR, R n and assume that ut, x satisfies the integral conditions 44, 45. Then ut, coincides with a trajectory of the semigroup P. The proof of this theorem is postponed to sections 3 and 4, where existence and uniqueness are proved. Before these technical details, we state the application of the above result to gas flow in section 2. Here we apply Theorem to establish the existence and uniqueness of the semigroup related to pipes with discontinuous cross sections. Furthermore, we show that our approach yields the same semigroup as the approach followed in [6] in the special case of two connected pipes. The technical details of section 2 can be found at the end of the paper in section 5. 2 Application to gas dynamics Theorem provides an existence and uniqueness result for pipes with Lipschitz continuous cross section where the equations governing the gas flow are given by ρ t + q x = a x ax q 2 q t + q ρ +p x = a x q 2 ax ρ e t + q e+p ρ x = a x ax q ρ e + p. 5 3
6 Here, as usual, ρ denotes the mass density, q the linear momentum, e is the energy density, a is the area of the cross section of the pipe and p is the pressure which is related to the conserved quantities ρ, q, e by the equations of state. In most situations, when two pipes of different size have to be connected, the length l of the adaptor is small compared to the length of the pipes. Therefore it is convenient to model these connections as pipes with a jump in the cross sectional area. These discontinuous cross sections however do not fulfill the requirements of Theorem. Nevertheless, we can use this Theorem to derive the existence of solutions to the discontinuous problem by a limit procedure. To this end, we approximate the discontinuous function ax = { a, x < a +, x > by a sequence a l C, R, R + with the following properties a, x < l 2 a l x = ϕ l x, x [ l 2, ] l 2 a +, x > l where ϕ l is any smooth monotone function which connects the two strictly positive constants a, a +. One possible choice of the approximations a l as well as the discontinuous pipe with cross section a are shown in figure. l a a + a a + Figure : Illustration of approximated and discontinuous cross-sectional area With the help of Theorem and the techniques used in its proof, we are now able to derive the following Theorem see also [7] for a similar result obtained with different methods. Theorem 2. If a l L = a+ a is sufficiently small, the semigroups P l related with the smooth section a l converge to a unique semigroup P. The limit semigroup satisfies and is uniquely identified by the integral estimates 44, 45 with U substituted by Ū see Section 5 for the point ξ =. More precisely let u : [, T] D be Lipschitz continuous as a map with values in LR, R n and assume that ut, x satisfies the integral conditions 44, 45 with U substituted by Ū for the point ξ =. Then ut, coincides with a trajectory of the semigroup P. Observe that the same Theorem holds for the 2 2 isentropic system see 4
7 Section 5 ρ t + q x = a x ax q q t + q 2 ρ +p x 8 = a x q 2 ax ρ. In [6] 2 2 homogeneous conservation laws at a junction are considered for given admissible junction conditions. The situation of a junctions with only two pipes with different cross sections can be modeled by our limit procedure or as in [6] with a suitable junction condition. If we define the function Ψ which describes the junction conditions as Ψ ρ, q, ρ 2, q 2 = ρ 2, q 2 Φa + a, ρ, q 9 then it fulfills the determinant condition in [6, Proposition 2.2] since it satisfies Lemma 3. Here Φa, u is the solution of the ordinary differential equation 62 in Section 5. With these junction conditions one can show that the semigroup obtained in [6] satisfies the same integral estimate see the following proposition as our limit semigroup hence they coincide. Proposition. The semigroup defined in [6] with the junction condition given by 5 satisfies the integral estimates 44, 45 with U substituted by Ū for the point ξ =. The proof is postponed to Section 5. Remark 2. Note that Proposition justifies the coupling condition 9 as well as the condition used in [9] to study the Riemann problem for the gas flow through a nozzle. 3 Existence of BV entropy solutions Throughout the next two sections, we follow the structure of []. We recall some definitions and notations in there, and also the results which do not depend on the L boundedness of the source term. We will prove only the results which in [] do depend on the L bound using our weaker hypotheses. 3. The non homogeneous Riemann-Solver Consider the stationary equations associated to, namely the system of ordinary differential equations: For any x o R, v Ω, consider the initial data fvx x = gx, vx. vx o = v. 5
8 As in [], we introduce a suitable approximation of the solutions to,. Thanks to 4, the map u fu is invertible inside some neighborhood of the origin; in this neighborhood, for small h >, we can define ] Φ h x o, u =f [fu + g x o + s, u ds. 2 This map gives an approximation of the flow of in the sense that f Φ h x o, u fu = g x o + s, u ds. 3 Throughout the paper we will use the Landau notation Ø to indicate any function whose absolute value remains uniformly bounded, the bound depending only on f and M L. Lemma. The function Φ h x o, u defined in 2 satisfies the following uniform with respect to x o R and to u in a suitable neighborhood of the origin estimates. Φ h x o, C 2 Ø, lim Id D uφ h x o, u = h lim sup Φ h x o, u u =, h x o R 4 Proof. The Lipschitz continuity of f and 3 imply Φ h x o, u u = Φh x o, u f fu Ø gx o + s, u ds Ø ωx o + s ds Ø ε h h. Next we compute D u Φ h x o, u =Df [fu + Dfu + gx o + s, u ds ] D u gx o + s, u ds 6
9 which together to the identity u = f fu implies D u Φ h x o, u Id = Du Φ h x o, u Df fu [ ] h Df fu + gx o + s, u ds Df fu Dfu + + Df fu Ø ε h h. D u gx o + s, u ds D u gx o + s, u ds Finally, denoting with D i the partial derivative with respect to the i component of the state vector and by Φ h,l the l component of the vector Φ h, we derive D i D j Φ h,l x o, u = h D k D k f l fu + gx o + s, u ds k,k + k D i f k u + D j f k u + D k f l fu + D j D i f k u + D i g k x o + s, u ds D j g k x o + s, u ds gx o + s, u ds D j D i g k x o + s, u ds so that D 2 Φ h x o, u Ø + Mx o + s ds Ø + M L Ø. For any x o R we consider the system, endowed with a Riemann initial datum: u, x = { ul if x < x o u r if x > x o. 5 If the two states u l, u r are sufficiently close, let Ψ be the unique entropic homogeneous Riemann solver given by the map u r = Ψßu l = ψ n σ n... ψ σ u l, 7
10 where ß = σ,..., σ n denotes the signed wave strengths vector in R n, []. Here ψ j, j =,...,n is the shock rarefaction curve of the j th family, parametrized as in [3] and related to the homogeneous system of conservation laws u t + fu x =. 6 Observe that, due to 4, all the simple waves appearing in the solution of 6, 5 propagate with non-zero speed. To take into account the effects of the source term, we consider a stationary discontinuity across the line x = x o, that is, a wave whose speed is equal to, the so called zero-wave. Now, given h >, we say that the particular Riemann solution: ut, x = { ul if x < x o u r if x > x o. t 7 is admissible if and only if u r = Φ h x o, u l, where Φ h is the map defined in 2. Roughly speaking, we require u l, u r to be approximately connected by a solution of the stationary equations. Definition. Given h > suitably small, x o R, we say that ut, x is a h Riemann solver for, 4, 5, if the following conditions hold a there exist two states u, u + which satisfy u + = Φ h x o, u ; b on the set {t, x < x o }, ut, x coincides with the solution to the homogeneous Riemann Problem 6 with initial values u l, u and, on the set {t, x > x o }, with the solution to the homogeneous Riemann Problem with initial values u +, u r ; c the Riemann Problem between u l and u is solved only by waves with negative speed i.e. of the families,..., p; d the Riemann Problem between u + and u r is solved only by waves with positive speed i.e. of the families p +,...,n. Lemma 2. Let x o R and u, u, u 2 be three states in a suitable neighborhood of the origin. For h suitably small, one has Φ h x o, u u = Ø Φ h x o, u 2 Φ h x o, u u 2 u = Ø u 2 u ωx o + s ds, 8 xo+h x o ωsds.9 Lemma 3. For any M > there exist δ, h >, depending only on M and the homogeneous system 6, such that the following holds. For all maps φ C 2 R n, R n satisfying φ C 2 M, φu u h, I Dφu h 8
11 u u - + u l u r Figure 2: Wave structure in an h Riemann solver. and for all u l B, δ, u r Bφ, δ there exist n + states w,...,w n+ and n wave sizes σ,..., σ n, depending smoothly on u l, u r, such that with previous notations: i w = u l, w n+ = u r ; ii w i = Ψ i σ i w i, i =,...,p; iii w p+ = φw p ; iv w i+ = Ψ i σ i w i, i = p +,...,n. The next lemma establishes existence and uniqueness for the h Riemann solvers see Fig.2. Lemma 4. There exist δ, h > such that the following holds: for any x o R, h [, h ], u l, u r B, δ, there exists a unique h Riemann solver in the sense of Definition. Proof. By Lemma if h > is chosen sufficiently small then for any h [, h ], x o R the map u Φ h x o, u meets the hypotheses of Lemma 3. Finally taking h eventually smaller we can obtain that there exists δ > such that B, δ B, δ BΦ hx o,, δ, for any h [, h ]. In the sequel, E stands for the implicit function given by Lemmas 3 and 4: ß =Eh, u l, u r ; x o, which plays the role of a wave size vector. We recall that, by Lemma 3, E is a C 2 function with respect to the variables u l, u r and its C 2 norm is bounded by a constant independent of h and x o. In contrast with the homogeneous case, the wave size ß in the h Riemann solver is not equivalent to the jump size u l u r ; an additional term appears coming from the Dirac source term see the special case u l = u r. 9
12 Lemma 5. Let δ, h be the constants in Lemma 4. For u l, u r B, δ, h [, h ], set ß = E[h, u l, u r ; x o ]. Then it holds: u l u r = Ø ß + ß = Ø u l u r + ωx o + s ds, ωx o + s ds Existence of a Lipschitz semigroup of BV entropy solutions Note that as shown in [] we can identify the sizes of the zero waves with the quantity σ = ωjh + sds. 2 With this definition all the Glimm interaction estimates continue to hold with constants that depend only on f and on M L, therefore all the wave front tracking algorithm can be carried out obtaining the existence of ε, h-approximate solutions as defined below. Definition 2. Given ǫ, h >, we say that a continuous map u ǫ,h : [, + L loc R, R n is an ǫ, h approximate solution of 2 if the following holds: As a function of two variables, u ǫ,h is piecewise constant with discontinuities occurring along finitely many straight lines in the x, t plane. Only finitely many wave-front interactions occur, each involving exactly two wave-fronts, and jumps can be of four types: shocks or contact discontinuities, rarefaction waves, non-physical waves and zero-waves: J = S R N P Z. Along each shock or contact discontinuity x α = x α t, α S, the values of u = u ǫ,h t, x α and u + = u ǫ,h t, x α + are related by u + = ψ kα σ α u for some k α {,..., n} and some wave-strength σ α. If the k α th family is genuinely nonlinear, then the Lax entropy admissibility condition σ α < also holds. Moreover, one has ẋ α λ kα u +, u ǫ where λ kα u +, u is the speed of the shock front or contact discontinuity prescribed by the classical Rankine-Hugoniot conditions. Along each rarefaction front x α = x α t, α R, one has u + = ψ kα σ α u, < σ α ǫ for some genuinely nonlinear family k α. Moreover, we have: ẋ α λ kα u + ǫ.
13 All non-physical fronts x = x α t, α N P travel at the same speed ẋ α = ˆλ > sup u,i λ i u. Their total strength remains uniformly small, namely: u ǫ,h t, x α + u ǫ,h t, x α ǫ, t >. α N P The zero-waves are located at every point x = jh, j hǫ, hǫ Z. Along a zero-wave located at x α = j α h, α Z, the values u = u ǫ,h t, x α and u + = u ǫ,h t, x α + satisfy u + = Φ h x α, u for all t > except at the interaction points. The total variation in space Tot.Var.u ǫ,h t, is uniformly bounded for all t. The total variation in time Tot.Var. { u ǫ,h, x; [, + } is uniformly bounded for x jh, j Z. Finally, we require that u ǫ,h,. u o LR ǫ. Keeping h > fixed, we are about to let first ǫ tend to zero. Hence we shall drop the superscript h for notational clarity. Theorem 3. Let u ǫ be a family of ǫ, h approximate solutions of 2. There exists a subsequence u ǫi converging as i + in L loc, + R to a function u which satisfies for any ϕ C c, + R: [uϕ t + fuϕ x ] dxdt R h + ϕt, jh g [jh + s, ut, jh ]ds dt =. 22 j Z Moreover Tot.Var. ut, is uniformly bounded and u satisfies the Lipschitz property ut, x ut, x dx C t t, t, t ; 23 R Now we are in position to prove [, Theorem 4] with our weaker hypotheses. As in [] we can apply Helly s compactness theorem to get a subsequence u hi converging to some function u in L loc whose total variation in space is uniformly bounded for all t. Moreover, working as in [2, Proposition 5.], one can prove that u hi t, converges in L to ut,, for all t. Theorem 4. Let u hi be a subsequence of solutions of equation 22 with uniformly bounded total variation converging as i + in L to some function u. Then u is a weak solution to the Cauchy problem 2. We omit the proofs of Theorem 3 and 4 since they are very similar to the proofs of [, Theorem 3 and 4]. We only observe that, in those proofs, the computations which rely on the L bound on the source term have to be substituted by the following estimates.
14 Concerning the proof of Theorem 3: g jh + s, u ε t, jh g jh + s, ut, jh ds gjh + s, C u ε t, jh ut, jh ds ε h u ε t, jh ut, jh. Concerning the proof of Theorem 4: and g jh + s, u h t, jh ds gjh + s, C ds ε h g jh + s, u h t, jh g jh + s, ut, jh + s ds gjh + s, C u h t, jh ut, jh + s ds ε h Tot.Var. { u h t,, [j h, j + h] } + j+h jh ωx u h t, x ut, x. We observe that all the computations done in [, Section 4] rely on the source g only through the amplitude of the zero waves and on the interaction estimates. Therefore the following two theorems still hold in the more general setting. Theorem 5. There exists δ > such that if ω LR is sufficiently small, then for any small h > there exist a non empty closed domain D h δ and a unique uniformly Lipschitz semigroup P h : [, + D h δ D h δ whose trajectories ut,. = P h t u o solve 22 and are obtained as limit of any sequence of ǫ, h approximate solutions as ǫ tends to zero with fixed h. In particular the semigroups P h satisfy for any u o, v o D h δ, t, s P h u o = u o, P h t P h s u o = P h s+t u o, 24 [ ] P h t u o Ps h v o LR L u o v o LR + t s 25 for some L >, independent on h. Theorem 6. If ω LR is sufficiently small, there exist a constant L >, a non empty closed domain D of integrable functions with small total variation and a semigroup P : [, + D D with the following properties i P u = u, u D; P t+s u = P t P s u, u D, t, s. 2
15 ii P s u P t v LR L s t + u v LR, u D, t, s. iii for all u o D, the function ut, = P t u o is a weak entropy solution of system. iv for some δ > and all h > small enough D D h δ. v There exists a sequence of semigroups P hi such that P hi t u converges in L to P t u as i + for any u D. Remark 3. Looking at [, 4.6] and the proof of [, Theorem 7] one realizes that the invariant domains D h δ and D depend on the particular source term gx, u. On the other hand estimate [, 4.4] shows that all these domains contain all integrable functions with sufficiently small total variation. Since the bounds Ø in Lemma 5 depend only on f and on M L, also the constant C in [, 4.4] depends only on f and on M L. Therefore there exists δ > depending only on f and on M L such that D h δ and D contain all integrable functions ux with Tot.Var. {u} δ. 4 Uniqueness of BV entropy solutions The proof of uniqueness in [] strongly depends on the boundedness of the source, therefore we have to consider it in a more careful way. 4. Some preliminary results As in [] we shall make use of the following technical lemmas whose proofs can be found in [3]. Lemma 6. Let a, b a possibly unbounded open interval, and let ˆλ be an upper bound for all wave speeds. If ū, v D h δ then for all t and h >, one has b ˆλt a+ˆλt Pt h ū x Pt h v x dx L b a ūx vx dx. 26 Lemma 7. Given any interval I = [a, b], define the interval of determinacy I t = [a + ˆλt, b ˆλt], t < b a 2ˆλ. 27 For every Lipschitz continuous map w : [, T] D h δ and h > : wt P h t w 28 LIt { } L t lim inf ws + P h ws LIs+ ds. 3
16 Remark 4. Lemmas 6, 7 hold also substituting P h with the operator P. In this case we have obviously to substitute the domains D h δ with the domain D of Theorem 6. Let now u l, u r be two nearby states and λ < ˆλ; we consider the function { ul if x < λt + x vt, x = o 29 u r if x λt + x o. Lemma 8. Call wt, x the self-similar solution given by the standard homogeneous Riemann Solver with the Riemann data 5. i In the general case, one has t + vt, x wt, x dx = Ø u l u r ; 3 ii Assuming the additional relations u r = R i σu l and λ = λ i u r for some σ >, i =,...,n one has the sharper estimate t + vt, x wt, x dx = Øσ 2 ; 3 iii Let u Ω and call λ <... < λ n the eigenvalues of the matrix A = fu. If for some i it holds A u r u l = λ i u r u l and λ = λ i in 29 then one has t + vt, x wt, x dx = Ø u l u r u l u + u u r ; 32 We now prove the next result which is directly related to our h-riemann solver. Lemma 9. Call wt, x the self-similar solution given by the h Riemann Solver in x o with the Riemann data 5. i In the general case one has t + vt, x wt, x dx = Ø u l u r + ii Assuming the additional relation u r = u l + [ f] u with λ = in 29 one has the sharper estimate t + = Ø vt, x wt, x dx gx o + s, u ds ωx o + s ds + u l u ωx o + s ds ; 33 ωx o + s ds. 34 4
17 Proof. Estimate i is a direct consequence of Lemma 5. Let us prove now ii. Since λ = we derive t + vt, x wt, x dx ˆλt = u l wt, x dx + u r wt, x dx 35 t ˆλt t [ p ] n = Ø σ ι + σ ι = Ø ß. ι= ι=p+ This leads to ß = E [h, u l, u r ; x o ] E [h, u l, Φ h x o, u l ; x o ] = Ø u r Φ h x o, u l. To estimate this last term, we define by, u = f fu + y and compute for some y, y 2 : u l + [ f] u y by 2, u l Ø y u u l + Ø y y 2 + u l + [ f] u l y 2 b y 2, u l. The function zy 2 = u l +[ f] u l y 2 b y 2, u l satisfies z =, D y2 z =, hence we have the estimate [ u l + [ f] u y by 2, u l Ø y u u l + y y 2 + y 2 2]. If in this last expression we substitute then, we get y = gx o + s, u ds, y 2 = gx o + s, u l ds u r Φ h x o, u l = O ωx o + s ds + u l u ωx o + s ds which proves Characterization of the trajectories of P In this section we are about to give necessary and sufficient conditions for a function ut, D to coincide with a semigroup s trajectory. To this end, we prove the uniqueness of the semigroup P and the convergence of all the sequence of semigroups P h towards P as h. 5
18 We begin by introducing some notations: given a BV function u = ux and a point ξ R, we denote by U u;ξ the solution of the homogeneous Riemann Problem 5 with data u l = lim x ξ ux, u r = lim x ξ+ ux, x o = ξ. 36 Moreover we define Uu;ξ as the solution of the linear hyperbolic Cauchy problem with constant coefficients w t + Ãw x = gx, w, x = ux, 37 with à = f uξ, gx = g x, uξ. We will need also the following approximations of Uu;ξ. Let v be a piecewise constant function. We will call w h the solution of the following Cauchy problem: w h t + Ãwh x = j Z δx jh gjh + s ds, w h, x = vx. Define u =uξ and let λ i = λ i u, r i = r i u, l i = l i u be respectively the i th eigenvalue, the i th right/left eigenvectors of the matrix Ã. As in [] w and w h have the following explicit representation wt, x = w h t, x = n x { l i, u x λ i t + λi i= n i= } l i, gx dx r i x λ it { l i, v x λ i t + λ i li, G h t, x } r i, 38 n where the function G h t, x = G h i t, xr i is defined by G h i t, x = j: jh x λ it,x i= j: jh x,x λ it l i, g jh + s ds if λ i > l i, g jh + s ds if λ i <. 39 Using 3 we can compute Gh i t, x x x λ it a+ˆλt a l i, gx dx = Ø ε h. 4 Hence, for any a, b R with a < b, we have the error estimate b ˆλt wt, x w h t, x [ ] b dx Ø ux vx dx + b a ε h. 4 6
19 From 38, 39, it is easy to see that w h t, x is piecewise constant with discontinuities occurring along finitely many lines on compact sets in the t, x plane for t. Only finitely many wave front interactions occur in a compact set, and jumps can be of two types: contact discontinuities or zero waves. The zero waves are located at the points jh, j Z and satisfy w h t, jh+ w h t, jh = [ f] u j+h jh g jh + sds. 42 Conversely a contact discontinuity of the i th family located at the point x α t satisfies ẋ α t = λ i u and w h t, x α t+ w h t, x α t = σr i u 43 for some σ R. Now, we can state the uniqeness result in our more general setting. Theorem 7. Let P : D [, + D be the semigroup of Theorem 6 and let ˆλ be an upper bound for all wave speeds. Then every trajectory ut, = P t u, u D, satisfies the following conditions at every τ. i For every ξ, one has lim ξ+ˆλ ξ ˆλ uτ +, x U dx uτ;ξ, x =. 44 ii There exists a constant C such that, for every a < ξ < b and < < b a, 2ˆλ one has b ˆλ uτ +, x U uτ;ξ, x dx 45 a+ˆλ [ b 2. C Tot.Var. {uτ; a, b} + ωx dx] Viceversa let u : [, T] D be Lipschitz continuous as a map with values in LR, R n and assume that the conditions i, ii hold at almost every time τ. Then ut, coincides with a trajectory of the semigroup P. Remark 5. The difference with respect to the result in [] is the presence of the integral in the right hand side of formula 45. If ω is in L, the integral can be bounded by Øb a and we recover the estimates in []. Note also that the quantity µ a, b = Tot.Var. {uτ; a, b} + a b a ωx dx is a uniformly bounded finite measure and this is what is needed for proving the sufficiency part of the above Theorem. 7
20 Proof. Part : Necessity Given a semigroup trajectory ut, = P t ū, ū D we now show that the conditions i, ii hold for every τ. As in [] we use the following notations. For fixed h,, ε > we define J t = Jt Jt o J t + with Jt = ξ 2 t + τ ˆλ, ξ t τˆλ ; [ ] J o t = J + t = ξ t τˆλ, ξ + t τˆλ ξ + t τˆλ, ξ + 2 t + τ ˆλ ; 46. Let U,ε uτ;ξ, x be the piecewise constant function obtained from U uτ;ξ, x dividing the centered rarefaction waves in equal parts and replacing them by rarefaction fans containing wave fronts whose strength is less than ε. Observe that: t + U,ε τ uτ;ξ, x U uτ;ξ, x dx = Øε. 47 Applying estimate 28 to the function U,ε uτ;ξ we obtain U,ε uτ;ξ, x P h U,ε uτ;ξ x dx 48 J τ+ τ+ U,ε uτ;ξ t τ + P hu,ε uτ;ξ t τ LJt+ L lim inf dt. The discontinuities of U,ε uτ;ξ do not cross the Dirac comb for almost all times t τ, τ +. Therefore we compute for such a time t: = J t+ U,ε uτ;ξ t τ +, x P h U,ε J t+ Jo t+ J+ t+ uτ;ξ t τ x dx 49 U,ε uτ;ξ t τ +, x P h U,ε uτ;ξ t τ x dx. Define W t the set of points in which U,ε uτ;ξt τ has a discontinuity while Z h is the set of points in which the zero waves are located. If is sufficiently small, the solutions of the Riemann problems arising at the discontinuities of U,ε uτ;ξ t τ do not interact, therefore Jt+ o = U,ε uτ;ξ t τ +, x P h U,ε + x Jt o Wt x Jt o Z h x+ˆλ U,ε x ˆλ uτ;ξ t τ x dx uτ;ξ t τ +, y P h U,ε uτ;ξ t τ y dy 8
21 Note that the shock are solved exactly both in U,ε uτ;ξ and in P h U,ε uτ;ξ therefore they make no contribution in the summation. To estimate the approximate rarefactions we use the estimate 3 hence x Jt o Wt Ø x+ˆλ U,ε x ˆλ x J o t W t rarefaction Øε uτ, ξ+ uτ, ξ uτ;ξ t τ +, x P h U,ε uτ;ξ t τ x dx { } σ 2 Øε Tot.Var. U,ε uτ;ξ t τ; J t 5 Concerning the zero waves, recall that t is chosen such that U,ε uτ;ξ is constant there, and P h is the exact solution of an h Riemann problem, hence we can apply 33 with u l = u r and obtain x J o t Z h Ø jh J o t x+ˆλ U,ε x ˆλ uτ;ξ t τ +, x P h U,ε uτ;ξ t τ x dx ωjh + s ds Ø J o t ωx dx + ε h 5 Finally using 5 and 5 we get in the end U,ε uτ;ξ t τ +, x P h U,ε uτ;ξ t τ x dx 52 Jt+ o { } = Ø ωx dx + ε h + ε. Moreover, following the same steps as before and using 3 and 33 with u l = u r we get U,ε J + uτ;ξ t τ +, x P h U,ε uτ;ξ t τ x dx 53 t+ { } = Ø ωx dx + ε h. Note that here there is no total variation of U,ε uτ;ξ since in J+ t it is constant. A similar estimate holds for the interval Jt+. Putting together 49, 52, 53, one has J t+ U,ε uτ;ξ t τ +, x P h U,ε uτ;ξ t τ x dx J t J + t = Ø ωx dx + ε h + ε. J τ 9
22 Hence, setting ṽ = U,ε uτ;ξ = U uτ;ξ by 48, we have U,ε uτ;ξ, x P h ṽ x dx = Ø ωx dx + ε h + ε. 54 J τ J τ+ Finally we take the sequence P hi converging to P. Using 26 we have P hi uτ P hi ṽ LJτ+ L uτ ṽ LJ τ 55 = L = ε, ξ+2ˆλ ξ 2ˆλ uτ, x ṽx dx where ε tends to zero as tends to zero due to the fact that uτ has right and left limit at any point: for any given ǫ > if is sufficiently small uτ, x ṽx = uτ, x uτ, ξ ǫ for x ξ 2ˆλ, ξ. Therefore by 47, 54, we derive: = ξ+ˆλ uτ +, x U dx uτ;ξ, x ξ ˆλ P uτ P hi uτ [ ] LR + ε + Ø ωx dx + ε hi. J τ The left hand side of the previous estimate does not depend on ε and h i, hence ξ+ˆλ ξ ˆλ uτ +, x U uτ;ξ, x dx = O ωx dx + ε. J τ Note that the intervals J τ depend on see 46. So taking the limit as in the previous estimate yields 44. To prove ii let > and a point τ, ξ be given together with an open interval a, b containing ξ. Fix ε > and choose a piecewise constant function v D satisfying vξ = uτ, ξ together with b a vx uτ, x dx ε, Tot.Var. { v; a, b} Tot.Var. {uτ; a, b} 56 Let now w h be defined by 38 u = vξ = uτ, ξ. From 4, 56 we have the estimate b ˆλ a+ˆλ Uuτ;ξ, x wh, x dx Ø ε + ε h b a. 57 2
23 Using 27, 28 we get b ˆλ L a+ˆλ τ+ τ w h, x P h wh x dx 58 w h t τ + P h w h t τ L Ĩ lim inf t+ dt where we have defined Ĩt+ = I t τ+. Let t τ, τ + be a time for which there is no interaction in w h ; in particular, discontinuities which travel with a non-zero velocity do not cross the Dirac comb this happens for almost all t. We observe that by the explicit formula 38: { } Tot.Var. w h t τ; Ĩt = Ø Tot.Var. { v; a, b} + w h t τ, x vξ = Ø Tot.Var. { v; a, b} + b a b a ωx dx + ε h 59 ωx dx + ε h. 6 As before for sufficiently small we can split homogeneous and zero waves w h t τ +, x P h wh t τ x dx 6 Ĩ t+ = x+ˆλ w h t τ +, x P h w h t τ x dx x Ĩt Wt + x Ĩt Z h x ˆλ The homogeneous waves in w h satisfy 43, with vξ in place of u, hence we can apply 32 which together with 59, 6 leads to x Ĩt Wt x+ˆλ x ˆλ w h t τ +, x P h w h t τ x dx Ø w h t τ, x Tot.Var. { v; a, b} + x Ĩt Wt Ø Tot.Var. Ø Tot.Var. { v; a, b} + { } w h t τ, Ĩt Tot.Var. { v; a, b} + b a ωx dx + ε h 2 b a b a ωx dx + ε h ωx dx + ε h where w h t τ, x denotes the jump of w h t τ at x. The zero waves in w h satisfy 42, hence we can apply 34 which together with 6 leads to 2
24 x Ĩt Z h Ø x Ĩt Z h x+ˆλ x ˆλ w h t τ +, x P h w h t τ x dx ωx + s ds Tot.Var. { v; a, b} + b a ωx dx + ε h b Ø ωx dx + ε h Tot.Var. { v; a, b} + ωx dx + ε h Ĩ t a b 2 Ø Tot.Var. { v; a, b} + ωx dx + ε h a Let now P hi be the subsequence converging to P. Since w h = v using 57, 58, 56, and the last estimates we get b ˆλ a+ˆλ uτ +, x Uuτ;ξ dx, x P uτ P hi uτ LR { ε + εhi b a +Ø uτ v LR + L + Tot.Var. { v; a, b} + b a ωx dx + ε hi 2 }. So for ε, h i we obtain the desired inequality. Part 2: Sufficiency By Remark 4 we can apply 28 to P and hence the proof for the homogeneous case presented in [3], which relies on the property recalled in Remark 5, can be followed exactly for our case, hence it will be not repeated here. Proof of Theorem It is now a direct consequence of Theorems 6 and 7. 5 Proofs related to Section 2 Consider the equation u t + fu x = a gu for some a BV. Equation 5 is comprised in this setting with the substitution a lna. For this kind of equations we consider the exact stationary solutions instead of approximated ones as in 2. Therfore call Φa, ū the solution of the following Cauchy problem: { d da ua = [D ufua] gua 62 u = ū 22
25 If a is sufficiently small, the map u Φa, u satisfies Lemma 3. We call a-riemann problem the Cauchy problem u t + fu x = a { gu a a, u, x =, u l if x < a +, u r if x > 63 its solution will be the function described in Definition using the map Φa + a, u instead of the Φ h in there. Observe that if a + = a the a-riemann solver coincides with the usual homogeneous Riemann solver. Definition 3. Given a function u BV and two states a, a +, we define Ū u t, x as the solution of the a-riemann solver 63 with u l = u and u r = u+. Proof of Theorem 2: Since a l L = a+ a, hypothesis P 2 is satisfied uniformly with respect to l, moreover the smallness of a + a ensures that the L norm of ω in P 3 is small. Therefore the hypotheses of Theorem are satisfied uniformly with respect to l. Let P l be the semigroup related with the smooth section a l. By Remark 3, if Tot.Var. {u} is sufficiently small, u belongs to the domain of P l for every l >. Since the total variation of Ptu l is uniformly bounded for a fixed initial data u, Helly s theorem guarantees that there is a converging subsequence Pt li u. By a diagonal argument one can show that there is a converging subsequence of semigroups converging to a limit semigroup P defined on an invariant domain see [, Proof of Theorem 7]. For the uniqueness we are left to prove the integral estimate 44 in the origin with U subsituted by Ū. Therefore we have to show that the quantity +ˆλ ˆλ uτ +, x Ū uτ, x dx 64 converges to zero as tends to zero. We will estimate 64 in several steps. First define v = Ū uτ, x and compute +ˆλ ˆλ P uτx P vx, x dx ǫ. 65 as in 55. Then we consider the approximating sequence P li corresponding to the source term a li and the semigroups P li,h which converge to P li in the sense of Theorem 6. Hence we have lim lim i h +ˆλ P li,h ˆλ vx P vx dx = 23
26 For notational convenience we skip the subscript i in l i. As in 47 we approximate rarefactions in Ū uτ introducing the function Ū,ε uτ. Then we define see Figure 3 Ū,ε,l,h uτ t τ, x = Ū,ε uτ t τ, x + l 2 for x < l/2 Ũx for l/2 x l/2 Ū,ε uτ t τ, x l 2 for x > l/2 where Ũx is piecewise constant with jumps in the points jh satisfying Ũjh+ = Ū,ε,l,h x = ˆλt x = l 2 + ˆλt t x uτ, = u l l 2 h l 2 uτ, + = u r Ũ Figure 3: Illustration of Ū,ε,l,h in the t, x plane Φjh, Ũjh. Furthermore Ũ l/2 = Ū,ε u;τt τ, and Φ is defined as in 2 using the source term gx, u = a l xgu. Observe that the jump between Ũl/2 and Ū,ε,l,h uτ t τ, l/2+ does not satisfy any jump condition, but as Ũx is an Euler approximation of the ordinary differential equation fu x = a l gu, this jump is of order ε h. Since Ū,ε uτ and Ū,ε,l,h uτ have uniformly bounded total variation we have the estimate +ˆλ ˆλ Ū,ε uτ, x Ū,ε,l,h uτ, x dx Ø l 24
27 the bound Ø not depending on h. We apply Lemma 7 on the remaining term +ˆλ P l,h ˆλ τ+ L τ lim inf vx Ū,ε,l,h uτ, x dx Ū,ε,l,h uτ t τ + P l,h Ū,ε,l,h uτ t τ LJt+ To estimate this last term we proceed as before. Observe that P l,h does not have zero waves outside the interval [ l 2 h, l 2 + h] since outside the interval [ l 2, l 2 ] the function a l is identically zero. If is small enough, the waves in P l,h Ū,ε,l,h uτ t τ do not interact, therefore the computation of the L norm in the previous integral, as before can be splitted in a summation on the points in which there are zero waves in P l,h or jumps in Ū,ε,l,h uτ t τ. Observe that the jumps of Ū,ε,l,h uτ t τ + in the interval l 2, + l 2, are defined exactly as the zero waves in P l,h so we have no contribution to the summation from this interval. Outside the interval [ l 2 h, l 2 + h], P h coincides with the homogeneous semigroup, hence we have only the second order contribution from the approximate rarefactions in Ū,ε,l,h uτ t τ as in 5. Furthermore we might have a zero wave in the interval [ l 2 h, l 2 ] and a discontinuity of Ū,ε,l,h uτ in the point x = l 2 of order ε h. Using 33 for the zero wave and 3 for the discontinuity since P h is equal to the homogeneous semigroup in x = l 2, we get lim inf Ū,ε,l,h uτ t τ + P l,h Ū,ε,l,h uτ t τ LJt+ Øε + ε h Which completes the proof if we let first ε tend to zero, then h tend to zero, then l tend to zero and finally tend to zero. As in the previous proof, the sufficiency part can be obtained following the proof for the homogeneous case presented in [3]. Proof of Proposition : Call S the semigroup defined in [6]. The estimates for this semigroup outside the origin are equal to the ones for the Standard Riemann Semigroup see [3]. Concerning the origin we first observe that the choice 9 implies that the solution to the Riemann problem in [6, Proposition 2.2] coincides with Ū uτ. We need to show that lim +ˆλ ˆλ uτ +, x Ū uτ, x dx =. 66 with Ū,ε uτ with ut, x = S t u o x. As before, we first approximate Ū uτ and uτ with Ū uτ = v then we apply Lemma 7 which holds also for the 25
28 semigroup S and compute +ˆλ L ˆλ τ+ τ dx S vx Ū,ε uτ, x lim inf Ū,ε uτ t τ + SŪ,ε uτ t τ LJ t+ The discontinuities of Ū,ε uτ are solved by S with exact shock or rarefaction for x and with the a Riemann solver in x = therefore the only difference between Ū,ε uτ t τ + and SŪ,ε uτ t τ are the rarefactions solved in an approximate way in the first function and in an exact way in the second. Recalling 3 we know that this error is of second order in the size of the rarefactions. To show that 66 holds, proceed as in 5. References [] Debora Amadori, Laurent Gosse, and Graziano Guerra. Global BV entropy solutions and uniqueness for hyperbolic systems of balance laws. Arch. Ration. Mech. Anal., 624: , 22. [2] Debora Amadori and Graziano Guerra. Uniqueness and continuous dependence for systems of balance laws with dissipation. Nonlinear Anal., 497, Ser. A: Theory Methods:987 4, 22. [3] Alberto Bressan. Hyperbolic systems of conservation laws, volume 2 of Oxford Lecture Series in Mathematics and its Applications. Oxford University Press, Oxford, 2. The one-dimensional Cauchy problem. [4] Rinaldo M. Colombo and Mauro Garavello. On the p-system at a junction. In Control methods in PDE-dynamical systems, volume 426 of Contemp. Math., pages Amer. Math. Soc., Providence, RI, 27. [5] Rinaldo M. Colombo, Graziano Guerra, Michael Herty, and Veronika Sachers. Modeling and optimal control of networks of pipes and canals. Preprint, 28. [6] Rinaldo M. Colombo, Michael Herty, and Veronika Sachers. On 2 2 conservation laws at a junction. SIAM J. Math. Anal., To appear. [7] Rinaldo M. Colombo and Francesca Marcellini. Smooth and discontinuous junctions in gas pipelines. In preparation, 28. [8] Rinaldo M. Colombo and Cristina Mauri. Euler system at a junction. Journal of Hyperbolic Differential Equations, To appear. [9] Paola Goatin and Philippe G. LeFloch. The Riemann problem for a class of resonant hyperbolic systems of balance laws. Ann. Inst. H. Poincaré Anal. Non Linéaire, 26:88 92, 24. [] P. D. Lax. Hyperbolic systems of conservation laws. II. Comm. Pure Appl. Math., : , 957. [] T.P. Liu. Quasilinear hyperbolic systems. Commun. Math. Phys., 68:4 72, 979. [2] A. I. Vol pert. Spaces bv and quasilinear equations. Mat. Sb. N.S., 73 5:255 32,
Hyperbolic 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 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 informationGlobal Existence of Large BV Solutions in a Model of Granular Flow
This article was downloaded by: [Pennsylvania State University] On: 08 February 2012, At: 09:55 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More 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 informationHurwitz generation of the universal covering of Alt(n)
Università di Milano Bicocca Quaderni di Matematica Hurwitz generation of the universal covering of Alt(n) M. A. Pellegrini, M. C. Tamburini Quaderno n. 1/2010 (arxiv:1001.5125v1) Stampato nel mese di
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationThe character table of a split extension of the Heisenberg group H 1 (q) by Sp(2, q), q odd
Università di Milano Bicocca Quaderni di Matematica The character table of a split extension of the Heisenberg group H (q by Sp(, q, q odd Marco Antonio Pellegrini Quaderno n. 7/8 (arxiv:math/85.48 Stampato
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationOn 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 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 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 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 informationInstability of Finite Difference Schemes for Hyperbolic Conservation Laws
Instability of Finite Difference Schemes for Hyperbolic Conservation Laws Alberto Bressan ( ), Paolo Baiti ( ) and Helge Kristian Jenssen ( ) ( ) Department of Mathematics, Penn State University, University
More informationOn finite time BV blow-up for the p-system
On finite time BV blow-up for the p-system Alberto Bressan ( ), Geng Chen ( ), and Qingtian Zhang ( ) (*) Department of Mathematics, Penn State University, (**) Department of Mathematics, University of
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationM. HERTY, CH. JÖRRES, AND B. PICCOLI
EXISTENCE OF SOLUTION TO SUPPLY CHAIN MODELS BASED ON PARTIAL DIFFERENTIAL EQUATION WITH DISCONTINUOUS FLUX FUNCTION M. HERTY, CH. JÖRRES, AND B. PICCOLI Abstract. We consider a recently [2] proposed model
More 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 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 informationBalance Laws as Quasidifferential Equations in Metric Spaces
Balance Laws as Quasidifferential Equations in Metric Spaces Rinaldo M. Colombo 1 Graziano Guerra 2 1 Department of Mathematics Brescia University 2 Department of Mathematics and Applications Milano -
More informationThe 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy
The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,
More informationOn a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws
On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws Stefano Bianchini and Alberto Bressan S.I.S.S.A., Via Beirut 4, Trieste 34014 Italy. E-mail addresses: bianchin@mis.mpg.de,
More informationShift Differentials of Maps in BV Spaces.
Shift Differentials of Maps in BV Spaces. Alberto Bressan and Marta Lewica SISSA Ref. 59/98/M (June, 998) Introduction Aim of this note is to provide a brief outline of the theory of shift-differentials,
More informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
More 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 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 informationExistence 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 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 informationNONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University. Email: wxs27@psu.edu November 5, 2015 Abstract We study two systems of conservation
More 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 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 information1. Introduction. This paper studies the initial value problem consisting of
ON 2 2 CONSERVATION LAWS AT A JUNCTION R.M. COLOMBO, M. HERTY, AND V. SACHERS Abstract. This paper deals with 2 2 conservation laws at a junction. For the Cauchy problem, existence, uniqueness and Lipschitz
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationEXISTENCE AND REGULARITY RESULTS FOR SOME NONLINEAR PARABOLIC EQUATIONS
EXISTECE AD REGULARITY RESULTS FOR SOME OLIEAR PARABOLIC EUATIOS Lucio BOCCARDO 1 Andrea DALL AGLIO 2 Thierry GALLOUËT3 Luigi ORSIA 1 Abstract We prove summability results for the solutions of nonlinear
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 informationSOURCE-DESTINATION FLOW ON A ROAD NETWORK
COMM. MATH. SCI. Vol. 3, No. 3, pp. 261 283 c 2005 International Press SOURCE-DESTINATION FLOW ON A ROAD NETWORK MAURO GARAVELLO AND BENEDETTO PICCOLI Abstract. We construct a model of traffic flow with
More informationDifferential Equations in Metric Spaces with Applications
Università di Milano Bicocca Quaderni di Matematica Differential Equations in Metric Spaces with Applications Rinaldo M. Colombo, Graziano Guerra Quaderno n. 18/27 arxiv:712.56v1 Stampato nel mese di dicembre
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
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 informationSTABILITY ESTIMATES FOR SCALAR CONSERVATION LAWS WITH MOVING FLUX CONSTRAINTS. Maria Laura Delle Monache. Paola Goatin
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX STABILITY ESTIMATES FO SCALA CONSEVATION LAWS WITH MOVING FLUX CONSTAINTS Maria Laura Delle Monache Department
More 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 informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University Abstract We study two systems of conservation laws for polymer flooding in secondary
More informationREGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS
SIAM J. MATH. ANAL. c 1988 Society for Industrial and Applied Mathematics Vol. 19, No. 4, pp. 1 XX, July 1988 003 REGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS BRADLEY J. LUCIER Abstract.
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 informationSingularity formation for compressible Euler equations
Singularity formation for compressible Euler equations Geng Chen Ronghua Pan Shengguo Zhu Abstract In this paper, for the p-system and full compressible Euler equations in one space dimension, we provide
More informationContractive Metrics for Nonsmooth Evolutions
Contractive Metrics for Nonsmooth Evolutions Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa 1682, USA bressan@mathpsuedu July 22, 212 Abstract Given an evolution
More informationSTRUCTURAL 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 informationNonclassical Shocks and the Cauchy Problem: General Conservation Laws
Contemporary Mathematics Contemporary Mathematics Volume 238, 1999 Volume 00, 1997 B 0-8218-1196-7-03536-2 Nonclassical Shocks and the Cauchy Problem: General Conservation Laws Paolo Baiti, Philippe G.
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 informationStructurally Stable Singularities for a Nonlinear Wave Equation
Structurally Stable Singularities for a Nonlinear Wave Equation Alberto Bressan, Tao Huang, and Fang Yu Department of Mathematics, Penn State University University Park, Pa. 1682, U.S.A. e-mails: bressan@math.psu.edu,
More informationPartial 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 information1. Introduction. Consider a strictly hyperbolic system of n conservation laws. 1 if i = j, 0 if i j.
SIAM J MATH ANAL Vol 36, No, pp 659 677 c 4 Society for Industrial and Applied Mathematics A SHARP DECAY ESTIMATE FOR POSITIVE NONLINEAR WAVES ALBERTO BRESSAN AND TONG YANG Abstract We consider a strictly
More informationSOLUTION OF AN INITIAL-VALUE PROBLEM FOR PARABOLIC EQUATIONS VIA MONOTONE OPERATOR METHODS
Electronic Journal of Differential Equations, Vol. 214 (214), No. 225, pp. 1 1. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SOLUTION OF AN INITIAL-VALUE
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 informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationDifferentiability with respect to initial data for a scalar conservation law
Differentiability with respect to initial data for a scalar conservation law François BOUCHUT François JAMES Abstract We linearize a scalar conservation law around an entropy initial datum. The resulting
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationError bounds for a Deterministic Version of the Glimm Scheme
Error bounds for a Deterministic Version of the Glimm Scheme Alberto Bressan and Andrea Marson S.I.S.S.A., Via Beirut 4, Trieste 34014 Italy. Abstract. Consider the hyperbolic system of conservation laws
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 informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationOn Asymptotic Variational Wave Equations
On Asymptotic Variational Wave Equations Alberto Bressan 1, Ping Zhang 2, and Yuxi Zheng 1 1 Department of Mathematics, Penn State University, PA 1682. E-mail: bressan@math.psu.edu; yzheng@math.psu.edu
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 informationEXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL
Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 05, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE OF NONTRIVIAL
More informationTraffic flow on networks: conservation laws models. Benedetto Piccoli I.A.C. C.N.R. Rome
Traffic flow on networks: conservation laws models Benedetto Piccoli I.A.C. C.N.R. Rome Conservation laws on networks u t + f(u) x=0 Dynamics at nodes? 1. The only conservation at nodes does not determine
More informationL 1 Stability for scalar balance laws. Control of the continuity equation with a non-local flow.
L 1 Stability for scalar balance laws. Control of the continuity equation with a non-local flow. Magali Mercier Institut Camille Jordan, Lyon Beijing, 16th June 2010 Pedestrian traffic We consider tu +
More informationUniversità di Milano Bicocca Quaderni di Matematica. A note on Schrödinger Newton systems with decaying electric potential. S.
Università di Milano Bicocca Quaderni di Matematica A note on Schrödinger Newton systems with decaying electric potential S. Secchi Quaderno n. 11/2009 arxiv:math.ap/0908.3768) Stampato nel mese di agosto
More informationRegularity and nonexistence results for anisotropic quasilinear elliptic equations in convex domains
Regularity and nonexistence results for anisotropic quasilinear elliptic equations in convex domains Ilaria FRAGALÀ Filippo GAZZOLA Dipartimento di Matematica del Politecnico - Piazza L. da Vinci - 20133
More informationScientiae Mathematicae Japonicae Online, Vol. 5, (2001), Ryo Ikehata Λ and Tokio Matsuyama y
Scientiae Mathematicae Japonicae Online, Vol. 5, (2), 7 26 7 L 2 -BEHAVIOUR OF SOLUTIONS TO THE LINEAR HEAT AND WAVE EQUATIONS IN EXTERIOR DOMAINS Ryo Ikehata Λ and Tokio Matsuyama y Received November
More informationRadon measure solutions for scalar. conservation laws. A. Terracina. A.Terracina La Sapienza, Università di Roma 06/09/2017
Radon measure A.Terracina La Sapienza, Università di Roma 06/09/2017 Collaboration Michiel Bertsch Flavia Smarrazzo Alberto Tesei Introduction Consider the following Cauchy problem { ut + ϕ(u) x = 0 in
More informationA maximum principle for optimally controlled systems of conservation laws
RENDICONTI del SEMINARIO MATEMATICO della UNIVERSITÀ DI PADOVA ALBERTO BRESSAN ANDREA MARSON A maximum principle for optimally controlled systems of conservation laws Rendiconti del Seminario Matematico
More informationConical Shock Waves for Isentropic Euler System
Conical Shock Waves for Isentropic Euler System Shuxing Chen Institute of Mathematical Research, Fudan University, Shanghai, China E-mail: sxchen@public8.sta.net.cn Dening Li Department of Mathematics,
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationAN OVERVIEW ON SOME RESULTS CONCERNING THE TRANSPORT EQUATION AND ITS APPLICATIONS TO CONSERVATION LAWS
AN OVERVIEW ON SOME RESULTS CONCERNING THE TRANSPORT EQUATION AND ITS APPLICATIONS TO CONSERVATION LAWS GIANLUCA CRIPPA AND LAURA V. SPINOLO Abstract. We provide an informal overview on the theory of transport
More informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationBlowup 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 informationThe Riemann problem for a class of resonant hyperbolic systems of balance laws
The Riemann problem for a class of resonant hyperbolic systems of balance laws Paola Goatin and Philippe G. LeFloch Abstract. We solve the Riemann problem for a class of resonant hyperbolic systems of
More information1 Basic Second-Order PDEs
Partial Differential Equations A. Visintin a.a. 2011-12 These pages are in progress. They contain: an abstract of the classes; notes on some (few) specific issues. These notes are far from providing a
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationHyperbolic Gradient Flow: Evolution of Graphs in R n+1
Hyperbolic Gradient Flow: Evolution of Graphs in R n+1 De-Xing Kong and Kefeng Liu Dedicated to Professor Yi-Bing Shen on the occasion of his 70th birthday Abstract In this paper we introduce a new geometric
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + FV εx, u = 0 is considered in R n. For small ε > 0 it is shown
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationASYMPTOTIC BEHAVIOR OF SOLUTIONS FOR PARABOLIC OPERATORS OF LERAY-LIONS TYPE AND MEASURE DATA
ASYMPTOTIC BEHAVIOR OF SOLUTIONS FOR PARABOLIC OPERATORS OF LERAY-LIONS TYPE AND MEASURE DATA FRANCESCO PETITTA Abstract. Let R N a bounded open set, N 2, and let p > 1; we study the asymptotic behavior
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More informationOptima and Equilibria for Traffic Flow on Networks with Backward Propagating Queues
Optima and Equilibria for Traffic Flow on Networks with Backward Propagating Queues Alberto Bressan and Khai T Nguyen Department of Mathematics, Penn State University University Park, PA 16802, USA e-mails:
More informationTwo dimensional exterior mixed problem for semilinear damped wave equations
J. Math. Anal. Appl. 31 (25) 366 377 www.elsevier.com/locate/jmaa Two dimensional exterior mixed problem for semilinear damped wave equations Ryo Ikehata 1 Department of Mathematics, Graduate School of
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 informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationSYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy
More informationExtremal Solutions of Differential Inclusions via Baire Category: a Dual Approach
Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Alberto Bressan Department of Mathematics, Penn State University University Park, Pa 1682, USA e-mail: bressan@mathpsuedu
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationA converse to the Ambrosetti-Prodi theorem
A converse to the Ambrosetti-Prodi theorem Marta Calanchi, Università di Milano with Carlos Tomei and André Zaccur (PUC-Rio, Brazil) Varese, RISM, September 2015 The Ambrosetti-Prodi Theorem 1/14 The Ambrosetti-Prodi
More informationGRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS
LE MATEMATICHE Vol. LI (1996) Fasc. II, pp. 335347 GRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS CARLO SBORDONE Dedicated to Professor Francesco Guglielmino on his 7th birthday W
More information