Numerical Solution to the Relativistic Shock Tube Problem

This website provides numerical solutions of the relativistic shock tube problem. The user can specify various parameters and choose between different numerical methods. The relativistic hydro code behind, was written by Moritz Fischer under supervision of Ahmad Hujeirat. This project emerged from his lecture on numerical techniques for modeling relativistic hydrodynamics in summer semester of 2017. The developed code solves the relativistic Euler equations and is based on the finite volume discretization. There are two schemes implemented: time-explicit using forward Euler and time-implicit using the Crank-Nicolson method to ensure second order temporal accuracy. Moreover flux limiters are implemented in the transport operators to achieve second order spatial accuracy.

Movies that demonstrate the accuracy of the time-implicit scheme can be found here.

To compute a shock profile specify the following parameters and push the "Start simulation" button. After waiting until the simulation has finished, you can take a look at the results.

Do you want to chose more parameters? Click here.