# Automatic control of finite element models for temperature-controlled radiofrequency ablation

- Dieter Haemmerich
^{1, 2}Email author and - John G Webster
^{3}

**4**:42

**DOI: **10.1186/1475-925X-4-42

© Haemmerich and Webster; licensee BioMed Central Ltd. 2005

**Received: **28 April 2005

**Accepted: **14 July 2005

**Published: **14 July 2005

## Abstract

### Background

The finite element method (FEM) has been used to simulate cardiac and hepatic radiofrequency (RF) ablation. The FEM allows modeling of complex geometries that cannot be solved by analytical methods or finite difference models. In both hepatic and cardiac RF ablation a common control mode is temperature-controlled mode. Commercial FEM packages don't support automating temperature control. Most researchers manually control the applied power by trial and error to keep the tip temperature of the electrodes constant.

### Methods

We implemented a PI controller in a control program written in C++. The program checks the tip temperature after each step and controls the applied voltage to keep temperature constant. We created a closed loop system consisting of a FEM model and the software controlling the applied voltage. The control parameters for the controller were optimized using a closed loop system simulation.

### Results

We present results of a temperature controlled 3-D FEM model of a RITA model 30 electrode. The control software effectively controlled applied voltage in the FEM model to obtain, and keep electrodes at target temperature of 100°C. The closed loop system simulation output closely correlated with the FEM model, and allowed us to optimize control parameters.

### Discussion

The closed loop control of the FEM model allowed us to implement temperature controlled RF ablation with minimal user input.

### Keywords

ablation liver ablation hepatic ablation radiofrequency ablation RF ablation finite element method## Background

*P*= constant), temperature-controlled mode (

*T*= constant) and impedance-controlled mode (

*Z*< constant) are commonly used. The electromagnetic energy is converted to heat by resistive heating. Tissue damage can occur at temperatures above 43°C with heating times of several hours; at 50°C cell necrosis occurs after ~3 min [12]. A commonly used mode is temperature-controlled ablation, where the tip temperature of the electrodes is kept at a predetermined value, usually around 100°C. The hepatic electrodes (Fig. 1) used for temperature-controlled ablation have temperature sensing elements (thermistors or thermocouples) embedded in the prong tips. The sensors report temperature back to the generator, which then applies an appropriate amount of power to the electrode to keep the temperature constant.

Researchers have been using the finite element method (FEM) to simulate both cardiac and hepatic RF ablation [3–11]. When using the FEM to model temperature-controlled ablation, the applied voltage has to be adjusted to keep the tip temperature constant. Previously, most researchers used manual adjustment of applied voltage and trial-and-error methods to perform temperature-controlled ablation [3, 5–8]. Implementation of temperature-controlled feedback in models is of importance to obtain results comparable to clinical devices that use this type of control. We modeled a clinically used hepatic ablation electrode (15 gauge, RITA medical systems model 30) as described previously [6]. We implemented a control algorithm for a PI controller – a commonly used controller type – in a C++ program to change the applied voltage between the time steps. We used a closed loop system simulation to optimize control parameters for the PI controller.

## Methods

### 1. Finite element method

RF ablation destroys tissue by thermal energy, which is converted from electric energy. The current flows from the conductive electrode through the tissue to a surface dispersive electrode. Tissue in close vicinity of the electrode tip is heated by resistive heating.

The heating of tissue during RF ablation is governed by the bioheat equation:

where *ρ* is the density (kg/m^{3}), *c* is the specific heat (J/(kg·K)), and *k* is the thermal conductivity (W/(m·K)). J is the current density (A/m^{2}) and E is the electric field intensity (V/m). *T*_{bl} is the temperature of blood, *ρ*_{bl} is the blood density (kg/m^{3}), *c*_{bl} is the specific heat of the blood (J/(kg·K)), and *w*_{bl} is the blood perfusion (1/s). *h*_{bl} is the convective heat transfer coefficient accounting for the blood perfusion. *Q*_{m} (W/m^{3}) is the energy generated by metabolic processes and was neglected since it is small compared to the other terms [13]. Equation (1) defines the solution in the spatial domain encompassing electrodes and tissue.

Tissue properties were assumed temperature independent and are described in more detail in [6]. We used the commercial software ABAQUS/Standard 6.3 (Hibbitt, Karlsson & Sorensen, Inc., Pawtucket, RI) for solving the coupled thermo-electrical analysis. All analysis was performed on a SUN Blade-1000 workstation equipped with 2.5 GB of RAM and 80 GB of hard disk space.

Fig. 1 shows the geometry of the RITA model 30, 4-prong electrode we used in our models. The electrode was placed within a cylinder (80 mm diameter, 50 mm length) of liver tissue. The outer surfaces of the cylinder were set to 37°C (thermal boundary condition), and 0 V (electrical boundary condition). We performed quasi-static analysis. Due to the symmetry of the arrangement, we could reduce computing time by only modeling a quarter of the cylinder. We used the same model as in a previous study [3]; the model consisted of ~35,000 tetrahedral elements and ~7,000 nodes. The node spacing was small next to the electrode (0.2 mm) and larger at the model boundary (2 mm). Perfusion was included in the model according to the Pennes model [14]. The blood perfusion *w* bl used in this model was 6.4·10^{-3} 1/s [15].

An input file was submitted to ABAQUS, in which geometry, material properties, step time and boundary conditions were specified. The applied voltage was one of the boundary conditions and remained constant during each step. ABAQUS created a results file in which temperatures of all nodes of the FEM model at the end of the step were written. We created a C++ program that reads the temperature at the electrode tip from the results file and creates a new input file where the applied voltage is set according to a control algorithm; the C++ program then calls ABAQUS and waits until the next step is completed. We implemented a PI controller in our control algorithm.

### 2. Closed loop system simulation

Since the FEM model takes several hours to complete, using the FEM model to optimize control parameters was not feasible. Therefore we used a closed loop simulation to adjust parameters of the controller, and then used these parameters in the FEM model. Before implementing the controller, we analyzed the dynamic system (i.e. the FEM model) to be controlled. This system consisted of the ablation electrode, tissue, and dispersive electrode. The input variable of the system was the voltage applied to the electrode. The output variable was the temperature measured at the tip of the electrode. Initially, we determined the step response by applying a constant voltage of 25 V for 180 s. We then approximated the transfer function of this system by a time-discrete transfer function of the following form:

*a*

_{0},

*a*

_{1},

*a*

_{2},

*b*

_{1}and

*b*

_{2}of (2). Fig. 2 shows the step responses of the original dynamic system (i.e. the FEM model) and of the approximation according to (2). The parameters used for the approximation in (2) were:

*a*

_{0}= 1.0,

*a*

_{1}= -0.959,

*a*

_{2}= 0.127,

*b*

_{1}= 1.734,

*b*

_{2}= -1.150. The sampling time of the approximated system was 10 s, since that was also the sampling time used later in the digital PI controller. This ensured a more accurate simulation model of the closed-loop control system.

*T*

_{s}is the desired set tip temperature and

*T*

_{t}is the current tip temperature. The input of the PI controller

*e*=

*T*

_{s}-

*T*

_{t}. The output of the PI controller

*u*(corresponds to the applied voltage) is fed into the dynamic system.

The transfer function of the PI controller is:

Equation (3) can be described in the time-domain by:

To implement this controller in software, we have to use the discrete time-domain version of (4),

The second term of (5) represents the approximation of the integral term in (4) by trapezoidal numerical integration. The behavior of the PI controller is determined by the two parameters *K*_{p} and *K*_{i}.

We simulated the behavior of the closed-loop system with the software ANA. From in-vivo experiments in pigs with the RITA 500 generator and the Model 30 electrode we determined that it takes between 1 and 2 min for the tip temperature to reach a target temperature of 100°C, which is the temperature used in clinical cases at our institution. The RITA 500 generator does not allow recording of the tip temperature signal, so we don't have exact data of the tip temperature over time. We empirically chose parameters for the PI controller to minimize overshoot and obtain similar temporal behavior as in the in-vivo experiments.

The parameters we used for the PI controller were *K*_{p} = 0.02, *K*_{i} = 0.0064.

### 3. Controlled FEM model

We implemented the PI controller with the parameters resulting from the control system simulation in a control program written in C++. The software control program determined the applied voltage and the step time. An input file for ABAQUS was created and ABAQUS solved the FEM model. ABAQUS created a result file, which included the temperatures at all nodes.

The tip temperature was read from the result file by the control program. The program then determined the applied voltage and step size for the next step, and ABAQUS was restarted with the modified input file. With decreasing step size (i.e. time that ABAQUS simulated the model with constant voltage) total computation time gets longer. As a compromise, we chose 10 s as the initial step size. Note that the FEM solver divided each step into smaller increments, starting at 0.05s. Also, the FEM solver performed convergence tests to ensure that the increment size was sufficiently small.

We simulated ablation for 12 min using the FEM model and the control program. The simulation took 250 min to run and was divided into 35 steps. The individual steps took between 6 and 11 min to run in ABAQUS. Without using the control program, manual interaction would be necessary after each step to manually change the input file for ABAQUS and apply a new input voltage and set the step size. In this case the operator would have to check the results file after each step (i.e. in 6 to 11 min intervals), modify the input file and restart ABAQUS.

## Results and Discussion

In hepatic RF ablation, ablation times clinically used go up to 35 min. Since the heat-up period is comparably much smaller (1 to 2 min.), it is not of essential importance that the temporal behavior of the control algorithm reproduces the control algorithm used in clinical devices during the heat-up period. From our experience, the temperature distribution in the FEM model reaches close to steady state at the end of the simulation due to the long simulation times. As long as the tip temperature is kept within a small range around the target temperature after the initial heat-up period, the model results (i.e. final temperature distribution) should not differ significantly. However, with knowledge of the actual control parameters and algorithms of commercial devices (e.g. obtained from measurements), accurate simulations of these devices is possible using our methods.

## Conclusion

We implemented a closed loop control system to a FEM model, to automate simulation of temperature-controlled RF ablation. We further used a closed loop control system simulation to optimize control parameters. Previously, researchers often applied constant power, or used time-consuming trial-and-error methods to determine required voltage. Furthermore, if control parameters and algorithms of commercial devices are known or can be measured, an accurate simulation of commercial devices is possible.

## Declarations

### Acknowledgements

This work was supported by NIH Grant HL56413.

## Authors’ Affiliations

## References

- Bosch FX, Ribes J, Borras J: Epidemiology of primary liver cancer.
*Semin Liver Dis*1999, 19: 271–85.View ArticleGoogle Scholar - Curley SA, Cusack JC Jr, Tanabe KK, Stoelzing O, Ellis LM: Advances in the treatment of liver tumors.
*Curr Probl Surg*2002, 39: 449–571. 10.1067/msg.2002.122810View ArticleGoogle Scholar - Haemmerich D, Staelin ST, Tungjitkusolmun S, Lee FT Jr, Mahvi DM, Webster JG: Hepatic bipolar radio-frequency ablation between separated multiprong electrodes.
*IEEE Trans Biomed Eng*2001, 48: 1145–52. 10.1109/10.951517View ArticleGoogle Scholar - Tungjitkusolmun S, Woo EJ, Cao H, Tsai JZ, Vorperian VR, Webster JG: Finite element analyses of uniform current density electrodes for radio-frequency cardiac ablation.
*IEEE Trans Biomed Eng*2000, 47: 32–40. 10.1109/10.817617View ArticleGoogle Scholar - Panescu D, Whayne JG, Fleischman SD, Mirotznik MS, Swanson DK, Webster JG: Three-dimensional finite element analysis of current density and temperature distributions during radio-frequency ablation.
*IEEE Trans Biomed Eng*1995, 42: 879–890. 10.1109/10.412649View ArticleGoogle Scholar - Tungjitkusolmun S, Staelin ST, Haemmerich D, Tsai JZ, Cao H, Webster JG, Lee FT, Mahvi DM, Vorperian VR: Three-dimensional finite-element analyses for radio-frequency hepatic tumor ablation.
*IEEE Trans Biomed Eng*2002, 49: 3–9. 10.1109/10.972834View ArticleGoogle Scholar - Jain MK, Wolf PD: Temperature-controlled and constant-power radio-frequency ablation: what affects lesion growth?
*IEEE Trans Biomed Eng*1999, 46: 1405–1412. 10.1109/10.804568View ArticleGoogle Scholar - Chang IA, Nguyen UD: Thermal modeling of lesion growth with radiofrequency ablation devices.
*Biomed Eng Online*2004, 3: 27. 10.1186/1475-925X-3-27View ArticleGoogle Scholar - Chang I: Finite element analysis of hepatic radiofrequency ablation probes using temperature-dependent electrical conductivity.
*Biomed Eng Online*2003, 2: 12. 10.1186/1475-925X-2-12View ArticleGoogle Scholar - Liu Z, Lobo SM, Humphries S, Horkan C, Solazzo SA, Hines-Peralta AU, Lenkinski RE, Goldberg SN: Radiofrequency tumor ablation: insight into improved efficacy using computer modeling.
*AJR Am J Roentgenol*2005, 184: 1347–52.View ArticleGoogle Scholar - Gopalakrishnan J: A mathematical model for irrigated epicardial radiofrequency ablation.
*Ann Biomed Eng*2002, 30: 884–93. 10.1114/1.1507845View ArticleGoogle Scholar - Graham SJ, Chen L, Leitch M, Peters RD, Bronskill MJ, Foster FS, Henkelman RM, Plewes DB: Quantifying tissue damage due to focused ultrasound heating observed by MRI.
*Magn Reson Med*1999, 41: 321–8. 10.1002/(SICI)1522-2594(199902)41:2<321::AID-MRM16>3.0.CO;2-9View ArticleGoogle Scholar - Chato JC: Heat transfer to blood vessels.
*J Biomech Eng*1980, 102: 110–8.View ArticleGoogle Scholar - Pennes HH: Analysis of tissue and arterial blood temperatures in the resting human forearm.
*J Appl Physiol*1948, 1: 93–122.Google Scholar - Ebbini ES, Umemura S, Ibbini M, Cain CA: A Cylindrical-section ultrasound phased-array applicator for hyperthermia cancer-therapy.
*IEEE Trans Ultrason Ferroelectr Freq Control*1988, 35: 561–572. 10.1109/58.8034View ArticleGoogle Scholar

## Copyright

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.