应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 是误差吗？？？？【FLUENT的一个问题】
181/1返回列表
查看: 3292  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kylafree

至尊木虫 (知名作家)

[交流] 是误差吗？？？？【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 ]
回复此楼

» 猜你喜欢

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

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-02-28 19:11:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

kylafree

至尊木虫 (知名作家)
kylafree: 回帖置顶 2012-02-29 10:13:04
FLUENT中，air的平均分子量为28.966，利用组分输运模型，
我自定义一种分子量为28.966的气体与air按摩尔分数为1：9混合，理应密度也为1.225kg/m3。采用双精度，按上述UDF测试其管内平均密度值为1.17667，UDM记录的管内每一个CELL的密度值都为1.176673870E+00; 与1.225相差挺有点大！！！

再自定义一种分子量为30的气体与air按摩尔分数为1：9混合，理应密度>1.225kg/m3。但按上述UDF测试其管内平均密度值为1.18087，UDM记录的管内每一个CELL的密度值都为1.180874250E+00 ; 与1.225相差也很大！！！

有谁能告诉，上述差值确实是由误差导致吗？？？还是我的测试有问题？？？上述UDF有问题吗？？？

[ Last edited by kylafree on 2012-2-29 at 10:14 ]
赞一下 回复此楼
5楼2012-02-28 20:18:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

heattransfet

新虫 (职业作家)
★
kylafree(金币+1):谢谢参与
xiegangmai(金币+1): 谢谢参与！ 2012-02-28 23:20:32
肯定是误差了，再说算密度也没必要用UDF吧，呵呵
赞一下 回复此楼
2楼2012-02-28 19:29:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
楼: Originally posted by heattransfet at 2012-02-28 19:29:16:
肯定是误差了，再说算密度也没必要用UDF吧，呵呵

误差咋这么大？？？难以接受啊！！！
另外我是在模拟的过程中发现FLUENTT算出的数据与直接用公式计算的数据不相符后，才考虑用UDF作一个测试。一试发现果然相关甚远啊！！！

[ Last edited by kylafree on 2012-2-28 at 19:37 ]
赞一下 回复此楼
3楼2012-02-28 19:34:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
哎，用双精度测试了一下，结果为1.225kg/m3。差别真大呀！！！！

但是用混合物作测试值又不对了

[ Last edited by kylafree on 2012-2-28 at 19:56 ]
赞一下 回复此楼
4楼2012-02-28 19:39:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

salas_1981

铜虫 (正式写手)
★
kylafree(金币+1):谢谢参与
其实，在我们专业看来
已经很准确了
你们是啥专业啊？
哈哈哈
赞一下 回复此楼
6楼2012-02-29 01:48:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
楼: Originally posted by salas_1981 at 2012-02-29 01:48:10:
其实，在我们专业看来
已经很准确了
你们是啥专业啊？
哈哈哈

第四楼的测试结果也很准确吗？？？？差好远啊！！！！
赞一下 回复此楼
9楼2012-02-29 08:42:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
★
kylafree(金币+1):谢谢参与
xiegangmai(金币+2): 鼓励讨论交流！ 2012-02-29 22:25:47
常态下实际空气密度就应该是1.17到1.18左右。
组分运输模型算出来的是对的。
那个1.225应该是某种状态下的值，可能是0℃吧，没验算过。
单精度和双精度，有这么大差别还是蛮正常的。
单精度下1000个0.001相加就不等于1。
赞一下 回复此楼
11楼2012-02-29 13:20:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
楼: Originally posted by 47425804 at 2012-02-29 13:20:07:
常态下实际空气密度就应该是1.17到1.18左右。
组分运输模型算出来的是对的。
那个1.225应该是某种状态下的值，可能是0℃吧，没验算过。
单精度和双精度，有这么大差别还是蛮正常的。
单精度下1000个0.001相加 ...

你好，FLUENT数据库中的空气密度是在一个大气压，15摄氏度时的数据(听别人说工程书上有介绍)，为1.225kg/m3   。这与用理想气体状态方程PV=nRT算出的值一致。
  而且请看baidu介绍如下：
“...我们一般采用的空气密度是指在0摄氏度、绝对标准指标下，密度为1.293kg/m3 　　...　通常情况下，即20摄氏度时，取1.205kg/m3。 ”

这些数据都比1.17大很多啊！！！！！
赞一下 回复此楼
12楼2012-02-29 14:20:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
xiegangmai(金币+1): 谢谢参与！ 2012-02-29 22:25:56
引用回帖:
12楼: Originally posted by kylafree at 2012-02-29 14:20:49:
你好，FLUENT数据库中的空气密度是在一个大气压，15摄氏度时的数据(听别人说工程书上有介绍)，为1.225kg/m3   。这与用理想气体状态方程PV=nRT算出的值一致。
  而且请看baidu介绍如下：
“...我们一般采用的空 ...

恩，但是还要考虑一些问题：
一个大气压到底是按101325Pa算的，还是按照100000Pa算。
海平面才是一个大气压，实际上都要小一些。
空气里还包含CO2等密度较大的气体，虽然含量少，但影响不可忽略。
组分运输模型里是按照N2和O2的混合物算的，不包含CO2。
赞一下 回复此楼
13楼2012-02-29 15:05:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
楼: Originally posted by 47425804 at 2012-02-29 15:05:38:
恩，但是还要考虑一些问题：
一个大气压到底是按101325Pa算的，还是按照100000Pa算。
海平面才是一个大气压，实际上都要小一些。
空气里还包含CO2等密度较大的气体，虽然含量少，但影响不可忽略。
组分运输模 ...

一个大气压按101325算的，在FLUENT材料面板中可以看到，air虽然是混合物，但FLUENT给出了它的平均分子量为28.966，密度为1.225kg/m3。
利用理想气体状态方程PV=nRT计算的话确实28.966的分子量与1.225kg/m3的密度刚好对应上。

P=101325 , R=8.31451(理想气体常数) , T=273.15+15 （15摄氏度）
取n=1 mol ，可算出此时的气体摩尔体积为V=23.64L；
现在取1m3=1000L空气，它的质量为（1000/23.64）*28.966=1.225kg，所以空气的密度为1.225kg/m3，对应的这么好！

但又有一点儿问题，我又算了一下其它几种气体，C2H2,C2H6,NH3,由分子量计算的密度与FLUENT数据库给的密度近似相等，有一定误差；而H2O蒸汽差的比较远（这是不是由于15度时水的饱和蒸汽压比较低的缘故？）

[ Last edited by kylafree on 2012-2-29 at 16:52 ]
赞一下 回复此楼
14楼2012-02-29 16:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
楼: Originally posted by 47425804 at 2012-02-29 15:05:38:
恩，但是还要考虑一些问题：
一个大气压到底是按101325Pa算的，还是按照100000Pa算。
海平面才是一个大气压，实际上都要小一些。
空气里还包含CO2等密度较大的气体，虽然含量少，但影响不可忽略。
组分运输模 ...

你认为我用UDF算得的密度值是可用可信的是吧，我还要用算的密度值来计算其它东西，由于发现它与FLUENT数据库给的数值对不上，所以才有此疑问而发帖
赞一下 回复此楼
15楼2012-02-29 16:55:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
★ ★
xiegangmai(金币+2): 鼓励讨论交流！ 2012-03-01 23:07:28
引用回帖:
14楼: Originally posted by kylafree at 2012-02-29 16:50:18:
一个大气压按101325算的，在FLUENT材料面板中可以看到，air虽然是混合物，但FLUENT给出了它的平均分子量为28.966，密度为1.225kg/m3。
利用理想气体状态方程PV=nRT计算的话确实28.966的分子量与1.225kg/m3 ...

这个恐怕还需要更专业的人士来解答了，气体密度到底还受到什么东西的影响。
可以考虑去查查物性手册，对比对比。

不过，你算平均密度的方法严格说来是不对的。
最正确应该是对质量和体积求和，然后相除得到平均密度。
赞一下 回复此楼
17楼2012-03-01 08:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)
引用回帖:
楼: 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);
}
赞一下 回复此楼
18楼2012-03-01 18:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
afg337楼
2012-02-29 06:51   回复  
kylafree(金币+1):谢谢参与
liuyufeng48楼
2012-02-29 07:30   回复  
kylafree(金币+1):谢谢参与
假大空10楼
2012-02-29 08:55   回复  
kylafree(金币+1):谢谢参与
zhoupeng8716楼
2012-02-29 17:18   回复  
kylafree(金币+1):谢谢参与
相关版块跳转 我要订阅楼主 kylafree 的主题更新
181/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +3 kulium 2026-08-21 4/200 2026-08-21 17:02 by gltch
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳今天，20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
关闭
请填处理意见
关闭