24小时热门版块排行榜

返回列表

xujinchao

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 913.3
帖子: 52
在线: 123.3小时
虫号: 1547085
注册: 2011-12-22
性别: GG
专业: 水力学与水信息学

[求助] 请问fluent两相流拽力模型，加载到三相流中，就出错，请问需要怎么设置或修改已有1人参与

谢谢大家

11.png

回复此楼

» 收录本帖的淘帖专辑推荐

模拟相关

» 猜你喜欢

2025冷门绝学什么时候出结果已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2015-08-11 16:03:05

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

xujinchao

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 913.3
帖子: 52
在线: 123.3小时
虫号: 1547085
注册: 2011-12-22
性别: GG
专业: 水力学与水信息学

/***************************************************************
UDF for customizing the Gidaspow drag law in Fluent
****************************************************************/

#include "udf.h"
#include "sg_mphase.h"

#define diam2 1.7e-5

DEFINE_EXCHANGE_PROPERTY(custom_drag_syam, 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;

thread_g = THREAD_SUB_THREAD(mix_thread, s_col);
thread_s = THREAD_SUB_THREAD(mix_thread, f_col);

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 reynolds 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.8*pow(void_g, 1.28);
else
  bfac = pow(void_g, 2.65);

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;
}