Bound-preserving OEDG schemes for Aw-Rascle-Zhang traffic models on networks

Wei Chen,Shumo Cui,Kailiang Wu,Tao Xiong
2024-09-25
Abstract:Physical solutions to the widely used Aw-Rascle-Zhang (ARZ) traffic model and the adapted pressure (AP) ARZ model should satisfy the positivity of density, the minimum and maximum principles with respect to the velocity $v$ and other Riemann invariants. Many numerical schemes suffer from instabilities caused by violating these bounds, and the only existing bound-preserving (BP) numerical scheme (for ARZ model) is random, only first-order accurate, and not strictly conservative. This paper introduces arbitrarily high-order provably BP DG schemes for these two models, preserving all the aforementioned bounds except the maximum principle of $v$, which has been rigorously proven to conflict with the consistency and conservation of numerical schemes. Although the maximum principle of $v$ is not directly enforced, we find that the strictly preserved maximum principle of another Riemann invariant $w$ actually enforces an alternative upper bound on $v$. At the core of this work, analyzing and rigorously proving the BP property is a particularly nontrivial task: the Lax-Friedrichs (LF) splitting property, usually expected for hyperbolic conservation laws and employed to construct BP schemes, does not hold for these two models. To overcome this challenge, we formulate a generalized version of the LF splitting property, and prove it via the geometric quasilinearization (GQL) approach [Kailiang Wu and Chi-Wang Shu, SIAM Review, 65: 1031-1073, 2023]. To suppress spurious oscillations in the DG solutions, we employ the oscillation-eliminating (OE) technique, recently proposed in [Manting Peng, Zheng Sun, and Kailiang Wu, Mathematics of Computation, in press], which is based on the solution operator of a novel damping equation. Several numerical examples are included to demonstrate the effectiveness, accuracy, and BP properties of our schemes, with applications to traffic simulations on road networks.
Numerical Analysis
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to solve the problem of numerical simulation instability in traffic flow models, especially for the Aw - Rascle - Zhang (ARZ) model and the Adaptive Pressure (AP) ARZ model. Specifically, these problems include: 1. **Maintaining physical boundary conditions**: Physical solutions should satisfy the positivity of density, the minimum and maximum principles of velocity \(v\), and the boundary conditions of other Riemann invariants. Many existing numerical methods generate instabilities due to violations of these boundary conditions. 2. **Limitations of existing BP schemes**: The only existing boundary - preserving numerical scheme (for the ARZ model) is stochastic, only first - order accurate and not strictly conservative. Therefore, it is necessary to develop high - order, strictly conservative numerical methods that can preserve boundaries. 3. **Inapplicability of Lax - Friedrichs (LF) splitting property**: For the ARZ and AP ARZ models, the LF splitting property, which is usually used to construct BP schemes, does not hold, which increases the difficulty of designing BP schemes. 4. **Conflict with the maximum principle of velocity \(v\)**: It has been found that the maximum principle of \(v\) conflicts with the consistency and conservativeness of numerical schemes. Although directly enforcing the maximum principle of \(v\) is not feasible, by strictly maintaining the maximum principle of another Riemann invariant \(w\), an upper bound on \(v\) can be indirectly imposed, thereby alleviating the velocity overshoot phenomenon. ### Main contributions of the paper To address the above challenges, the paper makes the following key contributions: 1. **Rigorous theoretical analysis**: Systematically prove the key properties of the invariant domain induced by constraint conditions, including convexity and the generalized LF splitting property. These theoretical efforts involve technical estimates and the Geometric Quasi - Linearization (GQL) method, which are novel and non - trivial. 2. **First proposed high - order BP DG scheme**: Based on the above theoretical analysis, a high - order BP DG scheme applicable to the original ARZ model and the AP ARZ model is developed, which can maintain the maximum / minimum principles of all Riemann invariants and the positivity of density. 3. **Effective velocity upper - bound control**: Although the maximum principle of \(v\) is not directly dealt with, by strictly maintaining the maximum principle of \(w\), an alternative upper bound on the value of \(v\) is effectively imposed, significantly alleviating the velocity overshoot phenomenon in numerical results, especially in regions close to the vacuum state. 4. **Local invariant domain estimation**: In addition to the global invariant domain, a method for estimating the local invariant domain is also provided. Numerical experiments show that the local BP - OEDG scheme is more robust than the global BP scheme, especially in solutions containing waves of different amplitudes. 5. **Oscillation suppression**: Incorporate a non - invasive scale - invariant OE process into the DG scheme to suppress numerical oscillations while maintaining the high - order accuracy of the DG scheme. Through these contributions, the paper provides a more stable and accurate solution for the numerical simulation of traffic flow models.