| 查看: 1364 | 回复: 2 | ||
哇塞骰铜虫 (小有名气)
|
[求助]
类似空话模型（cavitation）的udf，求大家帮忙看看对不对 已有1人参与
|
|
在用fluent做管内汽水两相流。实际系统比较复杂，所以入口条件不是特别明确。猜测进口是满管饱和水，流动过程中发生汽化，然后又气体又液化成水。但是因为整个流域温差很小，暂时不考虑汽化潜热所以汽化模型应该类似于空化。因为用vof模型来追踪气液相界，所以在使用vof的同时不能使用cavitation模型。就想着自己写个udf。 这个udf是按照沸腾的udf改编的，并将最后的能量源项去除不考虑，求大家帮忙看看有没有什么问题。相变参考量用的压力。主相是液相，次相是汽相。 编译没有问题，但是计算总是发散。求大家帮忙啦。感激不尽。 #include "udf.h" #include "sg_mphase.h" #define P_SAT 12826 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_P(cell, mix_th)<P_SAT) “混合区单元压力小于汽化压力，则液相发生汽化” { m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_P(cell, mix_th) -P_SAT)/P_SAT; dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_P(cell, mix_th) - P_SAT)/P_SAT; } else { “混合区单元压力大于汽化压力，则汽相液化” m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(P_SAT-C_P(cell,mix_th))/P_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_P(cell, mix_th)<T_SAT) “混合区单元压力小于汽化压力，则液相汽化” { m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_P(cell, mix_th) - P_SAT)/P_SAT; dS[eqn] = 0.; } else { “混合区单元压力大于汽化压力，则汽相液化” m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(P_SAT-C_P(cell,mix_th))/P_SAT; dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_P(cell, mix_th) - P_SAT)/P_SAT; } return m_dot_v; } |
回复此楼» 猜你喜欢
一志愿东华大学控制学硕320求调剂
已经有3人回复
-
考研化学学硕调剂，一志愿985
已经有7人回复
-
333求调剂
已经有6人回复
-
化学工程321分求调剂
已经有21人回复
-
0703化学调剂
已经有15人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有14人回复
-
326求调剂
已经有4人回复
-
311求调剂
已经有13人回复
-
311求调剂
已经有7人回复
-
材料与化工（0856）304求B区调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- udf不同壁面热流求助
已经有14人回复
- UDF并行编译求助-附程序段
已经有3人回复
- [已完结]UDF计算species 模型出错
已经有3人回复
- fluent14编译UDF时总是出错，求高手啊！！！难道没人遇到过吗？？？
已经有11人回复
- fluent 中二维模型中，材料属性随着温度变化的udf源代码？
已经有14人回复
- 急求：CSF表面张力模型的UDF代码
已经有6人回复
- fluent udf的应用
已经有20人回复
- 求助！请问怎么用UDF控制mixture模型第二相进口的volume fraction
已经有6人回复
- 全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题
已经有8人回复
- Fluent关于DPM模型UDF编程的求助
已经有4人回复
- 如何在udf里面 判断界面的位置啊（vof模型）
已经有9人回复
- 求助udf的输入输出
已经有16人回复
- vof模型算气液两相逆流 使用udf加动量源项 碰到问题 请高手解答
已经有14人回复
- 求助！！！我想要在fluent中针对一个面进行积分，那么UDF应该用什么宏函数？怎么用？
已经有5人回复
- 求助：UDF模型的编译问题
已经有15人回复
- udf中的平方怎么写
已经有5人回复
- 关于linux系统下使用UDF的问题
已经有7人回复
- 【求助】通气空泡多相流（气、水蒸汽、水）空化模型的UDF如何入手啊？？？？
已经有4人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
哇塞骰: 金币+10, ★★★★★最佳答案, 先给十个，谢谢回答 2014-04-12 10:02:56
感谢参与，应助指数 +1
哇塞骰: 金币+10, ★★★★★最佳答案, 先给十个，谢谢回答 2014-04-12 10:02:56
|
1，这个网上流传很广的udf例子，最多只能算个现象出来，这种方式算出来的相变速率（其实就是传热速率）与实际相差甚远，这跟他的一些常数取值有关（0.1， 这个值的取法目前看不出出处）。 2，如果你将你研究对象的相变类比成空化的话，建议参考空化的理论模型编写udf，帮助文档里对采用的空化模型的数学表达式有详细的描述。 3，空化模型有相应的使用范围，如空化体积分数等，超出一定限度不仅计算结果可能不合理甚至收敛也是很大的问题。 4，新版本里，vof模型下已经允许采用空化模型 |
2楼2014-04-11 22:35:53
哇塞骰
铜虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 134.8
- 散金: 25
- 帖子: 152
- 在线: 43.2小时
- 虫号: 2293475
- 注册: 2013-02-22
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2014-04-12 10:08:52
40