多核运行时出现这个错误，减少核数运行正常输出命令添加执行命令解决。

下面我们具体来进行操作，首先需要确保自己文件夹内的图片的名称是按照数字大小顺序排的，一般批量导出的图片都会自动命名，所以这里一般不需要其它操作，例如，我的文件夹下的图片是按照11.xxxx命名的，第一张图是11.0000，第二张是11.0001，以此类推。11.%04d.jpg 显然这里是图片名称的索引，注意，这里%04d代表四个位数的整数排序，前面的11.即代表前缀。这里的只采用了最基本的命令。

1、对于特征值运算eig()函数来说，GPU的加速效果是很明显的2、如果要求精度不高，采用单精度计算，加速效果更加明显。

长话短说，关于这个问题，有两个解决方案，一是将文件类型转换为VTK格式，二是采用github上的脚本，将OpenFOAM输出的位置文件转换成真实的颗粒的真实位置文件。

Re = 1000时，顶盖驱动的方腔流动速度云图和流线图如下：可以看该流动存在三个涡心，那么如何获得涡心的具体坐标呢，一开始我以为会有什么定量的计算方法，于是在网上找了一下，发现知乎上于是我去tecplot用户手册具体看了一下,发现这个命令在如下界面可以发现有两种计算方法，一种是基于速度梯度的，另一种是基于涡量方向的，点击Extract之后发现会报错，似乎只能在三维情况下使用。具体我还没有在三维流

但是其计算的结果是正确的，我开始就想当然的按着这个思路去做，但是那个代码是按照MRT来写的，我将其转换成LBGK时，没有执行第二点，于是计算结果不对，我以为是碰撞过程那个外力项Fa的问题，于是找来找去，一直debug不出来，后来我执行了第二步，但是我开始理解错了，我以为是程序的最后，在收敛之后在对速度上加上一项。这里面其实讲的十分清楚，但是我开始编程的时候却遇到了一些问题，本站里的那篇博文提供的代

在一个交流群里有人推荐simpleTex，免费好用，试了一下，果然还可以，最重要的是免费的。mathpix是识别公式将其转换为latex代码的非常好用的工具。

其实这里还有一点不是很理解，就是为什么温度梯度那一项经过平均后除以分母会等于1，如果是pure conduction的话，全场内任意处的温度梯度应该是一致的，因为温度会线性分布，但是现在是对流，全场温度梯度经过平均后会等于pure conduction的情况吗？在传热传质基础一书中，努塞尔数的定义如下，对流（convection）传热与纯热传导(pure conduction)传热之比。这里的下标

写的还挺详细的，不错。

虚拟机：链接：https://pan.baidu.com/s/1MPMTJQfl76mW0H5bbT_rAg提取码：rqli开机密码：530944988知乎博客：作者说明了如何关闭颗粒碰撞计算来达到提升计算速度。Github：这个是github上的CFDEM，在原来的基础上增加了许多求解器。官网：官网教程，但是比较鸡肋。