A globally hyperbolic moment system upto arbitrary order for the Wigner equation was derived in [6]. For numerically solving the high order hyperbolic moment system therein, we in this paper develop a preliminary numerical method for this system following the NRxx method recently proposed in [8], to validate the moment system of the Wigner equation. The method developed can keep both mass and momentum conserved, and the variation of the total energy under control though it is not strictly conservative. We systematically study the numerical convergence of the solution to the moment system both in the size of spatial mesh and in the order of the moment expansion, and the convergence of the numerical solution of the moment system to the numerical solution of the Wigner equation using the discrete velocity method. The numerical results indicate that the high order moment system in [6] is a valid model for the Wigner equation, and the proposed numerical method for the moment system is quite promising to carry out the simulation of the Wigner equation.

A new adaptive cell average spectral element method (SEM) is proposed to solve the time-dependent Wigner equation for transport in quantum devices. The proposed cell average SEM allows adaptive non-uniform meshes in phase spaces to reduce the high-dimensional computational cost of Wigner functions while preserving exactly the mass conservation for the numerical solutions. The key feature of the proposed method is an analytical relation between the cell averages of the Wigner function in the k-space (local electron density for finite range velocity) and the point values of the distribution, resulting in fast transforms between the local electron density and local fluxes of the discretized Wigner equation via the fast sine and cosine transforms. Numerical results with the proposed method are provided to demonstrate its high accuracy, conservation, convergence and a reduction of the cost using adaptive meshes.