A finite element model for radiofrequency ablation of the myocardium

A finite element model was developed to simulate the temperature distributions produced by radiofrequency catheter ablation. This model incorporated blood, myocardium and torso tissues. The Laplace equation was solved to determine the steady-state electric field. The heat generation in the tissues was then computed from the power density distribution and the bioheat equation was solved to determine the time-varying temperature distribution, taking into account the convective energy exchange at the blood-myocardium and torso-air interfaces. This model was used to predict the lesion depth and to evaluate the effects of electrode location, changes of the electrical and thermal conductivities, and the electrode radius on the thermally induced damage to the myocardium. Temperature distributions induced by radiofrequency ablation were found to be: i) not very sensitive to the reference electrode location, ii) more sensitive to electrical conductivity changes than to thermal conductivity changes, and iii) larger electrodes allow a current distribution at higher level of power with reducing the chance of impedance rise.