| 查看: 3948 | 回复: 4 | ||
klinedancer铜虫 (小有名气)
|
[求助]
fluent 自定义曳力系数
|
|
欧拉双流体模型自定义曳力系数，现在曳力系数模型是CD=CD0 *（-31.4565+32.4451void_g); （其中，CD为曳力系数，CD0为单颗粒曳力系数，void_g是空隙率，请问用udf怎么编写） DEFINE_EXCHANGE_PROPERTY(custom_drag_syam, cell, mix_thread, s_col, f_col) 宏返回值到底是什么东西？？是曳力，还是曳力系数，，还是什么别的东西？？求解 如果我的用DEFINE_EXCHANGE_PROPERTY(custom_drag_syam, cell, mix_thread, s_col, f_col)编写，是直接return CD=CD0 *（-31.4565+32.4451*void_g)吗？？ 我翻阅别人的如下 /* 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; } 我的函数返回值return应该是什么呢？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么用fluent直接得到的压力系数曲线是反的？
已经有8人回复
- 关于fluent模拟流化床的密度云图的问题
已经有8人回复
- FLUENT 中如何在VOF模型中嵌入曳力模型或者UDF自定义函数
已经有7人回复
- 怎样在Fluent里建立新的粒子材料
已经有6人回复
- 流化床气固模拟
已经有9人回复
- 细颗粒的流化床的带出速度计算
已经有5人回复
- 两相流欧拉模型模拟气液分层流动，第二相diameter如何确定
已经有10人回复
- fluent热辐射系数设置
已经有8人回复
- Fluent模拟流化床不同高度压力
已经有30人回复
- 你为何选择模拟？
已经有19人回复
- fluent如何设置不同直径的颗粒流？
已经有6人回复
- Fluent中进行多相流模拟时，只能进行瞬态模拟，不能进行稳态模拟吗？
已经有9人回复
- Fluent欧拉模型中相的体积分数的设定
已经有16人回复
- Fluent模拟流化床两相颗粒之间的曳力系数选择
已经有11人回复
- 关于气固两相流模拟中颗粒雷诺数的问题
已经有3人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- fluent中力矩系数的含义
已经有13人回复
- 曳力模型udf要并行计算，该怎么修改呢
已经有5人回复
- 曳力函数的的UDF
已经有15人回复
- 流化床曳力系数UDF
已经有9人回复
- 【讨论】Fluent边界条件及虚拟质量力
已经有6人回复
2楼2015-01-06 10:09:57
3楼2015-05-12 16:10:13
4楼2015-12-10 16:01:03
-183413
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3629.3
- 红花: 3
- 帖子: 224
- 在线: 139小时
- 虫号: 4037969
- 注册: 2015-08-26
- 性别: GG
- 专业: 化学反应工程
5楼2020-08-14 16:41:10
20