compressiblemultiphaseinterfoam计算的bug
-
我发现不论是compressiblemultiphaseinterfoam又或是compressibleinterfoam,计算类似于throttle带缩口问题,算到后面都会出错,原因是回流区的压力过低,虽然用pmin限制强行每次算完的p都大于0,但是不知道为什么气相密度会小于0,最后出错。所以我不知道如何才能避免这种问题,又或者本身求解器就是有问题的,无法解决。0_1504614071145_2017-09-05_20-19-18.png
-
@李东岳 ,大神能给点意见吗?
-
很难判断... :crying:
-
@李东岳 好吧,这个真的算是可压缩求解器的bug了,确实算会崩掉,但是不知道能不能通过调整什么能够避免这个问题。
-
@李东岳 东岳老师,我知道为什么会出现负密度了,原因是mixture.correctRho(p_rgh - p_rgh_0);,由于算缩口问题,迭代过程dp的负值太大了,导致密度就出现负值了,话说回来这个dp到底是什么原理呢,为什么还需要更新一次密度,如果直接按照状态方程求出每个相的rho不就行了吗?
-
修改一些mixture的类,闲置rho的最小值,应该可以跳过这个错误吧
-
@yuan_neu 能具体说下限制rho的最小值具体是什么意思吗?现在有个很大的问题是真的不知道如何去修改液体的状态方程。或者说更新rho密度是不完善的,会出现不合理的解。
-
这个问题使用fluent算貌似同样会出现这个问题,其实这时候已经出现了部分空化,处理不当产生负压。