Thermal ablation modeling via the bioheat equation and its numerical treatment