Modelling the thermalisation of non-equilibrium phonon distributions in bulk single crystals