A two-dimensional simulation of earth-contact heat transfer by applying the finite element method is presented. Relevant thermal properties for each material type are determined and employed in the numerical model. Initial conditions and boundary conditions are carefully defined to improve the accuracy of this numerical investigation. The results demonstrate rational connection between the measured data set and simulated result. This research provides a useful contribution and reference to the earth-contact simulation.