离散法求解群体平衡方程中的程序含义

[Caijinjin]

在离散法求解群体平衡方程中，程序文件中颗粒聚并生成项（birthByCoalescence）和文献中的给的式子形式不太一样，请问哪位大佬可以帮忙解释以下程序中为什么还要多乘以颗粒的直径？

方程里面η在程序里面应该对应gamma，指的是颗粒的分配体积，c对应程序里的coalescenceRate

[Caijinjin]

我自己暂时的理解是，程序中fjfj.phase()/fj.x()和fkfk.phase()/fk.x()指代的是j和k的颗粒数目，也就是右边公式中的Nj和Nk。但是其中的fj.phase和fk.phase是什么含义就有点不太明白了。还有就是聚并速率的量纲是m3/s的，程序里面为什么还要再乘以i颗粒体积fi.x()？这样子的话量纲就不是每秒了。