@上级

我在想能不能实现一个边界条件，先判断坡面上哪里有泥石流（边界面相分数判断），然后把边界面有相分数的的网格面提取出来，对这些网格面改变的值。

可以，就是要符合物理，这个模型估计要琢磨下

@五好青年

单核可以计算，并行报错

这个会发生，但是几率非常小，我搞openfoam这么多年，就发生过一次，我还去问openfoam官方了，他们也不知道咋回事。debug不太容易。非常的case dependent

@tens 在 OpenFOAM如何设置球形相域？ 中说：

@郑学习

球形用sphereToCell 在openfoam下\etc\caseDicts\preProcessing\setFieldsDict可以查到

感谢！跑通了！

@zhouxu 好的，谢谢，我尝试一下，因为刚使用对每个方程的计算不太了解

@cccrrryyy 好的，谢谢您。

@tens :136: 啊这，谢谢前辈:146:

@dy大世界 @cfdngu
蹲一个回复，关注2个问题
(1) 单核可以计算，并行报错
(2) 单核和并行计算，结果存在差距，差距是细微的可以忽略，还是跟核数有关？

@cfdngu
蹲一个回复，关注2个问题
(1) 单核可以计算，并行报错
(2) 单核和并行计算，结果存在差距，差距是细微的可以忽略，还是跟核数有关？

@白礼耕 我也遇到这个情况，只在特定网格，特定时间步下报错
换个情况可能就正常计算了
你最后知道原因了吗

@李东岳 是的，我apply固体域：
5ea19143-3aec-4a34-ab54-3761c2fd64bc-image.png
但是只有颜色变化：
f8539f48-419c-470a-9d9a-c14b025a53e9-image.png

我问了一下其他用fsiFoam的人，他是用tecplot打开的，也是只有流体域能变形。

@五好青年 是自己改的程序吗，我之前好像遇到过。
可能是某个变量在计算过程中随着时间变化但没有保存，重新开始计算后这个变量被初始化，初始化的值和当前时刻的差很多，导致发散

@bestucan 好的，谢谢老师，我试试

借鉴这个教程已解决
https://zhuanlan.zhihu.com/p/450526288

各位老师与大神们好，我想再补充一下我的问题。
首先对于model_.valid() ，BlendedInterfacialModel.H中的注释为“//- Model for region with no obvious dispersed phase”，我浅薄的理解为alpha_1=alpha_2=0.5；
而对于subtract，我初步理解为subtract=f1-f2，subtract=0时为false，非0时为true；
但按照以上两种理解，我不明白对于计算f2时subtract为何要在原f2上进行变号，也无法推导出类似M=f1M1+f2M2的blending公式形式。
我认为一定是我的理解存在问题，想请各位指点，非常感谢~:xiexie:

@乌云盖雪 复制粘贴？

@bestucan 非常感谢！

@星星星星晴 太感谢了！:xinxin2:

@tens 已成功，谢谢

@xpqiu 最后运维还是换成了OF5，好像还连外网了，要是当初就说能连外网，我觉得自己也能搞定，何必多花这一千块，前前后后折腾了一个月。

@bestucan 实际当然是有体积的，概化成无体积就是想简单点，网格逐列关闭也可以的。

实在不行就动网格

我先试试楼上李东岳老师的方法