| 查看: 3499 | 回复: 15 | |||
| 【悬赏金币】回答本帖问题，作者zhangxu86将赠送您 5 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
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; } |
回复此楼» 收录本帖的淘贴专辑推荐
 我的收藏 |
» 猜你喜欢
化学口 B06 面上开始会评了吗
已经有11人回复
-
今年国自必中，不中回家种田
已经有15人回复
-
Journal of Power Sources只有一个审稿人接受
已经有4人回复
-
普通老师、教研室主任、院长副院长都要做很多杂事，学科带头人、专业负责人最舒适
已经有12人回复
-
工材学部啥时候会评
已经有7人回复
-
现在申请CSC的国外访问学者，通过率是不是 很低，主要看重什么？？
已经有10人回复
-
博士有找工作搭子吗？
已经有19人回复
-
如何管理好研究生
已经有10人回复
-
让有意选自己当本科论文导师的同学联系自己VS看中某同学单独问TA是否愿意
已经有4人回复
-
孩子，我要告诉你
已经有30人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看
已经有4人回复
- UDF代码错误在哪里啊？
已经有6人回复
- 流化床曳力系数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人回复
- 【求助】UDF编译成功但导入边界条件时报错【已解决】
已经有23人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
- 【讨论】共享几个UDF
已经有28人回复
yangshuaibyx
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1740.5
- 散金: 106
- 帖子: 167
- 在线: 645.8小时
- 虫号: 2513219
- 注册: 2013-06-19
- 性别: GG
- 专业: 机械测试理论与技术
7楼2014-06-21 21:45:14
10