| 查看: 985 | 回复: 1 | ||||
w420049436新虫 (小有名气)
|
[求助]
做水冷凝结果不对，请教大神UDF哪里出错了？ 已有1人参与
|
|
我现在正在学习水冷凝的UDF，然后上网找了例子模拟，发现结果出错，但是不知道错在哪里。恳求各位大神指正 模拟的参数如下：水蒸汽以400度进入管道，管道壁面温度为283度，水蒸汽的速度为0.1m/s。管径为0.05m，长度为3m。 #include "udf.h" #include "sg_mphase.h" #define T_SAT 373 #define LAT_HT 2.257e6 DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_v; 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_v = 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_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.; } return m_dot_v; } DEFINE_SOURCE(lip_src, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_l; 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_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.; } 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.1*C_R(cell, sec_th)* fabs(C_T(cell, sec_th) - T_SAT)/T_SAT; } return m_dot_l; } 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; }  这个是计算66s后的气相体积分数，根本看不到有水滴行程  这个是计算66s后的温度 |
回复此楼» 收录本帖的淘帖专辑推荐
 博士科研新路历程 |
» 猜你喜欢
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一下大神们关于fluent中UDF的问题，谢谢！
已经有3人回复
- 利用DEFINE_EXCHANGE_PROPERTY编的蒸发和冷凝的udf，大侠们帮我看看对不对
已经有9人回复
- 有关水蒸汽冷凝的UDF
已经有8人回复
- 【求助】冷凝的UDF程序
已经有15人回复
xiaobingwill
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 47
- 红花: 1
- 帖子: 25
- 在线: 7.1小时
- 虫号: 1866037
- 注册: 2012-06-21
- 性别: GG
- 专业: 可再生与替代能源利用中的
2楼2014-10-14 01:30:03
10