Nanoimprint lithography(NIL) is becoming next generation lithography of significant interest due to its low cost and a potential patterning resolution of 10nm or less. Success of the NIL relies on the adequate conditions of pressure, temperature and time. To have the adequate conditions for NIL, one has to understand the polymer flowing behavior during the imprinting process. In this paper, an analytical approach of polymer flow in thermal NIL was performed based on the squeeze flow with partial slip boundary conditions. Velocity profiles and pressure distributions of the polymer flow were obtained and imprinting forces and residual thickness were predicted with the consideration of the slip velocity between the polymer and the mold/substrate. The results show that the consideration of the slip is very important for investigating the polymer flow in Thermal NIL.