This study presents a highly nonlinear dual-phase lag (DPL) bioheat model for simulating heat transfer in spherical biological tissues subjected to a focused laser heat source. The formulation integrates temperature-dependent metabolic heat generation and blood perfusion with an external laser heat source, yielding a strongly coupled nonlinear system of governing equations. The numerical solution is obtained through a hybrid framework that combines the implicit finite difference method with Newton-Raphson iterations, ensuring robust and efficient convergence in the strongly nonlinear regime. Results reveal that accounting for strong nonlinearity produces substantially higher peak temperatures and significantly modified heat penetration profiles compared to the linear counterpart. These findings are further supported by numerical simulations, experimental measurements, an alternative finite element approach, and comparisons made with existing literature, underscoring the necessity of nonlinear modeling for the accurate forecasting of thermal response and further enhancement of the accuracy and safety of laser biotherapy.
Khalid et al. (Mon,) studied this question.