| 查看: 977 | 回复: 1 | ||
zhao19861011银虫 (正式写手)
|
[求助]
fluent udf的相变问题，还请大神指导一下！！！
|
|
最近测试了一下udf发现气态源项的部分有可能存在不妥，我觉得是气体温度低于临界温度，气态转化为液态，由于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零；但是气体温度大于临界温度时，我就不知道怎么描述了，还请指点一下 以下为具体的udf #include "udf.h" #include "sg_mphase.h" DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) { Thread *mix_th, *pri_th; real m_dot_v; real T_SAT; mix_th =THREAD_SUPER_THREAD(sec_th); pri_th =THREAD_SUB_THREAD(mix_th,0); T_SAT=239.8+0.0000092568*C_P(cell, mix_th); if(C_T(cell, mix_th)<=T_SAT) { m_dot_v=10.0*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th)-T_SAT)/T_SAT; dS[eqn]=0.0; } else { m_dot_v=-10.0*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell, mix_th))/T_SAT; dS[eqn]=-10.0*C_R(cell, sec_th)*fabs(C_T(cell, mix_th)-T_SAT)/T_SAT; } return m_dot_v; } |
回复此楼» 猜你喜欢
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
2楼2017-05-22 09:48:25
