24小时热门版块排行榜

返回列表

jyl432

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7
帖子: 2
在线: 1.4小时
虫号: 1970048
注册: 2012-09-03

[交流] fluent计算多相流问题时，UDF能否提取单元内分相速度，以及单元平均速度已有4人参与

请问在vof两相流计算中，因为编写udf是要用到两相之间的速度差（Ug-Ul）。UDF中能否使用命令提取出cell中各分相的速度，或是两相间的滑移速度如何提取。

回复此楼

» 猜你喜欢

自荐读博已经有9人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

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

单机多核运行Fluent，加UDF,运算良好，但保存case,dat再次读取时提示错误！已经有7人回复
【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容) 已经有21人回复

1楼 2013-03-21 17:00:41

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

znon

金虫 (小有名气)

应助: 40 (小学生)
金币: 739.2
红花: 8
帖子: 164
在线: 74.7小时
虫号: 325496
注册: 2007-03-17
性别: GG
专业: 能源化工

★ ★ ★
小木虫: 金币+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;
}