| 查看: 3673 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
fluent计算多相流问题时，UDF能否提取单元内分相速度，以及单元平均速度 已有4人参与
|
|||
| 请问在vof两相流计算中，因为编写udf是要用到两相之间的速度差（Ug-Ul）。UDF中能否使用命令提取出cell中各分相的速度，或是两相间的滑移速度如何提取。 |
回复此楼» 猜你喜欢
什么时候开奖？
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
lsqtcc
金虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 448.1
- 散金: 6
- 帖子: 587
- 在线: 135.3小时
- 虫号: 1448221
- 注册: 2011-10-18
- 性别: GG
- 专业: 工程热物理相关交叉领域
5楼2015-03-17 23:19:37
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
xiegangmai: 金币+2, 谢谢参与！ 2013-03-24 14:32:09
小木虫: 金币+0.5, 给个红包，谢谢回帖
xiegangmai: 金币+2, 谢谢参与！ 2013-03-24 14:32:09
|
曳力系数的例子中就有，你看一下就知道该如何处理了。 DEFINE_EXCHANGE_PROPERTY(custom_drag,cell,mix_thread,s_col,f_col) { Thread *thread_g, *thread_s; real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y, rho_g, rho_s, mu_g, reyp, afac, bfac, void_g, vfac, fdrgs, taup, k_g_s; /* find the threads for the gas (primary) */ /* and solids (secondary phases) */ thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/* gas phase */ thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/ /* find phase velocities and properties*/ x_vel_g = C_U(cell, thread_g); y_vel_g = C_V(cell, thread_g); x_vel_s = C_U(cell, thread_s); y_vel_s = C_V(cell, thread_s); slip_x = x_vel_g - x_vel_s; slip_y = y_vel_g - y_vel_s; rho_g = C_R(cell, thread_g); rho_s = C_R(cell, thread_s); mu_g = C_MU_L(cell, thread_g); /*compute slip*/ abs_v = sqrt(slip_x*slip_x + slip_y*slip_y); /*compute Reynold's number*/ reyp = rho_g*abs_v*diam2/mu_g; /* compute particle relaxation time */ taup = rho_s*diam2*diam2/18./mu_g; void_g = C_VOF(cell, thread_g);/* gas vol frac*/ /*compute drag and return drag coeff, k_g_s*/ afac = pow(void_g,4.14); if(void_g<=0.85) bfac = 0.281632*pow(void_g, 1.28); else bfac = pow(void_g, 9.076960); vfac = 0.5*(afac-0.06*reyp+sqrt(0.0036*reyp*reyp+0.12*reyp*(2.*bfac- afac)+afac*afac)); fdrgs = void_g*(pow((0.63*sqrt(reyp)/ vfac+4.8*sqrt(vfac)/vfac),2))/24.0; k_g_s = (1.-void_g)*rho_s*fdrgs/taup; return k_g_s; } |
2楼2013-03-21 17:47:57
malu430
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1154.4
- 散金: 65
- 红花: 7
- 帖子: 235
- 在线: 133.1小时
- 虫号: 1277136
- 注册: 2011-04-25
- 性别: MM
- 专业: 多相流热物理学
3楼2013-03-22 14:46:19
4楼2013-03-22 20:08:14
