那位高手做过把simple算法由不可压推广到可压？

maximus

俺正在做可压，但原来不可压采用的simple算法，因此采用推广到可压的simple算法来计算可压，方成离散采用FVM（计算传热学里常用的控制容积积分法）。但现在碰到了几个问题，感觉算出来的温度场不对（速度场也不对了，但至少形状已经很好了），现在不知道是不是由于对源项里散度的离散有问题还是其他原因。
哪位高手做过的话，能否指点一下。很着急！！
谢谢