| 查看: 978 | 回复: 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; } |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有5人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有3人回复
-
售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人回复
2楼2017-05-22 09:48:25
