Question

Fluent 非稳态计算时 内部压力震荡 难以收敛

我是计算对冲式止回阀内流场，并进行受力分析。已经对阀进行了简化，画出的是流道部分，由于是旋转体，因此只画了旋转对称面。如阀体CAD图，图中标1的是阀座，可以轴向移动，标2的是阀头，也可以轴向移动。当流体从左侧进入时，一部分流体通过冲压管进入活塞体（图中2），一部分通过阀头和阀体环形部分向左流。由于外部节流作用，内部的压力大，推动阀头向右运动关闭阀门。 在用非稳态计算时，我发现里面的压力一会高，一会低，出现震荡。（下图也显示了，第一张图是稳定的初始流场，内部压力高；但是迭代了几步后，就出现内部压力低，而造成了阀头向左运动出现了负体积）。划分的是三角形网格，用的smooth和remesh（网格应该没问题，通过给阀头和阀体固定速度，他们可以运动）。压力松弛因子是0.01，动量松弛因子是0.5，湍流动能松弛因子也是0.5。我设置的是定压差，压力入口100000，压力出口96000。udf中看阀头受力的语句是 begin_f_loop(f,t) { F_AREA(A,f,t); f3=f3+F_P(f,t)*A[1]; } end_f_loop(f,t)
阀头因受力运动，因为设置了deform边的缘故，只能向右运动，而向左运动会出现负体积。按正常情况，阀头受合力是向右的，稳态计算下都是向右。　　我的时间步长都到e-5了，还是会出现这种震荡（就是活塞体里面压力一会儿变红，一会儿变蓝）。听说调弹性因子可以改进，但是好像也没成功。不知道问题出在哪儿？怎么改进？ 您曾经帮我解决过问题，所以还是想请您帮助一下，谢谢啦！





Answer 1

你这个udf求的是y方向的力，而不是x方向的力

追问

不好意思，在这上面打错了，我UDF里写的是A[0]。别的还有什么问题吗，会不会是有什么因子我没有调什么的

追答

那我觉得你应该考虑水的可压缩性

追问

如果考虑可压缩性，那么该怎么求解？还有会不会是网格质量有问题，以及求解器或者单精度双精度什么的？
谢谢啦

追答

考虑压缩性需要写水密度随压力变化的udf
基本上可以排除你说的这些原因