We demonstrate here, for the first time, the constitutive scaling approach applied to simulate a fully compressible, non-isothermal micro gas flow within a mainstream computational physics framework. First, the physics underlying these new constitutive-relation scaling models for rarefied gas flows at the microscale, in particular, the Knudsen layer, is discussed. Results for Couette-type flows in microchannels, including heat transfer effects due to rarefaction, are then reported and we show comparisons with both traditional Navier-Stokes-Fourier solutions and independent numerical studies. We discuss the limitations of the constitutive scaling process, such as the breakdown of the model as the Knudsen number increases and the influence of the wall interaction model on the numerical results. Advantages of the constitutive scaling technique are described, with particular reference to the practicality of using it for microscale engineering design.