【请教】气液两相流动+udf遇到的问题

andyecs

计算一个气液两相流动，不涉及换热，且流动过程中不考虑相变。采用mixture模型，速度进口（20m/S）,压力出口（500kpa），进口液相体积分数为0.1，出口液相体积分数未设置

物性中气、液相粘度为恒定值，液相密度为一定值，气相密度用udf定义，如下

//*****************************************//
//********vapor density********************//
//*****************************************//
#include "udf.h"
DEFINE_PROPERTY(cell_density,cell,thread)
{
real v_den;   //气相密度//
real v_p=C_P(cell,thread);  
v_den=24.62*v_p/500000;  //出口500kpa时气相密度为24.62，其它计算点气相密度近似可以认为和计算点的压力成正比//
return v_den;
}

用进口初始化，但每次算到第四步迭代就报错，不知道哪里出了问题。诚心请教，谢谢！

andyecs

  没人做过类似的算例么？

h04581803

density (as a function of pressure) for compressible liquid flows only (Mixture, cavitation models )
我在DEFINE-PROPERTY帮助看到的，希望可以帮助你，你做的问题应该是不可压缩流吧，那就不能这么定义了；

andyecs

感谢楼上。

问题解决了，调低松弛因子，可以得到收敛结果