向封闭腔充入可压缩油液的入口边界条件只能是设定为质量流速巴,我把其设成5kg/s才能达到8-10MPa左右的压力,我的密度udf是这样的
#include"udf.h"
#define modulus 7.0e8
#define row_ref 889.0
#define p_ref 101325
DEFINE_PROPERTY(vardensi,c,t)
{
real row;
real p,dp;
real p_operating;
p_operating=RP_Get_Real("operating-pressure");
p=C_P(c,t)+p_operating;
dp=p-p_ref;
row=row_ref*(dp+modulus)/modulus;
return row;
}
DEFINE_PROPERTY(soundspeed,c,t)
{
real a;
real p,dp;
real p_operating;
p_operating=RP_Get_Real("operating-pressure");
p=C_P(c,t)+p_operating;
dp=p-p_ref;
a=modulus*sqrt(1/(row_ref*(dp+modulus)));
return a;
}
现在我怀疑是p的问题,相关的C_P(c,t),p_operating是不是用在这里不妥阿? |