Papers › A Positivity-Preserving Limiting Strategy for Locally-Implicit Lax-Wendroff...
A Positivity-Preserving Limiting Strategy for Locally-Implicit Lax-Wendroff Discontinuous Galerkin Methods
Camille Felton, Mariana Harris, Caleb Logemann, Stefan Nelson, Ian Pelakh, James A. Rossmanith
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Nonlinear hyperbolic conservation laws admit singular solutions such as shockwaves (discontinuities in conserved variables), rarefaction waves (discontinuities in derivatives), and vacuum states (loss of strong hyperbolicity). When ostensibly high-order numerical methods are applied in such solution regimes, unphysical oscillations present themselves that can lead to large errors and a breakdown of the numerical simulation. In this work we develop a new Lax-Wendroff discontinuous Galerkin (LxW-DG) method with a limiting strategy that keeps the solution non-oscillatory and positivity-preserving for relevant variables, such as height in the shallow water equations and density and pressure in the compressible Euler equations. The proposed LxW-DG scheme updates the solution over each time-step with a locally-implicit predictor followed by an explicit corrector. The locally-implicit prediction phase is formulated in terms of primitive variables, which greatly simplifies the solver. The resulting system of nonlinear algebraic equations are approximately solved via a Picard iteration, where the number of iterations is equal to the order of accuracy of the method. The correction phase is an explicit evaluation formulated in terms of conservative variables in order to guarantee numerical conservation. In order to achieve full positivity-preservation, limiting is required in both the prediction and correction steps. The resulting scheme is applied to several standard test cases for the shallow water and compressible Euler equations. All of the presented examples are written in a freely available open-source Python code.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections