CFD中文网

谢谢李老师，那我是不是直接用turbulence->nut()就可以获得混合湍流运动粘度了？我是看到createFields.H定义了这样的指针，对应找到他的 incompressible::turbulenceModel 类，就发现mut和nut是一样的，对于这种多相流的，它实际用的是哪个类，在哪里体现了密度的影响？对openFOAM不太了解，谢谢老师指点
image.png

怎样输出与一个patch发生碰撞的颗粒的速度，位置信息

看起来像个bug。你可以试试双核的，然后手动把traslationRotation写进processor里面，然后看看报错不。不报错的话可能是bug

@lys 你可以看下你的solver中有没有添加storeGlobalPositions这个函数，我对你使用的这个solver不是很了解，拉格朗日是你自己添加进去的吗

@dxl 我印象中拉格朗日这边OF5以后发生了一些变化，不知道你用的什么版本。
但是归根到底就是在拉格朗日那边有一个判断。后续版本好像转为p.fraction了
如果dt太小，这拉格朗日的计算就跳过了

https://github.com/OpenFOAM/OpenFOAM-8/blob/master/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
line 334-348
// Avoid problems with extremely small timesteps
if (dt > rootVSmall)
{
// Update cell based properties
p.setCellValues(cloud, ttd);

p.calcDispersion(cloud, ttd, dt); if (cloud.solution().cellValueSourceCorrection()) { p.cellValueSourceCorrection(cloud, ttd, dt); } p.calc(cloud, ttd, dt); }

好久没碰了，特别细节忘记了。

@xpqiu
但是有些资料里面是有这样的设置的

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