CFD中文网

@李东岳 东岳老师好，分了32核，但是运行的时候也没有全部调用，不知道是不是线程多了内存需求也会变大？另外如果直接在linux系统跑的话，只有计算时间能缩短吗，内存占用这些是不是和虚拟机一样呢？

@ansel https://www.cnblogs.com/Xiwang-Sun/p/17463243.html
可以参考这个网址按这个画，一般来说网格大小尽量一致并且不要太尖锐，所以圆柱中间用一个长方体替代可能更加合适。

我的openfoam流体网格大概3000万左右，LIGGGHTS中预估颗粒数量也应该会几十万，之前试算过纯CFD的，已经很耗费计算资源了，我是在虚拟机里跑的，需要给虚拟机分配快200g内存和四线程并行计算，才能持续计算，因为流场很快稳定了，大概也就几天就算完了，但是后续CFDEM的话要持续加入颗粒，计算时间肯定不会短，想问各位老师有什么更好的解决办法？我的资源分配是不是合理呢？因为电脑就256g内存，分了200g去虚拟机其实已经很多了，而且我试了下线程越多好像对内存的占用也更多了。如果这些资源占用是合理的话，那可能只能去租服务器了

时隔一个月说一下结果，这样做不会报错，不过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并行跑了下，第一个时间步就报错了，我想可能是网格质量导致的问题。