| 查看: 428 | 回复: 2 | ||
[求助]
VOF气液逆流UDF请教
|
|
#include "udf.h" #include "sg_mphase.h" # define pi 4.*atan(1.) #define diam2 2.8 DEFINE_SOURCE(drag,c,t,dS,eqn) { Thread *thread_g, *thread_l; real x_vel_g, x_vel_l, y_vel_g, y_vel_l, abs_v, slip_x, slip_y, rho_g, rho_l, mu_g, reyp, fsc,flg void_g, drag ; /* 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_l = THREAD_SUB_THREAD(mix_thread, f_col);/* liquid phase*/ /* find phase velocities and properties*/ x_vel_g = C_U(cell, thread_g); y_vel_g = C_V(cell, thread_g); x_vel_l = C_U(cell, thread_l); y_vel_l = C_V(cell, thread_l); slip_x = x_vel_g - x_vel_s; slip_y = y_vel_g - y_vel_s; rho_g = C_R(cell, thread_g); rho_l = 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; if(y_vel_g=0) fsc = 16./reyp; else 1/sqrt(fsc) = 4.*log(reyp*sqrt(fsc))-0.4; flg = fsc(1.+24*pow(rho_l/rho_g,1/3)*/diam2); void_g = C_VOF(cell, thread_g);/* gas vol frac*/ if(void_g=0.5) drag = 0.5*rho_g*flg*diam2*fabs(abs_v)*abs_v; dS(eqn) = -1.*rho_g*flg*diam2*fabs(abs_v); else drag = 0 return drag } 为什么一直提示 38行 thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/*gas phase*/ thread_l = THREAD_SUB_THREAD(mix_thread, f_col);/* liquid phase*/ non-interger subscript expression :float 错误 麻烦大师帮忙看看这段  |
回复此楼» 猜你喜欢
基金申报
已经有5人回复
-
基金委咋了？2026年的指南还没有出来？
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
-
常年博士招收(双一流，工科)
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 多相流边界只允许液体通过，不允许气体通过
已经有6人回复
- vof模型中，相间作用力是怎么处理的，可以自己编udf吗？
已经有3人回复
- 如何在udf里面 判断界面的位置啊（vof模型）
已经有9人回复
- vof模型算气液两相逆流 使用udf加动量源项 碰到问题 请高手解答
已经有14人回复
- multiphase+相变UDF-by（vof-hrbeu)
已经有165人回复
- 【求助】用vof模型计算出的界面，有木有udf语句可以找得到
已经有6人回复
2楼2013-05-06 16:07:52
3楼2013-05-06 19:13:06