|
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
各位大虾,我做的是水遇高温气体蒸发成水蒸汽的问题,用的是mixture模型,要把水和水蒸汽之间的传质加入到质量源项,我在网上找到了一段有关水蒸发的udf,其中的温度项是用C_T(cell,mix_thread),但fluent 的udf 帮助中一个相似的例子中对应的温度项却是用C_T(cell,liq_thread)来表示的,请问各位路过的大虾,这个温度到底应该是用液相温度C_T(cell,liq_thread)表示,还是用混合相温度C_T(cell,mix_thread)?网上的udf如下:(请留意笑脸处)
DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相质量源项UDF
{
Thread *mix_th, *sec_th; ‘定义计算区线指针
real m_dot_l; ‘定义液相质量转移 kg/(m2.s)
mix_th = THREAD_SUPER_THREAD(pri_th); ‘指向混合区的主相即液相的指针
sec_th = THREAD_SUB_THREAD(mix_th, 1);指向单相控制区的气相的指针,气相为第二相
if(C_T(cell, mix_th) >=T_SAT) ‘如果液相单元的温度高于蒸发温度,液相向气相的质量质量转移
{
m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
‘ 定义源项对质量转移偏导 ,
}
else {
m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
fabs(T_SAT-C_T(cell,mix_th):) )/T_SAT;
‘如果指向混合区液相的单元温度小于蒸发温度,气相向液相的质量转移,液相得
dS[eqn] = 0.; ‘由于是气相向液相转移,所以液相的质量源项对质量转移的偏导为零
}
return m_dot_l;
}
fluent帮助中的udf:
if (C_T(cell, liq):) ) >= T_SAT)
{
m_lg = -0.1*C_VOF(cell,liq):) *C_R(cell,liq)*
fabs(C_T(cell,liq):) -T_SAT)/T_SAT;
[ 本帖最后由 h04581803 于 2009-3-16 07:41 编辑 ] |
评分
-
1
查看全部评分
-
|