pith. sign in

arxiv: 2606.31992 · v1 · pith:TZHZRKLXnew · submitted 2026-06-30 · 🧮 math.NA · astro-ph.IM· cs.NA· physics.comp-ph· physics.flu-dyn

GQL-Based Physical-Constraint-Preserving High-Order Finite Difference Schemes for Special Relativistic Hydrodynamics in Arbitrary Dimensions

Linfeng Xu , Shengrong Ding , Kailiang Wu This is my paper

Pith reviewed 2026-07-01 03:34 UTC · model grok-4.3

classification 🧮 math.NA astro-ph.IMcs.NAphysics.comp-phphysics.flu-dyn
keywords special relativistic hydrodynamicsphysical constraint preservingfinite difference schemesWENOflux limitinggeometric quasilinearizationhigh-order methods
0
0 comments X

The pith

Geometric quasilinearization turns nonlinear relativistic constraints into linear inequalities solvable by small eigenvalue problems for high-order PCP schemes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a flux-limiting framework that keeps rest-mass density, pressure, and velocity physical in high-order finite-difference simulations of special relativistic hydrodynamics. It starts from the geometric quasilinearization representation, which rewrites the nonlinear admissibility conditions as a set of linear inequalities on the conserved variables. A rational stereographic map of the normal vector then converts the worst-case limiter search into a generalized Rayleigh quotient that reduces to symmetric eigenvalue problems of size 2 by 2 in one dimension and (d+1) by (d+1) in d dimensions. Relaxed variants cut cost while preserving the guarantee. Benchmarks from one to three dimensions, including ultra-relativistic shocks and jets, show that the limiter maintains design order on smooth data and prevents unphysical states near discontinuities.

Core claim

By expressing the physical constraints of special relativistic hydrodynamics through geometric quasilinearization and solving the resulting limiter parameters via rational stereographic parameterization and small symmetric eigenvalue problems, high-order WENO finite-difference schemes can be made physical-constraint-preserving in arbitrary dimensions without iteration.

What carries the argument

Geometric quasilinearization (GQL) representation that converts nonlinear RHD constraints into linear inequalities, combined with rational stereographic parameterization of the GQL normal vector that reduces worst-case minimization to a generalized Rayleigh-quotient eigenvalue problem.

If this is right

  • The same GQL-based limiter can be inserted into any high-order finite-difference or finite-volume scheme for RHD without changing the underlying reconstruction or flux.
  • Design-order accuracy is retained for smooth solutions while discontinuities remain sharply captured and admissible.
  • The method extends directly from one to three spatial dimensions with only a modest increase in the size of the eigenvalue problems.
  • Relaxed variants lower computational cost in multidimensions while still guaranteeing physical states.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach may generalize to other hyperbolic systems whose admissibility sets admit a similar quasilinear inequality description.
  • Because the eigenvalue solves are small and independent per cell, the limiter is a natural candidate for GPU or distributed-memory parallelization.
  • If the GQL normal vector can be computed analytically from the equation of state, the entire limiter step becomes fully explicit.

Load-bearing premise

The geometric quasilinearization representation must exactly convert the nonlinear physical constraints into an equivalent family of linear inequalities.

What would settle it

A single run of the scheme on an ultra-relativistic Riemann problem in which the limiter is disabled and a cell is observed to reach negative density or superluminal velocity.

Figures

Figures reproduced from arXiv: 2606.31992 by Kailiang Wu, Linfeng Xu, Shengrong Ding.

Figure 1
Figure 1. Figure 1: Example 1: The numerical results obtained using WENO5 with the proposed PCP flux limiter. Left: 10 iso-surfaces of [PITH_FULL_IMAGE:figures/full_fig_p021_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Example 2: Numerical results for ρ, v1, p obtained using WENO5 with the relaxed PCP flux limiter at t = 0.45. (a) Example 2 (b) Example 3 [PITH_FULL_IMAGE:figures/full_fig_p023_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Spatiotemporal distribution of the proposed flux limiter coe [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Example 3: Numerical results for ρ, v1, p obtained using WENO5 or AiWENO5-Z with the proposed PCP flux limiter at t = 2 [PITH_FULL_IMAGE:figures/full_fig_p024_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Example 3: Plots of error between numerical and exact solution [PITH_FULL_IMAGE:figures/full_fig_p025_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: First 2D Riemann problem in Example 4: The contours of log [PITH_FULL_IMAGE:figures/full_fig_p027_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Same as Fig. 6 except for [PITH_FULL_IMAGE:figures/full_fig_p027_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Second 2D Riemann problem in Example 4: The contours of log [PITH_FULL_IMAGE:figures/full_fig_p028_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Example 4: Comparison of min i, j n θi+1/2, j , θi, j+1/2 o between different flux limiters over time. where the functions θˆ, w1, and w2 are given by θˆ(x, y) = 1 − ε 2 vtx 28π 2 e 1− x 2 1−v 2 c −y 2 , w1(x, y) = −y ˆf(x, y), w2(x, y) = x p 1 − v 2 c ˆf(x, y), with ˆf(x, y) = vut β 1 + β  x 2 1−v 2 c + y 2 , β = ε 2 vtx 10π 2 e 1− x 2 1−v 2 c −y 2 1.8 − ε 2 vtx 20π 2 e 1− x 2 1−v 2 c −y 2 . This setup … view at source ↗
Figure 10
Figure 10. Figure 10: Example 4: Contours of the logarithm log( [PITH_FULL_IMAGE:figures/full_fig_p029_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Example 5: Contour plots of log10(1 + |∇ρ|) at t = 19, obtained using WENO5 [PITH_FULL_IMAGE:figures/full_fig_p029_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Example 6: Top: Schlieren images of ln ρ obtained using WENO5 with our flux limiter (4.21); bottom: contours of the logarithm log(|ρex − ρre|), where ρex and ρre denote the density computed with the exact and relaxed flux limiters, respectively. 25 equally spaced contour levels from -11.5 to -0.5 are displayed. From left to right: configurations (i) at t = 30, (ii) at t = 25, and (iii) at t = 23. In the c… view at source ↗
Figure 13
Figure 13. Figure 13: Example 6: Top: Schlieren images of ln p obtained using WENO5 with our flux limiter (4.21); bottom: contours of the logarithm log(|pex − pre|), where pex and pre denote the pressure computed with the exact and relaxed flux limiters, respectively. 25 equally spaced contour levels from -11.5 to -1.2 are displayed. From left to right: configurations (i) at t = 30, (ii) at t = 25, and (iii) at t = 23. where P… view at source ↗
Figure 14
Figure 14. Figure 14: Example 7: Left: Schlieren images of ln ρ obtained using WENO5 with our flux limiter (4.21) at t = 100; bottom: Schlieren images of the logarithm log(|ρex − ρre|) at t = 100, where ρex and ρre denote the density computed with the exact and relaxed flux limiters, respectively [PITH_FULL_IMAGE:figures/full_fig_p034_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Example 7: Limiter activation frequency maps. Left: [PITH_FULL_IMAGE:figures/full_fig_p034_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Example 8: Schlieren images of ρ at t = 90, 180, 270, 360, 450 (from top to bottom). Left: on the slice z = 45; right: on the slice y = 45. The two orthogonal slices jointly demonstrate the rotational symmetry of the 3D flow field. Example 9. (3D Axisymmetric relativistic jets) This 3D example is based on the pressure-matched hot A1 model in the 2D case. The adiabatic index Γ = 4/3. The computational doma… view at source ↗
Figure 17
Figure 17. Figure 17: Example 8: 10 iso-surfaces of ρ equally spaced from 0.55 to 1.75, obtained using WENO5 with the relaxed PCP flux limiter (4.43). exactly across the x = 0 symmetry plane to display the symmetric half-domain [−7, 7]×[0, 7]×[0, 50]; this post-processing does not alter the numerical solution. Fig. 18a, 18b, and 18c present 15 iso-surfaces of ln ρ in the symmetric region [−7, 7]×[0, 7]×[0, 50] at t = 20, 40, a… view at source ↗
Figure 18
Figure 18. Figure 18: Example 9: Numerical results obtained using WENO5 with the relaxed PCP flux limiter (4.43). Left three subfigures: 15 [PITH_FULL_IMAGE:figures/full_fig_p036_18.png] view at source ↗
read the original abstract

High-order accurate simulations of special relativistic hydrodynamics (RHD) are prone to numerical breakdown if intrinsic physical constraints (positive rest-mass density/pressure and subluminal velocity) are violated near strong discontinuities. In this work, we develop a robust and efficient physical-constraint-preserving (PCP) flux-limiting framework for high-order schemes, using finite-difference WENO as a representative example. By leveraging the geometric quasilinearization (GQL) representation, which equivalently reformulates the nonlinear RHD constraints into a family of linear inequalities, we integrate a Zalesak-type Flux-Corrected Transport (FCT) update into a scalar-style limiter that acts directly on conservative variables. A critical innovation is the explicit, non-iterative determination of limiting parameters via a rational stereographic parameterization of the GQL normal vector. This technique transforms the required worst-case minimization over auxiliary variables into a generalized Rayleigh-quotient formulation, allowing the optimal parameters to be obtained by solving small symmetric eigenvalue problems ($2\times2$ in 1D; $(d+1)\times(d+1)$ in $d$ dimensions). Relaxed variants are further introduced to reduce computational costs in multidimensions while retaining the PCP guarantee. Extensive numerical benchmarks ranging from 1D to 3D, including ultra-relativistic Riemann problems and astrophysical jets, demonstrate that the proposed method robustly enforces physical admissibility, sharply resolves discontinuities, and maintains design-order accuracy for smooth solutions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript develops a physical-constraint-preserving (PCP) flux-limiting framework for high-order finite-difference WENO schemes for special relativistic hydrodynamics (RHD) in arbitrary dimensions. It employs geometric quasilinearization (GQL) to reformulate the nonlinear RHD admissibility constraints (positive density/pressure, subluminal velocity) as linear inequalities, integrates a Zalesak-type FCT update, and uses a rational stereographic parameterization of the GQL normal vector to convert the worst-case minimization into a generalized Rayleigh-quotient problem solved via small symmetric eigenvalue problems (2×2 in 1D; (d+1)×(d+1) in d dimensions). Relaxed variants reduce cost while retaining the PCP property. Extensive 1D–3D benchmarks, including ultra-relativistic Riemann problems and astrophysical jets, are presented to demonstrate enforcement of physical admissibility, sharp discontinuity capture, and retention of design-order accuracy on smooth solutions.

Significance. If the claimed exact algebraic equivalences in the GQL reformulation and stereographic parameterization hold without approximation, the work provides a non-iterative, dimension-independent PCP limiter that preserves high-order accuracy while guaranteeing admissibility. This is a meaningful advance for robust high-order RHD simulations in astrophysics, where standard methods frequently fail near strong shocks. The explicit reduction to small eigenproblems and the supporting numerical evidence constitute clear strengths.

minor comments (3)
  1. The description of the relaxed variants in multidimensions would benefit from an explicit statement of the trade-off between computational cost and the size of the admissible set (e.g., a short paragraph comparing flop counts or iteration counts for the full versus relaxed eigenproblems).
  2. Notation for the GQL normal vector and the auxiliary variables should be made uniform between the derivation of the linear inequalities and the subsequent Rayleigh-quotient formulation to improve readability.
  3. Figure captions for the 3D jet benchmarks should include a brief note on the grid resolution and the observed order of accuracy at the final time to allow direct comparison with the 1D/2D convergence tables.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work, the clear summary of its contributions, and the recommendation for minor revision. No specific major comments were listed in the report.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation chain begins from the stated GQL reformulation of nonlinear RHD constraints into linear inequalities, followed by an explicit rational stereographic parameterization that algebraically converts the min-max problem into a generalized Rayleigh quotient solved by small symmetric eigenproblems. These steps are presented as direct equivalences without any fitted parameter being relabeled as a prediction, without self-definitional loops, and without load-bearing uniqueness imported via self-citation. The numerical benchmarks function only as post-derivation validation. No enumerated circular pattern is exhibited by the quoted construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides insufficient detail to identify any free parameters, axioms, or invented entities; full text required for audit.

pith-pipeline@v0.9.1-grok · 5820 in / 1170 out tokens · 71804 ms · 2026-07-01T03:34:40.585080+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

63 extracted references · 3 canonical work pages

  1. [1]

    Abgrall, M

    R. Abgrall, M. Jiao, Y . Liu, and K. Wu. A novel and simple invariant-domain-preserving framework for PAMPA scheme: 1D case. SIAM J. Sci. Comput., 47(6):A3536–A3565, 2025

    2025

  2. [2]

    Abgrall, M

    R. Abgrall, M. Jiao, Y . Liu, and K. Wu. Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) scheme: one-dimensional case. Commun. Comput. Phys., 39(5):29–58, 2026

    2026

  3. [3]

    Abgrall, Y

    R. Abgrall, Y . Liu, and W. Boscheri. Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) on polyg- onal meshes. 2025. https: //arxiv.org/abs/2502.10069

    work page arXiv 2025

  4. [4]

    D. L. Book, J. P. Boris, and K. Hain. Flux-corrected transport II: Generalizations of the method.J. Comput. Phys., 18(3):248– 283, 1975

    1975

  5. [5]

    Borges, M

    R. Borges, M. Carmona, B. Costa, and W. S. Don. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. Journal of computational physics, 227(6):3191–3211, 2008

    2008

  6. [6]

    J. P. Boris and D. L. Book. Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works. J. Comput. Phys., 11(1):38–69, 1973

    1973

  7. [7]

    J. P. Boris and D. L. Book. Flux-corrected transport. III. Minimal-error FCT algorithms. J. Comput. Phys., 20(4):397–431, 1976

    1976

  8. [8]

    C. Cai, J. Qiu, and K. Wu. Provably convergent newton–raphson methods for recovering primitive variables with applications to physical-constraint-preserving hermite weno schemes for relativistic hydrodynamics.J. Comput. Phys., 498:112669, 2024

    2024

  9. [9]

    H. Cao, M. Peng, and K. Wu. Robust discontinuous Galerkin methods maintaining physical constraints for general relativistic hydrodynamics. J. Comput. Phys., 526:113770, 2025

    2025

  10. [10]

    Y . Chen, Y . Kuang, and H. Tang. Second-order accurate BGK schemes for the special relativistic hydrodynamics with the Synge equation of state. J. Comput. Phys., 442:110438, 2021

    2021

  11. [11]

    Chen and K

    Y . Chen and K. Wu. A physical-constraint-preserving finite volume WENO method for special relativistic hydrodynamics on unstructured meshes. J. Comput. Phys., 466:111398, 2022

    2022

  12. [12]

    S. Cui, K. Wu, and L. Xu. On local minimum entropy principle of high-order schemes for relativistic Euler equations. Math. Comp., in press, 2025. https://doi.org/10.1090/mcom/4139

    work page doi:10.1090/mcom/4139 2025

  13. [13]

    Duan and H

    J. Duan and H. Tang. High-order accurate entropy stable finite di fference schemes for one-and two-dimensional special relativistic hydrodynamics. Adv. Appl. Math. Mech., 12(1):1–29, 2020

    2020

  14. [14]

    Duan and H

    J. Duan and H. Tang. Entropy stable adaptive moving mesh schemes for 2D and 3D special relativistic hydrodynamics. J. Comput. Phys., 426:109949, 2021

    2021

  15. [15]

    Guermond, M

    J.-L. Guermond, M. Nazarov, B. Popov, and I. Tomas. Second-order invariant domain preserving approximation of the Euler equations using convex limiting. SIAM J. Sci. Comput., 40(5):A3211–A3239, 2018

    2018

  16. [16]

    Guermond and B

    J.-L. Guermond and B. Popov. Invariant domains and second-order continuous finite element approximation for scalar conservation equations. SIAM J. Numer. Anal., 55(6):3120–3146, 2017

    2017

  17. [17]

    Guermond, B

    J.-L. Guermond, B. Popov, and I. Tomas. Invariant domain preserving discretization-independent schemes and convex limiting for hyperbolic systems. Comput. Methods Appl. Mech. Eng., 347:143–175, 2019

    2019

  18. [18]

    He and H

    P. He and H. Tang. An adaptive moving mesh method for two-dimensional relativistic hydrodynamics. Commun. Comput. Phys., 11(1):114–146, 2012

    2012

  19. [19]

    X. Y . Hu, N. A. Adams, and C.-W. Shu. Positivity-preserving method for high-order conservative schemes solving com- pressible Euler equations. J. Comput. Phys., 242:169–180, 2013

    2013

  20. [20]

    P. A. Hughes, M. A. Miller, and G. C. Duncan. Three-dimensional hydrodynamic simulations of relativistic extragalactic jets. Astrophys. J., 572(2):713, 2002

    2002

  21. [21]

    Jiang and C.-W

    G.-S. Jiang and C.-W. Shu. E fficient implementation of weighted ENO schemes. J. Comput. Phys., 126(1):202–228, 1996. 42

    1996

  22. [22]

    Jiang and Z

    Y . Jiang and Z. Xu. Parametrized maximum principle preserving limiter for finite di fference WENO schemes solving convection-dominated diffusion equations. SIAM J. Sci. Comput., 35(6):A2524–A2553, 2013

    2013

  23. [23]

    D. Kuzmin. Positive finite element schemes based on the flux-corrected transport procedure. Computational Fluid and Solid Mechanics, Elsevier, pages 887–888, 2001

    2001

  24. [24]

    D. Kuzmin. Monolithic convex limiting for continuous finite element discretizations of hyperbolic conservation laws. Com- put. Methods Appl. Mech. Engrg., 361:112804, 2020

    2020

  25. [25]

    Kuzmin, R

    D. Kuzmin, R. L ¨ohner, and S. Turek. Flux-corrected transport: principles, algorithms, and applications . Springer Science & Business Media, 2012

    2012

  26. [26]

    Kuzmin, M

    D. Kuzmin, M. M ¨oller, and S. Turek. High-resolution FEM–FCT schemes for multidimensional conservation laws. Comput. Methods Appl. Mech. Engrg., 193(45-47):4915–4946, 2004

    2004

  27. [27]

    Liang and Z

    C. Liang and Z. Xu. Parametrized maximum principle preserving flux limiters for high order schemes solving multi- dimensional scalar hyperbolic conservation laws. J. Sci. Comput., 58(1):41–60, 2014

    2014

  28. [28]

    X.-D. Liu, S. Osher, and T. Chan. Weighted essentially non-oscillatory schemes. J. Comput. Phys., 115(1):200–212, 1994

    1994

  29. [29]

    Marquina, S

    A. Marquina, S. Serna, and J. M. Ib ´a˜nez. Capturing composite waves in non-convex special relativistic hydrodynamics. J. Sci. Comput., 81(3):2132–2161, 2019

    2019

  30. [30]

    J. M. Mart ´ı and E. M¨uller. Numerical hydrodynamics in special relativity. Living Rev. Relativ., 6(1):7, 2003

    2003

  31. [31]

    J. M. Mart ´ı and E. M ¨uller. Grid-based methods in relativistic hydrodynamics and magnetohydrodynamics. Living Rev. Comput. Astrophys., 1(1):3, 2015

    2015

  32. [32]

    M. M. May and R. H. White. Hydrodynamic calculations of general-relativistic collapse. Phys. Rev., 141(4):1232, 1966

    1966

  33. [33]

    M. Peng, Z. Sun, and K. Wu. OEDG: Oscillation-eliminating discontinuous Galerkin method for hyperbolic conservation laws. Mathematics of Computation, 94(353):1147–1198, 2025

    2025

  34. [34]

    Qin, C.-W

    T. Qin, C.-W. Shu, and Y . Yang. Bound-preserving discontinuous Galerkin methods for relativistic hydrodynamics. J. Comput. Phys., 315:323–347, 2016

    2016

  35. [35]

    Radice and L

    D. Radice and L. Rezzolla. Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and appli- cation to spherically symmetric spacetimes. Phys. Rev. D, 84(2):024010, 2011

    2011

  36. [36]

    Radice and L

    D. Radice and L. Rezzolla. THC: a new high-order finite-di fference high-resolution shock-capturing code for special- relativistic hydrodynamics. Astron. Astrophys., 547:A26, 2012

    2012

  37. [37]

    Rezzolla and O

    L. Rezzolla and O. Zanotti. Relativistic Hydrodynamics. Oxford University Press, 2013

    2013

  38. [38]

    Suresh and H

    A. Suresh and H. T. Huynh. Accurate monotonicity-preserving schemes with Runge–Kutta time stepping. J. Comput. Phys., 136(1):83–99, 1997

    1997

  39. [39]

    Wang and W

    B.-S. Wang and W. S. Don. A ffine-invariant WENO weights and operator. Applied Numerical Mathematics, 181:630–646, 2022

    2022

  40. [40]

    J. R. Wilson. Numerical study of fluid flow in a Kerr space. Astrophys. J., 173:431, 1972

    1972

  41. [41]

    K. Wu. Design of provably physical-constraint-preserving methods for general relativistic hydrodynamics. Phys. Rev. D, 95(10):103001, 2017

    2017

  42. [42]

    K. Wu. Positivity-preserving analysis of numerical schemes for ideal magnetohydrodynamics. SIAM J. Numer. Anal. , 56(4):2124–2147, 2018

    2018

  43. [43]

    K. Wu. Minimum principle on specific entropy and high-order accurate invariant region preserving numerical methods for relativistic hydrodynamics. SIAM J. Sci. Comput., 43(6):B1164–B1197, 2021

    2021

  44. [44]

    Wu and C.-W

    K. Wu and C.-W. Shu. A provably positive discontinuous Galerkin method for multidimensional ideal magnetohydrodynam- ics. SIAM J. Sci. Comput., 40(5):B1302–B1329, 2018

    2018

  45. [45]

    Wu and C.-W

    K. Wu and C.-W. Shu. Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes. Numer. Math., 142(4):995–1047, Aug 2019

    2019

  46. [46]

    Wu and C.-W

    K. Wu and C.-W. Shu. Entropy symmetrization and high-order accurate entropy stable numerical schemes for relativistic MHD equations. SIAM J. Sci. Comput., 42(4):A2230–A2261, 2020

    2020

  47. [47]

    Wu and C.-W

    K. Wu and C.-W. Shu. Provably physical-constraint-preserving discontinuous Galerkin methods for multidimensional rela- tivistic MHD equations. Numer. Math., 2021

    2021

  48. [48]

    Wu and C.-W

    K. Wu and C.-W. Shu. Geometric quasilinearization framework for analysis and design of bound-preserving schemes. SIAM Rev., 65(4):1031–1073, 2023

    2023

  49. [49]

    Wu and H

    K. Wu and H. Tang. High-order accurate physical-constraints-preserving finite di fference WENO schemes for special rela- tivistic hydrodynamics. J. Comput. Phys., 298:539–564, 2015

    2015

  50. [50]

    Wu and H

    K. Wu and H. Tang. Admissible states and physical-constraints-preserving schemes for relativistic magnetohydrodynamic equations. Math. Models Methods Appl. Sci., 27(10):1871–1928, 2017

    1928

  51. [51]

    Wu and H

    K. Wu and H. Tang. Physical-constraint-preserving central discontinuous Galerkin methods for special relativistic hydrody- namics with a general equation of state. Astrophys. J. Suppl. Ser., 228(1), 2017

    2017

  52. [52]

    K. Wu, X. Zhang, and C.-W. Shu. High order numerical methods preserving invariant domain for hyperbolic and related systems. SIAM Rev., in press, 2025. https://doi.org/10.48550/arXiv.2512.09116

    work page doi:10.48550/arxiv.2512.09116 2025

  53. [53]

    Xiong, J.-M

    T. Xiong, J.-M. Qiu, and Z. Xu. Parametrized positivity preserving flux limiters for the high order finite di fference WENO 43 scheme solving compressible Euler equations. J. Sci. Comput., 67(3):1066–1088, 2016

    2016

  54. [54]

    L. Xu, S. Ding, and K. Wu. High-order accurate entropy stable schemes for relativistic hydrodynamics with general synge- type equation of state. J. Sci. Comput., 98(2):43, 2024

    2024

  55. [55]

    Z. Xu. Parametrized maximum principle preserving flux limiters for high order schemes solving hyperbolic conservation laws: one-dimensional scalar problem. Math. Comp., 83(289):2213–2238, 2014

    2014

  56. [56]

    S. T. Zalesak. Fully multidimensional flux-corrected transport algorithms for fluids. J. Comput. Phys., 31(3):335–362, 1979

    1979

  57. [57]

    S. T. Zalesak. The design of Flux-Corrected Transport (FCT) algorithms for structured grids. Springer, 2012

    2012

  58. [58]

    Zhang and A

    W. Zhang and A. I. MacFadyen. RAM: A relativistic adaptive mesh refinement hydrodynamics code. Astrophys. J. Suppl. S., 164(1):255, 2006

    2006

  59. [59]

    X. Zhang. On positivity-preserving high order discontinuous Galerkin schemes for compressible Navier-Stokes equations. J. Comput. Phys., 328:301–343, 2017

    2017

  60. [60]

    Zhang and C.-W

    X. Zhang and C.-W. Shu. On maximum-principle-satisfying high order schemes for scalar conservation laws. J. Comput. Phys., 229(9):3091–3120, 2010

    2010

  61. [61]

    Zhang and C.-W

    X. Zhang and C.-W. Shu. On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes. J. Comput. Phys., 229(23):8918–8934, 2010

    2010

  62. [62]

    Zhang and C.-W

    X. Zhang and C.-W. Shu. Positivity-preserving high order finite difference WENO schemes for compressible Euler equations. J. Comput. Phys., 231(5):2245–2258, 2012

    2012

  63. [63]

    Zhao and H

    J. Zhao and H. Tang. Runge–Kutta discontinuous Galerkin methods with WENO limiter for the special relativistic hydrody- namics. J. Comput. Phys., 242:138–168, 2013

    2013