CONVERGENCE RATE OF APPROXIMATE SOLUTIONS TO CONSERVATION LAWS WITH INITIAL RAREFACTIONS
|
|
- Theodore Cameron
- 5 years ago
- Views:
Transcription
1 CONVERGENCE RATE OF APPROXIMATE SOLUTIONS TO CONSERVATION LAWS WITH INITIAL RAREFACTIONS HAIM NESSYAHU AND TAMIR TASSA Abstract. We address the question of local convergence rate of conservative Lip + -stable approximations, u (x, t, to the entropy solution, u(x, t, of a genuinely nonlinear conservation law. This question has been answered in the case of rarefaction free, i.e. Lip + -bounded, initial data. It has been shown that, by post-processing u, pointwise values of u and its derivatives may be recovered with an error as close to O( as desired, measures, in W,, the truncation error of the approximate solution u. In this paper we extend the previous results by including Lip + -unbounded initial data. Specifically, we show that for arbitrary L BV initial data, u and its derivatives may be recovered with an almost optimal, modulo a spurious log factor, error of O( ln. Our analysis relies on obtaining new Lip + -stability estimates for the speed, a(u, rather than for u itself. This enables us to establish an O( ln convergence rate in W,, which in turn, implies the above mentioned local convergence rate. We demonstrate our analysis for four types of approximate solutions: viscous parabolic regularizations, pseudo-viscosity approximations, the regularized Chapman-Enskog expansion and spectral-viscosity methods. Our approach does not depend on the geometry of the characteristics of the solution and, therefore, applies equally to finite-difference approximations of the conservation law. Key words. viscosity methods Conservation laws, Lip + -stability, W, -consistency, error estimates, parabolic regularizations, spectral. Intoduction. We study the convergence rate of approximate solutions of the single convex conservation law (. [u(x, t] + t x [f(u(x, t] = 0, t > 0, f α > 0, with compactly supported (or periodical initial condition (.2 u(x, t = 0 = u 0 (x, u 0 L BV. Our main focus in this paper is the extension of previous convergence results by allowing possibly Lip + - unbounded initial conditions, (.3, Lip + denotes the usual Lip + -semi-norm u 0 (x Lip +, ( + w(x w(y w(x Lip + ess sup, ( + max(, 0. x y x y It is well known that the solution of (. is not uniquely determined by the initial condition (.2 in the class of weak solutions. The unique physically relevant weak solution is the one which may be realized as a small viscosity solution of the parabolic regularization (.4 t [u (x, t] + x [f(u (x, t] = 2 x 2 [Q(u (x, t], Q 0, 0. We recall that these admissible, so-called entropy solutions, are characterized by their Lip + -stability [8]: (.5 a(u(, t Lip + a(u(, 0 Lip + + t a(u(, 0 Lip +, a( = f (. We therefore seek the convergence rate of conservative approximations to (., u (x, tdx = u 0 (xdx, t 0, x x School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv ISRAEL This research was supported in part by the Basic Research Foundation, Israel Academy of Sciences and Humanities.
2 which mimic this one sided Lipschitz stability of the exact entropy solution. This leads to Definition.. A family {u (x, t} >0 of approximate solutions of the conservation law (. is strongly Lip + -stable if (.6 a(u (, t Lip + a(u (, 0 Lip + + t a(u (, 0 Lip +, > 0. Our first convergence rate result is the content of the following theorem: Theorem.2. Let {u (x, t} >0 be a family of conservative and strongly Lip + -stable approximations to the entropy solution of (. (.2, u(x, t. Then, (i If u 0 Lip + <, the following error estimate holds (K and K 2 denote constants which depend on T : (.7a u (, T u(, T W, K u (, 0 u(, 0 W, + K 2 u t + f(u x L ([0,T ],W, (R x ; (ii If u 0 Lip + holds: (.7b = and the approximate solutions are also L -stable, the following error estimate ( u (, T u(, T W, O u (, 0 u(, 0 W,+ +O( u (, 0 BV + O( u(, 0 BV + O ( ln u t + f(u x L ([0,T ],W, (R x. Remarks.. An approximate solution operator, S (t, is considered L -stable, if for any two initial conditions, u 0 and v 0, (.8 S (tu 0 S (tv 0 L (R x Const t S (0u 0 S (0v 0 L (R x, t > The norm w(x, t W, is defined, when w(x, tdx = 0, as follows: R x w(x, t W, = w(x, t W, (R x w(ξ, tdξ L(R x. 3. The use of stability with respect to the Lip + -semi-norm in order to establish uniqueness for the Cauchy problem (. (.2, goes back to Oleinik [2] (see also Theorem.8 later on. Stability, in a similar sense, with respect to that semi-norm, was also used in [2] in order to obtain the total variation boundedness and entropy consistency of some finite difference approximations to (. and, consequently, their convergence to the entropy solution. However, this analysis lacks convergence rate estimates. The first to have used Lip + -stability in order to quantify the convergence rate, was Tadmor [8]. He used the Lip + -stability of both the entropy solution and its parabolic regularization, (.4, in order to quantify the convergence rate of the regularization. The same ideas were also used in [0, ] in the context of finite difference approximations. These works employed the Lip + -stability of the approximation itself, u (x, t, namely, an estimate of the sort (.9 u u (, 0 Lip + (, t Lip + + βt u (, 0 Lip +, 0 β α, in order to obtain convergence rate in the case of Lip + -bounded initial data. In fact, in that case, our first W, -error estimate, (.7a, holds even if the family of approximate solutions is merely Lip + -bounded, (.0 u (, t Lip + Const t, > 0, 2
3 and does not satisfy the strong Lip + -stability requirement (.6. However, estimates such as (.9 or (.0 are not sufficient in the case of Lip + -unbounded initial data and a stronger Lip + -stability, (.6, of a(u (x, t is required. As a counter-example we mention the Roe scheme (consult []: When u 0 Lip + < this scheme remains Lip + -bounded, (.0, and converges to the exact entropy solution. However, it is not strongly Lip + -stable and, therefore, it fails to converge to the entropy solution in case of Lip + -unbounded initial data (as demonstrated by the steady state solution obtained by this scheme for u 0 (x = sgn(x. The strong Lip + -stability, (.6, is indeed one of the main ingredients in establishing convergence rate estimates when initial rarefactions are present. Unfortunately, many well-known approximations of (. fail to satisfy this restricted condition. However, these approximations are still Lip + -stable in a weaker sense than that of Definition.. This weaker Lip + -stability proves sufficient in order to establish the same convergence rates as in Theorem.2. Definition.3. Let {u (x, t} >0 be a family of approximate solutions of (. and let W (t a(u (, t Lip +. Then this family is -weakly Lip + -stable if there exists a constant M such that whenever the following estimates hold for every T > 0: (. e W (0 M, 0 W (tdt O ( ; (.2 0 W e (τdτ t dt O ( ln. Remarks.. Any strongly Lip + -stable family of approximate solutions is also -weakly Lip + -stable (for any value of the constant M. 2. We henceforth refer by Lip + -stability to either weak or strong Lip + -stability. This notion of Lip + - stability is stronger than (.9, in view of the monotonicity of a(. The following theorem asserts that the convergence rate estimates, given in Theorem.2 for strongly Lip + -stable approximations, hold also for -weakly Lip + -stable ones. Theorem.4. Let {u (x, t} >0 be a family of conservative and Lip + -stable approximations to the entropy solution of (., u(x, t. Then, (i If u 0 Lip + <, error estimate (.7a holds; (ii If u 0 Lip + = and the approximate solutions are also L -stable, error estimate (.7b holds. In order to have convergence, the stability of the family of approximate solutions is not sufficient. The second crucial ingredient is consistency. if Definition.5. The family {u (x, t} >0 of approximate solutions is W, -consistent with (. (.2 (.3 u (, 0 u 0 ( W, { Const if u0 Lip + < Const 2 ln if u 0 Lip + = and (.4 u t + f(u x L ([0,T ],W, (R x Const T. 3
4 In view of Theorem.4 and Definition.5, we may now conclude the following convergence rate estimates. Corollary.6. (W, -Error Estimates. If the family {u (x, t} >0 of approximate solutions is conservative, W, -consistent with (. (.2, L -stable and Lip + -stable, then for every T > 0 there exists a constant C T such that (.5a u (, T u(, T W, C T, (.5b = { if u0 Lip + < ln if u 0 Lip + =. Remarks.. Error estimate (.5 suggests that whenever initial rarefactions are present, the convergence rate in W, is nearly O(. The ln term, which somewhat slows the rate of convergence, is a consequence of the initial rarefaction (as we show later on. 2. Error estimate (.5 relates to that of Harabetian in [3]. He has shown an O( ln convergence rate in L for the viscous parabolic regularizations, (.4, when the exact entropy solution amounts to a pure rarefaction wave. The W, error estimate (.5 may be translated, along the lines of [8, 0], into various global, as well as local, error estimates which we summarize as follows: Corollary.7. (Global and Local Error Estimates. Let {u (x, t} >0 be a family of conservative, W, -consistent, L -stable and Lip + -stable approximate solutions of the conservation law (. (.2. Then the following error estimates hold ( is as in (.5b: (E u (, T u(, T W s,p C T sp 2p, s p, p ; (E2 (u (, T φ δ (x u(x, T Const x,t p p+2, δ p+2, ( Const x,t = Const T + p! p x p u(, T L (x δ,x+δ and φ δ (x = δ φ ( x δ is any unit mass C 0 (, -mollifier, satisfying x k φ(xdx = 0 for k =, 2,..., p ; (E3 u (x, T u(x, T Const x,t 3, ( Const x,t = Const T + u x (, T L (x 3,x+ 3. Remark. A similar treatment enables the recovery of the derivatives of u(x, t as well, consult [8, 4]. 4
5 We would like to point out two straightforward consequences of Theorem.2, interesting for their own sake. The first is a simple proof of the uniqueness of Lip + -stable solutions to (. (.2, Theorem.8, and the second is the W, -stability of entropy solutions of (., Theorem.9. Theorem.8. Weak solutions of the convex conservation law (. which are Lip + -stable, (.5, are uniquely determined by their initial value. Theorem.9. Let u and v denote two entropy solutions of the conservation law (., subject to the L BV initial data u 0 and v 0, respectively. Then (.6 v(, t u(, t W, Const t v 0 u 0 η W,, η = if u 0 and v 0 are Lip + -bounded and η = 2 otherwise. This paper is organized as follows: After 2 in which we prove our main results, Theorems.2.9, the rest of the paper is devoted to applications to various types of approximations. In 3 we deal with the family of viscous parabolic regularizations, (.4. We prove that these approximations are L -contractive, W, -consistent and Lip + -stable, in order to conclude that they converge to the exact entropy solution and satisfy the convergence rate estimates (E (E3. We further show that if the viscosity coefficient satisfies (.7 ( Q a 0, then the resulting approximation is even strongly Lip + -stable. The most natural choice (already presented by Von-Neumann, Lax and Wendroff, [6] of a monotone regularization coefficient, Q(u, which satisfies (.7 is Q( = a(. Hence, we refer to regularizations which satisfy condition (.7 as speed-like. In 4 we apply our analysis to pseudo-viscosity approximations. These approximations are parabolic regularizations with a gradient dependent viscosity, u t + f(u x = Q(u, p x, p := u x, Q p 0. Such approximations, with Q = Q(p, were introduced by von Neumann and Richtmeyer in [9] and discussed later in [7]. We derive conditions on the pseudo-viscosity coefficient, Q, under which the resulting approximation is Lip + -stable and W, -consistent and, consequently, satisfies error estimates (E (E3. In 5 we discuss the regularized Chapman-Enskog expansion for hydrodynamics (consult [4, 7]. We focus our attention on Burgers equation and demonstrate our analysis in this case. Finally, in 6, we show how the Spectral Viscosity (SV method (consult [8, 9, 20] fits into our framework as well. In the course of the analysis performed there, we introduce an extension argument which removes the need for an a-priori L -bound. This argument may also be used for other approximate solutions of (. for which an a-priori L -bound is not known in advance. We close the Introduction by referring to the applicability of our framework to finite difference schemes, {v x } x>0. It is shown in [0, ] that finite difference schemes in viscosity form are conservative, BV - bounded and W, -consistent with (. (.2. Hence, so that our convergence rate results will apply to these schemes, all that remains to show is that they satisfy our strict notion of Lip + -stability, (.6 or (. (.2. However, the best Lip + -stability estimates which have been established for finite difference schemes are of the form (.9. Since we have not been able, so far, to sharpen those estimates, we do not include a treatment of these approximations in the present paper. 5
6 2. Proof of main results. We begin this section by proving our basic convergence rate estimates, as stated in Theorems.2 and.4 in the Introduction. Since Theorem.2 deals with strongly Lip + -stable approximations, which are, as noted before, weakly Lip + -stable as well, it suffices to prove Theorem.4. Proof (of Theorem.4. We deal with conservative approximations to (. which take the following form (2. t [u (x, t] + x [f(u (x, t] = r (x, t, t > 0, 0, r (x, t is the truncation error of the approximation, and we need to estimate, in W,, the error e (x, t u (x, t u(x, t. Step. We first assume that both the exact entropy solution, u(x, t, and its approximation, u (x, t, have a Lip + -bounded initial data, i.e., (2.2 L + 0 = max{ a(u(, 0 Lip +, a(u (, 0 Lip +} <. (2.3 Subtracting (. from (2. we arrive at the equation which governs the error e (x, t, t [e (x, t] + x [ā (x, te (x, t] = r (x, t, t > 0, ā (x, t = Note that the monotonicity of a( implies that 0 a (ξu (x, t + ( ξu(x, t dξ. (2.4 min{a(u, a(u } ā (x, t max{a(u, a(u }. Integration of (2.3 with respect to x yields (2.5 t [E (x, t] + ā (x, t x [E (x, t] = R (x, t, t > 0, E (x, t = x e (ξ, tdξ, R (x, t = x r (ξ, tdξ. Integration of (2.5 over R against sgn(e and rearranging, yield that ( d dt E (, t L ā (x, t (2.6 x E (x, t dx + R (, t L. x The main effort henceforth is concentrated on upper bounding the integral on the right hand side of (2.6. To this end we suggest to divide the real line into intervals, R = n I n (t, I n (t = [x n (t, x n+ (t, in such a way that neither sgn(e nor sgn(e change within the interior of these intervals (the implicit assumption of piecewise smoothness of the solution, as in [5], may be removed by considering a further vanishing parabolic regularization which is omitted. We use this division to define the following function: a(u(x, t if x I n (t and E (x, t 0 In(t (2.7 â (x, t =. a(u (x, t if x I n (t and E (x, t 0 6 In(t
7 We now claim (and prove later on that ( ā (x, t (2.8 x E (x, t dx x Integration by parts of the right hand side of (2.8 yields ( ā (x, t (2.9 x E (x, t dx x x ( â (x, t x E (x, t dx. x x [â (x, t] E (x, t dx. The following inequality (whose proof is postponed provides us an upper bound for the integral on the right hand side of (2.9: (2.0a x [â (x, t] E (x, t dx L (t E (, t L, (2.0b x { L L + (t = max 0 + tl + 0 Inserting (2.9 and (2.0a into (2.6, we arrive at the inequality }, W (t = a(u (, t Lip +. (2. which implies that d dt e (, t W, L (t e (, t W, + r (, t W,, (2.2 e (, T W, e 0 L (tdt e L (, 0 W, + e (τdτ t r (, t W,dt. 0 Since, by the definition of L + 0.3, that (2.3 and in (2.2, W (0 L + 0, we conclude, in view of Lip+ -stability (see Definition e W (tdt 0 Const, Const L+ 0 M, (2.4 0 e t W (τdτ dt Const 2, Const 2 lnm lnl + 0. Using (2.3, (2.4 and (2.0b in (2.2, proves the desired error estimate (.7a. Finally, in order to conclude Step, we return to justify (2.8 and (2.0: First, we prove (2.8 by showing that the inequality holds in each interval I n (t, i.e, ( ā (x, t ( x E (x, t dx â (x, t (2.5 x E (x, t dx. I n (t I n (t Suppose that E (, t 0 in I n (t. Then by definition (2.7, (2.6 â (x, t = a(u(x, t x I n (t. There are two possibilities to consider. If e (x, t 0 in I n (t then by (2.4 (2.7 ā (x, t a(u(x, t, x E (x, t = sgn(e (x, t e (x, t 0 x I n (t. 7
8 Therefore, (2.5 follows in this case by (2.6 and (2.7. If, on the other hand, e (x, t 0 in I n (t, then (2.8 ā (x, t a(u(x, t, x E (x, t 0 x I n (t and (2.5 follows in this case as well. The case E (, t In (t 0 is treated similarly. This concludes the proof of (2.8. Next, we prove inequality (2.0. In view of definitions (2.7 and (2.0b, we conclude, using the Lip + - stability of the exact solution, a(u(, t Lip + L+ 0 + tl + 0, that x [â (x, t] satisfies the following inequality in the sense of distributions: (2.9 x [â (x, t] L (t + [â (x n (t + 0, t â (x n (t 0, t]δ(x x n (t, the sum being taken over all division points x n (t â (, t experiences a jump discontinuity, namely sgn(e (, t changes. But, E (, t being a continuous primitive function vanishes at these points. Hence, integration of (2.9 against E (x, t proves (2.0a and completes Step. Step 2. Now we turn to the case of initial rarefactions and prove error estimate (.7b. To this end we introduce the function ψ δ ( = δ ψ( δ, δ > 0, which is the dilated mollifier of { x (2.20 ψ(x = 2 0 x >. 2 Clearly (2.2 ψ δ w w L O(δ w BV, and (2.22 ( ψ δ w Lip + O δ δ 0. With this in mind we return to the conservation law (. and its approximate solution (2. and define a new pair of solutions, u δ and u δ, corresponding to the mollified initial data: (2.23 (2.24 t [u δ(x, t] + x [f(u δ(x, t] = 0, u δ (, 0 = ψ δ u(, 0 ; t [u δ(x, t] + x [f(u δ(x, t] = rδ(x, t, u δ(, 0 = ψ δ u (, 0. We are now able to estimate the W, -error in (.7b by decomposing it as follows: (2.25 u (, T u(, T W, u (, T u δ(, T W, + u δ(, T u δ (, T W, + u δ (, T u(, T W,. Since for compactly supported functions, w W, supp{w} w L, we may bound the first term on the right hand side of (2.25, using (.8, (2.24 and (2.2, as follows (Ω T denotes the compact support at t = T : (2.26 u (, T u δ(, T W, Ω T u (, T u δ(, T L Note that in case u (, T is not compactly supported, the exponential decay which characterizes the tail of various viscosity-like approximations will suffice for our estimates. 8
9 Ω T C T u (, 0 u δ(, 0 L Ω T C T O(δ u (, 0 BV = O(δ u (, 0 BV. Similarly, the last term on the right hand side of (2.25, may be bounded by (2.27 u δ (, T u(, T W, O(δ u(, 0 BV. Hence, it remains only to deal with the term u δ (, T u δ(, T W,. This requires δ to be appropriately chosen so that (2.28 W δ (0 M, W δ (t = a(u δ(, t Lip + and, consequently, the Lip + -stability estimates x (. (.2 hold. If D denotes the largest positive jump in u (, 0 then the choice δ = 2D max[a (u (, 0]/M will do for (2.28. By doing so, we may conclude the -weak Lip + -stability estimates, (. (.2, for W δ (t: ( e W 0 δ (tdt O ; 0 W e t δ (τdτ dt O( ln. These estimates, together with error estimate (2.2 for e δ = u δ u δ, imply that (2.29 u δ(, T u δ (, T W, O ( u δ(, 0 u δ (, 0 W, + O( ln rδ L ([0,T ],W, (R x. Since ψ δ w W, w W,, estimate (2.29 implies that (2.30 u δ(, T u δ (, T W, O ( u (, 0 u(, 0 W, + O( ln r L ([0,T ],W, (R x. Therefore, since δ = O(, (.7b follows from (2.25, (2.26, (2.27 and (2.30 and the proof is thus concluded. Remark. Note that if the approximate solution smoothens the initial data so that ( u (, 0 Lip + O, e.g. the SV-method, there is no need to mollify the initial data of the approximation, as we did in (2.24. Hence, in this case, the error term (2.26 does not exist and, therefore, error estimate (.7b holds even if the approximate solution is not L -stable. We close this section with the proof of Theorems.8 and.9. Proof (of Theorem.8. Let u be the entropy solution of (. (.2 and v be another weak solution of (. (.2 which is also Lip + -stable x in the sense of (.5. Setting u = v, > 0, we have Hence, error estimate (.7b implies that u (, 0 u(, 0 = 0 and u t + f(u x = 0 > 0. v(, T u(, T W, = u (, T u(, T W, O( u 0 BV > 0. Letting 0, we conclude that u = v. 9
10 Proof (of Theorem.9. We set u = v for all > 0 and use error estimates (.7a and (.7b, given in Theorem (.. Since u is an exact entropy solution of (., the truncation error term on the right hand side of both estimates vanishes. In case that both u 0 and v 0 are Lip + -bounded, estimate (.7a holds and (.6 follows with Const t =K and η =. If either of the initial conditions is Lip + -unbounded, estimate (.7b holds and we conclude that v(, t u(, t W, O ( for all > 0. Taking = v 0 u 0 2 W,, proves (.6 with η = 2. v 0 u 0 W, + O( ( v 0 BV + u 0 BV 3. Viscous parabolic regularizations. We consider here viscous parabolic regularizations to (. of the form (.4. These regularizations are: Conservative; L -bounded, u (, t L u 0 L ; L -contractive and, therefore, thanks to translation invariance, BV -bounded (see Theorem 4., later on, for a proof of L -contraction in a more general setting; W, -consistent in the sense of Definition.5, since u (, 0 = u 0 ( and u t + f(u x W, = Q(u x L max Q (u u (, t BV O( ; u u 0 L Lip + -stable (Theorem 3.. In view of the above, error estimates (E (E3, x given in Corollary.7, apply to this family of approximate solutions. We are therefore left only with the task of proving Lip + -stability; this is done in the following theorem and lemma. (3. Theorem 3.. The (possibly degenerate parabolic regularization of (., is strongly Lip + -stable if t [u (x, t] + x [f(u (x, t] = 2 x 2 [Q(u (x, t], Q 0, 0, (3.2 ( Q a 0, and -weakly Lip + -stable otherwise. Proof. Let us first assume that Q is strictly positive so that the solution u is smooth. Multiplying (3. by a (u (x, t we get (3.3 By denoting (3.4 t [a(u ] + a(u x [a(u ] = a (u 2 x 2 [Q(u ]. w = w (x, t = a(u x the right hand side of (3.3 may be rewritten as follows: [ (3.5 a (u 2 x 2 [Q(u ] = Q (u w x + 0 = a (u u x, ( ] Q (u a (u (w 2.
11 Differentiation of (3.3 with respect to x and using identity (3.5 yields (3.6 w t + (w 2 + a(u w x = [ Q (u 2 w x 2 + Q (u a (u w w x + 2 ( Q (u a (u w w x + ( Q (u a (u ] (w 3 a (u Since u is smooth and compactly supported, w (, t attains its maximal value, say in x = x(t, and. (3.7 w (x(t, t 0, w (x(t, t = 0, x 2 w (x(t, t 0. x2 Hence, denoting (3.8 W (t = w (x(t, t = a(u (, t Lip +, we conclude by (3.6, (3.7 and the positivity of a and Q, that (3.9 (3.0 dw dt K α + (W 2 K(W 3, max u u 0 L [ (Q ] + (u a (u. In view of Lemma 3.2 below, inequality (3.9 implies -weak Lip + -stability. In case that condition (3.2 holds, K = 0 and inequality (3.9 amounts to Ricatti s inequality which implies strong Lip + -stability. dw dt + (W 2 0, If Q 0, equation (3. is degenerate and, therefore, admits non-smooth solutions. This case may be treated, as in [2], by introducing a further regularization. We replace Q( by the strictly monotone regularization term Q δ ( = Q( + δa(. Note that with this choice of Q δ, the value of K, (3.0, does not change. Hence, the corresponding solution, u δ, satisfies inequality (3.9 and by letting δ 0, we obtain the same inequality for the limit solution. Remark. The most common choice of a regularization coefficient is Q(u = u. For this special choice of Q(u, the speed-like condition (3.2 reads ( a 0, consult [6]. Lemma 3.2. Let y (t denote the solution of (3. (3.2 and c satisfies (3.3 dy dt + (y 2 = K(y 3, K > 0, t > 0, y (t = 0 = c K 0 < c c c <, 0.
12 Then, for any T > 0, (3.4 ( e 0 y (tdt O and (3.5 0 y e (τdτ t dt O( ln. The proof of this Lemma is postponed to the Appendix. Note that Lemma 3.2, together with (3.8 and inequality (3.9, show that the approximate solutions u (x, t are -weakly Lip + -stable with any constant M < /K (see Definition Pseudo-viscosity approximations. One of the methods for the approximation of phenomena governed by hyperbolic conservation laws is considering parabolic regularizations with a gradient dependent viscosity. These so-called pseudo-viscosity approximations take the form (4. u t + f(u x = Q(u, p x, p := u x, 0, (4.2 (4.3 u (x, 0 = u 0 (x, Q p 0. Note that this class of parabolic regularizations is wider than the class of viscous parabolic approximations, (3.. First, we note that these conservative approximations satisfy the maximum principle and, therefore, the solution remains uniformly bounded by u 0 L. Next, the following theorem (whose proof is postponed to the Appendix asserts that the solution operator of the pseudo-viscosity approximation is L -contractive. Therefore, thanks to translation invariance, the solution u remains BV -bounded. Theorem 4.. (L -Contraction. Let u and v be two solutions of (4., (4.3. Then (4.4 u (, t v (, t L u (, 0 v (, 0 L, t > 0. Finally, we address the question of Lip + -stability. We show that under suitable assumptions on the pseudo-viscosity coefficient, Q(u, p, the solution of (4. is weakly Lip + -stable. Theorem 4.2. (Lip + -Stability. Let Ω denote the domain in R 2, Ω = [inf u 0, sup u 0 ] [0,. Assume that the following hold for all (u, p Ω (M and M 2 denote some constants: (A Q p (u, p, Q up (u, p M ; (A2 Q uu (u, p M 2 p ; 2
13 (A3 Q pp (u, p 0. Then the solution of (4. (4.3 is -weakly Lip + -stable. Proof. We first deal with the uniformly parabolic case, Q p δ > 0. Let us denote w (x, t = x [a(u (x, t], W (t = max x w (x, t = a(u (, t Lip +. In view of Lemma 3.2, it suffices to show that there exists a constant K > 0, such that (4.5 d dt W (t + (W (t 2 K(W (t 3, t > 0. Multiplying (4. by a (u and differentiating with respect to x, we find that w = w (x, t satisfies w t + w 2 w + aw x = [Q 2 w uu a + 2Q up a (w x + A w 2 (w x + A w Q pp a + + Q u w x + Q p (w xx + 2A ww x + A w3 a ], a = a(u and A = A(u = /a (u. Let (x(t, t be a positive local maximum of w. Then w > 0 in that point and, since a α > 0, (., also p = u x > 0 there. Furthermore, w x = 0 and w xx 0 in that point. Therefore, in view of (4.3 and assumptions (A (A3, the above inequality implies that w t + w 2 Kw 3 in (x(t, t, for some constant K which depends on M, M 2, α and the uniform bounds on A and A. Therefore, (4.5 holds and that concludes the proof for the non-degenerate case. In the degenerate case, we replace Q(u, p by Q δ (u, p = Q(u, p+δp so that the resulting pseudo-viscous approximation will be uniformly parabolic, Q δ /p δ > 0, and admit a smooth solution, u δ. Note that Q δ, δ 0, still satisfies conditions (A (A3 with constants, say, M + and M 2. Therefore, inequality (4.5, with K independent of δ, holds for u δ, δ 0, and consequently it holds for u as well. Remark. Theorem 4.2 implies, in particular, the (-weak Lip + -stability of viscous parabolic regularizations, (3., stated earlier in Theorem 3.. These regularizations are identified by viscosity coefficients of the form (4.6 Q(u, p = q(u p, q(u 0. Such coefficients satisfy assumptions (A (A3, provided that q( is sufficiently smooth. We therefore conclude, in light of Theorems 4. and 4.2, that Theorem.4 applies to approximation (4. under assumptions (4.3 and (A (A3. Hence, if in addition, approximation (4. is W, -consistent with (., i.e., or simply, u t + f(u x W, (R x O(, (4.7 Q (u, u x L(R x Const, Corollary.7 may be applied and error estimates (E (E3 hold. We propose below a condition on Q(u, p which guarantees W, -consistency, (4.7. 3
14 Proposition 4.3. If there exists a constant C > 0, such that (4.8 Q(u, p C p (u, p [inf u 0, sup u 0 ] R, then equation (4. is W, -consistent with (.. Proof. Condition (4.8 implies that Therefore, (4.7 holds and the proof is concluded. Q (u, u x L(R x C u x L = C u BV C u 0 BV. An example of a family of pseudo-viscosity coefficients which satisfy all the above requirements, i.e., (4.3, (A (A3 and (4.8, is the following: (4.9 Q(u, p = Q q(u,β (u, p = q(u [ ( + p β ] sgn(p, q(u 0, 0 < β. Note that by letting β go to zero we obtain Q 0, which corresponds to the hyperbolic conservation law, while the other extreme case, β =, coincides with the standard viscous parabolic coefficient, (4.6. A special class of pseudo-viscosity approximations, (4., Q = Q(p, (4.0 u t + f(u x = Q(p x, Q 0, 0, was introduced by von Neumann and Richtmeyer in [9]. In [7] it is shown, by means of compensated compactness, that under further assumptions on the pseudo-viscosity coefficient, there exists a subsequence of weak solutions of (4.0, subject to the initial data (4.2, which converges in L p loc to the corresponding entropy solution of (., provided that u 0 W 2,. One of the additional restrictions assumed on Q in [7] is that it acts only on shock-waves and does not smear out rarefactions. Namely, (4. Q (p = 0 p 0 and Q (p > 0 p < 0. Note that restriction (4. guarantees Lip + -stability, since conditions (4.3 and (A (A3 are clearly satisfied in this case. An example of a family of such pseudo-viscosity coefficients which lead to W, -consistent approximations (in view of Proposition 4.3 is (4.2 Q β (p = [ Q,β (u, p ] = ( p β, 0 < β, Q q(u,β (u, p being defined in (4.9. The choice which corresponds to β =, Q (p = p, activates the regular parabolic regularization only on shock-waves and leaves rarefactions untouched. 5. The regularized Chapman-Enskog expansion. In this section we discuss the regularized Chapman- Enskog expansion for hydrodynamics, proposed by Rosenau [4]. This so-called R-C-E approximation is studied in [7], it is shown that it shares many of the properties of the viscosity approximation, e.g. existence of traveling waves, monotonicity, L -contraction and Lip + -stability. Let us briefly recall the main results of [7]. The R-C-E approximation is presented in the form (5. u t + f(u x = [Q m u x] x, 0, (5.2 u (, 0 = u 0 (, with the choice of the unit-mass viscosity kernel Q(x = 2 e x, Q m (x = ( x (5.3 m Q m. 4
15 This is a pseudo-local dissipative approximation of the conservation law, the viscosity coefficient is being activated by means of convolution rather than multiplication (compare (5. to (3.. When m 0, Q m tends to the Dirac measure and the R-C-E approximation, (5., turns into the viscous parabolic approximation (5.4 u t + f(u x = u xx. Equation (5. may be rewritten in the equivalent form u t + f(u x = m 2 [u Q m u ]. The solution of (5.4 remains as smooth as its initial data ([7, Theorem 2.] and, therefore, if the initial data are discontinuous, weak solutions must be admitted. Since such solutions are not uniquely determined by the initial data, (5.4 is augmented with a Kružkov-like [4] entropy condition ([7, (4.], (5.5 t u c + x {sgn(u c[f(u f(c]} m 2 { u c sgn(u c[q (u c]}, for all c R. In particular, by substituting c = + sup u or c = sup u, we obtain from (5.5 that u is, respectively, a supersolution or a subsolution of (5.4 and therefore a weak solution. Hence, u is considered an entropy solution of (5.4 if it satisfies inequality (5.5 in the sense of distributions for all c R. The above inequality, (5.5, implies L -contraction, and hence BV -boundedness, Since, by (5., we also have W, -consistency. u (, t v (, t L u (, 0 v (, 0 L, u (, t BV u 0 BV. u t + f(u x W, Q m u x L Q m L u 0 BV O(, Finally, we deal with the question of Lip + -stability. Adding the smoothing viscosity term δu,δ xx to (5.4 and differentiating with respect to x, we get that w u,δ x satisfies w t + a (u,δ w 2 + a(u,δ w x = m 2 [w Q m w] + δw xx. Letting δ 0, we get that W (t max x w(x, t is governed by the Ricatti differential inequality (5.6 W (t + αw 2 (t 0. Restricting our attention to Burgers equation, a(u = u, the R-C-E approximation turns to be strongly Lip + -stable, in virtue of (5.6. Therefore, we conclude, in view of Theorem (., that the R-C-E approximation converges to the entropy solution of Burgers equation and error estimates (E (E3 hold. This extends, for Burgers equation, the convergence rate result of ([7, Corollary 5.2] which was restricted to u 0 C. 5
16 6. The spectral viscosity method. The method of Spectral Viscosity (SV is used for the approximate solution of (. in the 2π-periodic case. The family of approximate solutions, {u N (x, t}, constructed by this method, consists of trigonometric polynomials, u N (x, t = N k= N ûk(te ikx, which approximate the spectral projection of the exact entropy solution, P N u. This method takes the following conservative form (consult [20]: (6. t u N(x, t + x P N (f(u N (x, t = N x Q N (x, t x u N (x, t, (6.2 u N (, 0 = P N u 0 (. The right hand side of (6. consists of a vanishing viscosity amplitude of size N 0 and a viscosity kernel, Q N (x, t = N k =m ˆQk N (te ikx, activated only on high wave numbers, k m N >>. As in [20] we deal with real viscosity kernels with increasing Fourier coefficients, ˆQk ˆQ k, which satisfy (6.3 ( 2q mn k ˆQ k (t, k m N, q = Const >.5, and the spectral viscosity parameters, N and m N, behave asymptotically as (6.4 N N θ log N, m N N θ 2q, 0 < θ <. The use of the projection P N on the initial data is problematic since even if u 0 has a bounded variation, P N u 0 BV may grow as much as O(log N. This may be avoided by taking, for instance, the spectrally accurate de la Vallee Poussin projection, N k N u N (x, 0 = V P N u 0 σ k û 0k e ikx 2 (6.5, σ k =, k= N 2 2k N k > N 2 which satisfies u N (, 0 BV = V P N u 0 BV 3 u 0 BV. This, according to the total-variation boundedness of the SV method (consult [20, Corollary 2.3], implies that (6.6 u N (, t BV Const T, t [0, T ]. Hence, we hereafter assume (6.5. At the end of this section we will deal with the case described in (6.2 of employing the regular spectral projection on the initial datat. The SV method smoothens the initial data by smearing its discontinuities: Since u 0 (x = BV, it follows that û 0k = O ( k. Hence x [u N(, 0] = N N ikσ k û 0k e ikx k û 0k O(N, and therefore k= N k= N k= û0ke ikx (6.7 u N (, 0 Lip + O(N <. We now turn to deal with the Lip + -stability of this approximation. To this end we rewrite (6., as in [20, (2.4], in the following form, (6.8a t u N(x, t + x f(u N(x, t = 6
17 2 = N x 2 u N(x, t N x R N(x, t x u N(x, t + E N, E N = x (I P Nf(u N is a spectrally small error term and (6.8b R N (x, t = N k= N ˆR k (te ikx, ˆRk (t = { k < mn ˆQ k (t k m N. Multiplying (6.8a by a (u N and differentiating with respect to x yields for w = x a(u N: w t + a(u N w x + w 2 = N [w xx + 2A (u N ww x + A w 3 ] (6.9 (u N a (u N ( N [a (u N A(u N x R N ( x u N w + a 2 (u N x 2 R N ] x u N + +a (u N A(u N we N + a (u N x E N. Here, as in 4, A( = /a (. As before, we find that W (t = max x w(x, t is governed by (6.0 (6. (6.2 and (6.3 d dt W (t + (W (t2 N K(W (t 3 + β N W (t + γ N, K max u u N L ( A + (u a, (u [ β N = M N x R N ] x u N L + E N L γ N = M 2 [ N 2 x 2 R N x u N L + ] x E N L ; M max a (ua(u u u N L ; M 2 max a (u. u u N L We now use estimates, obtained in [20], in order to estimate β N and γ N. First, we recall that [20, Lemma 3.] supplies us with a uniform bound for the spatial derivatives of R N : (6.4 s x s R N (, t L Const m s+ N log N, 0 s 2q. Using (6.4 with s =, 2 and the BV-boundedness (6.6, we conclude that (6.5 and (6.6 x R N x u N L x R N L u N BV Const m 2 N log N 2 x 2 R N x u N L 2 x 2 R N L u N BV Const m 3 N log N. Since E N L and x E N L are spectrally small, hence negligible, we conclude by (6.2, (6.3, (6.5, (6.6, (6.4 and (6.3 that (6.7 β N N θ( q 0, γ N N θ( 3 2q 0. 7
18 We may now state and prove the following weak Lip + -stability result: Theorem 6.. Consider the SV method (6., (6.3 (6.5, approximating the conservation law (. (.2. Assume that a = f satisfies ( a 0. Then the approximate solutions are -weakly Lip + -stable, with = N. Proof. Our assumption on a( implies that K, given in (6., equals zero. Hence, (6.0 reads in this case: (6.8 Solving (6.8 we get that d dt W (t (W (t2 + β N W (t + γ N. (6.9 W (t w + + w + w ηe (w + w t, (6.20 w ± = β N ± β 2 N + 4γ N 2, η = W (t = 0 w W (t = 0 w +. Note that w ± and η depend on N. Furthermore, by (6.20, (6.7 and (6.3 it follows that (6.2 Also, since by (6.7 (6.22 w ± = O ( γn N θ( 3 4q 2 0. N W (t = 0 N, we conclude by (6.20 and (6.2 that (6.23 η. N We claim that the weak Lip + -stability conditions, (. (.2, hold here with = N. Namely, (6.24 e W (tdt 0 O(N and ( W e (τdτ t dt O(log N. (6.26 In order to prove these two estimates we integrate (6.9 and find that t W (τdτ w (T t + log [ ηe (w + w ] T ηe (w + w t. Hence But since (6.27 [ ] T exp W (τdτ 0 e w T ηe(w + w T η η = w + w W (t = 0 w +, 8 = e w +T + e w T e(w + w T η.
19 we conclude that [ ] T exp W (τdτ 0 e w +T + e w T e(w + w T w + w (W (t = 0 w + and (6.24 follows by using (6.2 and (6.22. As for (6.25, inequality (6.26 implies (note that w 0: [ T ] T (6.28 exp W (τdτ dt 0 t ( ( T ηe (w+ w T + ηe(w+ w T ηe (w + w T log w + w η. First, we observe that (6.2 and (6.23 imply that (6.29 ηe (w+ w T N 0. Now, in order to estimate the second term on the right hand side of (6.28 we deal with each of its two multiplicands. Using (6.27, (6.2 and (6.22 we find that (6.30 ηe (w + w T w + w = e(w+ w T W (t = 0 w + + e(w+ w T w + w 0 + T = T. N Furthermore, by (6.27, (6.2 and (6.22, (6.3 ηe (w+ w T η = e (w+ w T + e(w+ w T w + w (W (t = 0 w + N. Hence, (6.28 (6.3 prove (6.25 and the proof is thus concluded. Corollary 6.2. Consider the SV method (6., (6.3 (6.5, approximating the conservation law (. (.2. Then (6.32 u N (, T u(, T W, { CT if u 0 Lip + < C T ln if u 0 Lip + = and ( a 0, with = N θ. Proof. The case of Lip + -bounded initial data is straightforward and we, therefore, concentrate on the case that u 0 Lip = and ( a 0. Since, by Theorem 6., we have N -weak Lip+ -stability in that case, and since N < = N θ, u N are also -weakly Lip + -stable. Hence, it remains to show -W, -consistency. W, -consistency with (., t u N + x f(u N K T N θ, L ([0,T ],W, (R x has already been shown in [20, (3.9b]. As for W, -consistency with the initial condition, we claim that (6.33a (6.33b u N (, 0 u(, 0 W, = V P N U 0 U 0 L K 0 N 2 log N, U 0 (x = x π 9 u 0 (ξdξ.
20 In order to prove (6.33, we recall that (consult [3, (2.2, (2.4, (2.5] (6.34 P N U 0 U 0 L Const log N N m U (m 0 L, m 0. Taking m = 2 in (6.34 we find that the initial error allowed by W, -consistency, is exhausted in this case: (6.35 P N U 0 U 0 L Const N 2 log N u 0 BV. We leave the reader to verify that (6.36 V P N U 0 P N U 0 L Const N 2 log N. Hence, (6.33 follows from (6.35, (6.36 and the proof is completed. We have thus far restricted our attention to the case ( a 0. In the general case, the cubic term on the right hand side of (6.0 does not vanish. Still, one can prove (along the lines of the proof of Lemma 3.2 weak Lip + -stability of order N = N θ log N, provided that W (t = 0 c N K for some c <. Alas, this condition does not hold in our case (consult (6.4 and (6.22. We, therefore, suggest to overcome this problem by considering a speed-like SV method, (6.37 t u N(x, t + x P N f(u N (x, t = N x Q N(x, t x a(u N (x, t with (6.3 (6.5 as before. This method, still conservative, differs from the regular SV method, (6., only in the spectral viscosity term on the right hand side, u N was replaced by a(u N. Remark (on an a-priori L bound. The question of uniform L -boundedness of this modified SV method may be tackled along the lines of [5]. However, we suggest here a simple argument which enables us to circumvent that question: Since the initial data are always assumed bounded, (.2, the exact entropy solution of (. (.2 will not be affected if we change the flux f outside the interval I 0 [min u 0, max u 0 ]. Therefore, we choose to smoothly extend f from I 0 to R, so that f, a = f,a,a, etc. remain uniformly bouded on R. By doing so we may conclude that f (i (u N, and by convexity, A (i (u N as well, i 0, are all uniformly bounded even if u N is not. Since our estimates depend only on f (i (u N L and A (i (u N L and never on the L -bound of u N itself, this argument is sufficient for our needs and no a-priori L -bound is required. We would like to comment that L -boundedness proofs for approximate solutions of (. (.2 may be sometimes tedious (as in our present case. Hence, it is sometimes customary to assume an a-priori L -bound, based, for instance, on numerical evidence. The above, to the best of our knowledge, innovative extension argument, may be applied to such approximations as well, so that assumptions, not fully justified, may be avoided. The convergence rate estimates for this modified SV method are given in the following theorem. Theorem 6.3. Consider the modified SV method (6.37, (6.3 (6.5, approximating the conservation law (. (.2. Then u N converges to the exact entropy solution u(x, t, as N, and for every T > 0 there exists a constant C T such that (6.38a u N (, T u(, T W, C T, (6.38b = { if u0 Lip + < ln if u 0 Lip + = and = N θ. 20
21 Proof. We first note that (6.37 is still L -stable (consult the proof of [20, Lemma 2.2] and hence (6.6 still holds. Therefore, (6.37 describes a family of conservative, L -stable and BV-bounded approximate solutions of (. (.2. Next, we address the question of weak Lip + -stability. We rewrite (6.37 as t u N + x f(u 2 N = N x 2 a(u N N x R N x a(u N + E N, R N and E N are as in (6.8. Multiplying by a (u N and differentiating with respect to x, we find that w = a(u N x satisfies (compare to (6.9: w t + a(u N w x + w 2 = N [a (u N w xx + a (u N A(u N ww x ] ( N [a (u N A(u N x R N ( x a(u N w + a 2 (u N x 2 R N ] x a(u N + We conclude that W (t = max x w(x, t satisfies +a (u N A(u N we N + a (u N x E N. d dt W (t + (W (t2 β N W (t + γ N, β N and γ N are not the same as in (6.2, (6.3 but still satisfy (6.7 (since a(u N BV remains uniformly bounded. This, according to the proof of Theorem 6., implies the N -weak Lip+ -stability of (6.37. Hence, by Theorem.2, error estimates (.7 hold with = N. Since it is easy to verify that our modified SV method is also W, -consistent of order N θ, error estimates (6.38 follow. Before concluding this section we consider the case of P N projecting the initial data, (6.2. We recall that the resulting approximation, u N, may not be bounded in BV and in fact u N BV may grow as much as O(log N. We note that this slightly changes our convergence rate results, stated in Corollary 6.2 and Theorem 6.3, so that (6.32 and (6.38 hold with = N θ log N, rather than = N θ. The first effect of replacing V P N by P N is that estimate (6.7 changes to β N N θ( q log N 0, γ N N θ( 3 2q log N 0, (consult (6.2, (6.3, (6.5 and (6.6. This, however, does not change the final result of -weak-lip + - stability with = N. Hence, by Theorem.2, error estimates (.7 still hold with = N. In view of (6.35 it remains only to consider the W, -consistency of u N (x, t with (.. Ignoring the spectrally small discretization error E N = x (I P Nf(u N, we obtain from (6. (the proof for (6.37 is similar that t u N (, t + x f(u N(, t W, N Q N (, t x u N(, t L. Using (6.4, Young inequality and the fact that Q N (, t L does not exceed O(log N (consult [20, (3.9b], we get (6.39 t u N(, t + x f(u N(, t W, N Q N (, t L u N (, t BV 2
22 Const N θ log N (log N2 = O(N θ log N. Hence, the order of W, -consistency reduced from O(N θ to O(N θ log N. Therefore, (.7, (6.35 and (6.39 imply an O(N θ log N convergence rate in W,. Appendix A. Appendix. Proof (of Lemma 3.2. By rescaling we may assume that K =. Since y (t is the solution of a perturbated Ricatti s equation, (3., we denote by y(t the solution of the regular Ricatti s equation, (A.a subject to the same initial condition, dy dt + y2 = 0, (A.b The solution of (A. is (A.2 y(0 = y (0 = c y(t =, 0 < c c c <. ( t + y, (0 while the solution of (3. (3.2 is given implicitly by (A.3 with (A.4 y (t = ( ( y t + D + ln, y D = ( y y (0 ln (0 y (0 First, we note that (3.2 and (3.3 imply that y (t is monotonically decreasing. Hence. (A.5 y (t y (0 t 0. Furthermore, since by (3. and (A.a (A.6 y (t y(t t 0, it follows, using (A.2 and monotonicity, (A.5, that (A.7 y (t y (T y(t = ( T + y t [0, T ]. (0 With the upper and lower bounds on y (t, (A.5 and (A.7, we may estimate the terms in (A.3 and (A.4. We start with the last term in the brackets in (A.3. Using (A.5 and (A.b it may be upper-bounded as follows, for all t 0: (A.8 ( y ln y ( ( ( = ln y ln y (0 ln c = O(. On the other hand, using (A.7 together with (A.b, we get a lower bound for this term: ( ( y T + y ln ln (0 (A.9 = O ( ln, 0 t T. y 22
23 Next, we estimate the constant D, given in (A.4. Using (A.b, (A.8 and (A.9 we obtain the following bounds: ( D T + c + ln y (0 (A.0 = O ( ln ; (A. D c + ln ( c = O(. (A.2 Hence we conclude by (A.3 and (A.8 (A. that y (t = [ t + O( + O ( ln ], 0 t T. With (A.2, estimates (3.4 and (3.5 may easily be verified. Indeed, ( e 0 y (tdt = T + O( + O ( ln O( + O ( ln O and and the proof is thus completed. 0 y e (τdτ t dt = 0 T + O( + O ( ln dt O ( ln, t + O( + O ( ln Proof (of Theorem 4.. Let u (x, t and v (x, t be two solutions of (4.. We assume that the regularization (4. is uniformly parabolic, Q p δ > 0, hence u and v are smooth. L -contraction for the degenerate case, Q p 0, easily follows by adding the term δp to the pseudo-viscosity coefficient Q(u, p and letting δ 0. As in [5], we divide the real line into intervals, R = n I n (t, I n (t = [x n (t, x n+ (t, so that ( n [u (, t v (A.3 (, t] 0 In(t and consequently (A.4 u (x n (t, t = v (x n (t, t. Using (A.3 and (A.4 we conclude that (A.5 d dt u (, t v (, t L (R = = d xn+(t ( n [u (x, t v (x, t]dx = dt x n (t n n xn+(t ( n [u t (x, t vt (x, t]dx. x n (t Using (4. and carrying out the integral on the right hand side of (A.5 we find that (A.6 d dt u (, t v (, t L(R = = n ( n[ ] xn+(t f(u + f(v + x n (t n 23 ( n[ ] xn+(t Q(u, u x Q(v, vx. x n (t
24 The first term on the right hand side of (A.6 vanishes in view of (A.4. Equality (A.4 also implies that the second term may be written as [ Q p (u, w [ ] ] x n+(t (A.7 ( n (u (x, t v (x, t, x n x n(t w is a mid-value between u x and vx. Since (A.3 implies that [ ( n (u (x, t v (x, t] x and x=x n+ (t [ ( n (u (x, t v (x, t] x x=x n(t and since Q p > 0, we conclude that (A.7 is non-positive. Therefore, (A.6 implies that d dt u (, t v (, t L (R 0, namely, the solution operator of (4. is L -contractive. 0 0 REFERENCES [] Y. Brenier, Roe s scheme and entropy solution for convex scalar conservation laws, INRIA Report 423, Le Chesnay, France 985. [2] Y. Brenier and S. Osher, The discrete one-sided Lipschitz condition for convex scalar conservation laws, SIAM J. Numer. Anal., 25 (988, pp [3] E. Harabetian, Rarefactions and large time behavior for parabolic equations and monotone schemes, Comm. Math. Phys., 4 (988, pp [4] S.N. Kružkov, First-order quasi-linear equations in several independent variables, Math. USSR Sb., 0 (970, pp [5] P.D. Lax, Hyperbolic Systems of conservation laws and the mathematical theory of shock waves, in Regional Conf. Series Lectures in Applied Math. Vol. (SIAM, Philadelphia, 972. [6] P. Le Floch and Z. Xin, Uniqueness via the adjoint problem for systems of conservation laws, to appear in Comm. on Pure and Appl. Math. [7] P. Marcati and R. Natalini, Convergence of the pseudo-viscosity approximation for conservation laws, to appear in Nolinear Analysis TMA. [8] Y. Maday and E. Tadmor, Analysis of the spectral viscosity method for periodic conservation laws, SIAM J. Numer. Anal., 26 (989, pp [9] J. von Neumann and R.D. Richtmeyer, A method for the numerical calculation of hydrodynamical shocks, J. Appl. Phys., 2 (950, pp [0] H. Nessyahu and E. Tadmor, The convergence rate of approximate solutions for nonlinear scalar conservation laws, SIAM J. Numer. Anal., 29 (992, pp [] H. Nessyahu, E. Tadmor and T. Tassa, The convergence rate of Godunov type schemes, TAU Report 5-92 (992, to appear in SIAM J. Numer. Anal. [2] O.A. Oleinik, Discontinuous solutions of nonlinear differential equations, Amer. Math. Soc. Trans. (2 26 (963, pp [3] A. Quarteroni, Some results of Bernstein and Jackson type for polynomial approximation in L P -spaces, Japan J. Appl. Math., (984, pp [4] P. Rosenau, Extending hydrodynamics via the regularization of the Chapman-Enskog expansion, Phys. Rev. A, 40 (989, pp [5] S. Schochet, The rate of convergence of spectral viscosity methods for periodic scalar conservation laws, SIAM J. Numer. Anal., 27 (990, pp [6] R.D. Richtmyer and K.W. Morton, Difference methods for initial-value problems, John Wiley & Sons, 967. [7] S. Schochet and E. Tadmor, The regularized Chapman-Enskog expansion for scalar conservation laws, Arch. Rational Mech. Anal., 9 (992, pp [8] E. Tadmor, Local error estimates for discontinuous solutions of nonlinear hyperbolic equations, SIAM J. Numer. Anal., 28 (99, pp [9] E. Tadmor, Convergence of spectral methods for nonlinear conservation laws, SIAM J. Numer. Anal., 26 (989, pp [20] E. Tadmor, Total variation and error estimates for spectral viscosity approximations, Math. Comp., 60 (993, pp [2] A.I. Volpert and S.I.Hudjaev, Cauchy s problem for degenerate second order quasi-linear parabolic equations, Math. USSR Sb., 7 (969, pp
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 informationSuper Viscosity and Spectral Approximations of Nonlinear Conservation Laws
umerical Methods for Fluid Dynamics. IV M. J. Baines and K. W. Morton, eds. Clarendon Press, Oxford, 1993, pp. 69-82 Super Viscosity and Spectral Approximations of onlinear Conservation Laws EITA TADMOR
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 informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
More 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 informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
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 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 information1. Introduction. We study the convergence of vanishing viscosity solutions governed by the single conservation law
SIAM J. NUMER. ANAL. Vol. 36, No. 6, pp. 1739 1758 c 1999 Society for Industrial and Applied Mathematics POINTWISE ERROR ESTIMATES FOR SCALAR CONSERVATION LAWS WITH PIECEWISE SMOOTH SOLUTIONS EITAN TADMOR
More informationConservation law equations : problem set
Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationExistence and Decay Rates of Solutions to the Generalized Burgers Equation
Existence and Decay Rates of Solutions to the Generalized Burgers Equation Jinghua Wang Institute of System Sciences, Academy of Mathematics and System Sciences Chinese Academy of Sciences, Beijing, 100080,
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 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 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 informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
More informationSHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS
SHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE Abstract. The present paper deals with the following hyperbolic elliptic coupled system, modelling
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 informationEffects of a Saturating Dissipation in Burgers-Type Equations
Effects of a Saturating Dissipation in Burgers-Type Equations A. KURGANOV AND P. ROSENAU Tel-Aviv University Abstract We propose and study a new variant of the Burgers equation with dissipation flues that
More informationLaplace 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 informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More 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 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 informationc1992 Society for Industrial and Applied Mathematics THE CONVERGENCE RATE OF APPROXIMATE SOLUTIONS FOR NONLINEAR SCALAR CONSERVATION LAWS
SIAM J. NUMER. ANAL. Vol. 9, No.6, pp. 505-59, December 99 c99 Society for Industrial and Applied Mathematics 00 THE CONVERGENCE RATE OF APPROXIMATE SOLUTIONS FOR NONLINEAR SCALAR CONSERVATION LAWS HAIM
More informationNonlinear Regularizing Effects for Hyperbolic Conservation Laws
Nonlinear Regularizing Effects for Hyperbolic Conservation Laws Ecole Polytechnique Centre de Mathématiques Laurent Schwartz Collège de France, séminaire EDP, 4 mai 2012 Motivation Consider Cauchy problem
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 informationOn 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 informationBose-Einstein Condensation and Global Dynamics of Solutions to a Hyperbolic Kompaneets Equation
Bose-Einstein Condensation and Global Dynamics of Solutions to a Hyperbolic Kompaneets Equation Joshua Ballew Abstract In this article, a simplified, hyperbolic model of the non-linear, degenerate parabolic
More informationWEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 94, pp. 1 14. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu WEAK ASYMPTOTIC
More 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 informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
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 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 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 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 informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
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 informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationGlobal unbounded solutions of the Fujita equation in the intermediate range
Global unbounded solutions of the Fujita equation in the intermediate range Peter Poláčik School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA Eiji Yanagida Department of Mathematics,
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 informationScaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations
Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.
More informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationThreshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations
Threshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455
More informationarxiv: v2 [math.fa] 17 May 2016
ESTIMATES ON SINGULAR VALUES OF FUNCTIONS OF PERTURBED OPERATORS arxiv:1605.03931v2 [math.fa] 17 May 2016 QINBO LIU DEPARTMENT OF MATHEMATICS, MICHIGAN STATE UNIVERSITY EAST LANSING, MI 48824, USA Abstract.
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
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 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 informationLARGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMERICS FOR 1-D SCALAR CONSERVATION LAWS
LAGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMEICS FO 1-D SCALA CONSEVATION LAWS L. I. IGNAT, A. POZO, E. ZUAZUA Abstract. In this paper we analyze the large time asymptotic behavior of the discrete
More 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 informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationFractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
Progress in Nonlinear Differential Equations and Their Applications, Vol. 63, 217 224 c 2005 Birkhäuser Verlag Basel/Switzerland Fractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationNECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES
NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES JUHA LEHRBÄCK Abstract. We establish necessary conditions for domains Ω R n which admit the pointwise (p, β)-hardy inequality u(x) Cd Ω(x)
More informationThe De Giorgi-Nash-Moser Estimates
The De Giorgi-Nash-Moser Estimates We are going to discuss the the equation Lu D i (a ij (x)d j u) = 0 in B 4 R n. (1) The a ij, with i, j {1,..., n}, are functions on the ball B 4. Here and in the following
More informationScaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations
Journal of Statistical Physics, Vol. 122, No. 2, January 2006 ( C 2006 ) DOI: 10.1007/s10955-005-8006-x Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations Jan
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 informationA RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS
A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS KATARINA JEGDIĆ, BARBARA LEE KEYFITZ, AND SUN CICA ČANIĆ We study a Riemann problem for the two-dimensional isentropic gas dynamics equations
More informationA COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 5, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A COUNTEREXAMPLE
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 informationConvergence of a first order scheme for a non local eikonal equation
Convergence of a first order scheme for a non local eikonal equation O. Alvarez, E. Carlini, R. Monneau, E. Rouy March 22, 2005 Abstract We prove the convergence of a first order finite difference scheme
More informationBanach Journal of Mathematical Analysis ISSN: (electronic)
Banach J. Math. Anal. 1 (2007), no. 1, 56 65 Banach Journal of Mathematical Analysis ISSN: 1735-8787 (electronic) http://www.math-analysis.org SOME REMARKS ON STABILITY AND SOLVABILITY OF LINEAR FUNCTIONAL
More informationLARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011
LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS S. G. Bobkov and F. L. Nazarov September 25, 20 Abstract We study large deviations of linear functionals on an isotropic
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationConvergence and sharp thresholds for propagation in nonlinear diffusion problems
J. Eur. Math. Soc. 12, 279 312 c European Mathematical Society 2010 DOI 10.4171/JEMS/198 Yihong Du Hiroshi Matano Convergence and sharp thresholds for propagation in nonlinear diffusion problems Received
More informationξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,
1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a
More informationOn Convergence of Minmod-Type Schemes
On Convergence of Minmod-Type Schemes Sergei Konyagin, Boan Popov, and Ognian Trifonov April 16, 2004 Abstract A class of non-oscillatory numerical methods for solving nonlinear scalar conservation laws
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationMAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS
MATHEMATICS OF COMPUTATION S 005-5780054-7 Article electronically published on May 6, 0 MAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS ORHAN MEHMETOGLU AND BOJAN POPOV Abstract.
More informationNumerische Mathematik
Numer. Math. (00) 87:70 735 Digital Object Identifier (DOI) 0.007/s0000000 Numerische Mathematik L -Stability and error estimates for approximate Hamilton-Jacobi solutions Chi-Tien Lin, Eitan Tadmor Department
More informationVariational formulation of entropy solutions for nonlinear conservation laws
Variational formulation of entropy solutions for nonlinear conservation laws Eitan Tadmor 1 Center for Scientific Computation and Mathematical Modeling CSCAMM) Department of Mathematics, Institute for
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More 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 informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
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 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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More 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 informationThreshold solutions and sharp transitions for nonautonomous parabolic equations on R N
Threshold solutions and sharp transitions for nonautonomous parabolic equations on R N P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract This paper is devoted to
More informationBoundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem
Boundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem Oleh Omel chenko, and Lutz Recke Department of Mathematics, Humboldt University of Berlin, Unter den Linden 6,
More informationP(E t, Ω)dt, (2) 4t has an advantage with respect. to the compactly supported mollifiers, i.e., the function W (t)f satisfies a semigroup law:
Introduction Functions of bounded variation, usually denoted by BV, have had and have an important role in several problems of calculus of variations. The main features that make BV functions suitable
More informationDUALITY SOLUTIONS FOR PRESSURELESS GASES, MONOTONE SCALAR CONSERVATION LAWS, AND UNIQUENESS. François Bouchut and François James
DUALITY SOLUTIONS FOR PRESSURELESS GASES, MONOTONE SCALAR CONSERVATION LAWS, AND UNIQUENESS François Bouchut and François James Comm. Partial Diff. Eq., 24 (1999), n 11-12, 2173-2189 Mathématiques, Applications
More informationContractive metrics for scalar conservation laws
Contractive metrics for scalar conservation laws François Bolley 1, Yann Brenier 2, Grégoire Loeper 34 Abstract We consider nondecreasing entropy solutions to 1-d scalar conservation laws and show that
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 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 informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationA LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION O. SAVIN. Introduction In this paper we study the geometry of the sections for solutions to the Monge- Ampere equation det D 2 u = f, u
More informationDynamics of Propagation and Interaction of Delta-Shock Waves in Conservation Law Systems
Dynamics of Propagation and Interaction of Delta-Shock Waves in Conservation Law Systems V. G. Danilov and V. M. Shelkovich Abstract. We introduce a new definition of a δ-shock wave type solution for a
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationEntropy-based moment closure for kinetic equations: Riemann problem and invariant regions
Entropy-based moment closure for kinetic equations: Riemann problem and invariant regions Jean-François Coulombel and Thierry Goudon CNRS & Université Lille, Laboratoire Paul Painlevé, UMR CNRS 854 Cité
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
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 information