流固耦合提高收敛性的经验总结

cainiao07

这两天做一个湍流模型与岩石的流固耦合，多次碰到不收敛，且问题为ADF 8000。
总结如下 ：

1- 报错显示能量方程计算不成功
   这有很多原因引起，包括边界条件、材料参数、模型的单位以及迭代的时间步不合理 。。。。。具体要根据你的模型去检查。   我把所有的都式过，最后发现有一个区域的热流通量存在极值，最后发现这个区域单元划分时有个边没有划分单元分数。

2- 报错显示单元重叠element overlapped
     这个后来发现是因为选用的FIBC单元+heat transfer组合引起的，查阅英文手册，发现这个单元在解决大模型（我的单元数是10万）时计算是不推荐的，ADINA推荐的是 FBIC-C单元，说这个单元室分离法求解耦合问题，更适合大模型。

3-警告 “Model may be unstable, ratio of diagonals >1.E11, please check your input data”，计算不到100步就终止。

开始以为是固相模型有问题，毕竟网页上关于固体力学求解碰到这个问题多。 后来我反复推敲固体模型，确实不存在问题。检查流体模型，一度怀疑是不是材料参数有问题，后来也不是。最后通过改变“载荷步+初始条件”，发现可以提高收敛性，效果很好。
   根据我的反复试用，小的载荷步是必须的，然后在初始条件中，尤其要注意调整 k, e （湍流是k-e模型）的值，这两个的初值对可压缩流体的湍流模型求解收敛十分重要。


谢谢，写的不好，请见谅。

弦犹在耳

支持原创高质量经验贴~

cik

学习了，谢谢楼主