The process of squeeze casting for metal matrix composites (MMCs) has been simulated numerically by using finite difference method. The governing equations to describe fluid flow through porous medium and heat transfer are applied to two dimensional model which is similar to a real system. A computational code has been developed to solve this problem. The influence on infiltration kinetics and solidification time of several parameters is investigated. Cooling curves and temperature distribution with time and position is also shown. The result can be used to design the squeeze casting for MMCs.