| 查看: 1002 | 回复: 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; } |
回复此楼» 猜你喜欢
294求调剂
已经有6人回复
-
322求调剂：一志愿湖南大学 材料与化工（085600），已过六级。
已经有14人回复
-
282求调剂
已经有6人回复
-
调剂
已经有9人回复
-
315求调剂
已经有13人回复
-
材料调剂
已经有8人回复
-
070300化学学硕311分求调剂
已经有6人回复
-
求调剂
已经有7人回复
-
材料调剂
已经有9人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有9人回复
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
