 Research
 Open Access
 Published:
Effect of variable heat transfer coefficient on tissue temperature next to a large vessel during radiofrequency tumor ablation
BioMedical Engineering OnLine volume 7, Article number: 21 (2008)
Abstract
Background
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 nonuniform and timevarying, and this feature makes the convective coefficient variable.
Methods
This paper presents a realistic timevarying model in which h is a function of the temperature distribution at the vessel wall. The finiteelement 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).
Results
When the ablation procedure takes longer than 1–2 min, the attained coagulation zone obtained with both timevarying 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.
Conclusion
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 timevarying h may be necessary.
Background
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 [1].
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 imageguidance 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 [2], at 50°C cell death occurs within 2–3 min [3]. 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 [4]. 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 [5].
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 [6]. 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 [7]. 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 [21]; in addition, there are no countercurrent vessels (i.e. venous and arterial vessels adjacent with opposing blood flow) in the liver as present in other tissues, and countercurrent vessels are an integral assumption of most heat transfer models. As a result, most computational models of hightemperature tumor ablation simulate microvascular perfusion using the commonly used Pennes formulation [8], and model large vessels separately since the Pennes model does not describe large vessel perfusion accurately, similar to a previous study [22].
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 [24]. 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 fullydeveloped 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 fullydeveloped. Thus, both the vessel wall temperature and h varies during the ablation procedure and cannot be considered constant [27]. 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 timevarying 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 timevarying h on the accuracy of the simulations of RF tumor ablation procedures.
Methods and models
We used the software FEMLAB (v. 3.2) to generate the finite element model. This software performs coupled electricalthermal field analysis, and provides all the elements needed to build the model, solve the problem and postprocess the results. The heating of tissue during RF ablation is modeled by the bioheat equation (1) [8].
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. $\overrightarrow{J}$(A*m^{2}) and $\overrightarrow{E}$(V*m^{1}) are current density and the electric field intensity and can be calculated with Laplace's equation, using a quasistatic 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 [28]. Microvascular perfusion was included in the model, considering the Pennes model [8]. 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 timevarying 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} [29], which is in the range of the perfusion for cirrhotic human liver tissue [30]. Equation 1 states that RF current flowing through the tissue is converted into thermal energy, which in turn cause tissue injury.
The thermal and electrical properties for the various materials of interest are listed in Table 1[6, 31–33]. The vessel wall was not modeled separately since it has similar properties as normal liver tissue [34]. Since we are mainly concerned with the impact of h on the thermal ablation zone, this study does not model the temperature dependence of the thermalelectrical parameters of the tissue.
A three dimensional view of the model and its physical domain are shown in Figure 1(a). Figure 1(b) shows the top view of the model. We set the temperature on the boundary of the model to 37°C and the voltage at outer surfaces at 0 V. We modeled a multitine electrode (15 gauge, model30, Rita Medical Systems, Inc., Fremont, CA).
The hepatic tissue where the electrode is embedded was modeled as a cylinder (100 mm diameter × 120 mm length). The initial temperature of tissueelectrode system and of the boundary of the model was set to 37°C. The voltages on the outer surfaces of the model were 0 V, assuming the Dirichlet boundary condition. The convective boundary condition was considered in the tissuevessel interface. Thus, the heat sink effect of the blood flow inside the vessel was modeled by a convective boundary condition at the tissuevessel interface. The value of convective coefficient was set in the model according to Equation 2, which assumes that the blood behaves as a Newtonian fluid in large vessels, the flow is laminar and the geometry of the vessel is linear [27]. Subsequently, with the use of Equation 2, the value of the convective coefficient is updated at each time step and is set as the new convective boundary condition on the vessel wall for the next step.
where h is the convective heat transfer coefficient (W*m^{2} *K^{1}) and accounts for the large vessel perfusion, ρ is the density of the blood (kg*m^{3}), c is the specific heat of the blood (J*kg^{1} *K^{1}), T _{0} is the inlet blood temperature (37°C), T _{1} is the temperature at the vessel outlet and T _{ s }is the maximum temperature of the heated region at the vessel wall. The parameters a and b are, respectively, the length of the heated region and the diameter of the vessel (Figure 2). u is the constant mean value of the blood velocity (0.2 m*s^{1}), considering the blood flow with laminar profile. Also, T _{ m }is 0:5·(T _{0} + T _{1}). The temperatures T _{0} and T _{1} are set at 37 and 37.05°C when the temperature of the vessel wall is 50°C. This choice is based on both theoretical and experimental work performed by our group [27, 35]. For other vessel wall temperatures, we used linear interpolation to evaluate the temperature T _{1}. The temperature T _{0} is kept constant at 37°C. The length of the heated region was calculated based on regions where the temperature in the vessel is above 37.01°C. Equation 2 considers the blood as a Newtonian fluid and one might be cautious about extending this equation to vessels smaller than 0.5 mm in diameter [27].
Initially, the profile temperature is uniform and there is no heated region in the vessel. Thus, the convective heat transfer coefficient is set initially at 0 W*m^{2} *K^{1}. Temperature controlled RF ablations were simulated by varying the voltage applied to the electrodes during the procedure. The voltage was controlled in order to keep the maximum tissue temperature at 90°C. We performed temperature control by means of a proportionalintegral controller as described previously [36]. The control algorithm for the finite element model is described by equations 3 and 4.
e = T _{ s } T _{ t }
(4)
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°Cisotherm 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,714element 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 builtin 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 timevarying behavior of h is taken into account and when h is constant and equal to the value of the timevarying h at the end of the 10 minutes ablation;

The behavior of the maximum temperature at the vessel wall with both constant and timevarying h.
Results
The behavior of the convective heat coefficient at the vessel wall
Figure 3(a) shows the evolution of the convective heat transfer coefficient at the vessel wall when the electrode is placed 1 mm away from the vessel. One can see in Figure 3(b) that h increases sharply at the beginning of the ablation procedure, followed by a sharp decrease within a second. The final value of h is indicated with a dotted line in Figure 3.
Figure 4(a) shows the evolution of the convective heat transfer coefficient on the vessel wall when the electrode is located 5 mm away from the vessel. Figure 4(b) shows that the maximum value of h occurs approximately at 0.5 s and decreases sharply thereafter. The final value of h is shown in Figure 4(a) and 4(b).
The evolution of the coagulation zone volume for varying h and for constant h cases
We analyzed the evolution of the coagulation zone volume (temperature greater than 50°C). First, we analyzed the coagulation zone development when varying h was taken into account. Second, we considered h as constant and equal to the final value of h in the first simulation. Figure 5(a) shows the coagulation zone volume development for both constant and timevarying h when the electrode is 1 mm away from the vessel. The absolute error is very small so that the curves are almost identical. Figure 5(b) shows the relative error (%) between the results obtained with the constant h and timevarying h. We also simulated the volume coagulation zone when the electrode is 5 mm away from the vessel. However, the differences in coagulation zones for constant h and timevarying h were below to 0.5% and therefore considered negligible.
The behavior of the maximum temperature at the vessel wall for varying h and for constant h
We analyzed the evolution of the temperature at the surface of the vessel. We determined the maximum temperature at the surface when the distance between the electrode and the vessel is 1 and 5 mm. Figure 6(a) shows the maximum temperature at the vessel wall when the electrode is 1 mm away from the vessel. Figure 6(b) shows the relative error (%). Figure 6(a) and 6(b) reveals that the difference in temperature is negligible after 80 s since the relative error approaches to 0%.
Figure 7(a) shows the temperature when the electrode is 5 mm from the vessel. Figure 7(b) shows the percentage error. One can see that the difference in temperature is nearly negligible during the whole procedure.
Discussion
The aim of this study was to evaluate the impact of the timevarying 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 [38]. 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 [39]. 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 timevarying h and constant h. In the first few seconds and when the electrode is close to the vessel, the timevarying 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 timevarying 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 timevarying h might be required but needs to be further investigated.
Conclusion
Previous studies considered a constant heat transfer coefficient throughout the ablation procedure. In this work, simulations were performed using a more realistic, timevarying 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 timevarying coefficient produces noticeable different results. However, this has no clinical impact for RF liver ablation procedure, which typically takes over 5 min.
References
 1.
Neeman Z, Wood BJ: Radiofrequency ablation beyond the liver. Tech Vasc Interv Radiol 2002, 5: 156–163. 10.1053/tvir.2002.36419
 2.
Dewhirst MW, Viglianti BL, LoraMichiels 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/0265673031000119006
 3.
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.
 4.
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.2322030184
 5.
Dewhirst MW, Vujaskovic Z, Jones E, Thrall D: Resetting the biologic rationale for thermal therapy. Int J Hyperthermia 2005, 21: 779–790. 10.1080/02656730500271668
 6.
Tungjiktkusolmun S, Staelin S, Haemmerich D, Tsai JZ, Cao H, Webster JG, LFT Jr, Mahvi DM, Voerperian V: ThreeDimensional FiniteElement Analyses for RadioFrequency Hepatic Tumor Ablation. IEEE Trans Biomed Eng 2002, 49: 3–9. 10.1109/10.972834
 7.
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.
 8.
Pennes HH: Analysis of tissue and arterial blood temperature in the resting human forearm. J Appl Physiol 1948, 1: 93–122.
 9.
Chen MM, Holmes KR: Microvascular contributions in tissue heat transfer. Ann New York Acad Sci 1980, 355: 137–150. 10.1111/j.17496632.1980.tb50742.x
 10.
Osman MM, Afify EM: Thermal modeling of the normal women's breast. J Biomech Eng 1984, 106: 123–130.
 11.
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.
 12.
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.2796092
 13.
Xuan Y, Roetzel W: Bioheat equation of the human thermal system. Chem Eng Technol 1997, 20: 268–276. 10.1002/ceat.270200407
 14.
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.2798007
 15.
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/02656730110081794
 16.
Deng ZS, Liu J: Blood perfusionbased model for characterizing the temperature fluctuation in living tissues. Physica A 2001, 300: 521–530. 10.1016/S03784371(01)003739
 17.
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/S07351933(01)00330X
 18.
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.1431263
 19.
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/00319155/48/11/307
 20.
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/00319155/50/15/010
 21.
Martini FH, Timmons MJ, McKinley MP: Human Anatomy. Upper Saddle River, New Jersey 07458: PrenticeHall Inc; 2000.
 22.
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.2895728
 23.
Stanczyk M, Leeuwen GMJV, V SAA: Discrete vessel heat transfer in perfused tissuemodel comparison. Phys Med Biol 2007, 52: 2379–2391. 10.1088/00319155/52/9/004
 24.
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/BF02348437
 25.
Gates G, Dore E: Streamline flow in the human portal vein. J Nucl Med 1973, 14: 79–83.
 26.
Garcier JM, Bousquet J, Alexandre M, Filaire M, Viallet JF, Vanneuville G, Boyer L: Visualisation of the portal flows by portoscanner. Surg Radiol Anat 2000, 22: 239–242. 10.1007/s0027600002394
 27.
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/00319155/48/24/010
 28.
Berjano EJ: Theoretical modeling for radiofrequency ablation: stateoftheart and challenges for the future. Biomed Eng Online 2006, 5: 24. 10.1186/1475925X524
 29.
Ebbini ES, Umemura SI, Ibbini M, Cain C: A cylindrical section ultrasound phasedarray applicator for hyperthermia cancer therapy. IEEE Trans. Ultrasonics, Ferroelectrics, and Frequency Control 1988, 35: 561–572. 10.1109/58.8034
 30.
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.
 31.
Valvano JW, Cochran JR, Diller KR: Thermal conductivity and diffusivity of biomaterials measured with selfheating thermistors. Int J Thermophys 1985, 6: 301–311. 10.1007/BF00522151
 32.
Panescu D, Whayne J, Fleschman SD, Mirotznik MS, Swanson DK, Webster JG: Threedimensional finiteelement analysis of current density and temperature distributions during radiofrequency ablation. IEEE Trans Biomed Eng 1995, 42: 879–890. 10.1109/10.412649
 33.
Tungjitkusolmun S, Cao H, Tsai J, Webster JG: Using Ansys for threedimensional electricalthermal models for radiofrequency catheter ablation. Proc 19th Ann Int Conf IEEE 1997, 1: 161–164.
 34.
Duck FA: Physical properties of tissue. London: Academic Press; 1990.
 35.
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.
 36.
Haemmerich D, Webster JG, Limanond P, Aziz D, Economou J, Busuttil R, Sayre J: Automatic control of finite element models for temperaturecontrolled radiofrequency ablation. Biomed Eng Online 2005, 4: 1–8. 10.1186/1475925X442
 37.
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/02656730601094415
 38.
Gray S: Gray's Anatomy. New York: Vintage Books; 1990.
 39.
Incropera FP, DeWitt DP: Fundamentals of Heat and Mass Transfer. New York, USA: John Willey & Sons, Inc; 1996.
Acknowledgements
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.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
IS and AFR conceived the study; IS and CSP carried out the computer simulation; AFR and DH helped with the interpretation of the results and with the draft of the manuscript. All authors have read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
dos Santos, I., Haemmerich, D., Pinheiro, C.d.S. et al. Effect of variable heat transfer coefficient on tissue temperature next to a large vessel during radiofrequency tumor ablation. BioMed Eng OnLine 7, 21 (2008) doi:10.1186/1475925X721
Received
Accepted
Published
DOI
Keywords
 Thermal Boundary Layer
 Ablation Procedure
 Convective Heat Transfer Coefficient
 Tumor Ablation
 Heated Region