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



  • 各位老师,我在做亚临界雷诺数下二维圆柱的涡激振动(y向),目前模拟的工况是Re=5000,圆柱直径0.0381m,速度0.132m/s,约化速度为8,质量比2.6,阻尼因子0.0036,振动方程通过编写四阶龙格库塔的UDF导入到fluent中,振动方程进行了无量纲化。
    微信截图_20200712212314.png
    模型如下,上游10D,下游20D,宽度上下各15D,圆柱附近采用四边形网格,其余为三角形网格。第一层网格高度0.0023D,时间步长0.0014s,湍流模型k-wSST。
    图片1.jpg
    在低雷诺数时基于编好的UDF进行过计算(此时的模型不一样,直径0.0016m),与实验数据接近,但在亚临界时,计算开始一段时间后圆柱位移一直增大,直到出现了负体积,此时的位移已接近10D。自己更改过时间步长,但是问题也是一样,想问下有老师做过类似的模拟吗?这个问题排查了很久还是没能解决,谢谢大家!
    编写的UDF如下。
    微信截图_20200712215323.png



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



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



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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2