24小时热门版块排行榜    

查看: 608  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了? +7 chengyan1220 2026-08-21 8/400 2026-08-21 15:19 by chengyan1220
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 13/650 2026-08-21 08:30 by 北京莱茵编辑
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见