A normal mode analysis for characterizing the temperature fluctuation in tissues was proposed based on the Penne’s bio-heat transfer equation. Closed-form analytical solutions to obtain the heating pattern due to the propagation of ultrasonic waves in tissue system are presented. The evaluation of temporal and spatial distributions of temperature is investigated with the effect of relaxation time. The derived method is evaluated with numerical simulations in 2D which are applied to tissue medium in simplified geometry.