We introduce a numerical method for solving Grad's moment equations or regularized moment equations for arbitrary order of moments. In our algorithm, we do not explicitly need the moment equations. Instead, we directly start from the Boltzmann equation and perform Grad's moment method. We define a conservative projection operator and propose a fast implementation which makes it convenient to add up two distributions and provides more efficient flux calculations compared with the classic method using explicit expressions of flux functions. For the collision term, the BGK model is adopted so that the production step can be done trivially based on the Hermite expansion. To regularize the system, we apply the Maxwell iteration to the infinite moment system and determine the magnitude of each moment with respect to the Knudsen number. After that, we obtain the approximation of high order moments and close the moment systems by dropping some high-order terms. Linearization is then performed to obtain a very simple regularization term, thus it is very convenient for numerical implementation.
view more