Effect of variable heat transfer coefficient on tissue temperature next to a large vessel during radiofrequency tumor ablation
© dos Santos et al; licensee BioMed Central Ltd. 2008
Received: 26 February 2008
Accepted: 11 July 2008
Published: 11 July 2008
One of the current shortcomings of radiofrequency (RF) tumor ablation is its limited performance in regions close to large blood vessels, resulting in high recurrence rates at these locations. Computer models have been used to determine tissue temperatures during tumor ablation procedures. To simulate large vessels, either constant wall temperature or constant convective heat transfer coefficient (h) have been assumed at the vessel surface to simulate convection. However, the actual distribution of the temperature on the vessel wall is non-uniform and time-varying, and this feature makes the convective coefficient variable.
This paper presents a realistic time-varying model in which h is a function of the temperature distribution at the vessel wall. The finite-element method (FEM) was employed in order to model RF hepatic ablation. Two geometrical configurations were investigated. The RF electrode was placed at distances of 1 and 5 mm from a large vessel (10 mm diameter).
When the ablation procedure takes longer than 1–2 min, the attained coagulation zone obtained with both time-varying h and constant h does not differ significantly. However, for short duration ablation (5–10 s) and when the electrode is 1 mm away from the vessel, the use of constant h can lead to errors as high as 20% in the estimation of the coagulation zone.
For tumor ablation procedures typically lasting at least 5 min, this study shows that modeling the heat sink effect of large vessels by applying constant h as a boundary condition will yield precise results while reducing computational complexity. However, for other thermal therapies with shorter treatment using a time-varying h may be necessary.
Radiofrequency (RF) tumor ablation is a treatment modality that uses radiofrequency electric current in an attempt to destroy cancer cells by localized heating of tumors. RF ablation is used for liver cancer, with increasing use in other organs such as kidney, lung, bone, and adrenal gland .
The heat due to RF heating causes tissue necrosis at predictable temperatures in relatively predictable volumes. During an RF ablation procedure, an active electrode is inserted percutaneously (i.e. through a small incision in the skin), during laparoscopy, or during open surgery with image-guidance into the tumors. Grounding pads are positioned at the patient's thighs or back muscles. RF energy is applied and current flows from the active electrode to the grounding pads. Thus, the patient becomes an element of the electrical circuit.
High RF current densities around the electrode results in resistive heating in surrounding tissue. While it takes several hours to induce cell necrosis at 43°C , at 50°C cell death occurs within 2–3 min . The tissue adjacent to the electrode is rapidly heated, and the remainder of the tissue is heated by thermal conduction, which is a slow process. Typically, a single ablation takes 12 – 35 min depending on device type . If the tissue temperature near the electrode is excessive, the tissue desiccates, and this process results in an electrically insulating layer preventing further energy deposition. As a result, the phenomenon limits the volume of the tissue that can be treated. A difference between tumor ablation and hyperthermia treatments (where lower temperatures of 43 – 45°C are used) is that during ablation most of the tissue volume that undergoes necrosis is at temperatures between 50 and 100°C. Therefore, the prevalent effect of cell death is coagulation necrosis conversely to hyperthermia treatments where a number of biologically complex mechanisms are in effect .
A shortcoming of current RF ablation devices is the limited performance adjacent to large blood vessels (diameter > 3 mm). When the active electrode is inserted near large vessels, the blood flow drags thermal energy away from the target tissue . This is a heat sink effect that can change both the shape and maximum volume that can be treated. In fact, the distance of the blood vessels from the tumor determines the location of the maximal tissue temperature. As a result, tumors in the vicinities of large vessels are associated with high recurrence rates . A considerable number of mathematical models have been suggested to describe heat transfer between tissue and vasculature [8–20] most of them are not directly applicable to liver tissue due to the specific blood supply of the liver. Most of the blood perfusing the liver (roughly 70%) is venous blood ; in addition, there are no counter-current vessels (i.e. venous and arterial vessels adjacent with opposing blood flow) in the liver as present in other tissues, and counter-current vessels are an integral assumption of most heat transfer models. As a result, most computational models of high-temperature tumor ablation simulate microvascular perfusion using the commonly used Pennes formulation , and model large vessels separately since the Pennes model does not describe large vessel perfusion accurately, similar to a previous study .
In order to estimate the heat sink effect of large vessels, many simulations performed hitherto assigned either constant temperature at the surface of the vessel [6, 23] or constant convective heat transfer coefficient, h, throughout the RF ablation procedure . In order to estimate the value of h, the simulations performed implied that the flow in the vessel is laminar, the thermal boundary layer is fully-developed and h remains constant throughout the procedure. Other works support that hypothesis that the flow is laminar in large vessels [25, 26]. However, since the heated region during the procedure varies in time and the vessel heated length is small, the thermal boundary layer is not fully-developed. Thus, both the vessel wall temperature and h varies during the ablation procedure and cannot be considered constant . Since the magnitude of h is a parameter that may significantly impact on the shape and size of the coagulation zone obtained during hepatic RF ablation, it is important to accurately model the time-varying behavior of h during the RF ablation of liver tumors in order to correctly determine the coagulation zone size. Thus, the overall objective of this work is to theoretically evaluate the impact of the time-varying h on the accuracy of the simulations of RF tumor ablation procedures.
Methods and models
where ρ, c and k are, respectively, the density (kg*m-3), the specific heat (J*kg-1*K-1) and the thermal conductivity (W*m-1 K-1) of the liver tissue. (A*m-2) and (V*m-1) are current density and the electric field intensity and can be calculated with Laplace's equation, using a quasi-static approximation. T is the temperature of the tissue, T bl is the temperature of the blood, ρ bl is the blood density, c bl is specific heat of the blood and w bl is blood perfusion (s-1). Q m (W*m-3) is the energy generated by metabolic processes and was neglected since it is small if compared to the other terms . Microvascular perfusion was included in the model, considering the Pennes model . Although the perfusion in normal and tumorous tissue can vary greatly, this paper focus on the temperature field close to large vessels and on the behavior of the time-varying h during RF ablation rather than simulating a clinical situation. The parameter w bl used in this model was 6.4.10-3 s-1 , which is in the range of the perfusion for cirrhotic human liver tissue . Equation 1 states that RF current flowing through the tissue is converted into thermal energy, which in turn cause tissue injury.
Thermal and electric properties of the materials.
Finite element region
where T t is the electrode tip temperature, T s is the set temperature measured at the tip of the electrode where maximum temperatures occur, and K p and K i are control parameters of the PI controller. The control parameters determine the behavior of the algorithm, e.g., response time, overshoot, swinging. The parameters we used for the PI controller were K p = 0:4 and K i = 0:04. These values were obtained based on tests performed with the dynamic system. In order to implement the time integral of the temperature difference between the electrode tip and the set temperature (equations 3 and 4), we added an ordinary differential equation in Femlab.
In this study, the maximum hepatic tissue temperature was kept at 90°C and the maximum power applied during the procedure was 28 W. We used the 50°C-isotherm to determine the coagulation zone boundary as has been done in previous models of RF tumor ablation [6, 24, 28, 37]. We used the temperature at the tip of the electrode for convergence test. We refined the mesh and run the analysis again. We searched the optimal mesh size in order to keep the simulation fast while keeping the accuracy better than 0.01°C at the tip of the electrode. We tested 2 mesh sizes: 65,759 and 13,714. We concluded that the 13,714-element model satisfied our requirement. The time step size was chosen at 0.05 s, so that the maximum temperature change was smaller than 0.01°C during each step. The simulations were performed on a PC with a 2.4 GHz PENTIUM Celeron CPU, with 1 GB of RAM and 30 GB of hard disk space. For post processing, we employed the built-in module in FEMLAB and MATLAB. The computation time required for 10 min simulation was about 3 hours.
We placed the electrode in two different configurations: (a) 1 mm away from the vessel (shown in Figure 1) and (b) 5 mm away from the vessel. For each of these two configurations, we performed three studies.
The behavior of the convective heat coefficient on the vessel;
The evolution of the coagulation zone volume when the time-varying behavior of h is taken into account and when h is constant and equal to the value of the time-varying h at the end of the 10 minutes ablation;
The behavior of the maximum temperature at the vessel wall with both constant and time-varying h.
The behavior of the convective heat coefficient at the vessel wall
The evolution of the coagulation zone volume for varying h and for constant h cases
The behavior of the maximum temperature at the vessel wall for varying h and for constant h
The aim of this study was to evaluate the impact of the time-varying behavior of the convective heat transfer coefficient h for the determination of the coagulation zone during RF tumor ablation procedures. We have simulated two different behaviors of the convective heat transfer coefficient: constant and variable. We also studied two configurations with the electrode either 1 mm, or 5 mm distant from a large vessel 10 mm in diameter . When the electrode is 1 mm away from the vessel, h initially rapidly increases up to 6000 W*m-2 *K-1 followed by a sharp decrease to approximately 554 W*m-2 *K-1 within 60 s (Figure 3). When the electrode is 5 mm away from the vessel, the convective coefficient rapidly increases to approximately 7000 W*m-2 *K-1 and then rapidly decreases to 336 W*m-2 *K-1 (Figure 4). This behavior occurs because the thermal energy reaches the vessel rapidly due to close proximity to the electrode. Thus, h initially increases due to the increase in temperature on the vessel wall, and after a while the temperature increase extends along the vessel and the convective heat transfer coefficient decreases. Hence, the maximum value of h for the first case occurs around 0.2 s, and for the second case it occurs around 0.5 s after start of the ablation. This behavior becomes more obvious by observing Equation 2. Initially, the vessel wall is not heated yet, and h is 0 W*m-2 *K-1 because there is no heat transfer to the vessel. Then the temperature at the vessel starts to increase and also the heated length grows larger typically up to 11 cm. Initially the heated length is very small compared to the temperature difference T 1 - T 0, thus h sharply increases. After a few seconds, the length starts to increase faster than the difference (T 1 - T 0) and h starts to decrease until both (T 1 - T 0) and h show little change because equilibrium between heat loss and heating is reached. Recall that the thermal boundary layer is very small at the entrance region. Thus, h is very high . It is noteworthy that if one wants to measure the dynamic behavior of the heat convection coefficient, the instrumentation must have a fast dynamic response and a very high static range as one can see from figures 3 and 4. We also analyzed the coagulation zone for the two cases of time-varying h and constant h. In the first few seconds and when the electrode is close to the vessel, the time-varying h has considerably higher values than the constant h and the relative error in coagulation zone volume due to the assumption of a constant h is very high (22%). For RF liver ablation this may not be important since initially the ablation volume is small (with a small absolute error), and the volume near the beginning of the treatment is not of high clinical value. As the time increases, the time-varying h sharply decreases, approaching its final value after 80 s for the 1 mm case and after 120 s for 5 mm case. Thus, the final coagulation zone volume is nearly identical in both cases. We also investigated the maximum temperature at the wall of the vessel, where tumor recurrence is possible due to insufficient temperatures. The maximum error in temperature is around 6% for the first 50–100 s, but then again is close to zero at the end of the ablation. Consequently, for tumor ablation procedures typically lasting at least 5 min, a very important result is the fact that modeling the heat sink effect of large vessels by applying constant h as a boundary condition will yield accurate results for RF ablation while reducing computational complexity. For other thermal therapies with shorter treatment times (e.g. 45 – 60 s long cardiac ablation close to coronary vessels) using a time-varying h might be required but needs to be further investigated.
Previous studies considered a constant heat transfer coefficient throughout the ablation procedure. In this work, simulations were performed using a more realistic, time-varying analytical expression of the convective heat transfer coefficient, which depends on the blood velocity and on the temperature distribution on the vessel wall. The simulations showed that the assumption of a constant convective coefficient leads to precise results when it is used for typical ablation procedures. Only during the first 1–2 min, a time-varying coefficient produces noticeable different results. However, this has no clinical impact for RF liver ablation procedure, which typically takes over 5 min.
We would like to thank the University of Brasilia, CNPq, CAPES and FAPDF for providing support for this research. Part of this work was conducted in a facility constructed with support from the National Institutes of Health, Grant Number C06 RR018823 from the Extramural Research Facilities Program of the National Center for Research Resources.
- Neeman Z, Wood BJ: Radiofrequency ablation beyond the liver. Tech Vasc Interv Radiol 2002, 5: 156–163. 10.1053/tvir.2002.36419View ArticleGoogle Scholar
- Dewhirst MW, Viglianti BL, Lora-Michiels M, Hanson M, Hoopes PJ: Basic principles of thermal dosimetry and thermal thresholds for tissue damage from hyperthermia. Int J Hyperthermia 2003,19(3):267–294. 10.1080/0265673031000119006View ArticleGoogle Scholar
- Graham SJ, Chen L, Leitch M, Peters RD, Bronskill MJ, Foster FS, Henkelman RM, Plewes DB: Quantifying tissue damage due to focused ultrasound heating observed by MRI. Int J Hyperthermia 1999,41(2):321–328.Google Scholar
- Pereira PL, Trubenbach J, Schenk M, Subke J, Kroeber S, Remy CT, Schmidt D, Brieger J, Claussen CD: Radiofrequency ablation: in vivo comparison of four commercially available devices in pig livers. Radiology 2004,232(2):482–90. 10.1148/radiol.2322030184View ArticleGoogle Scholar
- Dewhirst MW, Vujaskovic Z, Jones E, Thrall D: Re-setting the biologic rationale for thermal therapy. Int J Hyperthermia 2005, 21: 779–790. 10.1080/02656730500271668View ArticleGoogle Scholar
- Tungjiktkusolmun S, Staelin S, Haemmerich D, Tsai JZ, Cao H, Webster JG, LFT Jr, Mahvi DM, Voerperian V: Three-Dimensional Finite-Element Analyses for Radio-Frequency Hepatic Tumor Ablation. IEEE Trans Biomed Eng 2002, 49: 3–9. 10.1109/10.972834View ArticleGoogle Scholar
- Lu DS, Raman SS, Vodopich DJ, Wang W, Sayre J, Lassman C: Effect of vessel size on creation of hepatic radiofrequency lesions in pigs: assessment of the heat sink effect. Am J Roentgenol 2002, 178: 47–51.View ArticleGoogle Scholar
- Pennes HH: Analysis of tissue and arterial blood temperature in the resting human forearm. J Appl Physiol 1948, 1: 93–122.Google Scholar
- Chen MM, Holmes KR: Microvascular contributions in tissue heat transfer. Ann New York Acad Sci 1980, 355: 137–150. 10.1111/j.1749-6632.1980.tb50742.xView ArticleGoogle Scholar
- Osman MM, Afify EM: Thermal modeling of the normal women's breast. J Biomech Eng 1984, 106: 123–130.View ArticleGoogle Scholar
- Weinbaum S, Jiji LM: A new simplified bioheat equation for the effect of blood flow on local average tissue temperature. J Biomech Eng 1985, 107: 131–139.View ArticleGoogle Scholar
- Weinbaum S, Xu LX, Zhu L, Ekpene A: A new fundamental bioheat equation for muscle tissue: Part 1. Blood perfusion term. J Biomech Eng 1997, 119: 278–288. 10.1115/1.2796092View ArticleGoogle Scholar
- Xuan Y, Roetzel W: Bioheat equation of the human thermal system. Chem Eng Technol 1997, 20: 268–276. 10.1002/ceat.270200407View ArticleGoogle Scholar
- Roemer RB, Dutton AW: A generic tissue convective energy balance equation. Part 1. Theory and derivation. J Biomech Eng 1998, 120: 395–404. 10.1115/1.2798007View ArticleGoogle Scholar
- Wren J, Karlsson M, Loyd D: A hybrid equation for simulation of perfused tissue during thermal treatment. Int J Hyperthermia 2001, 17: 483–498. 10.1080/02656730110081794View ArticleGoogle Scholar
- Deng ZS, Liu J: Blood perfusion-based model for characterizing the temperature fluctuation in living tissues. Physica A 2001, 300: 521–530. 10.1016/S0378-4371(01)00373-9View ArticleGoogle Scholar
- Shih TC, Kou HS, Lin W: Effect of effective tissue conductivity on thermal dose distributions of living tissue with directional blood flow during thermal therapy. Int Commun Heat Mass Transfer 2002, 29: 115–126. 10.1016/S0735-1933(01)00330-XView ArticleGoogle Scholar
- Zhu L, Xu LX, He Q, Weinabum S: A new fundamental bioheat equation for muscle tissue: Part II. Temperature of SAV vessels. J Biomech Eng 2002, 124: 121–132. 10.1115/1.1431263View ArticleGoogle Scholar
- Khou HS, Shih TC, Lin WL: Effect of the directional blood flow on thermal dose distribution during thermal therapy: an application of a Green's function based on the porous model. Phys Med Biol 2003, 48: 1577–1589. 10.1088/0031-9155/48/11/307View ArticleGoogle Scholar
- Shrivastava D, Roemer R: An analytical study of Poisson conduction shape factors for two thermally significant vessels in a finite, heated tissue. Phys Med Biol 2005,50(15):3627–3641. 10.1088/0031-9155/50/15/010View ArticleGoogle Scholar
- Martini FH, Timmons MJ, McKinley MP: Human Anatomy. Upper Saddle River, New Jersey 07458: Prentice-Hall Inc; 2000.Google Scholar
- Rawnsley R, Roemer R, Dutton A: The simulation of discrete vessel effects in experimental hyperthermia. J Biomech Eng 1994,116(3):256–262. 10.1115/1.2895728View ArticleGoogle Scholar
- Stanczyk M, Leeuwen GMJV, V SAA: Discrete vessel heat transfer in perfused tissue-model comparison. Phys Med Biol 2007, 52: 2379–2391. 10.1088/0031-9155/52/9/004View ArticleGoogle Scholar
- Haemmerich D, Wright AE, Mahvi DM, Lee FT Jr, Webster JG: Hepatic bipolar radiofrequency ablation creates coagulation zones close to blood vessels: A finite element study. Med Biol Eng Comput 2003, 41: 317–323. 10.1007/BF02348437View ArticleGoogle Scholar
- Gates G, Dore E: Streamline flow in the human portal vein. J Nucl Med 1973, 14: 79–83.Google Scholar
- Garcier JM, Bousquet J, Alexandre M, Filaire M, Viallet JF, Vanneuville G, Boyer L: Visualisation of the portal flows by porto-scanner. Surg Radiol Anat 2000, 22: 239–242. 10.1007/s00276-000-0239-4View ArticleGoogle Scholar
- Consiglieri L, dos Santos I, Haemmerich D: Theoretical analysis of the heat convection coefficient in large vessels and the significance for thermal ablative therapies. Phys Med Biol 2003, 48: 4125–4134. 10.1088/0031-9155/48/24/010View ArticleGoogle Scholar
- Berjano EJ: Theoretical modeling for radiofrequency ablation: state-of-the-art and challenges for the future. Biomed Eng Online 2006, 5: 24. 10.1186/1475-925X-5-24View ArticleGoogle Scholar
- Ebbini ES, Umemura SI, Ibbini M, Cain C: A cylindrical- section ultrasound phased-array applicator for hyperthermia cancer therapy. IEEE Trans. Ultrasonics, Ferroelectrics, and Frequency Control 1988, 35: 561–572. 10.1109/58.8034View ArticleGoogle Scholar
- Van Beers BE, Leconte I, Smith AM, Jamart J, Horsmans Y: Hepatic perfusion parameters in chronic liver disease: dynamic CT measurements correlated with disease severity. AJR Am J Roentgenol 2001, 176: 667–673.View ArticleGoogle Scholar
- Valvano JW, Cochran JR, Diller KR: Thermal conductivity and diffusivity of biomaterials measured with self-heating thermistors. Int J Thermophys 1985, 6: 301–311. 10.1007/BF00522151View ArticleGoogle Scholar
- Panescu D, Whayne J, Fleschman SD, Mirotznik MS, Swanson DK, Webster JG: Three-dimensional finite-element analysis of current density and temperature distributions during radio-frequency ablation. IEEE Trans Biomed Eng 1995, 42: 879–890. 10.1109/10.412649View ArticleGoogle Scholar
- Tungjitkusolmun S, Cao H, Tsai J, Webster JG: Using Ansys for three-dimensional electrical-thermal models for radio-frequency catheter ablation. Proc 19th Ann Int Conf IEEE 1997, 1: 161–164.Google Scholar
- Duck FA: Physical properties of tissue. London: Academic Press; 1990.Google Scholar
- Miguel AF: An instrument to measure the convective heat coefficient on the hepatic artery and on the portal vein (in portuguese). In Master's thesis. University of Brasilia; 2006.Google Scholar
- Haemmerich D, Webster JG, Limanond P, Aziz D, Economou J, Busuttil R, Sayre J: Automatic control of finite element models for temperature-controlled radiofrequency ablation. Biomed Eng Online 2005, 4: 1–8. 10.1186/1475-925X-4-42View ArticleGoogle Scholar
- Liu Z, Ahmed M, Sabir A, Humphries S, Goldberg SN: Computer modeling of the effect of perfusion on heating patterns in radiofrequency tumor ablation. Int J Hyperthermia 2007, 23: 49–58. 10.1080/02656730601094415View ArticleGoogle Scholar
- Gray S: Gray's Anatomy. New York: Vintage Books; 1990.Google Scholar
- Incropera FP, DeWitt DP: Fundamentals of Heat and Mass Transfer. New York, USA: John Willey & Sons, Inc; 1996.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.