An automatic grid generation algorithm for complex computational geometry is proposed to deal with the moving phase change interface during a cryosurgerical process. When the algorithm is incorporated into the dual reciprocity boundary element method, the freezing process of the biological bodies embedded with tumor and subjected to surface cryo-probe is successfully simulated. Some of the heat transfer mechanisms are thus better understood.