| 查看: 1369 | 回复: 2 | |||
[求助]
相变udf求助
|
|
最近在做相变的模拟，udf使用的是论坛上大家通用的，我只是改了一下指针，结果就老是初始化报错，把udf放在下面，希望大家给看一下啊，谢谢各位，我个人感觉是温度的指针出了问题。 DEFINE_SOURCE(liq_src, cell, thread, dS, eqn) { /*Thread *mix_th, *sec_th;*/ Thread *mix_thread, *thread_gas, *thread_liq; real m_dot_l; real T_SAT; /*定义临界温度*/ real tem_pressure; /*定义单元压力*/ real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix; mix_thread= THREAD_SUPER_THREAD(thread); thread_gas= THREAD_SUB_THREAD(mix_thread, 1); thread_liq= THREAD_SUB_THREAD(mix_thread, 0); vof_gas=C_VOF(cell, thread_gas); vof_liq=C_VOF(cell, thread_liq); roul_gas=C_R(cell, thread_gas); roul_liq=C_R(cell, thread_liq); t_gas=C_T(cell, thread_gas); t_liq=C_T(cell, thread_liq); t_mix=C_T(cell, mix_thread); /*mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1);*/ tem_pressure=(C_P(cell,mix_thread)+101325)/1000000; /*读取单元压力*/ T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47; /*定义饱和温度与压力的关系*/ if(t_liq>=T_SAT) { m_dot_l = -0.1*vof_liq*roul_liq*fabs(t_liq - T_SAT)/T_SAT; dS[eqn] = -0.1*roul_liq*fabs(t_liq - T_SAT)/T_SAT; } { m_dot_l = 0.1*vof_gas*roul_gas*fabs(t_liq - T_SAT)/T_SAT; dS[eqn] =0; } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, thread, dS, eqn) { Thread *mix_thread, *thread_gas, *thread_liq; real m_dot_v; real T_SAT; /*定义临界温度*/ real tem_pressure; /*定义单元压力*/ real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix; mix_thread= THREAD_SUPER_THREAD(thread); thread_gas= THREAD_SUB_THREAD(mix_thread, 1); thread_liq= THREAD_SUB_THREAD(mix_thread, 0); vof_gas=C_VOF(cell, thread_gas); vof_liq=C_VOF(cell, thread_liq); roul_gas=C_R(cell, thread_gas); roul_liq=C_R(cell, thread_liq); t_gas=C_T(cell, thread_gas); t_liq=C_T(cell, thread_liq); t_mix=C_T(cell, thread); tem_pressure=(C_P(cell,mix_thread)+101325)/1000000; /*读取单元压力*/ T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47; /*定义饱和温度与压力的关系*/ if(t_gas>=T_SAT) { m_dot_v = 0.1*vof_liq*roul_liq* fabs(t_gas - T_SAT)/T_SAT; dS[eqn] = 0; /*1*/ } else { m_dot_v = -0.1*vof_liq*roul_liq*fabs(T_SAT-t_gas)/T_SAT; dS[eqn] = -0.1*roul_liq*fabs(t_gas - T_SAT)/T_SAT; } return m_dot_v; } DEFINE_SOURCE(enrg_src, cell, thread, dS, eqn) { Thread *mix_thread, *thread_gas, *thread_liq; real m_dot_l; real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix; real T_SAT; /*定义临界温度*/ real tem_pressure; /*定义单元压力*/ mix_thread= THREAD_SUPER_THREAD(thread); thread_gas= THREAD_SUB_THREAD(mix_thread, 1); thread_liq= THREAD_SUB_THREAD(mix_thread, 0); vof_gas=C_VOF(cell, thread_gas); vof_liq=C_VOF(cell, thread_liq); roul_gas=C_R(cell, thread_gas); roul_liq=C_R(cell, thread_liq); t_gas=C_R(cell, thread_gas); t_liq=C_R(cell, thread_liq); t_mix=C_T(cell, thread); tem_pressure=(C_P(cell,mix_thread)+101325)/1000000; /*读取单元压力*/ T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47; /*定义饱和温度与压力的关系*/ if(t_mix>=T_SAT) { m_dot_l = -0.1*vof_liq*roul_liq* fabs(t_mix - T_SAT)/T_SAT; dS[eqn] =-0.1*vof_liq*roul_liq/T_SAT; } else { m_dot_l = 0.1*vof_gas*roul_gas* fabs(t_mix - T_SAT)/T_SAT; dS[eqn] =0.1*vof_gas*roul_gas/T_SAT; /*1*/ } return LAT_HT*m_dot_l; } |
回复此楼» 收录本帖的淘帖专辑推荐
 博士科研新路历程 |
» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有17人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神帮检查下一个简单的UDF
已经有13人回复
- 有关fluent UDF的问题
已经有15人回复
- 关于温度的传递问题，能否有人帮我看一下这个UDF是哪里出了错误呢？
已经有10人回复
- Ansys Fluent 编译UDF及加载UDF出现错误
已经有14人回复
- 急，在FLUENT中如何同时调用两个UDF？
已经有11人回复
- Fluent中DPM模型的UDF中使用P_POS()指针，然后给我提示undeclared variable 怎么办?
已经有6人回复
- 请问谁做过气液相变的fluent模拟呀
已经有3人回复
- fluent中相变过程怎样设置固液密度
已经有10人回复
- fluent UDF 修改物性
已经有13人回复
- UDF中如何编写体积分数的梯度
已经有9人回复
- fluent中UDF编译的问题！
已经有6人回复
- 关于水蒸发udf
已经有19人回复
- 相变模型 UDF 求助
已经有9人回复
- Fluent中UDF的编译问题
已经有25人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- UDF应用之后的初始化错误
已经有6人回复
|
本帖内容被屏蔽 |
2楼2020-03-12 19:06:17
zhangcx1208
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 6889.3
- 散金: 51
- 红花: 5
- 帖子: 309
- 在线: 112.4小时
- 虫号: 468285
- 注册: 2007-11-27
- 专业: 流体力学
3楼2020-12-26 11:19:55
28