Efficient method to generate time evolution of the Wigner function for open quantum systems