Implementation of bivariate population balance equations in CFD codes for modelling nanoparticle formation in turbulent flames