@guanghuilineu 那看起来发散是个很正常的结果，如果相变区速度不为0时，也变速也变密度。
如果只变速，和周围空间交换能量，即动能转换为其他形式的能量，没问题。
如果速度和密度一起变，动能是依赖于质量的。这时候，质量也变了，不守恒了，就崩了。
消失的那部分质量，带有的动能如何传递出去。
压力就是传递能量的桥梁。在传递的时候一头突变了，就断了。可不可以先使相变区域停止流动再变密度。如果这样，不会出现动能突变的状况。或者在变密度操作后立刻更新一下速度；变前计算动能；变后用动能，根据变后的质量，计算新的速度。
或者最粗鲁的方法，用极小的时间步，这样突变也可以稳住。
