Numerical simulations of two-photon absorption time-resolved photoluminescence to extract the bulk lifetime of semiconductors under varying surface recombination velocities