Loading [MathJax]/jax/output/SVG/jax.js
Communication Special Issues

A role for the non-conserved N-terminal domain of the TATA-binding protein in the crosstalk between cell signaling pathways and steroid receptors

  • Received: 28 October 2014 Accepted: 13 March 2015 Published: 25 March 2015
  • Transcriptional induction by steroid receptors is coupled to cellular signal transduction pathways although,in general,the mechanisms governing these events are not well defined. Using TATA-binding protein (TBP) specificity mutants that recognize a TGTA box,we show that yeast TBP expressed in mammalian cells can support steroid-mediated gene induction to a similar degree as human TBP,however yeast TBP does not support the 8-Bromo-cAMP-mediated potentiation of glucocorticoid receptor (GR)-dependent transactivation. Chimeras between yeast and human TBP reveal that it is the non-conserved N-terminus of TBP that governs the potentiation of GR action. While the conserved core of TBP is sufficient for TATA-element binding and preinitiation complex formation,the role of the N-terminus has remained elusive. Our results suggest a role of the N-terminus of human TBP in coupling cell signaling events to steroid-mediated transcription,thereby establishing one of the few described functional roles of this polypeptide domain in a physiological process.

    Citation: James R. Lambert, Steven K. Nordeen. A role for the non-conserved N-terminal domain of the TATA-binding protein in the crosstalk between cell signaling pathways and steroid receptors[J]. AIMS Molecular Science, 2015, 2(2): 64-76. doi: 10.3934/molsci.2015.2.64

    Related Papers:

    [1] Abraham Sylla . Influence of a slow moving vehicle on traffic: Well-posedness and approximation for a mildly nonlocal model. Networks and Heterogeneous Media, 2021, 16(2): 221-256. doi: 10.3934/nhm.2021005
    [2] Paola Goatin, Chiara Daini, Maria Laura Delle Monache, Antonella Ferrara . Interacting moving bottlenecks in traffic flow. Networks and Heterogeneous Media, 2023, 18(2): 930-945. doi: 10.3934/nhm.2023040
    [3] Maria Laura Delle Monache, Paola Goatin . Stability estimates for scalar conservation laws with moving flux constraints. Networks and Heterogeneous Media, 2017, 12(2): 245-258. doi: 10.3934/nhm.2017010
    [4] Felisia Angela Chiarello, Giuseppe Maria Coclite . Nonlocal scalar conservation laws with discontinuous flux. Networks and Heterogeneous Media, 2023, 18(1): 380-398. doi: 10.3934/nhm.2023015
    [5] 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
    [6] Markus Musch, Ulrik Skre Fjordholm, Nils Henrik Risebro . Well-posedness theory for nonlinear scalar conservation laws on networks. Networks and Heterogeneous Media, 2022, 17(1): 101-128. doi: 10.3934/nhm.2021025
    [7] Gabriella Bretti, Roberto Natalini, Benedetto Piccoli . Numerical approximations of a traffic flow model on networks. Networks and Heterogeneous Media, 2006, 1(1): 57-84. doi: 10.3934/nhm.2006.1.57
    [8] Christophe Chalons, Paola Goatin, Nicolas Seguin . General constrained conservation laws. Application to pedestrian flow modeling. Networks and Heterogeneous Media, 2013, 8(2): 433-463. doi: 10.3934/nhm.2013.8.433
    [9] Jan Friedrich, Oliver Kolb, Simone Göttlich . A Godunov type scheme for a class of LWR traffic flow models with non-local flux. Networks and Heterogeneous Media, 2018, 13(4): 531-547. doi: 10.3934/nhm.2018024
    [10] Martin Gugat, Alexander Keimer, Günter Leugering, Zhiqiang Wang . Analysis of a system of nonlocal conservation laws for multi-commodity flow on networks. Networks and Heterogeneous Media, 2015, 10(4): 749-785. doi: 10.3934/nhm.2015.10.749
  • Transcriptional induction by steroid receptors is coupled to cellular signal transduction pathways although,in general,the mechanisms governing these events are not well defined. Using TATA-binding protein (TBP) specificity mutants that recognize a TGTA box,we show that yeast TBP expressed in mammalian cells can support steroid-mediated gene induction to a similar degree as human TBP,however yeast TBP does not support the 8-Bromo-cAMP-mediated potentiation of glucocorticoid receptor (GR)-dependent transactivation. Chimeras between yeast and human TBP reveal that it is the non-conserved N-terminus of TBP that governs the potentiation of GR action. While the conserved core of TBP is sufficient for TATA-element binding and preinitiation complex formation,the role of the N-terminus has remained elusive. Our results suggest a role of the N-terminus of human TBP in coupling cell signaling events to steroid-mediated transcription,thereby establishing one of the few described functional roles of this polypeptide domain in a physiological process.


    Delle Monache and Goatin developed in [20] a macroscopic model aiming at describing the situation in which a slow moving large vehicle – a bus for instance – reduces the road capacity and thus generates a moving bottleneck for the surrounding traffic flow. Their model is given by a Cauchy problem for Lightwill-Whitham-Richards scalar conservation law in one space dimension with local point constraint. The constraint is prescribed along the slow vehicle trajectory $ (y(t),t) $, the unknown $ y $ being coupled to the unknown $ \rho $ of the constrained LWR equation. Point constraints were introduced in [19,17] to account for localized in space phenomena that may occur at exits and which act as obstacles. The constraint in the model of [20] depends upon the slow vehicle speed $ \dot y $, where its position $ y $ verifies the following ODE

    $ ˙y(t)=ω(ρ(y(t)+,t)). $ (A)

    Above, $ \rho = \rho(x,t) \in [0,R] $ is the traffic density and $ \omega: [0,R] \to \mathbb{R}^{+} $ is a nonincreasing Lipschitz continuous function which links the traffic density to the slow vehicle velocity. Delle Monache and Goatin proved an existence result for their model in [20] with a wave-front tracking approach in the $ \mathbf{BV} $ framework. Adjustments to the result were recently brought by Liard and Piccoli in [28]. Despite the step forward made in [21], the uniqueness issue remained open for a time. Indeed, the appearance of the trace $ \rho(y(t)^{+},t) $ makes it fairly difficult to get a Lipschitz continuous dependency of the trajectory $ y = y(t) $ from the solution $ \rho = \rho(x,t) $. Nonetheless, a highly nontrivial uniqueness result was achieved by Liard and Piccoli in [27]. To describe the influence of a single vehicle on the traffic flow, the authors of [26] proposed a PDE-ODE coupled model without constraint on the flux for which they proposed in [9] two convergent schemes. In the present paper, we consider a modified model where the point constraint becomes nonlocal, making the velocity of the slow vehicle depend on the mean density evaluated in a small vicinity ahead the driver. More precisely, instead of A, we consider the relation

    $ ˙y(t)=ω(Rρ(x+y(t),t)μ(x)dx), $ (B)

    where $ \mu \in \mathbf{BV}( \mathbb{R} ; \mathbb{R}^{+}) $ is a weight function used to average the density. From the mathematical point of view, this choice makes the study of the new model easier. Indeed, the authors of [3,4,5] put forward techniques for full well-posedness analysis of similar models with nonlocal point constraints. From the modeling point of view, considering B makes sense for several reasons outlined in Section 3.5.

    The paper is organized as follows. Sections 2 and 3 are devoted to the proof of the well-posedness of the model. In Section 4 we introduce the numerical finite volume scheme and prove its convergence. An important step of the reasoning is to prove a $ \mathbf{BV} $ regularity for the approximate solutions. It serves both in the existence proof and it is central in the uniqueness argument. In that optic, the appendix is essential. Indeed, it is devoted to the proof of a $ \mathbf{BV} $ regularity for entropy solutions to a large class of limited flux models. Let us stress that we highlight the interest of the $ \mathbf{BV}_{\mathbf{loc}} $ discrete compactness technique of Towers [33] in the context of general discontinuous-flux problems. In the numerical section 5, first we perform numerical simulations to validate our model. Then we investigate both qualitatively and quantitatively the proximity between our model – in which we considered B – as $ \delta \to \mu_{0^+} $ and the model of [20] in which the authors considered A.

    Note that we find it convenient to study the problem in the bus frame, which means setting $ X = x-y(t) $ in the model of Delle Monache and Goatin in [20]. Keeping in mind what we said above about the nonlocal constraint, the problem we consider takes the following form:

    $ {tρ+x(F(˙y(t),ρ))=0R×(0,T)ρ(x,0)=ρ0(x+y0)xRF(˙y(t),ρ)|x=0Q(˙y(t))t(0,T)˙y(t)=ω(Rρ(x,t)μ(x)dx)t(0,T)y(0)=y0. $ (1)

    Above, $ \rho = \rho(x,t) $ denotes the traffic density, of which maximum attainable value is $ R > 0 $, and

    $ F(\dot y(t),\rho) = f(\rho) - \dot y(t)\rho $

    denotes the normal flux through the curve $ x = y(t) $. We assume that the flux function $ f : [0,R] \to \mathbb{R} $ is Lipschitz continuous and bell-shaped, which are commonly used assumptions in traffic dynamics:

    $ f(ρ)0, f(0)=f(R)=0, ! ¯ρ(0,R), f(ρ)(¯ρρ)>0 for a.e. ρ(0,R). $ (2)

    In [20], the authors chose the function $ {{Q(s) = \alpha \times \left(\frac{1-s}{2} \right)^2}} $ to prescribe the maximal flow allowed through a bottleneck located at $ x = 0 $. The parameter $ \alpha \in (0, 1) $ was giving the reduction rate of the road capacity due to the presence of the slow vehicle. We use the $ s $ variable to stress that the value of the constraint is a function of the speed of the slow vehicle. In the sequel the $ s $ variable will refer to quantities related to the slow vehicle velocity. Regarding the function $ Q $, we can allow for more general choices. Specifically,

    $ Q : \left[0,\|\omega\|_{\mathbf{L}^{{\infty}}} \right] \to \mathbb{R}^{+} $

    can be any Lipschitz continuous function. It is a well known fact that in general, the total variation of an entropy solution to a constraint Cauchy problem may increase (see [17,Section 2] for an example). However, this increase can be controlled if the constraint level does not reach the maximum level. A mild assumption on $ Q $ – see Assumption (17) below – will guarantee availability of $ \mathbf{BV} $ bounds, provided we suppose that

    $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}). $

    Throughout the paper, we denote by

    $ \Phi(a,b) = \text{sign}(a-b)(f(a) - f(b)) \quad \text{and} \quad \Phi_{\dot y(t)}(a,b) = \Phi(a,b) - \dot y(t) |a-b| $

    the entropy fluxes associated with the Kružkov entropy $ \rho \mapsto |\rho-\kappa| $, for all $ \kappa \in [0,R] $, see [25]. Following [20,17,6,15], we give the following definition of solution for Problem (1).

    Definition 2.1. A couple $ (\rho,y) $ with $ \rho \in \mathbf{L}^{{\infty}}( \mathbb{R} \times (0,T)) $ and $ y \in {\mathbf{W}^{{1},{\infty}}}((0,T)) $ is an admissible weak solution to (1) if

    (ⅰ) the following regularity is fulfilled:

    $ ρC([0,T];L1loc(R)); $ (3)

    (ⅱ) for all test functions $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times \mathbb{R}^{+}), \varphi \geq 0 $ and $ \kappa \in [0,R] $, the following entropy inequalities are verified for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττR|ρκ|tφ+Φ˙y(t)(ρ,κ)xφdxdt+R|ρ(x,τ)κ|φ(x,τ)dxR|ρ(x,τ)κ|φ(x,τ)dx+2ττR˙y(t)(κ,q(t))φ(0,t)dt0, $ (4)

    where

    $ \mathcal{R}_{\dot y(t)}(\kappa,q(t)) = F(\dot y(t),\kappa) - \min \left\{ F(\dot y(t),\kappa),q(t) \right\} \quad \text{and} \quad q(t) = Q(\dot y(t)); $

    (ⅲ) for all test functions $ \psi \in \mathbf{C}^{{\infty}}([0,T]), \psi \geq 0 $ and some given $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}) $ which verifies $ \varphi(0) = 1 $, the following weak constraint inequalities are verified for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττR+ρt(φψ)+F(˙y(t),ρ)x(φψ)dxdtR+ρ(x,τ)φ(x)ψ(τ)dx+R+ρ(x,τ)φ(x)ψ(τ)dxττq(t)ψ(t)dt; $ (5)

    (ⅳ) the following weak ODE formulation is verified for all $ t \in [0,T] $:

    $ y(t)=y0+t0ω(Rρ(x,s)μ(x)dx)ds. $ (6)

    Definition 2.2. We will call $ \mathbf{BV} $-regular solution any admissible weak solution $ (\rho,y) $ to the Problem (1) which also verifies

    $ \rho \in \mathbf{L}^{{\infty}}((0,T) ; \mathbf{BV}( \mathbb{R})). $

    Remark 1. It is more usual to formulate (4) with $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times [0,T)) $, $ \tau = 0 $ and $ \tau' = T $. The equivalence between the two formulations is due to the regularity (3).

    Remark 2. As it happens, the time-continuity regularity (3) is actually a consequence of inequalities (4). Indeed, we will use the result [12,Theorem 1.2] which states that if $ \Omega $ is an open subset of $ \mathbb{R} $ and if for all test functions $ \varphi \in \mathbf{C}_c^{{\infty}}(\Omega \times [0,T)), \varphi \geq 0 $ and $ \kappa \in [0,R] $, $ \rho $ satisfies the following entropy inequalities:

    $ \int_{0}^{T} \int_{\Omega} |\rho-\kappa| \partial_{t}\varphi + \Phi_{\dot y(t)}(\rho,\kappa) \partial_{x}\varphi \mathinner{\mathrm{d}{x}} \mathinner{\mathrm{d}{t}} + \int_{\Omega}|\rho_{0}(x)-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq 0, $

    then $ {{\rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}(\Omega))}} $. Moreover, since $ \rho $ is bounded and $ \overline{\Omega} \backslash \Omega $ has a Lebesgue measure $ 0 $, $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}(\overline{\Omega})) $. We will use this remark several times in the sequel of the paper, with $ \Omega = \mathbb{R}^{*} $.

    Remark 3. Any admissible weak solution $ (\rho,y) $ to Problem (1) is also a distributional solution to the conservation law in (1). Therefore, inequalities (5) imply the following ones for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττRρt(φψ)+F(˙y(t),ρ)x(φψ)dxdt+Rρ(x,τ)φ(x)ψ(τ)dxRρ(x,τ)φ(x)ψ(τ)dxττq(t)ψ(t)dt, $

    where $ \varphi $ and $ \psi $ are such as described in Definition 2.1 (ⅲ).

    The interest of weak formulations (5)-(6) for the flux constraint and for the ODE governing the slow vehicle lies in their stability with respect to $ \rho $. Formulation (4) – (6) is well suited for passage to the limit of a.e. convergent sequences of exact or approximate solutions.

    In this section, we prove stability with respect to the initial data and uniqueness for $ \mathbf{BV} $-regular solutions to Problem (1). We start with the

    Lemma 2.3. If $ (\rho,y) $ is an admissible weak solution to Problem (1), then $ \dot y \in {\mathbf{W}^{{1},{\infty}}}((0,T)) $. In particular, $ \dot y \in \mathbf{BV}([0,T]) $.

    Proof. Denote for all $ t \in [0,T] $,

    $ s(t) = \omega \left( \int_{ \mathbb{R}} \rho(x,t) \mu(x) \mathinner{\mathrm{d}{x}} \right). $

    Since $ \mu \in \mathbf{L}^{{1}}( \mathbb{R}) \cap \mathbf{L}^{{\infty}}( \mathbb{R}) $ and $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $, $ s $ is continuous on $ [0,T] $. By definition, $ y $ satisfies the weak ODE formulation (6). Consequently, for a.e. $ t \in (0,T) $, $ \dot y(t) = s(t) $. We are going to prove that $ s $ is Lipschitz continuous on $ [0,T] $, which will ensure that $ \dot y \in {\mathbf{W}^{{1},{\infty}}}((0,T)) $. Since $ \mu \in \mathbf{BV}( \mathbb{R}) $, there exists a sequence $ {{(\mu_{n})_{n \in \mathbb{N}} \subset \mathbf{BV}( \mathbb{R}) \cap \mathbf{C}_c^{{\infty}}( \mathbb{R})}} $ such that:

    $ \| \mu_{n} - \mu \|_{\mathbf{L}^{{1}}} {\ \underset{{n} \to {+\infty}}{\longrightarrow} \ } 0 \quad \text{and} \quad \mathbf{TV}(\mu_{n}) {\ \underset{{n} \to {+\infty}}{\longrightarrow} \ } \mathbf{TV}(\mu). $

    Introduce for all $ n \in \mathbb{N} $ and $ t \in [0,T] $, the function

    $ \xi_{n}(t) = \int_{ \mathbb{R}} \rho(x,t) \mu_{n}(x) \mathinner{\mathrm{d}{x}}. $

    Fix $ \psi \in \mathbf{C}_c^{{\infty}}((0,T)) $. Since $ \rho $ is a distributional solution to the conservation law in (1), we have for all $ n \in \mathbb{N} $,

    $ T0ξn(t)˙ψ(t)dt=T0Rρt(ψμn)dxdt=T0RF(˙y(t),ρ)x(ψμn)dxdt=T0(RF(˙y(t),ρ)μn(x)dx)ψ(t)dt, $

    which means that for all $ n \in \mathbb{N} $, $ \xi_{n} $ is differentiable in the weak sense, and that for a.e. $ t \in (0,T) $,

    $ \dot \xi_{n}(t) = \int_ \mathbb{R} F(\dot y(t),\rho) \mu_{n}'(x) \mathinner{\mathrm{d}{x}}. $

    In particular, since both the sequences $ (\|\mu_{n}\|_{\mathbf{L}^{{1}}})_{n} $ and $ ( \mathbf{TV}(\mu_{n}))_{n} $ are bounded – say by $ C > 0 $ – we also have for all $ n \in \mathbb{N} $,

    $ \|\xi_{n}\|_{\mathbf{L}^{{\infty}}} \leq RC \quad \text{and} \quad \|\dot \xi_{n}\|_{\mathbf{L}^{{\infty}}} \leq C(\|f\|_{\mathbf{L}^{{\infty}}} + \|\omega\|_{\mathbf{L}^{{\infty}}}R). $

    Therefore, the sequence $ (\xi_{n})_{n} $ is bounded in $ {\mathbf{W}^{{1},{\infty}}}((0,T)) $. Now, for all $ t,\tau \in [0,T] $ and $ n \in \mathbb{N} $, triangle inequality yields:

    $ |s(t)s(τ)|2ωLRμnμL1+ωL|R(ρ(x,t)ρ(x,τ))μn(x)dx|=2ωLRμnμL1+ωL|ξn(t)ξn(τ)|2ωLRμnμL1+CωL(fL+ωLR)K|tτ|. $

    Letting $ n \to + \infty $, we get that for all $ t,\tau \in [0,T] $, $ {{|s(t) - s(\tau)| \leq K|t - \tau|}} $, which proves that $ s $ is Lipschitz continuous on $ [0,T] $. The proof of the statement is completed.

    Before stating the uniqueness result, we make the following additional assumption:

    $ s[0,ωL], argmaxρ[0,R] F(s,ρ)>0. $ (7)

    This ensures that for all $ s \in [0,\|\omega\|_{\mathbf{L}^{{\infty}}}] $, the function $ F(s,\cdot) $ verifies the bell-shaped assumptions (53). For example, when considering the flux $ f(\rho) = \rho(R-\rho) $, (7) reduces to $ \|\omega\|_{\mathbf{L}^{{\infty}}} < R $, which only means that the maximum velocity of the slow vehicle is lesser than the maximum velocity of the cars.

    Theorem 2.4. Suppose that $ f $ satisfies (2) and (7). Fix $ \rho_{0}^{1}, \rho_{0}^{2} \in \mathbf{L}^{{1}}( \mathbb{R}; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $ and $ y_{0}^{1}, y_{0}^{2} \in \mathbb{R} $. We denote by $ (\rho^{1},y^{1}) $ a $ \mathbf{BV} $-regular solution to Problem (1) corresponding to initial data $ (\rho_{0}^{1}, y_{0}^{1}) $, and by $ (\rho^{2},y^{2}) $ an admissible weak solution with initial data $ (\rho_{0}^{2},y_{0}^{2}) $. Then there exist constants $ \alpha,\beta,\gamma > 0 $ such that

    $ for a.e. t(0,T), ρ1(t)ρ2(t)L1(|y10y20|TV(ρ10)+ρ10ρ20L1)exp(αt) $ (8)

    and

    $ t[0,T], |y1(t)y2(t)||y10y20|+(β|y10y20|+γρ10ρ20L1)(exp(αt)1). $ (9)

    In particular, Problem (1) admits at most one $ \mathbf{BV} $-regular solution.

    Proof. Since $ (\rho^1 ,y^1) $ is a $ \mathbf{BV} $-regular solution to Problem (1), there exists $ C \geq 0 $ such that

    $ \forall t \in [0,T], \ \mathbf{TV}(\rho^1 (t)) \leq C. $

    Lemma 2.3 ensures that $ \dot y^1, \dot y^2 \in \mathbf{BV}([0,T] ; \mathbb{R}^+) $. We can use result (54) to obtain that for a.e. $ t \in (0,T) $,

    $ ρ1(t)ρ2(t)L1|y10y20|TV(ρ10)+ρ10ρ20L1+(2QL+2R+C)t0|˙y1(s)˙y2(s)|ds. $ (10)

    Moreover, since for a.e. $ t \in (0,T) $,

    $ |\dot y^1 (t) - \dot y^2 (t)| \leq \|\omega'\|_{\mathbf{L}^{{\infty}}}\|\mu\|_{\mathbf{L}^{{\infty}}} \|\rho^1 (t) - \rho^2 (t)\|_{\mathbf{L}^{{1}}}, $

    Gronwall's lemma yields (8) with $ {{\alpha = \left( 2 \|Q'\|_{\mathbf{L}^{{\infty}}} + 2R + C \right) \|\omega'\|_{\mathbf{L}^{{\infty}}}\|\mu\|_{\mathbf{L}^{{\infty}}}}} $. Then for all $ t \in [0,T] $,

    $ |y1(t)y2(t)||y10y20|+t0|˙y1(s)˙y2(s)|ds|y10y20|+ωLμLt0ρ1(s)ρ2(s)L1ds|y10y20|+(β|y10y20|+γρ10ρ20L1)(exp(αt)1), $

    where

    $ \beta = \frac{ \mathbf{TV}(\rho_{0}^1)}{ 2 \|Q'\|_{\mathbf{L}^{{\infty}}} + 2R + C} \quad \text{and} \quad \gamma = \frac{1}{2 \|Q'\|_{\mathbf{L}^{{\infty}}} + 2R + C}. $

    The uniqueness of a $ \mathbf{BV} $-regular solution is then clear.

    Remark 4. Up to inequality (10), our proof was very much following the one of [21,Theorem 2.1]. However, the authors of [21] faced an issue to derive a Lipschitz stability estimate between the car densities and the slow vehicle velocities starting from

    $ |\omega\left(\rho^1 (0+,t)\right) - \omega\left(\rho^2 (0+,t)\right)|. $

    For us, due to the nonlocality of our problem, it was straightforward to obtain the bound

    $ \left|\omega \left(\int_{ \mathbb{R}} \rho^1 (x,t) \mu(x) \mathinner{\mathrm{d}{x}} \right) - \omega \left(\int_{ \mathbb{R}} \rho^2 (x,t) \mu(x) \mathinner{\mathrm{d}{x}} \right) \right| \leq \|\omega'\|_{\mathbf{L}^{{\infty}}}\|\mu\|_{\mathbf{L}^{{\infty}}} \|\rho^1 (t) - \rho^2 (t)\|_{\mathbf{L}^{{1}}}. $

    Remark 5. A noteworthy consequence of Theorem 2.4 is that existence of a $ \mathbf{BV} $-regular solution will ensure uniqueness of an admissible weak one.

    In [20], to prove existence for their problem, the authors took a wave-front tracking approach. We choose here to use a time-splitting technique. The main advantage of this technique is that it relies on a ready-to-use theory. More precisely, at each time step, we will deal with exact solutions to a conservation law with a flux constraint, which have now become standard, see [17,6,15].

    Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) $ and $ y_{0} \in \mathbb{R} $. Let $ \delta > 0 $ be a time step, $ N \in \mathbb{N} $ such that $ T \in [N \delta, (N+1)\delta) $ and denote for all $ n \in \{0, \ldots, N+1\} $, $ t^{n} = n \delta $. We initialize with

    $ \forall t \in \mathbb{R}, \ \rho^{0}(t) = \rho_{0}(\cdot +y_{0}) \quad \text{and} \quad \forall t \in [0,T], \ y^{0}(t) = y_{0}. $

    Fix $ n \in \{1, \ldots ,N+1\} $. First, we define for all $ t \in (t^{n-1},t^{n}] $,

    $ \sigma^{n}(t) = \omega \left( \int_{ \mathbb{R}} \rho^{n-1}(x,t - \delta) \mu(x) \mathinner{\mathrm{d}{x}} \right), \ s^{n} = \sigma^{n}(t^{n}) \ \text{and} \ q^{n} = Q(s^{n}). $

    Since both $ q^n $ and $ \rho^{n-1}(\cdot, t^{n-1}) $ are bounded, [6,Theorem 2.11] ensures the existence and uniqueness of a solution $ \rho^n \in \mathbf{L}^{{\infty}}( \mathbb{R} \times [t^{n-1},t^{n}]) $ to

    $ \left\{ tρ+x(F(sn,ρ))=0R×(tn1,tn)ρ(x,tn1)=ρn1(x,tn1)xRF(sn,ρ)|x=0qnt(tn1,tn), \right. $

    in the sense that $ \rho^n $ satisfies entropy/constraint inequalities analogous to (4)-(5) with suitable flux/constraint function and initial data, see Definition A.1. Taking also into account Remark 2, $ \rho^{n} \in \mathbf{C}^{{}}([t^{n-1},t^{n}] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $. We then define the following functions:

    $  ρδ(t)=ρ01R(t)+N+1n=1ρn(t)1(tn1,tn](t) σδ(t),qδ(t),sδ(t)=σn(t),qn,snif t(tn1,tn] yδ(t)=y0+t0σδ(u)du. $

    First, let us prove that $ (\rho_{\delta},y_{\delta}) $ solves an approximate version of Problem (1).

    Proposition 1. The couple $ (\rho_{\delta},y_{\delta}) $ is an admissible weak solution to

    $ {tρδ+x(F(sδ(t),ρδ))=0R×(0,T)ρδ(x,0)=ρ0(x+y0)xRF(sδ(t),ρδ)|x=0qδ(t)t(0,T)˙yδ(t)=ω(Rρδ(x,tδ)μ(x)dx)t(0,T)yδ(0)=y0, $ (11)

    in the sense that $ \rho_\delta \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $ and satisfies entropy/constraint inequalities analogous to (4)-(5) with flux $ F(s_{\delta}(\cdot), \cdot) $, constraint $ q_\delta $, and initial data $ \rho_0(\cdot + y_0) $; and $ y_{\delta} $ satisfies, instead of (6), the following weak ODE formulation:

    $ \forall t \in [0,T], \ y_{\delta}(t) = y_0 + \int_0^t \omega \left( \int_ \mathbb{R} \rho_\delta(x,s-\delta) \mu(x) \mathinner{\mathrm{d}{x}} \right) \mathinner{\mathrm{d}{s}}. $

    Proof. By construction, for all $ n \in \{1, \ldots, N+1\}, \ \rho^{n} \in \mathbf{C}^{{}}([t^{n-1},t^{n}] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $. Combining this with the "stop-and-restart" conditions $ \rho^{n}(\cdot,t^{n-1}) = \rho^{n-1}(\cdot,t^{n-1}) $, one ensures that $ \rho_{\delta} \in \mathbf{C}^{{}}([0, T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $. Let $ t \in [0,T] $ and $ n \in \{1, \ldots, N+1\} $ such that $ t \in [t^{n-1},t^{n}) $. Then,

    $ yδ(t)y0=n1k=1tktk1σk(s)ds+ttn1σn(s)ds=n1k=1tktk1ω(Rρk1(x,sδ)ρδ(x,sδ)μ(x)dx)ds+ttn1ω(Rρn1(x,sδ)ρδ(x,sδ)μ(x)dx)ds=t0ω(Rρδ(x,sδ)μ(x)dx)ds, $ (12)

    which proves that $ \dot y_\delta $ solves the ODE in (11) in the weak sense. Fix now $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times \mathbb{R}^{+}), \varphi \geq 0 $ and $ \kappa \in [0,R] $. By construction of the sequence $ ((\rho^{k},y^{k}))_{k} $, we have for all $ n,m \in \{ 0, \ldots, N+1\} $,

    $ tmtnR|ρδκ|tφ+Φsδ(t)(ρδ,κ)xφdxdt=mk=n+1tktk1R|ρkκ|tφ+Φsk(ρk,κ)xφdxdtmk=n+1{R|ρk(x,tk)κ|φ(x,tk)dxR|ρk(x,tk1)ρk1(x,tk1)κ|φ(x,tk1)dx2tktk1Rsk(κ,qk)φ(0,t)dt}=R|ρδ(x,tm)κ|φ(x,tm)dxR|ρδ(x,tn)κ|φ(x,tn)dx2tmtnRsδ(t)(κ,qδ(t))φ(0,t)dt. $

    It is then straightforward to prove that for all $ 0 \leq \tau < \tau' \leq T $,

    $ ττR|ρδκ|tφ+Φsδ(t)(ρδ,κ)xφ dxdt+R|ρδ(x,τ)κ|φ(x,τ)dxR|ρδ(x,τ)κ|φ(x,τ)dx+2ττRsδ(t)(κ,qδ(t))φ(0,t)dt0. $ (13)

    Proving that $ \rho_\delta $ satisfies constraint inequalities is very similar so we omit the details. One has to start from

    $ - \int_{\tau}^{\tau'} \int_{ \mathbb{R}^{+}} \rho_{\delta} \partial_{t} (\varphi \psi) + F(s_{\delta}(t),\rho_{\delta}) \partial_{x} (\varphi \psi) \ \mathinner{\mathrm{d}{x}}\mathinner{\mathrm{d}{t}} $

    and make use once again of the construction of the sequence $ ((\rho^{k},y^{k}))_{k} $ to obtain

    $ ττR+ρδt(φψ)+F(sδ(t),ρδ)x(φψ) dxdtR+ρδ(x,τ)φ(x)ψ(τ)dx+R+ρδ(x,τ)φ(x)ψ(τ)dxττqδ(t)ψ(t)dt. $ (14)

    This concludes the proof.

    Remark 6. Remark that we have for all $ \delta > 0 $,

    $ \|\sigma_{\delta}\|_{\mathbf{L}^{{\infty}}} \leq \|\omega\|_{\mathbf{L}^{{\infty}}} \quad \text{and} \quad \|y_{\delta}\|_{\mathbf{L}^{{\infty}}} \leq |y_{0}| + T \|\omega\|_{\mathbf{L}^{{\infty}}}. $

    This means that the sequence $ (y_{\delta})_{\delta} $ is bounded in $ {\mathbf{W}^{{1},{\infty}}}((0,T)) $. Then the compact embedding of $ {\mathbf{W}^{{1},{\infty}}}((0,T)) $ in $ \mathbf{C}^{{}}([0,T]) $ yields a subsequence of $ (y_{\delta})_{\delta} $, which we do not relabel, which converges uniformly on $ [0,T] $ to some $ y \in \mathbf{C}^{{}}([0,T]) $.

    At this point, we propose two ways to obtain compactness for the sequence $ (\rho_{\delta})_{\delta} $, which will lead to two existence results.

    Theorem 3.1. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0;R]) $ and $ y_{0} \in \mathbb{R} $. Suppose that $ f $ is Lipschitz continuous, satisfies (2)-(7) and the following nondegeneracy assumption

    $ for a.e. s(0,ωL), mes{ρ[0,R] | f(ρ)s=0}=0. $ (15)

    Then Problem (1) admits at least one admissible weak solution.

    Proof. Condition (15) combined with the obvious uniform $ \mathbf{L}^{{\infty}} $ bound

    $ \forall \delta > 0, \ \forall (x,t) \in \mathbb{R} \times [0,T], \ \rho_{\delta}(x,t) \in [0,R], $

    and the results proved by Panov in [30,31] ensure the existence of a subsequence – which we do not relabel – that converges in $ \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}^{*} \times (0,T)) $ to some $ \rho \in \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}^{*} \times (0,T)) $; and a further extraction yields the almost everywhere convergence on $ \mathbb{R} \times (0,T) $ and also the fact that $ \rho \in \mathbf{L}^{{\infty}}( \mathbb{R} \times (0,T) ; [0,R]) $. We now show that the couple $ (\rho,y) $ constructed above is an admissible weak solution to (1) in the sense of Definition 2.1.

    For all $ \delta > 0 $ and $ t \in [0,T] $,

    $ yδ(t)y0=t0ω(Rρδ(x,sδ)μ(x)dx)ds=tδδω(Rρδ(x,s)μ(x)dx)ds=t0ω(Rρδ(x,s)μ(x)dx)ds+(0δttδ)ω(Rρδ(x,s)μ(x)dx)ds. $

    The last term vanishes as $ \delta \to 0 $ since $ \omega $ is bounded. Then, Lebesgue theorem combined with the continuity of $ \omega $ gives, for all $ t \in [0,T] $,

    $ y_{\delta}(t) {\ \underset{{\delta} \to {0}}{\longrightarrow} \ } y_{0} + \int_{0}^{t} \omega \left( \int_{ \mathbb{R}} \rho(x,s) \mu(x) \mathinner{\mathrm{d}{x}} \right) \mathinner{\mathrm{d}{s}}. $

    This last quantity is also equal to $ y(t) $ due to the uniform convergence of $ (y_{\delta})_{\delta} $ to $ y $. This proves that $ y $ verifies (6). Now, we aim at passing to the limit in (13) and (14). With this in mind, we prove the a.e. convergence of the sequence $ (\sigma_{\delta})_{\delta} $ towards $ \dot y $. Since $ \mu \in \mathbf{BV}( \mathbb{R}) $, there exists a sequence of smooth functions $ (\mu_{n})_{n \in \mathbb{N}} \subset \mathbf{BV}( \mathbb{R}) \cap \mathbf{C}_c^{{\infty}}( \mathbb{R}) $ such that:

    $ \|\mu_n - \mu\|_{\mathbf{L}^{{1}}} {\ \underset{{n} \to {+\infty}}{\longrightarrow} \ } 0 \quad \text{and} \quad \mathbf{TV}(\mu_{n}) {\ \underset{{n} \to {+\infty}}{\longrightarrow} \ } \mathbf{TV}(\mu). $

    Introduce for every $ \delta > 0 $ and $ n \in \mathbb{N} $, the function

    $ \xi_{\delta}^{n}(t) = \int_{ \mathbb{R}} \rho_{\delta}(x,t) \mu_{n}(x) \mathinner{\mathrm{d}{x}}. $

    Since for all $ \delta > 0 $, $ \rho_{\delta} $ is a distributional solution to the conservation law in (11), one can show – following the proof of Lemma 2.3 for instance – that for every $ n \in \mathbb{N} $, $ \xi_{\delta}^{n} \in {\mathbf{W}^{{1},{\infty}}}((0,T)) $, and that for a.e. $ t \in (0,T) $,

    $ \dot \xi_{\delta}^{n}(t) = \int_{ \mathbb{R}} F(s_{\delta}(t),\rho_{\delta}) \mu_{n}'(x) \mathinner{\mathrm{d}{x}}. $

    Moreover, since both the sequences $ \left( \|\mu_{n}\|_{\mathbf{L}^{{1}}} \right)_{n} $ and $ ( \mathbf{TV}(\mu_{n}))_{n} $ are bounded, it is clear that $ \left(\xi_{\delta}^{n}\right)_{\delta,n} $ is uniformly bounded in $ {\mathbf{W}^{{1},{\infty}}}((0,T)) $, therefore so is $ \left(\omega (\xi_{\delta}^{n})\right)_{\delta,n} $. Consequently, for all $ n \in \mathbb{N} $, $ \delta > 0 $ and almost every $ t \in (0,T) $, triangle inequality yields:

    $ |σδ(t)ω(Rρ(x,t)μ(x)dx)|2ωLRμnμL1+δ supnN ω(ξnδ)W1,+ωL|R(ρδ(x,t)ρ(x,t))μ(x)dx|δ0n+0, $

    which proves that $ (\sigma_{\delta})_{\delta} $ converges a.e. on $ (0,T) $ to $ \dot y $. To prove the time-continuity regularity, we first apply inequality (13) with $ \tau = 0 $, $ \tau' = T $ (which is licit since $ \rho_\delta $ is continuous in time), $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}^{*} \times [0,T)), \varphi \geq 0 $ and $ \kappa \in [0,R] $:

    $ \int_{0}^{T} \int_{ \mathbb{R}} |\rho_{\delta}-\kappa| \partial_{t}\varphi + \Phi_{\sigma_{\delta}(t)}(\rho_{\delta},\kappa) \partial_{x} \varphi \ \mathinner{\mathrm{d}{x}}\mathinner{\mathrm{d}{t}} + \int_{ \mathbb{R}}|\rho_{0}(x+y_{0})-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq 0. $

    Then, we let $ \delta \to 0 $ to get

    $ \int_{0}^{T} \int_{ \mathbb{R}} |\rho-\kappa| \partial_{t}\varphi + \Phi_{\dot y(t)}(\rho,\kappa) \partial_{x} \varphi \mathinner{\mathrm{d}{x}}\mathinner{\mathrm{d}{t}} + \int_{ \mathbb{R}}|\rho_{0}(x+y_{0})-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq 0. $

    Consequently, $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $, see Remark 2. Finally, the a.e. convergences of $ (\sigma_{\delta})_{\delta} $ and $ (\rho_{\delta})_\delta $ to $ \dot y $ and $ \rho $, respectively, are enough to pass to the limit in (13). This ensures that for all test functions $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times \mathbb{R}^{+}), \varphi \geq 0 $ and $ \kappa \in [0,R] $, the following inequalities hold for a.e. $ 0 \leq \tau < \tau' \leq T $:

    $ ττR|ρκ|tφ+Φ˙y(t)(ρ,κ)xφ dxdt+R|ρ(x,τ)κ|φ(x,τ)dxR|ρ(x,τ)κ|φ(x,τ)dx+2ττR˙y(t)(κ,q(t))φ(0,t)dt0. $

    Observe that the expression in the left-hand side of the previous inequality is a continuous function of $ (\tau,\tau') $ which is almost everywhere greater than the continuous function $ 0 $. By continuity, this expression is everywhere greater than $ 0 $, which proves that $ \rho $ satisfies the entropy inequalities (4). Using similar arguments, we show that $ \rho $ satisfies the constraint inequalities (5). This proves the couple $ (\rho,y) $ is an admissible weak solution to Problem (1), and this concludes the proof.

    In this section, we proved an existence result for $ \mathbf{L}^{{\infty}} $ initial data, but we have no guarantee of uniqueness since a priori we have no information regarding the $ \mathbf{BV} $ regularity of such solutions.

    Assumption (15) ensures the compactness for sequences of entropy solutions to conservation laws with flux function $ F $. However, it prevents us from using flux functions with linear parts – which corresponds to constant traffic velocity for small densities – whereas such fundamental diagrams are often used in traffic modeling. The results of the next section will extend to this interesting case, under the extra $ \mathbf{BV} $ assumption on the data.

    To obtain compactness for $ (\rho_\delta)_\delta $, an alternative to the setting of Section 3.2 is to derive uniform $ \mathbf{BV} $ bounds.

    Theorem 3.2. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $ and $ y_{0} \in \mathbb{R} $. Suppose that $ f $ satisfies (2) and (7). Suppose also that

    $ s[0,ωL], F(s,)C1([0,R]{¯ρs}), $ (16)

    where $ \overline{\rho}_s = \underset{\rho \in [0,R]}{{\rm{argmax}}} \ F(s,\rho) $. Finally assume that $ Q $ satisfies the condition

    $ ε>0, s[0,ωL], Q(s)maxρ[0,R] F(s,ρ)ε. $ (17)

    Then Problem (1) admits a unique admissible weak solution, which is also $ \mathbf{BV} $-regular.

    Proof. Fix $ \delta > 0 $. Recall that $ (\rho_{\delta},y_{\delta}) $ is an admissible weak solution to (11). In particular, $ \rho_\delta $ is an admissible weak solution to the constrained conservation law in (11), in the sense of Definition A.1. It is clear from the splitting construction that for a.e. $ t \in (0,T) $,

    $ \sigma_{\delta}(t) = \omega \left( \int_ \mathbb{R} \rho_\delta (x,t-\delta) \mu(x) \mathinner{\mathrm{d}{x}} \right). $

    Following the steps of the proof of Lemma 2.3, we can show that for all $ \delta > 0 $, $ \sigma_{\delta} \in \mathbf{BV}([0,T]; \mathbb{R}^+) $. Even more than that, by doing so we show that the sequence $ \left( \mathbf{TV}(\sigma_{\delta}) \right)_{\delta} $ is bounded. Therefore, the sequence $ \left( \mathbf{TV}(s_{\delta}) \right)_{\delta} $ is bounded as well. Moreover, since $ Q $ verifies (17), all the hypotheses of Corollary 2 are fulfilled. Combining this with Remark 15, we get the existence of a constant $ \mathbf{C}^{{ \varepsilon}} = \mathbf{C}^{{ \varepsilon}} (\| \partial_s F \|_{\mathbf{L}^{{\infty}}}) $ such that for all $ t \in [0,T] $,

    $ TV(ρδ(t))TV(ρ0)+4R+Cε(TV(qδ)+TV(sδ))TV(ρ0)+4R+Cε(1+QL)TV(sδ). $ (18)

    Consequently for all $ t \in [0,T] $, the sequence $ (\rho_{\delta}(t))_{\delta} $ is bounded in $ \mathbf{BV}( \mathbb{R}) $. A classical analysis argument – see [24,Theorem A.8] – ensures the existence of $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $ such that

    $ \forall t \in [0,T], \ \rho_{\delta}(t) {\ \underset{{\delta} \to {0}}{\longrightarrow} \ } \rho(t) \ \text{in} \ \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}). $

    With this convergence, we can follow the proof of Theorem 3.1 to show that $ (\rho,y) $ is an admissible weak solution to (1). Then, when passing to the limit in (18), the lower semi-continuity of the $ \mathbf{BV} $ semi-norm ensures that $ (\rho,y) $ is also $ \mathbf{BV} $-regular. By Remark 5, it ensures uniqueness and concludes the proof.

    To end this section, we now study the stability of Problem (1) with respect to the weight function $ \mu $. More precisely, let $ \left(\mu^{\ell}\right)_{\ell} \subset \mathbf{BV}( \mathbb{R} ; \mathbb{R}^{+}) $ be a sequence of weight functions that converges to $ \mu $ in the weak $ \mathbf{L}^{{1}} $ sense:

    $ gL(R), Rg(x)μ(x)dx + Rg(x)μ(x)dx. $ (19)

    Let $ (y_{0}^{\ell})_{\ell} \subset \mathbb{R} $ be a sequence of real numbers that converges to some $ y_{0} $ and let $ (\rho_{0}^{\ell})_{\ell} \subset \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) $ be a sequence of initial data that converges to $ \rho_{0} $ in the strong $ \mathbf{L}^{{1}} $ sense. We suppose that the flux function $ f $ satisfies Assumptions (2)-(7)-(15). Theorem 3.1 allows us to define or all $ \ell \in \mathbb{N} $, the couple $ (\rho^{\ell},y^{\ell}) $ as an admissible weak solution to the problem

    $ \left\{ tρ+x(F(˙y(t),ρ))=0R×(0,T)ρ(x,0)=ρ0(x+y0)xRF(˙y(t),ρ)|x=0Q(˙y(t))t(0,T)˙y(t)=ω(Rρ(x,t)μ(x)dx)t(0,T)y(0)=y0. \right. $

    Remark 7. Using the same arguments as in Remark 6 and as in the proof of Theorem 3.1, we get that up to the extraction of a subsequence, $ (y^{\ell})_{\ell} $ converges uniformly on $ [0,T] $ to some $ y \in \mathbf{C}^{{}}([0,T]) $ and $ (\rho^{\ell})_{\ell} $ converges a.e. on $ \mathbb{R} \times (0,T) $ to some $ \rho \in \mathbf{L}^{{\infty}}( \mathbb{R} \times (0,T)) $.

    Theorem 3.3. The couple $ (\rho,y) $ constructed above is an admissible weak solution to Problem (1).

    Proof. The sequence $ (\mu^{\ell})_{\ell} $ converges in the weak $ \mathbf{L}^{{1}} $ sense and is bounded in $ \mathbf{L}^{{1}}( \mathbb{R}) $; by the Dunford-Pettis theorem, this sequence is equi-integrable:

    $ ε>0, α>0, AB(R), mes(A)<α N, Aμ(x)dxε $ (20)

    and

    $ ε>0, X>0, N, |x|Xμ(x)dxε. $ (21)

    Fix $ t \in (0,T) $ and $ \varepsilon > 0 $. Fix $ \alpha,X > 0 $ given by (20) and (21). Egoroff theorem yields the existence of a measurable subset $ E_{t} \subset [-X,X] $ such that

    $ \text{mes}([-X,X] \backslash E_{t}) < \alpha \ \text{and} \ \rho^{\ell}( \cdot, t) \longrightarrow \rho(\cdot,t) \ \text{uniformly on} \ E_{t}. $

    For a sufficiently large $ \ell \in \mathbb{N} $,

    $ |Rρ(x,t)μ(x)dxRρ(x,t)μ(x)dx||x|X|ρρ|μdx+|Et(ρρ)μdx|+|[X,X]Et(ρρ)μdx|+|RρμdxRρμdx|Rε+ρρL(Et)Etμ(x)dx+R[X,X]Etμ(x)dx+ε2(R+1)ε, $

    which proves that for a.e. $ t \in (0,T) $,

    $ Rρ(x,t)μ(x)dx + Rρ(x,t)μ(x)dx. $ (22)

    We get that $ y $ verifies the weak ODE formulation (6) by passing to the limit in

    $ y^{\ell}(t) = y_{0}^{\ell} + \int_{0}^{t} \omega \left( \int_{ \mathbb{R}} \rho^{\ell}(x,s) \mu^{\ell}(x) \mathinner{\mathrm{d}{x}} \right) \mathinner{\mathrm{d}{s}}. $

    By definition, for all $ \ell \in \mathbb{N} $, the couple $ (\rho^{\ell},y^{\ell}) $ satisfies the analogue of entropy/constraint inequalities (4)-(5) with suitable flux/constraint functions. Applying these inequalities with $ \tau = 0 $, $ \tau' = T $, $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}^{*} \times [0,T)), \varphi \geq 0 $ and $ \kappa \in [0,R] $, we get

    $ \int_{0}^{T} \int_{ \mathbb{R}} |\rho^{\ell}-\kappa| \partial_{t}\varphi + \Phi_{\dot y^{\ell}(t)}(\rho^{\ell},\kappa) \partial_{x}\varphi \mathinner{\mathrm{d}{x}} \mathinner{\mathrm{d}{t}} + \int_{ \mathbb{R}}|\rho_{0}^{\ell}(x+y_{0}^{\ell})-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq 0. $

    The continuity of $ \omega $ and the convergence (22) ensure that $ (\dot y^{\ell})_{\ell} $ converges a.e. to $ \dot y $. This combined with the a.e. convergence of $ (\rho^{\ell})_{\ell} $ to $ \rho $ and Riesz-Frechet-Kolmogorov theorem – $ \left( \rho_{0}^{\ell} \right)_{\ell} $ being strongly compact in $ \mathbf{L}^{{1}}( \mathbb{R}) $ – is enough to show that when letting $ \ell \to +\infty $ in the inequality above, we get, up to the extraction of a subsequence, that

    $ \int_{0}^{T} \int_{ \mathbb{R}} |\rho-\kappa| \partial_{t}\varphi + \Phi_{\dot y(t)}(\rho,\kappa) \partial_{x} \varphi \mathinner{\mathrm{d}{x}} \mathinner{\mathrm{d}{t}} + \int_{ \mathbb{R}}|\rho_{0}(x+y_{0})-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq 0. $

    Consequently $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $, see Remark 2. Finally, the combined a.e. convergences of $ (\dot y^{\ell})_{\ell} $ and $ (\rho^{\ell})_{\ell} $ to $ \dot y $ and $ \rho $, respectively, guarantee that $ (\rho,y) $ verifies inequalities (4)-(5) for almost every $ 0 \leq \tau < \tau' \leq T $. The same continuity argument we used in the proof Theorem 3.1 holds here to ensure that $ (\rho,y) $ actually satisfies the inequalities for all $ 0 \leq \tau < \tau' \leq T $. This concludes the proof of our stability claim.

    The last section concludes the theoretical analysis of Problem (1). The nonlocality in space of the constraint delivers an easy proof of stability with respect to the initial data in the $ \mathbf{BV} $ framework. Although a proof of existence using a fixed point theorem was possible (cf. [5]), we chose to propose a proof based on a time-splitting technique. The stability with respect to $ \mu $ is a noteworthy feature, which shows a certain sturdiness of the model. However, the case we had in mind – namely $ \mu \to \delta_{0^{+}} $ – is not reachable with the assumptions we used to prove the stability, especially (19). We will explore this singular limit numerically, after having built a robust convergent numerical scheme for Problem (1). Let us also underline that unlike in [27,28] where the authors required a particular form for the function $ \omega $ to prove well-posedness for their model, our result holds as long as $ \omega $ is Lipschitz continuous.

    As evoked earlier, the nonlocality in space of the constraint makes the mathematical study of the model easier. But in the modeling point of view, this choice also makes sense for several reasons. First of all, one can think that the velocity $ \dot y $ of the slow moving vehicle – unlike its acceleration – is a rather continuous value. Even if the driver of the slow vehicle suddenly applies the brakes, the vehicle will not decelerate instantaneously. Note that the LWR model allows for discontinuous averaged velocity of the agents, however while modeling the slow vehicle we are concerned with an individual agent and can model its behavior more precisely. Moreover, considering the mean value of the traffic density in a vicinity ahead of the driver could be seen at taking into account both the driver anticipation and a psychological effect. For example, if the driver sees – several dozens of meters ahead of him/her – a speed reduction on traffic, he/she will start to slow down. This observation can be related to the fact that, compared to the fluid mechanics models where the typical number of agents is governed by the Avogadro constant, in traffic models the number of agents is at least $ 10^{20} $ times less. Therefore, a mild nonlocality (evaluation of the downstream traffic flow via averaging over a handful of preceding cars) is a reasonable assumption in the macroscopic traffic models inspired by fluid mechanics. This point of view is exploited in the model of [16]. Note that it is feasible to substitute the basic LWR equation on $ \rho $ by the nonlocal LWR introduced in [16] in our nonlocal model for the slow vehicle. Such mildly nonlocal model remains close to the basic local model of [20]. It can be studied combining the techniques of [16] and the ones we developed in this section.

    In this section, we aim at constructing a finite volume scheme and at proving its convergence toward the $ \mathbf{BV} $-regular solution to (1). We will use the notations:

    $ a \vee b = \max\{a,b\} \quad \text{and} \quad a \wedge b = \min\{a,b\}. $

    Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) $ and $ y_{0} \in \mathbb{R} $.

    For a fixed spatial mesh size $ \Delta x $ and time mesh size $ \Delta t $, let $ x_{j} = j\Delta x $, $ t^{n} = n \Delta t $. We define the grid cells $ \mathbf{K}_{j+1/2} = (x_{j},x_{j+1}) $. Let $ N \in \mathbb{N} $ such that $ T \in [t^N, t^{N+1}) $. We write

    $ \mathbb{R} \times [0,T] \subset \bigcup\limits_{n = 0}^{N} \bigcup\limits_{j \in \mathbb{Z}} \mathcal{P}_{j+1/2}^{n}, \quad \mathcal{P}_{j+1/2}^{n} = \mathbf{K}_{j+1/2} \times [t^{n},t^{n+1}). $

    We choose to discretize the initial data $ \rho_{0}(\cdot + y_{0}) $ and the weight function $ \mu $ with $ \left(\rho_{j+1/2}^{0}\right)_{j \in \mathbb{Z}} $ and $ \left(\mu_{j+1/2}\right)_{j \in \mathbb{Z}} $ where for all $ j \in \mathbb{Z} $, $ \rho_{j+1/2}^{0} $ and $ \mu_{j+1/2} $ are their mean values on the cell $ \mathbf{K}_{j+1/2} $.

    Remark 8. Others choice could be made, for instance in the case $ \rho_{0} \in \mathbf{C}^{{}}( \mathbb{R}) $ such that $ {{\lim_{|x| \to +\infty} \rho_{0}(x)}} $ exists (in which case, the limit is zero due to the integrability assumption), the values $ \rho_{j+1/2}^{0} = \rho_{0}\left(x_{j+1/2} + y_{0}\right) $ can be used. The only requirements are

    $ \forall j \in \mathbb{Z}, \ \rho_{j+1/2}^{0} \in [0,R] \quad \text{and} \quad \rho_{\Delta}^0 = \sum\limits_{j \in \mathbb{Z}} \rho_{j+1/2}^{0} \mathbb{1}_{ \mathbf{K}_{j+1/2}} {\ \underset{{\Delta x} \to {0}}{\longrightarrow} \ } \rho_{0}(\cdot + y_{0}) \ \text{in} \ \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}). $

    Fix $ n \in \{0, \ldots, N-1\} $. At each time step we first define an approximate velocity of the slow vehicle $ s^{n+1} $ and a constraint level $ q^{n+1} $:

    $ sn+1=ω(jZρnj+1/2μj+1/2Δx),qn+1=Q(sn+1). $ (23)

    With these values, we update the approximate traffic density with the marching formula for all $ j \in \mathbb{Z} $:

    $ ρn+1j+1/2=ρnj+1/2ΔtΔx(Fn+1j+1(ρnj+1/2,ρnj+3/2)Fn+1j(ρnj1/2,ρnj+1/2)), $ (24)

    where, following the recipe of [6,15],

    $ Fn+1j(a,b)={Fn+1(a,b)if j0min{Fn+1(a,b),qn+1}if j=0, $ (25)

    $ \mathbf{F}^{n+1} $ being a monotone consistent and Lipschitz numerical flux associated to $ F(s^{n+1},\cdot) $. We will also use the notation

    $ ρn+1j+1/2=Hn+1j(ρnj1/2,ρnj+1/2,ρnj+3/2), $ (26)

    where $ \mathbf{H}_{j}^{n+1} $ is given by the expression in the right-hand side of (24). We then define the functions

    $  ρΔ(x,t)=Nn=0jZρnj+1/21Pnj+1/2(x,t) sΔ(t),qΔ(t)=sn+1,qn+1if t[tn,tn+1) yΔ(t)=y0+t0sΔ(u)du. $

    Let $ \Delta = (\Delta x, \Delta t) $. For our convergence analysis, we will assume that $ \Delta \to 0 $, with $ \lambda = \Delta t / \Delta x $ verifying the CFL condition

    $ λsups[0,ωL] (FsaL+FsbL)L1, $ (27)

    where $ \mathbf{F}^s = \mathbf{F}^s(a, b) $ is the numerical flux, associated to $ F(s,\cdot) $, we use in (24).

    Remark 9. When considering the Rusanov flux or the Godunov one, (27) is guaranteed when

    $ 2 \lambda ( \|f'\|_{\mathbf{L}^{{\infty}}} + \|\omega\|_{\mathbf{L}^{{\infty}}} ) \leq 1. $

    Proposition 2 ($ \mathbf{L}^{{\infty}} $ stability). The scheme (26) is

    (i) monotone: for all $ n \in \{0, \ldots, N-1\} $ and $ j \in \mathbb{Z} $, $ \mathbf{H}_{j}^{n+1} $ is nondecreasing with respect to its three arguments;

    (ii) stable:

    $ n{0,,N}, jZ,ρnj+1/2[0,R]. $ (28)

    Proof. (i) In the classical case – $ j \notin \{-1,0\} $ – we simply differentiate the Lipschitz function $ \mathbf{H}^{n+1}_{j} $ and make use of both the CFL condition (27) and the monotonicity of $ \mathbf{F}^{n+1} $. For $ j \in \{-1,0\} $, note that the authors of [6] pointed out (in Proposition 4.2) that the modification done in the numerical flux (25) does not change the monotonicity of the scheme.

    (ii) The $ \mathbf{L}^{{\infty}} $ stability is a consequence of the monotonicity and also of the fact that $ 0 $ and $ R $ are stationary solutions of the scheme. Indeed, as in [6,Proposition 4.2] for all $ n \in \{ 0, \ldots, N \} $ and $ j \in \mathbb{Z} $,

    $ \mathbf{H}_{j}^{n+1}(0,0,0) = 0, \quad \mathbf{H}_{j}^{n+1}(R,R,R) = R. $

    In order to show that the limit of $ (\rho_{\Delta})_{\Delta} $ – under the a.e. convergence up to a subsequence – is a solution of the conservation law in (1), we derive discrete entropy inequalities. These inequalities also contain terms that will help to pass to the limit in the constrained formulation of the conservation law, as soon as the sequence $ (q_{\Delta})_\Delta $ of constraints is proved convergent as well.

    Proposition 3 (Discrete entropy inequalities). The numerical scheme (26) fulfills the following inequalities for all $ n \in \{0, \ldots, N-1\} $, $ j \in \mathbb{Z} $ and $ \kappa \in [0,R] $:

    $ (|ρn+1j+1/2κ||ρnj+1/2κ|)Δx+(Φnj+1Φnj)ΔtRsn+1(κ,qn+1)Δt δj{1,0}+(Φn0Φnint)Δt(δj=1δj=0), $ (29)

    where $ \Phi_{j}^{n} $ and $ \Phi_{int}^{n} $ denote the numerical entropy fluxes:

    $ Φnj=Fn+1(ρnj1/2κ,ρnj+1/2κ)Fn+1(ρnj1/2κ,ρnj+1/2κ),Φnint=min{Fn+1(ρn1/2κ,ρn1/2κ),qn+1}min{Fn+1(ρn1/2κ,ρn1/2κ),qn+1} $

    and

    $ \mathcal{R}_{s^{n+1}}(\kappa,q^{n+1}) = F(s^{n+1},\kappa) - \min\{ F(s^{n+1},\kappa),q^{n+1}\}. $

    Proof. This result is a direct consequence of the scheme monotonicity. When the constraint does not enter the calculations i.e. $ j \notin \{-1,0\} $, the proof follows [23,Lemma 5.4]. The key point is not only the monotonicity, but also the fact that in the classical case, all the constants $ \kappa \in [0,R] $ are stationary solutions of the scheme. This observation does not hold when the constraint enters the calculations. For example if $ j = -1 $,

    $ \mathbf{H}_{-1}^{n+1}(\kappa,\kappa,\kappa) = \kappa + \lambda \mathcal{R}_{s^{n+1}}(\kappa,q^{n+1}). $

    Consequently, we have both

    $ \rho_{-1/2}^{n+1} \vee \kappa \leq \mathbf{H}_{-1}^{n+1}(\rho_{-3/2}^{n} \vee \kappa,\rho_{-1/2}^{n} \vee \kappa,\rho_{1/2}^{n} \vee \kappa) $

    and

    $ \rho_{-1/2}^{n+1} \wedge \kappa \geq \mathbf{H}_{-1}^{n+1}(\rho_{-3/2}^{n} \wedge \kappa,\rho_{-1/2}^{n} \wedge \kappa,\rho_{1/2}^{n} \wedge \kappa) - \lambda \mathcal{R}_{s^{n+1}}(\kappa,q^{n+1}). $

    By substracting these last two inequalities, we get

    $ |ρn+11/2κ|=ρn+11/2κρn+11/2κHn+11(ρn3/2κ,ρn1/2κ,ρn1/2κ)Hn+11(ρn3/2κ,ρn1/2κ,ρn1/2κ)+λRsn+1(κ,qn+1)=|ρn1/2κ|λ(min{Fn+1(ρn1/2κ,ρn1/2κ),qn+1}Fn+1(ρn1/2κ,ρn1/2κ))+λ(min{Fn+1(ρn1/2κ,ρn1/2κ),qn+1}Fn+1(ρn1/2κ,ρn1/2κ))+λRsn+1(κ,qn+1)=|ρn1/2κ|λ(Φn0Φn1)+λ(Φn0Φnint)+λRsn+1(κ,qn+1), $

    which is exactly (29) in the case $ j = -1 $. The case $ j = 0 $ is similar so we omit the details of the proof for this case.

    Starting from (24) and (29), we can obtain approximate versions of (4) and (5). Let us introduce the functions:

    $ \Phi_{\Delta}\left( \rho_{\Delta},\kappa \right) = \sum\limits_{n = 0}^{N} \sum\limits_{j \in \mathbb{Z}} \Phi_{j}^{n} \mathbb{1}_{ \mathcal{P}_{j+1/2}^{n}}; \; \mathbf{F}_{\Delta}(s_{\Delta},\rho_{\Delta}) = \sum\limits_{n = 0}^{N} \sum\limits_{j \in \mathbb{Z}} \mathbf{F}^{n+1}(\rho_{j-1/2}^{n},\rho_{j+1/2}^{n}) \mathbb{1}_{ \mathcal{P}_{j+1/2}^{n}}. $

    Proposition 4 (Approximate entropy/constraint inequalities). (i) Fix $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times \mathbb{R}^{+}), \varphi \geq 0 $ and $ \kappa \in [0,R] $. Then there exists a constant $ \mathbf{C}^{{\varphi}}_1 = \mathbf{C}^{{\varphi}}_1(R,T,\mathbf{L}^{{}}) $, nondecreasing with respect to its arguments, such that the following inequalities hold for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττR|ρΔκ|tφ+ΦΔ(ρΔ,κ)xφdxdt+R|ρΔ(x,τ)κ|φ(x,τ)dxR|ρΔ(x,τ)κ|φ(x,τ)dx+2ττRsΔ(t)(κ,qΔ(t))φ(0,t)dtCφ1(Δt+Δx). $ (30)

    (ii) Fix $ \psi \in \mathbf{C}^{{\infty}}([0,T]), \psi \geq 0 $ and $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}) $ such that $ \varphi(0) = 1 $. Then there exists a constant $ \mathbf{C}^{{\varphi, \psi}}_2 = \mathbf{C}^{{\varphi, \psi}}_2 (R,T,\mathbf{L}^{{}},\|Q\|_{\mathbf{L}^{{\infty}}}) $, nondecreasing with respect to its arguments, such that for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττR+ρΔt(φψ)+FΔ(sΔ,ρΔ)x(φψ)dxdtR+ρΔ(x,τ)φ(x)ψ(τ)dx+R+ρΔ(x,τ)φ(x)ψ(τ)dxττqΔ(t)ψ(t)dt+Cφ,ψ2(Δx+Δt). $ (31)

    Proof. Fix $ k,m \in \mathbb{N} $ such that $ \tau \in [t^{k},t^{k+1}) $ and $ \tau' \in [t^{m},t^{m+1}) $.

    (i) Define for all $ n \in \mathbb{N} $ and $ j \in \mathbb{Z} $, $ {{\varphi_{j+1/2}^{n} = \frac{1}{\Delta x \Delta t} \iint_{ \mathcal{P}_{j+1/2}^{n}} \varphi(x,t) \mathinner{\mathrm{d}{x}} \mathinner{\mathrm{d}{t}}}} $. Multiplying the discrete entropy inequalities (29) by $ \varphi_{j+1/2}^{n} $, then summing over $ n \in \{k,\ldots,m-1\} $ and $ j \in \mathbb{Z} $, one obtains after reorganization of the sums (using in particular the Abel/"summation-by-parts" procedure)

    $ A+B+C+D+E0, $ (32)

    with

    $ A=m1n=k+1jZ|ρnj+1/2κ|(φnj+1/2φn1j+1/2)ΔxB=m1n=kjZΦnj(φnj+1/2φnj1/2)ΔtC=jZ|ρkj+1/2κ|φkj+1/2ΔxjZ|ρmj+1/2κ|φm1j+1/2ΔxD=m1n=kRsn+1(κ,qn+1)(φn1/2+φn1/2)ΔtE=m1n=k(Φn0Φnint)(φn1/2φn1/2)Δt. $

    Inequality (30) follows from (32) with

    $ Cφ1=R(T maxt[0,T]2ttφ(,t)L1+4 maxt[0,T] tφ(,t)L1)+RL(T maxt[0,T] 2xxφ(,t)L1+2 maxt[0,T] xφ(,t)L1+4φL+2TxφL), $

    making use of the bounds:

    $ |AττR|ρΔκ|tφdxdt|R(T maxt[0,T] 2ttφ(,t)L1+2 maxt[0,T] tφ(,t)L1)Δt,|BττRΦΔ(ρΔ,κ)xφdxdt|RL(T maxt[0,T] 2xxφ(,t)L1Δx,+2 maxt[0,T] xφ(,t)L1Δt)|CR|ρΔ(x,τ)κ|φ(x,τ)dx+R|ρΔ(x,τ)κ|φ(x,τ)dx|2R maxt[0,T] tφ(,t)L1Δt,|D2ττRsΔ(t)(κ,qΔ(t))φ(0,t)dt|RL(4φLΔt+TxφLΔx)and|E|2RTLxφLΔx. $

    (ii) In this case, the constant $ \mathbf{C}^{{\varphi, \psi}}_2 $ reads

    $ Cφ,ψ2=RφL1(TψL+4ψL)+QLψL(2+TφL)+RLψL(2φL1+TφL1+TφL1). $

    Following the proof of (30), define for all $ n \in \mathbb{N} $ and $ j \in \mathbb{Z} $,

    $ \psi^{n} = \frac{1}{\Delta t} \int_{t^{n}}^{t^{n+1}} \psi(t) \mathinner{\mathrm{d}{t}}, \quad \text{and} \quad \varphi_{j+1/2} = \frac{1}{\Delta x} \int_{x_{j}}^{x_{j+1}} \varphi(x) \mathinner{\mathrm{d}{x}}, $

    multiply the scheme (24) by $ \varphi_{j+1/2} \psi^{n} $, then take the sum over $ n \in \{k,\ldots,m-1\} $ and $ j \geq 0 $. Since the proof is very similar to the one of (i), we omit the details.

    The final step is to obtain compactness for the sequences $ (\rho_{\Delta})_{\Delta} $ and $ (y_{\Delta})_{\Delta} $ in order to pass to the limit in (30)-(31). We start with $ (y_{\Delta})_{\Delta} $.

    Proposition 5. For all $ t \in [0,T] $,

    $ yΔ(t)=y0+t0ω(RρΔ(x,u)μ(x)dx)du. $ (33)

    Consequently, there exists $ y \in \mathbf{C}^{{}}([0,T]) $ such that up to an extraction, $ (y_{\Delta})_{\Delta} $ converges uniformly to $ y $ on $ [0,T] $.

    Proof. For all $ t \in [0,T] $, if $ t \in [t^{n},t^{n+1}) $ for some $ n \in \{0, \ldots, N\} $, then we can write

    $ yΔ(t)y0=n1k=0tk+1tksk+1du+ttnsn+1du=n1k=0tk+1tkω(jZRρkj+1/2μj+1/2Δx)du+ttnω(jZRρnj+1/2μj+1/2Δx)du=t0ω(RρΔ(x,u)μ(x)dx)du. $

    Let us also point out that from (23), we get that for all $ \Delta $ and almost every $ t \in (0,T) $,

    $ sΔ(t)=ω(RρΔ(x,t)μ(x)dx). $ (34)

    Combining (33) and (34), we obtain that for all $ \Delta $,

    $ \|\dot y_{\Delta}\|_{\mathbf{L}^{{\infty}}} = \|s_{\Delta}\|_{\mathbf{L}^{{\infty}}} \leq \|\omega\|_{\mathbf{L}^{{\infty}}} \quad \text{and} \quad \|y_{\Delta}\|_{\mathbf{L}^{{\infty}}} \leq |y_{0}| + T \|\omega\|_{\mathbf{L}^{{\infty}}}. $

    The sequence $ (y_{\Delta})_{\Delta} $ is therefore bounded in $ {\mathbf{W}^{{1},{\infty}}}((0,T)) $. Making use of the compact embedding of $ {\mathbf{W}^{{1},{\infty}}}((0,T)) $ in $ \mathbf{C}^{{}}([0,T]) $, we get the existence of $ y \in \mathbf{C}^{{}}([0,T]) $ such that up to the extraction of subsequence, $ (y_{\Delta})_{\Delta} $ converges uniformly to $ y $ on $ [0,T] $.

    The presence of a time dependent flux in the conservation law of (1) complicates the obtaining of compactness for $ (\rho_{\Delta})_{\Delta} $. In particular, the techniques used in [10,11] to derive localized $ \mathbf{BV} $ estimates don't apply here since our problem lacks time translation invariance. In the present situation, it would be possible to derive weak $ \mathbf{BV} $ estimates ([6,23]). We choose different options. Similarly to what we did in Section 3, we propose two ways to obtain compactness, which will lead to two convergence results.

    First, we choose to adapt techniques and results put forward by Towers in [33]. With this in mind, we suppose in this section that $ f \in \mathbf{C}^{{2}}([0,R]) $ and strictly concave. Therefore,

    $ α>0, ρ[0,R],f(ρ)α. $ (35)

    Though this assumption is stronger than the nondegeneracy one (15), since $ f $ is bell-shaped, these two assumptions are similar in their spirit. We will also assume, following [33], that

    $ the numerical flux chosen in (24) is eitherthe Engquist-Osher one or the Godunov one. $ (36)

    To be precise, the choice made for the numerical flux at the interface – i.e. when $ j = 0 $ in (25) – does not play any role. What is important is that away from the interface, one chooses either the Engquist-Osher flux or the Godunov one. We denote for all $ n \in \{0, \ldots, N+1\} $ and $ j \in \mathbb{Z} $,

    $ \mathbf{D}_{j}^{n} = \max \left\{ \rho_{j-1/2}^{n} - \rho_{j+1/2}^{n},0 \right\}. $

    We will also use the notation

    $ \hat{ \mathbb{Z}} = \mathbb{Z} \backslash \{-1,0,1\}. $

    In [33], the author dealt with a discontinuous in both time and space flux and the specific "vanishing viscosity" coupling at the interface. The discontinuity in space was localized along the curve $ \{x = 0\} $. Here, we deal with only a discontinuous in time flux, but we also have a flux constraint along the curve $ \{x = 0\} $ since we work in the bus frame. The applicability of the technique of [33] for our case with moving interface and flux-constrained interface coupling relies on the fact that one can derive a bound on $ \mathbf{D}_{j}^{n} $ as long as the "interface" does not enter the calculations for $ \mathbf{D}_{j}^{n} $ i.e. $ j \in \hat{\mathbb{Z}} $. This is what the following lemma points out under Assumptions (35)-(36). For readers' convenience and in order to highlight the generality of the technique of Towers [33], let us provide the key elements of the argumentation leading to compactness.

    Lemma 4.1. Let $ n \in \{0, \ldots, N-1\} $ and $ j \in \hat{ \mathbb{Z}} $. Then if $ {{a = \frac{\lambda \alpha}{4}}} $, we have

    $ Dn+1jmax{Dnj1,Dnj,Dnj+1}a(max{Dnj1,Dnj,Dnj+1})2 $ (37)

    and

    $ Dn+1j1min{|j|1,n+1}a. $ (38)

    Proof. (Sketched) Inequality (38) is an immediate consequence of inequality (37), see [33,Lemma 4.3]. Obtaining inequality (37) however, is less immediate. Let us give some details of the proof.

    First, note that by introducing the function $ \psi : z \mapsto z - az^{2} $, inequality (37) can be stated as:

    $ Dn+1jψ(max{Dnj1,Dnj,Dnj+1}). $ (39)

    Then, one can show – only using the monotonicity of both the scheme and the function $ \psi $ – that under the assumption inequality (39) holds when

    $  (ρnj+3/2ρnj+1/2), (ρnj1/2ρnj3/2)0, $ (40)

    it follows that inequality (39) holds for all cases. And finally in [33,Page 23], the author proves that if the flux considered is either the Engquist-Osher flux or the Godunov flux, then (40) holds.

    The following lemma is an immediate consequence of inequality (38).

    Lemma 4.2. Fix $ 0 < \varepsilon < X $. Let $ i,J \in \mathbb{N}^{*} $ such that $ \varepsilon \in \mathbf{K}_{i+1/2} $ and $ X \in \mathbf{K}_{J-1/2} $. Then if $ \Delta x / \varepsilon $ is sufficiently small, there exists a constant $ {{ \mathbf{B} = \mathbf{B}\left(R,X,\frac{1}{a},\frac{1}{ \varepsilon}\right)}} $, nondecreasing with respect to its arguments, such that for all $ n \geq i-1 $,

    $ J1j=i|ρnj+1/2ρnj1/2|,i1j=J+1|ρnj+1/2ρnj1/2|B $ (41)

    and

    $ J2j=i|ρn+1j+1/2ρnj+1/2|,i2j=J+1|ρn+1j+1/2ρnj+1/2|2λLB. $ (42)

    Proposition 6. There exists $ \rho \in \mathbf{L}^{{\infty}}( \mathbb{R} \times (0,T)) $ such that up to the extraction of a subsequence, $ (\rho_{\Delta})_{\Delta} $ converges almost everywhere to $ \rho $ in $ \mathbb{R} \times (0,T) $.

    Proof. Fix $ 0 < \varepsilon < X $ and $ t > \lambda \varepsilon $. Denote by $ {{\Omega(X, \varepsilon) = (-X,- \varepsilon) \cup ( \varepsilon, X)}} $. Introduce $ i,J,n \in \mathbb{N} $ such that $ \varepsilon \in \mathbf{K}_{i+1/2} $, $ X \in \mathbf{K}_{J-1/2} $ and $ t \in [t^{n},t^{n+1}) $. Remark that

    $ (n+1) \Delta t > t > \lambda \varepsilon \geq \lambda (i \ \Delta x) = i \Delta t, $

    i.e. $ n \geq i-1 $. Then if we suppose that $ \Delta x / \varepsilon $ is sufficiently small, we can use Lemma 4.2. From (41), we get

    $ TV(ρΔ(t)|Ω(X,ε))2B $ (43)

    and from (42), we deduce

    $ Ω(X,ε)|ρΔ(x,t+Δt)ρΔ(x,t)|dx4LBΔt. $ (44)

    Combining (43)-(44) and the $ \mathbf{L}^{{\infty}} $ bound (28), a functional analysis result ([24,Theorem A.8]) ensures the existence of a subsequence which converges almost everywhere to some $ \rho $ on $ \Omega(X, \varepsilon) \times (\lambda \varepsilon,T) $. By a standard diagonal process we can extract a further subsequence (which we do not relabel) such that $ (\rho_{\Delta})_{\Delta} $ converges almost everywhere to $ \rho $ on $ \mathbb{R} \times (0,T) $.

    Theorem 4.3. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) $ and $ y_{0} \in \mathbb{R} $. Suppose that $ f \in \mathbf{C}^{{2}} $ satisfies Assumptions (2)-(7)-(35). Suppose also that in (25), we use the Engquist-Osher flux or the Godunov one when $ j \neq 0 $ and any other monotone consistent and Lipschitz numerical flux when $ j = 0 $. Then under the CFL condition (27), the scheme (23) – (25) converges to an admissible weak solution to Problem (1).

    Proof. We have shown that – up to the extraction of a subsequence – $ y_{\Delta} $ converges uniformly on $ [0,T] $ to some $ y \in \mathbf{C}^{{}}([0,T]) $ and that $ \rho_{\Delta} $ converges a.e. on $ \mathbb{R} \times (0,T) $ to some $ \rho \in \mathbf{L}^{{\infty}}( \mathbb{R} \times (0,T)) $. We now prove that this couple $ (\rho,y) $ is an admissible weak solution to Problem (1) in the sense of Definition 2.1.

    Recall that for all $ \Delta $ and $ t \in [0,T] $,

    $ y_{\Delta}(t) = y_{0} + \int_{0}^{t} \omega \left( \int_{ \mathbb{R}} \rho_{\Delta}(x,u) \mu(x) \mathinner{\mathrm{d}{x}} \right) \mathinner{\mathrm{d}{u}}. $

    When letting $ \Delta \to 0 $, the dominated convergence theorem ensures that $ y $ satisfies (6). Apply inequality (30) with $ \tau = 0 $, $ \tau' = T $, $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}^{*} \times [0,T)), \varphi \geq 0 $ and $ \kappa \in [0,R] $ to obtain

    $ \int_{0}^{T} \int_{ \mathbb{R}} |\rho_{\Delta}-\kappa| \partial_{t}\varphi + \Phi_{\Delta}(\rho_{\Delta},\kappa) \partial_{x} \varphi \mathinner{\mathrm{d}{x}} \mathinner{\mathrm{d}{t}} + \int_{ \mathbb{R}}|\rho_{\Delta}^0-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq - \mathbf{C}^{{\varphi}}_1 (\Delta x + \Delta t). $

    Then the a.e. convergence of $ (s_{\Delta})_{\Delta} $ to $ \dot y $ – coming from (34) – and the a.e. convergence of $ (\rho_{\Delta})_{\Delta} $ to $ \rho $ ensure that when letting $ \Delta \to 0 $, we get

    $ \int_{0}^{T} \int_{ \mathbb{R}} |\rho-\kappa| \partial_{t}\varphi + \Phi_{\dot y(t)}(\rho,\kappa) \partial_{x} \varphi \mathinner{\mathrm{d}{x}} \mathinner{\mathrm{d}{t}} + \int_{ \mathbb{R}}|\rho_{0}(x+y_{0})-\kappa|\varphi(x,0) \mathinner{\mathrm{d}{x}} \geq 0, $

    and consequently $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $, see Remark 2. Now, we pass to the limit in (30) and (31) using the a.e. convergence of $ (s_{\Delta})_{\Delta} $ to $ \dot y $ and of $ (\rho_{\Delta})_{\Delta} $ to $ \rho $ as well as the continuity of $ Q $ and $ \omega $. Consequently, for all test functions $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times \mathbb{R}^{+}), \varphi \geq 0 $ and $ \kappa \in [0,R] $, the following inequalities hold for almost every $ 0 \leq \tau < \tau' \leq T $:

    $ ττR|ρκ|tφ+Φ˙y(t)(ρ,κ)xφdxdt+R|ρ(x,τ)κ|φ(x,τ)dxR|ρ(x,τ)κ|φ(x,τ)dx+2ττR˙y(t)(κ,q(t))φ(0,t)dt0. $

    To conclude, note that the expression in the left-hand side of the previous inequality is a continuous function of $ (\tau,\tau') $ which is almost everywhere greater than the continuous function $ 0 $. By continuity, this expression is everywhere greater than $ 0 $, which proves that $ \rho $ satisfies the entropy inequalities (4). Using similar arguments, one shows that $ \rho $ also satisfies the constraint inequalities (5). This shows that the couple $ (\rho,y) $ is an admissible weak solution to (1), and that concludes the proof of convergence.

    We proved than in the $ \mathbf{L}^{{\infty}} $ framework, the scheme converges to an admissible weak solution, but note that there is no guarantee of uniqueness in this construction. Also stress that we cannot extend this result to general consistent monotone numerical fluxes beyond hypothesis (36).

    The following result is the discrete version of Lemma 2.3 so it is consistent that the proof uses the discrete analogous arguments of the ones we used in the proof of Lemma 2.3.

    Lemma 4.4. Introduce for all $ \Delta > 0 $ the function $ \xi_{\Delta} $ defined for all $ t \in [0,T] $ by

    $ \xi_{\Delta}(t) = \int_{ \mathbb{R}} \rho_{\Delta}(x,t) \mu(x) \mathinner{\mathrm{d}{{x}}}. $

    Then $ \xi_{\Delta} $ has bounded variation and consequently, so does $ s_{\Delta} $.

    Proof. Since $ \mu \in \mathbf{BV}( \mathbb{R}) $, there exists a sequence of smooth functions $ (\mu_{\ell})_{\ell \in \mathbb{N}} \subset \mathbf{BV}( \mathbb{R}) \cap \mathbf{C}_c^{{\infty}}( \mathbb{R}) $ such that

    $ \| \mu_{\ell} - \mu \|_{\mathbf{L}^{{1}}} {\ \underset{{\ell} \to {+\infty}}{\longrightarrow} \ } 0 \quad \text{and} \quad \mathbf{TV}(\mu_{\ell}) {\ \underset{{\ell} \to {+\infty}}{\longrightarrow} \ } \mathbf{TV}(\mu). $

    Introduce for all $ \ell \in \mathbb{N} $ and $ t \in [0,T] $, the function $ {{\xi_{\Delta, \ell}(t) = \int_{ \mathbb{R}} \rho_{\Delta}(x,t) \mu_{\ell}(x) \mathinner{\mathrm{d}{x}}}} $ and let $ K > 0 $ such that

    $ \forall \ell \in \mathbb{N}, \quad \|\mu_{\ell}\|_{\mathbf{L}^{{1}}}, \mathbf{TV}(\mu_{\ell}) \leq K. $

    For all $ \ell \in \mathbb{N} $ and $ t,s \in [0,T] $, if $ t \in [t^{k},t^{k+1}) $ and $ s \in [t^{m},t^{m+1}) $, we have

    $ |ξΔ,(t)ξΔ,(s)|=|ξΔ,(tk)ξΔ,(tm)|=|RρΔ(x,tk)μ(x)dxRρΔ(x,tm)μ(x)dx|=|jZ(ρkj+1/2ρmj+1/2)μj+1/2Δx|,μj+1/2=1Δxxj+1xjμ(x)dx=|jZk1τ=m(ρτ+1j+1/2ρτj+1/2)μj+1/2Δx|=|k1τ=mjZ(Fτ+1j(ρτj1/2,ρτj+1/2)Fτ+1j+1(ρτj+1/2,ρτj+3/2))μj+1/2Δt|=|k1τ=mjZFτ+1j+1(ρτj+1/2,ρτj+3/2)(μj+3/2μj+1/2)Δt|RLk1τ=mTV(μ)ΔtRLK(|ts|+2Δt). $

    Consequently, for all $ \ell \in \mathbb{N} $, $ \Delta > 0 $ and $ t,\tau \in [0,T] $, the triangle inequality yields:

    $ \left|\xi_{\Delta}(t) - \xi_{\Delta}(\tau) \right| \leq 2R \|\mu-\mu_{\ell}\|_{\mathbf{L}^{{1}}} + R \mathbf{L}^{{}}K (|t-\tau| + 2\Delta t). $

    Letting $ \ell \to +\infty $, we get that for all $ \Delta > 0 $ and $ t,\tau \in [0,T] $,

    $ \left|\xi_{\Delta}(t) - \xi_{\Delta}(\tau) \right| \leq R \mathbf{L}^{{}}K (|t-\tau| + 2\Delta t), $

    which leads to

    $ \mathbf{TV}(\xi_{\Delta}) = \sum\limits_{k = 0}^{N} \left|\xi_{\Delta}(t^{k+1}) - \xi_{\Delta}(t^{k}) \right| \leq 3 R \mathbf{L}^{{}}K (T + \Delta t). $

    This proves that $ \xi_{\Delta} \in \mathbf{BV}([0,T]) $. Since $ \omega $ is Lipschitz continuous, $ s_{\Delta} $ also has bounded variation.

    Theorem 4.5. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $ and $ y_{0} \in \mathbb{R} $. Suppose that $ f $ satisfies (2)-(7)-(16) and that $ Q $ satisfies (17). Suppose also that in (25), we use the Godunov flux when $ j = 0 $ and any other monotone consistent and Lipschitz numerical flux when $ j \neq 0 $. Then under the CFL condition (27), the scheme (23) – (25) converges to a $ \mathbf{BV} $-regular solution to Problem (1).

    Proof. All the hypotheses of Lemma A.4 are fulfilled. Consequently, there exists a constant $ \mathbf{C}^{{ \varepsilon}} > 0 $ such that for all $ n \in \{0, \ldots, N-1\} $,

    $ TV(ρΔ(tn+1))TV(ρ0)+4R+Cε(nk=0|qk+1qk|+nk=0|sk+1sk|)TV(ρ0)+4R+Cε(1+QL)nk=0|sk+1sk|. $ (45)

    Making use of Lemma 4.4, we obtain that for all $ n \in \{0, \ldots, N\} $,

    $ nk=0|sk+1sk|=nk=0|sΔ(tk+1)sΔ(tk)|ωLnk=0|ξΔ(tk+1)ξΔ(tk)|3RLKωL(T+Δt). $

    where the constant $ K $ was introduced in the proof of Lemma 4.4. The two last inequalities imply that for all $ t \in [0,T] $, we have

    $ TV(ρΔ(t))TV(ρ0)+4R+3Cε(1+QL)ωLRLK(T+Δt). $ (46)

    Therefore, the sequence $ (\rho_{\Delta})_{\Delta} $ is uniformly in time bounded in $ \mathbf{BV}( \mathbb{R}) $. Using [22,Appendix], we get the existence of $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $ such that

    $ \forall t \in [0,T], \ \rho_{\Delta}(t) {\ \underset{{\Delta} \to {0}}{\longrightarrow} \ } \rho(t) \ \text{in} \ \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}). $

    Following the proof of Theorem 4.3, we show that $ (\rho,y) $ is an admissible weak solution. Then passing to the limit in (46), the lower semi-continuity of the $ \mathbf{BV} $ semi-norm ensures that $ (\rho,y) $ is also $ \mathbf{BV} $-regular.

    Remark 10. Note the complementarity of the hypotheses made in the above theorem with the ones of Theorem 4.3. Recall that in Theorem 4.3, we needed the Godunov flux only away from the interface.

    In this section we present some numerical tests performed with the scheme analyzed in Section 4. In all the simulations we take the uniformly concave flux $ f(\rho) = \rho(1-\rho) $ (the maximal car velocity and the maximal density are assumed to be equal to one). Following the hypotheses of Theorem 4.5, we choose the Godunov flux at the interface, and the Rusanov one away from the interface. We will use weight functions of the kind

    $ \mu_k (x) = 2^k \mathbb{1}_{\left[ 0 ; \frac{1}{2^k} \right]}(x), $

    for one (in Section 5.1) or several (in Section 5.3) values of $ k \in \mathbb{N}^* $.

    In this section, consider a two-lane road on which a bus travels with a speed given by the function

    $ \omega (\rho) = \left\{ α(β+ρ)2if0ρ0.61ρif0.6ρ1, \right. $

    where $ \alpha $ and $ \beta $ are chosen so that $ \omega(0) = 0.7 $ and $ \omega(0.6) = 0.4 $, as illustrated in Figure 1 (left). The set-up of the experiment is the following. Consider a domain of computation $ [0, 11] $, the weight function $ \mu_4 $ and the following data:

    $ \rho_{0}(x) = 0.5 \mathbb{1}_{[0.5 ; 1]} (x), \quad y_{0} = 1.5, \quad Q(s) = 0.75 \times \left( \frac{1-s}{2} \right)^2. $
    Figure 1.  Evolution in time of the subjective density $ \xi $ and the bus velocity $ \dot y $ (1280 cells).

    The idea behind the choice of $ Q $ is that in average (between the two lanes), the presence of the slow vehicle reduces by $ 25\% $ the maximum traffic flow. As we can see in Figure 1 (right), the slow vehicle nearly always travels at maximum velocity. It makes sense because even though we can see that cars are overtaking it (Figure 1, right and Figure 2), the density $ \xi $ ahead of it is never sufficiently important to make it go slower.

    Figure 2.  The numerical solution at different fixed times.

    Remark 11. The function $ \omega $ we chose above is not of the form as required in [27,28]. Once again, let us stress that the particular form $ \omega(\rho) = \min\left\{V_{\text{bus}}; 1-\rho \right\} $, where $ V_{\text{bus}} $ is the maximum bus velocity, is crucial for the well-posedness result of [27,28] to hold. Indeed, it is essential in the analysis of [27,28] that the velocity of the bus be constant (equal to $ V_{\text{bus}}) $ across the nonclassical shocks. Our nonlocal model is not bound to this restriction.

    We also perform a convergence analysis for this test. In Table 1, we computed the relative errors

    $ E_\rho^\Delta = \| \rho_{\Delta} - \rho_{\Delta/2} \|_{\mathbf{L}^{{1}}((0,T) ;\mathbf{L}^{{1}}( \mathbb{R}))} \quad \text{and} \quad E_y^\Delta = \|y_{\Delta} - y_{\Delta/2} \|_{\mathbf{L}^{{\infty}}}, $
    Table 1.  Measured errors at time $ T = 13 $.
    Number of cells $ E_{\rho}^\Delta \ (\times 10^{-2}) $ $ E_{y}^\Delta \ (\times 10^{-3}) $
    160 $ 24.053 $ $ 48.0643 $
    320 $ 15.731 $ $ 15.939 $
    640 $ 9.647 $ $ 7.698 $
    1280 $ 6.197 $ $ 3.715 $
    2560 $ 3.226 $ $ 1.777 $
    5120 $ 1.936 $ $ 0.889 $
    10240 $ 1.055 $ $ 0.443 $

     | Show Table
    DownLoad: CSV
    Figure 3.  Rates of convergence at time $ T = 13 $.

    for different number of space cells at the final time $ T = 13 $. We see that those ratio converge with convergence orders approximately equal to $ 0.76 $ for the car density and approximately equal to $ 1.1 $ for the slow moving vehicle position.

    Now we confront the numerical tests performed with our model with the tests done by the authors in [14] approximating the original problem of [20]. We deal with a road of length $ 1 $ parametrized by the interval $ [0,1] $ and choose the weight function $ \mu_3 $. Moreover,

    $ \omega(\rho) = \min\{0.3 ; \ 1 - \rho \} \quad \text{and} \quad Q(s) = 0.6 \times \left( \frac{1-s}{2} \right)^2. $

    First, consider the initial datum

    $ ρ0(x)={0.4ifx<0.50.5ifx>0.5 y0=0.5. $ (47)

    The numerical solution is composed of two classical shocks separated by a nonclassical discontinuity, as illustrated in Figure 4 (left).

    Figure 4.  Evolution in time of the numerical density corresponding to initial data (47) (left) and (48) (right), with 1280 cells.

    Next, we choose

    $ ρ0(x)={0.8ifx<0.50.5ifx>0.5 y0=0.5. $ (48)

    The values of the initial condition create a rarefaction wave followed by a nonclassical and classical shocks, as illustrated in Figure 4 (right).

    Finally, still following [14], we consider

    $ ρ0(x)={0.8ifx<0.50.4ifx>0.5 y0=0.4. $ (49)

    Here the solution is composed of a rarefaction wave followed by nonclassical and classical shocks on the density that are created when the slow vehicle approaches the rarefaction and initiates a moving bottleneck, as illustrated in Figure 5.

    Figure 5.  Evolution in time of the numerical density corresponding to initial data (49) (1280 cells).

    With these three tests, we can already see – in a qualitative way – the resemblance between the numerical approximations to the solutions to our model and the numerical approximations of [14]. One way to quantify their proximity is for example to evaluate the $ \mathbf{L}^{{1}} $ error between the car densities and the $ \mathbf{L}^{{\infty}} $ error between the bus positions. More precisely, denote by $ (\rho_{\Delta},y_{\Delta}) $ the approximation of the $ \mathbf{BV} $-regular solution to (1) obtained with the scheme (23) – (25), and denote by $ (\overline{\rho}_{\Delta},\overline{y}_{\Delta}) $ the couple obtained with this same scheme but

    $ \text{replacing} \quad s^{n+1} = \omega \left( \sum\limits_{j \in \mathbb{Z}} \rho_{j+1/2}^{n} \mu_{j+1/2} \Delta x \right) \quad \text{by} \quad s^{n+1} = \omega \left(\rho_{1/2}^{n} \right). $

    Let us precise that this is not the scheme the authors of [14] proposed. However, this scheme is consistent with the problem

    $ {tρ+x(F(˙y(t),ρ))=0R×(0,T)ρ(x,0)=ρ0(x+y0)xRF(˙y(t),ρ)|x=0Q(˙y(t))t(0,T)˙y(t)=ω(ρ(0+,t))t(0,T)y(0)=y0 $ (50)

    and behaves in a stable way in the calculations we performed. Therefore, the couple $ (\overline{\rho}_{\Delta},\overline{y}_{\Delta}) $ is expected to give a reasonable approximation of the solution to (50). With this in mind, for the case (49) and still with the weight function $ \mu_3 $, we computed in Table 2 the measured errors

    $ E_{\mathbf{L}^{{1}}}^\Delta = \|\rho_{\Delta} - \overline{\rho}_{\Delta} \|_{\mathbf{L}^{{1}}((0,T) ; \mathbf{L}^{{1}}( \mathbb{R}))} \quad \text{and} \quad E_{\mathbf{L}^{{\infty}}}^\Delta = \|y_{\Delta} - \overline{y}_{\Delta} \|_{\mathbf{L}^{{\infty}}}. $
    Table 2.  Measured errors at time $ T = 0.7245 $.
    Number of cells $ E_{\mathbf{L}^{{1}}}^\Delta \ (\times 10^{-4}) $ $ E_{\mathbf{L}^{{\infty}}}^\Delta \ (\times 10^{-3}) $
    160 $ 32.672 $ $ 18.519 $
    320 $ 14.236 $ $ 7.341 $
    640 $ 5.837 $ $ 3.701 $
    1280 $ 3.833 $ $ 4.879 $
    2560 $ 3.207 $ $ 6.405 $
    5120 $ 2.922 $ $ 7.144 $
    10240 $ 2.776 $ $ 7.501 $
    20480 $ 2.698 $ $ 7.674 $
    40960 $ 2.658 $ $ 7.759 $

     | Show Table
    DownLoad: CSV

    These calculations indicate that for a sufficiently large number of cells $ J \geq 40960 $,

    $ E_{\mathbf{L}^{{1}}}^\Delta \simeq 2.7 \times 10^{-4} \quad \text{and} \quad E_{\mathbf{L}^{{\infty}}}^\Delta \simeq 7.6 \times 10^{-3}. $

    This indicates the discrepancy between our nonlocal and the local model (50) of [20]. The idea is now to fix the number of cells $ J = 40960 $ and to make the length of the weight function support go to zero. In Table 3, we have computed, for different weight functions, the error between the approximations of the two models. This error corresponds, as in the above calculation, to the residual error observed starting from a sufficiently small $ \Delta x $.

    Table 3.  Measured errors at time $ T = 0.7245 $.
    weight function $ E_{\mathbf{L}^{{1}}}^\Delta $ $ E_{\mathbf{L}^{{\infty}}}^\Delta $
    $ \mu_1 $ $ 6.810 \times 10^{-3} $ $ 5.489 \times 10^{-2} $
    $ \mu_2 $ $ 1.105 \times 10^{-3} $ $ 1.972 \times 10^{-2} $
    $ \mu_3 $ $ 2.658 \times 10^{-4} $ $ 7.759 \times 10^{-3} $
    $ \mu_4 $ $ 9.232 \times 10^{-5} $ $ 2.913 \times 10^{-3} $
    $ \mu_5 $ $ 6.190 \times 10^{-5} $ $ 9.110 \times 10^{-4} $

     | Show Table
    DownLoad: CSV

    Remark 12. The previous simulations show a closeness between our model as $ \mu \to \delta_{0^{+}} $ and (50). Let us however point that the nonlocality in space for the slow vehicle introduces an undesirable artefact into the model. In the rarefaction regime one may observe that the large vehicle may move a bit faster that the surrounding flow. The situation where this effect becomes truly perceptible is when considering initial data of the type

    $ ρ0(x)={1ifx<xb0ifx>xb. $ (51)

    Indeed, for such data, there exists a small time interval $ [0, \delta] $ in which $ \dot y(t) > v(\rho(y(t)^+)) = 0 $, which would suggest that the slow vehicle moves forward while the cars in front of it do not. This time interval is in fact quite small due to the narrowness of the support of the weight function. The local model does not develop such phenomena. This qualitative artefact precludes us from giving a microscopic interpretation to the model, which main output is the global influence of the slow vehicle on the flow; however, let us stress that the phenomenon becomes quantitatively negligible for larger times. Indeed, Ole\u{ı}nik estimate on decay of positive waves ensures that data of the type (51) evolve into rarefaction waves and do not appear while driving: the classical LWR model precludes the formation of rarefaction waves focused at positive time. The modification of the classical LWR brought by the constraint may produce nonclassical waves at positive times; while these waves are downward jumps in density like in (51), they are situated precisely at the location of the constraint and not slightly behind it, like in (51).

    Even if we are unable, at this time, to rigorously link our problem (1) with $ \mu \to \delta_{0^{+}} $ and the original problem (50) of the authors in [20], this last experiment corroborates the conjecture that the local model (50) is the singular limit of our model in the case $ \omega $ is of the form $ \omega(\rho) = \min\left\{V_{\text{bus}} ; 1-\rho \right\} $. The other interesting question is whether the local model is well posed beyond this particular choice of $ \omega $.

    The author is most grateful to Boris Andreianov for his constant support and many enlightening discussions.

    We focus on the study of the following class of models:

    $ {tρ+x(F(s(t),ρ))=0R×(0,T)ρ(x,0)=ρ0(x)xRF(s(t),ρ)|x=0q(t)t(0,T), $ (52)

    where $ s \in \mathbf{BV}([0,T]; [0,\Sigma]) $ for some $ \Sigma > 0 $ and $ q \in \mathbf{BV}([0,T] ; \mathbb{R}^{+}) $. We suppose that $ F \in \mathbf{C}^{{1}}([0, \Sigma] \times [0,R]) $ and that for all $ s \in [0,\Sigma] $, $ F(s, \cdot) $ is bell-shaped i.e.

    $ s[0,Σ], F(s,0)=0, F(s,R)0and! ¯ρs(0,R), ρF(s,ρ)(¯ρsρ)>0 for a.e. ρ(0,R). $ (53)

    This framework covers the particular case when $ F $ takes the form:

    $ F(s(t),\rho) = f(\rho) - s(t) \rho, $

    with bell-shaped $ f : [0,R] \to \mathbb{R}^+ $, which our model (1) is based on. This class of models is well known, especially when the flux function is not time dependent, cf. [17,6]. In this appendix, we establish in passing the well-posedness of Problem (52), but our main interest lies in the $ \mathbf{BV} $ in space regularity of the solutions. More precisely, we aim at obtaining a bound on the total variation of the solutions to (52), using a finite volume approximation which allows for sharp control of the variation at the constraint. Note that the alternative offered by wave-front tracking would be cumbersome because of the explicit time-dependency in (52). In the general case, entropy solutions to limited flux problems like (52) do not belong to $ \mathbf{L}^{{\infty}}((0,T) ; \mathbf{BV}( \mathbb{R})) $, see [1]. We will show that it is the case under a mild assumption on the constraint function $ q $ – see Assumption (59) below – and provided that

    $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}). $

    Throughout the appendix, for all $ s \in [0,\Sigma] $ and $ a,b \in [0,R] $, we denote by

    $ \Phi_{s}(a,b) = \text{sign}(a-b)(F(s,a) - F(s,b)) $

    the classical Kružkov entropy flux associated with the Kružkov entropy $ \rho \mapsto |\rho-\kappa| $, for all $ \kappa \in [0,R] $, see [25].

    Let us first recall the following definition.

    Definition A.1. A bounded function $ \rho \in \mathbf{L}^{{\infty}}( \mathbb{R} \times (0,T)) $ is an admissible weak solution to (52) if

    (ⅰ) the following regularity is fulfilled: $ {{\rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}))}} $;

    (ⅱ) for all test functions $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R} \times \mathbb{R}^{+}), \varphi \geq 0 $ and $ \kappa \in [0,R] $, the following entropy inequalities are verified for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττR|ρκ|tφ+Φs(t)(ρ,κ)xφdxdt+R|ρ(x,τ)κ|φ(x,τ)dxR|ρ(x,τ)κ|φ(x,τ)dx+2ττRs(t)(κ,q(t))φ(0,t)dt0, $

    where

    $ \mathcal{R}_{s(t)}(\kappa,q(t)) = F(s(t),\kappa) - \min\left\{ F(s(t),\kappa),q(t) \right\}; $

    (ⅲ) for all test functions $ \psi \in \mathbf{C}^{{\infty}}([0,T]), \psi \geq 0 $ and some given $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}) $ which verifies $ \varphi(0) = 1 $, the following weak constraint inequalities are verified for all $ 0 \leq \tau < \tau' \leq T $:

    $ ττR+ρt(φψ)+F(s(t),ρ)x(φψ) dxdtR+ρ(x,τ)φ(x)ψ(τ)dx+R+ρ(x,τ)φ(x)ψ(τ)dxττq(t)ψ(t)dt. $

    Definition A.2. An admissible weak solution $ \rho $ will be called $ \mathbf{BV} $-regular if it verifies $ {{\rho \in \mathbf{L}^{{\infty}}((0,T) ; \mathbf{BV}( \mathbb{R}))}} $.

    As we pointed out before, this notion of solution is well suited for passage to the limit of a.e. convergent sequences of exact or approximate solutions. However, it is not so well-adapted to prove uniqueness. An equivalent notion of solution, based on explicit treatment of traces of $ \rho $ at the constraint, was introduced by the authors of [7]. This notion of solution leads to the following stability estimate.

    Theorem A.3. Fix $ s^1, s^2 \in \mathbf{BV}([0,T] ; [0,\Sigma]) $, $ \rho_{0}^1,\rho_{0}^2 \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $ and $ q^1,q^2 \in \mathbf{BV}([0,T] ; \mathbb{R}^+) $. Denote by $ \rho^1 $ a $ \mathbf{BV} $-regular solution to (52) with data $ \rho_0^1, q^1, s^1 $ and $ \rho^2 $ an admissible weak solution to (52) with data $ \rho_0^2, q^2, s^2 $. Suppose that the flux functions $ (t,\rho) \mapsto F(s^1(t),\rho),F(s^2(t),\rho) $ satisfy (53). Then for a.e. $ t \in (0,T) $, we have:

    $ ρ1(t)ρ2(t)L1ρ10ρ20L1+2t0|q1(τ)q2(τ)|dτ+2t0F(s1(τ),)F(s2(τ),)Ldτ+t0||ρF(s1(τ),)ρF(s2(τ),)||LTV(ρ1(τ))dτ. $ (54)

    In particular, Problem (52) admits at most one $ \mathbf{BV} $-regular solution.

    Proof. Since our interest to details lies rather on the numerical approximation point of view, we do not fully prove this statement but we give the essential steps leading to this stability result.

    Definition of solution. First, the authors of [7] introduce a subset of $ \mathbb{R}^2 $ called germ, which can be seen as the set of all the possible traces of a solution to (52). Then, they say that $ \rho $ is a solution to (52) if it satisfies entropy inequalities away from the interface – i.e. with $ \varphi \in \mathbf{C}_c^{{\infty}}( \mathbb{R}^* \times \mathbb{R}^+) $ in the entropy inequalities – and if the couple constituted of left-side and the right-side traces of $ \rho $ belongs to this so-called germ.

    Equivalence of the two definitions. The next step is to prove that this latter definition of solution is equivalent to Definition A.1. This part is done using good choices of test functions, see [7,Theorem 3.18] or [6,Proposition 2.5,Theorem 2.9].

    First stability estimate. One first shows that if $ s^1 = s^2 $, then for a.e. $ t \in (0,T) $, one has

    $ ρ1(t)ρ2(t)L1ρ10ρ20L1+2t0|q1(τ)q2(τ)|dτ. $ (55)

    The proof starts with the classical doubling of variables method of Kružkov [25,Theorem 1] and then uses the germ structure, what the authors of [7] called $ \mathbf{L}^{{1}} $-dissipativity, see [7,Definition 3.1] and [6,Lemma 2.7].

    Proof of estimate (54). The proof is based upon estimate (55) and elements borrowed from [8,18]. Most details can be found in the proof of [21,Theorem 2.1].

    Remark 13. Though the definition of solutions with the germ explicitly involves the traces of $ \rho $, we did not discuss the existence of such traces. A first way to ensure such existence is to deal with $ \mathbf{BV} $-regular solutions. That way, traces do exist and are to be understood in the sense of $ \mathbf{BV} $ functions. Outside the $ \mathbf{BV} $ framework, existence of strong traces for solutions to (52) is ensured provided an assumption on the fundamental diagram like (15), see [2,29]. Finally, if one does not want to impose such a condition on the flux, (which is our case in this appendix), one can follow what the authors of [7] proposed (in Section 2) and consider the "singular mapping traces."

    We now turn to the proof of the existence of $ \mathbf{BV} $-regular solutions by the means of a finite volume scheme.

    Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) $. For a fixed spatial mesh size $ \Delta x $ and time mesh size $ \Delta t $, let $ x_{j} = j\Delta x $, $ t^{n} = n \Delta t $. We define the grid cells $ \mathbf{K}_{j+1/2} = (x_{j},x_{j+1}) $ and $ N \in \mathbb{N}^{*} $ such that $ T \in [t^N, t^{N+1}) $. We write

    $ \mathbb{R} \times [0,T] \subset \bigcup\limits_{n = 0}^{N} \bigcup\limits_{j \in \mathbb{Z}} \mathcal{P}_{j+1/2}^{n}, \quad \mathcal{P}_{j+1/2}^{n} = \mathbf{K}_{j+1/2} \times [t^{n},t^{n+1}). $

    We choose to discretize the initial data $ \rho_{0} $ and the functions $ s,q $ with $ \left( \rho_{j+1/2}^{0}\right)_{j} $, $ (s^{n})_{n} $ and $ (q^{n})_{n} $ where for all $ j \in \mathbb{Z} $ and $ n \in \{0 ,\ldots, N\} $), $ \rho_{j+1/2}^{0} $, $ s^{n} $ and $ q^n $ are their mean values on each cell $ \mathbf{K}_{j+1/2} $ and $ [t^{n},t^{n+1}) $. Following [6], the marching formula of the scheme is the following: for all $ n \in \{0 ,\ldots, N-1\} $ and $ j \in \mathbb{Z} $:

    $ ρn+1j+1/2=ρnj+1/2ΔtΔx(Fnj+1(ρnj+1/2,ρnj+3/2)Fnj(ρnj1/2,ρnj+1/2)), $ (56)

    where

    $ Fnj(a,b)={Fn(a,b)if j0min{Fn(a,b),qn)}if j=0, $ (57)

    $ \mathbf{F}^{n} $ being a monotone consistent and Lipschitz numerical flux associated to $ F(s^{n}, \cdot) $. We then define

    $ \rho_{\Delta}(x,t) = \rho_{j+1/2}^{n} \ \text{if} \ (x,t) \in \mathcal{P}_{j+1/2}^{n} \quad \text{and} \quad s_{\Delta}(t), q_{\Delta}(t) = s^{n}, q^n \ \text{if} \ t \in [t^{n},t^{n+1}). $

    Let $ \Delta = (\Delta x, \Delta t) $. For the convergence analysis, we will assume that $ \Delta \to 0 $, with $ \lambda = \Delta t / \Delta x $, verifying the CFL condition

    $ λsups[0,Σ](FsaL+FsbL)L1, $ (58)

    where $ \mathbf{F}^s = \mathbf{F}^s(a, b) $ is the numerical flux – associated to $ F(s,\cdot) $ – we use in the scheme (56). From now, the analysis of the scheme follows the same path as in Section 4. In that order, we prove that the scheme (56)-(57) is $ \mathbf{L}^{{\infty}} $ stable, satisfies discrete entropy inequalities similar to (29) and approximate entropy/constraint inequalities similar to (30)-(31). Only the compactness for $ (\rho_{\Delta})_{\Delta} $ is left to obtain since the $ \mathbf{L}_{\mathbf{loc}}^{{1}} $ compactness for the sequences $ (s_{\Delta})_{\Delta} $ and $ (q_{\Delta})_{\Delta} $ is clear. One way to do so is to derive uniform $ \mathbf{BV} $ bounds.

    Lemma A.4. We suppose that $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $ and that $ q $ verifies the assumption

    $ ε>0, t[0,T], s[0,Σ], q(t)maxρ[0,R]F(s,ρ)ε:=qε(s). $ (59)

    Then there exists a constant $ \mathbf{C}^{{ \varepsilon}} = \mathbf{C}^{{ \varepsilon}}(\| \partial_s F \|_{\mathbf{L}^{{\infty}}}) $ nondecreasing with respect to its argument such that for all $ n \in \{0, \ldots, N-1\} $,

    $ TV(ρΔ(tn+1))TV(ρ0)+4R+Cε(nk=0|qk+1qk|+nk=0|sk+1sk|), $ (60)

    where $ \rho_{\Delta} = \left( \rho^{n}_{j+1/2} \right)_{n,j} $ is the finite volume approximation constructed with the scheme (56)-(57), using the Godunov numerical flux when $ j = 0 $ in (57).

    Proof. Fix $ n \in \{0, \ldots N-1\} $. With this set up we can follow the proofs of [13,Section 2] to obtain the following estimate:

    $ jZ|ρn+1j+1/2ρn+1j1/2|TV(ρ0)+4R+2nk=0|(ˆρsk+1(qk+1)ˆρsk(qk))(˘ρsk+1(qk+1)˘ρsk(qk))|, $

    where for all $ k \in \{0, \ldots, n \} $, the couple $ \left( \widehat{\rho}_{s^k}(q^{k}), \breve{\rho}_{s^k}(q^{k}) \right) \in [0,R]^{2} $ is uniquely defined by the conditions

    $ F(s^{k},\widehat{\rho}_{s^k}(q^{k})) = F(s^{k},\breve{\rho}_{s^k}(q^{k})) = q^{k} \quad \text{and} \quad \widehat{\rho}_{s^{k}}(q^{k}) > \breve{\rho}_{s^k}(q^{k}). $

    Denote by $ \Omega( \varepsilon) $ the open subset

    $ \Omega( \varepsilon) = \bigcup\limits_{s \in [0,\Sigma]} \Omega_{s}( \varepsilon) $

    where for all $ s \in [0,\Sigma] $, $ {{\Omega_{s}( \varepsilon) = \left( \breve{\rho}_{s}(q_ \varepsilon (s)), \widehat{\rho}_{s}(q_ \varepsilon (s)) \right)}} $. By Assumption (59), the continuous function $ (s,\rho) \mapsto | \partial_{\rho}F(s,\rho)| $ is positive on the compact subset $ [0,\Sigma] \times [0,R] \backslash \Omega( \varepsilon) $. Hence, it attains its minimal value $ \text{C}_{0} > 0 $. Consequently, for all $ s \in [0,\Sigma] $, if one denotes by $ I_{s} : [0,\breve{\rho}_{s}(q_ \varepsilon (s))] \to [0,q_{ \varepsilon}(s)] $ the increasing part of $ F(s,\cdot) $, this function carries out a $ \mathbf{C}^{{1}} $-diffeomorphism. Moreover,

    $ \forall q \in [0,q_{ \varepsilon}(s)], \ \left|(I_{s}^{-1})^{'}(q) \right| \leq \frac{1}{\text{C}_{0}}. $

    Then, for all $ k \in \{0, \ldots, n\} $,

    $ |˘ρsk+1(qk+1)˘ρsk(qk)|=|(I1sk+1)(qk+1)˘ρsk(qk)|1C0|qk+1qk|+|(I1sk+1)(qk)˘ρsk(qk)|=1C0|qk+1qk|+|(I1sk+1)(qk)(I1sk+1)  Isk+1(˘ρsk(qk))|1C0(|qk+1qk|+|qkIsk+1(˘ρsk(qk))|)=1C0(|qk+1qk|+|F(sk,˘ρsk(qk))F(sk+1,˘ρsk(qk))|)1C0(|qk+1qk|+sFL|sk+1sk|)1+sFLC0(|qk+1qk|+|sk+1sk|). $

    Using the same techniques, one can show that the same inequality holds when considering $ {{\left| \widehat{\rho}_{s^{k+1}}(q^{k+1}) - \widehat{\rho}_{s^{k}}(q^{k}) \right|}} $. Therefore, inequality (60) follows with

    $ \mathbf{C}^{{ \varepsilon}} = 4 \times \left( \frac{1 + \| \partial_s F \|_{\mathbf{L}^{{\infty}}}}{\text{C}_{0}} \right). $

    Remark 14. Recall we suppose that $ F : [0,\Sigma] \times [0,R] $ is continuously differentiable, but if we look in the details of the proof above, we actually need $ F = F(s,\rho) $ to be continuously differentiable with respect to $ s $ and

    $ \forall s \in [0,\Sigma], \ F(s,\cdot) \in \mathbf{C}^{{1}}([0,R] \backslash \{ \overline{\rho}_s\}), \quad \overline{\rho}_s = \underset{\rho \in [0,R]}{\text{argmax}} \ F(s,\rho). $

    Corollary 1. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $, $ s \in \mathbf{BV}([0,T], [0,\Sigma]) $ and $ q \in \mathbf{BV}([0,T], \mathbb{R}^+) $. Suppose that $ q $ verifies Assumption (59). Let $ \rho_{\Delta} = \left( \rho^{n}_{j+1/2} \right)_{n,j} $ be the finite volume approximate solution constructed with the scheme (56)-(57), using the Godunov numerical flux when $ j = 0 $ in (57), and any other monotone consistent and Lipschitz numerical flux when $ j \neq 0 $. Then there exists $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $ such that

    $ \forall t \in [0,T], \ \rho_{\Delta}(t) {\ \underset{{\Delta} \to {0}}{\longrightarrow} \ } \rho(t) \ {\rm{in}} \ \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R}). $

    Proof. Since $ s $ and $ q $ have bounded variation, inequality (60) leads to an uniform in time $ \mathbf{BV} $ bound for the sequence $ \left( \rho_{\Delta} \right)_\Delta $. Then the result from [22,Appendix] establish the compactness statement.

    Theorem A.5. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $, $ s \in \mathbf{BV}([0,T] ; [0,\Sigma]) $, $ F \in \mathbf{C}^{{1}}([0,\Sigma] \times [0,R]) $ verifying (53) and $ q \in \mathbf{BV}([0,T] ; \mathbb{R}^{+}) $. Suppose that in (57), we use the Godunov flux when $ j = 0 $ and any other monotone consistent and Lipschitz numerical flux when $ j \neq 0 $. Finally, suppose that $ q $ satisfies (59). Then under the CFL condition (58), the scheme (56)-(57) converges to an admissible weak solution $ \rho $, to (52), which is also $ \mathbf{BV} $-regular. More precisely, there exists a constant $ \mathbf{C}^{{ \varepsilon}} = \mathbf{C}^{{ \varepsilon}}(\| \partial_s F \|_{\mathbf{L}^{{\infty}}}) $ nondecreasing with respect to its argument such that

    $ t[0,T], TV(ρ(t))TV(ρ0)+4R+Cε(TV(q)+TV(s)). $ (61)

    Proof. From the scheme (56), one can derive approximate entropy/constraint inequalities analogous to (30)-(31) of Section 4. Let $ \rho $ be the limit to the finite volume scheme, the compactness of $ \left( \rho_{\Delta} \right)_\Delta $ coming from the last corollary. We already know that $ \rho \in \mathbf{C}^{{}}([0,T] ; \mathbf{L}_{\mathbf{loc}}^{{1}}( \mathbb{R})) $. By passing to the limit in the approximate entropy/constraint inequalities verified by $ \left( \rho_{\Delta} \right)_\Delta $ we get that $ \rho $ satisfies the entropy/constraint inequalities of Definition A.1. This shows that $ \rho $ is an admissible weak solution to Problem (52). Finally, from (60), the lower semi-continuity of the $ \mathbf{BV} $ semi-norm ensures that $ \rho \in \mathbf{L}^{{\infty}} ([0,T] ; \mathbf{BV}( \mathbb{R})) $ and verifies (61). This concludes the proof.

    Corollary 2. Fix $ \rho_{0} \in \mathbf{L}^{{1}}( \mathbb{R} ; [0,R]) \cap \mathbf{BV}( \mathbb{R}) $, $ s \in \mathbf{BV}([0,T] ; [0,\Sigma]) $, $ F \in \mathbf{C}^{{1}}([0,\Sigma] \times [0,R]) $ verifying (53) and $ q \in \mathbf{BV}([0,T] ; \mathbb{R}^{+}) $. Suppose that $ q $ satisfies Assumption (59). Then Problem (52) admits a unique $ \mathbf{BV} $-regular solution $ \rho $. Moreover, $ \rho $ satisfies the bound (61).

    Proof. Uniqueness comes from Theorem 54, the existence and the $ \mathbf{BV} $ bound comes from Theorem A.5.

    Remark 15. Under the hypotheses of Corollary 2, if we prove the existence of an other admissible weak solution $ \overline{\rho} $ to (11) (by another method, splitting for instance), then Theorem A.3 ensures that $ \overline{\rho} = \rho $.

    [1] Goodrich JA,Tjian R (1994) TBP-TAF complexes: selectivity factors for eukaryotic transcription. Curr Opin Cell Biol 6: 403-409. doi: 10.1016/0955-0674(94)90033-7
    [2] Buratowski S (1994) The basics of basal transcription by RNA polymerase II. Cell 77: 1-3. doi: 10.1016/0092-8674(94)90226-7
    [3] Evans RM (1998) The steroid and thyroid hormone receptor superfamily. Science 240: 889-895.
    [4] Dasgupta AB,Lonard DM,O'Malley BW (2014) Nuclear receptor coactivators: master regulators of human health and disease. Annu Rev Med 65: 279-292. doi: 10.1146/annurev-med-051812-145316
    [5] Xu J,Li Q (2003) Review of the in vivo Functions of the p160 Steroid Receptor Coactivator Family. Mol Endocrinol 17:1681-1692. doi: 10.1210/me.2003-0116
    [6] Smith CL,DeVera DG,Lamb DJ,et al. (2002) Genetic ablation of the steroid receptor coactivator-ubiquitin ligase,E6-AP,results in tissue-selective steroid hormone resistance and defects in reproduction. Mol Cell Biol 22: 525-535. doi: 10.1128/MCB.22.2.525-535.2002
    [7] Beitel LK,Elhaji YA,Lumbroso R,et al. (2002) Cloning and characterization of an androgen receptor N-terminal-interacting protein with ubiquitin-protein ligase activity. J Mol Endocrinol 29: 41-60. doi: 10.1677/jme.0.0290041
    [8] Khan SH,Ling J,Kumar R (2011) TBP binding-induced folding of the glucocorticoid receptor AF1 domain facilitates its interaction with steroid receptor coactivator-1. PLoS One 6: e21939
    [9] Kumar R,Moure CM,Khan SH,et al. (2013) Regulation of the structurally dynamic N-terminal domain of progesterone receptor by protein-induced folding. J Biol 288: 30285-3099.
    [10] Warnmark A,Wikstrom A,Wright AP,et al. (2001) The N-terminal regions of estrogen receptor alpha and beta are unstructured in vitro and show different TBP binding properties. J Biol Chem 276: 45939-45944. doi: 10.1074/jbc.M107875200
    [11] Keaveney M,Berkenstam A,Feigenbutz M,et al. (1993) Residues in the TATA-binding protein required to mediate a transcriptional response to retinoic acid in EC cells. Nature 365: 562-566. doi: 10.1038/365562a0
    [12] Pugh BF,Tjian R (1990) Mechanism of transcriptional activation by Sp1: evidence for coactivators. Cell 61: 1187-1197. doi: 10.1016/0092-8674(90)90683-6
    [13] Strubin M,Struhl K (1992) Yeast and human TFIID with altered DNA-binding specificity for TATA elements. Cell 68:721-730. doi: 10.1016/0092-8674(92)90147-5
    [14] Nordeen SK,Moyer ML,Bona BJ (1995) Modulation of glucocorticoid-regulated transcription by purines: novel characteristics and implications for tissue specificity of steroid responses. Endocrinology 136: 1120-1127.
    [15] Beck CA,Weigel NL,Moyer ML,et al. (1993) The progesterone antagonist RU486 acquires agonist activity upon stimulation of cAMP signaling pathways. Proc Natl Acad Sci U S A 90: 4441-4445. doi: 10.1073/pnas.90.10.4441
    [16] Nordeen SK,Bona BJ,Moyer ML (1993) Latent agonist activity of the steroid antagonist,RU486,is unmasked in cells treated with activators of protein kinase A. Mol Endocrinol 7: 731-742.
    [17] Moyer ML,Borror KC,Bona BJ,et al. (1993) Modulation of cell signaling pathways can enhance or impair glucocorticoid-induced gene expression without altering the state of receptor phosphorylation. J Biol Chem 268: 22933-22940.
    [18] Nordeen SK,Moyer ML,Bona BJ (1994) The coupling of multiple signal transduction pathways with steroid response mechanisms. Endocrinology 134: 1723-1732.
    [19] Nordeen SK (1994) Potentiation of glucocorticoid-mediated induction of the mouse mammary tumor virus promoter in mammary carcinoma cells by activation of protein kinase C. Biochem Biophys Res Commun 198: 1183-1188. doi: 10.1006/bbrc.1994.1167
    [20] Dutertre M,Smith CL (2003) Ligand-independent interactions of p160/steroid receptor coactivators and CREB-binding protein (CBP) with estrogen receptor-alpha: regulation by phosphorylation sites in the A/B region depends on other receptor domains. Mol Endocrinol 17: 1296-1314. doi: 10.1210/me.2001-0316
    [21] Wu RC,Smith CL,O'Malley BW (2005) Transcriptional regulation by steroid receptor coactivator phosphorylation. Endocr Rev 26: 393-399. doi: 10.1210/er.2004-0018
    [22] Kato S,Masuhiro Y,Watanabe M (2000) Molecular mechanism of a cross-talk between oestrogen and growth factor signaling pathways. Genes Cells 5: 593-601. doi: 10.1046/j.1365-2443.2000.00354.x
    [23] Kit S,Dubbs DR,Piekarksi LJ,et al. (1963) Deletion of thymidine kinase activity from L cells resistant to bromodeoxyuridine. Exp Cell Res 31: 297-312. doi: 10.1016/0014-4827(63)90007-7
    [24] Ralph P,Moore MA,Milsson K (1976) Lysozyme synthesis by established human and murine histiocytic lymphoma cell lines. J Exp Med 143: 1528-1533. doi: 10.1084/jem.143.6.1528
    [25] Gatlin J,Campbell LH,Schmidt MG,et al. (1995) Direct-rapid (DR) mutagenesis of large plasmids using PCR. Biotechniques 19: 559-564.
    [26] Nordeen SK (1988) Luciferase reporter gene vectors for analysis of promoters and enhancers. Biotechniques 6: 454-458.
    [27] Lieberman BA,Bona BJ,Edwards DP,et al. (1993) The constitution of a progesterone response element. Mol Endocrinol 7: 515-527.
    [28] Biggs JR,Ahn NG,Kraft AS (1998) Activation of the mitogen-activated protein kinase pathway in U937 leukemic cells induces phosphorylation of the amino terminus of the TATA-binding protein. Cell Growth Differ 9: 667-676.
    [29] Fujimoto N,Katzenellenbogen BS (1994) Alteration in the agonist/antagonist balance of antiestrogens by activation of protein kinase A signaling pathways in breast cancer cells: antiestrogen selectivity and promoter dependence. Mol Endocrinol 8: 296-304.
    [30] Keenan CR,Mok JS,Harris T,et al (2014) Bronchial epithelial cells are rendered insensitive to glucocorticoid transactivation by transforming growth factor-β1. Respiratory Res 15: 55
    [31] Piovan E,Yu J,Tosello V,et al. (2013) Direct reversal of glucocorticoid resistance by AKT inhibition in acute lymphoblastic leukemia. Cancer Cell 24: 766-776. doi: 10.1016/j.ccr.2013.10.022
    [32] Beck IM,Vanden Berghe W,Vermeulen L (2009) Crosstalk in inflammation: the interplay of glucocorticoid receptor-based mechanisms and kinases and phosphatases. Endocr Rev 30: 830-882. doi: 10.1210/er.2009-0013
    [33] Kumar R,Calhoun WJ (2008) Differential regulation of the transcriptional activity of the glucocorticoid receptor through site-specific phosphorylation. Biologics 2: 845-854.
    [34] Ward RD,Weigel NL (2009) Steroid receptor phosphorylation: Assigning function to site-specific phosphorylation. Biofactors 35: 528-536. doi: 10.1002/biof.66
    [35] Chen W,Dang T,Blind RD,et al. (2008) Glucocorticoid receptor phosphorylation differentially affects target gene expression. Mol Endocrinol 22: 1754-66. doi: 10.1210/me.2007-0219
    [36] Espinas ML,Roux J,Pictet R,et al. (1995) Glucocorticoids and protein kinase A coordinately modulate transcription factor recruitment at a glucocorticoid-responsive unit. Mol Cell Biol 15: 5346-5354.
    [37] Rowan BG,Garrison N,Weigel NL,et al. (2000) 8-Bromo-cyclic AMP induces phosphorylation of two sites in SRC-1 that facilitate ligand-independent activation of the chicken progesterone receptor and are critical for functional cooperation between SRC-1 and CREB binding protein. Mol Cell Biol 20: 8720-8730. doi: 10.1128/MCB.20.23.8720-8730.2000
    [38] Johnson AB,O'Malley BW (2012) Steroid receptor coactivators 1,2,and 3: critical regulators of nuclear receptor activity and steroid receptor modulator (SRM)-based cancer therapy. Mol Cell Endocrinol 348: 430-439. doi: 10.1016/j.mce.2011.04.021
    [39] Das D,Scovell WM (2001) The binding interaction of HMG-1 with the TATA-binding protein/TATA complex. J Biol Chem 276: 32597-32605. doi: 10.1074/jbc.M011792200
    [40] Boonyaratanakornkit V,Melvin V,Prendergast P,et al. (1998) High-mobility group chromatin proteins 1 and 2 functionally interact with steroid hormone receptors to enhance their DNA binding in vitro and transcriptional activity in mammalian cells. Mol Cell Biol 18: 4471-4487.
    [41] Melvin VS,Edwards DP (1999) Coregulatory proteins in steroid hormone receptor action: the role of chromatin high mobility group proteins HMG-1 and -2. Steroids 64: 576-586. doi: 10.1016/S0039-128X(99)00036-7
    [42] Hernandez N (1993) TBP,a universal eukaryotic transcription factor? Genes Dev 7:1291-1308. doi: 10.1101/gad.7.7b.1291
    [43] Mittal V,Hernandez N (1997) Role for the amino-terminal region of human TBP in U6 snRNA transcription. Science 275: 1136-1140. doi: 10.1126/science.275.5303.1136
    [44] Mittal V,Ma B,Hernandez N (1999) SNAP(c): a core promoter factor with a built-in DNA-binding damper that is deactivated by the Oct-1 POU domain. Genes Dev 13: 1807-1821. doi: 10.1101/gad.13.14.1807
    [45] Zhao X,Herr W (2002) A regulated two-step mechanism of TBP binding to DNA: a solvent-exposed surface of TBP inhibits TATA box recognition. Cell 108: 615-627. doi: 10.1016/S0092-8674(02)00648-7
    [46] Um M,Yamauchi J,Kato S,et al. (2001) Heterozygous disruption of the TATA-binding protein gene in DT40 cells causes reduced cdc25B phosphatase expression and delayed mitosis. Mol Cell Biol 21: 2435-2448. doi: 10.1128/MCB.21.7.2435-2448.2001
    [47] Schmidt EE,Bondareva AA,Radke JR,et al. (2003) Fundamental cellular processes do not require vertebrate-specific sequences within the TATA-binding protein. J Biol Chem 278: 6168-6174. doi: 10.1074/jbc.M211205200
    [48] Hobbs NK,Bondareva AA,Barnett S,et al. (2002) Removing the vertebrate-specific TBP N terminus disrupts placental beta2m-dependent interactions with the maternal immune system. Cell 110: 43-54. doi: 10.1016/S0092-8674(02)00806-1
  • This article has been cited by:

    1. Paola Goatin, Chiara Daini, Maria Laura Delle Monache, Antonella Ferrara, Interacting moving bottlenecks in traffic flow, 2023, 18, 1556-1801, 930, 10.3934/nhm.2023040
    2. Abraham Sylla, A LWR model with constraints at moving interfaces, 2022, 56, 2822-7840, 1081, 10.1051/m2an/2022030
    3. Boris Andreianov, Abraham Sylla, Existence analysis and numerical approximation for a second-order model of traffic with orderliness marker, 2022, 32, 0218-2025, 1295, 10.1142/S0218202522500294
    4. Boris Andreianov, Carlotta Donadello, Massimiliano D. Rosini, Entropy solutions for a two-phase transition model for vehicular traffic with metastable phase and time depending point constraint on the density flow, 2021, 28, 1021-9722, 10.1007/s00030-021-00689-5
    5. Boris Andreianov, Abraham Sylla, Finite volume approximation and well-posedness of conservation laws with moving interfaces under abstract coupling conditions, 2023, 30, 1021-9722, 10.1007/s00030-023-00857-9
    6. Boris Andreianov, Massimiliano D. Rosini, Graziano Stivaletta, On existence, stability and many-particle approximation of solutions of 1D Hughes' model with linear costs, 2023, 369, 00220396, 253, 10.1016/j.jde.2023.06.004
    7. Boris Andreianov, Theo Girard, Existence of Solutions for a Class of One-Dimensional Models of Pedestrian Evacuations, 2024, 56, 0036-1410, 3386, 10.1137/23M1550256
    8. Malaya Mohanty, Madhu Lisha Pattanaik, Satya Ranjan Samal, Shivam Shivam, Assessing the Impact of Three-Wheelers on Traffic Flow in India: A Case Study Using Ann, 2024, 13354205, 10.26552/com.C.2024.033
  • Reader Comments
  • © 2015 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(4115) PDF downloads(965) Cited by(0)

Figures and Tables

Figures(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog