| 查看: 1236 | 回复: 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的不一样 求大侠指点下，已经想了好几天了还没有头绪 |
回复此楼» 猜你喜欢
一志愿同济大学323分（080500）求调剂
已经有3人回复
-
330分求调剂
已经有7人回复
-
336求调剂
已经有3人回复
-
285求调剂
已经有5人回复
-
调剂推荐
已经有9人回复
-
材料调剂
已经有6人回复
-
材料专硕322分
已经有3人回复
-
321求调剂
已经有5人回复
-
301求调剂
已经有6人回复
-
085600，321分求调剂
已经有13人回复
2楼2015-04-10 01:32:56
3楼2015-04-10 08:59:07
L161019 
主管区长
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
4楼2015-04-10 10:52:14
5楼2015-04-10 13:11:03
L161019
超级版主
★
小木虫: 金币+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
