24小时热门版块排行榜    

查看: 816  |  回复: 9
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 气液两相流 的 13 个金币

气液两相流

铜虫 (小有名气)

[交流] 【求助】模拟过程中出错

error:floating point error :invalid number
error:object:()
请大侠指点下这个错误是什么原因造成的?我在用欧拉欧拉模型做气液两相流动的非稳定计算时遇到的。

[ Last edited by adu886886 on 2010-3-15 at 08:40 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuming524

专家顾问 (著名写手)

石化蛀虫


adu886886(金币+1):谢谢交流,欢迎常来仿真编程版 2010-03-15 10:51
浮点错误
2楼2010-03-15 10:19:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuhailangyu

荣誉版主 (著名写手)

小木虫蛀虫防疫站站长

优秀版主

这个问题已经有很多帖子了,不是很容易就能说清了就解决的

因为原因是多方面的,你可以搜一下原来的帖子,学习一下

字面意思大家都懂,呵呵
有能力有脾气的人--怀才不遇;有能力没脾气的人--春风得意;没能力有脾气的人--一事无成;没能力没脾气的人--贵人相助。
3楼2010-03-15 10:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuming524

专家顾问 (著名写手)

石化蛀虫

★ ★
kuhailangyu(金币+2):从哪弄这么大一段,厉害,谢谢呵呵 2010-03-15 13:39
气液两相流(金币+5): 2010-03-16 10:11
原理】     
从数值计算方面看,计算机所执行的运算在计算机内是以浮点数(floating point number)来表示的。那些由于用户的非法数值计算或者所用计算机的限制所引起的错误称为floating point error。
1)非法运算:最简单的例子是使用Newton Raphson方法来求解f(x)=0的根时,如果执行第N次迭代时有,x=x(N),f’(x(N))=0,那么根据公式x(N+1)=x(N)-f(x(N))/ f’(x(N))进行下一次迭代时就会出现被0除的错误。
2)上溢或下溢:这种错误是数据太大或太小造成的,数据太大称为上溢,太小称为下溢。这样的数据在计算机中不能被处理器的算术运算单元进行计算。
3)舍入错误:当对数据进行舍入时,一些重的数字会被丢失并且不可再恢复。例如,如果对0.1进行舍入取整,得到的值为0,如果再对它又进行计算就会导致错误。
【避免方法】
计算和迭代:设一个比较小的时间步长会比较好的。或者改成小的欠松驰因子也会比较好。经验来看,把欠松驰因子设为默认值的1/3;降低欠松驰因子或使用耦合隐式求解;改变欠松驰因子,如果是非稳态问题可能是时间步长太大;改善solver-control-limits比例或许会有帮助;你需要降低Courant数;如果仍然有错误,不选择compute from初始化求解域,然后单击init。再选择你想从哪个面初始化并迭代,这样应该会起作用。另外一个原因可能是courant数太大,就样就是说两次迭代之间的时间步太大并且计算结果变化也较大(残差高)。
网格问题:缩放网格时就会发生这个错误。在Gambit中,所有的尺寸都是以mm为单位,在fluent按scale按钮把它转换成m,然后迭代几百次时就会发生这种错误。但是当不把网格缩放到m时,让它和在Gambit中一样,迭代就会成功;我认为你应当检查网格,网格数太多了,使用较少的网格问题就会解决;网格太多,计算机资源不够用,使使比较粗的网格。
边界条件:一个wall边界条件来代迭axis边界条件,结果fluent拒绝计算并告诉floating point error。你的边界条件不能代表真实的物理现象;错误的边界条件定义可能会导致floating point error。例如把内边界设成interior;使用对称边界条件模拟2D区间时也遇到这种问题,把symmetry设为axe symmetric,就发生了floating point error;检查你设的湍流参数,减小湍流强度,先进行50次迭代。
多处理器问题,进行多处理器模拟时也遇到相似的问题。问题的解决方法是在单个处理器上运行,这样就运算得很好。
错误迭代:
以错误的条件来初始化,在开始迭代时就会发生floating point error。

以上别人的观点,个人认为计算策略也很重要,先稳态计算看看~~
4楼2010-03-15 11:09:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wgf332

铁杆木虫 (著名写手)


kuhailangyu(金币+1):谢谢参与 2010-03-15 13:39
气液两相流(金币+1): 2010-03-19 09:59
1:可能是网格问题,调整网格,检查网格质量,基本上要控制在0.8一下。
2:重新初始化,改变初始化方式
3:试试双精度
5楼2010-03-15 11:37:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

气液两相流

铜虫 (小有名气)

引用回帖:
Originally posted by wuming524 at 2010-03-15 11:09:44:
原理】     
从数值计算方面看,计算机所执行的运算在计算机内是以浮点数(floating point number)来表示的。那些由于用户的非法数值计算或者所用计算机的限制所引起的错误称为floating point error。
1)非法 ...

谢谢你的回复。
6楼2010-03-16 11:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

气液两相流

铜虫 (小有名气)

引用回帖:
Originally posted by wuming524 at 2010-03-15 10:19:52:
浮点错误

谢谢你的回复
7楼2010-03-16 11:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

气液两相流

铜虫 (小有名气)

引用回帖:
Originally posted by kuhailangyu at 2010-03-15 10:32:44:
这个问题已经有很多帖子了,不是很容易就能说清了就解决的

因为原因是多方面的,你可以搜一下原来的帖子,学习一下

字面意思大家都懂,呵呵

谢谢你的回复
8楼2010-03-16 11:07:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

气液两相流

铜虫 (小有名气)

引用回帖:
Originally posted by wgf332 at 2010-03-15 11:37:43:
1:可能是网格问题,调整网格,检查网格质量,基本上要控制在0.8一下。
2:重新初始化,改变初始化方式
3:试试双精度

什么是双精度?
9楼2010-03-19 09:59:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyinuniverse58

木虫 (正式写手)


nono2009(金币+1):鼓励应助! 2010-03-25 17:07
气液两相流(金币+1): 2010-03-26 08:33
引用回帖:
Originally posted by 气液两相流 at 2010-03-19 09:59:20:

什么是双精度?

双精度就是你在打开fluent时系统让你做的选择,不要直接点确定,选择2ddp就是2维双精度了,其他以此类推。
10楼2010-03-25 14:42:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 气液两相流 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 17:30 by 4FFAWE8HcgUD
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 17:22 by 4FFAWE8HcgUD
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 10/500 2026-08-29 16:49 by zICmwzsBXjbN
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见