NEW COMPUTATIONAL APPROACHES FOR BIOPHYSICAL HEAT TRANSFER IN TISSUE UNDER ULTRASONIC WAVES: THE VARIATIONAL ITERATION AND CHEBYSCHEV SPECTRAL SIMULATIONS

Abstract

A mathematical and numerical study is presented for simulating temperature distribution in a two-dimensional tissue medium using Pennes bioheat transfer equation, when the tissue is subjected to ultrasonic waves. Following nondimensionalization of the governing partial differential equation, a novel variational iteration method (VIM) solution is developed. This excellent technique introduced by He [Variational iteration method — a kind of non-linear analytical technique: Some examples, Int J Non-Linear Mech.34:699–708, 1999] employs Lagrange multipliers which can be identified optimally via variational theory. The space and time distributions of temperature are studied and solutions visualized via Mathematica. The influence of thermal conductivity and relaxation time are also examined. Excellent stability and convergence characteristics of VIM are demonstrated. Validation is achieved with a Chebyschev spectral collocation method (CSCM). The present work demonstrates the excellent potential of this powerful semi-numerical method in nonlinear biological heat transfer and furthermore provides an alternative strategy to conventional finite element and finite difference computational simulations. The model finds applications in minimally-invasive spinal laser treatments, glaucoma therapy in ophthalmology and thermoradiotherapy for malignant tumors.