CFD中文网

我之前使用相同网格，
fsiProperties:

couplingReuse 1;
couplingScheme Aitken;

时非常好用。
不知道为什么改成了:
fsiProperties:

couplingReuse 1;
couplingScheme IQN-ILS;

就出现了以上错误。

我在使用foam-extend-4.0与第三方软件fsiFoam进行流固耦合并行计算，突然就报错了：

Time = 0.1067, iteration: 2
Maximal accumulated displacement of interface points: 0.122364
GAMG:  Solving for cellMotionUx, Initial residual = 0.992404, Final residual = 0.000953657, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0.999421, Final residual = 0.000628662, No Iterations 1
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] Some processors detect different empty (2-D) directions.  Probably using empty patches on a bad parallel decomposition.
Please check your geometry and empty patches
[1]
[1]     From function void polyMesh::calcDirections() const
[1]     in file meshes/polyMesh/polyMesh.C at line 101.
[1]
FOAM parallel run aborting
[1]
 --------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

跪求大佬指点迷津！

@cuimai 个人感觉solids4Foam有点不稳定，我用它自带算例也跑不通。cfdonline也有很多人反应有bug

@ice_flow 您好！在starccm中创建切割体网格如图所示。

@ice_flow 您好！个人感觉这点误差是可以接受的。还有向您请教一下，openfoam识别多面体网格么？是否可以在paraview中展现多面体网格

跟进一下，这是不同边界条件下的压力云图。可以直观的看到入口处压力场出现震荡。

我也尝试从数值角度说一下我的看法。计算求解的是NS方程，对于入口处速度和压力均为定值的初始条件，可能会导致结果不满足NS方程。造成所谓的过度约束。而入口出口分别指定速度和压力为固定值，则会避免过度约束的现象发生。

欢迎各位大佬批评指正！

@李东岳 谢谢李老师！

圆柱绕流是各种教程中最常见的算例。其中对于出口和入口的边界条件常见为：入口（速度指定值，压力zeroGradient），出口（速度zeroGradient，压力指定值p=0）

但为什么这样呢？将入口速度和压力均设为指定值（U=1，p=0）,出口速度和压力均设为零梯度，这样会有什么不妥呢？

我看资料有说入口这样的边界条件会过定义，导致方程无法求解。但具体原因和细节不甚了解。

请大佬赐教！（最好可以从数值计算和使用经验的角度讲解）提前拜谢了！

@Mina-Lee 不是。类似于channel的边界不是刚性而是弹性的，channel中有刚体运动（计划用重叠网格实现）。

@Mina-Lee 谢谢大佬！

顶，欢迎更多大佬指点迷津

@HITSC30 谢谢

各位大佬好！小弟想研究流固耦合问题。计划用OpenFOAM算流体，preCICE做媒介，结构求解器哪个好一些？calculix还是solids4Foam。solids4Foam自带的流固耦合我品鉴过了，在fe41基础上添加的solids4Foam。感觉流固耦合不如fsiFoam计算效率高。calculix没有涉及。跪求大佬指点迷津，建议从上手难度和稳定性的角度给予更合适的建议。提前拜谢了！

@yuan您好！我觉得是坐标值而非绝对值。

@李东岳 好的，李老师

@李东岳 李老师您好！这个链接失效了，麻烦补一下，谢谢

cfd-online也有人提出了相同的问题，但Jasak表示没有问题，或许真的是安装版本的问题，请问用fe40的兄弟们有遇到相同问题吗？

如果可以正常运行，可以分享一下fe40安装包吗？

@李东岳 麻了，哭。。。我多么希望是我自己的问题

@李东岳 我要用到基于fe40的fsiFoam。这么说我的代码没有问题么？

如图所示，第一个时间步网格会出现些许运动，之后就停止了。