请教各位使用过，浸入边界法IBM，解决流固耦合fsi问题的大佬。

最近在用foam-extend中的浸入边界法，算圆柱绕流（刚性问题） 如果我想考虑圆柱的流固耦合fsi，请问设置上的区别在哪呢？
即：浸入边界法，做流固耦合问题，和刚性问题，设置上区别在哪呢？
万分感谢！

7ffe581e-acb0-42ea-abb7-64507e4652c1-image.png

@oitocfd
您好，请问foam-extend，浸入边界法，做流固耦合问题，和刚性问题，
设置上区别在哪呢?

@micro 自挖自填。
大家可以参考东岳老师关于icoFoam解析的文章，说的很正确的，一个时间步骤动量方程的系数是不变的，因为动量方程每个时间步中仅被组建了一次。但是每个时间步中，U可是被更新了好几次（在压力泊松循环中），两者并不矛盾！另外，我们同样可知，一个时间步中，我们反复求解压力泊松方程（一般是3次），压力泊松方程的系数同样是不变的，变化的是右端项（因为右端项是用速度组建的）。

在设置探针的时候，有个fixedLocation的条目可以设置为false，这样的话对于动网格探针应该是会随着网格拓扑更新移动的。

@xpqiu 谢谢您的回复！关于2，我设置的phi的tolerance是10^-6，solver是GAMG。把nNonOrthogonalCorrectors设为10之后，最后几次计算的No. Iteration 就随计算次数减小了，最后为0，我想应该是残差已经降到最低。谢谢！

@cresendo 那就是纯对流方程了

@xpqiu 收到，还是非常感谢！！解释的已经很清晰了！

@五好青年 tutorials里面的wingMotion

目前也不确定是openMPI的问题还是OpenFOAM太挑剔，所以各位使用集群计算的老师同学们，能不能看看你们集群的hostname是不是有下划线的，能让我定位一下问题的位置。:xiezuoye:

@李东岳 了解了，谢谢东岳老师

实在没空弄了 我看看有没有人给你弄 https://www.cfd-china.com/topic/4861/领取任务免费赠送cfd中文网注册资格

@wying 要非常高的精度，比如普方法、有限差分那种7 8 9阶，有限体积法大体都不行

如题，根据probes的字母意思解释，是Set of locations to sample（设置位置集），我看了下probes.H和probes.C的源程序，没有看懂，probes.C的部分代码如下：

const label celli = mesh.findCell(location); elementList_[probei] = celli; if (celli != -1) { const labelList& cellFaces = mesh.cells()[celli]; const vector& cellCentre = mesh.cellCentres()[celli]; scalar minDistance = great; label minFaceID = -1; forAll(cellFaces, i) { label facei = cellFaces[i]; vector dist = mesh.faceCentres()[facei] - cellCentre;//是指面心和体心的插值么？ if (mag(dist) < minDistance) { minDistance = mag(dist); minFaceID = facei; } } faceList_[probei] = minFaceID; }

@李东岳 是这样啊，我还以为是软件系统文件缺少呢，感谢感谢

在OpenFOAM文件家下打开终端输入

find src applications -name "*.L" -type f | xargs sed -i -e 's=\(YY\_FLEX\_SUBMINOR\_VERSION\)=YY_FLEX_MINOR_VERSION < 6 \&\& \1='

@cfdngu 好的，非常感谢您快速的解答，我再尝试尝试。

@samuel-tu 主要这文章看看也很有帮助，它的角度偏代码整体架构而不是力学，对理解openfoam的设计还挺好的。

@小狗狗 谢谢您，我学习一下

@Yongbo :high: 谢谢

@bestucan 后来我也是这么解决的。。相当于直接安了新的Paraview