CFD中文网

大家好，我将文献中的公式植入到foam—extend4.1的force.C文件中发现算出来的数据与文献的数据相差较大，下面是我的代码和公式，希望有大佬帮我看一下我编的是否正确。ps：我用的是浸没边界法，不知道这个速度这样获取可以不。
be1e4612d51ceea1f453b513f11f757.png
1.png
2.png ```

在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中实现，希望有大佬能够实现并给大伙解惑。