MC simulation of double-resonant-phonon depopulation THz QCLs for high operating temperatures