| 查看: 1065 | 回复: 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; } |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有16人回复
-
帮忙看看fileCode
已经有9人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有4人回复
-
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
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
100