Thermal-electrical finite element modelling for radio frequency cardiac ablation: effects of changes in myocardial properties

Finite element (FE) analysis has been utilised as a numerical tool to determine the temperature distribution in studies of radio frequency (RF) cardiac ablation. However, none of the previous FE analyses clarified such computational aspects as software requirements, computation time or convergence test. In addition, myocardial properties included in the previous models vary greatly. A process of FE modelling of a system that included blood, myocardium, and an ablation catheter with a thermistor embedded at the tip is described. The bio-heat equation is solved to determine the temperature distribution in myocardium using a commercial software application (ABAQUS). A Cauchy convergence test (epsilon = 0.1 degreesC) was performed and it is concluded that the optimal number of elements for the proposed system is 24610. The effects of changes in myocardial properties (+/- 50% electric conductivity, + 100%/-50% thermal conductivity, and + 100%/-50% specific heat capacity) in both power-controlled (PCRFA) and temperature-controlled RF ablation (TCRFA) were studied. Changes in myocardial properties affect the results of the FE analyses of PCRFA more than those of TCRFA, and the maximum changes in lesion volumes were -58.6% (-50% electric conductivity), -60.7% + 100% thermal conductivity), and +43.2% (-50% specific heat).