Discontinuous Galerkin finite element scheme for a conserved higher-order traffic flow model by exploring Riemann solvers
Introduction
A Riemann solver, or numerical flux, is sufficient for designing a first-order accurate scheme for hyperbolic conservation laws. The Godunov, Engquist-Osher (EO), and Lax-Friedrichs (LF) are well-known numerical fluxes for solving scalar equations, such as the Lighthill-Whitham-Richards (LWR) model [16], [21] and the Burgers equation, which give rise, respectively, to the Godunov, EO, and LF schemes [4], [6], [23], [24] (see also [2], [3], [5], [8], [9], [11], [15], [26], [29], [35] for studies of the LWR or Multi-class LWR models). To solve the conserved higher-order (CHO) model [33], which is an extension of the LWR model, the aforementioned fluxes were extended with a so-called traffic flow (TF) flux to design a number of first-order accurate schemes [30]. Although these extended fluxes were able to reproduce stable and convergent numerical solutions, considerable numerical viscosities or errors were observed.
This sequence simplifies the expression of these fluxes and the computation becomes more efficient; however, the shock profile of the reproduced wide moving jam in the CHO model becomes smoother and the backward moving wave becomes faster, which suggests an increase in numerical viscosity. It seems that the improvement in efficiency from computing a simpler numerical flux is counteracted by an increase in numerical viscosity. If a first-order accurate scheme is adopted, the grid needs to be refined to improve the solution. Generally, improving the scheme’s accuracy should suppress the redundant numerical viscosity in a simpler flux. More precisely, numerical fluxes (either simple or complicated) are expected to produce similarly accurate solutions in a higher-order scheme. In this case, using a simpler flux will improve the scheme’s efficiency.
Following this insight, this paper designs a discontinuous Galerkin (DG) scheme for solving the CHO model by exploiting the numerical fluxes extended from the Godunov, EO, LF and TF fluxes [19], [30]. The DG scheme is found to be more accurate than the first-order scheme in that it is able to suppress differences between the explored solvers when generating numerical viscosities. Moreover, we find that the performance of the second-order DG scheme is similar to that of the fifth-order WENO schemes, a similarity that is due to the existence of a shock profile in the solution, although the latter scheme is expected to be more accurate in generating smooth solutions. These findings indicate that the second-order DG scheme is adequate for the model, especially when a simpler flux with a lower computational cost is adopted. To enhance the argument, an analytical solution for a wide moving jam is constructed for comparison with the numerical solutions.
The remainder of this paper is organized as follows. Section 2 discusses the CHO model, together with its solution for a wide moving jam. In Section 3, the extended Godunov, EO, LF and TF numerical fluxes are derived by applying the Riemann solution to the model’s homogeneous system. In Section 4, the DG schemes are constructed using these fluxes. In Section 5, numerical solutions to the wide moving jam are derived through simulation, and the solutions produced by the first-order accurate scheme, the fifth-order accurate WENO scheme and the second-order accurate DG scheme are compared with the analytical solution. Section 6 concludes the paper.
Section snippets
Model equations
Taken as a continuum, the mass conservation of traffic flow is described through the following partial differential equation:where and are the density and velocity in location x at time t, respectively, and denotes the flow at . Assume that there is a determined equilibrium velocity-density relationship . Then, the substitution of v in Eq. (1) by this equation gives rise to the following well-known LWR model [16], [21]:
Riemann solvers for the homogeneous system
System (3), (4) is rewritten in the following vector form:with and . To derive a numerical flux, the Riemann problem in the homogeneous systemis of interest. The problem is actually set with the initial valueswhere and .
The discontinuous Galerkin finite element scheme
By dividing the computational interval into cells, , with , we multiply system (12) by a test function and integrate the resultant equations over . By applying the integration of the parts to the second term, we have
Numerical simulation
In system (3), (4), we generally set the relaxation time as , and the velocity-density relationships aswhere , and veh./m. Similar settings were chosen in [19], [30]. To ensure the numerical stability of a scheme, the time step should satisfy the following CFL condition:where .
The computational interval is set as , with
Conclusions
The CHO model is composed of the mass and pseudo-mass conservations, which are favorably consistent with the LWR model. Underlying this consistency is the finding that the Riemann solver for the model can be expressed by the classical Godunov flux, which is based on a Riemann invariant that remains unchanged through propagation from the initial state on the upstream side to the interface. Accordingly, the classical EO and LF fluxes, together with a TF flux, are extended, allowing us to obtain a
Acknowledgements
This study was jointly supported by grants from the National Natural Science Foundation of China (11072141, 11272199), the National Basic Research Program of China (2012CB725404), the Shanghai Program for Innovative Research Team in Universities, and a National Research Foundation of Korea grant funded by the Korean government (MSIP) (NRF-2010-0029446).
References (35)
- et al.
A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier–Stokes equations
J. Comput. Phys.
(1997) - et al.
On the distribution schemes for determining flows through a merge
Transport. Res. Part B
(2003) - et al.
Semi discrete discontinuous Galerkin methods and stage-exceeding-order, strong-stability-preserving Runge–Kutta time discretizations
J. Comput. Phys.
(2007) - et al.
The Aw-Rascle and Zhang’s model: vacuum problems, existence and regularity of the solutions of the Riemann problem
Transport. Res. Part B
(2007) An improved macroscopic model of traffic flow: derivation and links with the Lighthill–Whitham model
Math. Comput. Modell.
(2002)- et al.
Riemann solver for a kinematic wave traffic model with discontinuous flux
J. Comput. Phys.
(2013) - et al.
A conserved higher-order anisotropic traffic flow model: description of equilibrium and non-equilibrium flows
Transport. Res. Part B: Method.
(2009) - et al.
A semi-discrete model and its approach to a solution for wide moving jam in traffic flow
Physica A
(2012) - et al.
Hyperbolicity analysis of polydisperse sedimentation models via a secular equation for the flux Jacobian
SIAM J. Appl. Math.
(2010) - et al.
A family of numerical schemes for kinematic flows with discontinuous flux
J. Eng. Math.
(2008)