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

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如下。