Implementation of the population balance equation in CFD codes for modelling soot formation in turbulent flames