找回密码
 注册
查看: 2193|回复: 8

向高手请教一个发散的问题

[复制链接]
发表于 2010-1-5 20:03:09 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
计算的是一个三维火车的模型,采用耦合隐式求解,残差图如下,想问一下有哪位能看出来发散原因是什么?或者可能的原因!谢谢!希望大家多给建议!

残差图

残差图
 楼主| 发表于 2010-1-5 20:54:52 | 显示全部楼层

回复 1# heal2005 的帖子

怎么没人说话,没人知道吗!提供点建议也可以。
发表于 2010-1-5 21:15:38 | 显示全部楼层
你看了流场了吗? 不一定是发散。能不能接着往下算?
 楼主| 发表于 2010-1-5 22:10:04 | 显示全部楼层

回复 3# 通流 的帖子

接着算和这差不多,下降后又继续上升,没算完,我看了压强场太大了,应该是错了!
发表于 2010-1-6 15:45:40 | 显示全部楼层
调一下能量的收敛因子,调到0.9试试
发表于 2010-1-6 16:53:45 | 显示全部楼层
对于隐式格式求解,我是不太懂。我瞎猜一下,可能是网格有问题。

如果,流动是定常的,隐式格式应该是很好的,因为时间步长可以取得很大。但如果是非定常的,那么大时间步长的后果,我搞不清。我用CFD一般先用pressure based(应该就是SIMPLE的算法)。这东西收敛快。如果是高速流,当然的用density based(时间推进)做。如果开始时收敛不了,就先改用非定常算一会。基本的规律是,你的数值模拟越接近实际情况,就越稳定可靠。实际流动都是非定常的,非定常计算也最保险。代价就是计算时间。不过现在计算已经变得便宜了,就让它先算一晚上也没什么吗。
 楼主| 发表于 2010-1-6 18:55:58 | 显示全部楼层

回复 6# 通流 的帖子

谢谢版主的分析,那我就试一下你的方法,我这个是高速流,采用的是定常可压模型,我也怀疑是网格问题,所以把网格细化了一下,效果很好,可以收敛了!不过结果不是非常好,所以会采用你的建议再试一下!
发表于 2010-1-7 21:56:46 | 显示全部楼层
版主讲解的很深呢,学习了。另外,想问一下,如果定常求解器求解非定常的问题,是不是也容易出现上面的情况。谢了
发表于 2010-1-9 11:10:31 | 显示全部楼层
有时候,残差降不下去的原因就是局部的流动是非定常的,所以显得不收敛。这时候,就靠用CFD的人的判断了。如果认为局部的问题不会影响整个答案,那么就可以不理局部的问题。如果,认为局部的问题很重要,那么,就得做非定常计算。大部分的计算,局部的问题确实不影响大局。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表