CFD中文网

时隔一个月说一下结果，这样做不会报错，不过liggghts运算显示的是1×1×4（假设分成四个processor），但是of里的确是scotch划分的，不知道liggghts有没有自动识别

LIGGGHTS中stl网格导入报错，复杂的stl网格导入liggghts中总是报这样的错，看起来是网格质量问题，有没有老师遇到过，用什么软件去修复可以成功呢？

@李东岳 是的，我也不知道liggghts能不能自动识别，想看看论坛里有没有老师恰好做过这部分内容，谢谢东岳老师

我的网格模型非常复杂，因此在of中用scotch方法划分计算域，但是liggghts并行计算中in脚本需要写processors * * * ，其中 * 代表三个维度的划分个数，这种情况应该如何确定*值呢？

@chapaofan 好的，谢谢

@chapaofan 您好，想问个题外话，CFDEM并行运算相较于单核运算会快很多吗，我有试过并行，但是相较于单线程运算好像计算效率并没有明显提升

图像和代码里都显示没有加边界层，另外可以尝试边缘处的加密等级提高一点试试。

@李东岳 我又试了几次，总是跑一会就库朗数太大发散了，我觉得可能是我shm中minvol的体积改小了的缘故，又或者是时间步长太大了，反复调的心累。不过还是非常感谢东岳老师的耐心回答！

@李东岳 icofoam，我试着算过，都报错了

@李东岳 东岳老师，非常感谢您的耐心回复。我刚才又重新跑了一下，这次600多w网格，网格看起来挺不错的，但依旧有skewness的错误，我把全过程的log都记了下来，但是log里显示shm过程中好像没有skewness大于4的值 ，为什么checkmesh中会有呢？
log.CFDrun

@李东岳 第一次的是100多万网格，第二次就是1600w网格了，但是一旦加上feature的细化，网格质量就有问题，不过我再试一下看看，谢谢东岳老师。

如果不考虑feature的细化level的话，全都是六面体网格网格质量是很好的，我在想可能真的是stl文件的原因吧，毕竟过于复杂，如果没什么好的办法的话，那只能不考虑其feature的细化了，只用第一次第二次那样的网格。

@李东岳 东岳老师好，这是某种多孔隙岩石，我想研究流体的其内部的渗透，stl文件是从ct处理软件中导出来的，我后面又试了些shm的参数，目前网格质量问题主要是skewness会有偏大的值，导致后续计算不收敛。

@李东岳 东岳老师好，我进行了三次对比模拟，主要改变的是shm文件中的两个参数。结果如下
第一次：

第二次：

第三次 ：

其实我想要的是第三次的这种网格结果，网格边缘更加平滑，但可惜网格质量不ok。请问一下老师，这种应该往哪个方向去调整？显示贴合与隐式贴合我也都试过，没有太大区别。我上传了blockmesh和snappyhexmesh文件，以及三次的log文件。
blockMeshDict
snappyHexMeshDict
log.shm1
log.shm2
log.shm3

@李东岳 谢谢东岳老师，我再重新画网格看看，是blockmesh里画粗一点，snappyhexmesh细化等级高一些好还是blockmesh里画细一点，snappyhexmesh细化等级低一点好呢？不知道有没有普遍的结论，我都试试看看。

第二个问题解决了，现在还是网格质量的问题，有没有用snappyhexmesh画过比较复杂的网格的老师可以指点一下。
我用icoFoam并行跑了下，第一个时间步就报错了，我想可能是网格质量导致的问题。

我在进行snappyhexmesh并行划分网格遇到两个问题：
1.首先是网格看起来大部分位置很好，但是checkmesh后有错误（我的stl文件比较复杂），我试过不对surface进行精细化，网格checkmesh后没有问题，但是就会缺少大量的网格细节，效果不好。问问各位老师有什么比较好的办法可以得到质量好的网格。


2.进行CFD计算时报错：processor进程文件夹中的0文件下的各种物理量都没有stl文件的信息（边界条件等等）
前者是0文件夹下p文件，后者是processor进程文件夹的p文件，缺失reef3cm（stl网格面）的信息