In this paper, we apply a hybrid Nonequilibrium Green’s Functions and relaxation rate approximation approach to investigate frequency multiplication in the GHz to THz range delivering results in good agreement with experimental data for different input frequencies. A discussion of the current limitation of housing designs for the superlattices is given, based a simple model for the local electric field and in comparison with measured input powers delivered by a backward wave oscillator.