@李东岳 对的，大概有十几公里，那谢谢东岳老师，我再使劲努努力钻研一下。个人觉得根据文献植入算法太难了，虽然能看懂离散和求解过程，但是想在openfoam里面实现觉得难度挺大的，当时老师您带我们植入的单相可压缩的，目前这是两相可压缩，我努力努力，再次感谢东岳老师的持续关心:xinxin:

问题已解决，勾选paraView中的“Use VYkPolyhedron”选项后交错就会消失

@bestucan 我看了目前用到这个 -I. 的地方，包括reactingParcelFoam跟sprayFoam，都有子文件夹，比如sprayFoam下的egineFoam，应该跟compile子求解器有关

@李东岳 东岳老师是想诱惑我再报一次您的辅导班:141:

@cresendo 我也是看了亿遍~~~~~

在opemfoam中有一个内聚力，可以使粒子之间相互黏结。它是通过粒子之间的碰撞的能量转换实现的。其中通过一个叫cohesionEnergyDensity的系数来决定它的强度。这个密度越大，黏结越明显。希望我说的这个是你想问的。

@史浩 VOF和LES都是高度网格依赖类求解器。主要是LES只能捕获大于网格尺度的涡，因此，网格越小，LES越能捕获到更小的涡。类似的，VOF捕获到的界面只能大于网格尺度，如果计算域存在气泡，如果网格越小，VOF可以捕获到的气泡越小。

@李东岳 实在是有点困难~请问老师如何让两个场（如图片中的a和vDotAlpha）对应网格的值相乘，这不是内积也不是外积，用forAll循环的话，会报assignment of read-only location错误。

forAll(a,celli) { vDotcAlphal[celli] = a[celli]*vDotcAlphal0[celli]; vDotvAlphal[celli] = a[celli]*vDotvAlphal0[celli]; }

图片2.png

麻烦帮忙改正一下。想要规定thetaBoundary为一个常数值，然后对theta（相分数）进行对比，总是实现不了。谢谢。

const Field<scalar>& theta();
const Field<scalar>& thetaBoundary();
thetaBoundary = {0.5};
if (theta < thetaBoundary)

算出来的等值面结果如图
2bfde5ee-55f9-497e-9109-0ab0d48a0789-image.png

Hi，东岳大佬，请教下openfoam里可以实现分区网格计算，比如两个分区，并设置interface boundary condition来coupling这两个区域吗？如果可以，麻烦给点思路，谢谢。
我暂时想到的只有fluid-solid interaction，即FSI toolbox in Foam-extend，改写solid domain为fluid-fluid interaction，来实现，不过有些繁琐了。

BR，
Shiyu

@李东岳 谢谢东岳老师，我还是在win下安装paraview吧~

@lys 用的已有模型，专门针对非牛顿流体的曳力模型还没用过，希望能进一步交流啊

有用过Foam-extend 4.0 里的interIbFoam的吗？这个求解器是不是不能加湍流模型啊？一加湍流模型边界条件就报错

谢谢李老师，确实是alphatWallFunction出了问题。重新修改了这个函数后，就可以了。再次感谢。:xinxin:

@zhe 因为最后一步你重新给p.U()速度了啊。。

@荷兰后腰阿贾克斯 嗯嗯，谢谢老哥了:xinxin3:

Thanks a lot. I wonder if the partial elimination algorithm mentioned in Passalaqcua's paper was implemented in the current version?

@我是河滩 在 指定cellzone区域添加动量源项的问题 中说：

遍历了整个网格区域进行搜索

前辈您好，您是如何遍历整个网格区域并获取网格坐标的呢，是利用mesh.C().component()语句吗:141:

@李东岳 感谢老师关心！已经在家啦，趁着有空补补基础hhh