24小时热门版块排行榜    

查看: 4186  |  回复: 29

kylafree

至尊木虫 (知名作家)

组分方程残差降到e-6  连续性残差及速度的残差都降到了e-4以下了,并且流场变量的变化很小了, 不过没作网格独立性检验。
11楼2011-10-10 20:32:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shookware

铜虫 (初入文坛)


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

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

shookware

铜虫 (初入文坛)


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

kylafree

至尊木虫 (知名作家)

引用回帖:
12楼: Originally posted by shookware at 2011-10-11 09:56:59:
使用的是RANS么?
计算特征长度和特征速度分别多少?
时间空间格式选取怎么选的

没有使用RANS,仅是层流模拟;
真正的模型非常不复杂且非常不规则,是一个动物的鼻腔内部,之前为简单说明问题说成了是一个三维圆柱。特征长度取0.001m吧,特征速度0.6m/s;
插值格式用的默认格式,其中空间格式用二阶格式算过,但差别不大。
14楼2011-10-12 10:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

引用回帖:
13楼: Originally posted by shookware at 2011-10-11 12:50:58:
今天特地翻了一下J.Blazek的书,上面6.1.4明确说 对于定常问题,收敛后不影响计算精度的。 建议尽可能扩大时间步长。

你的意思是指,既然我这个问题本身可以算到稳态,那么用非定常求解器求解时,时间步长可以随便取?

还是仅指要等到算到稳态了之后,时间步长便可以随便取呢?但这样好象就没什么意义了,反正都已经算了稳态了 再算下去也还是这个样子。
15楼2011-10-12 10:40:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

还有一个问题,按照稳定性要求推算出一个时间步长,但这个时间步长非常之小,比方时间步长为1e-5,如果要算1s的工况的话,就得算10 0000个时间步,即使每个时间步只迭代1次的话,总共也得迭代10 0000次。对于网格数量比较多,而电脑又不给力的人来说,这岂不相当于说根本就做不了非稳态模型了吗?
16楼2011-10-12 11:01:31
已阅   回复此楼   关注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都成上万了……
用层流模型算有问题吧。。
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了……

» 本帖已获得的红花(最新10朵)

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台以上的节点机算一个礼拜才能出结果,还是对非常简单的一个三维问题…… 恩 还是结构化网格。。
19楼2011-10-12 13:17:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

引用回帖:
17楼: Originally posted by shookware at 2011-10-12 13:11:01:
复杂外形流动的层流流场本来就不好算……
我没算错的话,按照你给参数对应的Re都成上万了……
用层流模型算有问题吧。。

Re数才几十啊,
取进气气体为空气,其密度为 1.225kg/m3 , 动力粘度为1.8e-5 ,
流速0.6 , 特征长度0.001  则 Re=   1.225*0.6*0.001/1.8e-5 = 41
20楼2011-10-12 14:53:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kylafree 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见