Research article Special Issues

Nonlocal scalar conservation laws with discontinuous flux

  • We prove the well-posedness of entropy weak solutions for a class of space-discontinuous scalar conservation laws with nonlocal flux. We approximate the problem adding a viscosity term and we provide L and BV estimates for the approximate solutions. We use the doubling of variable technique to prove the stability with respect to the initial data from the entropy condition.

    Citation: Felisia Angela Chiarello, Giuseppe Maria Coclite. Nonlocal scalar conservation laws with discontinuous flux[J]. Networks and Heterogeneous Media, 2023, 18(1): 380-398. doi: 10.3934/nhm.2023015

    Related Papers:

    [1] Darko Mitrovic . Existence and stability of a multidimensional scalar conservation law with discontinuous flux. Networks and Heterogeneous Media, 2010, 5(1): 163-188. doi: 10.3934/nhm.2010.5.163
    [2] Raimund Bürger, Harold Deivi Contreras, Luis Miguel Villada . A Hilliges-Weidlich-type scheme for a one-dimensional scalar conservation law with nonlocal flux. Networks and Heterogeneous Media, 2023, 18(2): 664-693. doi: 10.3934/nhm.2023029
    [3] Mauro Garavello, Roberto Natalini, Benedetto Piccoli, Andrea Terracina . Conservation laws with discontinuous flux. Networks and Heterogeneous Media, 2007, 2(1): 159-179. doi: 10.3934/nhm.2007.2.159
    [4] . Adimurthi, Siddhartha Mishra, G.D. Veerappa Gowda . Existence and stability of entropy solutions for a conservation law with discontinuous non-convex fluxes. Networks and Heterogeneous Media, 2007, 2(1): 127-157. doi: 10.3934/nhm.2007.2.127
    [5] Shyam Sundar Ghoshal . BV regularity near the interface for nonuniform convex discontinuous flux. Networks and Heterogeneous Media, 2016, 11(2): 331-348. doi: 10.3934/nhm.2016.11.331
    [6] Clément Cancès . On the effects of discontinuous capillarities for immiscible two-phase flows in porous media made of several rock-types. Networks and Heterogeneous Media, 2010, 5(3): 635-647. doi: 10.3934/nhm.2010.5.635
    [7] Jan Friedrich, Sanjibanee Sudha, Samala Rathan . Numerical schemes for a class of nonlocal conservation laws: a general approach. Networks and Heterogeneous Media, 2023, 18(3): 1335-1354. doi: 10.3934/nhm.2023058
    [8] Shyam Sundar Ghoshal, John D. Towers, Ganesh Vaidya . BV regularity of the adapted entropy solutions for conservation laws with infinitely many spatial discontinuities. Networks and Heterogeneous Media, 2024, 19(1): 196-213. doi: 10.3934/nhm.2024009
    [9] Raimund Bürger, Stefan Diehl, María Carmen Martí . A conservation law with multiply discontinuous flux modelling a flotation column. Networks and Heterogeneous Media, 2018, 13(2): 339-371. doi: 10.3934/nhm.2018015
    [10] Boris Andreianov, Kenneth H. Karlsen, Nils H. Risebro . On vanishing viscosity approximation of conservation laws with discontinuous flux. Networks and Heterogeneous Media, 2010, 5(3): 617-633. doi: 10.3934/nhm.2010.5.617
  • We prove the well-posedness of entropy weak solutions for a class of space-discontinuous scalar conservation laws with nonlocal flux. We approximate the problem adding a viscosity term and we provide L and BV estimates for the approximate solutions. We use the doubling of variable technique to prove the stability with respect to the initial data from the entropy condition.



    The aim of this paper is to study nonlocal conservation laws characterized by a flux discontinuous in space. In particular, the nonlocality consists in the fact that the velocity function depends on a convolution term that averages the solution in space. It is worth pointing out that the discontinuity appears in the flux through a multiplicative way. We will focus on the following equation,

    tρ+x(ρ(1wηρ)v(x))=0,(t,x)(0,)×R,

    where the function v=v(x) is defined as follows:

    v(x)={vl,if x<0,vr,if x>0.

    The idea comes from the work in [27] in which traveling waves are studied for a nonlocal scalar space discontinuous traffic model that describes the beaviour of drivers on two consecutive roads with different speed limits. Indeed, in recent years nonlocal conservation laws have been provided to describe several phenomena, for example: flux granular flows [2], sedimentation [6], supply chains [19], conveyor belts [18], structured populations dynamics [26] and traffic flows [7,9,10,16,28]. For these reasons, we believe the matter of discontinuous nonlocal conservation laws mathematically challenging and interesting while applicable to different real-life scenarios. Here, we prove the wellposedness of a nonlocal space discontinuous problem and our approach is based on a viscous regularizing approximation of the problem and standard compactness estimates. To our knowledge, these are the first results regarding discontinuous nonlocal problem using the vanishing viscosity technique. In particular, we have been inspired by the adaptation of the classical vanishing viscosity argument for scalar conservation laws [23] to erosion models [15], scalar equations with discontinous fluxes [8,24,25] and triangular systems [14]. This technique is based on the approximation of the solution of the starting problem through a sequence of smooth solutions of the corresponding viscous parabolic problem. The convergence to a solution of the starting problem is obtained proving compactness estimates on the sequence of smooth solutions. The existence of the approximate smooth solutions is proved through a fixed point theorem. In [9,10,27] conservation laws with nonlocal flux have been applied to the traffic flow setting. In particular, in [9,10] the authors study conservation laws with continuous flux functions and the well-posedness is obtained approximating the problem through an adapted numerical scheme and proving standard compactness estimates on the sequence of approximate solutions. In [27], travelling waves for a space-discontinuous traffic model describing two roads with rough conditions are studied. In the present work we do not need to apply an appropriate numerical discretisation of our problem due to the vanishing viscosity technique. We would like to count other more recent, noteworthy and interesting results about discontinuous nonlocal problems in [21] obtained with the fixed-point theorem technique. Our aim is to study a nonlocal equation in which the space-discontinuity occurs in the multiplicative term. It is not straightforward to deal with more general flux functions in the nonlocal setting satisfying the 'crossing condition' as in the paper [20]. Indeed, considering two different nonlocal flux functions for x<0 and x>0 would imply that the crossing point is not fixed but it changes position in time and this makes harder the analytical study. The paper is organized as follows. In Section 2, we describe the main results in this paper. In Section 3, we prove the existence of weak solutions of our problem, approximating it through a viscous problem and giving L and BV bounds. Finally, in Section 4, we show the uniqueness of entropy solutions, deriving an L1 contraction property using a doubling of variables argument.

    We consider the following scalar conservation equation with discontinuous nonlocal flux coupled with an initial datum

    {tρ+xf(t,x,ρ)=0,(t,x)(0,)×R,ρ(0,x)=ρ0(x),xR, (2.1)

    where

    f(t,x,ρ)=ρ(1wηρ)v(x),(wηρ)(t,x)=x+ηxρ(t,y)wη(yx)dy,η>0,

    and the velocity function v=v(x) is defined as follows

    v(x)={vl,if x<0,vr,if x>0.

    In this context ρ represents the unknown function, wη is a non-increasing kernel function whose length of the support is η. The equation in (2.1) is the space discontinuous version of the one in [7], where a nonlocal traffic model is presented.

    On wη,v,ρ0 we shall assume that

    0<vl<vr; (2.2)
    wηC2([0,η]),wη(η)=wη(η)=0,wη0wη,wηL1(0,η)=1; (2.3)
    0ρ01,ρ0L1(R)BV(R). (2.4)

    Assumption (2.3) implies that, if ρ is continuous,

    x(wηρ)(t,x)=(wηρ)(t,x)wη(0)ρ(t,x),x(wηρ)(t,x)=(wηρ)(t,x)wη(0)ρ(t,x). (2.5)

    Remark 2.1. The assumption Eq (2.3) does not allow the usual choices of kernels in traffic literature, such as: the constant and the linear decreasing kernels. Our kernels are like restrictions to [0,η] of cut-off functions that are equal to 1 in [η/2,η/2] and vanish outside [η,η] or, for example,

    wη(x):={e1x2η2x[0,η),0otherwise, (2.6)

    observing that limxηwη(x)=0.

    Remark 2.2. It is interesting to notice that if Eq (2.2) does not hold, namely vl>vr, we cannot say even in the local case that

    ρL((0,)×R)ρ0L(R).

    Let us consider this very easy example in the classical local case

    {tρ+xf(ρ)=0,(t,x)(0,)×(,0),tρ+xg(ρ)=0,(t,x)(0,)×(0,),ρ(0,x)=ρ0(x),xR, (2.7)

    where

    f(ρ)=2(ρ(1ρ)),g(ρ)=ρ(1ρ),ρ0(x)={0.25,if x<0,0.77,if x>0.

    The entropy weak solution to the above Cauchy problem is

    ρ={ρl=0.25,if x<f(ρ)f(ρl)ρρlt,ρ=0.9,if f(ρ)f(ρl)ρρlt<x<0,ρ+=ρr=0.77,if x>0. (2.8)
    Figure 1.  Fundamental diagrams relative to Eq (2.7).

    A complete description of conservation laws with discontinuous flux can be found in [17,22].

    We use the following definitions of solution.

    Definition 2.1. We say that a function ρ:[0,)×RR is a weak solution of Eq (2.1) if

    0ρ1,ρ(t,)L1(R)ρ0L1(R), (2.9)

    for almost every t>0 and for every test function φC1c(R2)

    0R(ρtφ+f(t,x,ρ)xφ)dtdx+Rρ0(x)φ(0,x)dx=0.

    Definition 2.2. A function ρ(L1L)(R+×R;[0,ρmax]) is an entropy weak solution of Eq (2.1), if

    (1) for all κR, and any test function φC1c(R2;R+) which vanishes for x0,

    +0R+|ρκ|φt+|ρκ|(1wηρ)vrφxsgn(ρκ)κx(wηρ)vrφdxdt+R+|ρ0(x)κ|φ(0,x)dx0;

    (2) for all κR, and any test function φC1c(R2;R+) which vanishes for x0,

    +0R|ρκ|φt+|ρκ|(1wηρ)vlφxsgn(ρκ)κx(wηρ)vlφdxdt+R|ρ0(x)κ|φ(0,x)dx0;

    (3) for all κR, and any test function φC1c(R2;R+)

    +0R|ρκ|φt+|ρκ|(1wηρ)v(x)φx+0Rsgn(ρκ)κx((wηρ)v(x))φdxdt+R|ρ0(x)κ|φ(0,x)dx++0|(vrvl)κ(1wηρ)|φ(t,0)dt0;

    (4) the traces are such that the jump

    |ρlρr| (2.10)

    is the smallest possible that satisfies the Rankine-Hugoniot condition

    f(t,0+,ρr)=f(t,0,ρl) i.e. vlρl=vrρr,

    where we denoted with

    f(t,0±,ρr,l)=limx0±v(x)ρ(t,x)(1x+ηxwη(yx)ρ(y,t)dy).

    Remark 2.3. We would like to underline that the existence of strong right and left traces, respectively ρr and ρl, is ensured by the genuine non-linearity of our flux function, as it is proved in [1,4].

    The main result of this paper is the following.

    Theorem 2.1. Assume Eq (2.2), Eq (2.3), and Eq (2.4). Then, the initial value problem in Eq (2.1) possesses an unique entropy solution u in the sense of Definition 2.2. Moreover, if u and v are two entropy solutions of Eq (2.1) in the sense of Definition 2.2, the following inequality holds

    u(t,)v(t,)L1(R)eKtu(0,)v(0,)L1(R), (2.11)

    for some suitable constant K>0.

    Our existence argument is based on passing to the limit in a vanishing viscosity approximation of Eq (2.1). We have been inspired by the viscous approximation in [11,Theorem 3.1], the sign of the term in absolute value |1wηρε| follows from Lemma 3.1 below.

    Fix a small number ε>0 and let ρε=ρε(t,x) be the unique classical solution of the following problem

    {tρε+(1wηρε)vε(x)xρε+ρε|1wηρε|vε(x)+ρε(wηρε)vε(x)+ρ2εwη(0)vε(x)=ε2xxρε,(t,x)(0,)×R,ρε(0,x)=ρ0,ε(x),xR, (3.1)

    where ρ0,ε and vε are C(R) approximations of ρ0 and v such that

    ρ0,ερ0,a.e. and in Lp(R),1p<,0ρ0,ε1,ρ0,εL1(R)ρ0L1(R),xρ0,εL1(R)C0vlvεvr,vε0,vε(x)={vlif x<ε,vrif x>ε, (3.2)

    for every ε>0 and some positive constant C0 independent on ε. The well-posedness of Eq (3.1) can be obtained following the same arguments of [11,12,13].

    Let us prove some a priori estimates on ρε denoting with C0 the constants which depend only on the initial data, and C(T) the constants which depend also on T.

    Lemma 3.1 (L estimate). Let ρε be a solution of (3.1). We have that

    0ρε1,

    for every ε>0.

    Proof. Thanks to Eq (3.2), 0 is a subsolution of Eq (3.1), due to the Maximum Principle for parabolic equations we have that

    ρε0. (3.3)

    We have to prove

    ρε1. (3.4)

    Assume by contradiction that Eq (3.4) does not hold.

    Let us define the function r(t,x)=eλtρε(t,x). We can choose λ>0 so small that

    rL((0,)×R)>1. (3.5)

    Thanks to Eq (3.1), r solves the equation

    tr+λr+(1wηρε)vε(x)xrε2xxr=r(wηρε)vε(x)r|1wηρε|vε(x)eλtr2wη(0)vε(x). (3.6)

    Since

    (wηρε(t,))(x)=x+ηxwη(yx)(ρε(t,y)ρεL((0,)×R))dyρεL((0,)×R)wη(0),

    we can write

    tr+λr+(1wηρε)vε(x)xrε2xxr=r(wη(ρερεL((0,)×R)))vε(x)+rwη(0)ρεL((0,)×R)vε(x)r|1wηρε|vε(x)eλtr2wη(0)vε(x)=r(wη(ρερεL((0,)×R)))vε(x)r|1wηρε|vε(x)+r(ρεL((0,)×R)ρε)wη(0)vε(x)0. (3.7)

    Let (ˉt,ˉx) be such that

    rL((0,)×R)=r(ˉt,ˉx).

    Since, thanks to Eq (3.5),

    r(0,)L(R)1<r(ˉt,ˉx),

    we must have

    ˉt>0.

    Therefore we can evaluate Eq (3.6) in (ˉt,ˉx) and gain

    0<λrL((0,)×R)0.

    Since, this cannot be, Eq (3.4) is proved.

    Using Eq (2.3) and Lemma 3.1, we know that

    0wηρε1 (3.8)

    and then we can rewrite Eq (3.1) as follows

    {tρε+x(ρε(1wηρε)vε(x))=ε2xxρε,(t,x)(0,)×R,ρε(0,x)=ρ0,ε(x),xR. (3.9)

    Lemma 3.2 (L1 estimate). Let ρε be a solution of Eq (3.1). We have that

    ρε(t,)L1(R)ρ0L1(R), (3.10)
    (wηρε)(t,)L1(R)ρ0L1(R), (3.11)
    x(wηρε)(t,)L1(R)2wη(0)ρ0L1(R), (3.12)

    for every t0 and ε>0.

    Proof. We have

    ddtRρεdx=Rtρεdx=εR2xxρεdxRx(ρε(1wηρε)vε(x))dx=0.

    Therefore,

    ρε(t,)L1(R)=ρ0,εL1(R),

    and Eq (3.10) follows from Eq (3.2).

    Using Eqs (2.3), (2.5), (3.8), and Lemma 3.1

    R(wηρε)(t,x)dx=Rx+ηxwη(yx)ρε(t,y)dydx=Rη0wη(y)ρε(t,y+x)dydx=wηL1(R)ρε(t,)L1(R)=ρε(t,)L1(R),R|x(wηρε)(t,x)|dxRx+ηx|wη(yx)|ρε(t,y)dxdy+wη(0)Rρεdx=Rη0wη(y)ρε(t,y+x)dxdy+wη(0)Rρεdx=2wη(0)ρε(t,)L1(R).

    Therefore, Eq (3.10), Eq (3.11), and Eq (3.12) follow from Eq (3.2).

    Lemma 3.3 (BV estimate in x). Let ρε be a solution of Eq (3.1). We have that

    xρε(t,)L1((,2δ)(2δ,))Cδ,

    for every t0 and ε,δ>0 where Cδ is a constant depending on δ but not on ε.

    Proof. Let us consider the function

    χ(x)={1,x(,2δ)(2δ,+),0,x(δ,δ),

    such that

    χC(R),0χ(x)1,χ(x)0 for x[0,+),χ(x)0 for x(,0].

    It is not restrictive to assume ε<δ. In such a way we have that the supports of χ and vε are disjoint. Finally, we observe that

    χχ,χχL(R).

    Differentiating the equation in (3.9) w.r.t. the space variable

    2txρε+x((1wηρε)vε(x)xρε)+x(ρε(1wηρε)vε(x))+x(ρε(wηρε)vε(x))+wη(0)x(ρ2εvε(x))=ε3xxxρε.

    Using [5,Lemma 2] and Lemmas 3.1, and 3.2

    ddtR|χ(x)xρε|dx=Rχ(x)2txρεsgnxρεdx=εRχ(x)3xxxρεsgnxρεdxRχ(x)x((1wηρε)vε(x)xρε)sgnxρεdxRχ(x)x(ρε|1wηρε|vε(x))sgnxρεdx=0Rχ(x)x(ρε(wηρε)vε(x))sgnxρεdxwη(0)Rχ(x)x(ρ2εvε(x))sgnxρεdx=εRχ(x)(2xxρε)2δ{xρε=0}0εRχ(x)2xxρεsgnxρε=x|xρε|dx+Rχ(x)(1wηρε)vε(x)xρε2xxρεδ{xρε=0}dx=0+Rχ(x)(1wηρε)vε(x)|xρε|dxRχ(x)(wηρε)vε(x)|xρε|dxRχ(x)(ρε(wηρε)vε(x)sgnxρεdx=0+Rχ(x)ρε(wηρε)vε(x)sgnxρεdx+wη(0)Rχ(x)ρ2εvε(x)sgnxρεdx2wη(0)Rχ(x)ρεvε(x)|xρε|dxwη(0)Rχ(x)ρ2εvε(x)sgnxρεdx=0cRχ(x)|xρε|dx+cRρεdxcRχ(x)|xρε|dx+cρ0L1(R),

    where δ{xρε=0} is the Dirac delta concentrated on the set {xρε=0} and c is a constant that depends on δ and does not depend on ε. Thanks to the Gronwall Lemma we get

    χxρε(t,)L1(R)ectχxρ0,εL1(R)+c(ect1),

    and using (3.2) we get the claim.

    Lemma 3.4 (Compactness). There exists a function ρ:[0,)×RR and a subsequence {εk}k(0,),εk0, such that

    0ρ1,ρBV((0,)×((,δ)(δ,))),δ>0,ρεkρa.e. and in Lploc((0,)×R),1p<.

    Proof. Thanks to Lemma 3.3 the sequence {ρεχIδ}ε,δ>0 of approximate solutions to Eq (2.1) constructed by vanishing viscosity has uniformly bounded variation on each interval of the type Iδ=(,δ)(δ,+),δ>0. Moreover, thanks to Lemma 3.1 the Lnorm of the sequence {ρεχIδ}ε,δ>0 is bounded by 1. Thus, applying Helly's Theorem and by a diagonal procedure, we can extract a subsequence {ρεkχIδk}kN that converges to a function ρ:[0,)×RR that satisfies the following conditions: ρBV((0,)×((,δ)(δ,))) and 0ρ1,

    ρεkχIδkρ a.e. and in Lploc((0,)×R),1p<.

    Thus, we obtain the compactness of the sequence {ρεk}kN a.e. in (0,)×R and for this reason we get the claim. It is worth remarking that being δ as small as we want we get the convergence on the whole space R.

    We are now ready to complete the proof of Theorem 2.1.

    Proof of Theorem 2.1. The existence of entropy solutions follows using the same arguments of [3] and Lemma 3.4. In particular, the nature of entropy solution of our limit function is related to the equivalence between [3,Definition 3] and [3,Definition 4] based on the germs theory, being our solution obtained through the vanishing viscosity technique. Moreover, one can observe that the points 1 and 2 of Definition 2.2 are directly satisfied multiplying equation (3.9) times the sgn(ρk), integrating with respect to time and space, and passing to the limit as ε0. The sketch of this proof is the following: we start from an L1 contraction property proved using the doubling of variables technique. After that we choose appropriate test functions in order to deal with the discontinuity in 0. We apply some limit procedures on the test functions and the classical Rankine-Hugoniot condition. At the end the Gronwall's inequality gives us the statement.

    Let us prove the inequality Eq (2.11). In Lemma 4.1 we prove the following inequality through the doubling of variables technique. For any two entropy solutions u and v we derive the L1 contraction property:

    R+×R(|uv|ϕt+sgn(uv)(f(t,x,u)f(t,x,v))ϕx)dxdtKR+×R|uv|ϕdxdt, (4.1)

    for any 0ϕCc(R+×R{0}). We remove the assumption in Eq (4.1) that ϕ vanishes near 0, by introducing the following Lipschitz function for h>0

    μh(x)={1h(x+2h),x[2h,h],1,x[h,h],1h(2hx),x[h,2h],0,|x|2h.

    Now we can define Ψh(x)=1μh(x), noticing that Ψh1 in L1 as h0. Moreover, Ψh vanishes in a neighborhood of 0. For any 0ΦCc(R+×R), we can check that ϕ=ΦΨh is an admissible test function for Eq (4.1). Using ϕ in Eq (4.1) and integrating by parts we get

    R+×R(|uv|ΦtΨh+sgn(uv)(f(t,x,u)f(t,x,v))ΦxΨh)dxdtR+×Rsgn(uv)(f(t,x,u)f(t,x,v))Φ(t,x)Ψh(x)dxdtJ(h)KR+×R|uv|ΦΨhdxdt.

    Sending h0 we end up with

    R+×R(|uv|Φt+sgn(uv)(f(t,x,u)f(t,x,v))Φx)dxdtKR+×R|uv|Φdxdt+limh0J(h).

    We can write

    limh0J(h)=limh01h+02hhsgn(uv)(f(t,x,u)f(t,x,v))Φ(t,x)dxdtlimh01h+0h2hsgn(uv)(f(t,x,u)f(t,x,v))Φ(t,x)dxdt=+0[sgn(uv)(f(t,x,u)f(t,x,v))]x=0+x=0Φ(t,0)dt,

    where we indicate the limits from the right and left at x=0. The aim is to prove that the limit limh0J(h)0, then it is sufficient to prove

    S:=[sgn(uv)(f(t,x,u)f(t,x,v))]x=0+x=00.

    In particular, denoting the right and left traces of u and v with u± and v±, we can write

    S=vrsgn(u+v+)(u+(1η0u(t,y)wη(y)dy)v+(1η0v(t,y)wη(y)dy))vlsgn(uv)(u(1η0u(t,y)wη(y)dy)v(1η0v(t,y)wη(y)dy))=vrsgn(u+v+)(v+u+)η0u(t,y)wη(y)dyvrsgn(u+v+)v+η0(u(t,y)v(t,y))wη(y)dy+vr|u+v+|vlsgn(uv)(vu)η0u(t,y)wη(y)dy+vlsgn(uv)vη0(u(t,y)v(t,y))wη(y)dyvl|uv|=(vr|u+v+|vl|uv|)=0(1η0u(y,t)wη(y)dy)+(vrv+vlv)=0sgn(uv)η0(v(t,y)u(t,y))wη(y)dy.

    A simple application of the Rankine-Hugoniot condition yields S=0, being u+=vlvru and v+=vlvrv.

    In this way we know that (4.1) holds for any 0ϕCc(R+×R). For r>1, let γr:RR be a C function which takes values in [0,1] and satisfies

    γr(x)={1,|x|r,0,|x|r+1.

    Fix s0 and s such that 0<s0<s. For any τ>0 and k>0 with 0<s0+τ<s+k, let βτ,k:[0,+]R be a Lipschitz function that is linear on [s0,s0+τ[[s,s+k] and satisfies

    βτ,k(t)={0,t[0,s0][s+k,+],1,t[s0+τ,s].

    We can take the admissible test function via a standard regularization argument ϕ=γr(x)βτ,k(t). Using this test function in Eq (4.1) we obtain

    1ks+ksR|u(t,x)v(t,x)|γr(x)dxdt1τs0+ks0R|u(t,x)v(t,x)|γr(x)dxdtKs0+ks0R|uv|γr(x)dxdt+γrs+ks0r|x|r+1sgn(uv)(f(t,x,u)f(t,x,v))dxdt.

    Sending s00, we get

    1ks+ksrr|u(t,x)v(t,x)|γr(x)dxdtrr|u0(x)v0(x)|dx+1ττ0rr|v(t,x)v0(x)|dxdt+1ττ0rr|u(t,x)u0(x)|dxdt+Kt+τ0R|uv|γr(x)dxdt+o(1r).

    Observe that the second and the third terms on the right-hand side of the inequality tends to zero as τ0 following the same argument in [20,Lemma B.1], because our initial condition is satisfied in the "weak" sense of the definition of our entropy condition. Sending τ0 and r, we have

    1ks+ksR|u(t,x)v(t,x)|dxdtR|u0(x)v0(x)|dx+Ks+k0R|u(t,x)v(t,x)|dxdt.

    Sending k0 and an application of Gronwall's inequality gives us the statement.

    Lemma 4.1 (A Kružkov-type integral inequality). For any two entropy solutions u=u(t,x) and v=v(t,x) the integral inequality of Eq (4.1) holds for any 0ϕCc(R+×R{0}).

    Proof. The proof follows [20]. Let 0ϕCc((R+×R{0})2),ϕ=ϕ(t,x,s,y),u=u(t,x) and v=v(s,y). From the definition of entropy solution for u=u(t,x) with κ=v(s,y) we get

    R+×R(|uv|ϕt+sgn(uv)(f(t,x,u)f(t,x,v))ϕx)dtdx+R+×R{0}sgn(uv)f(t,x,v)xϕdtdx0.

    Integrating over (s,y)R+×R, we find

    (R+×R)2(|uv|ϕt+sgn(uv)(f(t,x,u)f(t,x,v))ϕx)dtdxdsdy+(R+×R{0})2sgn(uv)f(t,x,v)xϕdtdxdsdy0. (4.2)

    Similarly, for the entropy solution v=v(s,y) with α(y)=u(t,x)

    (R+×R)2(|vu|ϕs+sgn(vu)(f(s,y,v)f(s,y,u))ϕx)dtdxdsdy+(R+×R{0})2sgn(uv)f(t,x,v)xϕdtdxdsdy0. (4.3)

    Note that we can write, for each (t,x)R+×R{0},

    sgn(uv)(f(t,x,u)f(t,x,v))ϕxsgn(uv)f(t,x,v)xϕ=sgn(uv)(f(t,x,u)f(s,y,v))ϕxsgn(uv)[(f(t,x,v)f(s,y,v))ϕ]x,

    so that

    (R+×R)2sgn(uv)(f(t,x,u)f(t,x,v))ϕxdtdxdsdy+(R+×R{0})2sgn(uv)f(t,x,v)xϕdtdxdsdy=(R+×R)2sgn(uv)(f(t,x,u)f(s,y,v))ϕxdtdxdsdy+(R+×R{0})2sgn(uv)[(f(t,x,v)f(s,y,v))ϕ]xdtdxdsdy.

    Similarly, writing, for each (y,s)R+×R{0}

    sgn(vu)(f(s,y,v)f(s,y,u))ϕysgn(vu)f(s,y,u)yϕ=sgn(uv)(f(s,y,v)f(s,y,u))ϕysgn(uv)[(f(t,x,u)f(s,y,u))ϕ]x,

    so that

    (R+×R)2sgn(uv)(f(s,y,v)f(s,y,u))ϕydtdxdsdy+(R+×R{0})2sgn(uv)f(s,y,u)yϕdtdxdsdy=(R+×R)2sgn(uv)(f(t,x,v)f(s,y,u))ϕxdtdxdsdy+(R+×R{0})2sgn(uv)[(f(t,x,u)f(s,y,u))ϕ]ydtdxdsdy.

    Let us introduce the notations

    t+s=t+s,x+y=x+y,2x+y=(x+y)2=2x+2xy+2y.

    Adding Eq (4.2) and Eq (4.3) we obtain

    (R+×R)2(|uv|t+sϕ+sgn(uv)(f(t,x,u)f(s,y,v))x+yϕ)dtdxdsdy+R+×R{0}sgn(uv)(x[(f(t,x,v)f(s,y,v))ϕ]+y[(f(t,x,u)f(s,y,u))ϕ])dtdxdsdy0. (4.4)

    We introduce a non-negative function δCc(R), satisfying δ(σ)=δ(σ),δ(σ)=0 for |σ|1, and Rδ(σ)dσ=1. For u>0 and zR, let δp(z)=1pδ(zp). We take our test function ϕ=ϕ(t,x,s,y) to be of the form

    Φ(t,x,s,y)=ϕ(t+s2,x+y2)δp(xy2)δp(ts2),

    where 0ϕCc(R+×R{0}) satisfies

    ϕ(t,x)=0,(t,x)[0,T]×[h,h],

    for small h>0. By making sure that

    p<h,

    one can check that Φ belongs to Cc((R+×R{0})2). We have

    t+sΦ(t,x,s,y)=t+sϕ(t+s2,x+y2)δp(xy2)δp(ts2),x+yΦ(t,x,s,y)=x+yϕ(t+s2,x+y2)δp(xy2)δp(ts2),

    Using Φ as test function in Eq (4.4)

    (R+×R)2(I1(t,x,s,y)+I2(t,x,s,y))δp(xy2)δp(ts2)dtdxdsdy(R+×R{0})2(I3(t,x,s,y)+I4(t,x,s,y)+I5(t,x,s,y))dtdxdsdy,

    where

    I1=|u(t,x)v(s,y)|t+sϕ(t+s2,x+y2),I2=sgn(u(t,x)v(s,y))(f(t,x,u)f(s,y,v))x+yϕ(t+s2,x+y2),I3=sgn(u(t,x)v(s,y))(xf(t,x,v)yf(s,y,u))ϕ(t+s2,x+y2,)δp(xy2)δp(ts2),I4=sgn(u(t,x)v(s,y))δp(xy2)δp(ts2)[xϕ(t+s2,x+y2,)(f(t,x,v)f(s,y,v))yϕ(t+s2,x+y2)(f(t,x,u)f(s,y,u))],I5=(F(x,u(t,x),v(s,y))F(y,u(t,x),v(s,y)))ϕ(t+s2,x+y2)xδp(xy2)δp(ts2),

    where F(x,u,c):=sgn(uc)(f(t,x,u)f(t,x,c)).

    We now use the change of variables

    ˜x=x+y2,˜t=t+s2,z=xy2,τ=ts2,

    which maps (R+×R)2 in ΩR4 and (R+×R{0})2 in Ω0R4, where

    Ω={(˜x,˜t,z,τ)R4:0<˜t±τ<T},Ω0={(˜x,˜t,z,τ)Ω:˜x±z0},

    resepectively. With this changes of variables,

    t+sϕ(t+s2,x+y2)=˜tϕ(˜t,˜x),x+yϕ(t+s2,x+y2)=˜xϕ(˜t,˜x).

    Now we can write

    Ω(I1(˜t,˜x,τ,z)+I2(˜t,˜x,τ,z))δp(z)δp(τ)d˜td˜xdτdzΩ0(I3(˜t,˜x,τ,z)+I4(˜t,˜x,τ,z)+I5(˜t,˜x,τ,z))d˜td˜xdτdz,

    where

    I1(˜t,˜x,τ,z)=|u(˜t+τ,˜x+z)v(˜tτ,˜xz)|˜tϕ(˜t,˜x),I2(˜t,˜x,τ,z)=sgn(u(˜t+τ,˜x+z)v(˜tτ,˜xz))(f(˜t+τ,˜x+z,u)f(˜tτ,˜xz,v))˜xϕ(˜t,˜x),I3(˜t,˜x,τ,z)=sgn(u(˜t+τ,˜x+z)v(˜tτ,˜xz))(˜x+zf(˜t+τ,˜x+z,v)˜xzf(˜tτ,˜xz,u))ϕ(˜t,˜x)δp(z)δp(τ),I4(˜t,˜x,τ,z)=sgn(u(˜t+τ,˜x+z)v(˜tτ,˜xz))˜xϕ(˜t,˜x)δp(z)δp(τ)[(f(˜t+τ,˜x+z,v)f(˜tτ,˜xz,v))+(f(˜t+τ,˜x+z,u)f(˜tτ,˜xz,u))],I5(˜t,˜x,τ,z)=(F(˜x+z,u(˜t+τ,˜x+z),v(˜tτ,˜xz))F(˜xz,u(˜t+τ,˜x+z),v(˜tτ,˜xz)))ϕ(˜t,˜x)zδp(z)δp(τ).

    Employing Lebesgue's differentiation theorem, to obtain the following limits

    limp0ΩI1(˜t,˜x,τ,z)δp(z)δp(τ)d˜td˜xdτdz=R+×R|u(t,x)v(t,x)|tϕ(t,x)dtdx,limp0ΩI2(˜t,˜x,τ,z)δp(z)δp(τ)d˜td˜xdτdz=R+×Rsgn(u(t,x)v(t,x))(f(t,x,u)f(t,x,v))xϕ(t,x)dtdx.

    Let us consider the term I3. Note that I3(˜t,˜x,τ,z)=0, if ˜x[h,h], since then ϕ(˜t,˜x)=0 for any ˜t, or if |z|p. On the other hand, if ˜x[h,h], then ˜x±z<0 or ˜x±z>0, at least when |z|<p and p<h. Defining U(t,x)=1wηu and V(t,x)=1wηv, and sending p0:

    limp0Ω0I3(˜t,˜x,τ,z)d˜td˜xdτdz=R+×R{0}sgn(u(t,x)v(t,x))v(x)(vxVuxU)ϕ(t,x)dtdxvrxVR+×R{0}|uv|ϕ(t,x)dtdx+vrR+×R{0}|ρ||xVxU|dtdxK1R+×R{0}|uv|ϕ(t,x)dtdx.

    In fact,

    |xVxU|ωηu(t,)v(t,)L1+wη(0)(|uv|(t,x+η)+|uv|(t,x)).

    The term I4 converges to zero as p0. Finally, the term I5

    limp0Ω0I5(˜t,˜x,τ,z)d˜td˜xdτdzK2R+×R{0}|uv|ϕ(t,x)dtdx.

    In this paper we proved the well-posedness of a Cauchy problem characterized by a nonlocal conservation law with space-discontinuous flux using the vanishing viscosity technique. This kind of equations can be applied to describe different real phenomena, such as: traffic flow, sedimentation, conveyor belts and others. It is worth noticing that the discontinuity appears in a multiplicative way. For this reason, one can think to consider more general nonlocal flux functions satisfying proper 'crossing conditions' in a future work and to study nonlocal-to-local limit in this space-discontinuous setting.

    The authors are members of the Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). GMC has been partially supported by the Research Project of National Relevance "Multiscale Innovative Materials and Structures" granted by the Italian Ministry of Education, University and Research (MIUR Prin 2017, project code 2017J4EAYB and the Italian Ministry of Education, University and Research under the Programme Department of Excellence Legge 232/2016 (Grant No. CUP - D94I18000260001).

    The authors declare there is no conflict of interest.



    [1] J. Aleksić, D. Mitrović, Strong traces for averaged solutions of heterogeneous ultra-parabolic transport equations, J. Hyperbolic Differ. Equ., 10 (2013), 659–676. https://doi.org/10.1142/S0219891613500239 doi: 10.1142/S0219891613500239
    [2] S. W. Amadori Debora, An integro-differential conservation law arising in a model of granular flow, J. Hyperbolic Differ. Equ., 9 (2012), 105–131.
    [3] B. Andreianov, K. H. Karlsen, N. H. Risebro, On vanishing viscosity approximation of conservation laws with discontinuous flux, Netw. Heterog. Media, 5 (2010), 617–633. https://doi.org/10.3934/nhm.2010.5.617 doi: 10.3934/nhm.2010.5.617
    [4] B. Andreianov, D. Mitrović, Entropy conditions for scalar conservation laws with discontinuous flux revisited, Ann. Inst. H. Poincaré Anal. Non Linéaire, 32 (2015), 1307–1335. https://doi.org/10.1016/j.anihpc.2014.08.002 doi: 10.1016/j.anihpc.2014.08.002
    [5] C. Bardos, A. Y. le Roux, J. C. Nédélec, First order quasilinear equations with boundary conditions, Comm. Partial Differential Equations, 4 (1979), 1017–1034. https://doi.org/10.1080/03605307908820117 doi: 10.1080/03605307908820117
    [6] F. Betancourt, R. Bürger, K. H. Karlsen, E. M. Tory, On nonlocal conservation laws modelling sedimentation, Nonlinearity, 24 (2011), 855–885.
    [7] S. Blandin, P. Goatin, Well-posedness of a conservation law with non-local flux arising in traffic flow modeling, Numer. Math., 132 (2016), 217–241. http://dx.doi.org/10.1007/s00211-015-0717-6 doi: 10.1007/s00211-015-0717-6
    [8] A. Bressan, G. Guerra, W. Shen, Vanishing viscosity solutions for conservation laws with regulated flux, J. Differential Equations, 266 (2019), 312–351. https://doi.org/10.1016/j.jde.2018.07.044 doi: 10.1016/j.jde.2018.07.044
    [9] F. A. Chiarello, J. Friedrich, P. Goatin, S. Göttlich, O. Kolb, A non-local traffic flow model for 1-to-1 junctions, Eur. J. Appl. Math., 31 (2020), 1029–1049.
    [10] F. A. Chiarello, P. Goatin, Global entropy weak solutions for general non-local traffic flow models with anisotropic kernel, ESAIM: M2AN, 52 (2018), 163–180. https://doi.org/10.1051/m2an/2017066 doi: 10.1051/m2an/2017066
    [11] G. M. Coclite, M. M. Coclite, Conservation laws with singular nonlocal sources, J. Differential Equations, 250 (2011), 3831–3858. https://doi.org/10.1016/j.jde.2010.12.001 doi: 10.1016/j.jde.2010.12.001
    [12] G. M. Coclite, L. di Ruvo, Well-posedness of the Ostrovsky-Hunter equation under the combined effects of dissipation and short-wave dispersion, J. Evol. Equ., 16 (2016), 365–389. https://doi.org/10.1007/s00028-015-0306-2 doi: 10.1007/s00028-015-0306-2
    [13] G. M. Coclite, H. Holden, K. H. Karlsen, Wellposedness for a parabolic-elliptic system, Discrete Contin. Dyn. Syst., 13 (2005), 659–682. https://doi.org/10.3934/dcds.2005.13.659 doi: 10.3934/dcds.2005.13.659
    [14] G. M. Coclite, K. H. Karlsen, S. Mishra, N. H. Risebro, Convergence of vanishing viscosity approximations of 2×2 triangular systems of multi-dimensional conservation laws, Boll. Unione Mat. Ital., 2 (2009), 275–284.
    [15] G. M. Coclite, E. Jannelli, Well-posedness for a slow erosion model, J. Math. Anal. Appl., 456 (2017), 337–355. https://doi.org/10.1016/j.jmaa.2017.07.006 doi: 10.1016/j.jmaa.2017.07.006
    [16] R. M. Colombo, M. Garavello and M. Lécureux-Mercier, A class of nonlocal models for pedestrian traffic, Math. Models Methods Appl. Sci., 22 (2012), 1150023.
    [17] M. Garavello, R. Natalini, B. Piccoli and A. Terracina, Conservation laws with discontinuous flux, Netw. Heterog. Media, 2 (2007), 159–179. https://doi.org/10.3934/nhm.2007.2.159. doi: 10.3934/nhm.2007.2.159
    [18] S. Göttlich, S. Hoher, P. Schindler, V. Schleper, A. Verl, Modeling, simulation and validation of material flow on conveyor belts, Appl. Math. Model., 38 (2014), 3295–3313.
    [19] M. Gröschel, A. Keimer, G. Leugering, Z. Wang, Regularity theory and adjoint-based optimality conditions for a nonlinear transport equation with nonlocal velocity, SIAM J. Control Optim., 52 (2014), 2141–2163.
    [20] K. H. Karlsen, N. H. Risebro, J. D. Towers, L1 stability for entropy solutions of nonlinear degenerate parabolic convection-diffusion equations with discontinuous coefficients, Skr. K. Nor. Vidensk. Selsk., 3 (2003), 1–49.
    [21] A. Keimer, L. Pflug, Discontinuous nonlocal conservation laws and related discontinuous ODEs Existence, Uniqueness, Stability and Regularity, arXiv: 2110.10503, [Preprint], (2021), [cited 2022 Dec 28 ]. Available from: https://arXiv.org/pdf/2110.10503.pdf
    [22] C. Klingenberg, N. H. Risebro, Convex conservation laws with discontinuous coefficients. Existence, uniqueness and asymptotic behavior, Comm. Partial Differential Equations, 20 (1995), 1959–1990. https://doi.org/10.1080/03605309508821159 doi: 10.1080/03605309508821159
    [23] S. N. Kružkov, First order quasilinear equations with several independent variables, Mat. Sb. (N.S.), 81 (1970), 228–255.
    [24] E. Y. Panov, Erratum to: Existence and strong pre-compactness properties for entropy solutions of a first-order quasilinear equation with discontinuous flux, Arch. Ration. Mech. Anal., 196 (2010), 1077–1078. https://doi.org/10.1007/s00205-010-0303-0 doi: 10.1007/s00205-010-0303-0
    [25] E. Y. Panov, Existence and strong pre-compactness properties for entropy solutions of a first-order quasilinear equation with discontinuous flux, Arch. Ration. Mech. Anal., 195 (2010), 643–673. https://doi.org/10.1007/s00205-009-0217-x doi: 10.1007/s00205-009-0217-x
    [26] B. Perthame, Transport equations in biology, Basel: Springer, 2007.
    [27] W. Shen, Traveling waves for conservation laws with nonlocal flux for traffic flow on rough roads, Netw. Heterog. Media, 14 (2019), 709–732. https://doi.org/10.3934/nhm.2019028 doi: 10.3934/nhm.2019028
    [28] A. Sopasakis, M. A. Katsoulakis, Stochastic modeling and simulation of traffic flow: asymmetric single exclusion process with Arrhenius look-ahead dynamics, SIAM J. Appl. Math., 66 (2006), 921–944.
  • This article has been cited by:

    1. Alexander Keimer, Lukas Pflug, 2023, 15708659, 10.1016/bs.hna.2022.11.001
    2. F. A. Chiarello, J. Friedrich, S. Göttlich, A non-local traffic flow model for 1-to-1 junctions with buffer, 2024, 19, 1556-1801, 405, 10.3934/nhm.2024018
    3. Veerappa Gowda G. D., Sudarshan Kumar Kenettinkara, Nikhil Manoj, Convergence of a second-order scheme for non-local conservation laws, 2023, 57, 2822-7840, 3439, 10.1051/m2an/2023080
    4. Alexander Keimer, Lukas Pflug, Discontinuous nonlocal conservation laws and related discontinuous ODEs – Existence, Uniqueness, Stability and Regularity, 2023, 361, 1778-3569, 1723, 10.5802/crmath.490
    5. F. A. Chiarello, H. D. Contreras, L. M. Villada, Existence of entropy weak solutions for 1D non-local traffic models with space-discontinuous flux, 2023, 141, 0022-0833, 10.1007/s10665-023-10284-5
    6. Aekta Aggarwal, Ganesh Vaidya, Convergence of the numerical approximations and well-posedness: Nonlocal conservation laws with rough flux, 2024, 0025-5718, 10.1090/mcom/3976
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2093) PDF downloads(99) Cited by(6)

Figures and Tables

Figures(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog