超声波雾化vof to dpm残差值不断变大

老师，做非稳态超声波雾化，气体为压缩空气考虑了气体的可压缩性。先算了流场残差正常，加水后残差值就不断变大，之前计算不断变大后出站粘性比受限，紧接着是温度受限，最后库朗数过大发散。 调整过松弛因子，治标不治本。 采用了网格自适应，算法用的simple. 打开了warped face geadient correction，也只是延缓了发散。 是不是有哪些问题没注意到？