compressiblemultiphaseinterfoam计算的bug

[mohui]

我发现不论是compressiblemultiphaseinterfoam又或是compressibleinterfoam，计算类似于throttle带缩口问题，算到后面都会出错，原因是回流区的压力过低，虽然用pmin限制强行每次算完的p都大于0，但是不知道为什么气相密度会小于0，最后出错。所以我不知道如何才能避免这种问题，又或者本身求解器就是有问题的，无法解决。0_1504614071145_2017-09-05_20-19-18.png

[mohui]

@李东岳 ，大神能给点意见吗？

[李东岳]

很难判断... :crying:

[mohui]

@李东岳 好吧，这个真的算是可压缩求解器的bug了，确实算会崩掉，但是不知道能不能通过调整什么能够避免这个问题。

[mohui]

@李东岳 东岳老师，我知道为什么会出现负密度了，原因是mixture.correctRho(p_rgh - p_rgh_0);，由于算缩口问题，迭代过程dp的负值太大了，导致密度就出现负值了，话说回来这个dp到底是什么原理呢，为什么还需要更新一次密度，如果直接按照状态方程求出每个相的rho不就行了吗？

[yuan_neu]

修改一些mixture的类，闲置rho的最小值，应该可以跳过这个错误吧

[mohui]

@yuan_neu 能具体说下限制rho的最小值具体是什么意思吗？现在有个很大的问题是真的不知道如何去修改液体的状态方程。或者说更新rho密度是不完善的，会出现不合理的解。

[youmengtian]

这个问题使用fluent算貌似同样会出现这个问题，其实这时候已经出现了部分空化，处理不当产生负压。