| 查看: 1298 | 回复: 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的不一样 求大侠指点下，已经想了好几天了还没有头绪 |
回复此楼» 猜你喜欢
药学硕士，第一、第二作者已发表6 篇 SCI，药理方向及相关方向2026年/2027年博士申请
已经有6人回复
-
一篇MDPI论文改变了学习工作和生活
已经有5人回复
-
26年博士申请自荐-电催化
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有6人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有4人回复
-
考博
已经有5人回复
-
26年申博自荐-计算机视觉
已经有4人回复
-
药化及相关博士的申请
已经有3人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
6楼2015-04-10 13:34:56
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