Solving GPBE

Recently, I am trying to understand the state of the art for solving GPBE. However, I have no experience in this continum approach. So, I wrote a miniparagraph and posted it here asking for OpenFOAMers to help discuss remaining issues and miss understandings of the literature. Thanks in advance.
 The generalized population balance equation (GPBE) (Marchisio and Fox, 2013) describes the evolvement of particle number density function in dimensions of time, space, velocity and particle properties, such as size. However, it is infeasible to directly solve this integraldifferential equation due to its high dimensionality and the unclosed term such as particle collisions.
 To reduce the transport equation into only time and space dimensions, both the velocities and particle properties should be integrated out. This results in at least four dimensions of moments. High order multivariable momentinversion algorithm used to find the weights and abscissas is very challenging (Fox, 2008). One approach to solve this problem is assuming that the velocity depends on particle size (Buffo et al., 2013). Then, for a specific particle size, the velocity distribution is integrated out and result in a macro scale averaged particle velocity conditioned on particle size. This macroscale velocity is then solved using hydrodynamic models such as twofluid model or multifluid model (Fan et al., 2004; Buffo and Marchisio Daniele, 2014). Note that, the momentum equation of the hydrodynamic model used in the twofluid model or multifluid model is derived from Boltzmann equation which doesn’t include the effect of particle size evolution. Solving GPBE without these assumptions is rarely seen due to its high dimensionality.

It looks correct at least from my side

AFAIK, you are the only one working in this area. Any update?