Benchmark solution for the prediction of temperature distributions during radiofrequency ablation of cardiac tissue

Several studies on radiofrequency (RF) ablation are aimed at accurately predicting tissue temperature distributions by numerical solution of the bioheat equation. This paper describes the development of a solution that can serve as a benchmark for subsequent numerical solutions. The solution was obtained using integral transforms and evaluated using a C program. Temperature profiles were generated at various times and for different convection coefficients. In addition, a numerical model was developed using the same assumptions made in obtaining the benchmark solution. Comparison of surface and axial temperature profiles shows that the two solutions match very closely, cross validating the numerical methods used in evaluating both solutions.