Consistent estimation of the bispectral density function of a harmonizable process