- Open Access
In vitro calibration of a system for measurement of in vivo convective heat transfer coefficient in animals
BioMedical Engineering OnLine volume 5, Article number: 57 (2006)
We need a sensor to measure the convective heat transfer coefficient during ablation of the heart or liver.
We built a minimally invasive instrument to measure the in vivo convective heat transfer coefficient, h in animals, using a Wheatstone-bridge circuit, similar to a hot-wire anemometer circuit. One arm is connected to a steerable catheter sensor whose tip is a 1.9 mm × 3.2 mm thin film resistive temperature detector (RTD) sensor. We used a circulation system to simulate different flow rates at 39°C for in vitro experiments using distilled water, tap water and saline. We heated the sensor approximately 5°C above the fluid temperature. We measured the power consumed by the sensor and the resistance of the sensor during the experiments and analyzed these data to determine the value of the convective heat transfer coefficient at various flow rates.
From 0 to 5 L/min, experimental values of h in W/(m2·K) were for distilled water 5100 to 13000, for tap water 5500 to 12300, and for saline 5400 to 13600. Theoretical values were 1900 to 10700.
We believe this system is the smallest, most accurate method of minimally invasive measurement of in vivo h in animals and provides the least disturbance of flow.
Radiofrequency (RF) catheter ablation has been used to treat many types of cardiac arrhythmias such as atrial tachyarrhythmias (atrial tachycardia, atrial flutter, atrial fibrillation), atrioventricular nodal re-entrant tachycardia, Wolf-Parkinson-White syndrome, symptomatic supraventricular and ventricular tachycardia (fascicular VT, bundle branch re-entrant VT, idiopathic VT, ischemic VT) [1–5] with high success rate because of its controllability, high efficacy, low complications and minimal invasiveness. The technique of radiofrequency catheter ablation is to deliver high-frequency alternating electric current from 350 kHz to 1 MHz through the electrode catheter to generate a thermal lesion in myocardial tissue . The tissue in direct contact with the catheter is heated by resistive (ohmic) heating. Resistive heating in tissue is proportional to the power density . The thermal energy from the directly contacted tissue is then transferred to its vicinity by means of conduction and forms an RF lesion. Flowing blood near the tissue transfers heat away by convection and is a major cause of heat loss from the RF lesion. Convective heat transferred by the epicardial coronary artery is considered a minor heat loss in RF catheter ablation [6, 7].
Hot wire/film anemometer systems have been used to measure the blood velocity in vivo. Nerem et al.  used a hot film probe to study the velocity distribution in the aortas of dogs. Paulsen and Nissen  developed a safety system for a hot-film anemometer for blood-velocity measurement in humans. Yamaguchi et al.  performed turbulence measurements in the center of the canine ascending aorta using a hot-film anemometer. Paulsen et al.  analyzed the dynamic properties of a hot-film anemometer system for blood velocity measurements in humans.
In order to improve the electrodes and success rate of RF catheter ablation, researchers have used finite element method (FEM) modelling to simulate the ablation [7, 12, 13]. The value of the endocardial convective heat transfer coefficient (h) is essential for the simulation of heat loss from the endocardium to the blood pool in the cardiac chambers [7, 13]. Researchers have used values of h ranging from 44 to 6090 W/(m2·K) for various locations in the cardiac chambers [7, 13–15], however, none of those values came from in vivo measurement in animals. Because we did not find in vivo measurements of h in animals, we built our instrument, which is a Wheatstone-bridge circuit connected to a thin film resistive temperature detector (RTD) sensor, similar to a hot-wire anemometer circuit, and used it to measure the endocardial h first in vitro and then in vivo .
Absolute accuracy in measuring heat transfer coefficient h might be advantageous, but will be very difficult to attain given the real geometry of the sensor and the operating conditions envisioned for its use. Of equal value, and far more realizable, is a sensor that can provide repeatable, reasonably accurate estimates of local h, exhibiting the variations in time and space expected in the in vivo application. We will demonstrate the sensor described provides repeatable results in in vitro experiments, which are consistent with an estimated lower bound on h from heat transfer correlations on an idealized geometry.
In order to obtain the value of h inside the cardiac chambers, we performed in vitro experiments to ensure the capability of our measuring system, to test for any leakage current and to calibrate the system. The in vitro experiment setting was consisted of a circulation system (Fig. 1), the sensor (Figs. 2 and 3), the measuring circuit (Fig. 4) and the data acquisition program. The circulation system simulated flow rates of blood from 0 to 5 L/min to provide different flow rates for three different fluids (distilled water, tap water and saline). The catheter sensor, whose tip is a Pt thin film resistive temperature detector (RTD) sensor, was placed in the fluid flow, and formed one arm of the Wheatstone bridge circuit (the measuring circuit). The electric power consumed by the sensor and the temperature of the sensor were measured and saved using a data acquisition program and were then analyzed to yield h (using equation (16)).
A. Circulation system
Fig. 1 shows that the circulation system was composed of a fluid bath, a centrifugal pump, a flow meter and a fluid container. We measured flow with a rotameter Model 7200–0061 (King Instrument Co., Garden Grove, CA), which has a full scale of 2 GPM (8 L/min) and 3% accuracy. The fluid bath maintained the temperature of the fluid at 39°C (swine body temperature). From the fluid bath, the pump pumped the fluid to the flow meter. A valve at the flow meter was adjusted to vary the flow rate from 0 to 5 L/min. The fluid then flowed to the container through a 20 mm diameter, 500 mm long, solid PVC tube with the heated sensor on the axis in the tube.
Since the tube was more than 20 times longer than the diameter of the tube, this suggests that the expected laminar flow was fully developed at the measuring point. As the fluid flowed past the sensor, it dissipated heat from the sensor and cooled the sensor by means of convection. The higher the flow rate, the faster the heat was dissipated.
Because this system will be used to measure the endocardial convective heat transfer coefficient, h, in vivo, it is important that the measuring system (especially, the catheter sensor) will not create any leakage current while performing the h measurements. Since the electrical conductivity (σ) of blood is around 0.60 to 0.67 S/m at 20 to 25°C , we must perform the in vitro tests in fluid that has electrical conductivity similar to that of blood. Therefore, three different types of fluid were used: distilled water, tap water and saline. Distilled water (σ ≈ 10-4 to 10-3 S/m ), which has a very low electrical conductivity, was first used for the in vitro experiments (as control experiments) in order to avoid having any leakage current during the experiment. Tap water (σ ≈ 0.15 S/m ), which has an electrical conductivity higher than that of distilled water but lower than that of blood, and saline (σ ≈ 0.82 S/m ), which has electrical conductivity slightly higher than that of blood, were later used in the experiments in order to determine if any leakage current was created from the catheter sensor as well as to observe the differences of h obtained from these fluids with different electrical conductivity.
B. Measuring circuit
Fig. 4 shows our measuring system (similar to a hot wire/film anemometer system) using a resistive temperature detector (RTD) (Fig. 2), R s, forming one arm of a Wheatstone bridge. This circuit maintained the resistance of the sensor (R s) at a constant value, hence holding the temperature of the sensor constant (constant temperature mode). During the measurement, the sensor was heated approximately 5°C above the temperature of the flowing fluid (distilled water, tap water and saline). The electric power consumed by the sensor and the resistance of the sensor during the experiment were recorded and analyzed to yield the value of h.
The Wheatstone bridge had a resistive temperature detector (RTD) sensor in one arm of the bridge, and three wire-wound resistors (1 W), and a precision potentiometer in the other arms. A power operational amplifier supplied the power to the circuit and the sensor. The resistance of the sensor was maintained constant by the bridge at R S = R R (G - 2)/(G + 2) where R R is the total resistance of the upper right arm of the bridge, and G is the gain of the op amp . The sensor was electrically heated about 5°C above the fluid temperature because if using this system in vivo, higher temperatures may cause blood coagulation on the surface of the sensor and may damage the surrounding cells . Therefore we restricted the temperature difference to 5°C also for the in vitro tests. Dissipation of the heat occurred because the flowing fluid carried the warmed fluid away from the surface of the heated sensor. Increased velocity of the fluid reduced the fluid temperature next to the sensor and therefore the ohmic resistance of the sensor. Restoration of the sensor to its original working temperature was achieved by feedback controlled by the op amp. The sensor current increased to increase the power and the temperature of the sensor.
The voltage across the sensor (V 1) and the voltage across R 1 and R S (V A) (see Fig. 4) were recorded by the data acquisition program through the analog-to-digital converter (ADC) (12 bit, 100 kS/s, 8 analog inputs). V A and V 1 were used to determine the current that flowed through the sensor. V 1 was also used to determine the resistance and the temperature of the sensor. Once we knew the value of the current, the resistance and the temperature of the sensor, we calculated h using Newton's law of cooling :
q = hA(T s - T ∞) = Q h (1)
where: q is the heat flow (W), Q h is the electric power consumed by heating the sensor (W), h is the convective heat transfer coefficient (W/(m2·K)), A is the sensor area (m2), T S is the temperature of the heated sensor (K), and T ∞ is the temperature of the bulk fluid (K). However, since the sensor we used was not a bare Pt thin film, we calculated a correction.
C. Catheter sensor design
To measure in vivo h, we required a catheter that could be externally steered and placed against cardiac chamber and vessel walls. The sensor (model TFD, Omega Company) was at the tip of a cardiac ablation catheter and Loctite sealed the electric connection, covered the backside of the probe and rounded the sharp edges of the sensor. Fig. 2 shows the structure of the thin film Pt sensor. Its size was 1.9 × 9.4mm with a sensing area of 1.9 × 3.2mm (after covering the less-temperature sensitive area and the rough edges with Loctite). Fig. 3 shows the structure of the catheter with the sensor at the tip.
D. Temperature vs. resistance for the sensor
When electrically connecting the bare sensor to a catheter, the overall resistance of the catheter sensor was greater than that of the bare sensor because of the added resistance of the lead wire. We used an adjustable temperature water bath and a digital multimeter (model# HP34401A) to plot the resistance versus the temperature shown in Fig. 5. The average resistance difference between the bare sensor and the sensor with the lead wire was about 8.51 Ω. The resistance differences of the lead wire in air (25°C, room temperature) and in heated water (up to 45°C) were less than 0.1 Ω. Using this information, we used the DIN EN 60751 , resistance vs. temperature table. The resistance of the catheter sensor was R c = R s + 8.51 Ω, where R s was the resistance of the bare sensor.
Theoretical calculation of forced convection at the center of a 20 mm diameter tube
It was important to consider the closest theoretical prediction possible in order to guide our experimental testing. Therefore, we calculated the value of forced h for laminar flow, constant heat flux , using:
where: Nu is the Nusselt number (dimensionless), k is thermal conductivity (W/(m·K)), L is the length of the heated area (parallel to the direction of the flow), (For our catheter sensor, the length of the exposed, heated area was 3.2 × 10-3 m).
The average Nusselt number (Nu) for laminar flow, constant heat flux from a flat plate heated from x = 0 to x = L, and Re < 5 × 105 can be estimated from the following equation :
= 0.664(Re1/2)(Pr1/3) (3)
where: Re is the Reynolds number (dimensionless), Pr is the Prandtl number (dimensionless), α is the Womersley number (dimensionless), St is Strouhal number (dimensionless), ρ is the density of fluid (kg/m3), μ is dynamic viscosity (N·s/m2), ν is the kinematic viscosity = μ/ρ (m2/s), C P is specific heat at constant pressure (J/(kg·K)), u is the velocity of the flow at the sensor location (m/s)
F = Flow, D is the diameter of the tube (m), f is a characteristic frequency of the waveform (in this case, the heart rate) (Hz), A t is the cross-sectional area of the 20 mm diameter tube = π(0.01)2 = 3.1416 × 10-4m2.
For the calculation of , Re and h, the values of Pr, ρ, μ, C P, and k at 37°C and 39°C are listed in Table 1. For the in vitro experiments, the unsteadiness of the flow was insignificant (with very small value of the Womersley number) and could be neglected. However, in physiological situations, a larger value of the Womersley number would be obtained, thus equations (4) and (5) must be adjusted according to the unsteadiness of the flow due to the heart rate.
Theoretical calculation of free convection for a vertical plane
When placing a heated sensor in a still liquid, heat generated by the sensor is dissipated by free or natural heat convection. The average value of free or natural convection, caused by the fluid movement resulting from the change of the fluid density due to the heating process, can be calculated using :
where the subscript f indicates that the properties in the dimensionless groups are evaluated at the film temperature (T f), where:
where: T ∞ is the fluid temperature (°C), T w is the wall temperature (of the heated plane) (°C), and for constant heat flux surface:
where: Grf is the Grashof number at the film temperature, Prf is the Prandtl number at the film temperature, g is acceleration of gravity (m/s2), β is the temperature coefficient of thermal conductivity (1/K), x is the length of the exposed area of the heated sensor (= 3.2 × 10- 3 m), ΔT is temperature difference of the fluid and the wall (T ∞ - T w), C is a constant, and can be evaluated by the value of GrfPrf, m is a constant, and can also be evaluated by the value of GrfPrf, h is the local free convective heat transfer coefficient (W/(m2·K)).
For free convection from isothermal vertical planes, the values of local h and Nuf can be analyzed according to Table 2. Since the value of GrxPrf is lower than 104, we used Fig. 7-7 of  (free convection heat transfer from vertical isothermal plate, Nu vs. GrPr) to determine the value of without using the value of C and m.
The average convective heat transfer coefficient () for laminar region can be evaluated by :
In a water bath with constant temperature of 39°C, a sensor was placed vertically. The sensor was heated to a constant 5°C above the water temperature. The free h was h local = 1500 W/(m2·K) and = 1870 W/(m2·K).
Table 3 shows the values of u, Re, Nu and . Note that at zero flow rate, is calculated using free convection equations.
We calculated a correction for the sensor to obtain a more accurate value of h since a glass layer covers the top of the Pt thin film. A thermal circuit analogous to an electric circuit was used for this correction. We treated the heat transfer-rate (q) as a flow. We calculated the thermal resistances from the thermal conductivity, convective heat transfer coefficient, and thickness of the material and the area of the material. The temperature difference is analogous to the potential difference. The Fourier equation  may be written as:
When the sensor is heated, the heat from the Pt thin film conducts through this thin layer of glass, which is on the top of the Pt film. The equivalent thermal resistance of the glass layer is R G. The equivalent thermal resistance of the ceramic layer is R C and the equivalent thermal resistance of the Loctite layer is R L. Convection occurs on the surface of this layer, and is included in this thermal circuit as R conv. Fig. 6 shows the one-dimensional structure of the catheter sensor and the equivalent circuit of this system. The value of R G, R C, R L and R conv can be calculated using:
where: d G is the thickness of the glass layer = 0.027 mm, k G is the thermal conductivity of the glass = 1.38 W/(m·K), d C is the thickness of the ceramic substrate = 0.45 mm, k C is the thermal conductivity of the ceramic = 6.06 W/(m·K), d L is the thickness of the Loctite = 1.00 mm, k L is the thermal conductivity of the Loctite = 0.55 W/(m·K), A is the surface area of the glass layer = the surface area of the ceramic substrate = 3.2 mm × 1.9 mm = 6.08 × 10-6 m2
Heat from the layer of thin film Pt also conducts through the ceramic substrate and Loctite layer as shown in Fig. 6. The combination of the thermal resistance of the Loctite layer (R L) and the thermal resistance of the ceramic substrate (R C) is parallel to the thermal resistance of the glass layer. Thus for the catheter sensor, we calculate h, using:
where A is the exposed surface area of the catheter sensor = 1.9 mm × 3.2 mm
In vitro experimental results
Fig. 7 shows the value of h from 16 experiments (for each flow rate) at various flow rates in three different types of media: distilled water, tap water and saline. Table 4 shows the average values of h and the standard deviations of the values of h from each flow rate. All of the in vitro experiments were performed at 39°C similar to the swine body temperature.
The experimental results show that from 0 to 5 L/min, values of h in W/(m2·K) were for distilled water 5100 to 13000, for tap water 5500 to 12300, and for saline 5400 to 13600 (as shown in Fig. 7). At low flow rates (from 1 to 3 L/min), distilled water yielded the highest value of h (8700 W/(m2·K) at 1 L/min, which was about 13% higher than the lowest value of h obtained from tap water). At the higher flow rates (from 4 L/min or higher), saline yielded the highest value of h (13600 W/(m2·K) at 5 L/min, which was about 9% higher than the lowest value of h obtained from tap water). No significant differences of h among these three media were found at any flow rate.
It is important to consider the closest theoretical prediction possible in order to guide our experimental testing. Therefore, we compared the experimental results with the theoretical values and found that at all flow rates, the h results from the in vitro experiments of those three types of media yielded a higher value of approximately 1500 W/(m2·K) above the theoretical value.
We have built an instrument for the measurement of h and tested it in distilled water, tap water and saline. The measured h varied significantly with flow rate, as theory suggests it should, but did not follow the square root dependence on Reynolds number that was expected. The present data varied approximately with the 0.3 power, whereas theory suggests 0.5.
Because of the small size of the sensor, the Grashof or Rayleigh Number and the Reynolds number (for the forced convection flows) were small. Most of the correlations were more reliable at the upper end of the Grashof or Reynolds number ranges, hence the experimental data agreed better with the theoretical calculation at higher flows.
At the lower flow rate in stagnant liquid (at low Grashof Number), the experimental flow regime is uncertain. The natural convection correlation assumes a flat vertical surface, no edge effects, constant fluid properties, and no forced convection (fluid movement) at all, which are difficult to replicate. Most other experiments in natural convection attempt to achieve relatively large Grashof numbers, where the flow is strongly driven by free convection with negligible small amounts of forced convection. However, in our case, Gr was small, so normally negligible small forced convection had a big effect and could not be neglected.
The data changed slightly with the type of medium. All three data sets were within ± 7% of the mean line (excluding the free convection results, which were not expected to lie on the same lines as forced convection). No leakage current was detected in all in vitro experiments, which suggested that this system should be safe to perform the measurement in vivo.
We did not measure h at the tube wall in vitro because the tube wall was rigid and the sensor would have been forced away from the wall into flow streamlines different from those at the tube wall. This is the reason we measured h on the axis, with known flow. Once we calibrated h on the axis, we later used the sensor in vivo against the soft tissue of the endocardium or vessels to measure h. The sensor was small enough so that it indented the soft tissue to place the sensor face close to the wall contours and achieve minimal disruption of the flow streamlines at the endocardium or vessel wall. This enabled measurement of the value of h at the endocardium or vessel wall, which was less than the value of h away from the wall.
We believe this system is the smallest, most accurate method of minimally invasive measurement of in vivo h in animals and provides the least disturbance of flow.
Chen SA, Chiang CE, Tai CT, Lee H, Chang MS: Future ablation concepts of tachyarrhythmias. J Cardiovasc Electrophysiol 1995, 6: 852–862.
Haissaguerre M, Jais PD, Shah DC, Gencel L, Pradeau V, Garrigues S, Chouairi S, Hocini M, Le Metayer P, Roudaut R, Clementy J: Right and left atrial radiofrequency catheter therapy of paroxysmal atrial fibrillation. J Cardiovasc Electrophysiol 1996, 7: 1132–1144.
Lau CP, Tai YT, Lee PWH: The effects of radiofrequency ablation versus medical therapy on the quality-of-life and exercise capacity in patients with accessory pathway-medical supraventricular tachycardia: a treatment comparison study. PACE 1995, 18: 424–432.
Mirotznik MS, Demazumder DT, Jones JR, Schwartzman DS: Heating distribution of multipolar radiofrequency ablation catheters. Proc 19th Int Conf IEEE Eng Med Biol Soc 1997, 1: 157–160.
Stuart TP, Nicholson IA, Nunn GR, Rees A, Trieu L, Daly MPJ, Wallace EM, Ross DL: Effect of atrial radiofrequency ablation designed to cure atrial fibrillation on atrial mechanical function. J Cardiovasc Electrophysiol 2000, 11: 77–82.
Haines DE: The biophysics of radiofrequency catheter ablation in the heart: the importance of temperature monitoring. PACE 1993, 16: 586–591.
Tungjitkusolmun S, Vorperian VR, Bhavaraju N, Cao HJ, Tsai JZ, Webster JG: Guidelines for predicting lesion size at common endocardial locations during radio-frequency ablation. IEEE Trans Biomed Eng 2001, 48: 194–201.
Nerem RM, Seed WA, Wood NB: An experimental study of the velocity distribution and transition to turbulence in the aorta. J Fluid Mech 1972, 52: 137–160.
Paulsen PK, Nissen T: Patient safety unit for a hot-film anemometer, used for blood-velocity determination in humans. Med Biol Eng Comput 1982, 20: 625–627.
Yamaguchi T, Kikkawa S, Yoshikawa T, Tanishita K, Sugawara M: Measurement of turbulence intensity in the center of the canine ascending aorta with a hot-film anemometer. J Biomed Eng 1983, 105: 177–187.
Paulsen PK, Hasenkam JM, Nygaard H, Gormsen J: Analysis of the dynamic properties of a hot-film anemometer system for blood velocity measurements in humans. Med Biol Eng Comput 1987, 25: 195–200.
Cao H, Speidel MA, Tsai JZ, Van Lysel MS, Vorperian VR, Webster JG: FEM analysis of predicting electrode-myocardium contact from RF cardiac catheter ablation system impedance. IEEE Trans Biomed Eng 2002, 49: 520–526.
Jain MK, Wolf PD: Finite element analysis predicts dose-response relationship for constant power and temperature controlled radiofrequency ablation. Proc 19th Int Conf-IEEE Eng Med Biol Soc 1997, 1: 165–168.
Bhavaraju NC: Heat transfer modeling during cardiac ablation in swine myocardium. PhD thesis. University of Texas at Austin, Department of Biomedical Engineering; 1999.
Tungjitkusolmun S, Woo EJ, Cao H, Tsai JZ, Vorperian VR, Webster JG: Thermal-electric finite element modeling for radio frequency cardiac ablation: effects of changes in myocardial properties. Med Biol Eng Comput 2000, 38: 562–568.
Tangwongsan C, Will JA, Webster JG, Meredith KL Jr, Mahvi DM: In vivo measurement of swine endocardial convective heat transfer coefficient. IEEE Trans Biomed Eng 2004, 51: 1478–1486.
Bronzino JD: The Biomedical Engineering Handbook. Volume I. 1st edition. New York, CRC Press; 1995.
Hayt WH Jr: Engineering electromagnetics. 5th edition. New York, McGraw-Hill international edition; 1989.
Lide DR: Handbook of Chemistry and Physics. 81st edition. New York, CRC Press; 2000–2001.
Lenntech TDS and electrical conductivity [http://www.lenntech.com/tds-ec_engels.htm]
Freire RCS, Deep GS, Oliveira A: Effect of operational amplifier parameters on the performance of feedback structures with thermoresistive sensors. IEEE Instrum Meas Technol Conf 1997, 898–903.
Holman JP: Heat Transfer. 6th edition. New York, McGraw-Hill; 1986.
RTD table for DIN EN60751 for class A and class B In Handbook and Encyclopedia: Temperature. MM, omega.com; 2000:z-251-z-255.
This work was supported by the National Institute of Health (NIH) under grant HL56413. We thank Prof. Robert J. Moffat, emeritus professor from Stanford University, Department of Mechanical Engineering, for his valuable suggestions and comments. We thank Dr. James A. Will from Department of Animal Health and Biomedical Science, University of Wisconsin-Madison for helping us in many ways. We thank Dr. Dorin Panescu of the Boston Scientific Corporation, who supplied the probe.
The author(s) declare that they have no competing interests.
CT carried out the design of the test equipment, ran tests, and drafted the manuscript. LC ran tests. JGW conceived of the study, and participated in its design and coordination and helped to draft the manuscript. PVF contributed to the heat transfer sections. All authors read and approved the final manuscript.
Chanchana Tangwongsan, Louay Chachati, John G Webster and Patrick V Farrell contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Tangwongsan, C., Chachati, L., Webster, J.G. et al. In vitro calibration of a system for measurement of in vivo convective heat transfer coefficient in animals. BioMed Eng OnLine 5, 57 (2006). https://doi.org/10.1186/1475-925X-5-57
- Free Convection
- Forced Convection
- Convective Heat Transfer Coefficient
- Constant Heat Flux
- Wheatstone Bridge