2
depends on both the density and the spatial derivative of the density (see also [11], [20], [21], [15]),
a feature that is rarely studied in the literature where typically the diffusion coefficient exclusively
depends on the density, see for instance [11], [25], [26], [27], [29], [30], [33] and references therein.
The nonlinear heat equation under consideration is inspired by a second order macroscopic traffic
flow model that was recently derived in [16] by the design of cruise controllers for autonomous
vehicles and has several fluid-like characteristics. Macroscopic traffic flow models have been widely
used as they are able to capture the collective behavior of vehicles in a traffic stream and can express
certain relationships between traffic flow, density and mean speed. Macroscopic traffic flow models
can be in general distinguished between first-order models, governed by the continuity equation, and
second-order models, which include an additional differential equation for the speed (see for instance
[2], [12], [34]).
In the present work, it is first shown that the solutions of the corresponding second-order
macroscopic model can be approximated by the solutions of this particular nonlinear heat equation.
Furthermore, inspired by the mechanical energy of the original second-order system, certain
functionals are defined, providing entropy-like conditions that characterize physically meaningful
solutions, see [9], [31]. In particular, we define the concept of a weak solution for the problem at
hand that is more demanding than the typical definition of weak solution (see [24], [31]), as it requires
that certain entropy conditions hold and, in addition, the mass is conserved. Finally, it is shown that,
under certain regularity assumptions, this type of weak solutions is equivalent to classical solutions,
under the assumption of certain entropy-like conditions (see Proposition 1).
The study of the nonlinear heat equation is inevitably related to numerical techniques in order to
obtain an approximation of its solution. More specifically, several finite-difference methods have
been proposed to study a variety of nonlinear heat equations, see for instance [3], [4], [6], [7], [8],
[10], [13], [14], and references therein. In this paper, we propose a new first-order explicit finite-
difference scheme that, in addition to being conservative (Proposition 2), it also respects the
corresponding entropy conditions (Proposition 3 and Proposition 5). Moreover, a condition on the
step size is provided that is a nonlinear version of the Courant-Friedrichs-Levy condition (Proposition
4). Finally, a link between the numerical solution produced by the proposed numerical scheme and
the corresponding weak solution is provided (Proposition 6). To demonstrate the properties of the
solutions of the nonlinear heat equation, we present two examples: (i) an academic example, where
the efficiency of the proposed numerical scheme is illustrated; and (ii) a realistic traffic simulation
scenario showing that the mean flow produced by automated vehicles is much higher than that of the
LWR (see [22], [28]).
Therefore, the contribution of the paper is summarized as follows:
• The study of a novel nonlinear heat equation that arises in traffic flow theory for automated
vehicles;
• The construction of a first-order numerical scheme that can be applied to nonlinear heat
equations with an explicit condition for numerical stability;
• The study of the connection of the weak solution for the nonlinear heat equation with the
solution produced by the proposed numerical scheme.
The structure of the paper is as follows. Section 2 is devoted to the presentation of a second-order
traffic flow model and the approximation of its solution by the solution of a particular nonlinear heat
equation. Section 3 presents the proposed numerical scheme and its properties. Section 4 presents the
two aforementioned examples. Finally, some concluding remarks are given in Section 5. All proofs
are provided in the Appendix.