@疏影横斜水清浅 网格太大导致的精度误差，不应该跟软件有关系，只跟你想要研究的问题、你关注的现象的尺度以及数值格式这些有关系。说白了，只要案例设置的好，精度误差不会跟别的软件有明显区别。

你提到的计算大尺度的开源软件和程序是指哪些？我这边做风场模拟，对这些还比较感兴趣。

@李东岳 辛苦李老师了，方便的话可以试试看。后面那种场景，主要是公司内集群是所有专业混用，很容易出现这种情况。甚至单单只针对OpenFOAM来说，暂时没限制节点的使用，所以会导致一个节点上能同时跑好几个计算或者计算的一部分，分别是隶属于不同作业的。

@李东岳 老师，我自己觉得是可以，只是我没想好具体的实际方法，这也是为什么我选择在均匀的网格下弄。 我对于我现在均匀网格的就这么搞得 我自定义了我的固体活动范围里的网格编号和关系 均匀网格的连接关系很简单 我甚至连数组都不用开个去记录 我只是为了并行 预处理了固体活动范围内每个网格所在核心的编号和所处分块区域内我自定义编号的cell在对应核心对应分解域内的真实编号 计算时候根据连接关系，o(1)的运算时间取出临近网格的自定义标号 然后根据预处理记录的它在所属分块内的真实编号来计算 。

只是对于非结构网格，不均匀，有变形这种，我具体不太知道openfoam里编号的规律，也不知道老师您这里具体的应用场景是什么。我个人觉得，也是在进入时间步循环前，搞个预处理，把非结构网格的cell中心坐标提前循环一遍预处理，映射给自己自定义的一套均匀网格，然后一个数组记录，数组下标是自定义映射均匀网格的编号，该数组下标里数组对应内容是落在这个自定义网格内的非结构网格的真实编号，均匀网格的关系很简单，可以简单查找临近关系。而且这个映射网格的作用纯粹就是，映射个大致的位置关系。但是这个映射的均匀网格的尺寸要小于实际非结构最小cell的尺寸，大致我是这么想的。但是这样我不知道提前记录的内存开销有多大，我觉得甚至还不如，自己知道每个非结构网格要查找的周围感兴趣网格范围，每个模拟开始进入时间步前，提前预处理输出给文件打个表 （因为是预处理，所以无所谓，各种findcell也花不太久，因为只需要执行一次并记录），然后运算时候，对于某个非结构网格，直接读取预处理打表的那个文件，再取出周围符合关系的网格编号。这里读取文件，和字符串处理花一点运算时间。

可以。我记得湍流方程右边都是只考虑网格点的作用，不需要考虑边界。因此你的G声明了内部场，没有问题。

@李东岳 谢谢！: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），这种用粒子包去代替粒子的方法我没了解过，但肯定已经不是物理粒子了。所以要看哪个层面，哪种模型的碰撞聚合了。