| 查看: 1234 | 回复: 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的不一样 求大侠指点下，已经想了好几天了还没有头绪 |
回复此楼» 猜你喜欢
271求调剂
已经有12人回复
-
322求调剂
已经有3人回复
-
环境工程 085701，267求调剂
已经有15人回复
-
材料0856 英一数二 323 求调剂
已经有3人回复
-
301求调剂
已经有3人回复
-
348求调剂
已经有3人回复
-
350求调剂
已经有3人回复
-
309分085801求调剂
已经有7人回复
-
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
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
