| 查看: 1075 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xujinchao铜虫 (小有名气)
|
[求助]
请问fluent两相流拽力模型，加载到三相流中，就出错，请问需要怎么设置或修改 已有1人参与
|
||
谢谢大家 11.png |
回复此楼» 收录本帖的淘帖专辑推荐
 模拟相关 |
» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有13人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组-博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
xujinchao
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 913.3
- 帖子: 52
- 在线: 123.3小时
- 虫号: 1547085
- 注册: 2011-12-22
- 性别: GG
- 专业: 水力学与水信息学
3楼2015-08-13 14:48:13
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; } |
2楼2015-08-12 23:22:34
4楼2015-12-07 22:09:48