| 查看: 1904 | 回复: 4 | |||
[交流]
fluent问题急求大神帮助！！！ 已有4人参与
|
|
最近在模拟饱和蒸汽倾斜管内凝结换热的例子，可是液相体积分数一直是在10e-7左右徘徊，udf就是用的网上搜到的蒸发冷凝模型，另外还有个质量输运方程，具体见下方，然后我设置气相为主相液相为次相，用的mix模型，有大神能帮我解答一下吗，要是有热心人能帮帮忙模拟一下最好了，小弟感激不尽，因为是初学者，所以好多地方一头雾水，谢谢大家了！管子模型是0.2m×10m，2维模型，倾角60°，管内压力15kpa恒定，蒸汽温度327K，壁温293K，入口质量边界条件0.1kg/s，出口压力边界条件出口混合物温度325K #include "udf.h" #include "sg_mphase.h" #define T_SAT 327 #define LAT_HT 2.e3 DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_l; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th,1); if(C_T(cell, mix_th)>=T_SAT) { m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; } else { m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.; } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_v; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); if(C_T(cell, mix_th)>=T_SAT) { m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dS[eqn] = 0.; } else { m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; } return m_dot_v; } DEFINE_MASS_TRANSFER(lic_gas_source,cell,thread,from_index,from_species_index,to_index,to_species_index) { Thread *gas, *liq; real m_1g; gas=THREAD_SUB_THREAD(thread,from_index); liq=THREAD_SUB_THREAD(thread,to_index); m_1g=0.; if(C_T(cell,liq)>=T_SAT) m_1g=-0.1*C_VOF(cell,liq)*C_R(cell,liq)*fabs(C_T(cell,liq)-T_SAT)/T_SAT; else{ if(m_1g==0.) m_1g=0.1*C_VOF(cell,gas)*C_R (cell,gas)* fabs(T_SAT-C_T(cell,gas))/T_SAT; } return (m_1g); } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) { Thread *pri_th, *sec_th; real m_dot; pri_th = THREAD_SUB_THREAD(mix_th, 0); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)>=T_SAT) { m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; } else { m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT; } return LAT_HT*m_dot; } [ Last edited by li2hebe on 2014-2-2 at 17:00 ] |
回复此楼» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 崩溃了 求大神指导
已经有6人回复
- fluent封闭流场的遇到的一些问题，麻烦大家来探讨一下，求大神指点。多谢各位
已经有9人回复
- 求两个值中最大的值但编写的代码不知道哪里错了，求大神帮忙。
已经有10人回复
- 求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）
已经有6人回复
- fluent模拟冷热流体换热出现问题，求大神指导！！！
已经有3人回复
- 求助各位大神fluent的安装问题
已经有9人回复
- 关于fluentUDF的一个小问题！
已经有5人回复
- 有人用fluent计算过自由风扇的流场么？求大神指导
已经有15人回复
- 全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题
已经有8人回复
- Fluent在稳态求解器的时候还能用动网格吗？
已经有4人回复
- Fluent欧拉两相流，内部动态边界条件设定，求大神！
已经有8人回复
- Ｆｌｕｅｎｔ的动网格出现负体积问题求助？
已经有18人回复
- 几道政治问题，求大神解答
已经有6人回复
- fluent安装新问题，本人安装过几次fluent，遇到一个新问题
已经有4人回复
- 导入UDF时出这种错误！！
已经有4人回复
- 接触fluent时间很短。。弱弱的问几个问题 希望有大神解答
已经有13人回复
- GView5.08 CConnectionGFCHK:: 错误
已经有4人回复
2楼2014-02-20 00:14:45
3楼2014-02-20 08:47:10
4楼2014-05-18 20:39:51
5楼2018-06-05 10:10:31
