24小时热门版块排行榜    

查看: 490  |  回复: 3

1991xiao11

铜虫 (初入文坛)

[求助] 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));}有没有大神可以解释一下,十分感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

f5287028

铁虫 (初入文坛)

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rayamay

新虫 (初入文坛)

宣告好像要放前面??

#include "udf.h"
real t;
DEFINE_PROFILE(cylinder_velocity)
{ cell_t f;
     begin_f_loop(f,therad)     
     {
       t=RP_Get_Real("iteration";   /*check*/
       f=-0.1256*sin(0.0628*t);
      }
    end_f_loop(f,therad)

}
3楼2015-04-08 15:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rayamay

新虫 (初入文坛)

抱歉沒注意到還有別的
先COMPILE看看能不能用

#include "udf.h"
real t;
cell_t f;
DEFINE_PROFILE(cylinder_velocity)
{
     begin_f_loop(f,therad)     
     {
       t=RP_Get_Real("iteration";   /*check*/
       f=-0.1256*sin(0.0628*t);
      }
    end_f_loop(f,therad)
}
4楼2015-04-08 15:21:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1991xiao11 的主题更新
信息提示
请填处理意见