mfix中输出DEM颗粒的固相速度到网格

基于mfix-19.1.2版本

方法一：直接输出差值网格固相速度

注：这种方式只适用于garg 2012颗粒差值格式
在DEM中，默认是无法输出固相速度的网格值的：

但是通过搜索des文件夹下V_s关键字可以看到，实际上网格固相速度是会通过差值计算得到的：

并且在计算DEM颗粒脉动速度的地方也用到了网格固相速度：

说明网格固相速度只是默认没有写出，因此只需要人为写出一下即可。

首先找到网格输出的代码位置：

由于MMAX代表TFM的固相数，因此在循环范围部分需要修改，默认只写出TFM的网格固相速度。修改如下：

最后在dat文件里手动添加一下VTK_U_S(VTK, PHASE)这个关键字（且确保颗粒的差值格式为garg 2012）例如：

结果如下：

方法二：保存网格内颗粒的平均速度到网格

网格颗粒速度还可以仿照DES_GRANULAR_TEMPERATURE函数内计算网格拟颗粒温度的方式，对网格内颗粒速度进行统计平均。
按照J. Musser的建议，可以把需要后处理的变量存到ReactionRates这个数组中，虽然这个数组本来是用来后处理反应速率的，不过由于可以直接在gui里设置该数组的个数以及输出频率，而且该变量仅仅用来保存输出，计算中未被使用，因此可以灵活使用，来保存其他需要保存的变量。

例如这里我要保存网格内颗粒Y方向的平均速度V_s，可以先设置一个ReactionRates数组，以备后面使用：