In the present study, a numerical alogrithm is developed in order to simulate the non-isothermal filling and compression stages of the injection-compression molding process. The generialized Hele-Shaw flow model combined with the control volume/finite element method is implemented to predict the melt front advancement and the distributions of pressure, temperature and flow velocity dynamically during the molding process. From on the simulated results, it was found that the compression speed is the most important factor that affect the molding pressure. Under proper compression speed, the predicted molding pressure is lower than that required by the purely conventional injection molding and the part residual stress can be minimized. It was also found that part temperature is more uniform at the end of filling process by using injection-compression molding. If the switch time from injection to compression is too low, it may also cause high molding pressure.