| 查看: 950 | 回复: 3 | ||
[求助]
流化床曳力模型UDF的问题 已有1人参与
|
|
请问曳力模型这个函数中的0.281632和9.076960是如何计算的出的呢？如果是非球形颗粒应该改哪个地方？ #include "udf.h" #define pi 4.*atan(1.) #define diam2 3.e-4 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; } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 重金请教，fluent边界条件设置中液体出口设置问题！！！！
已经有8人回复
- 关于DPM模型模拟液滴的若干问题
已经有16人回复
- Fluent中DPM模型的UDF中使用P_POS()指针，然后给我提示undeclared variable 怎么办?
已经有6人回复
- fluent模拟流化床的固相颗粒泄露问题
已经有3人回复
- Fluent的UDF是否可以自定义辐射模型或湍流模型
已经有11人回复
- 用VOF模型如何模拟不用流速的气液两相流
已经有4人回复
- fluent 中二维模型中，材料属性随着温度变化的udf源代码？
已经有14人回复
- 单机多核并行计算下UDF的问题
已经有9人回复
- 急求：CSF表面张力模型的UDF代码
已经有6人回复
- Fluent模拟流化床两相颗粒之间的曳力系数选择
已经有11人回复
- 一个UDF并行报错
已经有17人回复
- 气液固三相反应器的模拟
已经有21人回复
- Fluent关于DPM模型UDF编程的求助
已经有4人回复
- 欧拉模型和VOF之间的差别
已经有11人回复
- 欧拉双流体模型中如何模拟出布风板？
已经有10人回复
- UDF 改为并行计算
已经有9人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- 曳力函数的的UDF
已经有15人回复
- 流化床曳力系数UDF
已经有9人回复
- 求助：UDF模型的编译问题
已经有15人回复
- fluent中颗粒速度求助
已经有19人回复
- 【求助】通气空泡多相流（气、水蒸汽、水）空化模型的UDF如何入手啊？？？？
已经有4人回复
- 【求助】请问fluent能模拟超重力的情况吗？
已经有4人回复
- 【求助】fluent模拟两段流化床采用欧拉和DPM模型问题
已经有11人回复
- 【求助】Fluent 蒸汽凝结模型求助
已经有3人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
welson
银虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1371.8
- 散金: 22
- 红花: 3
- 帖子: 123
- 在线: 182.9小时
- 虫号: 1579962
- 注册: 2012-01-16
- 性别: GG
- 专业: 传热传质学
2楼2014-05-25 15:21:45
3楼2014-05-25 17:35:43
welson
银虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1371.8
- 散金: 22
- 红花: 3
- 帖子: 123
- 在线: 182.9小时
- 虫号: 1579962
- 注册: 2012-01-16
- 性别: GG
- 专业: 传热传质学
4楼2014-05-29 16:40:08
10