在帮助文档给出的RK方程UDRGM案例基础上对其中声速计算的更改？

mengdaolemeng

本人想要计算含杂质超临界CO2减压波速度，想用fluent对激波管进行模拟，计算过程总是一开始就发散，觉得可能是模拟的泄放过程相态要经过两相区，而软件内自带的真实气体状态方程不做修改的话都无法做到，因此开始尝试自己编UDF，想在官方案例的基础上进行修改，使得两相区内声速的计算改为如图片所示；请问如何在udf中如何实现？