| 查看: 623 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
UDF
|
||
|
B：在Fluent中使用Herschel-Bulkley粘性模型： #include "udf.h” real T,vis, s_mag, s_mag_c, sigma_y,n,k; real C_1 = 1.0; real C_2 = 1.0; real C_3 = 1.0; real C_4 = 1.0; int ia ; DEFINE_PROPERTY(hb_viscosity,c,t) { T=C_T(c, t); s_mag = CELL_STRAIN_RATE_MAG(c,t); if (ia==0.0) { C_1 = RP_Get_Real("c_1"  ;C_2 = RP_Get_Real("c_2" ;C_3 = RP_Get_Real("c_3" ;C_4 = RP_Get_Real("c_4" ;ia = 1;} k= C_1 ; n= C_2 ; sigma_y = C_3 ; s_mag_c = C_4 ; if (s_mag < s_mag_c) {vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c)*pow(s_mag_c,(n-1));} else { vis = sigma_y / s_mag + k*pow(s_mag, (n-1));} return vis; } 中{vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c)*pow(s_mag_c,(n-1));}有没有大神可以解释一下，十分感谢！ |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
3楼2015-04-08 15:19:10
|
fluent动网格使用udf，圆柱面边界的速度随时间变化，写了程序，总是编写错误，求大神给看看，其实就是在x轴的速度是sint函数（做的二维流场，左边来流，右边出口，中间圆柱在振动），别的也没啥，c基础太差，求指导，新手不能发帖，希望没有打扰到您，万分感激 #include "udf.h" DEFINE_PROFILE(cylinder_velocity) { cell_t f; begin_f_loop(f,therad) { real t=RP_Get_Real("iteration" ;f=-0.1256*sin(0.0628*t); } end_f_loop(f,therad) } |
2楼2015-04-08 10:19:30
4楼2015-04-08 15:21:17
10