在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

最近对两个概念非常混淆，即摩擦系数和切线回弹系数，这两个系数是否就是同一个概念？

比如在以下的颗粒pairCollision参数设置中，法向碰撞的回弹系数根据alpha设定。但是切向的回弹系数，是不是就是这个mu?

pairSpringSliderDashpotCoeffs { useEquivalentSize no; alpha 0.01; b 1.5; mu 0.1; cohesionEnergyDensity 0; collisionResolutionSteps 1; };

@波流力 上传了所有案例文件，需要的请自行下载。见桥梁水动力：

https://mp.weixin.qq.com/s?__biz=MzU5NTgwODcyMg==&mid=2247483815&idx=1&sn=a0ef629459736bde6b94ebf0d31f68b1&chksm=fe6d191ac91a900cd2a2c2d45a6aafb711da8aab3f3fc7ad7b1d40969df1cc1f0f741080bc10&token=1944678271&lang=zh_CN#rd

压力速度迭代感觉国际上这面玩的少了 就刚出那时候搞出来一些SIMPLEC之类的

@同学博 好的 十分感谢！

各位朋友们好！

我最近在看lagrangian库中的solidParticle库代码，对于里面的 类模板的使用 不太懂，solidParticleCloud.C部分代码如下：

void Foam::solidParticleCloud::move(const dimensionedVector& g) //在solidParticleCloud类中定义一个move函数 { const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho"); const volVectorField& U = mesh_.lookupObject<const volVectorField>("U"); const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu"); interpolationCellPoint<scalar> rhoInterp(rho); interpolationCellPoint<vector> UInterp(U); interpolationCellPoint<scalar> nuInterp(nu); solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); //在solidParticleCloud类下面定义了一个trackingData类的对象td？？ Cloud<solidParticle>::move(td, mesh_.time().deltaTValue()); }

代码中的这句话：

solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value());

目前自己很困惑，为什么在solidParticleCloud类的函数move中可以 定义一个 solidParticle类中的 trackingData类的 对象td呢？

自己查找了solidParticleCloud.H文件发现 solidParticleCloud的基类是 Cloud<solidParticle>类，这个Cloud<solidParticle>类 中solidParticle是实参，对Cloud <ParticleType> 类模板进行了实例化.

我怀疑难道是因为Cloud<solidParticle>类中的 实参是solidParticle，所以其子类solidParticleCloud中函数 就可以调用/定义 solidParticle类中的trackingData类？

请哪位高手能够指点一下呢？感激不尽:xinxin: