An approach based on a discontinuous Galerkin discretization is proposed for the Bhatnagar-Gross-Krook model kinetic equation. This approach allows for a high order polynomial approximation of molecular velocity distribution function both in spatial and velocity variables. It is applied to model normal shock wave and heat transfer problems. Convergence of solutions with respect to the number of...