Improving accurate numerical modeling for relativistic simulations is a critical area of research. A key advancement in this field was the adaptation of modern high-resolution shock-capturing (HRSC) schemes, originally designed for non-relativistic flows, to incorporate relativistic effects. However, simply extending these Newtonian HRSC schemes to relativistic
shock-capturing has proven challenging. Several studies have shown that achieving grid-convergent numerical solutions often requires impractically fine grid resolutions. This persistent challenge demands further in-depth
investigation. In this presentation, we’ll share our latest progress in enhancing shock-capturing numerical methods for relativistic flows. We’ll start by highlighting common numerical issues encountered when solving 1D shock-tube problems using existing HRSC approaches, especially when relativistic tangential velocities are present. To address these deficiencies, we’ll introduce our new hybrid method. This approach significantly improves the simulation of relativistic shock-tube solutions, delivering sharply resolved shock and contact waves.