CFD中文网

我的猜测（不一定正确）：
在初始化的时候，seffield和toposet，已经指定了重叠网格的前景网格和背景网格，相当于背景网格指定了zone的ID
如果在计算过程中，用SnappyHexMesh细化了网格，那么此时的重叠网格zoneID就和初始时间步的zoneID不一致了，出现了找不到的情况。
所以，重叠网格应该是不支持在计算过程中使用SnappyHexMesh吧。。

@Joker 不好意思，我没用过这个ex.....。iostream 是负责读取信息和输出信息的，从你的报错信息来看，似乎是文件格式存在问题，比如，可能你这个需要指定额外的参数，但是你没给定。

@Joker 问题应该不是在于扩散方案的选择，而是在于你的**“设置”**错误。扩散方案选的不好，顶多会导致动网格拉伸严重。。。不会出现只有一层网格运动的情况。

@Joker 这应该是你的设置问题，我测试过这个fixedpoint，没问题的，这个功能只是让一部分区域的扩散系数设置为0，你这个应该是动网格设置有问题。

https://www.cfd-online.com/Forums/openfoam-meshing/134902-how-combine-dynamic-mesh-static-mesh.html 给你个好东东

@Gunther 嗷嗷，那可以做，祝好~

@Gunther 你要用precice做什么案例？不会是浮式风机吧。

推荐你使用calculix或者dealii。

calculix其实就是低版本的abaqus，支持abaqus的inp文件，只不过关键词少几个，并且有前处理器prepomax可以使用，油管上有很多视频可以自学，很多关键词可以参考abaqus的技术文件。

preCICE和solids4foam的区别可以参考下图：

稳定性方面，不太清楚solid4foam的算法，可以明确的是，这个软件是基于有限体积法的。有限体积法求解结构动力学表现的咋样不太清楚，但是应该是比较适合多物理场耦合。

还可以考虑GPU运行，不使用cpu运行，具体可以看一下这个WX公众号的内容：TFluid多物理场仿真

@尚善若水 可能不太行，我这个所有三方库、所有软件都用新C++重新编译的。

我之前在超算里安装precice+ccx+of，安装了一个月才装上，同样也是没有root权限，找了一大堆人，花了七八张毛爷爷才搞上。
我遇到的问题，OF编译就会报错，超算里所有的C++都试过了，都报错。最后找人查到了原因，超算的C++缺mpicxx环境，重新装的C++才搞上。

@msm 不好意思，根据你所提供的网格信息，无法判断是哪里的问题。
你提供的图片没看出来是哪个位置的。

1、可以参考fluent一样，把边界层提取出来，给边界层也指定位移。
2、也可以使用二次逆距离扩散系数试一下，你现在用的是一次逆距离扩散系数。
3、你需要关注一下，是不是因为单个时间步长的位移过大引起的网格问题，如果不是的话，可以调整动网格方法。

你应该放几个动网格的运动图片上来，并且描述一下你的动网格方法。

@李东岳 DEM软件是没有网格的概念的，因为DEM软件就是对每个颗粒求解牛顿第二定律，计算得到每个颗粒的运动。

所以OF耦合liggghts，OF是提取liggghts颗粒的位置、速度......，然后OF识别颗粒所在网格，最后计算拖曳力、体积分数......。反之，根据OF传递进来对每个颗粒的力，liggghts计算牛顿第二定律，得到颗粒位置。

换而言之，CFD-DEM耦合，大多数的活都需要CFD软件干。

所以我觉得，要做这个耦合，OF的代码修改量可能会远远大于liggghts的代码修改量。

@capillaryFix 我之前想干的事，是把openfoam和liggghts通过precice耦合起来，但是这个工作量太大了，一个人很难完成。
现在的OF和liggghts是通过侵入式的耦合，导致了升级OF版本，就需要重新修改求解器。

这是我之前的设想，这个设想，不需要修改liggghts源代码，但是，因为现在苦于毕业没文章，所以这个想法只能暂时搁置了。

您好，请问您这个问题解决了吗？能分享一下解决办法吗？

您好，请问您这个问题解决了吗？能分享一下解决办法吗？

https://github.com/louisgag/openFoam-Overset-SpeedUp

@RolandLent precice中没有计算OF力的能力，力都是在OF里计算的。OF中的力也是通过force函数计算的。
pecice只干三件事:
1、引用of中的力变量，提取流固耦合交界面的力变量；引用位移变量。2、precice内部运算。3、覆盖引用的位移变量。

因为你编译的openfoam-adapter里，已经对流固耦合交界面的force变量引用了，OF在运行完成后，adapter会引用force变量，precice运行完毕后，会覆盖这个force变量。
precice就是给openfoam施加了一个狄利克雷位移边界条件-理解了这句话就全明白了。