The validity of variable sphere molecular model for the Monte Carlo simulation of rarefied gas flow is examined by the simulation of thermal creep flow in a two-dimensional cavity. Macroscopic profiles in the cavity, such as density, velocity, and temperature profiles are quite in good agreement with those results of inverse power law potential.