This paper considers the three-phase lag (TPL) bioheat model to study the phase change phenomena in skin tissue during cryosurgery. The considered TPL model is based on the model of thermo-elasticity, i.e., the combination of the rate of thermal conductivity and new phase lag due to thermal displacement. An effective heat capacity-based numerical algorithm is established to solve the nonlinear governing equation for biological tissue freezing. Space and time derivatives appearing in the mathematical model are approximated using the radial basis function (RBF) and finite difference method (FDM), respectively. The impact of three nonclassical models, single-phase lag (SPL), dual-phase lag (DPL), and TPL, on the freezing process is studied. The effects of phase lags involved in the models on freezing are also part of this study.