| 查看: 948 | 回复: 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; } |
回复此楼» 猜你喜欢
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有27人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
-
什么时候出结果，有咨询渠道？？？
已经有3人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
2楼2017-05-22 09:48:25
