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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 非稳态计算的一个问题
111/2返回列表12下一页
查看: 4275  |  回复: 29
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kylafree

至尊木虫 (知名作家)

[交流] 非稳态计算的一个问题 已有8人参与

非稳态计算时，
采用了一个稍小的时间步长计算，最后一直算到了稳态；
同样这个工况，采用稍大一点时间步长计算，最后也算到了稳态；
但两种情况算出来的结果却相差很大。比较困惑。哪个结果更可靠一点呢？
回复此楼

» 收录本帖的淘帖专辑推荐

fluent帖子收集 模拟相关

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-28 22:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xiegangmai(金币+1): 谢谢参与！ 2011-10-10 14:37:17
引用回帖:
7楼: Originally posted by kylafree at 2011-10-10 08:24:56:
还真没遇到过提示，不管采用多大的时间步长都没遇到过提示

不知道你具体的计算背景是什么，毕竟我不知道Fluent具体的源代码是啥，你选取的具体解算器是什么结构的。 但是从一般的经验上来看，时间步长小一些的精度高一些。
赞一下(2人) 回复此楼
8楼2011-10-10 09:44:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
臭水沟(金币+2): 谢谢交流~~ 2011-10-10 18:01:32
引用回帖:
9楼: Originally posted by kylafree at 2011-10-10 14:06:32:
简单介绍下背景： 一个三维圆管，一个进口，一个出口，壁面；
仅有流动和组分输运，其中管壁面对其中一种组分有吸收。

你是按URANS算的？
这样时间步长小一点好一些

还有一种可能是 因为你不是自己写的程序，而是用的fluent
由于使用了组分方程 所以本身方程实际上并没收敛……
ps：这是因为Fluent的解算器本身的限制，带组分方程的时候，求解是分离求解的（即便是耦合求解器，也是把组分方程啥的都扔外头），然后组分方程没解耦，然后有一些数值扰动一直在流场里面耗散不掉~~
当然，这是我的一种猜想

你收敛的判断依据是什么？
赞一下(3人) 回复此楼
10楼2011-10-10 17:19:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
11楼: Originally posted by kylafree at 2011-10-10 20:32:01:
组分方程残差降到e-6  连续性残差及速度的残差都降到了e-4以下了，并且流场变量的变化很小了， 不过没作网格独立性检验。

使用的是RANS么？
计算特征长度和特征速度分别多少？
时间空间格式选取怎么选的
赞一下(1人) 回复此楼
12楼2011-10-11 09:56:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
今天特地翻了一下J.Blazek的书，上面6.1.4明确说 对于定常问题，收敛后不影响计算精度的。 建议尽可能扩大时间步长。
赞一下(1人) 回复此楼
13楼2011-10-11 12:50:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
14楼: Originally posted by kylafree at 2011-10-12 10:35:10:
没有使用RANS，仅是层流模拟；
真正的模型非常不复杂且非常不规则，是一个动物的鼻腔内部，之前为简单说明问题说成了是一个三维圆柱。特征长度取0.001m吧，特征速度0.6m/s；
插值格式用的默认格式，其中空间格 ...

复杂外形流动的层流流场本来就不好算……
我没算错的话，按照你给参数对应的Re都成上万了……
用层流模型算有问题吧。。
赞一下(1人) 回复此楼
17楼2011-10-12 13:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
15楼: Originally posted by kylafree at 2011-10-12 10:40:52:
你的意思是指，既然我这个问题本身可以算到稳态，那么用非定常求解器求解时，时间步长可以随便取？

还是仅指要等到算到稳态了之后，时间步长便可以随便取呢？但这样好象就没什么意义了，反正都已经算了稳态了 ...

用非定常方程计算定常问题本来时间步长都是按照当地时间法取的
你随便找任何一本专业点的CFD教材都是这么说的……
事实上跟你算过程中除了最后完全收敛以外的中间任何结果都是伪的……
没任何物理意义

建议你直接导出整个流场的不同时刻 仔细比流场的所有细节…… 这样比较好

按照你给的Re，用层流模型算，这都方法步骤快成做DNS了……
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

fatsen
18楼2011-10-12 13:13:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by kylafree at 2011-10-12 11:01:31:
还有一个问题，按照稳定性要求推算出一个时间步长，但这个时间步长非常之小，比方时间步长为1e-5，如果要算1s的工况的话，就得算10 0000个时间步，即使每个时间步只迭代1次的话，总共也得迭代10 0000次。对于网格 ...

所以实际上大多数工程算法都是用隐式格式，采用双时间（引入伪时间）法来求解的
本来最早引入算法的时候就是求解的定常问题的解，然后有人说你这东西好是好，就是算不了非定常问题，结果这哥们就鼓捣出来了双时间法解决了这个问题……

做流动稳定性数值模拟的时候就是按照你说的用显格式R-K这么一点一点的推进算的，所以这东西才这么不好算。 我以前算一个算例需要20台以上的节点机算一个礼拜才能出结果，还是对非常简单的一个三维问题…… 恩 还是结构化网格。。
赞一下(1人) 回复此楼
19楼2011-10-12 13:17:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
20楼: Originally posted by kylafree at 2011-10-12 14:53:39:
Re数才几十啊，
取进气气体为空气，其密度为 1.225kg/m3 , 动力粘度为1.8e-5 ,
流速0.6 ， 特征长度0.001  则 Re=   1.225*0.6*0.001/1.8e-5 = 41

我以为介质是水了……
赞一下(1人) 回复此楼
23楼2011-10-12 17:38:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xiegangmai(金币+2): 谢谢参与！ 2011-10-12 20:52:14
引用回帖:
22楼: Originally posted by kylafree at 2011-10-12 15:33:03:
对阁下讲的前两段很感兴趣，能推荐一本CFD教材吗

你用Fluent 看Fluent的理论手册就行。 大多数求解问题都说的很清楚了，专业的CFD教材里面都是推公式…… 你入门更复杂 -_-#
入门教材一般p解算器的看陶文栓的那本数值传热学就行
        密度解算器的入门的可以看看anderson的计算流体入门与应用
        上次在书店貌似看到这本的中文译文版本了，但是不知道翻译的咋样
赞一下(2人) 回复此楼
24楼2011-10-12 17:42:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)
引用回帖:
21楼: Originally posted by kylafree at 2011-10-12 15:30:34:
fluent中默认的就是隐式格式，且只能选隐式格式，不能改变。
双时间法，挺高深的内容啊，现在不懂，得花时间看看...

这里的显式说的是时间层上的 Fluent这方面了解的不多，Fastran用的时候是支持R-K的显式格式的。
赞一下 回复此楼
25楼2011-10-12 17:43:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kylafree 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭