We present an extended ensemble Monte Carlo approach, allowing for the self-consistent modeling of terahertz difference frequency generation in quantum cascade lasers. Our simulations are validated against available experimental data for a current room temperature design. Tera-hertz output powers in the mW range are predicted for ideal light extraction.