Traveling waves for degenerate diffusive equations on networks

TL;DR

This paper investigates the existence of traveling wave solutions for degenerate diffusive equations on star graphs, with applications to traffic flow modeling at crossroads, providing algebraic conditions and detailed examples.

Contribution

It establishes a necessary and sufficient algebraic condition for traveling wave existence on networks with degenerate diffusion, extending previous results to more general settings.

Findings

Derived algebraic condition for traveling wave existence

Analyzed quadratic and logarithmic flux functions

Applied results to traffic flow models

Abstract

In this paper we consider a scalar parabolic equation on a star graph; the model is quite general but what we have in mind is the description of traffic flows at a crossroad. In particular, we do not necessarily require the continuity of the unknown function at the node of the graph and, moreover, the diffusivity can be degenerate. Our main result concerns a necessary and sufficient algebraic condition for the existence of traveling waves in the graph. We also study in great detail some examples corresponding to quadratic and logarithmic flux functions, for different diffusivities, to which our results apply.

Equations

$$\displaystyle\rho_{h,t}+f_{h}(\rho_{h})_{x}=\mathopen{}\mathclose{{}\left(D_{h}(\rho_{h})\rho_{h,x}}\right)_{x},$$

$$F_h({\rho }_h,{\rho }_{h,x})\doteq f_h({\rho }_h)-D_h({\rho }_h){\rho }_{h,x}$$

$$\displaystyle F_{j}\mathopen{}\mathclose{{}\left(\rho_{j}(t,0^{+}),\rho_{j,x}(t,0^{+})}\right)=\sum_{i\in\mathsf{I}}\alpha_{i,j}\,F_{i}\mathopen{}\mathclose{{}\left(\rho_{i}(t,0^{-}),\rho_{i,x}(t,0^{-})}\right)$$

$$\sum _{j\in \mathsf{J}}{\alpha }_{i,j}=1,$$

$$\displaystyle\sum_{j\in\mathsf{J}}F_{j}\mathopen{}\mathclose{{}\left(\rho_{j}(t,0^{+}),\rho_{j,x}(t,0^{+})}\right)=\sum_{i\in\mathsf{I}}F_{i}\mathopen{}\mathclose{{}\left(\rho_{i}(t,0^{-}),\rho_{i,x}(t,0^{-})}\right)$$

$${\rho }_i(t,0^{-})={\rho }_j(t,0^{+}),$$

$$\displaystyle\rho_{h,t}+f_{h}(\rho_{h})_{x}=\mathopen{}\mathclose{{}\left(D_{h}(\rho_{h})\rho_{h,x}}\right)_{x},$$

$$\mathopen{}\mathclose{{}\left[F_{h}(\varphi_{h},\varphi_{h}^{\prime})-c_{h}\varphi_{h}}\right]^{\prime}=0,$$

$$\displaystyle\mathopen{}\mathclose{{}\left(D_{h}(\varphi_{h})\varphi_{h}^{\prime}}\right)^{\prime}-g_{h}^{\prime}(\varphi_{h})\,\varphi_{h}^{\prime}=0,$$

$$g_h(\rho )\doteq f_h(\rho )-c_h\rho$$

$$\int_{\mathbb{R}}\mathopen{}\mathclose{{}\left[D_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi)}\right)\varphi_{h}^{\prime}(\xi)-g_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi)}\right)}\right]\psi^{\prime}(\xi)\,{\rm d}\xi=0,$$

$${\phi }_h(\pm \infty )={\ell }_h^{\pm },$$

$$I_{h}\doteq\mathopen{}\mathclose{{}\left\{\xi\in\mathbb{R}\,:\,\ell_{h}^{-}<\varphi_{h}(\xi)<\ell_{h}^{+}}\right\}.$$

$$0\le {\ell }_h^{-}<{\ell }_h^{+}\le 1$$

$$\underset{\xi \downarrow {\nu }_h^{-}}{\lim }{\phi }_h^{\prime }(\xi )=\{\begin{array}{ll}\frac{{\ell }_h^{+}{f}_h^{\prime }(0)-f_h({\ell }_h^{+})}{{\ell }_h^{+}{D}_h^{\prime }(0)}& if{D}_h^{\prime }(0)>0,\\ \infty & if{D}_h^{\prime }(0)=0,\end{array}$$

$$\displaystyle\lim_{\xi\downarrow\nu_{h}^{-}}D_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi)}\right)\varphi_{h}^{\prime}(\xi)=0.$$

$$\displaystyle\lim_{\xi\uparrow\nu_{h}^{+}}\varphi_{h}^{\prime}(\xi)=\begin{cases}\frac{\mathopen{}\mathclose{{}\left(1-\ell_{h}^{-}}\right)f_{h}^{\prime}(1)+f_{h}(\ell_{h}^{-})}{\mathopen{}\mathclose{{}\left(1-\ell_{h}^{-}}\right)D_{h}^{\prime}(1)}&\hbox{if }D_{h}^{\prime}(1)<0,\\ \infty&\hbox{if }D_{h}^{\prime}(1)=0,\end{cases}$$

$$\displaystyle\lim_{\xi\uparrow\nu_{h}^{+}}D_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi)}\right)\varphi_{h}^{\prime}(\xi)=0.$$

$$\underset{\xi \to \pm \infty }{\lim }{\phi }_h^{\prime }(\xi )=0.$$

$$g_h({\ell }_h^{+})=g_h({\ell }_h^{-})=-\frac{f_h({\ell }_h^{+}){\ell }_h^{-}-f_h({\ell }_h^{-}){\ell }_h^{+}}{{\ell }_h^{+}-{\ell }_h^{-}}$$

$$I_h=({\nu }_h^{-},{\nu }_h^{+}).$$

$${\omega }_h\doteq \min \{c_h^{-1}{\nu }_h^{-},c_h^{-1}{\nu }_h^{+}\}.$$

$$c_{h}\varphi_{h}(\xi)-F_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi),\varphi_{h}^{\prime}(\xi)}\right)=c_{h}\varphi_{h}(\xi_{-})-F_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi_{-}),\varphi_{h}^{\prime}(\xi_{-})}\right).$$

$$\displaystyle F_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi),\varphi_{h}^{\prime}(\xi)}\right)=c_{h}\varphi_{h}(\xi)+g_{h}(\ell_{h}^{\pm}),$$

$$\displaystyle D_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi)}\right)\,\varphi_{h}^{\prime}(\xi)=g_{h}\mathopen{}\mathclose{{}\left(\varphi_{h}(\xi)}\right)-g_{h}(\ell_{h}^{\pm}),$$

$$\displaystyle c_{j}\,\varphi_{j}(c_{j}t)+g_{j}(\ell_{j}^{\pm})=\sum_{i\in\mathsf{I}}\alpha_{i,j}\mathopen{}\mathclose{{}\left[c_{i}\,\varphi_{i}(c_{i}t)+g_{i}(\ell_{i}^{\pm})}\right],$$

$$\displaystyle F_{j}\bigl{(}\varphi_{j}(-c_{j}t),\varphi_{j}^{\prime}(-c_{j}t)\bigr{)}=\sum_{i\in\mathsf{I}}\alpha_{i,j}\,F_{i}\bigl{(}\varphi_{i}(-c_{i}t),\varphi_{i}^{\prime}(-c_{i}t)\bigr{)},$$

$$\displaystyle\max\mathopen{}\mathclose{{}\left\{f_{j}(\ell_{j}^{-}),f_{j}(\ell_{j}^{+})}\right\}$$

$$\displaystyle\min\mathopen{}\mathclose{{}\left\{f_{j}(\ell_{j}^{-}),f_{j}(\ell_{j}^{+})}\right\}$$

$${{\rm Y}}_j(t)=\sum _{i\in \mathsf{I}}{\alpha }_{i,j}{{\rm Y}}_i(t),$$

$$\displaystyle\lim_{t\to\infty}\Upsilon_{h}(t)=\max\mathopen{}\mathclose{{}\left\{f_{h}(\ell_{h}^{-}),f_{h}(\ell_{h}^{+})}\right\},$$

Full text

AMS Subject Classification: 35K65; 35C07, 35K55, 35K57

Keywords: Parabolic equations, wavefront solutions, traffic flow on networks.

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, 23, 29]. 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

[TABLE]

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})\geq 0$ 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, 24]. We refer to [3, 15, 17, 19, 21, 26] 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, 25] 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 [20]; in the case of a second-order model without diffusion but including a relaxation term, we refer to [9, 27]; 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 [29, 30] for the semilinear diffusive case and to [18] 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, 25]; 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.

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

[TABLE]

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, 24]. 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 [22]. The simplest model for the diffusivity is then $D_{h}(\rho_{h})=-\delta_{h}\rho_{h}\,v_{h}^{\prime}(\rho_{h})$, where $\delta_{h}$ is an anticipation length [3, 20].

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

[TABLE]

and require

[TABLE]

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

[TABLE]

Conditions (2.2) and (2.3) imply

[TABLE]

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

[TABLE]

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.

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

[TABLE]

where we keep for future reference the index $h$. Equation (3.1) has no source terms (differently from [29, 30]) 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 [18, 28] because we are considering non-constant profiles. The profile must satisfy the equation

[TABLE]

namely,

[TABLE]

in the weak sense, where

[TABLE]

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}^{\prime}\in\mathbf{L^{1}_{loc}}(\mathbb{R};\mathbb{R})$ and

[TABLE]

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

[TABLE]

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

[TABLE]

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

[TABLE]

We have that \varphi_{h}\in\mathbf{C}^{2}\mathopen{}\mathclose{{}\left(I_{h};(\ell_{h}^{-},\ell_{h}^{+})}\right) is unique (up to shifts) and $\varphi_{h}^{\prime}(\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\leq\nu_{h}^{-}$. In this case*

[TABLE]

• (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\geq\nu_{h}^{+}$. In this case*

[TABLE]

• (iii)

In all the other cases $I_{h}=\mathbb{R}$ and

[TABLE]

We observe that for $c_{h}$ given by (3.7), we deduce by (f) that $g_{h}(\rho)\geq 0$ for all $\rho\in[\ell_{h}^{-},\ell^{+}_{h}]$, see Figure 2. Moreover, we have

[TABLE]

and no $\rho\neq\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 (i) or (ii) 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 (i) or at $\nu_{h}^{+}$ in case (ii), in the sense that $\varphi_{h}^{\prime}$ 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):

[TABLE]

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

[TABLE]

Lemma 3.5**.**

Let $\rho_{h}$ be a traveling-wave solution of (3.1); then we have the following.

1. (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$). 2. (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}^{\prime}(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 [math] or $1$, say [math]; 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}\neq 0$, 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}^{\prime}(\xi)$ is singular at $\xi=\nu_{h}^{-}$ and $\mathbf{C^{1}}$ elsewhere, whence the regularity of $\xi\mapsto\varphi_{h}^{\prime}(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}^{\prime}(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

[TABLE]

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

[TABLE]

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

[TABLE]

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 $\eqref{eq:model}_{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 if and only if $\varphi_{h}$ is a solution to (3.5)-(3.16) for any $h\in\mathsf{H}$ and

[TABLE]

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

Proof.

By plugging $\rho_{h}(t,x)=\varphi_{h}(x-c_{h}t)$ in (2.2) and recalling that by Theorem 3.2 the profiles are continuous in $\mathbb{R}$, we obtain

[TABLE]

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

[TABLE]

Proof.

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

[TABLE]

where the map $t\mapsto\Upsilon_{h}(t)\doteq c_{h}\,\varphi_{h}(c_{h}t)+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_{h}t)-\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,

[TABLE]

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:

1. (i)

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

$c_{i}=0$* for all $i\in\mathsf{I}$;* 3. (iii)

$c_{j}=0$* for all $j\in\mathsf{J}$.*

Proof.

By subtracting (4.3) to (4.2) we obtain

[TABLE]

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

[TABLE]

Of course, by Lemma 4.4, $c_{\rm i}\neq 0$ for some ${\rm i}\in\mathsf{I}$ if and only if $c_{j}\neq 0$ 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

[TABLE]

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}\geq\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}\neq 0$ for some ${\rm i}\in\mathsf{I}$; then $c_{j}\neq 0$ by Lemma 4.4, i.e., $f_{j}(\ell_{j}^{-})\neq 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

[TABLE]

Proposition 4.6**.**

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

[TABLE]

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

[TABLE]

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

[TABLE]

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). ∎

Clearly, if both $D_{h}(0)\neq 0$ and $D_{h}(1)\neq 0$ 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\neq\ell_{h}^{-}<\ell_{h}^{+}\neq 1$ 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}^{-})\neq f_{\rm i}(\ell_{\rm i}^{+})$ and $f_{j}(\ell_{j}^{-})\neq f_{j}(\ell_{j}^{+})$ for every $j\in\mathsf{J}$. We can therefore introduce the following notation:

[TABLE]

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

[TABLE]

We notice that $\mathsf{I_{0}^{c}}\neq\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)

[TABLE]

Finally, by (f) we deduce that

[TABLE]

Proposition 4.8**.**

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

[TABLE]

Proof.

By Proposition 4.2 it is sufficient to prove that by (4.4) condition (4.1) is equivalent to (4.10). 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_{j}t$, condition (4.1) is

[TABLE]

that is equivalent to (4.10). ∎

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

[TABLE]

Proposition 4.8 shows how each outgoing profile $\varphi_{j}$ can be expressed by (4.10) 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.10).

Lemma 4.9**.**

Let $\varphi_{i}$, for $i\in\mathsf{I}$, be the profiles provided by Theorem 3.2 and assume that $\mathsf{I_{0}^{c}}\neq\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}\neq 0$. Then $l_{j}^{-}<l_{j}^{+}$. Moreover, denote by $\ell_{j}(\xi)$ the right-hand side of (4.10); 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)

[TABLE]

By definition of $\ell_{j}$ we have \ell_{j}^{\prime}(\xi)=\sum_{i\in\mathsf{I_{0}^{c}}}\alpha_{i,j}\,c_{i,j}^{2}\,\varphi_{i}^{\prime}\mathopen{}\mathclose{{}\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

[TABLE]

We notice that Proposition 4.8 exploits condition (2.2) through its expression (4.1) for the profiles; the diffusivities $D_{h}$ are not involved in (4.10). Indeed, Proposition 4.8 imposes strong necessary conditions on the diffusivities as we discuss now as a preparation to (4.16).

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.10**.**

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

[TABLE]

Proof.

Let us introduce the following conditions:

1. $(B)^{\prime}$

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; 2. $(B)^{\prime\prime}$

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)^{\prime}$ and $(B)^{\prime\prime}$. 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)^{\prime}$ and $(B)^{\prime\prime}$ holds. To complete the proof it is therefore sufficient to show that (B), $(B)^{\prime}$ and $(B)^{\prime\prime}$ are equivalent. By Lemma 3.5(b) and (4.4), the conditions (B), $(B)^{\prime}$ and $(B)^{\prime\prime}$ are respectively equivalent to

1. (I)

$\omega_{h}$ is finite for every $h\in\mathsf{I_{0}^{c}}\cup\mathsf{J}$, 2. (II)

for some $i\in\mathsf{I_{0}^{c}}$ we have that $\omega_{i}$ is finite, 3. (III)

for some $j\in\mathsf{J}$ we have that $\omega_{j}$ is finite,

where $\omega_{h}$ is defined in (3.14). Differentiating (4.10) gives

[TABLE]

More precisely, by Lemma 3.5, formula (4.13) holds for \xi\in\mathbb{R}\setminus\mathopen{}\mathclose{{}\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}^{\prime}(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.13) implies (4.12). By (4.12) we have that the above statements (I), (II) and (III) are equivalent and then also (B), $(B)^{\prime}$ and $(B)^{\prime\prime}$ are so. ∎

As for Lemma 4.4, we notice that Lemma 4.10 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.12).

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.10 and, hence, omitted.

Corollary 4.11**.**

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}^{\prime}(\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

[TABLE]

In fact, when $\ell$ is replaced by $\varphi_{h}(\xi)$, then $\gamma_{h}(\ell)=\varphi_{h}^{\prime}(\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.14), 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

[TABLE]

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:

1. (i)

$\mathsf{I_{0}^{c}}\neq\emptyset$; 2. (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}^{-})\neq f_{j}(\ell_{j}^{+})$; 3. (iii)

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

[TABLE]

where $\varphi_{1},\ldots,\varphi_{m}$ are solutions to (3.5)-(3.16) and, for $k_{j}$ as in (4.9),

[TABLE]

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.8 the profiles $\varphi_{h}$ satisfy (3.5)-(3.16) and (4.10). 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}}\neq\emptyset$ and $f_{j}(\ell_{j}^{-})\neq f_{j}(\ell_{j}^{+})$ for all $j\in\mathsf{J}$. By (4.15) with $h=j$ we have

[TABLE]

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

[TABLE]

for $\xi\in\mathbb{R}$ in the non-degenerate case and for $\xi\in\mathbb{R}\setminus\{\omega\}$ with $\omega$ given by (4.12) in the degenerate case. Identity (4.16) follows because $\ell_{j}\equiv\varphi_{j}$ by (4.10) and by comparing (4.18), (4.19).

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.10). We know by assumption that $\mathsf{I_{0}^{c}}\neq\emptyset$, $\ell_{j}^{\pm}\in[0,1]$ satisfy (4.7) and $c_{j}\neq 0$; we can apply therefore Lemma 4.9 and deduce that $\ell_{j}^{-}<\ell_{j}^{+}$, $\varphi_{j}$ is non-decreasing and satisfies (3.5) with $h=j$. By Proposition 4.8, what remains to prove is that $\varphi_{j}$ satisfies (3.16). But by (4.10) we deduce (4.19) for a.e. $\xi\in\mathbb{R}$, because $\varphi_{1},\ldots,\varphi_{m}$ satisfy (3.16) and hence, recalling the extension (4.14), also (4.15); then by (4.16) we conclude that $\varphi_{j}$ satisfies (4.15) 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**.**

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.16), 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 [18, 29, 30] 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 $\eqref{eq:model}_{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}$.

1. (i)

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

[TABLE] 2. (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

[TABLE] 3. (iii)

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

[TABLE]

Proof.

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

1. (i)

Condition (2.5) and (5.1) are clearly equivalent. 2. (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}\neq 0$ 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}\neq 0$ for all $h\in\mathsf{H}$. By differentiating (5.1) with respect to $t$ we deduce

[TABLE]

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}^{\prime}(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}\mathopen{}\mathclose{{}\left(\Phi(\xi)}\right)\Phi^{\prime}(\xi)=c_{h}\mathopen{}\mathclose{{}\left[g_{h}(\Phi(\xi))-g_{h}(\ell^{\pm})}\right] for all $h\in\mathsf{H}$, whence (5.4) by the extension (4.14). 3. (iii)

To deduce $\eqref{e:4eqn}_{1}$, we differentiate (4.1) and then exploit (5.6). Formula $\eqref{e:4eqn}_{2}$ follows by summing $\eqref{e:4eqn}_{1}$ with respect to $j$ and by (2.3). By $\eqref{e:4eqn}_{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 $\eqref{e:4eqn}_{4}$. Finally, (5.3) and $\eqref{e:4eqn}_{4}$ imply $\eqref{e:4eqn}_{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}^{+})\neq\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\neq\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 $\eqref{e:proportionality}_{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}$

[TABLE]

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 $\eqref{e:4eqn}_{4}$ because

[TABLE]

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 $\eqref{e:4eqn}_{4}$ and (4.11) imply $\eqref{e:4eqn}_{3}$, namely $k_{j}=0$. By Lemma 5.1 (i) and Proposition 4.8 it remains to prove (5.1) and (4.10). 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.14) and (5.9) we have

[TABLE]

Therefore we conclude that (5.1) holds. Finally, (4.10) follows immediately from (5.1), $\eqref{e:4eqn}_{3}$ and $\eqref{e:4eqn}_{4}$. ∎

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

[TABLE]

for some constants $v_{h},\delta_{h}>0$. Denote

[TABLE]

We notice that now we have

[TABLE]

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

[TABLE]

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^{-})\neq 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_{j}v_{j}\delta_{j}^{-1}=c_{i}v_{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 $\eqref{e:Mugnolo}_{2}$ is equivalent to $\eqref{e:4eqn}_{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]

[TABLE]

with no further mention. The case when only $\eqref{e:proportionality}_{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

[TABLE]

In particular, by $\eqref{e:Gspeedg}_{2}$

[TABLE]

Moreover, $g_{h}(\ell)=v_{h}\,\ell\,[\ell_{h}^{+}+\ell_{h}^{-}-\ell]$ implies

[TABLE]

and therefore (3.16) becomes

[TABLE]

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

[TABLE]

Proposition 6.1**.**

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

[TABLE]

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), $\eqref{eq:ch}_{1}$ and (4.8) are satisfied if and only if for any $h\in\mathsf{H}$ and $j\in\mathsf{J}$

[TABLE]

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\emph{S}\doteq\bigcap_{h\in\mathsf{H}}(\ell_{h}^{-},\ell_{h}^{+})\neq\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}\neq 0$ for $h\in\mathsf{H}$, namely, the traveling wave is completely non-stationary. By (4.7), (6.1) and (6.2), from (4.9) we deduce

[TABLE]

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

[TABLE]

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:

1. (i)

$\ell_{1}^{+}+\ell_{1}^{-}\neq 1$; 2. (ii)

$\ell_{1}^{\pm}\in\bigcap_{j\in\mathsf{J}}\mathcal{L}_{j}^{c}$; 3. (iii)

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

[TABLE]

where $k_{j}$ is defined in (6.5) with $\ell_{j}^{\pm}$ being solutions to

[TABLE]

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}^{-})\neq 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,

[TABLE]

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

[TABLE]

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) are then real numbers and $\ell_{j}^{\pm}$ given in (6.8) 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}^{+})\neq f_{j}(\ell_{j}^{-})$ because $f_{1}(\ell_{1}^{+})\neq 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}^{-}\neq 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.16), namely,

[TABLE]

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

[TABLE]

We point out that the above expression of $\ell_{j}$ is deduced from (4.17) by applying (6.5); moreover (6.9) is deduced from (4.16) by applying (6.3). Recall that both fractions in (6.9) are meant as in (4.14). Since

[TABLE]

we have that (6.9) is equivalent to

[TABLE]

To conclude now that the above condition is equivalent to (6.6) it is sufficient to recall that by Lemma 4.9 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.14) and (6.24).

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

[TABLE]

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

[TABLE]

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

[TABLE]

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

[TABLE]

solves (6.12) 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.13) to problem (6.12).

Proposition 6.4**.**

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

[TABLE]

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

[TABLE]

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

[TABLE]

Proof.

By Theorem 3.2, any solution to (6.12) has the form (6.15) 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.14)-(6.16). Straightforward computations show that in the present case (4.1) can be written as

[TABLE]

where $\zeta_{h}(t)\doteq\exp z_{h}(t)$, for $z_{h}(t)\doteq\frac{v_{h}}{\delta_{h}}\,(\ell_{h}^{+}-\ell_{h}^{-})(c_{h}t+\sigma_{h})$, $h\in\mathsf{H}$. By Proposition 4.6 we have

[TABLE]

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

[TABLE]

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

[TABLE]

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

[TABLE]

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

[TABLE]

In both cases we proved that (4.1) is equivalent to (6.14)-(6.16); this concludes the proof. ∎

Remark 6.5**.**

Consider conditions $\eqref{eq:AmonAmarth}_{1}$, $\eqref{eq:AmonAmarth}_{2}$ and (6.14). Any two of them implies the third one.

Proposition 6.6**.**

Assume (6.11). Problem (2.1)-(2.2) admits a (completely) non-stationary traveling wave satisfying (2.5) if and only if (6.10) 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.10) implies (6.14), by Proposition 6.4 any (completely) non-stationary traveling wave $\rho$ has a profile of the form (6.15)-(6.16).

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

[TABLE]

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

[TABLE]

If $\ell_{h}^{-}=0$, then the function

[TABLE]

solves (6.19) because $\ell_{h}^{-}<\ell_{h}^{+}$; by (3.6) we have I_{h}=\mathopen{}\mathclose{{}\left(-\frac{\delta_{h}}{v_{h}}\ln 2,\infty}\right). If $\ell_{h}^{-}>0$, then $I_{h}=\mathbb{R}$ and the function $\psi_{h}$ implicitly given by

[TABLE]

solves (6.19) 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

[TABLE]

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

[TABLE]

Lemma 6.7**.**

Assume (6.18). If $\ell_{1}^{+}+\ell_{1}^{-}\neq 1$, then condition (6.6) is equivalent to

[TABLE]

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.24)1, the latter is equivalent to (6.24)2 by (6.5)4, because $\kappa_{j}=0$. ∎

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

[TABLE]

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

[TABLE]

Proposition 6.8**.**

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

[TABLE]

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

[TABLE]

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

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

[TABLE]

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.24) by Lemma 6.7 and (6.7) as (6.25). Then, (6.24) and (6.25) with $\ell_{h}^{-}=0$, $h\in\mathsf{H}$, reduce to the relation among the end states

[TABLE]

By (6.29) we obtain $v_{1,j}^{2}/\delta_{1,j}=\alpha_{1,j}\,v_{1,j}\,\ell_{1}^{+}/\ell_{j}^{+}$ and then $\eqref{e:values}_{1}$; by plugging $\eqref{e:values}_{1}$ into $\eqref{eq:Deftones}_{2}$ we get $\eqref{e:values}_{2}$ and then the claim.

Assume there is a degenerate non-stationary traveling wave; then $\ell_{1}^{+}$ satisfies (6.28)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}\neq v^{3}_{1,j}$ for every $j\in\mathsf{J}$. The former case is equivalent to (6.10). In the latter case we can explicitly compute $\ell_{1}^{+}$ by $\eqref{e:values}_{2}$ for any $j\in\mathsf{J}$ and impose the constraint $0<\ell_{1}^{+}<1$, namely,

[TABLE]

A direct computation shows that this is equivalent to (6.26). In conclusion, either condition (6.10) or (6.26) is necessary for the existence of a non-stationary traveling wave with $\ell_{1}^{-}=0$.

Conversely, assume condition (6.10). In this case $\alpha_{1,j}\,\delta_{1,j}^{2}=\alpha_{1,j}\,v_{1,j}^{4}=v_{1,j}^{3}$. Then (6.28)2 is trivially satisfied for every $\ell_{1}^{+}\in(0,1)$ and from (6.28)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.26). In this case the values for $\ell_{1}^{+}$ and $\ell_{j}^{+}$ in (6.27) are well defined since $v^{3}_{1,j}\neq\alpha_{1,j}\delta_{1,j}$ and they are the unique solution to (6.28). In particular, condition (ii) in Proposition 6.2 is automatically satisfied. By the estimates in (6.26) 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.27). Furthermore, by Lemma 5.1, condition (2.5) implies (5.3), which is precluded by (6.26). Hence, the traveling wave does not satisfy (2.5).

At last, by Theorem 3.2, any solution to (6.19) has the form (6.22). By (4.13), that in the present case becomes

[TABLE]

and the regularity of $\psi_{h}$ defined in (6.20), we have

[TABLE]

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

The following result treats the non-degenerate case.

Proposition 6.9**.**

Assume (6.18). Problem (2.1)-(2.2) admits a non-degenerate (completely) non-stationary traveling wave if and only if condition (6.10) 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.22) with $\psi_{h}$ implicitly defined by (6.21) and $\sigma_{h}\in\mathbb{R}$, $h\in\mathsf{H}$, satisfying (6.23).

Proof.

Assume that there is a non-degenerate non-stationary traveling wave; then $\ell_{h}^{-}\neq 0$ and $1\neq\ell_{h}^{+}+\ell_{h}^{-}$, $h\in\mathsf{H}$. Moreover, by Proposition 6.2, condition ($\mathcal{T}_{q}$) is satisfied, where (6.6) becomes (6.24) by Lemma 6.7 and (6.7) is (6.25). When dividing (6.25) by (6.24)2 we obtain

[TABLE]

By adding the above relations we have $\ell_{1}^{-}+\ell_{1}^{+}=\ell_{j}^{-}+\ell_{j}^{+}$, hence

[TABLE]

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

Conversely, assume (6.10). Then (6.24) and (6.25) write

[TABLE]

The same computations as before give that if we impose $\ell_{h}^{-}\neq 0$ and $1\neq\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.19) has the form (6.22). Fix $j\in\mathsf{J}$. By (2.5) we have (5.1), namely

[TABLE]

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

[TABLE]

By dividing the above equalities and taking the logarithm we get

[TABLE]

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

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

[TABLE]

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

[TABLE]

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

[TABLE]

Therefore (3.16) becomes

[TABLE]

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

[TABLE]

Proposition 7.1**.**

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

[TABLE]

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}$

[TABLE]

Hence \ell_{1}^{+}=f_{r}^{-1}\mathopen{}\mathclose{{}\left(-\ell_{1}^{-}\ln(\ell_{1}^{-})}\right) 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\emph{S}\doteq\bigcap_{h\in\mathsf{H}}(\ell_{h}^{-},\ell_{h}^{+})\neq\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

[TABLE]

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

[TABLE]

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:

1. (i)

$\ell_{1}^{-}\ln(\ell_{1}^{-})\neq\ell_{1}^{+}\ln(\ell_{1}^{+})$; 2. (ii)

$\ell_{1}^{\pm}\in\bigcap_{j\in\mathsf{J}}\mathcal{L}_{j}^{c}$; 3. (iii)

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

[TABLE]

where $g_{h}$ is given in (7.2), $c_{1,j}$ in (7.4), $A_{1,j}$ in (4.9)2 and $k_{j}$ in (4.9)3, with $\ell_{j}^{\pm}$ being solutions to

[TABLE]

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.16) is equivalent to

[TABLE]

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.9)2 and $k_{j}$ in (4.9)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.9, 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:

1. (i)

$\ell_{1}^{-}=0$; 2. (ii)

$\ell_{j}^{-}=0$* for all $j\in\mathsf{J}$;* 3. (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

[TABLE]

Therefore, by (7.4), (4.9)2 and (4.9)3 we have that

[TABLE]

Condition (7.5) can be written as

[TABLE]

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

[TABLE]

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

[TABLE]

By Lemma 7.3 we have either $\ell_{h}^{-}=0$, $h\in\mathsf{H}$, or $\ell_{h}^{-}\neq 0$, $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.10) holds true or

[TABLE]

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

[TABLE]

and do not satisfy (5.3).

Proof.

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

[TABLE]

Hence (7.6) can be written as

[TABLE]

and therefore (7.5) becomes

[TABLE]

namely

[TABLE]

System (7.10)-(7.11) admits a solution if and only if either (6.10) 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.10). 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}\neq\ell^{+}_{h}$. We introduce the change of variable

[TABLE]

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

[TABLE]

where

[TABLE]

Furthermore, equation (A.2) has a wavefront solution $\psi_{h}$ from $1$ to [math] 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

[TABLE]

and $\varphi_{h}$ is obtained by $\psi_{h}$ by the change of variable (A.1), i.e.

[TABLE]

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

[TABLE]

By the definition of $G_{h}$ we have

[TABLE]

Then, inequality (A.4) is equivalent to

[TABLE]

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 [math] 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

[TABLE]

where

[TABLE]

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

[TABLE]

which implies $\psi_{h}^{\prime}<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}^{\prime}>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}^{\prime}$ 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

[TABLE]

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

[TABLE]

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

[TABLE]

By applying (A.3) we conclude that $\varphi_{h}(\xi)=\ell_{h}^{-}$ for $\xi\leq\nu_{h}^{-}$ and the estimates in (3.8) are satisfied. Furthermore, by the change of variables (A.1), we obtain that

[TABLE]

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

[TABLE]

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

[TABLE]

By applying (A.3) we conclude that $\varphi_{h}(\xi)=\ell_{h}^{+}$ for $\xi\geq\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 the Gruppo Nazionale per l’Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). 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 project Balance Laws in the Modeling of Physical, Biological and Industrial Processes of GNAMPA.