Numerical approximation of the Vlasov-Poisson-Fokker-Planck system in two dimensions