DQMOM might be prone to numerical issues as the quadrature is dynamically updated for every cell and sometimes the problem gets so stiff. For that reaon I highly recommend that you continue experimenting first with inhomogeneous class method. If there you still have issues then you should revisit the way you are providing the bings and the closure laws.