The Dual Reciprocity Boundary Element Method (DRBEM) is extended to simulatethe thermal wave propagation in biological tissues. The higher the thermal relaxation timeis, the stronger the thermal wave effect will be. Under changing heat source, bioheat trans-fer has distinct wave characters. The thermal wave propagation in biological tissues obeysthe superposition and resolution principle of ordinary wave. Reflected by a rigid wall’ (thefirst boundary condition), the thermal wave will show a phase jumping phenomenon. TheDRBEM is an efficiellt pure boundary iotegral method without domain integral for solvingthermal wave problems. Thermal wave and their refiection, phase jumping, superposition,resolution can be correctly located and sharply captured. There are no the oscillatory behav-ior in the wave front and wave peak region, which is presented in reported finite differencesolution with TVD high accuracy scheme.
