The goal of this research is developing a unified numerical method for simulating continuum and transitional flow. To achieve our ultimate goal, first, hyperbolic-relaxation equations are introduced, then a new discretization method is developed. The method is based on Huynh’s upwind moment scheme, with implicit treatment of the source term. Our previous linear method is generalized to 1-D nonl...