We describe a numerical technique for the solution of macroscopic traffic flow models on networks of roads. On individual roads, we consider the standard Lighthill-Whitham-Richards model which is discretized using the discontinuous Galerkin method along with suitable limiters.
In order to solve traffic flows on networks, we construct suitable numerical fluxes at junctions based on preferences of the drivers. We prove that our semi-discrete DG solution is L2 stable on several types of networks.
We present numerical experiments, including a junction with complicated traffic light patterns with multiple phases.