24小时热门版块排行榜

返回列表

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

kylafree

至尊木虫 (知名作家)

应助: 166 (高中生)
金币: 17194
帖子: 6163
在线: 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 ]

回复此楼

» 猜你喜欢

青椒八年已不青，大家都被折磨成啥样了？已经有11人回复
限项规定已经有9人回复
免疫学博士有名额，速联系已经有4人回复
交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

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

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

正交试验，空白列作误差列，极值较大，那所做数据有效吗？已经有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 (高中生)
金币: 17194
帖子: 6163
在线: 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);
}