| 查看: 1392 | 回复: 5 | ||||
| 当前主题已经存档。 | ||||
[交流]
【资源】水蒸发udf源文件 已有1人参与
|
||||
|
水蒸发udf源文件 #include "udf.h" #include "sg_mphase.h" #define T_SAT 373 #define LAT_HT 2256254 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); 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.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_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,mix_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,mix_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; } [ Last edited by kuhailangyu on 2009-11-14 at 21:55 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 杂记 |
» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有3人回复
-
关于代码变化问题，想知道的进来
已经有16人回复
-
fileCode有新解读？
已经有18人回复
-
这样的filecode谁见过
已经有15人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2009-11-14 20:35:01
3楼2010-01-18 10:23:51
4楼2010-01-18 12:55:38
5楼2010-01-27 15:22:13
6楼2010-04-15 10:02:20
