Research
Overview
HAMILTON-JACOBI EQUATIONS AND SPATIAL DISCONTINUITIES
Hamilton-Jacobi equations on junctions
We consider first order Hamilton-Jacobi (HJ) equations posed on a "junction", that is to say the union of a finite number of half-lines with a unique common point. For this continuous HJ problem, we propose a finite difference scheme and prove two main results. As a first result, we show bounds on the discrete gradient and time derivative of the numerical solution. Our second result is the convergence (for a subsequence) of the numerical solution towards a viscosity solution of the continuous HJ problem, as the mesh size goes to zero. When the solution of the continuous HJ problem is unique, we recover the full convergence of the numerical solution. We apply this scheme to compute the densities of cars for a traffic model. We recover the well-known Godunov scheme outside the junction point and we give a numerical illustration.
References:
Costeseque, G., Lebacque, J.P., Monneau, R., (2015). A convergent scheme for Hamilton-Jacobi equations on a junction: application to traffic, Numerische Mathematik, 129 (3), 405-447.
Costeseque, G., Lebacque, J.P., (2012). Intersection modeling using a convergent scheme based on Hamilton-Jacobi equation, Proceedings of 15th EWGT, Social & Behavioral Sciences, Volume 54, 736-748.
Costeseque, G., Lebacque, J.P., (2014). Discussion about traffic junction modelling: conservation laws VS Hamilton-Jacobi equations, Discrete Cont. Dyn. Syst. S, 7 (3), 411-433.
Note: This work has been extended in:
Guerand, J., & Koumaiha, M. (2019). Error estimates for a finite difference scheme associated with Hamilton–Jacobi equations on a junction. Numerische Mathematik, 142(3), 525-575.
Hamilton-Jacobi equations and source terms
This paper revisits the variational theory (VT) of traffic flow, now under the presence of continuum lateral inflows and outflows to the freeway say Eulerian source terms. It is found that a VT solution can be easily exhibited only in Eulerian coordinates when source terms are exogenous meaning that they only depend on time and space, but not when they are a function of traffic conditions, as per a merge model. In discrete time, however, these dependencies become exogenous, which allowed us to propose improved numerical solution methods. In Lagrangian and vehicle number-space coordinates, VT solutions may not exist even if source terms are exogenous.
Reference:
Laval, J. A., Costeseque, G., Chilukuri, B.R., (2016). The impact of source terms in the variational representation of traffic flow, Transportation Research part B: Methodological, 94, 204-216.
Multilanes section near a diverge
The aim of this paper is to propose a new event-based mesoscopic approach to model multi-class traffic flow on multi-lane road sections. The mesoscopic model was first proposed by Leclercq and Bécarie (2012) and turns out to be equivalent to the resolution of the seminal LWR model in lagrangian-space coordinates n − x. It is fully consistent at a macroscopic scale with the LWR model while keeping track of individual vehicles. Our model is built on Hamilton-Jacobi equations which have been proven to provide an efficient framework in traffic flow modeling for exact numerical methods at a low computational cost. The paper revisits the multi-class problem with a continuous moving bottleneck approach (instead of a sequential one), introducing a capacity drop parameter for multi-lane sections. It also overhauls the Daganzo diverge model with a relaxed FIFO assumption.
Reference:
Costeseque, G., Duret, A., (2016). Mesoscopic multiclass traffic flow modeling on multi-lane sections, Proceedings of the Transportation Research Board 95th Annual Meeting (2016), Washington DC, USA.
SCALING IN TRAFFIC MODELS
From micro to macro without reaction time
We set a rigorous link between microscopic and macroscopic traffic models by homogenization of a system of Ordinary Differential Equations (ODEs) that describe the dynamics of vehicles. At the microscopic level, each vehicle behaves according to a given speed-spacing function depending on its leading vehicle. After a proper asymptotic rescaling that can be seen as unzooming, we show that solutions of these systems of ODEs converge to the solution of a macroscopic traffic flow model, based on Hamilton-Jacobi equation on one dimensional domain.
Reference:
Costeseque, G., (2011). Modélisation du trafic routier : passage du microscopique au macroscopique, Master Thesis, ED MEGA, Université Lyon.
Disclaimer: This work is mainly based on Monneau's lectures. It has been extended in:
Firozaly, J., (2016). Homogenization of a 1D pursuit law with delay and a counter-example, preprint.
From micro to macro with reaction time
We derive first order continuum traffic flow models from a microscopic delayed follow-the-leader model. This micro model is based on an optimal velocity function and a reaction time parameter. The corresponding macroscopic formulations in Eulerian or Lagrangian coordinates result in first order convection-diffusion equations. More precisely, the convection is described by the optimal velocity while the diffusion term depends on the reaction time. A linear stability analysis for homogeneous solutions of both continuous and discrete models are provided. The behavior of the novel model is illustrated thanks to numerical simulations. Transitions to collision-free self-sustained stop-and-go dynamics are obtained if the reaction time is sufficiently large. The results show that the dynamics of the microscopic model can be well captured by the macroscopic equations. For non–zero reaction times we observe a scattered fundamental diagram. The scattering width is compared to real data.
Reference:
Tordeux, A., Costeseque, G., Herty, M., Seyfried, A., (2018). From traffic and pedestrian follow-the-leader models with reaction time to first order convection-diffusion flow models, SIAM Journal on Applied Mathematics, 78, 63-79.
From micro to macro with multi-anticipation
We deal with a macroscopic model of multi-anticipative car-following behaviour i.e. driving behaviour taking into account several vehicles ahead. Some empirical studies have suggested that drivers not only react to the closest leader vehicle but also anticipate on traffic conditions further ahead. Using a recent mathematical result of homogenization for a general class of car-following models (and also available for multi-anticipative models), we deeply investigate the effects of multi-anticipation at the microscopic level on the macroscopic traffic flow. To investigate multi-anticipation behaviour may be fundamental to understand better cooperative traffic flow dynamics.
Reference:
Costeseque, G., Lebacque, J.P., (2014). Multi-anticipative car-following behaviour: macroscopic modeling, Proceedings of the Traffic and Granular Flow'13 conference (TGF'13), Julich, Germany.
From macro to very large areas
This paper is concerned with the numerical homogenization of Hamilton-Jacobi equations posed on a periodic network embedded in R^2. For this continuous HJ problem, we propose a finite difference scheme which was previously designed for HJ equations on a single junction. We provide some qualitative properties of the effective Hamiltonian. We also apply a derived scheme to compute the densities of vehicles for a traffic flow model on the network and we recover some macroscopic features obtained in traffic flow theory. We finally provide numerical computations of the effective Hamiltonian in the case of two consecutive traffic lights on R.
Reference:
Costeseque, G., Ghorbel, A., (2017). Numerical homogenization of Hamilton-Jacobi equations on networks and applications to traffic, working paper.
MACRO SECOND ORDER MODELS
Variational formulation for second order models
This paper deals with numerical methods providing semi-analytic solutions to a wide class of macroscopic traffic flow models for piecewise affine initial and boundary conditions. In a very recent paper, a variational principle has been proved for models of the Generic Second Order Modeling (GSOM) family, yielding an adequate framework for effective numerical methods. Any model of the GSOM family can be recast into its Lagrangian form as a Hamilton-Jacobi equation (HJ) for which the solution is interpreted as the position of vehicles. This solution can be computed thanks to Lax-Hopf like formulas and a generalization of the inf-morphism property. The efficiency of this computational method is illustrated through a numerical example and finally a discussion about future developments is provided.
Reference:
Costeseque, G., Lebacque, J.P., (2014). A variational formulation for higher order macroscopic traffic flow models: numerical investigation, Transportation Research part B: Methodological, 70, 112-133.
Bounded acceleration: queue estimation on arterials
We consider a macroscopic traffic model able to reproduce the boundedness of the vehicles acceleration. It is based on a modified Lighthill-Whitham-Richards (LWR) model recast as a Hamilton-Jacobi Partial Differential Equation for which there exists explicit solution methods. We developed an optimization-based framework which is able to be extended to large scale networks in order to estimate queue lengths on arterial road networks, taking into account data from classical density or flow sensors and mobile sensors. We validated our results on real data extracted from the NGSIM Lankershim Boulevard dataset, comparing the estimation result from LWR model with and without the bounded acceleration constraint. Comparing to previous method without considering bounded acceleration, our method improves the queue lengths estimation precision under most cases.
References:
Canepa, E.S., Costeseque, G., Claudel, C.G., (2016). A queue estimation method aware of bounded acceleration: application to arterials, Proceedings of the TRISTAN conference (2016), Aruba, the Netherlands.
Qian, K., Costeseque, G., Canepa, E.S. and Claudel, C.G. (2019). Optimization-based queue length estimation on arterial networks using a macroscopic traffic model enforcing vehicles bounded acceleration, submitted, 49 pages.
Software:
Coupled PDE-ODE model for modeling bounded acceleration
This paper details a new macroscopic traffic flow model accounting for the boundedness of traffic acceleration, which is required for physical realism. Our approach relies on the coupling between a scalar conservation law, which refers to the seminal LWR model, and a system of Ordinary Differential Equations describing the trajectory of accelerating vehicles, which we treat as moving constraints. We propose a Wave Front Tracking Algorithm to construct approximate solutions. We use this algorithm to prove the existence of solutions to the associated Cauchy Problem, and provide some numerical examples illustrating the solution behaviour.
References:
Laurent-Brouty, N., Costeseque, G., Goatin, G. (2018). A coupled PDE-ODE model for bounded acceleration in macroscopic traffic flow models, 15th IFAC Symposium on Control in Transportation Systems (CTS 2018), Savona, Italy.
Laurent-Brouty, N., Costeseque, G. and Goatin, P. (2021). A macroscopic traffic flow model accounting for bounded acceleration, SIAM Journal on Applied Mathematics, 81(1), 173-189.
SECOND ORDER MODELS AND SPATIAL DISCONTINUITIES
Aw-Rascle-Zhang traffic flow model on junctions
This article deals with macroscopic traffic flow models on a road network. More precisely, we consider coupling conditions at junctions for the Aw-Rascle-Zhang (ARZ) second order model consisting of a hyperbolic system of two conservation laws. These coupling conditions conserve both the number of vehicles and the mixing of Lagrangian attributes of traffic through the junction. The proposed Riemann solver is based on assignment coefficients, multi-objective optimization of fluxes and priority parameters. We prove that this Riemann solver is well posed in the case of special junctions, including 1-to-2 diverge and 2-to-1 merge.
Reference:
Kolb, O., Costeseque, G., Goatin, G., Göttlich, S. (2018). Pareto-optimal coupling conditions for the Aw-Rascle-Zhang traffic flow model at junctions, SIAM Journal on Applied Mathematics, 78(4), 1981-2002.
Lagrangian GSOM traffic flow models on junctions
This paper is concerned with a generic class of second order models (GSOM) on networks. Such higher order models can be more easily solved in a Lagrangian framework whose coordinates move with the traffic. The difficulty is to deal with Eulerian-fixed discontinuities such as junctions. The aim of this work is twofold: first, to propose adapted junction models for second order macroscopic traffic flow models and second, to solve the resulting model in a moving framework. Our numerical methodology is based on a finite difference scheme that bridges the gap with a microscopic description of the traffic.
References:
Costeseque, G., Lebacque, J.P., Khelifi, A., (2015). Lagrangian GSOM traffic flow models on junctions, Proceedings of the 8th Vienna International Conference on Mathematical Modelling, Vienna, Austria.
Khelifi, A., Haj-Salem, H., Lebacque, J. P., Lioris, J., Costeseque, G., Nabli, L. (2016). Contribution à la modélisation du GSOM Lagrangien. Academic Journal of Civil Engineering, 34(1), 1167-1174.
Costeseque, G., Lebacque, J.P., Khelifi, A., (2017). Higher order macroscopic traffic flow models on networks: a Lagrangian perspective, working paper.