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.
