 # The Open Biomedical Engineering Journal

ISSN: 1874-1207 ― Volume 13, 2019

# Mathematical Models Based on Transfer Functions to Estimate Tissue Temperature During RF Cardiac Ablation in Real Time

Jose Alba-Martínez1, Macarena Trujillo2, Ramon Blasco-Gimenez3, Enrique Berjano1, *
1 Biomedical Synergy, Electronic Engineering Department, Universitat Politècnica de Valencia, Spain
2 Instituto Universitario de Matemática Pura y Aplicada, Universitat Politècnica de Valencia, Spain
3 Dpto. Ingeniería de Sistemas y Automática, Universitat Politècnica de Valencia, Spain

## Abstract

Radiofrequency cardiac ablation (RFCA) has been used to treat certain types of cardiac arrhythmias by producing a thermal lesion. Even though a tissue temperature higher than 50ºC is required to destroy the target, thermal mapping is not currently used during RFCA. Our aim was thus to develop mathematical models capable of estimating tissue temperature from tissue characteristics acquired or estimated at the beginning of the procedure (electrical conductivity, thermal conductivity, specific heat and density) and the applied voltage at any time. Biological tissue was considered as a system with an input (applied voltage) and output (tissue temperature), and so the mathematical models were based on transfer functions relating these variables. We used theoretical models based on finite element method to verify the mathematical models. Firstly, we solved finite element models to identify the transfer functions between the temperature at a depth of 4 mm and a constant applied voltage using a 7Fr and 4 mm electrode. The results showed that the relationships can be expressed as first-order transfer functions. Changes in electrical conductivity only affected the static gain of the system, while specific heat variations produced a change in the dynamic system response. In contrast, variations in thermal conductivity modified both the static gain and the dynamic system response. Finally, to assess the performance of the transfer functions obtained, we conducted a new set of computer simulations using a controlled temperature protocol and considering the temperature dependence of the thermal and electrical conductivities, i.e. conditions closer to those found in clinical use. The results showed that the difference between the values estimated from transfer functions and the temperatures obtained from finite element models was less than 4ºC, which suggests that the proposed method could be used to estimate tissue temperature in real time.

Keywords: Cardiac ablation, closed loop control, finite element method, radiofrequency ablation, temperature controlled ablation, theoretical model..

### Article Information

#### Identifiers and Pagination:

Year: 2012
Volume: 6
First Page: 16
Last Page: 22
Publisher Id: TOBEJ-6-16
DOI: 10.2174/1874120701206010016

#### Article History:

Acceptance Date: 3/2/2012
Electronic publication date: 8/3/2012
Collection year: 2012

#### Article Metrics:

0

##### Total Statistics:

Full-Text HTML Views: 730
Abstract HTML Views: 561

##### Unique Statistics:

Full-Text HTML Views: 461
Abstract HTML Views: 394
Geographical View

open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted, non-commercial use, distribution and repr oduction in any medium, provided the work is properly cited.

* Address correspondence to this autor at the Biomedical Synergy, Electronic Engineering Department, Universitat Politècnica de Valencia, Spain; Tel: 34-963877607; Fax: 34-963877609; Email: eberjano@eln.upv.es

## INTRODUCTION

Radiofrequency cardiac ablation (RFCA) is currently used to treat some types of cardiac arrhythmia. This technique uses radiofrequency (RF) current (≈500 kHz) to produce a thermal lesion and hence tissue necrosis in the target zone responsible for the arrhythmia. Electrical current is delivered to the tissue through a small active electrode placed at the tip of a percutaneous catheter and a large dispersive electrode located on the patient’s back. There are two modes of delivering radiofrequency energy: constant voltage and temperature-controlled. The constant voltage mode means that the applied voltage is constant; under this condition high temperatures (approximately 100ºC) can be reached in the tissue around the electrode tip in few seconds and lesion size is limited by the charring of tissue around the electrode. The temperature-controlled mode, which can make a larger lesion than the constant voltage approach [1OG Anfinsen, H Aass, E Kongsgaard, A Foerster, H Scott, and JP Amlie, "“Temperature-controlled radiofrequency catheter ablation with a 10-mm tip electrode creates larger lesions without charring in the porcine heart”", J. Interv. Card. Electrophysiol Dec, vol. 3 , no. 4 , pp. 343-351, 1999.], is the one most frequently used in clinical practice [1OG Anfinsen, H Aass, E Kongsgaard, A Foerster, H Scott, and JP Amlie, "“Temperature-controlled radiofrequency catheter ablation with a 10-mm tip electrode creates larger lesions without charring in the porcine heart”", J. Interv. Card. Electrophysiol Dec, vol. 3 , no. 4 , pp. 343-351, 1999., 2HH Petersen, X Chen, A Pietersen, JH Svendsen, and S Haunso, "“ Tissue temperatures and lesion size during irrigated tip catheter radiofrequency ablation: An in vitro comparison of temperaturecontrolled irrigated tip ablation, power-controlled irrigated tip ablation, and standard temperature-controlled ablation”", Pacing Clin. Electrophysiol, vol. 23, no. 1 , pp. 8-17, 2000.] and consists of modulating the applied voltage to keep the temperature constant at the active electrode. The temperature sensor embedded in the electrode tip measures the approximate temperature at the tissue-electrode interface. Even though a tissue temperature higher than 50ºC is required to destroy the target, thermal mapping is not currently used during RFCA. This is an important issue, since tissue temperatures above 100ºC are associated with charring and can cause a thrombus [3F Gaita, D Caponi, M Pianelli, M Scaglione, E Toso, F Cesarani, C Boffano, G Gandini, MC Valentini, R De Ponti, F Halimi, and JF Leclercq, "“Radiofrequency catheter ablation of atrial fibrillation: a cause of silent thromboembolism? Magnetic resonance imaging assessment of cerebral thromboembolism in patients undergoing ablation of atrial fibrillation”", Circulation, vol. 26 122(), no. 17 , pp. 1667-1673, 2010.]. Safer ablations could thus be achieved by knowing the temperature in the tissue in order to ensure the creation of the thermal lesion and to avoid overheating. Our aim was thus to develop mathematical models to estimate tissue temperature in real time from the tissue characteristics (electrical conductivity (σ), thermal conductivity (k), specific heat (c) and density (ρ)) measured or estimated at the beginning of the procedure and the applied voltage at any time.

Previous work has been conducted on estimating lesion progress on-line by studying impedance progress [4WM Hartung, ME Burton, AG Deam, PF Walter, K McTeague, and JJ Langberg, "“Estimation of temperature during radiofrequency catheter ablation using impedance measurements”", Pacing Clin Electrophysiol Nov , vol. 18 , no. 11 , pp. 2017-2021, 1995.-6WC Ko, SK Huang, JL Lin, WY Shau, LP Lai, and PH Chen, "“ New method for predicting efficiency of heating by measuring bioimpedance during radiofrequency catheter ablation in humans”", J. Cardiovasc Electrophysiol Jul, vol. 12 , no. 7 , pp. 819-823, 2001.]. Numerical modeling (e.g. based on the Finite Element Method (FEM)) has also been proposed to estimate tissue temperature during RFCA [7YC Lai, YB Choy, D Haemmerich, VR Vorperian, and JG Webster, "“Lesion size estimator of cardiac radiofrequency ablation at different common locations with different tip temperatures”", IEEE Trans. Biomed. Eng, vol. 51, pp. 1859-1864, 2004.]. This method calculates the temperature distribution by means of solving numerically partial differential equations problems [8EJ Berja, "“Theoretical modeling for radiofrequency ablation: state-of-the-art and challenges for the future”", Biomed. Eng. Online Apr, vol. 18 , no. 5 , p. 24, 2006.] based on physical laws and hence the models using this method can be considered as structural models. We previously proposed and experimentally verified autoregressive moving average models with exogenous input (ARMAX) aimed at estimating tissue temperature during ablation [9R Blasco-Gimenez, JL Lequerica, M Herrero, F Hornero, and EJ Berja, "“Black-box modeling to estimate tissue temperature during radiofrequency catheter cardiac ablation: Feasibility study on an agar phantom model”", Physiol. Meas Apr, vol. 31 , no. 4 , pp. 581-594, 2010.]. These were clearly non-structural models since they were not based on physical laws. Our present proposal is based on mathematical models to estimate tissue temperature at a depth of 4 mm (T4) by measuring the applied voltage (U=V2) and initial value of the tissue characteristics, (see Fig. 1 ). These models are neither directly based on physical laws (e.g. FEM based) nor are they totally non-structural (e.g. ARMAX models), but are based on transfer functions which consider the biological tissue as a system with inputs (applied voltage and initial value of tissue characteristics) and output (tissue temperature). They would be relatively simple to implement while retaining all the important characteristics of more complex models and hence could be used for real time estimation of tissue response or as a theoretical basis to obtain more sophisticated temperature controllers. In this study, we set up the models using theoretical models based on FEM, and then verified the feasibility of the proposal by other FEM models. Fig. (1)Diagram of the proposed mathematical model to estimate tissue temperature. The initial value of the tissue values characteristics, (electrical conductivity (σ), thermal conductivity (k), specific heat (c)) and the applied voltage (U=V2) are inputs of the mathematical model. The output is the estimated temperature at a depth of 4 mm (T4). Fig. (2)A: Two-dimensional theoretical model. The active electrode tip is inserted into cardiac tissue to depth P=1.25 mm. Dimensions R, Z and L are obtained by means of convergence tests (see text). T4 shows the temperature estimation point (4 mm deep). B: biological tissue is modeled as a dynamic system with an input (applied voltage) U=V2and output tissue temperature (at depth of 4 mm), which corresponds with temperature at T4 (see Fig. 1A ), thus obtaining a first order transfer function G(s). Fig. (3)Temperature evolution obtained at 4 mm insertion depth (T4) when tissue characteristics were changed A: electrical conductivity; B: thermal conductivity; and C: specific heat. Fig. (4)Effect of the variation invarying tissue characteristics (electrical and thermal conductivities and specific heat) on gain Kv (left) and constant time τ (right) of the transfer functions obtained to implement the mathematical models. Fig. (5)Temperature evolution from FEM based models (solid line) and temperature estimated from mathematical model (dashed line). Ten cases were considered with changes in the values of tissue characteristics (see Table 2). Fig. (6)Error evolution in ºC. The error was obtained from the difference between the estimated tissue temperature (computed from the mathematical model) and the real temperature (obtained by the theoretical model using COMSOL Multiphysics). Ten cases were considered with changing in the values of tissue characteristics (see Table 2).

## METHODS

### Description of the Theoretical Model

We considered an active 7Fr and 4 mm long platinum-iridium electrode as used in clinical procedures. This electrode was positioned perpendicular to the tissue, which implied a rotational symmetry axis and allowed a two-dimensional model to be considered. Fig. (2A ) shows the geometry and dimensions of the model, which includes the active electrode with a fragment of the polyurethane probe, a thermistor embedded in the active electrode tip, a section of the thermistor coating material and a fragment of cardiac tissue. We considered the electrode to be inserted into the cardiac tissue to a depth of 1.25 mm [10D Schutt, EJ Berjano, and D Haemmerich, "“Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: a computational modeling study”", Int. J. Hyperthermia Mar, vol. 25 , no. 2 , pp. 99-107, 2009.] and the study focused on the temperature evolution at a depth of 4 mm (T4) (see Fig. 2A ). From a mathematical point of view, the model was based on a coupled electric-thermal problem which was solved using the Bioheat Equation [8EJ Berja, "“Theoretical modeling for radiofrequency ablation: state-of-the-art and challenges for the future”", Biomed. Eng. Online Apr, vol. 18 , no. 5 , p. 24, 2006.]:

(1)
where T is temperature, t is time, q is heat source produced by the RF power, Qp is heat loss by blood perfusion and Qm is metabolic heat generation. The terms Qp and Qm are insignificant in RFCA and were thus not considered. The heat source q (Joule losses) is given by q=JE where J is the current density and E is the electric field strength. The electrical problem was solved using the Laplace equation ∇(∇V)=0 in which V is voltage. The electric field was calculated by means of E and J using Ohm's law (JE), where σ is electrical conductivity. We used a quasi-static approach due to the frequencies used in RF ablation (≈500 kHz) and for the geometric area of interest the tissues can be considered as purely resistive [11JD Doss, "“Calculation of electric fields in conductive media”", Med. Phys Jul-Aug, vol. 9 , no. 4 , pp. 566-573, 1982.]. The blood circulating inside the cardiac chamber was not considered in the model domain. However, the blood-tissue interface and blood-electrode interface were really model boundaries (see Fig. 2A ).

The electrical boundary conditions were: voltage zero density in the transversal direction to the symmetry axis, blood-tissue interface and blood-electrode interface. A boundary condition of voltage zero was set at the dispersive electrode. The voltage at the active electrode tip was 15 V. This value was chosen in order to keep Tmax below 100ºC. The thermal boundary conditions were: null thermal flux in the transversal direction to the symmetry axis and constant temperature at the dispersive electrode and at the outer end of the plastic probe. The circulating blood inside the cardiac chamber produces a cooling effect on the tissue and on the electrode surface, which was modeled by means of two thermal convection coefficients htissue and helec which represent the cooling effect at the blood-tissue interface and at the blood-electrode interface respectively (see Fig. 2A ). We considered a medium blood flow rate, with values for helec of 3636 W/m2K and for htissue of 708 W/m2K [10D Schutt, EJ Berjano, and D Haemmerich, "“Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: a computational modeling study”", Int. J. Hyperthermia Mar, vol. 25 , no. 2 , pp. 99-107, 2009.] and an initial temperature of 37ºC. Table 1 shows the values of the physical characteristics of the materials [10D Schutt, EJ Berjano, and D Haemmerich, "“Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: a computational modeling study”", Int. J. Hyperthermia Mar, vol. 25 , no. 2 , pp. 99-107, 2009., 12D Panescu, JG Whayne, SD Fleischman, MS Mirotznik, DK Swanson, and JG Webster, "“Three-dimensional finite element analysis of current density and temperature distributions during radiofrequency ablation”", IEEE Trans. Biomed. Eng Sep, vol. 42 , no. 9 , pp. 879-890, 1995.]. Electrical and thermal conductivity variations with temperature were not considered in these FEM models in order to obtain a relatively simple mathematical model which could be used as a first step towards a more complex (non-linear model). Table 1

Characteristics of the Materials Used in the Model [10D Schutt, EJ Berjano, and D Haemmerich, "“Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: a computational modeling study”", Int. J. Hyperthermia Mar, vol. 25 , no. 2 , pp. 99-107, 2009., 12D Schutt, EJ Berjano, and D Haemmerich, "“Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: a computational modeling study”", Int. J. Hyperthermia Mar, vol. 25 , no. 2 , pp. 99-107, 2009.]

We used COMSOL Multiphysics software (COMSOL Inc., Burlington, MA, USA) to implement the numerical solution based on the FEM. In order to avoid boundary effects, the model dimensions R, Z and L (see Fig. 2A ) were estimated by means of convergence tests in which the value of the maximum temperature reached in the tissue (Tmax) after 120 s of RF heating was used as a control parameter. Parameter values were increased by 1 mm in each simulation. If there was a difference of less than 0.5% between Tmax and the same parameter in the previous simulation, the former values were considered as appropriate. The mesh was initially heterogeneous and a finer mesh was used around the electrode-tissue interface, where the highest gradient was expected. The mesh was then refined to determine adequate spatial resolution and the time-step was reduced to determine adequate temporal resolution.

### Identification of the Transfer Function

In the context of our study, the biological tissue was considered to be a dynamic system with applied voltage (V) as input signal and the temperature obtained at a depth of 4 mm (T4) as an output signal (see Fig. 2B ). It is assumed that the dynamic evolution of the system output (T4) with respect to the system input (U) can be expressed mathematically in terms of a differential equation. If the differential equation is linear, the mathematical input/output relationship can be expressed as a transfer function.

The first step was to obtain an accurate identification of the system response of the modeled tissue, which involves obtaining the transfer function G(s). This was carried out firstly by applying a constant voltage of 15 V for 300 s. The The model is assumed to be U(t)=V2(t) in which V was the applied voltage. We checked that the active electrode tip temperature stabilized around 55°C. We then obtained the evolution of T4, and MATLAB was then used to estimate the model response.

A first-order transfer function was obtained from the simulated evolution. The estimated transfer functions G(s) in Laplace-transform quantities had the form:

(2)
where ΔT4 is the temperature increase (above a base temperature of 37º), U(s) is the square of the voltage applied to the electrode, KV is the steady state gain of the system expressed in ºC/V2, and τ is the time constant. A first order transfer function was used since it is the simplest transfer function that leads to results reasonably similar to those obtained from finite element simulations.

From a practical point of view and using the above parameters for the first model obtained, the ablation final temperature could be calculated as follows:

(3)
where ΔT4final =V2 KV, being V the applied voltage and Tinitial = 37°C.

### Description of the Analyzed Cases

In order to study the influence of the tissue characteristics on the mathematical behavior of the ablation process we carried out several computer simulations with varying values of σ, k and c. We considered a variation from −75% to 100%, in steps of 25%, beginning with the value shown in Table 1. We show an example:

(4)
where k=0.531 W/mK and k-25%= 0.398 W/mK is the thermal conductivity variation. A first-order transfer function (2) was identified for each case, with corresponding values of KV and τ. This group of transfer functions was obtained in order to study the relationship between the transfer function gain (KV) and time constant (τ) with the variations in tissue specific heat and electrical and thermal conductivity. From the series of the KV and τ values, a general mathematical relationship for the variation of KV and τ for different values of σ, k and c was obtained.

### Building the Mathematical Model

We divided the study into two parts: firstly, we assessed the influence of the values of σ, k and c independently, i.e. for each variable separately in order to study the relationship between gain and time constant in each tissue characteristic. From these relationships we were able to obtain a mathematical model.

Secondly, in order to obtain better-fitting parameters in the first-order transfer function of the mathematical model, we considered simultaneous variations around the ends of σ, k and c and then used the linear minimum square fit to obtain more accurate parameters.

### Assessment of the Performance of the Models Under Different Conditions

Having built the mathematical model based on a first-order transfer function with fixed applied voltage and temperature non dependent on the electrical and thermal conductivity, in order to assess the accuracy of the models we carried out new computer simulations using a temperature-controlled protocol (variable applied voltage to maintain target temperature in the temperature sensor) and considering the temperature to be dependent on electrical and thermal conductivity, mimicking conditions similar to those found in clinical use.

In the models, we used random values around the center of variations (see Table 2). In addition, for a better approximation to clinical reality we considered the simulations to have a controlled temperature protocol and also that temperature was dependent on tissue thermal and electrical conductivity, with the aim of keeping the temperature in the active electrode approximately constant. Therefore, we used a PI controller with a target temperature of 55ºC. The values of the Kp and Ki of the PI controller were 4.78 and 3.39, respectively [13J Alba, M Trujillo, R Blasco-Gimenez, and E Berja, "“Could be advantageous to tune the temperature controller during radiofrequency ablation? A feasibility study using theoretical models”", Int. J. Hyperthermia, vol. 27 , no. 6 , pp. 539-548, 2011.]. The change in cardiac tissue electrical conductivity was taken as +1.5%/ºC [10D Schutt, EJ Berjano, and D Haemmerich, "“Effect of electrode thermal conductivity in cardiac radiofrequency catheter ablation: a computational modeling study”", Int. J. Hyperthermia Mar, vol. 25 , no. 2 , pp. 99-107, 2009.] and in thermal conductivity as +0.001195K-1 [14NC Bhavaraju, H Cao, DY Yuan, JW Valvano, and G Webster John, "“Measurement of directional thermal properties of biomaterials”", IEEE Trans. Biomed. Eng, vol. 48 , no. 2 , pp. 61-267, 2001.]. The temperature evolution at a depth of 4 mm was considered as the real temperature, in order to be comparable to that estimated by the mathematical models. From the transfer functions temperature progress was estimated by means of:

(5)
where T4(t) is the inverse Laplace transform of G(s)V2(s), and is obtained by numerical integration. The error was calculated as the difference between the estimated temperature (obtained by the mathematical model) and the real temperature (obtained by COMSOL).

## RESULTS

### Study of Temperature Evolution at a Depth of 4 mm

Fig. (3A ) shows the temperature evolution at a depth of 4 mm (T4) when electrical conductivity was varied. No changes were observed in the time needed to reach Tfinal during these variations. The final temperature was seen to be directly proportional to the variation of electrical conductivity. Fig. (3B ) shows temperature evolution with variations in thermal conductivity. The value of T4 was found to be inversely proportional to the variations of thermal conductivity, i.e. high values for T4 were obtained with low thermal conductivity values and vice versa. Changes were observed in the time needed to reach Tfinal, high thermal conductivity values produced faster heating than low values. Finally, Fig. (3C ) shows temperature evolution for changes in specific heat. In this case, T4 was also inversely proportional to the variation of specific heat. Changes were observed in the time needed to reach Tfinal, low specific heat values produced faster heating than high values.

### Static and Dynamic Characterization

Fig. (4 ) shows the values of the gain (KV) and time constant (τ) against variations in tissue characteristics. Fig. (4A ) shows a linear variation of the gain with respect to electrical conductivity. Gain was practically constant when specific heat was varied, and was a hyperbolic function of thermal conductivity. Fig. (4B ) shows the evolution of τ against variations in tissue characteristics. Note that the value of the time constant was not affected by changes in electric conductivity. Conversely, the time constant was linearly dependent on specific heat variations and inversely proportional to thermal conductivity.

From this behavior we tried to obtain transfer functions for each tissue characteristic:

### Fit for Variations in Electrical Conductivity

In this case, the relation was constant with time constant and linear with the gain was linear. We therefore proposed the following relationship:

(6)
where the gain can be represented by Kv=Ass and the time constant τ=τ0 where As=0.07632 ºCm/s obtained by using a linear minimum square fit on the obtained data and t is at 52.7 s.

### Fit for Variations in Thermal Conductivity

Gain and the time constant depended on thermal conductivity as a hyperbolic + constant function. Thermal conductivity variations were therefore considered as:

(7)
where Ak=0.01900 ºCW/mK, Bk=0.005059ºC, Dk=25.2171 and Ek=5.423 were obtained by using a linear minimum square fit on the obtained data.

### Fit for Variations in Specific Heat

The time constant showed a linear dependence on thermal conductivity and gain was almost constant. Specific heat variations were therefore considered as:

(8)
where the gain can be represented by Kvc=0.041ºC and Dc=0.01478 and Ec=5.009 are obtained by using a linear minimum square fit.

### General Mathematical Model

We know that the time constants in thermal systems can be expressed as τ=c/k where τ is the time constant, c is the specific heat and k is thermal conductivity. By taking this into account, we combined equations (6), (7) and (8) to propose the following mathematical model:

(9)
where A, B, D and E could be obtained by using a linear minimum square fit on the obtained data. However, to obtain a better fit we considered simultaneous variations of σ, k and c. The resulting values were A=0.0342ºCW/sK, B=0.0114ºCm/s, D=0.0083 and E=3.8908.

### Assessment of the performance of the models under different conditions

Fig. (5 ) shows the temperature evolutions obtained from FEM and from the mathematical model, both of which correspond to tissue temperature at a depth of 4 mm. Ten different tissue characteristic values were considered (see Table 2). The mathematical models underestimated the temperature in all cases. Fig. (6 ) shows the error evolution in each case, which was obtained from the difference between the estimated tissue temperature and the real temperature (obtained from FEM-based models). It can be seen that the differences in temperature are less than −4ºC.

## DISCUSSION

With the aim of developing mathematical models to estimate in real time tissue temperature from the tissue characteristics measured or estimated at the beginning of the procedure (electrical conductivity, thermal conductivity, specific heat and density) and applied voltage at a given time, we built a mathematical model based on a first-order transfer function. The temperatures estimated from the mathematical models proved to be reasonably similar to those obtained from numerical models, since they always kept lower than 5ºC, which is around the range of lethal isotherm found in experiments using cardiac tissue [15M Wood, S Goldberg, M Lau, A Goel, D Alexander, F Han, and S Feinstein, "“Direct measurement of the lethal isotherm for radiofrequency ablation of myocardial tissue”", Circ. Arrhythm. Electrophysiol Jun, vol. 4 , no. 3 , pp. 373-378, 2011.]. The maximum error of 4ºC was obtained from the maximum difference between the tissue temperature estimated by the mathematical model and tissue temperature obtained from COMSOL Multiphysics. In all the cases the temperature was underestimated. The temperature dependence of the thermal and electric conductivity was not initially considered in the model and may be the likely reason for the consistent underestimation of tissue temperature.

In order to obtain the input parameters for the mathematical models the initial values of the tissue characteristics (σ, k and c) must be known and can be obtained from experimental measurements, as described in previous studies. For example, techniques based on automatic swept-frequency network and impedance analysers have been proposed to measure electrical conductivity [16S Gabriel, RW Lau, and C Gabriel, "The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz Phys”", Med. Biol, vol. 41, pp. 2251-2269, 1996.], and the empirical relation between the power delivered to a heated thermistor by a heating pulse and the temperature rise recorded by the sensing thermistor placed at a distance from the sample could be used to measure thermal conductivity [14NC Bhavaraju, H Cao, DY Yuan, JW Valvano, and G Webster John, "“Measurement of directional thermal properties of biomaterials”", IEEE Trans. Biomed. Eng, vol. 48 , no. 2 , pp. 61-267, 2001.]. Further studies will be undertaken in order to obtain a mathematical model that takes into account temperature-dependent parameters and new methods of estimating initial tissue characteristics.

Clearly, tissue temperature estimation is currently a challenge. Although previous studies estimated temperature by means of impedance measurements [4WM Hartung, ME Burton, AG Deam, PF Walter, K McTeague, and JJ Langberg, "“Estimation of temperature during radiofrequency catheter ablation using impedance measurements”", Pacing Clin Electrophysiol Nov , vol. 18 , no. 11 , pp. 2017-2021, 1995.-6WC Ko, SK Huang, JL Lin, WY Shau, LP Lai, and PH Chen, "“ New method for predicting efficiency of heating by measuring bioimpedance during radiofrequency catheter ablation in humans”", J. Cardiovasc Electrophysiol Jul, vol. 12 , no. 7 , pp. 819-823, 2001.] and numerical modeling [7YC Lai, YB Choy, D Haemmerich, VR Vorperian, and JG Webster, "“Lesion size estimator of cardiac radiofrequency ablation at different common locations with different tip temperatures”", IEEE Trans. Biomed. Eng, vol. 51, pp. 1859-1864, 2004.], here we present a new method based on simple mathematical models. Our findings suggest that if the applied voltage (U=V2) and the initial value of the tissue characteristics (σ, k and c) are known, mathematical models based on transfer functions could be used to estimate tissue temperature at specific points with reasonable accuracy. This study could lead to the further development of more detailed mathematical models by considering additional factors, such as insertion depth and circulation blood flow.

## CONCLUSIONS

The results suggest that it is possible to estimate temperature at a specific point (at a depth of 4 mm in this case) from the evolution of the applied voltage when the initial value of the tissue characteristics is known. The proposed mathematical model agreed reasonably accurately (around 4ºC) with the results given by FEM.

## ACKNOWLEDGMENTS

This work received financial support from the Spanish “Plan Nacional de I+D+I del Ministerio de Ciencia e Innovación” Grant No. TEC2008-01369/TEC and FEDER Project MTM2010-14909.

## REFERENCES

### Endorsements

"Open access will revolutionize 21st century knowledge work and accelerate the diffusion of ideas and evidence that support just in time learning and the evolution of thinking in a number of disciplines."

Daniel Pesut
(Indiana University School of Nursing, USA)

"It is important that students and researchers from all over the world can have easy access to relevant, high-standard and timely scientific information. This is exactly what Open Access Journals provide and this is the reason why I support this endeavor."

Jacques Descotes
(Centre Antipoison-Centre de Pharmacovigilance, France)

"Publishing research articles is the key for future scientific progress. Open Access publishing is therefore of utmost importance for wider dissemination of information, and will help serving the best interest of the scientific community."

Patrice Talaga
(UCB S.A., Belgium)

"Open access journals are a novel concept in the medical literature. They offer accessible information to a wide variety of individuals, including physicians, medical students, clinical investigators, and the general public. They are an outstanding source of medical and scientific information."

Jeffrey M. Weinberg
(St. Luke's-Roosevelt Hospital Center, USA)

"Open access journals are extremely useful for graduate students, investigators and all other interested persons to read important scientific articles and subscribe scientific journals. Indeed, the research articles span a wide range of area and of high quality. This is specially a must for researchers belonging to institutions with limited library facility and funding to subscribe scientific journals."

Debomoy K. Lahiri
(Indiana University School of Medicine, USA)

"Open access journals represent a major break-through in publishing. They provide easy access to the latest research on a wide variety of issues. Relevant and timely articles are made available in a fraction of the time taken by more conventional publishers. Articles are of uniformly high quality and written by the world's leading authorities."

Robert Looney

"Open access journals have transformed the way scientific data is published and disseminated: particularly, whilst ensuring a high quality standard and transparency in the editorial process, they have increased the access to the scientific literature by those researchers that have limited library support or that are working on small budgets."

Richard Reithinger
(Westat, USA)

"Not only do open access journals greatly improve the access to high quality information for scientists in the developing world, it also provides extra exposure for our papers."

J. Ferwerda
(University of Oxford, UK)

"Open Access 'Chemistry' Journals allow the dissemination of knowledge at your finger tips without paying for the scientific content."

Sean L. Kitson
(Almac Sciences, Northern Ireland)

"In principle, all scientific journals should have open access, as should be science itself. Open access journals are very helpful for students, researchers and the general public including people from institutions which do not have library or cannot afford to subscribe scientific journals. The articles are high standard and cover a wide area."

Hubert Wolterbeek
(Delft University of Technology, The Netherlands)

"The widest possible diffusion of information is critical for the advancement of science. In this perspective, open access journals are instrumental in fostering researches and achievements."

Alessandro Laviano
(Sapienza - University of Rome, Italy)

"Open access journals are very useful for all scientists as they can have quick information in the different fields of science."

Philippe Hernigou
(Paris University, France)

"There are many scientists who can not afford the rather expensive subscriptions to scientific journals. Open access journals offer a good alternative for free access to good quality scientific information."

Fidel Toldrá
(Instituto de Agroquimica y Tecnologia de Alimentos, Spain)

"Open access journals have become a fundamental tool for students, researchers, patients and the general public. Many people from institutions which do not have library or cannot afford to subscribe scientific journals benefit of them on a daily basis. The articles are among the best and cover most scientific areas."

M. Bendandi
(University Clinic of Navarre, Spain)

"These journals provide researchers with a platform for rapid, open access scientific communication. The articles are of high quality and broad scope."

Peter Chiba
(University of Vienna, Austria)

"Open access journals are probably one of the most important contributions to promote and diffuse science worldwide."

Jaime Sampaio
(University of Trás-os-Montes e Alto Douro, Portugal)

"Open access journals make up a new and rather revolutionary way to scientific publication. This option opens several quite interesting possibilities to disseminate openly and freely new knowledge and even to facilitate interpersonal communication among scientists."

Eduardo A. Castro
(INIFTA, Argentina)

"Open access journals are freely available online throughout the world, for you to read, download, copy, distribute, and use. The articles published in the open access journals are high quality and cover a wide range of fields."

Kenji Hashimoto
(Chiba University, Japan)

"Open Access journals offer an innovative and efficient way of publication for academics and professionals in a wide range of disciplines. The papers published are of high quality after rigorous peer review and they are Indexed in: major international databases. I read Open Access journals to keep abreast of the recent development in my field of study."

Daniel Shek
(Chinese University of Hong Kong, Hong Kong)

"It is a modern trend for publishers to establish open access journals. Researchers, faculty members, and students will be greatly benefited by the new journals of Bentham Science Publishers Ltd. in this category."

Jih Ru Hwu
(National Central University, Taiwan)

Webmaster Contact: info@benthamopen.net