The interaction between thrust plume and satellite base region was investigated by using direct simulate Monte-Carlo calculations. For the accurate simulation of N2 and H2 collisions and rotation-translation transition, a variable soft-sphere model and a recent rotational relaxation model of N2 and H2 were used. For the investigation of the interaction between thrust plume and base region, the number density distribution for each species, translational and rotational temperature distributions, heat flux, and pressure were examined by direct simulation of Monte-Carlo calculations. It was found that most of the surface properties are affected by H2 collisions and a strong non-equilibrium state is observed at the base region. It was demonstrated that an accurate model is needed to simulate H2 collisions and the rotation-translation transition. The results by the present calculation are more accurate than previous direct simulation Monte-Carlo calculations because more accurate rotational relaxation models were used in simulating the inelastic collisions.