| 查看: 1125 | 回复: 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; } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
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
