| 查看: 1250 | 回复: 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的不一样 求大侠指点下，已经想了好几天了还没有头绪 |
回复此楼
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
4楼2015-04-10 10:52:14
2楼2015-04-10 01:32:56
3楼2015-04-10 08:59:07
5楼2015-04-10 13:11:03
