24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kylafree

至尊木虫 (知名作家)

应助: 166 (高中生)
金币: 16959
帖子: 6120
在线: 1161.6小时
虫号: 1061079

[交流] 是误差吗？？？？【FLUENT的一个问题】

FLUENT的材料库中，可看到air的密度为1.225kg/m3，
现在我在一个圆柱管模型中，采用air作为流体，初始化后，管内全为air；
执行以下UDF，应可算出管内air的平均密度，实际算出的值为1.2353kg/m3。
这两个值不一样，是误差吗？好象误差挺大啊！！！！

采用双精度测试了一下，结果为1.225kg/m3。这个数值倒是准确。

但是用混合物作测试值又不对了，请见下面一楼：

[ Last edited by kylafree on 2012-2-28 at 19:56 ]

#include "udf.h"
DEFINE_ON_DEMAND(st2)
{
cell_t c;
Thread *t = Lookup_Thread(Get_Domain(1), 8);
real s=0.0;
int n=0;
begin_c_loop(c,t)
{
   n=n+1;
   s=s+C_R(c,t);
   C_UDMI(c,t,0)=C_R(c,t);
}
end_c_loop(c,t)
s=s/n;
Message("\n s n =%g %d\n", s,n);
}

[ Last edited by kylafree on 2012-2-29 at 10:16 ]

回复此楼

» 猜你喜欢

寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

正交试验，空白列作误差列，极值较大，那所做数据有效吗？已经有13人回复
PL谱和透射吸收谱表征材料禁带宽度的时候，误差多大？已经有8人回复
Design expert 设计实验，需要自己做平行吗？已经有14人回复
正交试验设计问题已经有11人回复
元素分析误差分析已经有13人回复
帮忙看下图上“工”是什么意思？是最大、最小值吗？origin怎么弄上去的？谢谢！已经有11人回复
ORIGIN里做误差棒已经有12人回复
【求助】配标曲药物各浓度时是直接从母液稀释还是从母液逐各稀释已经有11人回复
【求助】利用GPC测试PLA-b-PEG的分子量，为甚麽误差很大已经有9人回复
Origin怎样做误差线啊。怎么算？已经有3人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2012-02-28 19:11:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

应助: 166 (高中生)
金币: 16959
帖子: 6120
在线: 1161.6小时
虫号: 1061079

引用回帖:

楼: Originally posted by 47425804 at 2012-03-01 08:34:08:
这个恐怕还需要更专业的人士来解答了，气体密度到底还受到什么东西的影响。
可以考虑去查查物性手册，对比对比。

不过，你算平均密度的方法严格说来是不对的。
最正确应该是对质量和体积求和，然后相除得到 ...

按你说的求平均密度。空气分子量为28.966，自定义一种分子量为28.966的气体，作测试，代码如下，单精度，直接对密度平均，值为1.18422，密度*体积求和后再除总体积，值为1.18864。还是与1.225相差很多。

#include "udf.h"
DEFINE_ON_DEMAND(st2)
{
cell_t c;
Thread *t = Lookup_Thread(Get_Domain(1), 8);
real m,s=0.0,q=0.0,r=0.0;
int n=0;
begin_c_loop(c,t)
{
   n=n+1;
   s=s+C_R(c,t);
   C_UDMI(c,t,0)=C_R(c,t);
   r=r+C_VOLUME(c,t);
   q=q+C_VOLUME(c,t)*C_R(c,t);
}
end_c_loop(c,t)
s=s/n;
m=q/r;
Message("\n s m n =%g %g %d\n", s,m,n);
}