Magnetic　hyperthermia　regulates　the　therapeutic　temperature　within　a　specific　range　to　damage　malignant　cells　after　exposing　the　magnetic　nanoparticles　inside　tumor　tissue　to　an　alternating　magnetic　field.　The　therapeutic　temperature　of　living　tissues　can　be　generally　predicted　using　Pennes＇　bio-heat　equation　after　ignoring　both　the　inhomogeneity　of　biological　structure　and　the　microstructural　responses.　Although　various　of　the　bio-heat　transfer　models　proposed　in　literature　fix　these　shortages,　there　is　still　a　lack　of　a　comprehensive　report　on　investigating　the　discrepancy　for　different　models　when　applied　in　the　magnetic　hyperthermia　context.　This　study　compares　four　different　bio-heat　equations　in　terms　of　the　therapeutic　temperature　distribution　and　the　heat-induced　damage　situation　for　a　proposed　geometric　model,　which　is　established　based　on　computed　tomography　images　of　a　tumor　bearing　mouse.　The　therapeutic　temperature　is　also　used　as　an　index　to　evaluate　the　effect　of　two　key　relaxation　times　for　the　phase　lag　behavior　on　bio-heat　transfer.　Moreover,　this　work　evaluates　the　effects　of　two　blood　perfusion　rates　on　both　the　treatment　temperature　and　the　cumulative　equivalent　heating　minutes　at　43　degrees　C.　Numerical　analysis　results　reveal　that　relaxation　times　for　phase-lag　behavior　as　well　as　the　porosity　for　living　tissues　directly　affect　the　therapeutic　temperature　variation　and　ultimately　the　thermal　damage　for　the　malignant　tissue　during　magnetic　hyperthermia.　The　dual-phase-lag　equation　can　be　converted　into　Pennes＇　equation　and　simple-phase-lag　equation　when　relaxation　times　meet　specific　conditions　during　the　process　of　heat　transfer.　In　addition,　different　blood　perfusion　rates　can　result　in　an　amplitude　discrepancy　for　treatment　temperature,　but　this　parameter　does　not　change　the　characteristics　of　thermal　propagation　during　therapy.