24小时热门版块排行榜    

查看: 1546  |  回复: 7

逆袭defy

新虫 (初入文坛)

[交流] 利用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的不一样
求大侠指点下,已经想了好几天了还没有头绪
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fujixing

新虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
请问report的平均速度选择的是Area Aver还是Volume Aver

[ 发自小木虫客户端 ]
2楼2015-04-10 01:32:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

逆袭defy

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by fujixing at 2015-04-10 01:32:56
请问report的平均速度选择的是Area Aver还是Volume Aver

我用的是面积的平均,两种面积平均都试了下,而且看图显示的也和report的差不多,就是和我直接udf统计出来的不一样,请问您知道怎么回事吗
3楼2015-04-10 08:59:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

L161019

新虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

4楼2015-04-10 10:52:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

逆袭defy

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by L161019 at 2015-04-10 10:52:14
是不是出口有回流,你取了绝对值!

如果我不取绝对值的话,那样就更小了,比report出来的差距更大

而且我用这种方法统计速度进口的速度特比我实际输入的参数要小一点,是不是我定位到边界的cell thread有问题?
5楼2015-04-10 13:11:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

L161019

新虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

6楼2015-04-10 13:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

逆袭defy

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by L161019 at 2015-04-10 13:34:56
CHUKOUID你给了值了咯...

gas_thread= Lookup_Thread(solid_domain, CHUKOUID);
                chukou_thread = THREAD_T0(gas_thread);

THREAD_T0
感觉下面那句定位到邻近的cell 层这里是不是有问题

7楼2015-04-10 13:37:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

panjunxiu

新虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
我也遇到了,Udf加载的速度与实际monitor的不一样,请问楼主解决了没
8楼2018-11-01 19:58:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 逆袭defy 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 06:50 by xSrQNSSjcgoX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:38 by xSrQNSSjcgoX
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:43 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 10/500 2026-11-15 00:39 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:51 by Equinoxhua
[考博] 考博求助 10+3 社会实践对别 2026-11-08 3/150 2026-11-13 14:04 by 沐小虫252
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[考博] 27博士申请(市政/环境) +3 ladf 2026-11-12 3/150 2026-11-13 08:19 by Equinoxhua
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 9/450 2026-11-12 11:22 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 07:09 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
信息提示
请填处理意见