Novel model for integrated demand-responsive transit service considering rail transit schedule

1.   Introduction

Partial differential equations on networks have been considered in the last years by several authors, in particular in the parabolic case; we quote for instance [8,10,11,16,24,30]. According to the modeling in consideration and to the type of equations on the edges of the underlying graph, different conditions at the nodes are imposed. In most of the cases, precise results of existence of solutions are given, even for rather complicated networks.

In this paper, the main example we have in mind comes from traffic modeling, where the network is constituted by a crossroad connecting $m$ incoming roads with $n$ outgoing roads; the traffic in each road is modeled by the scalar diffusive equation

where $t$ denotes time and $x$ the position along the road. In this case $\rho_h$ is a vehicle density; about the diffusivity $D_h(\rho_h)\ge0$ we do not exclude that it may vanish at some points. System (1.1) is completed by a condition of flux conservation at the crossroad, which implies the conservation of the total number of cars. Such a model is derived from the famous Lighthill-Whitham-Richards equation [17,25]. We refer to [3,15,17,20,22,27] for several motivations about the introduction of (possibly degenerate) diffusion in traffic flows and in the close field of crowds dynamics. We also refer to the recent books [10,11,26] for more information on the related hyperbolic modeling.

We focus on a special class of solutions to (1.1), namely, traveling waves. In the case of a single road, traveling waves are considered, for instance, in [21]; in the case of a second-order model without diffusion but including a relaxation term, we refer to [9,28]; for a possibly degenerate diffusion function and in presence of a source term, detailed results are given in [6,7]. In the case of a network, the papers dealing with this subject, to the best of our knowledge, are limited to [30,31] for the semilinear diffusive case and to [19] for the case of a dispersive equation. In these papers, as in most modeling of diffusive or dispersive partial differential equations on networks, both the continuity of the unknown functions and the Kirchhoff condition (or variants of it) are imposed at the nodes. We emphasize that while the classical Kirchhoff condition implies the conservation of the flow and then that of the mass, some variants of this condition are dissipative and, then, imply none of the conservations above. While these assumptions are natural when dealing with heat or fluid flows, they are much less justified in the case of traffic modeling, where the density must be allowed to jump at the node while the conservation of the mass must always hold. Moreover, they impose rather strong conditions on the existence of the profiles, which often amount to proportionality assumptions on the parameters in play.

In this paper we only require the conservation of the (parabolic) flux at the node, as in [4]; differently from that paper and the other ones quoted above, we do not impose the continuity condition. A strong motivation for dropping this condition comes from the hyperbolic modeling [1,10,11,26]; nevertheless, we show how our results simplify when such a condition is required. In particular, in Sections 6 and 7 we provide explicit conditions for traveling wave solutions which do not satisfy the continuity condition; in some other cases, such a condition is indeed always satisfied. Our main results are essentially of algebraic nature and concern conditions about the end states, flux functions, diffusivities and other parameters which give rise to a traveling wave moving in the network.

Here follows a plan of the paper. In Section 2 we introduce the model and give some basic definitions; for simplicity we only focus on the case of a star graph. Section 3 deals with a general existence result in the case of a single equation; its proof is provided in Appendix A. Section 4 contains our main theoretical results about traveling waves in a network. In that section we characterize both stationary/non-stationary and degenerate/non-degenerate waves; in particular, Theorem 4.12 contains an important necessary and sufficient condition that we exploit in the following sections. Section 5 focus on the continuity condition; in this case the conditions for the existence of traveling wave solutions are much stricter than in the previous case. Detailed applications of these results are provided in Sections 6 for quadratic fluxes and in Section 7 for logarithmic fluxes; in particular, in subsection 6.2 and in the whole Section 7 the diffusivity is as in [3]. For simplicity, we only deal there with the case of a single ingoing road but we consider both constant and degenerate diffusivities.

From a theoretical point of view, the extension of our work to more general networks is possible by applying the results in Section 4 to each node iteratively. However, a quick look at Sections 6 and 7 shows that the corresponding results become quickly very technical. In particular, in the general case the whole problem must be recast by using the graph notation, see for instance [18].

2.   The model

In terms of graph theory, we consider a semi-infinite star-graph with $m$ incoming and $n$ outgoing edges; this means that the incidence vector $d\in{\mathbb R}^{m+n}$ has components $d_i=1$ for $i\in\mathsf{I} \doteq \{ 1, \ldots, m \}$ and $d_j=-1$ for $j\in\mathsf{J} \doteq \{ m+1, \ldots, m+n \}$. We also denote $\mathsf{H} \doteq \{ 1,\ldots, m+n \}$ and refer to Figure 1. For simplicity, having in mind the example in the Introduction, we always refer to the graph as the network, to the node as the crossroad and to the edges as the roads. Then, incoming roads are parametrized by $x \in {\mathbb R}_-\doteq(-\infty,0]$ and numbered by the index $i$, outgoing roads by $x \in {\mathbb R}_+\doteq[0,\infty)$ and $j$; the crossroad is located at $x=0$ for both parameterizations. We denote the generic road by $\Omega_h$ for $h \in \mathsf{H}$; then $\Omega_i \doteq {\mathbb R}_-$ for $i \in \mathsf{I}$ and $\Omega_j \doteq {\mathbb R}_+$ for $j \in \mathsf{J}$. The network is defined as $\mathcal{N}\doteq \prod_{h\in\mathsf{H}}\Omega_h$.

Following the above analogy, we understand the unknown functions $\rho_h$ as vehicular densities in the roads $\Omega_h$, $h\in\mathsf{H}$; $\rho_h$ ranges in $[0,\overline{\rho}_h]$, where $\overline{\rho}_h$ is the maximal density in the road $\Omega_h$. Without loss of generality we assume that $\overline{\rho}_h=1$ for every $h\in\mathsf{H}$; the general case is easily recovered by a change of variables and modifying (2.2)-(2.3) below for a multiplicative constant. With a slight abuse of notation we denote $\rho \doteq (\rho_1, \ldots, \rho_{m+n}): {\mathbb R} \times \mathcal{N} \to [0,1]^{m+n}$ understanding that $\rho(t,x_1,\ldots, x_{m+n})=(\rho_1(t,x_1),\ldots,\rho_{m+n}(t,x_{m+n}))$.

For each road we assign the functions $f_h$, the hyperbolic flux, and $D_h$, the diffusivity; we assume for every $h \in \mathsf{H}$

(f) $f_h\in\mathbf{C}^1([0,1];{\mathbb R}_+)$ is strictly concave with $f_h(0)=f_h(1)=0$;

(D) $D_h\in\mathbf{C}^1([0,1];{\mathbb R}_+)$ and $D_h(\rho)>0$ for any $\rho \in (0,1)$.

We emphasize that in (D) we can possibly have either $D_h(0)=0$ or $D_h(1)=0$, or even both possibilities at the same time. The evolution of the flow is described by the equations

Assumption (f) is standard when dealing with traffic flows [2]. More precisely, in that case $f_h(\rho_h)=\rho_h \, V_h(\rho_h)$, where $V_h$ is the velocity. Then, assumption (f) is satisfied if, for instance, $V_h\in\mathbf{C}^2([0,1];{\mathbb R}_+)$ is either linear or strictly concave, decreasing and satisfying $V_h(1) = 0$, see [17,25]. The prototype of such a velocity satisfying (f) is $V_h(\rho) = v_h (1-\rho)$ with $v_h > 0$, which was introduced in [14]; another example is given in [23]. The simplest model for the diffusivity is then $D_h(\rho_h) = -\delta_h\rho_h \, V_h'(\rho_h)$, where $\delta_h$ is an anticipation length [3,21].

The coupling among the differential equations in (2.1) occurs by means of suitable conditions at the crossroad. In this paper, having in mind the previous example, we impose a condition on the conservation of the total flow at the crossroad, see [4,5]; in turn, this implies the conservation of the mass. More precisely, we define the parabolic flux by

and require

for given constant coefficients $\alpha_{i,j} \in (0, 1]$ satisfying

Conditions (2.2) and (2.3) imply

which is the conservation of the total flow at the crossroad. Conditions (2.2) and (2.3) deserve some comments. First, by no means they imply

Condition (2.5) is largely used, together with some Kirchhoff conditions, when dealing with parabolic equations in networks and takes the name of continuity condition. Second, above we assumed $\alpha_{i,j}>0$ for every $i$ and $j$. The case when $\alpha_{i,j}=0$ for some $i$ and $j$ would take into account the possibility that some outgoing $j$ roads are not allowed to vehicles coming from some incoming $i$ roads; this could be the case, for instance, if only trucks are allowed in road $i$ but only cars are allowed in road $j$. For simplicity, we do not consider this possibility. Third, we notice that assumption (2.2) destroys the symmetry of condition (2.4); indeed, with reference to the example of traffic flow, the loss of symmetry is due to the fact that all velocities $V_h$ are positive.

We point out that condition (2.3) is almost never explicitly exploited in Section 4 and in most of the following: our results hold for every choice of the coefficients $\alpha_{i,j}>0$. Clearly, condition (2.3) must be imposed to have mathematically (see [11,Lemma 5.1.9]) and physically meaningful solutions; however, this requirement only adds algebraic conditions on the parameters we are dealing with (the end states, the ratios of the traveling-wave speeds, the proportionality coefficients -see (6.11), for instance -and so on) which do not affect the final results.

Then, we are faced with the system of equations (2.1) that are coupled through (2.2), with the $\alpha_{i,j}$ satisfying (2.3). Solutions to (2.1)-(2.2) are meant in the weak sense, namely $\rho_h\in\mathbf{C^1}({\mathbb R}\times\Omega_h;[0,1])$ a.e.; see also [2,11] for an analogous definition in the hyperbolic case. We do not impose any initial condition because we only consider traveling waves, which are introduced in the next sections.

3.   Traveling waves for a single equation

In this section we briefly remind some definitions and results about traveling waves [12] for the single equation

where we keep for future reference the index $h$. Equation (3.1) has no source terms (differently from [30,31]) and then any constant is a solution; for simplicity we discard constant solutions in the following analysis.

Definition 3.1. A weak solution $\rho_h(t,x)$ to (3.1) is a traveling-wave solution of (3.1) if $\rho_h(t,x) = \varphi_h(x-c_h t)$ for $(t,x) \in \mathbb{R} \times \Omega_h$, for a non-constant profile $\varphi_h:{\mathbb R}\to[0,1]$ and speed $c_h \in {\mathbb R}$.

This definition coincides with that given in [19,29] because we are considering non-constant profiles. The profile must satisfy the equation

namely,

in the weak sense, where

is the reduced flux, see Figure 2.

This means that $\varphi_h \in \mathbf{C^{0}}(\mathbb{R}; [0,1])$, $D_h(\varphi_h) \, \varphi_h'\in \mathbf{L^{1}_{loc}}(\mathbb{R};\mathbb{R})$ and

for every $\psi\in \mathbf{C_c^{\infty}}(\mathbb{R};\mathbb{R})$. Equation (3.3) is coupled with the limit conditions

for $\ell_h^\pm\in[0,1]$. Clearly, solutions to (3.3) -(3.5) are determined up to a shift. We define

The existence of profiles is a well-established result [12]; nevertheless, we state for completeness the following theorem, where we point out the qualitative properties of these fronts. The proof is deferred to Appendix A.

Theorem 3.2. Assume (f) and (D). Equation (3.1) admits a traveling-wave solution $\rho_h$ with profile $\varphi_h$ satisfying (3.5) if and only if

We have that $\varphi_h \in \mathbf{C}^2\left(I_h; (\ell_h^-, \ell_h^+)\right)$ is unique (up to shifts) and $\varphi_h'(\xi)>0$ for $\xi \in I_h$; moreover, the following holds true.

(i) $D_h(0)=0=\ell_h^-$ if and only if there exists $\nu_h^- \in \mathbb{R}$ such that $I_h \subseteq (\nu_h^-,\infty)$ and $\varphi_h(\xi)=0$ for $\xi \le \nu_h^-$. In this case

(ii) $D_h(1)=0=1-\ell_h^+$ if and only if there exists $\nu_h^+ \in \mathbb{R}$ such that $I_h \subseteq (-\infty,\nu_h^+)$ and $\varphi_h(\xi)=1$ for $\xi \ge \nu_h^+$. In this case

(iii) In all the other cases $I_h=\mathbb{R}$ and

We observe that for $c_h$ given by (3.7), we deduce by (f) that $g_h(\rho) \ge 0$ for all $\rho \in[\ell_h^-,\ell^+_h]$, see Figure 2. Moreover, we have

and no $\rho\ne\ell_h^\pm$ makes $g_h(\rho)$ equal to that value.

Theorem 3.2 motivates the following definition.

Definition 3.3. A traveling-wave solution $\rho_h$ is stationary if $c_h=0$. It is degenerate if at least one of conditions (ⅰ) or (ⅱ) of Theorem 3.2 holds.

Remark 3.4. A consequence of assumption (f) is that if $\rho_h$ is degenerate, then the profile $\varphi_h$ is singular either at $\nu_h^-$ in case (ⅰ) or at $\nu_h^+$ in case (ⅱ), in the sense that $\varphi_h'$ cannot be extended to the whole of ${\mathbb R}$ as a continuous function.

In case (i) (or (ii)) of Theorem 3.2 does not hold we define $\nu_h^- \doteq -\infty$ (respectively, $\nu_h^+ \doteq \infty$). In this way the interval $(\nu_h^-,\nu_h^+)$ is always defined and coincides with the interval $I_h$ defined in (3.6):

The interval $I_h$ is bounded if and only if both (i) and (ii) hold; in this case $\rho_h$ is both degenerate and stationary. As a consequence, if $\rho_h$ is non-stationary then $I_h$ is unbounded and coincides either with a half line (if $\rho_h$ is degenerate) or with ${\mathbb R}$ (if $\rho_h$ is non-degenerate). At last, $\rho_h$ is degenerate if and only if either $\nu_h^-$ or $\nu_h^+$ is finite.

In the case of non-stationary traveling-wave solutions $\rho_h$ we use the notation

Lemma 3.5. Let $\rho_h$ be a traveling-wave solution of (3.1); then we have the following.

(a) If $\rho_h$ is stationary, then it is degenerate if and only if $D_h(0)D_h(1)=0$ and $\ell_h^- = 0$ (hence $\ell_h^+=1$).

(b) If $\rho_h$ is non-stationary, then it is degenerate if and only if one of the following equivalent statements hold:

• either $D_h(0) = 0 = \ell_h^-$ or $D_h(1) = 0 = 1-\ell_h^+$, but not both;

• $\omega_h$ is finite.

In this case the function $\xi \mapsto \varphi_h'(c_h\xi)$ is singular at $\xi = \omega_h$ and $\mathbf{C^{1}}$ elsewhere.

Proof. We recall that $\rho_h$ is degenerate if and only if either $D_h(0) = 0 = \ell_h^-$ or $D_h(1) = 0 = 1-\ell_h^+$. This means that at least one of the end states must be $0$ or $1$, say $0$; but then $c_h=0$ if and only if the other end state is $1$. This proves (a) and the first part of (b).

Now, we prove the second part of (b). Since $c_h \ne 0$, $\rho_h$ is degenerate if and only if exactly one between (i) and (ii) of Theorem 3.2 occurs, namely, exactly one between $\nu_h^-$ and $\nu_h^+$ is finite. If $\nu_h^-$ is finite and $\nu_h^+=\infty$, then $c_h=f_h(\ell_h^+)/\ell_h^+>0$ and $\omega_h = c_h^{-1}\nu_h^-$ is finite. By Remark 3.4, we know that $\xi\mapsto\varphi_h'(\xi)$ is singular at $\xi=\nu_h^-$ and $\mathbf{C^{1}}$ elsewhere, whence the regularity of $\xi \mapsto \varphi_h'(c_h\xi)$. Analogously, if $\nu_h^+$ is finite and $\nu_h^-=-\infty$, then $c_h=-f_h(\ell_h^-)/(1-\ell_h^-)<0$ and $\omega_h = c_h^{-1}\nu_h^+$ is finite. The statement about the smoothness of $\xi \mapsto \varphi_h'(c_h\xi)$ is proved as above.

Finally, the converse is straightforward. In fact, if $\omega_h$ is finite, then either $\omega_h = c_h^{-1}\nu_h^-$ and $\nu_h^-$ is finite, or $\omega_h = c_h^{-1}\nu_h^+$ and $\nu_h^+$ is finite; in both cases $\rho_h$ is degenerate.

Because of the smoothness properties of the profile proved in Theorem 3.2, we can integrate equation (3.2) in $(\xi_-,\xi)\subset I_h$ and we obtain

If $\xi_-\downarrow\nu_h^-$ in the previous expression, by applying (3.9) or (3.12) we deduce

We observe that (3.15) is trivially satisfied in case (i) when $\xi < \nu_h^-$ and in case (ii) when $\xi > \nu_h^+$; moreover, by a continuity argument, we deduce from (3.9) and (3.11) that (3.15) is satisfied in case (i) at $\xi = \nu_h^-$ and in case (ii) at $\xi = \nu_h^+$, respectively. In conclusion, we have that (3.15) holds in the whole ${\mathbb R}$, namely

4.   Traveling waves in a network

In this section we consider the traveling-wave solutions of problem (2.1)-(2.2) in the network $\mathcal{N}$. We first introduce the definition of traveling-wave solution in $\mathcal{N}$.

Definition 4.1. For any $h \in \mathsf{H}$, let $\rho_h$ be a traveling-wave solution of $(2.1)_h$ in the sense of Definition 3.1 and set $\rho \doteq (\rho_1, \ldots, \rho_{m+n})$. With reference to Definition 3.3, we say that:

• $\rho$ is stationary if each component $\rho_h$ is stationary;

• $\rho$ is completely non-stationary if none of its components is stationary;

• $\rho$ is degenerate if at least one component $\rho_h$ is degenerate;

• $\rho$ is completely degenerate if each of its components is degenerate.

Finally, we say that $\rho$ is a traveling-wave solution of problem (2.1)-(2.2) in the network $\mathcal{N}$ if (2.2) holds.

For brevity, from now on we simply write "traveling wave" for "traveling-wave solution". In analogy to the notation above, we say that $\varphi \doteq (\varphi_1, \ldots, \varphi_{m+n})$ is a profile for $\rho$ if $\varphi_h$ is a profile corresponding to $\rho_h$ for every $h \in \mathsf{H}$.

For clarity of exposition, we collect our general results for stationary and non-stationary traveling waves in the following subsections.

4.1.   General results

In this subsection, as well as in the following ones, we always assume (f) and (D) without explicitly mentioning it. Moreover, by Definition 4.1 and Theorem 3.2, the end states and the speeds of the profiles must satisfy (3.7) for every $h\in\mathsf{H}$; both conditions in (3.7) are tacitly assumed as well.

Proposition 4.2. The function $\varphi$ is the profile of a traveling wave for (2.1)-(2.2) if and only if $\varphi_h$ is a solution to (3.5)-(3.16) for any $h \in \mathsf{H}$ and

In (4.1) any combination of the signs $\pm$ is allowed.

Proof. By plugging $\rho_h(t,x)=\varphi_h(x-c_ht)$ in (2.2) and recalling that by Theorem 3.2 the profiles are continuous in ${\mathbb R}$, we obtain

which is equivalent to (4.1) by (3.16). At last, we can clearly choose any combination of signs in (4.1) because of (3.13).

Differently from what specified in Proposition 4.2, in the following the choice of the signs "$\pm$" follows the usual rules, i.e., top with top and bottom with bottom.

Lemma 4.3. Assume that problem (2.1)-(2.2) admits a traveling wave. Then for any $j\in\mathsf{J}$ we have

Proof. Fix $j\in\mathsf{J}$. We notice that (4.1) is equivalent to

where the map $t \mapsto \Upsilon_h(t) \doteq c_h \, \varphi_h(c_ht) + g_h(\ell_h^-)$ is non-decreasing because the profiles are so, by Theorem 3.2. Since we can write $\Upsilon_h(t) = f_h(\ell_h^-) + c_h [\varphi_h(c_ht) - \ell_h^-]$, we see that $\Upsilon_h$ ranges between $f_h(\ell_h^-)$ and $f_h(\ell_h^+)$ because of (3.7) and the fact that $\xi\mapsto\varphi_h(\xi)$ takes values in $[\ell_h^-,\ell_h^+]$. As a consequence,

Hence, by passing to the limit for $t\to\pm\infty$ in (4.1) we obtain (4.2) and (4.3), respectively.

Lemma 4.4. Assume that problem (2.1)-(2.2) admits a traveling wave. The traveling wave is stationary if and only if one of the following equivalent statements hold:

(i) there exists ${\rm j} \in \mathsf{J}$ such that $c_{\rm j} = 0$;

(ii) $c_i = 0$ for all $i \in \mathsf{I}$;

(iii) $c_j = 0$ for all $j \in \mathsf{J}$.

Proof. By subtracting (4.3) to (4.2) we obtain

Since $c_h=0$ if and only if $f_h(\ell_h^-) = f_h(\ell_h^+)$, from the above equation we immediately deduce that $(i)$, $(ii)$ and $(iii)$ are equivalent. By the equivalence of (ii) and (iii), a traveling wave is stationary if and only if one of the statements above holds.

Lemma 4.4 shows that either a traveling wave is stationary, and then $c_h=0$ for every $h\in\mathsf{H}$, or it is non-stationary, and then

Of course, by Lemma 4.4, $c_{\rm i}\ne0$ for some ${\rm i}\in\mathsf{I}$ if and only if $c_j\ne0$ for every $j\in\mathsf{J}$.

Proposition 4.5. Fix $\ell_i^\pm \in [0,1]$ with $\ell_i^- < \ell_i^+$, $i \in \mathsf{I}$. Then for any $j\in\mathsf{J}$ there exist $\ell_j^\pm \in [0,1]$ with $\ell_j^- < \ell_j^+$ and satisfying (4.2)-(4.3) if and only if

In this case, the end states $\ell_j^\pm$ are uniquely determined if and only if $c_i=0$ for every $i \in \mathsf{I}$.

Proof. Assume that there exist $\ell_j^\pm \in [0,1]$, with $\ell_j^- < \ell_j^+$, which satisfy (4.2)-(4.3). Then clearly we have $\max_{[0,1]}f_j \ge \sum_{i \in \mathsf{I}} \alpha_{i,j} \, \max\{f_i(\ell_i^-), f_i(\ell_i^+)\}$. If $c_i=0$ for every $i\in\mathsf{I}$, then we have $c_j=0$ for every $j\in\mathsf{J}$ by Lemma 4.4; the equality $\max_{[0,1]}f_j= f(\ell_j^-)=f(\ell_j^+)$ would imply $\ell_j^-=\ell_j^+$ because of (f), a contradiction, and then $\max_{[0,1]}f_j> f(\ell_j^-)=f(\ell_j^+)$. This proves (4.5).

Conversely, assume (4.5). If $c_i=0$ for every $i\in\mathsf{I}$, then $\ell_j^-<\ell_j^+$ are uniquely determined because of the strict concavity of $f_j$. Assume, on the contrary, that $c_{\rm i}\ne0$ for some ${\rm i} \in \mathsf{I}$; then $c_j\ne0$ by Lemma 4.4, i.e., $f_j(\ell_j^-)\ne f_j(\ell_j^+)$. Thus (4.2)-(4.3) determine exactly four possible choices of end states $\ell_j^\pm$ with $\ell_j^- < \ell_j^+$, see Figure 3.

By Proposition 4.5 and Lemma 4.4 we deduce that the end states $\ell_j^\pm$ are uniquely determined in terms of the end states $\ell_i^\pm$ if and only if the traveling wave is stationary and the first condition in (4.5) holds.

We now give an algebraic result about determining the end states of the outgoing profiles in terms of the end states of the ingoing ones. We introduce

Proposition 4.6. Assume that problem (2.1)-(2.2) admits a traveling wave. Then for any $j\in\mathsf{J}$ we have

Moreover, (4.7) is equivalent to (4.2)-(4.3).

Proof. Fix $j\in\mathsf{J}$. By Lemma 4.3 it is sufficient to prove that (4.7) is equivalent to (4.2)-(4.3). If $c_j>0$, and then $f_j(\ell_j^+) > f_j(\ell_j^-)$, by (4.6) we have

and therefore (4.7) is equivalent to (4.2)-(4.3). The case $c_j<0$ is analogous. If $c_j=0$, then $f_j(\ell_j^+) = f_j(\ell_j^-)$ and by Lemma 4.4 we have $f_i(\ell_i^+) = f_i(\ell_i^-)$ for any $i \in \mathsf{I}$. In this case formulas (4.2)-(4.3) reduce to a single equation, which coincides with (4.7).

4.2.   The stationary case

In this short subsection we briefly consider stationary traveling waves.

Theorem 4.7. Problem (2.1)-(2.2) admits infinitely many stationary traveling waves; such waves are characterized by the conditions on the end states

Proof. Clearly, (4.8) is trivially satisfied if $\ell_h^-=0$ and $\ell_h^+ = 1$ for all $h \in \mathsf{H}$. We claim that there exist infinitely many choices of $\ell_1^\pm,\ldots,\ell_{m+n}^\pm$ satisfying (4.8). To prove the claim, we choose $\ell_i^\pm \in [0,1]$, with $\ell_i^- < \ell_i^+$, such that $f_i(\ell_i^-) = f_i(\ell_i^+)$ are sufficiently small to satisfy the first condition in (4.5) for all $j \in \mathsf{J}$. Then, by a continuity argument, we can choose $\ell_j^\pm \in [0,1]$ so that $\ell_j^- < \ell_j^+$ and $f_j(\ell_j^-) = f_j(\ell_j^+) = \sum_{i \in \mathsf{I}} \alpha_{i,j} \, f_i(\ell_i^-)$. This proves the claim.

With this choice of the end states, by Theorem 3.2 we deduce the existence of a stationary traveling wave in each road satisfying (2.1). At last we notice that, in the stationary case, condition (4.1) is equivalent to the latter condition in (4.8).

Remark 4.8. The previous result does not use condition (2.3); in case it holds, then (4.8) implies

Conversely, consider the stationary case and assume conditions (4.8), (4.9); if moreover $f_i(\ell_i^-)\ne0$ and either $\sum_{j \in \mathsf{J}} \alpha_{i,j} \ge 1$ or $\sum_{j \in \mathsf{J}} \alpha_{i,j} \le 1$ for every $i\in\mathsf{I}$, then (2.3) holds. Indeed, (4.8)$_2$ and (4.9) imply

namely $\sum_{i \in \mathsf{I}} f_i(\ell_i^-) \bigl[1- \sum_{j\in\mathsf{J}}\alpha_{i,j}\bigr]=0$, whence (2.3).

Clearly, if both $D_h(0)\ne0$ and $D_h(1)\ne0$ for every $h\in\mathsf{H}$, then problem (2.1)-(2.2) admits no degenerate traveling wave. However, even in the general case, the proof of Theorem 4.7 shows that (2.1)-(2.2) admits infinitely many non-degenerate stationary traveling waves: just choose $0\ne\ell_h^-<\ell_h^+\ne1$ satisfying (4.8). Moreover, if there exists ${\rm h}\in\mathsf{H}$ such that either $D_{\rm h}(0)=0$ or $D_{\rm h}(1)=0$, then (2.1)-(2.2) admits also infinitely many degenerate stationary traveling waves: just choose $\ell_{\rm h}^-=0 = 1 - \ell_{\rm h}^+$ and determine the other end states by (4.8).

4.3.   The non-stationary case

In this subsection we consider non-stationary traveling waves. By Lemma 4.4 this is equivalent to consider the scenario in (4.4): there exists ${\rm i}\in\mathsf{I}$ such that $f_{\rm i}(\ell_{\rm i}^-) \ne f_{\rm i}(\ell_{\rm i}^+)$ and $f_j(\ell_j^-) \ne f_j(\ell_j^+)$ for every $j\in\mathsf{J}$. We can therefore introduce the following notation:

where $L_{i,j}$ is defined in (4.6) and

We notice that $\mathsf{I_0^c}\ne\emptyset$ by (4.4) and that both $\mathsf{I_0}$ and $\mathsf{I_0^c}$ depend on the end states $\ell_i^\pm$, $i\in\mathsf{I}$, indeed. Moreover, $k_j$ is well defined because by (4.7)

Finally, by (f) we deduce that

Proposition 4.9. The function $\varphi$ is the profile of a non-stationary traveling wave for (2.1)-(2.2) if and only if $\varphi_h$ is a solution to (3.5)-(3.16) for any $h \in \mathsf{H}$ and

Proof. By Proposition 4.2 it is sufficient to prove that by (4.4) condition (4.1) is equivalent to (4.11). By (3.13) we have $g_i(\ell_i^+) = g_i(\ell_i^-) = g_i(L_{i,j}^+) = g_i(L_{i,j}^-)$ and then by (4.7) we have $\kappa_j = g_j(\ell_j^\pm) - \sum_{i \in \mathsf{I}} \alpha_{i,j} g_i(L_{i,j}^\pm)$. Hence, by (4.4), with the change of variable $\xi = c_jt$, condition (4.1) is

that is equivalent to (4.11).

We observe that $k_j$ and (4.11) can be written in a little bit more explicit form by avoiding the use of $L_{i,j}^\pm$ as follows

Proposition 4.9 shows how each outgoing profile $\varphi_j$ can be expressed by (4.11) in terms of the ingoing profiles $\varphi_i$, $i\in\mathsf{I}$. We know a priori that $\varphi_j$ is increasing and its end states are contained in the interval $[0,1]$. Now, we prove a sort of converse implication, which shows that these properties of the profile $\varphi_j$ are enjoined by the function defined by the right-hand side of (4.11).

Lemma 4.10. Let $\varphi_i$, for $i\in\mathsf{I}$, be the profiles provided by Theorem 3.2 and assume that $\mathsf{I_0^c}\ne\emptyset$; fix $j\in\mathsf{J}$ and consider any $l_j^\pm\in[0,1]$ satisfying (4.7) and such that, for the corresponding $c_j$, it holds $c_j\ne0$. Then $l_j^- < l_j^+$. Moreover, denote by $\ell_j(\xi)$ the right-hand side of (4.11); then $\xi\mapsto\ell_j(\xi)$ is non-decreasing and $\ell_j(\pm\infty)=l_j^\pm$.

Proof. Since by Theorem 3.2 we know that $\ell_i^-<\ell_i^+$, then by (4.7)

By definition of $\ell_j$ we have \(\ell_j'(\xi) = \sum_{i \in \mathsf{I_0^c}} \alpha_{i, j} \, c_{i, j}^2 \, \varphi_i'\left(c_{i, j} \, \xi\right) \) for a.e. $\xi\in{\mathbb R}$, hence $\xi\mapsto\ell_j(\xi)$ is non-decreasing since all profiles $\varphi_i$ do. Moreover, $\ell_j(\pm\infty) = l_j^\pm$ because by the definitions of $\ell_j$ and $\kappa_j$ we have

We notice that Proposition 4.9 exploits condition (2.2) through its expression (4.1) for the profiles; the diffusivities $D_h$ are not involved in (4.11). Indeed, Proposition 4.9 imposes strong necessary conditions on the diffusivities as we discuss now as a preparation to (4.17).

We notice that if both $\nu_h^-$ and $\nu_h^+$ are finite, then $\ell_h^-=0=1-\ell_h^+$ and consequently $c_h=0$; therefore either $\nu_h^-$ or $\nu_h^+$ (possibly both) is infinite for any $h\in \mathsf{I_0^c} \cup \mathsf{J}$.

The following result is similar to Lemma 4.4.

Lemma 4.11. Problem (2.1)-(2.2) admits a degenerate non-stationary traveling wave $\rho$ if and only if at least one of the following conditions holds:

(A) for some $i \in \mathsf{I_0}$ we have $D_i(0) D_i(1)= 0$ and $\ell_i^- = 0$ (hence $\ell_i^+=1$);

(B) for every $h \in \mathsf{I_0^c}\cup\mathsf{J}$ we have either $D_h(0) = 0 = \ell_h^-$ or $D_h(1) = 0 = \ell_h^+-1$, but not both. In this case we have

Proof. Let us introduce the following conditions:

$(B)'$ for some $i \in \mathsf{I_0^c}$ we have either $D_i(0) = 0 = \ell_i^-$ or $D_i(1) = 0 = \ell_i^+-1$, but not both;

$(B)''$ for some $j \in \mathsf{J}$ we have either $D_j(0) = 0 = \ell_j^-$ or $D_j(1) = 0 = \ell_j^+-1$, but not both.

Clearly $(B)$ implies both $(B)'$ and $(B)''$. Moreover, by Lemma 3.5 and (4.4), problem (2.1)-(2.2) admits a degenerate non-stationary traveling wave $\rho$ if and only if at least one of the conditions (A), $(B)'$ and $(B)''$ holds. To complete the proof it is therefore sufficient to show that (B), $(B)'$ and $(B)''$ are equivalent. By Lemma 3.5(b) and (4.4), the conditions (B), $(B)'$ and $(B)''$ are respectively equivalent to

(Ⅰ) $\omega_h$ is finite for every $h \in \mathsf{I_0^c} \cup \mathsf{J}$,

(Ⅱ) for some $i\in\mathsf{I_0^c}$ we have that $\omega_i$ is finite,

(Ⅲ) for some $j\in\mathsf{J}$ we have that $\omega_j$ is finite,

where $\omega_h$ is defined in (3.14). Differentiating (4.11) gives

More precisely, by Lemma 3.5, formula (4.14) holds for $\xi\in{\mathbb R}\setminus\left(\{\omega_j\}\cup\bigcup_{i\in\mathsf{I_0^c}}\omega_i\right)$; moreover, by the same lemma we know that $\xi\mapsto\varphi_h'(c_h\xi)$ is singular at $\xi = \omega_h$ and $\mathbf{C^{1}}$ elsewhere, for $h\in\mathsf{I_0^c}\cup\mathsf{J}$. Hence, (4.14) implies (4.13). By (4.13) we have that the above statements (Ⅰ), (Ⅱ) and (Ⅲ) are equivalent and then also (B), $(B)'$ and $(B)''$ are so.

As for Lemma 4.4, we notice that Lemma 4.11 implies that a non-stationary traveling wave $\rho$ is either non-degenerate, and then $\rho_h$ is non-degenerate for every $h \in \mathsf{H}$, or $\rho$ is degenerate, and then either there exists ${\rm i} \in \mathsf{I_0}$ such that $\rho_{\rm i}$ is degenerate, or $\rho_h$ is degenerate for all $h \in \mathsf{I_0^c} \cup \mathsf{J}$. In both cases a non-stationary traveling wave $\rho$ satisfies (4.13).

When modeling traffic flows it is natural to use different diffusivities, which however share some common properties. For instance, this led to consider in [3,7] the following subcase of (D):

(D1) $D_h$ satisfies (D) and $D_h(0)=0$, $D_h(1)>0$, for every $h\in\mathsf{H}$.

The proof of the following result is an immediate consequence of Lemma 4.11 and, hence, omitted.

Corollary 1. Assume that problem (2.1)-(2.2) has a non-stationary traveling wave $\rho$ and (D1) holds. Then $\rho$ is degenerate if and only at least one of the following conditions holds:

(A) for some $i \in \mathsf{I_0}$ we have $\ell_i^- = 0$ (hence $\ell_i^+ = 1$);

(B) for every $h \in \mathsf{I_0^c}\cup\mathsf{J}$ we have $\ell_h^- = 0$ (hence $\ell_h^+ \neq 1$).

The case when $D_h$ satisfies (D) and $D_h(0)=0=D_h(1)$ for every $h\in\mathsf{H}$, see [6,7], can be dealt analogously.

The next result is the most important of this paper; there, we give necessary and sufficient conditions for the existence of non-stationary traveling waves in a network. About its statement, let us recall Theorem 3.2: we have $\varphi_h'(\xi)=0$ in case (i) if $\xi<\nu_h^-$ or in case (ii) if $\xi>\nu_h^+$. Since $\varphi_h$ satisfies equation (3.16), we are led to extend the quotient $\ell\mapsto \frac{g_h(\ell)-g_h(\ell_h^-)}{D_h(\ell)}$ to the whole of ${\mathbb R}$ by defining

In fact, when $\ell$ is replaced by $\varphi_h(\xi)$, then $\gamma_h(\ell)=\varphi_h'(\xi)$ for $\xi \in {\mathbb R} \setminus \{\nu_h^-, \nu_h^+\}$. We remark that condition $D_h(\ell)=0$ occurs at most when either $\ell=0$ or $\ell=1$. To avoid the introduction of the new notation (4.15), in the following we simply keep on writing $\frac{g_h(\ell)-g_h(\ell_h^-)}{D_h(\ell)}$ for $\gamma_h(\ell)$. As a consequence, any non-stationary traveling wave of problem (2.1)-(2.2) satisfies

Theorem 4.12. Assume conditions (f) and (D). Problem (2.1)-(2.2) admits a non-stationary traveling wave if and only if the following condition holds.

($\mathcal{T}$) There exist $\ell_1^\pm,\ldots,\ell_m^\pm\in [0,1]$ with $\ell_i^- < \ell_i^+$, $i\in\mathsf{I}$, such that:

(i) $\mathsf{I_0^c}\ne\emptyset$;

(ii) for any $j \in \mathsf{J}$ there exist $\ell_j^\pm\in[0,1]$ satisfying (4.7) and such that $f_j(\ell_j^-) \ne f_j(\ell_j^+)$;

(iii) for any $j \in \mathsf{J}$ we have

where $\varphi_1,\ldots,\varphi_m$ are solutions to (3.5)-(3.16) and, for $k_j$ as in (4.10),

Proof. First, assume that problem (2.1)-(2.2) admits a non-stationary traveling wave $\rho$ with profiles $\varphi_h$, end states $\ell_h^\pm$ and speeds $c_h$, for $h\in\mathsf{H}$. By Theorem 3.2 we have that $\ell_h^\pm$ and $c_h$ satisfy (3.7). By Proposition 4.9 the profiles $\varphi_h$ satisfy (3.5)-(3.16) and (4.11). The end states $\ell_j^\pm$, $j\in\mathsf{J}$, satisfy (4.7) by Proposition 4.6. Since $\rho$ is non-stationary we are in the scenario given by (4.4): $\mathsf{I_0^c}\ne\emptyset$ and $f_j(\ell_j^-) \ne f_j(\ell_j^+)$ for all $j \in \mathsf{J}$. By (4.16) with $h = j$ we have

for $\xi \in {\mathbb R}$ in the non-degenerate case and for $\xi \in{\mathbb R}\setminus\{ \omega\}$ with $\omega$ given by (4.13) in the degenerate case. On the other hand, by differentiating (4.11) and applying (4.16) with $h = i$ we deduce

for $\xi \in {\mathbb R}$ in the non-degenerate case and for $\xi \in{\mathbb R}\setminus\{ \omega\}$ with $\omega$ given by (4.13) in the degenerate case. Identity (4.17) follows because $\ell_j \equiv \varphi_j$ by (4.11) and by comparing (4.19), (4.20).

Conversely, assume that condition ($\mathcal{T}$) holds. We remark that the existence of $\varphi_i$, $i \in \mathsf{I}$, is assured by Theorem 3.2. Fix $j\in\mathsf{J}$. By defining $\varphi_j \doteq \ell_j$ we obtain (4.11). We know by assumption that $\mathsf{I_0^c} \ne \emptyset$, $\ell_j^\pm \in [0,1]$ satisfy (4.7) and $c_j \ne 0$; we can apply therefore Lemma 4.10 and deduce that $\ell_j^-<\ell_j^+$, $\varphi_j$ is non-decreasing and satisfies (3.5) with $h=j$. By Proposition 4.9, what remains to prove is that $\varphi_j$ satisfies (3.16). But by (4.11) we deduce (4.20) for a.e. $\xi \in {\mathbb R}$, because $\varphi_1,\ldots,\varphi_m$ satisfy (3.16) and hence, recalling the extension (4.15), also (4.16); then by (4.17) we conclude that $\varphi_j$ satisfies (4.16) for a.e. $\xi \in {\mathbb R}$ and then (3.16) for a.e. $\xi \in {\mathbb R}$. Finally, (3.16) holds by the regularity ensured by Theorem 3.2 for the profiles.

Remark 4.13. As previously in Theorem 4.7, also in Theorem 4.12 we do not use condition (2.3). We observe that (2.3) together with (4.7), or equivalently (4.2)-(4.3) by Proposition 4.6, imply

Conversely, consider the non-stationary case and assume conditions (4.2), (4.21); if moreover $\max\{f_i(\ell_i^-), f_i(\ell_i^+)\}\ne0$ and either $\sum_{j \in \mathsf{J}} \alpha_{i,j} \ge 1$ or $\sum_{j \in \mathsf{J}} \alpha_{i,j} \le 1$ for every $i\in\mathsf{I}$, then (2.3) holds. Indeed, (4.2) and (4.21) imply

namely $\sum_{i \in \mathsf{I}} \max\left\{f_i(\ell_i^-), f_i(\ell_i^+)\right \}\bigl[1- \sum_{j\in\mathsf{J}}\alpha_{i,j}\bigr]=0$, whence (2.3).

Remark 4.14. Fix $\ell_i^\pm \in [0,1]$, $i \in \mathsf{I}$, so that $\ell_i^- < \ell_i^+$ and (4.5) holds. We know by Proposition 4.5 that for every $j\in\mathsf{J}$ there exists $(\ell_j^-,\ell_j^+) \in [0,1]^2$, with $\ell_j^- < \ell_j^+$, that satisfies (4.7), but it is not unique. If beside (4.7) we impose also (4.17), then we may have three possible scenarios: such $(\ell_j^-,\ell_j^+)$ either does not exist, or it exists and is unique, or else it exists but is not unique. We refer to Subsections 6.1 and 6.2 for further discussion.

5.   The continuity condition

In this section we discuss the case when solutions to (2.1)-(2.2) are also required to satisfy the continuity condition (2.5); this makes the analysis much easier because (2.5) implies several strong conditions.

First, we provide the main results about traveling waves satisfying condition (2.5). We point out that some of the consequences below have already been pointed out in [19,30,31] in the case that some Kirchhoff conditions replace the conservation of the total flow (2.2). In order to emphasize the consequences of the continuity condition (2.5), the first two parts of the following lemma do not assume that also condition (2.2) holds.

Lemma 5.1. For any $h \in \mathsf{H}$, let $\rho_h$ be a traveling wave of $(2.1)_h$ in the sense of Definition 3.1 and set $\rho \doteq (\rho_1, \ldots, \rho_{m+n})$; then the following holds for every $(i,j) \in \mathsf{I}\times\mathsf{J}$ and $h\in\mathsf{H}$.

(i) $\rho$ satisfies (2.5) if and only if

(ii) If $\rho$ satisfies (2.5), then either it is stationary (hence (5.1) reduces to $\varphi_j(0) = \varphi_i(0)$), or it is completely non-stationary and the speeds $c_h$ have the same sign (hence $c_{i,j}>0$). In the latter case, $\rho$ is either non-degenerate or completely degenerate; moreover

(iii) If $\rho$ is non-stationary and satisfies both (2.2) and (2.5), then

Proof. We split the proof according to the items in the statement.

(i) Condition (2.5) and (5.1) are clearly equivalent.

(ii) Since we are discarding constant profiles, by (5.1) we have that either $c_h=0$ for all $h\in\mathsf{H}$ or $c_h\ne0$ for all $h\in\mathsf{H}$. The stationary case is trivial; therefore we consider below only the non-stationary case and assume that $c_h\ne0$ for all $h\in\mathsf{H}$. By differentiating (5.1) with respect to $t$ we deduce

Then (5.6) implies that either $\rho$ is non-degenerate or it is completely degenerate. Moreover (5.6) implies (5.2) because, by Lemma 3.5, we have that $\rho_h$ is degenerate if and only if the map $\xi \mapsto \varphi_h'(c_h \xi)$ is singular at $\xi = \omega_h \in {\mathbb R}$ and $\mathbf{C^{1}}$ elsewhere. By taking $t \in \mathcal{I}$ in (5.6) we deduce that $c_i$ and $c_j$ have the same sign. As a consequence we have $L_{i,j}^\pm=\ell_i^\pm$ and then $\ell_i^\pm = \ell_j^\pm$ by letting $t\to\pm\infty$ in (5.1). By (5.1), (3.16) and (5.3) we have $D_h\left(\Phi(\xi)\right) \Phi'(\xi) = c_h \left(g_h(\Phi(\xi)) - g_h(\ell^\pm)\right)$ for all $h\in \mathsf{H}$, whence (5.4) by the extension (4.15).

(iii) To deduce $(5.5)_1$, we differentiate (4.1) and then exploit (5.6). Formula $(5.5)_2$ follows by summing $(5.5)_1$ with respect to $j$ and by (2.3). By $(5.5)_1$ we have $\sum_{i \in \mathsf{I}} A_{i,j} = \sum_{i \in \mathsf{I}} \alpha_{i,j} c_i c_j^{-1}=1$, which proves $(5.5)_4$. Finally, (5.3) and $(5.5)_4$ imply $(5.5)_3$.

In the following proposition we deal with stationary traveling waves satisfying condition (2.5).

Proposition 5.2. Problem (2.1)-(2.2) admits infinitely many stationary traveling waves satisfying (2.5); their end states $\ell_h^\pm$ satisfy (4.8) and are such that$S \doteq \bigcap_{h\in\mathsf{H}}(\ell_h^-,\ell_h^+)\ne\emptyset$.

Proof. By (5.1) condition (2.5) holds in the stationary case if and only if $\varphi_i(0)=\varphi_j(0)$ for $(i,j)\in\mathsf{I}\times\mathsf{J}$. Recalling the proof of Theorem 4.7, it is sufficient to take $\ell_h^\pm \in [0,1]$ satisfying (4.8) and such that $S \ne \emptyset$, $\ell_0 \in S$ and the unique solution $\varphi_h$ to (3.5)-(3.16) such that $\varphi_h(0) = \ell_0$. There are infinitely many of such profiles because of the arbitrariness of $\ell_h^\pm$.

We point out that condition $S=\emptyset$ can occur if the functions $f_h$ assume their maximum values at different points. This is not the case when the following condition $(5.10)_1$ is assumed.

The following result is analogous to Theorem 4.12 in the case (2.5) holds.

Theorem 5.3. Assume conditions (f) and (D). Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave satisfying (2.5) if and only if the following condition holds.

($\mathcal{T}_c$) There exist $\ell^\pm\in [0,1]$ with $\ell^- < \ell^+$, such that for any $h\in\mathsf{H}$, $i \in \mathsf{I}$ and $j \in \mathsf{J}$

where $c_h$ is given by (3.7), $\varphi_h$ is a solution to (3.16) such that$\varphi_h(\pm\infty) = \ell^\pm$ and $\varphi_1(0) = \ldots = \varphi_{m+n}(0)$.

Proof. Assume that condition ($\mathcal{T}_c$) holds; the other implication is obvious. We remark that the existence of $\varphi_1,\ldots,\varphi_{m+n}$ is assured by Theorem 3.2; indeed, for any $\ell_0 \in (\ell^-,\ell^+)$, up to shifts it is always possible to assume that $\varphi_h(0) = \ell_0$, $h \in \mathsf{H}$. By (5.8) we have $(5.5)_4$ because

By (5.7) we have that $\mathsf{I_0^c} = \mathsf{I}$, $\mathsf{I_0} = \emptyset$ and $\rho$ corresponding to the profile $\varphi \doteq (\varphi_1, \ldots, \varphi_{m+n})$ is completely non-stationary. Then $(5.5)_4$ and (4.12) imply $(5.5)_3$, namely $k_j = 0$. By Lemma 5.1 (i) and Proposition 4.9 it remains to prove (5.1) and (4.11). We start with (5.1). Clearly (5.1) holds for $t=0$ because $\varphi_h(0) = \ell_0$, $h\in\mathsf{H}$. Then by the extension (4.15) and (5.9) we have

Therefore we conclude that (5.1) holds. Finally, (4.11) follows immediately from (5.1), $(5.5)_3$ and $(5.5)_4$.

Consider in particular the case when the functions $f$ and $D$ satisfy (f) and (D), respectively, and assume that

for some constants $v_h,\delta_h>0$. Denote

We notice that now we have

In the following proposition we apply Theorem 5.3 when (5.10) is assumed; in this case conditions (5.8) and (5.9) no longer depend on the end states and the statement is somewhat simplified.

Proposition 5.4. Assume (5.10) with $f$ and $D$ satisfying (f) and (D), respectively. Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave satisfying (2.5) if and only if for every $i\in\mathsf{I}$ and $j\in\mathsf{J}$ we have

Proof. We only need to translate condition $(\mathcal{T}_c)$ to the current case. Let $\ell^\pm \in [0,1]$ with $\ell^-<\ell^+$ and $f(\ell^-)\ne f(\ell^+)$. By (5.10) it is obvious that (5.7) is satisfied. If $\ell^-=0$ or $\ell^+=1$ condition (5.8) is satisfied by (f). In all the other cases (5.8) is equivalent to \(\sum_{i \in \mathsf{I}} \alpha_{i, j} v_{i, j} = 1 \) by (5.10). Similarly, condition (5.9) reduces to $c_jv_j\delta_j^{-1}=c_iv_i\delta_i^{-1}$ and hence, by (5.12), it is equivalent to $v_{i,j}^2=\delta_{i,j}$.

Remark that by (5.12) condition $(5.13)_2$ is equivalent to $(5.5)_4$.

6.   Application to the case of a quadratic flux, $m=1$

In this section we assume (5.10) for some constants $v_h,\delta_h>0$, $D$ satisfying (D) and the quadratic flux [14]

with no further mention. The case when only $(5.10)_1$ holds is doable and follows with slight modifications. We use the notation introduced in (5.11).

For simplicity, in the whole section we focus on the case $m=1$, see Figure 4, even without explicitly mentioning it. Then $\mathsf{I}=\{1\}$, $\mathsf{J}=\{2,\ldots,n+1\}$, $\mathsf{H}=\{1,2,\ldots, n+1\}$. The general case $m>1$ offers no further difficulties than heavier calculations.

In this case, condition (3.7) becomes

In particular, by $(6.1)_2$

Moreover, $g_h(\ell) = v_h \, \ell \, [\ell_h^+ + \ell_h^ - - \ell]$ implies

and therefore (3.16) becomes

We first consider stationary traveling waves and specify Theorem 4.7 and Proposition 5.2 in the current framework. We define the intervals

Proposition 6.1. Problem (2.1)-(2.2) admits infinitely many stationary traveling waves; their end states are characterized by the conditions

Moreover, up to shifts, any stationary traveling wave satisfies (2.5).

Proof. The first part of the proposition follows from Theorem 4.7. Indeed, conditions (6.2), $(3.7)_1$ and (4.8) are satisfied if and only if for any $h\in\mathsf{H}$ and $j\in\mathsf{J}$

then it is sufficient to compute $\ell_j^\pm$ and to observe that the definition of $\mathcal{L}_j^0$ guarantees that they are real numbers.

The latter part of the proposition is deduced by Proposition 5.2 because $1/2 \in s \doteq \bigcap_{h \in \mathsf{H}}(\ell_h^-, \ell_h^+) \ne \emptyset$.

In the following we treat the existence of non-stationary traveling waves. Since $m=1$, by Lemma 4.4 this is equivalent to assume $c_h\ne 0$ for $h \in \mathsf{H}$, namely, the traveling wave is completely non-stationary. By (4.7), (6.1) and (6.2), from (4.10) we deduce

The following result translates Theorem 4.12 to the present case. We define the intervals

Proposition 6.2. Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave if and only if the following condition holds.

($\mathcal{T}_q$) There exist $\ell_1^{\pm}\in [0,1]$ with $\ell_1^-<\ell_1^+$ such that:

(i) $\ell_1^++\ell_1^-\ne 1$;

(ii) $\ell_1^\pm \in \bigcap_{j \in \mathsf{J}}\mathcal{L}_j^c$;

(iii) for any $j \in \mathsf{J}$ we have

where $k_j$ is defined in (6.5) with $\ell_j^\pm$ being solutions to

Proof. The proof consists in showing that, in the present case, condition ($\mathcal{T}$) of Theorem 4.12 is equivalent to ($\mathcal{T}_q$).

$\bullet$ The first item of ($\mathcal{T}$) is clearly equivalent to the first item of ($\mathcal{T}_q$).

$\bullet$ We prove now that the second item of ($\mathcal{T}$) is equivalent to the second item of ($\mathcal{T}_q$).

''$\Rightarrow$'' Assume that for any $j \in \mathsf{J}$ there exist $\ell_j^\pm\in[0,1]$ satisfying (4.7) and such that $f_j(\ell_j^-) \ne f_j(\ell_j^+)$. Fix $j \in \mathsf{J}$. Clearly (4.7) is equivalent to (6.7).

If we denote $z_{1,j}^\pm \doteq 4\alpha_{1,j}v_{1,j}L_{1,j}^\pm(1-L_{1,j}^\pm)$, then the $\ell_j^\pm$-solutions to (6.7) are, see Figure 3,

The square roots in (6.8)-(6.9) are real numbers if and only if $z_{1,j}^\pm\le1$, namely,

It is easy to see that the above estimate is equivalent to require $\ell_1^\pm \in \mathcal{L}_j^c$.

''$\Leftarrow$'' Assume that $\ell_1^\pm \in \bigcap_{j \in \mathsf{J}}\mathcal{L}_j^c$ and fix $j \in \mathsf{J}$. The square roots in (6.8)-(6.9) are then real numbers and $\ell_j^\pm$ given in (6.8)-(6.9) satisfy (6.7), namely (4.7). Obviously $\ell_j^\pm$ belong to $[0,1]$. Finally, since $\ell_j^\pm$ are solutions to (6.7), it is easy to see that $f_j(\ell_j^+) \ne f_j(\ell_j^-)$ because $f_1(\ell_1^+) \ne f_1(\ell_1^-)$.

$\bullet$ We prove now that ($\mathcal{T}$) implies the last item of ($\mathcal{T}_q$). Since the first two items in ($\mathcal{T}$) are equivalent to the first two items in ($\mathcal{T}_q$), we can assume that $\ell_1^++\ell_1^-\ne 1$, $\ell_1^\pm \in \bigcap_{j \in \mathsf{J}}\mathcal{L}_j^c$ and that for any $j \in \mathsf{J}$ we have (4.17), namely,

for a.e. $\xi \in {\mathbb R}$, where $\varphi_1$ is a solution to (3.5)-(6.4) and

We point out that the above expression of $\ell_j$ is deduced from (4.18) by applying (6.5); moreover (6.10) is deduced from (4.17) by applying (6.3). Recall that both fractions in (6.10) are meant as in (4.15). Since

we have that (6.10) is equivalent to

To conclude now that the above condition is equivalent to (6.6) it is sufficient to recall that by Lemma 4.10 the continuous function $\xi \mapsto \ell_j(\xi)$ is increasing and $\ell_j(\pm\infty) = \ell_j^\pm$ and that $\ell_j(\xi) = A_{1,j} \, \varphi_1(c_{1,j} \, \xi) - k_j$ by (6.5).

$\bullet$ Finally, to prove that ($\mathcal{T}_q$) implies the last item of ($\mathcal{T}$) it is enough to trace backwards the proof of the previous item.

We notice that if $D$ is a polynomial with degree $\mathfrak{d}$, then (6.6) is equivalent to $\mathfrak{d} + 1$ conditions on the parameters, see for instance (6.15) and (6.25).

Remark 6.3. We point out that by Proposition 5.4 we have that problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave satisfying (2.5) if and only if

The special cases of constant or linear diffusivities are treated in the following subsections.

6.1.   The case of constant diffusivities

In this subsection we assume

and in this case problem (3.5)-(6.4) reduces to

For any $h\in\mathsf{H}$, the function

solves (6.13) because $\ell_h^- <\ell_h^+$; all the other solutions are of the form $\varphi_h(\xi)=\psi_h(\xi+\sigma_h)$ for $\sigma_h \in {\mathbb R}$. Notice that $\psi_h(0)=(\ell_h^+ + \ell_h^-)/2$.

We rewrite Proposition 6.2 in the current setting; we emphasize that the shifts appear below because in this case we have the explicit solution (6.14) to problem (6.13).

Proposition 6.4. Assume (6.12). Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave if and only if

In this case any non-stationary traveling wave $\rho$ has a profile $\varphi$ of the form

with $\ell_h^\pm$ satisfying (i), (ii) and (6.7) in Proposition 6.2 and $\sigma_h\in {\mathbb R}$, $h\in\mathsf{H}$, such that

Proof. By Theorem 3.2, any solution to (6.13) has the form (6.16) with $\sigma_h\in {\mathbb R}$, $h\in\mathsf{H}$. Therefore, by Proposition 4.2 it only remains to prove that (4.1) is equivalent to (6.15)-(6.17). Straightforward computations show that in the present case (4.1) can be written as

where $\zeta_h(t) \doteq \exp z_h(t)$, for $z_h(t) \doteq \frac{v_h}{\delta_h} \, (\ell_h^+-\ell_h^-) (c_ht+\sigma_h)$, $h \in \mathsf{H}$. By Proposition 4.6 we have

$\bullet$ In the former case, identity (6.18) is equivalent to

Since by assumption $f_j(\ell_j^+) \neq f_j(\ell_j^-)$, it must be $\zeta_j \equiv \zeta_1$, i.e., $z_j(t) = z_1(t)$, namely

$\bullet$ In the latter case, identity (6.18) is equivalent to

Since by assumption $f_j(\ell_j^+) \neq f_j(\ell_j^-)$, it must be $\zeta_j \, \zeta_1 \equiv 1$, i.e. $z_j(t) = -z_1(t)$, namely

In both cases we proved that (4.1) is equivalent to(6.15)-(6.17); this concludes the proof.

Remark 6.5. Consider conditions $(6.11)_1$, $(6.11)_2$ and (6.15). Any two of them implies the third one.

Proposition 6.6. Assume (6.12). Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave satisfying (2.5) if and only if (6.11) holds true. In this case a non-stationary traveling wave satisfies (2.5) if and only if its end states satisfy (5.3).

Proof. The first part of the statement is just Remark 6.3. In this case, since (6.11) implies (6.15), by Proposition 6.4 any (completely) non-stationary traveling wave $\rho$ has a profile of the form (6.16)-(6.17).

The second part of the statement characterizes the end states. If a non-stationary traveling wave $\rho$ satisfies (2.5), then (5.3) holds because of Lemma 5.1. Conversely, if the end states of $\rho$ satisfy (5.3), then long but straightforward computations show that (5.1) holds true, and therefore $\rho$ satisfies (2.5).

6.2.   The case of linear diffusivities

In this subsection we assume

We notice that $D$ degenerates at $0$ and this makes the subject more interesting. In this case problem (3.5)-(6.4) reduces to

If $\ell_h^-=0$, then the function

solves (6.20) because $\ell_h^- <\ell_h^+$; by (3.6) we have $I_h=\left(-\frac{\delta_h}{v_h}\ln2,\infty\right)$. If $\ell_h^->0$, then $I_h={\mathbb R}$ and the function $\psi_h$ implicitly given by

solves (6.20) because $\ell_h^- <\ell_h^+$. Notice that in both cases $\psi_h(0)=(\ell_h^+ + \ell_h^-)/2$ and all the other solutions are of the form $\varphi_h(\xi)=\psi_h(\xi+\sigma_h)$ for $\sigma_h \in {\mathbb R}$. Hence, any non-stationary traveling wave $\rho$ has a profile $\varphi$ of the form

In the sequel we prove that the shifts $\sigma_h$, $h \in \mathsf{H}$, satisfy (6.17), or equivalently

Lemma 6.7. Assume (6.19). If $\ell_1^+ + \ell_1^-\ne1$, then condition (6.6) is equivalent to

Proof. In the present case, condition (6.6) becomes $(c_{1,j}v_{1,j} - \delta_{1,j}) \, \ell - \delta_{1,j} \, k_j = 0$ for $\ell \in (\ell_j^-,\ell_j^+)$: it is satisfied if and only if both $c_{1,j}v_{1,j} = \delta_{1,j}$ and $k_j=0$. The former is equivalent to (6.25)$_1$, the latter is equivalent to (6.25)$_2$ by (6.5)$_4$, because $\kappa_j=0$.

We observe that (6.25)$_1$ and (6.5)$_1$ imply that $c_{1,j} = \delta_{1,j}/v_{1,j} > 0$; therefore (6.7) becomes

As a consequence $\rho$ is either non-degenerate or completely degenerate.

Now, we discuss (completely) non-stationary traveling waves by considering separately the (completely) degenerate and non-degenerate case. We denote

Proposition 6.8. Assume (6.19). Problem (2.1)-(2.2) admits a traveling wave that is both (completely) degenerate and (completely) non-stationary if and only if either (6.11) holds true or

In the first case, problem (2.1)-(2.2) has infinitely many of such waves; each of them satisfies (5.3) and (up to shifts) (2.5).

In the second case, problem (2.1)-(2.2) has a unique (up to shifts) such wave, which does not satisfy (for no shifts) (2.5). Its end states do not satisfy (5.3) and are

In both cases, any degenerate non-stationary traveling wave $\rho$ has a profile $\varphi$ of the form (6.23) with$\psi_h$ defined by (6.21) and $\sigma_h\in {\mathbb R}$, $h\in\mathsf{H}$, satisfying (6.24).

Proof. We claim that the existence of a degenerate non-stationary traveling wave is equivalent to the existence of $\ell_h^+\in (0,1), \, h \in \mathsf{H}$, such that

In fact, by Proposition 6.2 the existence of a non-stationary traveling wave is equivalent to condition ($\mathcal{T}_q$), where (6.6) can be written as (6.25) by Lemma 6.7 and (6.7) as (6.26). Then, (6.25) and (6.26) with $\ell_h^-=0$, $h \in \mathsf{H}$, reduce to the relation among the end states

By (6.30) we obtain $v_{1,j}^2/\delta_{1,j} = \alpha_{1,j} \, v_{1,j} \, \ell_1^+/\ell_j^+$ and then $(6.29)_1$; by plugging $(6.29)_1$ into $(6.30)_2$ we get $(6.29)_2$ and then the claim.

Assume there is a degenerate non-stationary traveling wave; then $\ell_1^+$ satisfies (6.29)$_2$. As a consequence, we have either $\alpha_{1,j}\delta^2_{1,j}- v^3_{1,j}=v_{1,j}^2-\delta_{1,j}=0$ or $\alpha_{1,j}\delta^2_{1,j}\ne v^3_{1,j}$ for every $j \in \mathsf{J}$. The former case is equivalent to (6.11). In the latter case we can explicitly compute $\ell_1^+$ by $(6.29)_2$ for any $j \in \mathsf{J}$ and impose the constraint $0<\ell_1^+<1$, namely,

A direct computation shows that this is equivalent to (6.27). In conclusion, either condition (6.11) or (6.27) is necessary for the existence of a non-stationary traveling wave with $\ell_1^-=0$.

Conversely, assume condition (6.11). In this case $\alpha_{1,j} \, \delta_{1,j}^2 =\alpha_{1,j} \, v_{1,j}^4=v_{1,j}^3$. Then (6.29)$_2$ is trivially satisfied for every $\ell_1^+\in (0,1)$ and from (6.29)$_1$ we deduce $\ell_1^+=\ell_j^+$. Hence, there is an infinite family of non-stationary traveling waves parameterized by $\ell_1^+ \in (0,1)$ and satisfying (5.3); as a consequence, they do not coincide up to shifts and they all satisfy (up to shifts) the continuity condition (2.5).

Assume now condition (6.27). In this case the values for $\ell_1^+$ and $\ell_j^+$ in (6.28) are well defined since $v^3_{1,j}\ne \alpha_{1,j}\delta_{1,j}$ and they are the unique solution to (6.29). In particular, condition (ii) in Proposition 6.2 is automatically satisfied. By the estimates in (6.27) we have $\ell_1^+,\ell_j^+ \in (0,1)$ for $j \in \mathsf{J}$. Hence, there is a unique (up to shifts) degenerate non-stationary traveling wave and its end states satisfy (6.28). Furthermore, by Lemma 5.1, condition (2.5) implies (5.3), which is precluded by (6.27). Hence, the traveling wave does not satisfy (2.5).

At last, by Theorem 3.2, any solution to (6.20) has the form (6.23). By (4.14), that in the present case becomes

and the regularity of $\psi_h$ defined in (6.21), we have

which is equivalent to (6.17) because $c_{1,j} = \delta_{1,j}/v_{1,j}$.

The following result treats the non-degenerate case.

Proposition 6.9. Assume (6.19). Problem (2.1)-(2.2) admits a non-degenerate (completely) non-stationary traveling wave if and only if condition (6.11) is satisfied. In this case any non-degenerate non-stationary traveling wave satisfies (up to shifts) (2.5); moreover, it has a profile $\varphi$ of the form (6.23) with $\psi_h$ implicitly defined by (6.22) and$\sigma_h\in {\mathbb R}$, $h\in\mathsf{H}$, satisfying (6.24).

Proof. Assume that there is a non-degenerate non-stationary traveling wave; then $\ell_h^- \ne 0$ and $1 \ne \ell_h^+ + \ell_h^-$, $h \in \mathsf{H}$. Moreover, by Proposition~6.2, condition ($\mathcal{T}_q$) is satisfied, where (6.6) becomes (6.25) by Lemma 6.7 and (6.7) is (6.26). When dividing (6.26) by (6.25)$_2$ we obtain

By adding the above relations we have $\ell_1^- + \ell_1^+ = \ell_j^- + \ell_j^+$, hence

It is now easy to conclude that (5.3) is satisfied and then also (6.11) holds true by (6.25). At last, the traveling wave satisfies (up to shifts) (2.5) by Remark 6.3.

Conversely, assume (6.11). Then (6.25) and (6.26) write

The same computations as before give that if we impose $\ell_h^- \ne 0$ and $1 \ne \ell_h^+ + \ell_h^-$, $h \in \mathsf{H}$, then the above conditions are equivalent to (5.3); the existence of infinitely many non-degenerate non-stationary traveling waves satisfying (2.5) easily follows.

At last, by Theorem 3.2, any solution to (6.20) has the form (6.23). Fix $j \in \mathsf{J}$. By (2.5) we have (5.1), namely

This identity together with (6.22) and (5.3) imply

By dividing the above equalities and taking the logarithm we get

Since $\ell^- \ne \ell^+$ and $c_{1,j} = \delta_{1,j}/v_{1,j}$, the above equality is equivalently to (6.24).

7.   Application to the case of a logarithmic flux, $m=1$

In this section we assume (5.10) for some constants $v_h,\delta_h>0$, $D \doteq 1$ and the logarithmic flux [13] defined by

for $\rho\in(0,1]$ with $f(0)=0$ by continuity; in the following we simply write $\rho \ln (\rho)$ for $\rho \in [0,1]$. We use the notation introduced in (5.11); then, in the present case the diffusivity $D_h$ coincides with the anticipation length $\delta_h$ of [3], see Section 2. As in Section 6, we focus on the case $m=1$ and do not mention in the following these assumptions on $f_h$, $D_h$ and $m$.

Condition (3.7) becomes

Moreover we have, for $h \in \mathsf{H}$,

Therefore (3.16) becomes

for $\xi \in {\mathbb R}$. Let $f_\ell^{-1}:[0,e^{-1}] \to [0,e^{-1}]$ and $f_r^{-1}:[0,e^{-1}] \to [e^{-1},1]$ be the inverse functions of the restrictions $f_\ell$ and $f_r$ of $f$ to $[0, e^{-1}]$ and $[e^{-1}, 1]$, respectively.

We first consider the case of stationary waves. We define the intervals

Proposition 7.1. Problem (2.1)-(2.2) admits infinitely many stationary traveling waves; their end states are characterized by the conditions

Moreover, up to shifts, any stationary traveling wave satisfies (2.5).

Proof. The first part of the proposition follows from Theorem 4.7. Indeed, conditions (3.7)$_1$ and (4.8) are satisfied if and only if for any $h\in\mathsf{H}$ and $j\in\mathsf{J}$

Hence $\ell_1^+=f_r^{-1}(-\ell_1^-\ln(\ell_1^-))$ and it is sufficient to determine $\ell_j^\pm$. Observe that the definition of $\mathcal{L}_j^0$ guarantees that they can be uniquely computed. At last, the latter part of the proposition follows by the proof of Proposition 5.2 since $e^{-1} \in S \doteq \bigcap_{h\in\mathsf{H}}(\ell_h^-,\ell_h^+) \ne \emptyset$.

In the following we discuss the existence of non-stationary traveling waves. Since $m=1$, by Lemma 4.4 this is equivalent to assume that the traveling wave is completely non-stationary. By (7.1)$_2$ we deduce

The following result translates Theorem 4.12 to the current framework. We define the intervals

Proposition 7.2. Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave if and only if the following condition holds.

($\mathcal{T}_l$) There exist $\ell_1^{\pm}\in [0,1]$ with $\ell_1^-<\ell_1^+$ such that:

(i) $\ell_1^-\ln(\ell_1^-)\ne \ell_1^+ \ln(\ell_1^+)$;

(ii) $\ell_1^\pm \in \bigcap_{j \in \mathsf{J}}\mathcal{L}_j^c$;

(iii) for any $j \in \mathsf{J}$ we have

where $g_h$ is given in (7.2), $c_{1,j}$ in (7.4), $A_{1,j}$ in (4.10)$_2$ and $k_j$ in (4.10)$_3$, with $\ell_j^\pm$ being solutions to

Proof. The proof consists in showing that, in the present case, ($\mathcal{T}$) of Theorem 4.12 is equivalent to ($\mathcal{T}_l$). The first two items in ($\mathcal{T}$) and ($\mathcal{T}_l$) are clearly equivalent. It remains to discuss the third one. Condition (4.17) is equivalent to

where $\varphi_1$ is a solution to (3.5)-(7.3) and $\ell_j(\xi) \doteq A_{1,j}\varphi_1(c_{1,j}\xi)-k_j$ for $c_{1,j}$ in (7.4), $A_{1,j}$ in (4.10)$_2$ and $k_j$ in (4.10)$_3$. By Theorem 3.2, $\varphi_1$ is strictly increasing and so is the function $\ell_j$. Put $\ell \doteq \ell_j(c_j\xi)$. Hence $\ell \in (\ell_j^-, \, \ell_j^+)$, by Lemma 4.10, and then (7.7) is equivalent to (7.5).

In the following we focus on the case of (completely) non-stationary traveling waves with $\ell_h^-=0$ for some $h \in \mathsf{H}$.

Lemma 7.3. Assume that problem (2.1)-(2.2) admits a traveling wave. The following statements are equivalent:

(i) $\ell_1^- = 0$;

(ii) $\ell_j^- = 0$ for all $j \in \mathsf{J}$;

(iii) there exists ${\rm j} \in \mathsf{J}$ such that $\ell_{\rm j}^- = 0$.

Proof. First, we prove that (i) implies (ii). Fix $j \in \mathsf{J}$. Since $\ell_1^-=0$, then condition (7.6) implies that either $\ell_j^-=0$ or $\ell_j^+=1$, for $j \in \mathsf{J}$. Assume by contradiction that $\ell_j^+=1$. Since $c_{1,j}<0$, condition (7.6) becomes

Therefore, by (7.4), (4.10)$_2$ and (4.10)$_3$ we have that

Condition (7.5) can be written as

for $\ell \in (\ell_j^-, 1)$. By differentiating the above equation three times we obtain

This is a contradiction because the two sides have opposite sign. This proves (ii).

Since the implication (ii) $\Rightarrow$ (iii) is obvious, it remains to show that (iii) $\Rightarrow$ (i). Let $\ell_{\rm j}^-=0$ for some ${\rm j} \in \mathsf{J}$. By (7.6) it follows that either $\ell_1^-=0$ or $\ell_1^+=1$. In the latter case by arguing as above it is easy to obtain a contradiction and then (iii) follows.

At last, we give a result which is similar to the one given in Proposition 6.8. We denote

By Lemma 7.3 we have either $\ell_h^-=0$, $h \in \mathsf{H}$, or $\ell_h^-\ne0$, $h \in \mathsf{H}$. Below we consider the first case.

Proposition 7.4. Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave with $\ell_h^-=0$, $h \in \mathsf{H}$, if and only if either (6.11) holds true or

In the first case, problem (2.1)-(2.2) has infinitely many of such waves; each of them satisfies (5.3) and (up to shifts) (2.5).

In the second case, problem (2.1)-(2.2) has a unique (up to shifts) such wave and such wave, which does not satisfy (for no shifts) (2.5). Its end states are

and do not satisfy (5.3).

Proof. Fix $j \in \mathsf{J}$. Since $c_{1,j} > 0$, by (7.4) the formulas in (4.10) and (7.2) become

Hence (7.6) can be written as

and therefore (7.5) becomes

namely

System (7.10)-(7.11) admits a solution if and only if either (6.11) or (7.8) holds true. In the former case, (7.10)-(7.11) has infinitely many solutions and they satisfy (5.3); in the latter, the unique solution of (7.10)-(7.11) is (7.9)$_{2,3}$. We examine separately these cases.

Assume (6.11). In this case condition ($\mathcal{T}_l$) of Proposition 7.2 with $\ell_1^-=0=\ell_j^-$ is equivalent to $\ell_1^+=\ell_j^+ \in (0,1), \, j \in \mathsf{J}$, and then there are infinitely many traveling waves. They all satisfy (2.5) by Remark 6.3.

Assume (7.8). In this case condition ($\mathcal{T}_l$) of Proposition 7.2 with $\ell_1^-=0=\ell_j^-$ is equivalent to $\ell_h^+ \in (0,1)$, $h \in \mathsf{H}$, satisfying (7.10)-(7.11), namely to (7.8)-(7.9). In particular, (7.8)$_1$, (7.9) imply that $\ell_j^+$ and $\ell_1^+$ are distinct, namely they do not satisfy (5.3). Moreover, by Remark 6.3 the traveling wave does not satisfies (2.5).

At last, the reverse implications are direct consequences of previous discussion about the solutions of (7.10)-(7.11) and then the proof is complete.

Appendix A.   Proof of Theorem 3.2

Let $\ell^\pm_h \in [0, 1]$ with $\ell^{-}_h\ne \ell^{+}_h$. We introduce the change of variable

which implies $\rho_h = \ell^{+}_h - (\ell^{+}_h-\ell^{-}_h) \, r_h$, $\rho_{h,t} = -(\ell^{+}_h - \ell^{-}_h) \, r_{h,t}$ and $\rho_{h,x} = -(\ell^{+}_h - \ell^{-}_h) \, r_{h,x}$. Consequently, equation (2.1) can be written

where

Furthermore, equation (A.2) has a wavefront solution $\psi_h$ from $1$ to $0$ with wave speed $\theta_h$ if and only if equation (2.1) has a wavefront solution $\varphi_h$ from $\ell^{-}_h$ to $\ell^{+}_h$ with the same speed. Notice that $\psi_h$ satisfies the equation

and $\varphi_h$ is obtained by $\psi_h$ by the change of variable (A.1), i.e.

We discuss now the existence of a wavefront solution $r_h(t,x)=\psi_h(x-\theta_h t+\sigma_h)=\psi_h(\xi)$ of (A.2). In order to make use of [12,Theorem 9.1], we only need to show that

By the definition of $G_h$ we have

Then, inequality (A.4) is equivalent to

if and only if $\ell_h^-<\ell_h^+$. By the strict concavity of $f_h$ the last inequality is satisfied and then, by [12,Theorem 9.1], we deduce the existence of wavefront solutions $\psi_h$ from $1$ to $0$ for (A.2). The wave speed, in this case, is $\theta_h \doteq G_h(1)$. Furthermore, the profile $\psi_h$ is unique up to shifts and, if $\psi_h(0)\doteq\nu$ for some $0<\nu< 1$, then

where

Notice that, by differentiating (A.5) in the interval $(\nu_h^-, \nu_h^+)$, we obtain that

which implies $\psi_h'<0$ in $(\nu_h^-, \nu_h^+)$ because of (A.4).

Consider now $\varphi_h$ defined in (A.3); it satisfies (3.3) with $I_h=(\nu_h^-, \nu_h^+)$ and $\varphi_h'>0$ in $I_h$. Also condition (3.7) is true and $\varphi_h \in \mathbf{C^{2}}(I_h, (\ell^{-}_h, \ell^{+}_h)$ by the regularity of $D_h$ and $f_h$.

Now it remains to consider the boundary conditions of $\varphi_h'$ at the extrema of $I_h$ in the different cases. We have the following.

(i) Assume $\ell_h^-=0=D_h(0)$. We show that

To prove (A.7), notice that $E_h(1)=D_h(0)=0$ and that $-G_{h}(s)+ sG_{h}(1) \to 0$ as $s\to 1^-$. In addition, by means of the strict concavity of $f_h$ we obtain that

and then, by applying de l'Hospital Theorem we prove condition (A.7). Moreover, by condition (A.6), we get

By applying (A.3) we conclude that $\varphi_h(\xi)=\ell_h^-$ for $\xi \le \nu_h^-$ and the estimates in (3.8) are satisfied. Furthermore, by the change of variables (A.1), we obtain that

and hence, by (A.6), we deduce (3.9).

(ii) Assume $1-\ell_h^+=0=D_h(1)$. With a similar reasoning as in (i) we prove that $\nu_h^+>-\infty$. In fact, $E_h(0)=D_h(1)=0$ and $-s \, G_h(s)+s \, G_h(1) \to 0$ as $s \to 0^+$. Moreover

and again, by applying de l'Hospital Theorem we prove that $\nu_h^+>-\infty$. Moreover, by the estimate (A.6), we have that

By applying (A.3) we conclude that $\varphi_h(\xi)=\ell_h^+$ for $\xi \ge \nu_h^+$ and the estimates in (3.10) are satisfied; by (A.1) and (A.6) we derive (3.11).

(iii) In all the other cases it is easy to show that $I_h=\mathbb{R}$ and again the slope condition (3.12) can be obtained by the estimate (A.6).

Acknowledgments

The authors are members of INdAM-GNAMPA. They were supported by the Project Macroscopic models of traffic flows: qualitative analysis and implementation, sponsored by the University of Modena and Reggio Emilia. The first author was also supported by the INdAM -GNAMPA Project 2016 ''Balance Laws in the Modeling of Physical, Biological and Industrial Processes''. The last author was also supported by the INdAM -GNAMPA Project 2017 ''Nonlocal hyperbolic equations''.