| 查看: 1088 | 回复: 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; } |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有12人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
mycc
专家顾问 (正式写手)
-
专家经验: +181 - 应助: 230 (大学生)
- 金币: 2589.6
- 红花: 35
- 帖子: 531
- 在线: 58.8小时
- 虫号: 881523
- 注册: 2009-10-23
- 性别: GG
- 专业: 流体力学
- 管辖: 仿真模拟
3楼2017-07-03 10:52:08
2楼2017-07-02 22:48:43
30