24小时热门版块排行榜    

查看: 760  |  回复: 3

sxfisgod

金虫 (初入文坛)

[求助] 流化床曳力模型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;
}
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

welson

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
fluent 15里面symalal o'brian para可以算c和d。
也有引文介绍c和d是怎么算的。

比较复杂,具体我也没细看。
2楼2014-05-25 15:21:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxfisgod

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by welson at 2014-05-25 15:21:45
fluent 15里面symalal o'brian para可以算c和d。
也有引文介绍c和d是怎么算的。

比较复杂,具体我也没细看。

我用的还是6.3额
3楼2014-05-25 17:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

welson

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sxfisgod: 金币+50, ★★★★★最佳答案 2014-05-29 21:11:27
引用回帖:
3楼: Originally posted by sxfisgod at 2014-05-25 17:35:43
我用的还是6.3额...

我贴出来给你看看。

这个文献“The Derivation of a Drag Coefficient Formula from Velocity-Voidage Correlations”上面好像也有。
流化床曳力模型UDF的问题
2.png


流化床曳力模型UDF的问题-1
22.png

4楼2014-05-29 16:40:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sxfisgod 的主题更新
信息提示
请填处理意见