| 查看: 1003 | 回复: 2 | |||
hjnature铜虫 (初入文坛)
|
[求助]
在源项中调用温度梯度 已有2人参与
|
|
想在两相流源项的计算中使用温度梯度，但是为什么udf能顾编译通过，但是算了一步之后计算就报错呢，源程序如下，我是fluent新手，希望各位大神解答， #include "udf.h" #include "sg_mphase.h" #define T_SAT 373 #define L 230000.0 #define d 1000 #define k_L 0.6791 DEFINE_SOURCE(liq_src, cell, thread, dS, eqn) { real void_s, tg,m_dot_l=0; void_s=C_VOF(cell,thread); if (NULL != THREAD_STORAGE(thread,SV_T_G)) { tg=C_T_G(cell,thread)[0]; if(C_T(cell,thread)>=T_SAT&&void_s>0&&void_s<1) { m_dot_l =-k_L*tg*d/L; dS[eqn] =0; } else if(C_T(cell,thread)<T_SAT&&void_s>0&&void_s<1) { m_dot_l =k_L*tg*d/L; dS[eqn] =0.; } else{ m_dot_l =0; } } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, thread, dS, eqn) { real void_s, tg,m_dot_v=0; void_s=C_VOF(cell, thread); if (NULL != THREAD_STORAGE(thread,SV_T_G)) { tg=C_T_G(cell,thread)[0]; if(C_T(cell, thread)>=T_SAT&&void_s>0&&void_s<1) { m_dot_v =k_L*tg*d/L; dS[eqn] = 0.; } else if(C_T(cell, thread)<T_SAT&&void_s>0&&void_s<1) { m_dot_v = -k_L*tg*d/L; dS[eqn] = 0; } else { m_dot_v =0; } } return m_dot_v; } |
回复此楼» 猜你喜欢
复试调剂
已经有5人回复
-
一志愿北京化工085600 310分求调剂
已经有11人回复
-
一志愿北京2，材料与化工308求调剂
已经有10人回复
-
0817化学工程与技术求调剂，一志愿中海洋319
已经有9人回复
-
一志愿电子科技大学085600材料与化工 329分求调剂
已经有6人回复
-
304求调剂（085602，过四级，一志愿985）
已经有7人回复
-
求调剂
已经有3人回复
-
一志愿北京化工大学，初试成绩350求调剂
已经有14人回复
-
323分（计算机视觉和大模型项目）能直接上手
已经有3人回复
-
311分 22408 求调剂
已经有3人回复
2楼2017-07-02 22:48:43
mycc
专家顾问 (正式写手)
-
专家经验: +181 - 应助: 230 (大学生)
- 金币: 2589.6
- 红花: 35
- 帖子: 531
- 在线: 58.8小时
- 虫号: 881523
- 注册: 2009-10-23
- 性别: GG
- 专业: 流体力学
- 管辖: 仿真模拟
3楼2017-07-03 10:52:08
