| 查看: 1405 | 回复: 7 | |||
[交流]
利用udf统计出来的速度和直接report出来的速度差很多 已有3人参与
|
 我利用udf读取边界的速度，代码如下， DEFINE_ADJUST(gasification, domain) { cell_t c; Thread *gas_thread; Thread *jinkou_thread; Thread *chukou_thread; Domain *mixure_domain,*solid_domain; real zhiliangfenshu1=0,tiji1=0; real zhiliangfenshu2=0,tiji2=0; real chukousuduzonhe=0; real chukousuduzonhe1=0; real cellshuliang=0; solid_domain=Get_Domain(3); gas_thread= Lookup_Thread(solid_domain, CHUKOUID); chukou_thread = THREAD_T0(gas_thread); begin_c_loop(c,chukou_thread) { chukouwenduzonghe+=C_T(c,chukou_thread)*C_VOLUME(c,chukou_thread); chukousuduzonhe+=fabs(C_U(c,chukou_thread)*C_VOLUME(c,chukou_thread));/*注意这里用的是x方向的速度作为出口速度，根据mesh的具体坐标而定*/ chukousuduzonhe1+=fabs(C_U(c,chukou_thread)); cellshuliang+=1; tiji1+=C_VOLUME(c,chukou_thread); } end_c_loop(c,chukou_thread) chukousudu1=chukousuduzonhe1/cellshuliang; chukousudu=chukousuduzonhe/tiji1; Message ("出口固体速度: %g\n",chukousudu); Message ("出口固体速度平均: %g\n",chukousudu1); } 我用两种方法求速度，一个是把 速度*体积 加到一起然后除以总体积，另外一种是把所有的数量加一起，然后除以cell总数量，两种方法求出来的边界的速度很接近，但是和利用flunet直接report出来的速度差的很多，，但是这两种方法求出来的温度和质量分数和report出来的都是一样的，只有速度和report的不一样 求大侠指点下，已经想了好几天了还没有头绪 |
回复此楼» 猜你喜欢
确定了，国自然21号放榜
已经有4人回复
-
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
2楼2015-04-10 01:32:56
3楼2015-04-10 08:59:07
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
4楼2015-04-10 10:52:14
5楼2015-04-10 13:11:03
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
6楼2015-04-10 13:34:56
7楼2015-04-10 13:37:34
8楼2018-11-01 19:58:44
100