Computational analysis of nonlinear bioheat transfer equation of multilayered brain tissue with porosity during tumor thermal therapy

Abstract The present model incorporates a bioheat transport framework within multilayered brain tissue, accounting for blood perfusion and thermal conductivity as temperature-dependent variables, while also considering porosity during magnetic fluid hyperthermia. Additionally, the skin’s surface is subjected to a nonlinear boundary condition of heat-flux kind. The treatment’s temperature level is crucial for tumor eradication, which is forecasted by resolving the current model utilizing the Crank–Nicolson scheme. To validate the current scheme, we compared with the Laplace transform solution in specific situation against both types of boundary conditions. The analysis indicates that augmenting the nonlinear thermal conductivity coefficient ( $$\alpha _{2}$$ α 2 ) and the blood perfusion coefficient ( $$\beta _2$$ β 2 ) reduces the temperature profile. Moreover, augmenting the porosity parameter ( $$\mathfrak {e}_i$$ e i ) substantially reduces tissue temperature. Furthermore, augmenting the magnetic parameters $$H_0$$ H 0 and f boosts the heating profile specifically within the tumor. An elevation in the convective heat transfer coefficient ( $$h_c$$ h c ) under Robin conditions reduces outer-layer temperatures, whereas an increased ambient temperature ( $$T_{00}$$ T 00 ) elevates the boundary temperature. The Neumann condition reduces outer-layer temperatures by increasing heat flux ( $$q_o$$ q o ). A comparative investigation demonstrates that the nonlinear model produces lower peripheral tissue temperatures and more effectively protects healthy tissue than the linear model, suggesting its superior appropriateness for therapeutic heating applications.