如何解决OpenFoam升阻力系数剧烈震荡的问题？
-
大家好，我在使用foam-extend4.1的pimpleDyMIbFoam求解柔性翼型波动的算例时，发现计算出来的升阻力系数剧烈震荡，加密网格，缩小时间步长，更改离散格式（试了一部分），更改fvsolution里的设置等方法都试过了，发现算出来的力系数依旧存在震荡，有大佬知道怎么解决吗？
是不是应该要换一个求解器呢？我打算将pimpleDyMIbFoam改成icoDyMIbFoam用来求解看看会不会有改善。
现在我对比pimpleDyMIbFoam和icoDyMFoam，把pimpleDyMIbFoam涉及浸没边界部分的程序移植到icoDyMFoam中，做一个新的求解器icoDyMIbFoam。但是当我使用icoDyMIbFoam求解时会出现pcorr求解奇异的问题。心累
请问各位大佬有什么好的解决方案吗？
这个是我用icoDyMIbFoam计算fe41里的某个算例的报错
-
fe4.1-IB求解器，印象中测试圆柱绕流的力系数也存在较大震荡，不知道是不是求解器没做好。可以试试fe-5.0会不会好用一点
-
@luofq-sysu fe5.0已经试过了，还是会存在震荡
-
@GrpOne123 在 如何解决OpenFoam升阻力系数剧烈震荡的问题？ 中说：
是不是应该要换一个求解器呢？我打算将pimpleDyMIbFoam改成icoDyMIbFoam用来求解看看会不会有改善
请教一下：这边移植到icoDyMlbFoam的用意是什么呢？因为不加湍流模型嘛
-
@luofq-sysu 是的，我这个算例Re=100，用的是层流，想看看用icoFoam能不能算出来。关于pcorr方程奇异的问题，你有遇到过吗？
-
@GrpOne123 不好意思，我是很久之前测试的，只记得容易出bug，然后力也似乎算不准。然后就没继续了
刚看了一下，我当时用pimpleDyMIb算了一个旋转圆柱类型的工况，是能并行的。pcor方程的问题没有遇到过或者没印象了。 -
@luofq-sysu 哦哦哦，好的谢谢
-
在foam-extend4.1的ibm里折腾了两三个月，发现用foam-extend4.1的浸没边界法计算升阻力系数产生震荡是其本身的问题。具体原因可以参考这篇论文Sources of spurious force oscillations from an immersed boundary
method for moving-body problems
这篇论文提到了三个解决方法：1、在力源项后面再增加一个质量源项以抑制动边界的力震荡 2、增大时间步长 3、将动边界所处位置的网格加密
目前我只是用了后面两种方法可以有效的减少其震荡，但震荡还是会存在，如果想获得一条平滑曲线，可能需要与第一种方法结合。第一种方法涉及到一些算法问题，鉴于目前的水平有限还不会在openfoam中实现，希望有大佬能够实现并给大伙解惑。