NUMERICAL SCHEMES FOR COMPUTING DISCONTINUOUS SOLUTIONS OF THE DEGASPERIS-PROCESI EQUATION
|
|
- Ella Wright
- 6 years ago
- Views:
Transcription
1 Dept. of Math./CMA Univ. of Oslo Pure Mathematics No. 6 ISSN March 006 NUMEICAL SCHEMES FO COMPUTING DISCONTINUOUS SOLUTIONS OF THE DEGASPEIS-POCESI EQUATION G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO Abstract. ecent work [4] has shown that the Degasperis-Procesi equation is well-posed in the class of (discontinuous) entropy solutions. In the present paper we construct numerical schemes and prove that they converge to entropy solutions. Additionally, we provide several numerical examples accentuating that discontinuous (shock) solutions form independently of the smoothness of the initial data. Our focus on discontinuous solutions contrasts notably with the existing literature on the Degasperis-Procesi equation, which seems to emphasize similarities with the Camassa-Holm equation (bi-hamiltonian structure, integrabillity, peakon solutions, H 1 as the relevant functional space). 1. Introduction In this paper we present and analyze several numerical schemes for capturing discontinuous solutions of the Degasperis-Procesi equation [8] (1.1) t u 3 txxu + 4u x u = 3 x u xxu + u 3 xxxu, (x, t) (0, T ). We are interested in the Cauchy problem for this equation where an initial function u 0 is prescribed at time t = 0: u t=0 = u 0. The Degasperis-Procesi equation is a special case of the more general equation (1.) t u 3 txxu + 4 x f(u) = 3 xxxf(u), where (1.1) is recovered by choosing f(u) = 1 u. We shall concentrate on this equation, which we coin the generalized Degasperis-Procesi equation, and show that our numerical schemes converge if (1.3) f(u) κu for all u, for some constant κ > 0. About the motivation for considering (1.) instead of (1.1), besides the increased generality, we remark at this stage only that (1.) will play a central role in the analysis of numerical schemes for (1.1) as well. Degasperis and Procesi [8] studied the following six parameter family of third order dispersive equations t u + c 0 x u + γ 3 xxxu α 3 txxu = x ( c1 u + c ( x u) + c 3 u xxu ). Using the method of asymptotic integrability, they found that three equations of this family to be asymptotically integrable up to third order: the KdV equation Date: March 9, Mathematics Subect Classification. 35G5, 35L05, 35L65, 74S0, 65M1. Key words and phrases. Shallow water equation, integrable equation, hyperbolic equation, discontinuous solution, weak solution, entropy condition, numerical scheme, convergence. K. H. Karlsen is supported by an Outstanding Young Investigators Award. The authors are grateful to Hans Lundmark for interesting discussions. 1
2 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO (α = c = c 3 = 0), the Camassa-Holm equation (c 1 = (c 3 )/(α ), c = c 3 ) and the equation t u + x u + 6u x u + xxxu 3 α ( txxu xu xxu + 3 ) u 3 xxxu = 0. By rescaling, shifting the dependent variable, and applying a Galilean boost, this transforms to (1.1), see [7, 6] for details. Formally, equation (1.) is equivalent to the hyperbolic-elliptic system (1.4) t u + x f(u) + x P = 0, P xxp = 3f(u), which also provides the starting point for defining weak solutions of (1.). The Helmholtz operator (1 xx) 1 can be written as a convolution ( ) xx (g)(x) = e x y g(y) dy, and consequently P (x, t) = 3 e x y f(u(y, t)) dy. The Degasperis-Procesi equation (1.1) has a form similar to the Camassa-Holm shallow water wave equation (see, e.g., [, 5, 10, 19] and the references therein for more information about this equation), and many authors have emphasized that (1.1) share several properties with the Camassa-Holm equation such as bi- Hamiltonian structure, integrabillity, exact solutions that are a superposition of multipeakons, and H 1 as the relevant functional space for well-posedness. Let us also mention that a simplified model for radiating gases is given by a hyperbolicelliptic system that bears some similarities with (1.4), see, e.g., [1, 13, 18]. Degasperis, Holm, and Hone [7] proved the exact integrability of (1.1) and found exact non-smooth multipeakon solutions. Lundmark and Szmigielski [15, 16] used an inverse scattering approach to determine an explicit formula for the general n- peakon solution of the Degasperis-Procesi equation (1.1). Mustafa [17] showed that smooth solutions to (1.1) have the infinite speed of propagation property. Yin studied the globall well-posedness of the Degasperis-Procesi equation (1.1) in a series of papers [0, 1,, 3]. The (weak) solutions encompassed by Yin s well-posedness theory are all H 1 regular (typical examples are the peakons), a fact that is reminiscent of the Camassa-Holm equation (see, e.g., [5]). We refer to [10] for an overview of the Camassa-Holm, Degasperis-Procesi, and other related equations. There the authors present many numerical examples and include some comments on the relevance of these equations as shallow water models. In a different direction, two of the authors of the present paper promoted recently [4] the view that the Degasperis-Procesi equation could admit discontinuous (shock wave) solutions and that a well-posedness theory should rely on functional spaces containing discontinuous functions. Indeed, [4] proved the existence, uniqueness, and L 1 stability of so-called entropy solutions which belong to the space BV of functions of bounded variation. Lundmark [14] recently found some explicit shock solutions of the Degasperis- Procesi equation (1.1) that are entropy solutions (in the sense of [4]). A simple example of such a discontinuous solution to (1.1) is provided by { u(x, t) = 1+t ex t, x < t, t e (x t), x > t,
3 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 3 { 1 + e x, x < 0, which corresponds to the initial data u 0 (x) = 1 e x, x > 0. Since entropy solutions are of importance to us, let us recall their definition. Definition 1.1 (Entropy solution). A function u(x, t) is an entropy solution of the generalized Degasperis-Procesi equation (1.) with initial data u t=0 = u 0 L () provided the following two conditions hold: (1.5) i) u L ( (0, T )), and iii) for any convex C entropy η : with corresponding entropy flux q : defined by q (u) = η (u) f (u) there holds t η(u) + x q(u) + η (u) x P u 0 in D ( [0, T )), that is, φ Cc ( [0, T )), φ 0, (η(u) t φ + q(u) x φ η (u) x P u φ) dx dt + T 0 where P u (x, t) = 3 e x y f(u(y, t)) dy. η(u 0 (x))φ(x, 0) dx 0, From (1.3) and part i) of Definition 1.1, it is clear that x P u L ( + ) and thus the entropy formulation (1.6) makes sense. Another remark is that it is sufficient to verify (1.5) for the Kružkov entropies/entropy fluxes η(u) = u c, q(u) = sign (u c) (f(u) f(c)), c. Using the Kružkov entropies it can be seen that the following weak formulation (1.6) is a consequence of the entropy formulation (1.5): φ Cc ( [0, T )), T (1.6) (u t φ + f(u) x φ x P u φ) dx dt + u 0 (x)φ(x, 0) dx = 0. 0 egarding the L 1 stability (and thus uniqueness) of entropy solutions to the generalized Degasperis-Procesi equation, the following result is proved in [4]: Theorem 1.1 (L 1 stability). Suppose f C 1 (). Let u, v C([0, T ]; L 1 ()) be two entropy solutions of the generalized Degasperis-Procesi equation (1.) with initial data u 0, v 0 L 1 () L (), respectively. Then for any t (0, T ) (1.7) u(t, ) v(t, ) L 1 () ect u 0 v 0 L 1 (), where (1.8) C = 3 f L ([ M,M]), M = max [ u L ((0,T ) ), v L ((0,T ) ) Consequently, there exists at most one entropy solution to (1.) with initial data u t=0 = u 0 L 1 () L (). The existence of entropy solutions to (1.) is proved in [4] under the assumption that u 0 L 1 () BV (). These entropy solutions u satisfy u L (0, T ; L ()), u L (0, T ; L 1 () BV ()), where the latter L 1 BV bound, which depends on T, is a consequence of the former L bound. The L bound required by Definition 1.1 is a consequence of ].
4 4 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO the BV bound. The L bound, which is independent of T, is at the heart of the matter in [4], and formally comes from the conservation law (1.9) t (( xxv) + 5( x v) + 4v ) ( ) + x 3 u3 + 4vw + x v x w 4u v = 0, where v := (4 xx) 1 (u), w := (1 xx) 1 (u ). A consequence of this equation is that v(t, ) H () and thereby u(t, ) L (), for any t 0. Let us now turn to the topic of the present paper, namely numerical schemes for capturing discontinuous entropy solutions of the generalized Degasperis-Procesi equation (1.), which is to be interpreted in the sense of (1.4). First of all, the form of the equation (1.4) resembles a model of two phase flow in porous media. In a simplified form this model reads } t u + div (vf (u)) = 0 (1.10) (x, t) U (0, T ], div (Λ(u) P ) = 0 where U is a bounded domain in or 3, and Λ is a given function of u, and v = Λ(u) P. Here u [0, 1] denotes the saturation of one of the phases, and the flux function F is an increasing function F : [0, 1] [0, 1] with F (0) = 0 and F (1) = 1. The variable P represents the pressure. Since (1.10) is assumed to hold in a bounded domain, it is necessary to prescribe boundary values, at least for the total velocity v, as well as initial values for u. We remark that the well-posedness of the boundary/initial value problem for (1.10) is open, although results exist when capillary pressure forces (a degenerate diffusion operator is added to the equations), see [3]. Numerically, (1.10) is often solved using an operator splitting strategy in which the elliptic equation for P (the pressure equation) is solved to give a pressure which is used for the hyperbolic equation (the saturation equation), see for example [1] and the references therein. The rationale behind this is that the time scales in the pressure and saturation equations are sufficiently different, so that (presumably) the total velocity will change much slower than the saturation. This type of operator splitting inspired the construction of the numerical schemes presented in this paper. Informally, this splitting approximation is constructed as follows. Let S g t denote the solution operator for the balance law (1.11) t v + x f(v) + g(x) = 0, v(x, 0) = v 0 (x), i.e., we formally write the entropy solution v(x, t) as v(x, t) = S g t (u 0 )(x). Then we can define a sequence of functions {u n } n 0 by u 0 := u 0 and (1.1) u n := S gn 1 t (u n 1 ), g n 1 := x P n 1, P n 1 P n 1 xx = 3f ( u n 1). The resulting approximation u t = u t (x, t) is defined as u t (, t) = u n for n t t < (n + 1) t. Here t is a small parameter (the time step ). The hope is that u t u as t 0 in some reasonable topology, and that u is the entropy solution to (1.4). In this paper we establish this convergence. Additionally, we establish convergence of fully discrete numerical schemes constructed using the splitting approach and monotone difference schemes to discretize (1.11) and the equation in (1.1) for P n.
5 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 5 To make the topology for the conservation law and the elliptic equation for P compatible, we assume that f is globally Lipschitz continuous in the convergence analysis of the numerical schemes. It is difficult at the discrete level to establish (1.9), which would have provided us with a uniform L estimate on the numerical approximations. This in turn would have implied uniform BV, L 1, L estimates, and compactness of the numerical approximations would then have followed. On the other hand, if the nonlinearity f is globally Lipschitz continuous, then BV, L 1, L estimates can be derived directly without going via an L estimate. Thanks to the uniqueness of entropy solutions (Theorem 1.1), we can without loss of generality assume that f is globally Lipschitz continuous. If not, we can replace the function f by a globally Lipschitz continuous function, call it f L, such that f L coincides with f on [ L, L] for L > 0 and f L (u) κu for all u with κ given in (1.3). Let u L be the entropy solution of (1.) with f replaced by f L and initial data u L t=0 = u 0. It follows from the results in [4] that the L 1, L, L, BV bounds on u L coincides with those attached to u. Choosing L > M we conclude by the uniqueness of entropy solutions that the functions u and u L must coincide. In other words, thanks to the well-posedness analysis in [4] and without any loss of generality, we can restrict ourselves to the case where f in (1.) is globally Lipschitz continuous. The rest of this paper is organized as follows: In Section we establish some preliminary estimates for the balance equation where P is given by the initial data only. Furthermore, as a warm-up, we establish the convergence of a monotone difference scheme for this balance law. In Section 3 we define the operator splitting method. First we establish convergence of a semi-discrete splitting approximation, and then we proceed to show the convergence of fully discrete schemes. Finally, in Section 4 we present some numerical results. These examples show that shock solutions form independently of the smoothness of the initial data, and hence that the entropy solution framework is reasonable for this type of equations.. Preliminary estimates for balance laws As a starting point we consider the Cauchy problem for the conservation law with a source (which in this case is determined by the initial condition) t u + x f(u) + x P = 0, x, 0 < t < T, (.1) xxp + P = 3f(u 0 ), x, u(x, 0) = u 0 (x), x, where we shall assume that the initial condition satisfies (.) u 0 L 1 () BV (). As explained in the introduction, it is sufficient to consider nonlinearities f that are globally Lipschitz continuous. More precisely, we assume (.3) f C 1 (), f(0) = 0, f(ξ) f(ξ ) L ξ ξ, for all ξ, ξ and some positive constant L > Estimates on the exact solutions. First we collect some estimates satisfied by the source P in (.1). Lemma.1 (Estimates on P ). Assume that (.) and (.3) are satisfied. Then (.4) P L (), xp L () 3 L u 0 L 1 (),
6 6 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO (.5) P L 1 (), xp L 1 () 3L u 0 L 1 (), (.6) xxp 6L u L 1 () 0 L 1 (), (.7) xx P 3 L () L u 0 L 1 () + 3L u 0 BV (). Proof. The proof follows from the observation that e x / is the Green s function of the differential operator 1 xx, and thus we have that P (x) = 3 e x y f(u 0 (y)) dy, x, x P (x) = 3 e x y sign (y x) f(u 0 (y)) dy, x. Since e x 1, e x y dy =, for all x, by (.3), P (x), x P (x) 3 f(u 0 (y)) dy 3 L u 0 (y) dy, x, P (x) dx, x P (x) dx 3 f(u 0 (y)) dy 3L u 0 (y) dy, which give (.4), (.5), respectively. estimates. Using xxp = P 3f(u 0 ) gives the final Lemma.1 gives us the necessary control to prove some estimates for the entropy solution u of the conservation law (.1) with source x P. We collect these in the following lemma, whose proof is standard and therefore omitted (see [9]). Lemma. (Estimates on u). Assume that (.) and (.3) hold. Then for each t (0, T ) (.8) u(, t) L 1 () (1 + 3Lt) u 0 L 1 (), (.9) u(, t) L () u(, t) BV () u 0 BV () + 6Lt u 0 L 1 (). Additionally, there is a constant C T, independent of the discretization parameters, such that u(, t) u(, s) L 1 () C T t s, t, s [0, T ]... Finite difference schemes. We shall analyze the following finite difference discretization of (.1): (.10) U n+1 U n U 0 = 1 + D F (U n ; ) + DP = 0, Z, n {0,..., N 1}, ( 1 ) (+ 1 ) u 0 (x) dx, Z, D + D P + P = 3f ( U 0 ), Z, where, > 0 denote the temporal and spatial discretization parameters, respectively, while N is the smallest integer such that N t T. Moreover,
7 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 7 D +, D, D denote the forward, backward, and central spatial difference operators, respectively: (.11) D + ξ := ξ +1 ξ Dξ := D +ξ + D ξ, D ξ := ξ ξ 1, = ξ +1 ξ 1. The ratio / = λ is assumed to be constant. We use the following notation for the numerical flux function F : F (U; ) = F (U p, U p+1,..., U +q ), for some integers p, q 0. The numerical flux is assumed to be consistent, i.e., F (u,..., u) = f(u), and such that for the conservation law without source the scheme would be TVD. This means that (.1) ξ ξ 1 λ(f (ξ; ) + F (ξ; 1)) {ξ } bv, = for sufficiently small λ, and where {ξ } Z is any real sequence. Throughout this paper we use the notations bv {ξ } = l ξ ξ 1, {ξ } = l ξ, {ξ } 1 = sup ξ. ecall that for (.1) to hold, we must assume that the CFL-condition t λ holds, for a specific constant λ determined by f and u 0 [9]. In this paper, we assume that the TVD property (.1) is valid. There are numerous examples of numerical fluxes, see for example [9]. For the numerical examples in this paper we use the Engquist-Osher flux (.13) F (u, v) = 1 ( v ) f(u) + f(v) f (w) dw. Lemma.3 (Estimates on P ). Assume that (.) and (.3) hold. Then there exists a positive constant C 1, independent of f, u 0,,, such that (.14) {P } l, {D P } l, {D + P } l, {DP } l C 1 L u 0 L 1 (), {P } l1, {D P } l1, {D + P } l1, {DP } l1 C 1 L u 0, L1() u (.15) (.16) {D + D P } l1 C 1 L u 0. L1() Proof. The discrete equation for P can be solved explicitly: (.17) P = 3h e κ i f(ui 0 ), Z, where (.18) h := (1 + 1 ) e κ 1 (, κ := ln () 1 + () () ).
8 8 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO These formulas can be derived by applying the discrete Fourier transform to the equation determining P. Observe that e κ 1 (.19) h = O(),, 1 e κ 1 + O(), Then we can proceed as in the continuous case: (.0) D P = P P 1 = 3h = 3h = 3h e κ i e κ i +1 f(ui 0 ) e κ(i ) e κ(i +1) f(ui 0 ) + 3h i= κ(i ) 1 e κ e f(u i 0 ) + 3h i= 1 1 e κ 1. e κ(i ) e κ(i +1) f(ui 0 ) κ(i ) 1 eκ e f(u i 0 ). Using (.3) and (.19), for each Z and some constants c 1, c, c 3 > 0, P 3h e κ i f(ui 0 ) 3h f(ui 0 ) 3hL Ui 0 c 1 L D P 3h 3hc i= U 0 i c 1 L u 0 L 1 (), e κ(i ) 1 e κ 3c 3 L f(u 0 i ) + 3h e κ i f(u 0 i ) 3hc L U 0 i 3c 3 L u 0 L 1 (). 1 The remaining estimates are proved in the same way. e κ(i ) 1 eκ f(u 0 i ) U 0 i Now we have enough estimates on the discrete source term DP estimates on U n implying compactness. to establish Lemma.4 (Estimate on U n ). Assume that (.), (.3), and (.1) hold. Then (.1) {U n } l 1 (1 + C 1 Ln) u 0 L 1 (), (.) {U n } l {U n } bv {U 0 } bv + C 1 Ln u 0 L 1 (), for each n {0,..., N}, where C 1 is given in Lemma.3. Additionally, there exists a constant C T,L, independent of the discretization parameters, such that (.3) {U n } {U m } l 1 C T,L t n m, n, m {0,..., N}. Proof. This lemma is proved in the same ways as the corresponding result for a conservation law with a source [9].
9 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 9 Define the piecewise constant function u by (.4) u (x, t) := I n (x, t)u n, n 0, Z where I n(x) = I (x)i n (t), and I n is the characteristic function of the interval [n, (n + 1)) and I of the interval [( 1/), ( + 1/)). Similarly we define the piecewise linear function P by x (.5) P (x) = DP I (y) dy. Z Theorem.1 (Convergence). Assume that (.), (.3), and (.1) hold. Then as 0 (.6) u u in L p ( [0, T ]) p < and a.e. in [0, T ], where u is the unique entropy solution to (.1) in the sense of Kružkov (see [9, Section.1]). We also have that (.7) DP x P a.e. in + as 0, where P xxp = 3f(u 0 ) Proof. Due to (.1) and (.), u (, t) L 1 () (1 + C 1Lt) u 0 L 1 (), u (, t) L () u (, t) BV () {U 0 } bv + LtC 1 u 0 L 1 (), for any t (0, T ). Additionally, due to (.3), there is a constant C T,L, independent of the discretization parameters, such that u (, t) u (, t) L 1 () C T,L t s, t, s [0, T ]. Using a standard argument (see, e.g., [9]), there exists a limit function u and a sequence { l } l, with l 0 as l, such that (.8) (.9) u l u in L p ( [0, T ]) p < and a.e. in [0, T ], u(, t) L 1 (), u(, t) L (), u(, t) BV () C T,L, for each t (0, T ), where C T,L is a finite constant. In addition, u C([0, T ]; L 1 ()). Next we prove (.7) which will ensure that we have the right source term in the limit in the entropy formulation of the conservation law. Observe that (see (.17)) where P l = G l (3f(u l )), P = G (3f(u)), h x l i, G l (x) = h e κ G(x) = e x l, x, where denotes convolution in the x variable and [z] denotes the the integer part of z. By (.8), (.9), and the dominated convergence theorem, we have only to prove that (.30) G l (x) G(x), for each x as l, Due to (.18) and (.19) we have h 1 (.31) lim = lim l l l l + 1 = 1 e κ. l
10 10 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO Fixing x, due to (.18), (.3) i x h lim l e κ l = lim (e κ ) l = lim l h x i l (1 + ( l) = lim (1 + l ) l = lim (1 + l ) l + ) l 4 + (l ) h x i l x l = e x. h x i l Since (.31) and (.3) give (.30), (.7) is proved. Finally, we prove that u is an entropy solution to (.1). To this end, we use a standard argument. Let k N. Define the discrete entropy flux (.33) Q k (ξ; ) := Q k (ξ p,..., ξ +q ) where = F (ξ k; ) F (ξ k; ) = F (ξ p k,..., ξ +q k) F (ξ p k,..., ξ +q k), {ξ }, We begin by verifying that (.34) 1 a b := max{a, b}, a b := min{a, b}, a, b. ( U n+1 k U n k ) + D Q k (U n ; ) + DP sign ( U n+1 k ) 0. To show this we distinguish two cases. If U n+1 k, from (.10) and the monotonicity of the scheme (see [9, Theorem 3.6]) we get U n+1 k = U n+1 k = U n k λ ( F (U n ; ) F (U n ; 1) ) DP U n k λ ( Qk (U n ; ) Q k (U n ; 1) ) DP = U n k λ ( Qk (U n ; ) Q k (U n ; 1) ) DP sign ( U n+1 k ). If U n+1 < k, from (.10) and the monotonicity of the scheme (see [9, Theorem 3.6]) we get U n+1 k U n k ( λ Qk (U n ; ) Q k (U n ; 1) ) DP sign ( U n+1 k ), which proves (.34). Fix φ Cc ( [0, )), φ 0, and multiply (.34) by φ(, n) to get n+1 U k U n (.35) φ(, n) k n=0 = + + n=0 = n=0 = φ(, n)d Q k (U n ; ) φ(, n)dp sign ( U n+1 k ) 0.
11 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 11 Summing by parts yields n=0 = = n=1 = n+1 U k U n φ(, n) k = n=0 = hence (.35) becomes (.36) We claim that (.37) where φ(, n) φ(, (n 1)) U n k φ(, 0) U 0 k, φ(, n)d Q k (U n ; ) = = n=1 = + n=0 = Q k (U n ; )D + φ(, n), φ(, n) φ(, (n 1)) U n k n=0 = n=0 = = Q k (U n ; )D + φ(, n) φ(, n)dp sign ( U n+1 k ) φ(, 0) U 0 k 0. Qk (U n ; ) q k (U n ) 0, as 0 for each n N, (.38) q k (ξ) := sign (ξ k) (f(ξ) f(k)), ξ, is the flux associated to the entropy ξ ξ k. By the definition of the discrete entropy flux Q k, we get q k (U n ) = Q k (U n,..., U n ), and hence, from (.), there is a constant C L,T, independent of the discretization parameters, such that = Qk (U n ; ) q k (U n ) CT,L. Let us recall the following well-known fact (see for example [11]). Let Ω d be a bounded open set, g L 1 (Ω), and suppose that g ν (x) g(x) a.e. in Ω as ν. Then there exists an at most countable set Θ such that for any k \ Θ sign (g ν (x) k) sign (g(x) k) a.e. in Ω.
12 1 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO Using this fact, (.9), (.37), (.7) and the dominated convergence theorem, and sending l in (.36), we get (.39) ( u k t φ + q k (u) x φ sign (u k) x P φ) dx dt 0 u 0 (x) k φ(0, x) dx 0. for any k \ Θ, where Θ is an at most countable set Θ. To conclude from this that (.39) holds for all k we repeat the argument in [11]. Consequently, u is the Kružkov entropy solution to (.1). Finally, convergence of the whole sequence {u } >0 follows from the uniqueness of the entropy solution to (.1) [9, Section.4]. 3. Numerical schemes based on operator splitting In this section we construct numerical schemes for the generalized Degaperis- Procesi equation (1.). As explained in the introduction, there is no loss of generality in assuming that nonlinearity f is globally Lipschitz continuous. The idea is to fix P for a small interval t, and let u solve the conservation law with source (.1) for t [0, t). Then we update P and proceed to the next interval. We start by doing this for the semi-discrete scheme, where we use exact solutions of the conservation law and of the elliptic equation The semi-discrete approximation. We consider the problem t u + x f(u) + x P = 0, x, 0 < t < T, (3.1) P xxp = 3f(u), x, 0 < t < T, u(x, 0) = u 0 (x), x, where we assume that (.) and (.3) hold. We fix a positive number t and set N to be the smallest integer such that N t T. Then we define two sequences {u n } N n=1 and {vn } N n=1 by t v n + x f(v n ) + x P n 1 = 0, x, 0 < t < t, (3.) P n 1 xxp n 1 = 3f(u n 1 ), x, v n (x, 0) = u n 1 (x), x, and (3.3) u n ( ) := v n (, t). We initialize by setting u 0 = u 0. In the following we show that this sequence converges to an entropy solution of (3.1) as t 0, see Definition 1.1. To do this we establish the usual estimates. Lemma 3.1. Assume that (.) and (.3) hold. Then (3.4) (3.5) u n L 1 () e 3LT u 0 L 1 (), u n L () un BV () u 0 BV () + 6LT e 3LT u 0 L 1 (). for each n {0,..., N}. Additionally, there exists a constant C T,L, independent of the discretization parameters, such that (3.6) u n u m L 1 () C T,L t n m, n, m {0,..., N}.
13 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 13 Proof. Due to (.8), hence, from (3.3), v n (, t) L 1 () (1 + 3Lt) u n 1 L 1 (), u n L 1 () (1 + 3L t) u n 1 L 1 (), which implies (3.4). We continue by proving (3.5). By (.9), 0 t t, v n (, t) BV () u n 1 BV () + 6Lt u n 1 L 1 (), 0 t t, hence, from (3.3), Therefore, u n BV () u n 1 BV () + 6L t u n 1 L 1 (). n 1 u n BV () u 0 BV () + 6L t u k L 1 () k=0 u 0 BV () + 6LT e 3LT u 0 L 1 (), which is (3.5). Finally, (3.6) is a consequence of (3.5) and the fact that the source x P n 1 in (3.) is bounded in L, see (.4). For the convergence analysis, we introduce an approximation defined for all t by (3.7) u t (, t) := u n, n t t < (n + 1) t, n {0,..., N 1}. Theorem 3.1. Assume that (.) and (.3) hold. Then as t 0 (3.8) u t u in L p ( [0, T ]) p < and a.e. in [0, T ], where u is the unique entropy solution to (3.1) in the sense of Definition 1.1. Proof. By (3.4) and (3.5), (3.9) u t (, t) L 1 () e3lt u 0 L 1 (), u t (, t) L () u t(, t) BV () u 0 BV () + 6LT e 3LT u 0 L 1 (), for any 0 t T. Additionally, due to (3.6), there exists a constant C T,L, independent of the discretization parameters, such that (3.10) u t (, t) u t (, s) L 1 () C T,L t s, t, s [0, T ]. Using a standard argument (see, e.g., [9]), there exists a limit function u and a sequence { t l } l, with t l 0 as l, such that (3.11) (3.1) u tl u in L p ( [0, T ]) p < and a.e. in [0, T ], u(, t) L 1 (), u(, t) L (), u(, t) BV () C T,L, for each t (0, T ), where C T,L is a finite constant. In addition, u C([0, T ]; L 1 ()). To continue, we introduce an auxiliary function v t (x, t) defined by (3.13) v t (, t) := v n (, t n t), (n 1) t < t n t, n {1,..., N}, where v n solves (3.). Using (3.4), (3.5), and (3.10), there exists a constant C T,L, independent of the discretization parameters, such that for all t, s [0, T ] v t (, t) L 1 () C T,L,
14 14 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO (3.14) (3.15) v t (, t) L () v t(, t) BV () C L,T, v t (, t) v t (, t) L 1 () C T,L t s. We claim that (3.16) u t (, t) v t (, t) L 1 () 0, as t 0, for each 0 t T. To this end, let t [(n 1) t, n t). Then, since v n (, 0) = u n 1 and thanks to (3.10) and (3.15), u t (, t) v t (, t) L 1 () = u n 1 v n (, t n t) L 1 () = v n (, 0) v n (, t n t) L 1 () C T,L (t n t) C T,L t, which proves (3.16). Thus v t it is the unique entropy solution to t v t + x f(v t ) + x P t = 0, x, 0 < t < T, (3.17) P t xxp t = 3f(u t ), x, 0 < t < T, v t (x, 0) = u 0 (x), x. Hence, for any convex entropy/entropy flux pair (η, q) and for any nonnegative test function φ Cc ( [0, T )) we have, referring to (3.11) for the converging subsequence { t l }, T (3.18) (η(v tl ) t φ + q(v tl ) x φ η (v tl ) x P tl φ) dx dt 0 + η(u 0 (x))φ(0, x) dx 0. By (3.11) and (3.16), (3.19) v tl u a.e. in [0, T ] as l. Due to (3.11), (3.1), and the fact that P tl (x, t) = 3 e x y f(u tl (y, t)) dy, we conclude (3.0) x P tl x P a.e. in [0, T ] as l, where P xxp = 3f(u). Using (3.9), (3.14), (3.19), (3.0) and the dominated convergence theorem, sending l in (3.18) yields that the limit u satisfies T (η(u) t φ + q(u) x φ η (u) x P u φ) dx dt + η(u 0 (x))φ(0, x) dx 0. 0 Then u is the entropy solution to (3.1) in the sense of Definition 1.1. Finally, convergence of the whole sequence {u t } t>0 follows from the uniqueness of the entropy solution to (3.1) (see Theorem 1.1).
15 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION Fully discrete schemes. In this subsection we introduce and analyze fully discrete schemes for (3.). To this end, we introduce a time step t > 0 and a local time step < t and let N( t) = N N, M( t, ) = M N be the smallest integers satisfying (3.1) M = t, N t = T. Define (3.) U 0 := 1 and for n {1,..., N} (3.3) V n,m+1 V n,m D + D P n 1 V n,0 = U n 1, (+ 1 ) where Z, m {0,..., M 1}, and ( 1 ) u 0 (x)dx, Z, + D F (V n,m ; ) + DP n 1 = 0, + P n 1 = 3f ( U n 1 ), (3.4) U n := V n,m, Z, where we use the notations introduced in Section.. As before, we assume that F is a consistent numerical flux such that the scheme for the conservation law (without P ) would be monotone. ecall that this implies that a CFL-condition must hold. Lemma 3.. Assume that (.), (.3), and (.1) hold. Then { U n } (3.5) e C1LT u 0 l 1 L 1 (), { U n } { U n } { U 0 } (3.6) + C 1 LT e C1LT u 0 l L 1 (). bv bv for each n {0,..., N}. Additionally, there exists a constant C T,L, independent of the discretization parameters, such that (3.7) {U n } {U m } l 1 C T,L t n m, n, m {0,..., N}. Proof. We begin by proving (3.5). By (.1) we have (3.8) { V n,m } (1 + C 1 Lm) { U n 1 } l 1 l 1, 0 m M 1. Hence, from (3.1), (3.4), { U n } (1 + C 1 LM) { U n 1 } = (1 + C 1 L t) { U n 1 } l 1 l 1 Thus (3.5) holds. We continue by proving (3.6). Due to (.) we have (3.9) { V n,m } { U n 1 } + C 1 Lm { U n 1 }, bv bv so, from (3.1) and (3.4), we obtain { U n } { U n 1 } + C bv 1 LM { U n 1 } bv = { U n 1 } + C 1 L t { U n 1 }. bv Now (3.6) follows as in the proof of Lemma 3.1. l 1 l 1 l 1 l 1.
16 16 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO The proof of (3.7) is a consequence of the following calculation: m V n,i V n,i+1 i=0 = λ m i=0 = + m ( F = V n,i p DP n,..., V n,i +q ) ( F V n,i n,i p 1,..., V+q 1) { } c 1 sup V n,i i + m { DP n 1 } bv l 1 ( { } ) c 1 U 0 + C 1 LT e C1LT u 0 L 1 () + mc 1 Le C1LT u 0 L 1 (), BV for any n {1,..., N}, m {0,..., m 1}. Set = ( t, ), and define u (x, t) =,n I (x)i n (t)u n, where I n now denotes the characteristic function of the interval [n t, (n+1) t) and I (x) is still the characteristic function of the interval [( 1/), ( + 1/)). Theorem 3.. Assume (.), (.3), and (.1) hold. Then as t, 0, with t 0, (3.30) u u in L p ( [0, T ]) p < and a.e. in [0, T ], where u is the unique entropy solution to (3.1) in the sense of Definition 1.1. Proof. Due to (3.5), (3.6), u (, t) L 1 () ec1lt u 0 L 1 (), u (, t) L () u (, t) BV () { U 0 } + C 1 LT e C1LT u 0 L 1 (), bv for any 0 t < T. Additionally, due to (.3), there exists a constant C T,L, independent of the discretization parameters, such that u (, t) u (, t) L 1 () C T,L t s, t, s [0, T ]. Using a standard argument, there exists a limit function u and two sequences { t l } l, { l } l, l := ( t l, l ) 0 such that (3.31) (3.3) u l u in L p ( [0, T ]) p < and a.e. in [0, T ], u(, t) L 1 (), u(, t) L (), u(, t) BV () C T,L, for each t (0, T ), where C T,L is a finite constant. In addition, u C([0, T ]; L 1 ()). Define v (x, t) = I n (x, t)j m (t n t)v n,m n,m 0, Z with J m denoting the characteristic interval of [m, (m + 1)). Here we recall that V n,m solves (3.3). Using (3.5), (3.6), (3.7), (3.8), (3.9) we see that (3.33) v (, t) L 1 () ec1lt u 0 L 1 (),
17 (3.34) NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 17 v (, t) BV () for any 0 t < T. We claim that { U 0 } + C 1 LT e C1LT u 0 L 1 (), bv (3.35) u (, t) v (, t) L 1 () 0, as 0, for each 0 t T. Let 0 t T and n {0,..., N 1}, m {0,..., M 1}, such that n t + m t < n t + (m + 1). For the sake of simplicity we consider only the case m 0, n t + m < t < n t + (m + 1). Now we find that (see the proof of (3.)) u (, t) v (, t) L 1 () = V n+1,0 V n+1,m+1 = m i=0 = V n+1,i V n+1,i+1 C T,L τ where C T,L > 0 independent of the discretization parameters, and (3.35) thus holds. We continue by showing that u is an entropy solution to (3.1). Let k be any real number and define Q k, q k as in the proof of Theorem.1, see (.33) and (.38). Arguing as for (.34) we get V n,m+1 k V n,m k ( ) + D Q k (V n,m ; ) + DP n 1 sign V n,m+1 k 0. For a fixed test function φ Cc ( [0, T )), φ 0, we compute n,m+1 V k V n,m k (3.36) φ(, t n,m ) m=0 n=0 = + + m=0 n=0 = m=0 n=0 = φ(, t n,m )D Q k (V n,m ; ) ( ) φ(, t n,m )DP n 1 sign V n,m+1 k 0, where t n,m := n t + m. Summing by parts and using (3.1) produce n,m+1 V k V n,m k φ(, t n,m ) m=0 n=0 = = m=1 n=0 = + n=0 = n=0 = m=1 n=0 = φ(, t n,m ) φ(, t n,m 1 ) V n,m k n,m V k φ(, t n, ) n,0 V k φ(, t n,0 ) φ(, t n,m ) φ(, t n,m 1 ) V n,m k
18 18 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO = = + n=0 = n=0 = m=1 n=0 = + n=0 = n=0 = m=1 n=0 = + + n=1 = = = m=0 n=0 = = n+1,0 V k φ(, t n, ) n,0 V k φ(, t n,0 ) φ(, t n,m ) φ(, t n,m 1 ) V n,m k n+1,0 V k φ(, t n+1, 1 ) n,0 V k φ(, t n,0 ) φ(, t n,m ) φ(, t n,m 1 ) V n,m k φ(, n t ) φ(, n t) V n,0 k φ(, N t ) V N,0 k φ(, 0) U 0 k, φ(, t n,m )D Q k (V n,m ; ) m=0 n=0 = Q k (V n,m ; )D + φ(, t n,m ), hence, due to the identity N t = T (see (3.1)), (3.36) becomes (3.37) m=1 n=0 = + + n=1 = = φ(, t n,m ) φ(, t n,m 1 ) V n,m k φ(, n t ) φ(, n t) V n,0 k φ(, T ) V N,0 k φ(, 0) U 0 k =
19 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 19 + m=0 n=0 = m=0 n=0 = Q k (V n,m ; )D + φ(, t n,m ) ( ) φ(, t n,m )DP n sign V n,m+1 k 0. Since supp (φ) is a compact subset of [0, T ), we have (, T ) supp (φ) for small enough, and since V N,0 is bounded we find that (3.38) Furthermore (3.39) = n=1 = φ(, T ) V N,0 k 0, as t, 0. φ(, n t ) φ(, n t) t T { } V n,0 k,n l V n,0 k sup t φ(, t) L 1 () 0. 0 t T Next, thanks to the bv estimate (3.9) and the consistency of the numerical flux, (3.40) m=0 n=0 = As before, we have that Qk (V n,m ; ) q k (V n,m ) 0, as t, 0. (3.41) DP tl, l x P a.e. in [0, T ] as l, where P xxp = 3f(u), and x P tl, l (x, t) = DP n I n (y, t) dy. Using (3.31)-(3.35), (3.38)-(3.41) and the dominated convergence theorem when sending l in (3.37), we find that the limit u satisfies the entropy inequality (1.5) for all k \ Θ, where Θ is an at most countable set Θ (see the proof of Theorem.1). To conclude from this that (.39) holds for all k, i.e., that u is an entropy solution to (3.1) in the sense of Definition 1.1, we repeat the argument in [11]. Again, convergence of the whole sequence {u t } t>0 follows from the uniqueness of entropy solutions. 4. Numerical examples In this section we test our scheme on two examples. In order to compare our results with an exact solution we use the shock-peakons described by Lundmark in [14]. These are exact entropy solutions to the Degasperis-Procesi equation (1.1) and are given by the formulas (4.1) u(x, t) = N {m k (t) s k (t)sign (x x k (t))} e x xk(t), k=1
20 0 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO where the coefficients x k, m k, and s k satisfy the following system of ordinary differential equations: (4.) (4.3) (4.4) d dt x k(t) = u (x k, t), d [ dt m k(t) = s k (t)u (x k, t) m k (t) d dt s k(t) = s k (t)) N ] (s (t) sign (x k (t) x (t)) m (t)) e x k(t) x (t), =1 N (s (t) sign (x k (t) x (t)) m (t)) e x k(t) x (t). =1 Unless N = 1 it is not possible to find explicit solutions to this system of equations, and we used Matlab s routine ode3s (after some trials) to integrate the system numerically. The result of this was then used as an exact solution with which we compared the approximations calculated with our numerical scheme. If N = 1 we have three different types of traveling waves, the peakon the antipeakon u(x, t) = e x t, u(x, t) = e x t, and the shockpeakon sign (x) e x u(x, t) =. 1 + t We have tested the scheme using the Engquist-Osher numerical flux function (.13). For our test we used N = 3 and the initial values x 1 (0) = 5, x (0) = 0, x 3 (0) = 5, m 1 (0) = 1, m (0) = 0, m 3 (0) = 1, s 1 (0) = 0, s (0) = 1, s 3 (0) = 0, which gives the initial function (4.5) u 0 (x) = e x+5 + sign (x) e x e x 5. The exact solution will be a peakon, shockpeakon, antipeakon collision at t = 5 and x = 0. For t > 5 the solution continues as a decaying shockpeakon. We used M uniformly spaced grid points in the interval x [ 10, 10]. The relationship between and t were such that (M/4) = t, and was determined according to the CFL-condition = 0.5 max { U n In Figure 1 we show the approximate solution obtained using 8 gridpoints in [ 10, 10] for t [0, 8]. Since we have an (almost) exact solution to compare with, we tabulated the relative l 1 error, defined by U n (4.6) err = max un n, where u n denotes the exact solution at (x, t n ). In Table 1 we report these errors. u n }.
21 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION !0.!0.4!0.6!0.8 0!8!6!4! Figure 1. The approximation to with initial data (4.5) and 56 gridpoints in [ 10, 10]. M err Table 1. elative l 1 errors. Here = 0 M and err is defined in (4.6). The initial data are given by (4.5). This computations confirm the convergence established in the previous section, and give a numerical convergence rate of about 1/3. We remark that the observed numerical convergence is improved if = t, but we have not been able to prove that the resulting scheme converges. For the solutions of the Degasperis-Procesi equation given by (4.1) shocks cannot form from continuous initial data. If s k = 0 for all k, then the third equation gives s k (t) = 0, and thus s k remains zero unless m i blows up. We observed such singular behavior numerically for peakon-antipeakon collisions, but it was not possible to integrate the system past the collision time. Despite this we observed that shock formation was generic if u 0(x) < 0 and u 0 (x) < 0 for some x. To illustrate this we show a computed example where shocks form as they would in the conservation law. In this case the initial function is given by (4.7) u 0 (x) = e 0.5x sin(πx), for x [, ], where we assume that u 0 is extended periodically outside this interval. In Figure we show the computed solution using 8 gridpoints in the
22 G. M. COCLITE, K. H. KALSEN, AND N. H. ISEBO interval [, ] for t [0, 1.3]. We see that the N-waves familiar from Burgers equation do indeed form, even though the initial data is continuous !1 0.4! 0.!3 0!1.5!1! Figure. An approximate solution with initial data (4.7). eferences [1] F. Bratvedt, K. Bratvedt, C. F. Buchholz, T. Gimse, H. Holden, L. Holden, and N. H. isebro. FONTLINE and FONTSIM: two full scale, two-phase, black oil reservoir simulators based on front tracking. Surveys Math. Indust., 3(3):185 15, []. Camassa and D. D. Holm. An integrable shallow water equation with peaked solitons. Phys. ev. Lett., 71(11): , [3] Z. Chen. Degenerate two-phase incompressible flow. I. Existence, uniqueness and regularity of a weak solution. J. Differential Equations, 171():03 3, 001. [4] G. M. Coclite and K. H. Karlsen. On the well-posedness of the Degasperis-Procesi equation. J. Funct. Anal., 33(1):60 91, 006. [5] A. Constantin and L. Molinet. Global weak solutions for a shallow water equation. Comm. Math. Phys., 11(1):45 61, 000. [6] A. Degasperis, D. D. Holm, and A. N. W. Hone. Integrable and non-integrable equations with peakons. In Nonlinear physics: theory and experiment, II (Gallipoli, 00), pages World Sci. Publishing, iver Edge, NJ, 003. [7] A. Degasperis, D. D. Holm, and A. N. I. Khon. A new integrable equation with peakon solutions. Teoret. Mat. Fiz., 133(): , 00. [8] A. Degasperis and M. Procesi. Asymptotic integrability. In Symmetry and perturbation theory (ome, 1998), pages World Sci. Publishing, iver Edge, NJ, [9] H. Holden and N. H. isebro. Front tracking for hyperbolic conservation laws, volume 15 of Applied Mathematical Sciences. Springer-Verlag, New York, 00. [10] D. D. Holm and M. F. Staley. Wave structure and nonlinear balances in a family of evolutionary PDEs. SIAM J. Appl. Dyn. Syst., (3): (electronic), 003. [11] K. H. Karlsen, N. H. isebro, and J. D. Towers. Front tracking for scalar balance equations. J. Hyperbolic Differ. Equ., 1(1): , 004.
23 NUMEICAL SCHEMES FO THE DEGASPEIS-POCESI EQUATION 3 [1] C. Lattanzio and P. Marcati. Global well-posedness and relaxation limits of a model for radiating gas. J. Differential Equations, 190(): , 003. [13] H. Liu and E. Tadmor. Critical thresholds in a convolution model for nonlinear conservation laws. SIAM J. Math. Anal., 33(4): (electronic), 001. [14] H. Lundmark. Formation and dynamics of shock waves in the Degasperis Procesi equation. Preprint, 006. [15] H. Lundmark and J. Szmigielski. Multi-peakon solutions of the Degasperis-Procesi equation. Inverse Problems, 19(6): , 003. [16] H. Lundmark and J. Szmigielski. Degasperis-Procesi peakons and the discrete cubic string. IMP Int. Math. es. Pap., ():53 116, 005. [17] O. G. Mustafa. A note on the Degasperis-Procesi equation. J. Nonlinear Math. Phys., 1(1):10 14, 005. [18] D. Serre. L 1 -stability of constants in a model for radiating gases. Commun. Math. Sci., 1(1):197 05, 003. [19] Z. Xin and P. Zhang. On the weak solutions to a shallow water equation. Comm. Pure Appl. Math., 53(11): , 000. [0] Z. Yin. Global existence for a new periodic integrable equation. J. Math. Anal. Appl., 83(1):19 139, 003. [1] Z. Yin. On the Cauchy problem for an integrable equation with peakon solutions. Illinois J. Math., 47(3): , 003. [] Z. Yin. Global solutions to a new integrable equation with peakons. Indiana Univ. Math. J., 53(4): , 004. [3] Z. Yin. Global weak solutions for a new periodic integrable equation with peakon solutions. J. Funct. Anal., 1(1):18 194, 004. (Giuseppe Maria Coclite) Department of Mathematics, University of Bari Via E. Orabona 4, 7015 Bari, Italy address: coclitegm@dm.uniba.it (Kenneth Hvistendahl Karlsen) Centre of Mathematics for Applications (CMA), University of Oslo P.O. Box 1053, Blindern N 0316 Oslo, Norway address: kennethk@math.uio.no UL: (Nils Henrik isebro) Centre of Mathematics for Applications (CMA), University of Oslo P.O. Box 1053, Blindern, N 0316 Oslo, Norway address: nilshr@math.uio.no UL:
Proceedings of WASCOM 2005, pages World Sci. Publishing, 2006
Proceedings of WASCOM 2005, pages 128 133. World Sci. Publishing, 2006 A SEMIGROUP OF SOLUTIONS FOR THE DEGASPERIS-PROCESI EQUATION G. M. COCLITE AND K. H. KARLSEN Centre of Mathematics for Applications,
More informationNumerical methods for conservation laws with a stochastically driven flux
Numerical methods for conservation laws with a stochastically driven flux Håkon Hoel, Kenneth Karlsen, Nils Henrik Risebro, Erlend Briseid Storrøsten Department of Mathematics, University of Oslo, Norway
More informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More informationAbstract. A front tracking method is used to construct weak solutions to
A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar
More informationVISCOSITY SOLUTIONS OF HAMILTON JACOBI EQUATIONS WITH DISCONTINUOUS COEFFICIENTS
Dept. of Math. Univ. of Oslo Pure Mathematics No. 2 ISSN 86 2439 January 25 VISCOSITY SOLUTIONS OF HAMILTON JACOBI EQUATIONS WITH DISCONTINUOUS COEFFICIENTS GIUSEPPE MARIA COCLITE AND NILS HENRIK RISEBRO
More 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 informationDerivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations
Derivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations H. A. Erbay Department of Natural and Mathematical Sciences, Faculty of Engineering, Ozyegin University, Cekmekoy 34794,
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 informationJournal of Hyperbolic Differential Equations, Vol. 1, No. 1 (2004),
FONT TACKING FO SCALA BALANCE EQUATIONS K. H. KALSEN, N. H. ISEBO, AND J. D. TOWES Journal of Hyperbolic Differential Equations, Vol. 1, No. 1 (2004), 115-148 Abstract. We propose and prove convergence
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 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 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 informationON TRAVELING WAVE SOLUTIONS OF THE θ-equation OF DISPERSIVE TYPE
ON TRAVELING WAVE SOLUTIONS OF THE -EQUATION OF DISPERSIVE TYPE TAE GAB HA AND HAILIANG LIU Abstract. Traveling wave solutions to a class of dispersive models, u t u txx + uu x = uu xxx +(1 )u x u xx,
More informationWave Breaking Phenomena and Stability of Peakons for the Degasperis-Procesi
Wave Breaking Phenomena and Stability of Peakons for the Degasperis-Procesi Yue Liu Technical Report 2008-10 http://www.uta.edu/math/preprint/ Wave Breaking Phenomena and Stability of Peakons for the Degasperis-Procesi
More informationVariational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized
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 informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More informationPERSISTENCE PROPERTIES AND UNIQUE CONTINUATION OF SOLUTIONS OF THE CAMASSA-HOLM EQUATION
PERSISTENCE PROPERTIES AND UNIQUE CONTINUATION OF SOLUTIONS OF THE CAMASSA-HOLM EQUATION A. ALEXANDROU HIMONAS, GERARD MISIO LEK, GUSTAVO PONCE, AND YONG ZHOU Abstract. It is shown that a strong solution
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 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 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 informationConvergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates
Convergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates Xin Wen and Shi Jin Abstract We study the L 1 -error estimates
More informationGLOBAL WELL-POSEDNESS FOR NONLINEAR NONLOCAL CAUCHY PROBLEMS ARISING IN ELASTICITY
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 55, pp. 1 7. ISSN: 1072-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu GLOBAL WELL-POSEDNESS FO NONLINEA NONLOCAL
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 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 informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
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 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 informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationElastic Collisions Among Peakon Solutions for the Camassa-Holm Equation
Elastic Collisions Among Peakon Solutions for the Camassa-Holm Equation Alina Chertock, Jian-Guo Liu and Terrance Pendleton Abstract The purpose of this paper is to study the dynamics of the interaction
More informationCONVERGENCE OF THE LAX-FRIEDRICHS SCHEME AND STABILITY FOR CONSERVATION LAWS WITH A DISCONTINUOUS SPACE-TIME DEPENDENT FLUX
Chin. Ann. Math. 25B:3(24),287 318. CONVERGENCE OF THE LAX-FRIEDRICHS SCHEME AND STABILITY FOR CONSERVATION LAWS WITH A DISCONTINUOUS SPACE-TIME DEPENDENT FLUX K. H. KARLSEN J. D. TOWERS Abstract The authors
More informationGLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC SYSTEM
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 35, Number 4, April 7, Pages 7 7 S -99396)8773-9 Article electronically published on September 8, 6 GLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC
More informationBlow-up for a Nonlocal Nonlinear Diffusion Equation with Source
Revista Colombiana de Matemáticas Volumen 46(2121, páginas 1-13 Blow-up for a Nonlocal Nonlinear Diffusion Equation with Source Explosión para una ecuación no lineal de difusión no local con fuente Mauricio
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 informationANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL
Electronic Journal of Differential Equations, Vol. 217 (217), No. 289, pp. 1 6. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS
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 informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationWELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE
WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper
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 informationA Class of Shock Wave Solutions of the Periodic Degasperis-Procesi Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(2010) No.3,pp.367-373 A Class of Shock Wave Solutions of the Periodic Degasperis-Procesi Equation Caixia Shen
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 informationA CAUCHY-KOVALEVSKY THEOREM FOR NONLINEAR AND NONLOCAL EQUATIONS. In memory of M. Salah Baouendi
A CAUCHY-KOVALEVSKY THEOREM FOR NONLINEAR AND NONLOCAL EQUATIONS RAFAEL F. BAROSTICHI, A. ALEXANDROU HIMONAS* & GERSON PETRONILHO Abstract. For a generalized Camassa-Holm equation it is shown that the
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 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 informationExponential Energy Decay for the Kadomtsev-Petviashvili (KP-II) equation
São Paulo Journal of Mathematical Sciences 5, (11), 135 148 Exponential Energy Decay for the Kadomtsev-Petviashvili (KP-II) equation Diogo A. Gomes Department of Mathematics, CAMGSD, IST 149 1 Av. Rovisco
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
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 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 informationUNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS
Electronic Journal of Differential Equations, Vol. 16 16, No. 7, pp. 1 11. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIFORM DECAY OF SOLUTIONS
More informationThe Riemann problem for the Leray-Burgers equation
The iemann problem for the Leray-Burgers equation H. S. Bhat. C. Fetecau December 29, 28 Abstract For iemann data consisting of a single decreasing jump, we find that the Leray regularization captures
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 informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationMATH 6337: Homework 8 Solutions
6.1. MATH 6337: Homework 8 Solutions (a) Let be a measurable subset of 2 such that for almost every x, {y : (x, y) } has -measure zero. Show that has measure zero and that for almost every y, {x : (x,
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More informationarxiv: v2 [math.ap] 23 Feb 2012
Analytic solutions and Singularity formation for the Peakon b Family equations arxiv:.4983v [math.ap] 3 Feb G. M. Coclite, F. Gargano, and V. Sciacca Abstract. Using the Abstract Cauchy-Kowalewski Theorem
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationCUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION
CUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION HANS CHRISTIANSON Abstract. This paper shows how abstract resolvent estimates imply local smoothing for solutions to the Schrödinger equation.
More informationGlobal conservative solutions of the Camassa-Holm equation
Global conservative solutions of the Camassa-Holm equation Alberto Bressan Deptartment of Mathematics, Pennsylvania State University, University Park 168, U.S.A. e-mail: bressan@math.psu.edu and Adrian
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 informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationarxiv: v1 [math-ph] 17 Sep 2008
arxiv:080986v [math-ph] 7 Sep 008 New solutions for the modified generalized Degasperis Procesi equation Alvaro H Salas Department of Mathematics Universidad de Caldas Manizales Colombia Universidad Nacional
More informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
More informationAsymptotic behavior of the degenerate p Laplacian equation on bounded domains
Asymptotic behavior of the degenerate p Laplacian equation on bounded domains Diana Stan Instituto de Ciencias Matematicas (CSIC), Madrid, Spain UAM, September 19, 2011 Diana Stan (ICMAT & UAM) Nonlinear
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
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 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 informationHARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx.
Electronic Journal of Differential Equations, Vol. 003(003), No. 3, pp. 1 8. ISSN: 107-6691. UL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) HADY INEQUALITIES
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More 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 informationEXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS
Electronic Journal of Differential Equations, Vol. 2008(2008), No. 98, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) EXISTENCE
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 informationPublications of Giuseppe Maria Coclite
Publications of Giuseppe Maria Coclite Thesis [1] G. M. Coclite. Metodi Variazionali Applicati allo Studio delle Equazioni di Schrödinger-Maxwell. Graduate thesis, University of Bari (1999), supervisors:
More informationGIOVANNI COMI AND MONICA TORRES
ONE-SIDED APPROXIMATION OF SETS OF FINITE PERIMETER GIOVANNI COMI AND MONICA TORRES Abstract. In this note we present a new proof of a one-sided approximation of sets of finite perimeter introduced in
More informationDETERMINATION OF THE BLOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION
DETERMINATION OF THE LOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION y FRANK MERLE and HATEM ZAAG Abstract. In this paper, we find the optimal blow-up rate for the semilinear wave equation with a power nonlinearity.
More informationUNIQUENESS OF SELF-SIMILAR VERY SINGULAR SOLUTION FOR NON-NEWTONIAN POLYTROPIC FILTRATION EQUATIONS WITH GRADIENT ABSORPTION
Electronic Journal of Differential Equations, Vol. 2015 2015), No. 83, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIQUENESS OF SELF-SIMILAR
More informationSkr. K. Nor. Vidensk. Selsk., 2003, no. 3, 49 pp
STABILITY FOR ENTROPY SOLUTIONS OF NONLINEAR DEGENERATE PARABOLIC CONVECTION-DIFFUSION EQUATIONS WITH DISCONTINUOUS COEFFICIENTS L 1 K. H. KARLSEN, N. H. RISEBRO, AND J. D. TOWERS Skr. K. Nor. Vidensk.
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationOn uniqueness of weak solutions to transport equation with non-smooth velocity field
On uniqueness of weak solutions to transport equation with non-smooth velocity field Paolo Bonicatto Abstract Given a bounded, autonomous vector field b: R d R d, we study the uniqueness of bounded solutions
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationOSCILLATION-INDUCED BLOW-UP TO THE MODIFIED CAMASSA HOLM EQUATION WITH LINEAR DISPERSION
OSCILLATION-INDUCED BLOW-UP TO THE MODIFIED CAMASSA HOLM EQUATION WITH LINEAR DISPERSION ROBIN MING CHEN, YUE LIU, CHANGZHENG QU, AND SHUANGHU ZHANG Abstract. In this paper, we provide a blow-up mechanism
More informationA study of the modelling error in two operator splitting algorithms for porous media flow
Computational Geosciences 2 (1998) 23 36 23 A study of the modelling error in two operator splitting algorithms for porous media flow K. Brusdal a,h.k.dahle a, K. Hvistendahl Karlsen a and T. Mannseth
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 informationCONVERGENCE OF RELAXATION SCHEMES FOR HYPERBOLIC CONSERVATION LAWS WITH STIFF SOURCE TERMS
MATHEMATICS OF COMPUTATION Volume 68, Number 227, Pages 955 970 S 0025-5718(99)01089-3 Article electronically published on February 10, 1999 CONVERGENCE OF RELAXATION SCHEMES FOR HYPERBOLIC CONSERVATION
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 253 Monotonicity of Implicit Finite Difference Methods for Hyperbolic Conservation Laws Michael
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
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 informationOn universality of critical behaviour in Hamiltonian PDEs
Riemann - Hilbert Problems, Integrability and Asymptotics Trieste, September 23, 2005 On universality of critical behaviour in Hamiltonian PDEs Boris DUBROVIN SISSA (Trieste) 1 Main subject: Hamiltonian
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 informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationThe Cauchy problem for nonlocal evolution equations
The Cauchy problem for nonlocal evolution equations Alex Himonas Department of Mathematics Mini-workshop non-local dispersive equations NTNU, Trondheim, Norway September 21 23, 2015 Abstract We shall consider
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 informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationOn the distributional divergence of vector fields vanishing at infinity
Proceedings of the Royal Society of Edinburgh, 141A, 65 76, 2011 On the distributional divergence of vector fields vanishing at infinity Thierry De Pauw Institut de Recherches en Mathématiques et Physique,
More informationhal , version 1-22 Nov 2009
Author manuscript, published in "Kinet. Relat. Models 1, 3 8) 355-368" PROPAGATION OF GEVREY REGULARITY FOR SOLUTIONS OF LANDAU EQUATIONS HUA CHEN, WEI-XI LI AND CHAO-JIANG XU Abstract. By using the energy-type
More information