Bound Preserving Lax-Wendroff Flux Reconstruction Method for Special Relativistic Hydrodynamics

Lax-Wendroff flux reconstruction (LWFR) schemes have high order of accuracy in both space and time despite having a single internal time step. Here, we design a Jacobian-free LWFR type scheme to solve the special relativistic hydrodynamics equations on Cartesian grids. We then blend the scheme with a first-order finite volume scheme to control the oscillations near discontinuities. We also use a scaling limiter to preserve the physical admissibility of the solution after ensuring the scheme is admissible in means. A particular focus is given to designing a discontinuity indicator model to detect the local non-smoothness in the solution of the highly non-linear relativistic hydrodynamics equations. Finally, we present numerical results for a wide range of test cases to show the robustness and efficiency of the proposed scheme.