The comparison of lesion outline and temperature field determined by different ways in atrial radiofrequency ablation
© The Author(s) 2016
Published: 28 December 2016
The aim of this study is to research the lesion outline and temperature field in different ways in atrial radiofrequency ablation by using finite element method.
This study used the method which considered the thermal dosage to determine the boundary between viable and dead tissue, and compared to the 50 °C isotherm results in analyzing lesion outline. Besides, we used Hyperbolic equation which considered the relaxation time to calculate the temperature field and contrasted it with Pennes’ bioheat transfer equation.
As the result of the comparison of the lesion outline, when the ablation time was 120 s, the isotherm of the thermal dosage was larger than the 50 °C isotherm and with the increasing of the voltage the gap increased. When the ablation voltage was 30 V, the 50 °C isotherm was larger than the thermal dosage isotherm when the ablation time was less than 160 s. The isotherms overlapped when the time was 160 s. And when the ablation time was more than 160 s, the 50 °C isotherm was less than the thermal dosage isotherm. As to the temperature field, when the ablation voltage was 30 V with the ablation time 120 s the highest temperature decided by Hyperbolic was 0.761 °C higher. The highest temperature changed with relaxation time. In most cases, the highest temperature of the Hyperbolic was higher otherwise the relaxation time was 30–40 s.
It is better to use CEM43 °C to estimate the lesion outline when the ablative time within 160 s. For temperature distribution, the Hyperbolic reflects the influence of heat transmission speed, so the result is more close to the actual situation.
KeywordsAtrial fibrillation Radiofrequency ablation Finite element method Thermal dosage Hyperbolic equation
Atrial fibrillation (AF) is the most common arrhythmia cardiac symptoms, the incidence of it increases with age. It also contacts with some other diseases, such as stroke and heart failure which can degrade the quality of life and increase rates of death . Anti-arrhythmic drugs and surgical operation are the two major therapeutic methods in restoring and maintaining sinus rhythm. However, the incompletely effective and the dangerous side effects of the anti-arrhythmic drugs limit their long-term use. Also, the complexity of the surgical operation limit the popularization of this method . While radiofrequency ablation (RFA) which used electric current to cut the accessory pathways of abnormal tissue and targeted at certain points which produce cardiac arrhythmia. In recent years, due to its advantages of safety, minimally invasive and so on, it was widely used in treatment atrial fibrillation .
In the process of RFA the electromagnetic energy is converted to heat. Tissue temperature above 50 °C are reserved for direct treatment, and the therapy is termed ablation [4, 5]. Dewhirst found that the cell survival/CEM43 relationship closely aligns with isothermal exposure of tissue to temperatures of 50 °C . During RFA 50 °C isotherm which only consider temperature is regarded as the boundary of necrotic cells and survival cells [7–11]. However, many studies show tissue damage is concerned with both temperature and time. On the basis of the conception, thermal dosage which consider both temperature and time was put forward. 43 °C is used as the benchmark temperature how long the tissue to absorb heat in order to keep 43 °C is expressed as cumulative equivalent minutes at 43 °C (CEM43 °C) . So this study created finite element method (FEM) models of cardiac and determined the temperature field in the tissue solving by the Pennes’ bioheat transfer equation [12, 13]. And comparison of the lesion range which determined by the 50 °C isotherm with the lesion range which determined by thermal dosage.
In most of the simulations, Pennes’ bioheat transfer equation which is based on Fourier’s heat theory determined the temperature field. But Fourier’s heat conduction theory is the law of macro-continuity. This theory does not have time item, which implies the speed of heat is infinite. So it also applies in the immediate energy diffusion at the infinite propagation speed in the medium. In most situations, Pennes’ bioheat transfer equation can meet the conditions, but when involved in very low or high temperature, very high heat flux or very short heating duration Fourier’s heat conduction theory breaks down. The reason is the wave nature of heating processes becomes pronounced [14–19]. Hyperbolic equation is based on No-Fourier’s heat theory which thinks about heat wave propaganda speed (relaxation time). In order to consider the condition of the short heat duration time, this studies took Hyperbolic equation into account and contrasted the temperature field which determined by the two equations.
This paper uses the finite element method of cardiac radiofrequency ablation to determine the damage area size. The simulation research, comprehensive analysis of the differences of the temperature field determined by both thermal dose method and 50 °C isotherm. The thermodynamics equation Pennes’ bioheat transfer equation and Hyperbolic equation was also considered. The model was referenced in the paper .
Electromagnetic and Pennes’ bioheat transfer equation
In the other situation, the temperature distribution in the cardiac model was obtained by solving the Hyperbolic equation which considered relaxation time .
Definition of the boundary between the viable and dead tissue
According to the temperature distribution, there are two main transformation algorithms to estimate the ablation region, and to calculate the transverse width, longitudinal length, and the area of the lesion. In this model, we used 50 °C isotherm and CEM43 °C to estimate the thermal damage region.
For tissue temperature above 50 °C are reserved for direct treatment, and the therapy is termed ablation [4, 5]. Dewhirst found that the cell survival/CEM43 relationship closely aligns with isothermal exposure of tissue to temperatures of 50 °C . In RFA, 50 °C is regarded as the boundary of necrotic cells and survival cells, and it is easy and convenient for 50 °C isotherm threshold. So in this model, we use it to analyze.
As there is no critical CEM43 °C for cardiac tissue, we estimate CEM43 °C for myocardium to be 128 min (7680 s) to analyze .
The comparison of the lesion outline decided by the CEM43 °C and decided by 50 °C isotherm were selected for result illustration. In addition, the distribution of temperature field calculated by the Hyperbolic equation and the Pennes’ bioheat transfer equation were both researched.
The comparison of the lesion outline in different voltages
We can see from Fig. 3 that the lesion outline determined by the CEM43 °C was a little smaller than which decided by the 50 °C isotherm. And with the increasing of the voltage, the lesion range and the highest temperature of the two methods increasing. Under the same ablation time, the highest temperature rose about 11 °C with each additional 5 V.
The comparison of the lesion outline at different times
The lesion outline of both methods and the highest temperature increased with the increasing of the time. At the same voltage, the highest temperature rose about 1 °C with every additional 80 s.
The comparison of lesion size determined by isotherm and thermal dosage
The comparison of lesion size determined by the 50 °C isotherm and the CEM43 °C at t = 60 s
Lesion size of 50 °C (mm2)
Lesion size of CEM43 (mm2)
Overestimate rate (%)
The comparison of lesion size determined by the 50 °C isotherm and the CEM43 °C at t = 120 s
Lesion size of 50 °C (mm2)
Lesion size of CEM43 (mm2)
Overestimate rate (%)
The comparison of the temperature field determined by Pennes’ bioheat transfer equation and Hyperbolic equation
Due to the blood flow, the temperature in the blood can hardly raise, the temperature field looks like a semicircle in the myocardium. In above figures, it can be known the temperature increase with both the increase of voltage and the time, and the temperature increase with voltage is more obviously than the increase with time.
In this study, we used 50 °C isotherm and thermal dosage to analyze the lesion outline, and both of them increased with the increase of the voltage and time. From results we can see that the lesion outline of the thermal dosage is smaller than the lesions outline of the 50 °C isotherm in all the time when the voltage increased. The lesion range of CEM43 is smaller than the lesion range of the 50 °C isotherm when the time within 160 s.
In cardiac ablation, the duration time is 60 s to 120 s, and the voltage is 30–60 V in clinical currently. Since we compared the lesion size when the time was 60 and 120 s and the voltage 30–60 V, the size of 50 °C isotherm were larger than the size of the CEM43. Some studies researched in ex vivo liver  and found the two methods give similar results, but the critical isotherm was chosen; and others found that the isotherm overestimate 4.8% for the final lesion diameter than the critical CEM43 °C  since the RF has been shut down. Therefore the thermal dose gets more precise data.
We contrasted the temperature field determined by Pennes’ bioheat transfer equation and Hyperbolic equation. And in this study, we assumed the relaxation time was 16 s to analyze at first. And the results showed the shape of the temperature distribution similarly. We also chose 10 values from the relaxation time from 10 to 50 s to simulate, and we found that only the maximum temperature which determined by the Hyperbolic equation is little higher in addition to the relaxation time between 30 and 40 s. The temperature determined by the Hyperbolic equation was smaller than which determined by the Pennes equation in the start and with the growth of the time it become larger than the temperature of Pennes. The highest temperature determined by the Hyperbolic equation changes up and down in a small range which is smaller than 2 °C.
In pre-ablation, due to the theory of the Hyperbolic considers the limitation of the transfer of the energy heat, at the beginning of the ablation, the heat that produced by the myocardium cannot deliver in the layer and transfer through the surface blood convection. Thus the lesion area was smaller than the Pennes transfer method. As the ablation time was over the relaxation time, the electric field produced more heat, and the heat is trapped inside the myocardium, so the area was bigger than the Pennes in post ablation.
Compared to 50 °C isotherm, CEM43 °C considers the temperature history, it can be more accuracy to describe the change of temperature. And from above we can know that the 50 °C isotherm overestimate the thermal dose when the ablative time within 160 s. Now the cardiac ablation time is 60–120 s, so use thermal dose to analyze can get more accurate data and its necessary for treatment.
Pennes’ bioheat transfer equation does not consider the speed of the heat while Hyperbolic equation take it into account, and the Hyperbolic method reflects the influence of heat propagation speed. In this study, we contrasted the maximum temperature in different relaxation time and found that the highest temperature changing over it. The result of the Hyperbolic is more close to the actual situation of cardiac radiofrequency ablation, so the latter can describe temperature more exactly in cardiac ablation.
There are still some limits in this research. As lack of a time-varying blood perfusion data in clinic, we did not consider this condition into the simulation. In future, we’ll take this condition into account and to get more accuracy data to compare the results.
ZT was responsible for the design, data collection and overall investigation. ZT, QN and XHN were responsible for computational modeling and data analysis part. TD and RRW were responsible for the statistical analysis part. All authors (1) have made substantial contributions to conception and design, or acquisition of data, or analysis and interpretation of data; (2) have been involved in drafting the manuscript or revising it critically for important intellectual content; and (3) have given final approval of the version to be published. Each author has participated sufficiently in the work to take public responsibility for appropriate portions of the content. All authors read and approved the final manuscript.
The work was supported by Natural Science Foundation of Beijing (3162006), National Natural Science Foundation of China (No. 31070754), Education Project Scientific and Technological Program of Beijing Municipal Commission (KM201410005028), the Importation and Development of High-Caliber Talents Project of Beijing Municipal Institutions and Basic Research Foundation of Beijing University of Technology (X4015999201401).
Other than the grants listed in the acknowledgement section, the authors declare that they have no other competing interests.
About this supplement
This article has been published as part of BioMedical Engineering OnLine Volume 15 Supplement 2, 2016. Computational and experimental methods for biological research: cardiovascular diseases and beyond. The full contents of the supplement are available online http://biomedical-engineering-online.biomedcentral.com/articles/supplements/volume-15-supplement-2.
Availability of data and materials
All data are fully available without restriction.
Publication of this article was funded by Natural Science Foundation of Beijing (3162006).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Camm AJ, Kirchhof P, Lip GY, et al. Guidelines for the management of atrial fibrillation: the task force for the management of atrial fibrillation of the european society of cardiology (ESC). Eur Heart J. 2010;31:2369–429.View ArticleGoogle Scholar
- Dong J, Liu X, Ma C. Catheter ablation of atrial fibrillation. J Clin Cardiovasc Dis. 2010;26(9):641–3.Google Scholar
- Berjano E. Theoretical modeling for radiofrequency ablation: state-of-the-art and challenges for the future. Biomed Eng Online. 2006;5:24.View ArticleGoogle Scholar
- Eyring H. The energy of activation for bimolecular reactions involving hydrogen and the halogens, according to the quantum mechanics. J Am Chem Soc. 1931;53(2):2537–49.View ArticleGoogle Scholar
- Eyring H. The activated complex in chemical reactions. J Chem Phys. 1935;3(2):107–15.View ArticleGoogle Scholar
- Viglianti BL, Dewhirst MW, Abraham JP, et al. Rationalization of thermal injury quantification methods: application to skin burns. Burns J Int Soc Burn Inj. 2014;40(5):896–902.View ArticleGoogle Scholar
- Haemmerich D, Tungjitkusolmun S, Staelin T, et al. Finite-element analysis of hepatic multiple probe radio-frequency ablation. IEEE Trans Biomed Eng. 2002;49(8):836–42.View ArticleGoogle Scholar
- Jain MK, Wolf PD. A three-dimensional finite element model of radiofrequency ablation with blood flow and its experimental validation. Ann Biomed Eng. 2000;28(9):1075–84.View ArticleGoogle Scholar
- Tungjitkusolmun S, Staelin T, Haemmerich D, et al. Three-dimensional finite-element analyses for radio-frequency hepatic tumor ablation. IEEE Trans Biomed Eng. 2002;49(1):3–9.View ArticleGoogle Scholar
- Tungjitkusolmun S, Haemmerich D, Cao H, et al. Modeling bipolar phase-shifted multielectrode catheter ablation. IEEE Trans Biomed Eng. 2002;49(1):10–7.View ArticleGoogle Scholar
- Haemmerich D, Webster J G, Mahvi D M: Thermal dose versus isotherm as lesion boundary estimator for cardiac and hepatic radio-frequency ablation. Engineering in medicine and biology society, 2003. In: Proceedings of the 25th annual international conference of the IEEE, vol 1. New York: IEEE; 2003. p. 134–7.Google Scholar
- Diederich CJ. Thermal ablation and high-temperature thermal therapy: overview of technology and clinical implementation. Int J Hyperth. 2005;21(8):745–53.View ArticleGoogle Scholar
- Zhang H, Nan Q, Liu Y. Thermal distribution of microwave antenna for atrial fibrillation catheter ablation. Int J Hyperth. 2013;29(6):582–9.View ArticleGoogle Scholar
- Vernotte P. Les Paradoxes de la Théorie Continue de l’équation de lachaleur. Comptes Rendus. 1958:246;3154–5.MATHGoogle Scholar
- Vedavarz A, Mitra K, Kumar S. Hyperbolic temperature profiles for laser surface interactions. J Appl Phys. 1994;76(9):5014–21.View ArticleGoogle Scholar
- Mitra K, Kumar S, Vedavarz A, et al. Experimental evidence of Hyperbolic heat conduction in processed meat. J Heat Transfer. 1995;117(3):568–73.View ArticleGoogle Scholar
- Shih TC, Kou HS, Liauh CT, et al. The impact of thermal wave characteristics on thermal dose distribution during thermal therapy: a numerical study. Med Phys. 2005;32(9):3029–36.View ArticleGoogle Scholar
- Sieniutycz S. The variational principles of classical type for non-coupled non-stationary irreversible transport processes with convective motion and relaxation. Int J Heat Mass Transf. 1977;20(11):1221–31.View ArticleMATHGoogle Scholar
- Kaminski W. Hyperbolic heat conduction equation for materials with a nonhomogeneous inner structure. J Heat Transf. 1990;112(3):555–60.View ArticleGoogle Scholar
- Liu P. Finite element method simulation study of lesion size in cardiac catheter radiofrequency ablation for atrial fibrillation. Hangzhou: Zhejiang University; 2012.Google Scholar
- Guo X, Zhai F, Nan Q. The temperature field simulation of radiofrequency catheter-based renal sympathetic denervation for resistant hypertension. Bio Med Mater Eng. 2014;24(1):315–21.Google Scholar
- Liu P, Liu J, Duan H. Finite element modeling simulation of cardiac radiofrequency catheter ablation based on heart chamber anatomy. Chin J Biomed Eng. 2012;31(4):532–8.Google Scholar
- Panescu D, Whayne JG, Fleischman SD, et al. Three-dimensional finite element analysis of current density and temperature distributions during radio-frequency ablation. IEEE Trans Biomed Eng. 1995;42(9):879–90.View ArticleGoogle Scholar
- Tungjitkusolmun S, Cao H, Tsai J Z, et al: Using ANSYS for three-dimensional electrical-thermal models for radio-frequency catheter ablation. Engineering in medicine and biology society, 1997. In: Proceedings of the international conference of the IEEE; 1997. p.161–4.Google Scholar
- Incropera FP. Fundamentals of heat and mass transfer. New York: Wiley; 1985. p. 139–62.Google Scholar
- Sparrow EM, Abraham JP. A simulation of gas-based, endometrial-ablation therapy for the treatment of menorrhagia. Ann Biomed Eng. 2008;36(1):171–83.View ArticleGoogle Scholar
- Abraham JP, Sparrow EM, Ramadhyani S. Numerical simulation of a BPH thermal therapy—a case study involving TUMT. J Biomech Eng. 2007;129:548–57.View ArticleGoogle Scholar
- Abraham JP, Sparrow EM. A thermal ablation model including liquid-to-vapor phase change, necrosis-dependent perfusion, and moisture-dependent properties. Int J Heat Mass Transf. 2007;50:2537–44.View ArticleMATHGoogle Scholar
- Sapareto SA, Dewey WC. Thermal dose determination in cancer therapy. Int J Radiat Oncol Biol Phys. 1984;10(6):787–800.View ArticleGoogle Scholar
- .Li F, Liu L, Gong X, et al: Research on thermal dosage for high intensity focused ultrasound in bovine liver based on the magnetic resonance thermometry. International conference on biomedical engineering and informatics, BMEI 2009, October 17–19, 2009, Tianjin; 2009. p. 1–3.Google Scholar
- Graham SJ, Chen L, Leitch M, et al. Quantifying tissue damage due to focused ultrasound heating observed by MRI. Magn Reson Med. 1999;41(2):321–8.View ArticleGoogle Scholar