亚临界雷诺数下二维圆柱涡激振动出现负体积

[ycx1042341217]

各位老师，我在做亚临界雷诺数下二维圆柱的涡激振动（y向），目前模拟的工况是Re=5000，圆柱直径0.0381m，速度0.132m/s，约化速度为8，质量比2.6，阻尼因子0.0036，振动方程通过编写四阶龙格库塔的UDF导入到fluent中，振动方程进行了无量纲化。

模型如下，上游10D，下游20D，宽度上下各15D，圆柱附近采用四边形网格，其余为三角形网格。第一层网格高度0.0023D，时间步长0.0014s，湍流模型k-wSST。

在低雷诺数时基于编好的UDF进行过计算（此时的模型不一样，直径0.0016m），与实验数据接近，但在亚临界时，计算开始一段时间后圆柱位移一直增大，直到出现了负体积，此时的位移已接近10D。自己更改过时间步长，但是问题也是一样，想问下有老师做过类似的模拟吗？这个问题排查了很久还是没能解决，谢谢大家！
编写的UDF如下。

[yxg2018]

你的动网格更新怎么设置的

[ycx1042341217]

@yxg2018 我现在动网格是采用弹簧更新和网格重构，将圆柱边界设置为刚体运动

[yxg2018]

@ycx1042341217 网格重构应该是只更新三角形网格，你可以将圆柱和边界层网格一起做刚体运动

[ycx1042341217]

@yxg2018 谢谢，我将网格按照单元类型分离并将边界层网格也设定为刚体运动区域，计算依旧会报错，提示Divergence detected in AMG solver: x-momentum，我将迭代次数增加并改小了动量的松弛因子，计算还是发散。