This work presents a bivariate extension of the moment projection method (BVMPM) for solving the two-dimensional population balance equations involving particle inception, growth, shrinkage, coagulation and fragmentation. A two–dimensional Blumstein and Wheeler algorithm is proposed to generate a set of weighted particles that approximate the number density function. With this algorithm, the number of the smallest particles can be directly tracked, closing the shrinkage and fragmentation moment source terms. The performance of BVMPM has been tested against the hybrid method of moments (HMOM) and the stochastic method. Results suggest that BVMPM can achieve higher accuracy than HMOM in treating shrinkage and fragmentation processes where the number of the smallest particles plays an important role.