@李东岳 谢谢！:140: :140:

@李东岳 李老师你好，那麻烦你推荐几篇非牛顿流体+颗粒流方面的论文啊，多谢了

@李东岳 您能简单推测大概是什么地方出现了问题吗？

具体研究体系是气液固催化反应，催化剂粉末与液体作为拟均相混合，气液不互溶，固体只存在于液相，流动过程中夹杂催化剂的液体会被切碎成小液滴。本人小白看了论坛上很多帖子，目前感觉这个MPPICInterFoam最符合我的体系，不知道有没有用过的前辈可以指点一二

@李东岳 好的，李老师，我去调试看看

最近在做风浪相互作用的模拟时遇到一个问题：在气液剪切过程中，湍流的研究是非常重要的，我在加入湍流模型之后，计算结果耗散特别严重，但是没有气体作用（nowind）时，SSTk-w湍流模型能够较准确的模拟，从下面的水位图中也可以看出来。并且湍动能k和湍流粘度nut在出口outlet逐渐向计算域内部扩散导致波高衰减（见下图），但是不知道具体原因。单纯的平底地形结果也没问题，但是遇到到这种斜坡地形就有衰减。
目前我已经尝试过of自带的几乎所有RANS类的湍流模型，计算结果均耗散严重（以SSTk-w为例），也尝试过对边界条件和壁面边界层的各种处理方式，还是效果不理想，想请教一下各位大佬，想解决这个问题，我还应该从哪些地方入手？
二维模型如图：
模型图.png
水位图：
水位图.png
不同时刻湍动能k分布：
k-12s.png
k-30s.png
k-60s.png
不同时刻湍流粘度nut分布：
nut-12s.png
nut-30s.png
nut-60s.png

@李东岳 是的李老师，我理解成了这个边界用在入口时相当于固定值，由inletValue指定，用在出口相当于零梯度，应该是我理解错了。。。

@veen 您好，关于“Primary job terminated normally...”的报错您后面还有研究过吗？我最近也遇到了相同的报错信息，不知道该怎么处理

你可能需要仿照MRFZone类写一个平动的类了，但是相对来说要更加容易些，主要工作是修正面通量，推荐用三角分解计算。还要注意下include的边界面上速度修正。

@我是河滩 请问您这个问题解决了吗？我最近也遇到了同样的问题

或者你应该在intU[1]=...之前加个if(cellID > -1).

想问一下，在绘制二维截面的流线图时，在paraview中以列表形式导出了点的数据，存在以下疑问：
①列表中各列数据代表什么含义，比如normals,integrationtime,point ID?
②各点为什么会有block number之分?block number是什么含义？
③能否通过列表数据做出流线图？
④保存数据时，为什么每个时刻的列表文件会被保存成2个csv文件？
以上问题求指教
3e5026fe-ea05-4101-86ee-b5ad8e43bfad-1666618809(1).png

刚接触OpenFOAM，安装的是OpenFOAM-v2206，了解到需要用pisofoam或者pimpleFoam求解器，貌似在tutorials里面没有现成的NACA翼型低速绕流大涡模拟不可压并行算例。

@cyberk 抱歉，这个问题我没有遇到过，但我看了你之前用ICEM画的网格，这个六叶Ruston桨的搅拌槽的网格不太复杂的呀，这个你可以试试只画1/6个区域，然后对称过去，画完整体之后再去openfoam中用topoSet划分搅拌的区域

@李东岳
谢谢老师，虽然只是一个小建议和思路，却给我节约了大量的时间，目前已经在OF7中调通了mutiphaseInterFoam求解器。:xiexie:

@hotairballon1997

在您回帖的链接中好像没有看到关于这个问题的解答，

对呀，要分类讨论，就像你说的，是用 PBM，还是拉格朗日。
我发的那个链接里 MPPIC 和 DPMFoam 中用到的粒子都是粒子包（particle parcel），这种用粒子包去代替粒子的方法我没了解过，但肯定已经不是物理粒子了。所以要看哪个层面，哪种模型的碰撞聚合了。

@李东岳 感谢东岳前辈

@李东岳 已经解决，多谢李老师

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

@李东岳 非常感谢李老师，终于搞定，能输出数据了，但还需要进一步验证是否正确。差点放弃想升级新版本of。

volScalarField& summass = const_cast<volScalarField&>(p.mesh().objectRegistry::lookupObject<volScalarField>("summass")); const label cellIp = p.cell(); scalar& PPC = summass.internalField() [cellIp]; PPC += p.nParticle()*p.mass();

希望对以后需要的朋友有点借鉴作用。特别感谢 @星星星星晴 不厌其烦的解答:146: