我们使用基于二次冲击近似的Riemann求解器，在大冲击波的存在下保持稳定，在米尔坐标中构造了新的Godunov型相对论流体力学代码。通过比较各种问题的数值计算和解析解，例如冲击管，物质在真空中的膨胀，Landau–Khalatnikov解以及在Bjorken流和Gubser流附近的波动传播，我们检查了数值算法的正确性。我们在具有高能重离子碰撞的初始条件的纵向流体动力膨胀的测试问题中研究代码的能量和动量守恒性质。我们还将讨论物质在真空中膨胀和保持特性的测试问题中的数值粘度。此外，我们讨论了Milne坐标中相对论性数值流体动力学的源项如何影响数值稳定性。