| 查看: 2622 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
流化床曳力系数UDF
|
||
| 简单的液固两相流化床，用欧拉双流体模型，两相间的曳力系数UDF怎么编写，曳力系数公式已经有了，UDF不会编，有没有高手能给个类似的.c文件 |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- [求助]光吸收系数
已经有13人回复
- 请教指定浓度的法向梯度为0的UDF如何写
已经有5人回复
- 不同扫描速率下的循环伏安的图求扩散系数
已经有13人回复
- UDF源项的一个疑问
已经有14人回复
- 曳力函数的的UDF
已经有15人回复
- 氧化铝陶瓷为什么有较低的摩擦系数？
已经有3人回复
- 吸附的动力学方程拟合
已经有18人回复
- udf 重力方向变化
已经有17人回复
- 【求助】udf读取混合物密度失败、及平均流速读取？
已经有8人回复
- 【求助】Warning: incorrect cg motion UDF value on zone 35
已经有5人回复
- 不同温度下丙酮的表面张力系数
已经有6人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
3楼2011-10-18 15:13:21
2楼2011-10-18 15:10:47
4楼2011-10-18 19:52:00
【答案】应助回帖
★
mokalin(金币+5): 谢谢 2011-10-19 13:47:24
spc08(金币+1): 感谢应助 2011-10-19 15:18:29
mokalin(金币+5): 谢谢 2011-10-19 13:47:24
spc08(金币+1): 感谢应助 2011-10-19 15:18:29
|
就是fluent自带的UDF手册，讲DEFINE_EXCHANGE_PROPERTY宏的那一节，有自定义曳力模型的例子。 我也给你贴上来吧 /****************************************************************** UDF for customizing the default Syamlal drag law in ANSYS FLUENT *******************************************************************/ #include "udf.h" #define pi 4.*atan(1.) #define diam2 3.e-4 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; } |
5楼2011-10-19 13:21:49