@Imech 在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：

@学流体的小明 老哥，这个问题你解决了吗，我最近算完美滑移边界发现算的结果与参考对不上

最近还在搞，你是什么问题？算例是什么？

@Zhy2022 您好，frictionalStressModel和granularPressureModel里面的模型我都试了一下,并且有的跟fluent选的一样，但还是这个问题。我能否把算例分享给您，你帮我看一下，我的模型网格是20倍粒径，单核跑起来很快，不知道怎么分享给您合适？还有一个问题就是，我如果不把物性参数改成water，用官方算例air，结果就是正常的，能达到紧密堆积0.6.不确定是不是我物性参数改的有问题？faac38a7-acef-406c-8ae4-47b7d5bbfd64-image.png

@队长别开枪 感谢大佬回复，这就去学习一下！

compressibleInterFoam自带的tutorial里面就有

为什么在controlDict中把writeInterval调大，计算就发散？？？为什么会有这么诡异的bug
writeInterval不就是控制写入频率的么，跟计算有关系吗？
有没有大佬遇见过
7f83b6e1550cdbc87dc895138ddbea5.png 1748766340208.png

@队长别开枪 好的，感谢！

@李东岳 好的李老师，搞定了

@WangGH 直接wmake就行，在那个C++代码所在的文件夹中

感觉这个很流弊，我一直没尝试。之前让学生研究过步骤，也是半途而废了。

感谢李老师指点！

还有怎么在两个分开的流体域之间建立cyclic边界条件，谢谢！

@上级 我想到一个特别的方法，writeCellCentres输出网格中心坐标，然后写个Python/Matlab等脚本，判断一下中心点在不在你想要的范围内，然后给赋值
（感觉听起来难，操作起来应该挺简单）

@李东岳 :zoule: :ok:

@heyi999 feeling_coolabc@163.com

@李东岳 谢谢东岳老师

@李东岳好的，谢谢老师！

@biaowai 是的，所以怀疑是不是设置出了一些问题

@fan820 上述问题已解决只需要在''Surfaces to Plot''处右键选择Boundary cell faces，内部面就可以显示出来:146:

@学流体的小明 在 setFields是否可以并行计算提高效率 中说：

确实很慢。我想手动划分计算域，使用setFields也很慢。看别人的思路，先topoSet标记区域，然后再对区域进行setFields，也感觉挺慢的，时间上差不多吧？但是topoSet之后可以重复使用这个被标记的区域。
另一个思路是自己写程序进行赋值。我推测setFields和topoSet很慢，是因为每个赋值块它都要遍历一次网格，设置几千个球状的体积分数，就要遍历几千次网格？自己写程序的话，遍历一次就够了。

感谢您提到的思路，我还没有尝试过。之所以提到这个问题，一方面是因为计算速度，另一方面是因为一次性setFields几千个球体经常会在中间报错（各种错误），导致我一直失败。昨天晚上，我尝试了将几千个球体分成四份setFieldDict文件，大概每次set 1k多球体，这样可以保证不报错。实际体验来看，1s可以设置1个球体，4k个球体的体量大概1h左右就可以setFields完成，不报错的话就可以接受这样的速度。

@李东岳 谢谢李老师！