Numerical models simulating the process of NAPL from contaminated soil or groundwater through steam injection can be a useful tool for designing and evaluating the cleanup strategy under various field conditions. One and two dimensional numerical analyses were conducted based on the governing equations describing the NAPL removal as a non-isothermal, multi-phase and multi component process. Relatively good agreements were obtained between the numerical results and the observations from one-dimensional laboratory experiment, except some discrepancy due to experimental difficulties. Simulation effectively identified the steam displacement process of xylene floating on the water table and TCE sinking on the aquifer bottom in a two-dimensional analysis. Overall, simulation models have a high potential in the design/appraisal of a system for field application of the technique as well as in the examination of complex processes such as vaporization which is hard to identify experimentally.