Abstract: Numerical simulations with the Galerkin’s finite element method are made to study steady-state flow under different simulation conditions at a vertical well site of soil vapor extraction (SVE). Flow-field simulation shows that radial pressure gradient goes down fast and air flowrate is nearly directly proportional to the well vacuum and the soil permeability coefficient. It also indicates that stop-off coating will play an important role in SVE. This axisymmetric 2D flow model can predict vacuum distribution, which gives a much better fit to field data by contrast to 1D model. The numerical modeling process for SVE remediation is achieved after a steady alternant direction implicit (ADI) finite difference format is developed. In this paper, as an example, the 60 days remediation under 4 different conditions is simulated, which proves that the clear-up result is mainly governed by soil permeability. This work will be used as a tool for engineering design and prediction for SVE remediation.