For numerical simulation of the rarefied gas flows, the classical Navier-Stokes-Fourier (NSF) equations and Fourier’s law of heat conduction based on the continuum flow assumption are no longer valid, while solving the Boltzmann equation based on gas kinetic theory requires enormous computational resources. To balance computational accuracy and cost in the numerical simulation, a hybrid algorithm of the discrete velocity method (DVM) and moment method (MM) is developed for rarefied gas flows in the early transition regime to extend the validity range of the moment method. Implicit DVM in combination with Maxwell’s wall boundary condition is employed in the near-wall region, while the moment equations are used to describe the bulk flow field.