| 查看: 4617 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
zhangxu86铜虫 (小有名气)
|
[求助]
曳力函数的的UDF 已有2人参与
|
||
|
请问各位大侠，我的模型使用的是欧拉模型，物质是沙和水，请问一下，编写沙水两者之间的曳力时，是使用液液交换系数还是使用液体-固体交换系数呢？下面的是我看到一篇流化床立体上的程序，不知道为什么计算两次K和f呢？上面一段程序使用的就是Syamlal-O‘Brien，下面一段就不明白了。。。。。 /* This is a routine for customizing default Syamlal drag law in Fluent 6. The default drag law uses 0.8 (for void<=0.85) and 2.65 (void>0.85) for bfac. This is for a min fluid vel of 25 cm/s. The current drag law has been tuned for a min fluid vel of 8 cm/s and uses 0.28 and 9.07 for these parameters. */ #include "udf.h" #include "sg_mphase.h" # define pi 4.*atan(1.) #define diam2 3.e-4 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; /* find the threads for the gas (primary) and solids (secondary phases). These phases appear in columns 2 and 1 in the Interphase panel respectively*/ 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 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.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; } DEFINE_EXCHANGE_PROPERTY(custom_drag_ihme, 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, cd, eg, void_g, k_g_s; /* find the threads for the gas (primary) and solids (secondary phases). These phases appear in columns 2 and 1 in the Interphase panel respectively*/ 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 reynolds number*/ reyp = rho_g*abs_v*diam2/mu_g; cd = (24./(reyp+SMALL)) + 5.48*pow((reyp+SMALL),-0.573) + 0.36; void_g = C_VOF(cell, thread_g);/* gas vol frac*/ eg = pow(void_g,-2.65); k_g_s = (3./4.)*(cd*(1.-void_g)*abs_v*rho_g*eg)/diam2; return k_g_s; } |
回复此楼» 收录本帖的淘帖专辑推荐
 我的收藏 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看
已经有4人回复
- 流化床曳力系数UDF
已经有9人回复
- udf profile update interval 是什么意思？
已经有5人回复
- 使用UDF后，初始化错误
已经有11人回复
- 重金！关于udf的问题！
已经有5人回复
- udf 重力方向变化
已经有17人回复
- UDF 编译出现不认识函数名称，显示变量问题
已经有5人回复
- 【讨论】Fluent边界条件及虚拟质量力
已经有6人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
- 【求助】密度函数的UDF
已经有6人回复
- 【求助】Warning: incorrect cg motion UDF value on zone 35
已经有5人回复
- 【讨论】请问大家用UDF的时候，有没有碰到如下代码，很重要的
已经有7人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
liubao2012
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2126.3
- 散金: 10
- 红花: 1
- 帖子: 208
- 在线: 123.9小时
- 虫号: 1886137
- 注册: 2012-07-10
- 性别: GG
- 专业: 化学反应工程
9楼2015-12-28 22:14:11
2楼2011-10-25 22:52:01
zhangxu86
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 160
- 帖子: 115
- 在线: 33.5小时
- 虫号: 1279329
- 注册: 2011-04-27
- 性别: MM
- 专业: 水文、水资源
3楼2011-10-26 09:12:50
4楼2011-10-26 10:13:08
