| 查看: 1334 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
------------------------------------------- 已有4人参与
|
||||
|
----------------------------- [ Last edited by lovewb on 2010-7-2 at 21:14 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 fluent仿真问题。 |
» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
sosososo1114
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2377.2
- 散金: 96
- 红花: 4
- 帖子: 488
- 在线: 57.2小时
- 虫号: 939993
- 注册: 2010-01-10
- 专业: 红外物理与技术
2楼2010-05-08 15:43:42
sosososo1114
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2377.2
- 散金: 96
- 红花: 4
- 帖子: 488
- 在线: 57.2小时
- 虫号: 939993
- 注册: 2010-01-10
- 专业: 红外物理与技术
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+1):这个UDF本站也有哦，呵呵 2010-05-08 22:10:33
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+1):这个UDF本站也有哦，呵呵 2010-05-08 22:10:33
|
人品再爆发一下咯，本人前段时间也在做蒸发冷凝，所以从网上搜了点资料，这是煮开水的UDF，特此声明：作者不是我，这个UDF是转载的，但是绝对能用！ #include "udf.h" #include "sg_mphase.h" #define T_SAT 373 #define LAT_HT 3.e3 real a; DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_l; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); a=C_T(cell, mix_th); 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, pri_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_R(cell, pri_th)* fabs(C_T(cell, pri_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; } DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_v; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); if(C_T(cell, mix_th)>=T_SAT){ m_dot_v = 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.; } else { m_dot_v = -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.1*C_R(cell, sec_th)* fabs(C_T(cell, sec_th) - T_SAT)/T_SAT; } return m_dot_v; } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) { Thread *pri_th, *sec_th; real m_dot; pri_th = THREAD_SUB_THREAD(mix_th, 0); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)>=T_SAT){ m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; } else { m_dot = 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.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;} return LAT_HT*m_dot; } |
3楼2010-05-08 15:56:46
sosososo1114
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2377.2
- 散金: 96
- 红花: 4
- 帖子: 488
- 在线: 57.2小时
- 虫号: 939993
- 注册: 2010-01-10
- 专业: 红外物理与技术
9楼2010-05-09 13:23:20