| 查看: 714 | 回复: 0 | |||
xxfreedomxx新虫 (初入文坛)
|
[求助]
求助各位大神，我的wen-yu曳力UDF不能收敛
|
|
求助各位大神，如下是我的UDF，我发现自己写的别的曳力公式和这个wen-yu的UDF是一个情况，我做的是固液搅拌体系，500个时间步之后第二项速度在10-3以上就收敛不了了，慢慢增大，调节收敛因子，时间步长，湍流模型，都不行，单纯用自带模型就没有问题，选择的时候加入了brucato的模型进行修正。在这里卡了两个星期了。 #include "udf.h" #define pi 4.*atan(1.) DEFINE_EXCHANGE_PROPERTY(s_l_drag_IL_sl, cell, mix_thread, s_col, f_col) { Thread *thread_l, *thread_s; real abs_v, slip_x, slip_y, slip_z; real reyp, void_s, rho_l, rho_s, mu_l, fdrag, k_l_s, void_l, diam2, Res; /* find the threads for the liquid(primary) */ /* and solids (secondary phases) */ thread_l = THREAD_SUB_THREAD(mix_thread, f_col); /* liquid phase */ thread_s = THREAD_SUB_THREAD(mix_thread, s_col); /* solid phase */ /* find phase velocities and properties*/ rho_l = C_R(cell, thread_l); mu_l = C_MU_L(cell, thread_l); rho_s = C_R(cell, thread_s); diam2 = C_PHASE_DIAMETER(cell, thread_s); slip_x = C_U(cell, thread_l) - C_U(cell, thread_s); slip_y = C_V(cell, thread_l) - C_V(cell, thread_s); slip_z = C_W(cell, thread_l) - C_W(cell, thread_s); /*compute slip*/ abs_v = sqrt(slip_x*slip_x + slip_y*slip_y + slip_z*slip_z); /*compute Reynold's number*/ reyp = rho_l*abs_v*diam2/mu_l; /* compute particle relaxation time */ void_s = C_VOF(cell, thread_s);/* solid vol frac*/ void_l = C_VOF(cell, thread_l);/* liquid vol frac*/ Res=void_l*reyp; if(Res<=1000){ fdrag = 24.*(1. + 0.15*pow(Res, 0.687)/Res; } else(Res>1000){ fdrag = 0.44; } k_l_s = 0.75*void_s*void_l*mu_l*fdrag*abs_v*pow(void_l, -2.65)/diam2; return k_l_s; }  |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
30