Numerical schemes for solving the non-stationary Boltzmann-Poisson system for two-dimensional semiconductor devices