The objective of this study was to compare three different heat transfer models for radiofrequency ablation of in vivo liver tissue using a cooled electrode and three different voltage levels. The comparison was between the simplest but less realistic Pennes' equation and two porous media-based models, i.e. the Local Thermal Non-Equilibrium (LTNE) equations and Local Thermal Equilibrium (LTE) equation, both modified to take into account two-phase water vaporization (tissue and blood). Different blood volume fractions in liver were considered and the blood velocity was modeled to simulate a vascular network. Governing equations with the appropriate boundary conditions were solved with Comsol Multiphysics finite-element code. The results in terms of coagulation transverse diameters and temperature distributions at the end of the application showed significant differences, especially between Pennes and the modified LTNE and LTE models. The new modified porous media-based models covered the ranges found in the few in vivo experimental studies in the literature and they were closer to the published results with similar in vivo protocol. The outcomes highlight the importance of considering the three models in the future in order to improve thermal ablation protocols and devices and adapt the model to different organs and patient profiles.