| 查看: 2464 | 回复: 12 | ||
wzhhgd3210银虫 (小有名气)
|
[求助]
UDS及其后处理求助！ 已有2人参与
|
|
|
楼主刚开始接触CFD，现在模拟室内空调通风，需要利用每个节点上已计算出的温度和风速计算求解该点的PMV值并输出云图（该值可以看做室内温度和风速的函数），看到paper上都说用UDF做，我在网上找了一个开源的修改了一下，可以interpret但是迭代计算就会报错Fluent received fatal signal(ACCESS_VIOLATION)。 代码如下： #include "udf.h" enum { pmv }; DEFINE_ADJUST(pmv_myudf,domain) { Thread *t; cell_t c; thread_loop_c(t,domain) { begin_c_loop(c,t) { double t_myudf=C_T(c,t)-273.15; double speed_u_myudf=C_U(c,t); double speed_v_myudf=C_V(c,t); double speed_w_myudf=C_W(c,t); double mh2o_myudf=C_YI(c,t,0); double clo=0.55; double icl_myudf=0.08525; double fcl_myudf; double tcl_myudf; double tcl1_myudf; double tcl2_myudf; double temperary1_myudf; double temperary2_myudf; double hc_myudf; double tr_myudf=t_myudf; double M_myudf=58.0; double W_myudf=0.0; double p_myudf=101325; double pa_myudf; double a_myudf; double b_myudf; double c_myudf; double d_myudf; double e_myudf; double f_myudf; double L_myudf; double speed_myudf=sqrt(pow(speed_u_myudf,2.0)+pow(speed_v_myudf,2.0)+pow(speed_w_myudf,2.0)); pa_myudf=29*mh2o_myudf/(18+11*mh2o_myudf)*p_myudf; if (icl_myudf<0.078) fcl_myudf=1.00+1.290*icl_myudf; else fcl_myudf=1.05+0.645*icl_myudf; tcl1_myudf=40; hc_myudf=1; tcl2_myudf=35.7-0.025*(M_myudf-W_myudf)-icl_myudf*(0.0000000396*fcl_myudf*(pow((tcl1_myudf+273),4)-pow((tr_myudf+273),4))+fcl_myudf*hc_myudf*(tcl1_myudf-t_myudf)); while (tcl1_myudf-tcl2_myudf>0.000001) { tcl1_myudf=tcl2_myudf; temperary1_myudf=2.38*pow((tcl1_myudf-t_myudf),0.25); temperary2_myudf=12.1*pow(speed_myudf,0.5); if (temperary1_myudf<temperary2_myudf) hc_myudf=temperary2_myudf; else hc_myudf=temperary1_myudf; tcl2_myudf=35.7-0.025*(M_myudf-W_myudf)-icl_myudf*(0.0000000396*fcl_myudf*(pow((tcl1_myudf+273),4)-pow((tr_myudf+273),4))+fcl_myudf*hc_myudf*(tcl1_myudf-t_myudf)); } tcl_myudf=tcl2_myudf; a_myudf=0.0000000396*fcl_myudf*(pow((tcl_myudf+273),4.0)-pow((tr_myudf+273),4.0)); b_myudf=fcl_myudf*hc_myudf*(tcl_myudf-t_myudf); c_myudf=0.00305*(5733-6.99*(M_myudf-W_myudf)-pa_myudf); d_myudf=0.42*(M_myudf-W_myudf-58.15); e_myudf=0.000017*M_myudf*(5867-pa_myudf); f_myudf=0.0014*M_myudf*(34-t_myudf); L_myudf=M_myudf-W_myudf-(a_myudf+b_myudf+c_myudf+d_myudf+e_myudf+f_myudf); C_UDSI(c,t,pmv)=(0.303*exp(-0.036*M_myudf)+0.028)*L_myudf; } end_c_loop(c,t) } } 请各位大侠帮忙看一下问题在哪，另外楼主是小白，还想请教一下UDS的具体操作：先interpret代码，然后在Funtion Hooks里的Adjust选项里选中我的UDF激活，再把User-Difined Scalars里面Number添加1，设置边界条件迭代。这样的操作就可以了么？还是说Memory那里也要改？先谢谢各位大侠了~ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
2楼2015-01-03 11:40:53
nhnt
铜虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 456.8
- 散金: 105
- 帖子: 87
- 在线: 60.1小时
- 虫号: 3129029
- 注册: 2014-04-11
- 性别: GG
- 专业: 传热传质学
5楼2015-06-27 16:20:12
3楼2015-01-03 14:26:15
nhnt
铜虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 456.8
- 散金: 105
- 帖子: 87
- 在线: 60.1小时
- 虫号: 3129029
- 注册: 2014-04-11
- 性别: GG
- 专业: 传热传质学
4楼2015-06-26 11:10:33
wangyan2387
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 108
- 帖子: 105
- 在线: 67.9小时
- 虫号: 2679233
- 注册: 2013-09-25
- 性别: MM
- 专业: 工程热物理相关交叉领域
6楼2016-08-31 17:12:08
wangyan2387
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 108
- 帖子: 105
- 在线: 67.9小时
- 虫号: 2679233
- 注册: 2013-09-25
- 性别: MM
- 专业: 工程热物理相关交叉领域
7楼2016-08-31 17:12:59
|
本帖内容被屏蔽 |
8楼2016-09-09 16:50:45
9楼2017-07-04 15:34:22
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
10楼2017-07-06 00:10:59
