第一步，在Define>Custom Field Function Calculator 里编辑一个所求量的算法。
第二步，在File>Export>During Calculation>Solution data里选择 File Type。如果有专业软件，选择相应格式，如果用Excel处理，选择ASCII码。 Location 选择Cell center或node 看你自己需求。
第三步， Frequency里选择保存频率。是时间步长的倍数。
第四步，在surface里选择需要的表面，你的话应该选择outlet
第五步， quantities里选择刚才定义的算法方程。
第六步， OK
横坐标小r是指画图点的位置，R是反应器半径，比值是指该点位于反应器的径向位置，比如，0 就是值反应器中心，-1，就是左壁面，1 就是右边壁面，能迅速定位到径向位置。
@walten 东岳老师，我试出来了。还是用的您给的公式，
只不过这里的αd不是指床层的固含率，而是指 所有填充颗粒占整个反应器的体积分数，这样的话文献中的初始固含率0.6经过换算应该为0.033。然后代入N公式，最后计算出来的值与文献基本一致。
感谢东岳老师百忙之中的指点，我终于解开了这个谜团了。。。
另：我又发现了您分享文献里的几个笔误。
1：公式94下边第一句话，standard 写成了dtandard。
2：公式96分母右边括号多了一半。
@东岳 东岳老师，C语言中的语法是 log 是以e为底，log10 是以10为底。。但是平时书写log是以10为底，以e为底是不是习惯写成ln。所以看到您文章里的log我就想当然的用了10为底的log，excel中log也是默认以10为底。。所以会有以上的计算误差。现在剩下的问题就是N的计算了。。
@东岳 东岳老师，我把计算excel表中A3和A4 两个数值的公式中的log改成了ln，得到的结果比较一致了，但是N仍然不对。您看下。
@东岳 是的，0阶矩是颗粒总数量，1阶是颗粒总长度，m1/m0应该是平均粒径，所以我计算时是不是哪里出问题了？我也还在排查，而且N的计算值也不对。
