CFD中文网

我最近也碰到这个问题，单核计算时对某个特定的patch上的面心可以赋值，但并行分块之后也是运行不了。

@mollychen95 请问你解决这个问题了吗？我也遇到只有一条edge运动的情况，似乎动网格没有扩散出去

@wh3296 请问您解决了这个问题了吗？我现在想要获取整个patch的位移，这是我的代码，但是计算出来的位移是0

        forAllConstIter(labelHashSet, patchSet_, iter)
        {
            label patchi = iter.key();
	    vectorField displacement = mesh.Cf().boundaryField()[patchi] - mesh.Cf().oldTime().boundaryField()[patchi];
	    vectorField Uw = displacement/t.deltaT().value();
            vectorField pf = Sfb[patchi]*(p.boundaryField()[patchi] - pRef);
            
            vectorField vf = Sfb[patchi] & devRhoReffb[patchi];                
            fm.first().first() += rho(p)*sum(pf);
            
            fm.first().second() += sum(vf);
            
            vectorField CPs = (((pf + vf)&IE_) * (Uw&IE_))*IE_;

            fm.second().first() += sum(CPs);

            fm.second().second() = IE_;
        }

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

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


```

@wuyukai 请问你这个问题解决了吗？我遇到同样的问题

@luofq-sysu 哦哦哦，好的谢谢

@luofq-sysu 是的，我这个算例Re=100，用的是层流，想看看用icoFoam能不能算出来。关于pcorr方程奇异的问题，你有遇到过吗？

@luofq-sysu fe5.0已经试过了，还是会存在震荡

大家好，我在使用foam-extend4.1的pimpleDyMIbFoam求解柔性翼型波动的算例时，发现计算出来的升阻力系数剧烈震荡，加密网格，缩小时间步长，更改离散格式（试了一部分），更改fvsolution里的设置等方法都试过了，发现算出来的力系数依旧存在震荡，有大佬知道怎么解决吗？
是不是应该要换一个求解器呢？我打算将pimpleDyMIbFoam改成icoDyMIbFoam用来求解看看会不会有改善。
现在我对比pimpleDyMIbFoam和icoDyMFoam，把pimpleDyMIbFoam涉及浸没边界部分的程序移植到icoDyMFoam中，做一个新的求解器icoDyMIbFoam。但是当我使用icoDyMIbFoam求解时会出现pcorr求解奇异的问题。心累
请问各位大佬有什么好的解决方案吗？

这个是我用icoDyMIbFoam计算fe41里的某个算例的报错