  In this paper, we propose a method to simulate the microflows with Shakhov
model using the NRxx method developed in [4, 5, 6]. The equation under
consideration is the Boltzmann equation with force terms and the Shakhov model
is adopted to achieve the correct Prandtl number. As the focus of this paper,
we derive a uniform framework for different order moment systems on the wall
boundary conditions, which is a major difficulty in the moment methods.
Numerical examples for both steady and unsteady problems are presented to show
the convergence in the number of moments.
