关于QMOM矩初值的给定

[walten]

@东岳 感谢感谢！这个问题困扰我好久了，就是不知道初始矩是怎么算出来的。。

[李东岳]

@walten 看附录 https://www.jianguoyun.com/p/DXZ5oIoQ9s3ZBhiHkaAC

[walten]

@东岳 东岳老师，您17年AIChE文章里的公式和您刚分享的文章里的公式不太一样，不知道哪个是笔误。

[李东岳]

93忘记了个根号，一会我更新下

[李东岳]

更新了，你那面算出来了么

[walten]

@东岳

东岳老师您好，您分享的文章中N值按理说就是文献中矩阵标准化之前的M0值，即单位体积的粒子数量，我按照您文章中的算法计算了N 与文献中M0差距了5个数量级。感觉M0还是没有算对。然后对计算值和文献值的矩进行标准化后，也还是存在差别。标黄的部分是B列的计算公式，您看下那里出问题了。。。感激不尽

[李东岳]

@walten 粒子直径单位都是米么？

[walten]

@东岳 颗粒平均直径 d10 我用的米的单位，即颗粒粒径是235微米。

[walten]

@东岳

文献中的一些物性参数。

[李东岳]

你的d10 = 0.000235，你的计算的$m_1/m_0=0.26$？这两个值应该相等

[walten]

@东岳 是的，0阶矩是颗粒总数量，1阶是颗粒总长度，m1/m0应该是平均粒径，所以我计算时是不是哪里出问题了？我也还在排查，而且N的计算值也不对。

[walten]

@东岳 东岳老师，我把计算excel表中A3和A4 两个数值的公式中的log改成了ln，得到的结果比较一致了，但是N仍然不对。您看下。

[walten]

@东岳 东岳老师，C语言中的语法是 log 是以e为底，log10 是以10为底。。但是平时书写log是以10为底，以e为底是不是习惯写成ln。所以看到您文章里的log我就想当然的用了10为底的log，excel中log也是默认以10为底。。所以会有以上的计算误差。现在剩下的问题就是N的计算了。。

[walten]

@walten 东岳老师，我试出来了。还是用的您给的公式，
只不过这里的αd不是指床层的固含率，而是指 所有填充颗粒占整个反应器的体积分数，这样的话文献中的初始固含率0.6经过换算应该为0.033。然后代入N公式，最后计算出来的值与文献基本一致。

感谢东岳老师百忙之中的指点，我终于解开了这个谜团了。。。
另：我又发现了您分享文献里的几个笔误。
1：公式94下边第一句话，standard 写成了dtandard。
2：公式96分母右边括号多了一半。

[李东岳]

@walten 试出来就好！谢谢，我回去更新下

[walten]

@东岳 就是还是那个问题，书写习惯中log是不是一般是指以10为底的对数。如果这样的话，您文章里应该写成ln，避免引起误会。

[李东岳]

@walten 感谢，已更新，也可以看看前面的内容，矩方法这面算法的问题非常多，尤其是realizablitity的问题，Fluent并没有处理realizablitity和boundedness的问题，不知道你后续会不会发现

[walten]

@东岳 我现在接触的还很浅。。以后碰到这方面问题再来请教~~