| 查看: 1738 | 回复: 5 | ||
[求助]
曳力模型udf要并行计算，该怎么修改呢
|
|
这是我单机进行计算时用的udf代码： #include "udf.h" #define pi 4.*atan(1.) #define diam2 3.e-4 DEFINE_EXCHANGE_PROPERTY(custom_drag3,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; } 现在需要换到服务器上进行并行计算，需要加一些语句，该怎么加呢？哪位知道的帮帮忙啊 [ 来自科研家族 流体数值模拟 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 还原力测定时出现沉淀是怎么回事
已经有6人回复
- 图谱修改
已经有3人回复
- 用ANSYS画2D模型时如何在里面画导流板（即隔板）
已经有5人回复
- 关于决策树分类模型
已经有5人回复
- 大家使用哪些免费软件建立模型呢？
已经有11人回复
- 高效液相Purge时压力高怎么办
已经有33人回复
- journal of membrane science 修改稿求助
已经有14人回复
- PK/PD模型是否必须在病理动物上做??
已经有6人回复
- 制作的纳米颗粒已经溶在水中了，怎么提取做XRD呢
已经有34人回复
- 论文被拒了。。两天后，编辑又来信，让重新修改，然后再以信文章投稿。。
已经有36人回复
- 曳力函数的的UDF
已经有15人回复
- 流化床曳力系数UDF
已经有9人回复
- 求助：UDF模型的编译问题
已经有15人回复
- 【求助】用vof模型计算出的界面，有木有udf语句可以找得到
已经有6人回复
- 【求助】单机多核UDF怎么并行计算
已经有15人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
2楼2012-03-30 17:17:50
3楼2012-03-30 20:39:13
4楼2012-03-30 20:50:57
5楼2012-04-11 16:53:35
6楼2015-10-18 20:16:07
10