The aim of the study was to design a novel radiofrequency (RF) electrode for larger and rounder ablation volumes and its ability to achieve the complete ablation of liver tumors larger than 3 cm in diameter using finite element method. A new RF expandable electrode comprising three parts (i.e., insulated shaft, changing shaft, and hooks) was designed. Two modes of this new electrode, such as monopolar expandable electrode (MEE) and hybrid expandable electrode (HEE), and a commercial expandable electrode (CEE) were investigated using liver tissue with (scenario I) and without (scenario II) a liver tumor. A temperature-controlled radiofrequency ablation (RFA) protocol with a target temperature of 95 and an ablation time of 15 min was used in the study. Both the volume and shape of the ablation zone were examined for all RF electrodes in scenario I. Then, the RF electrode with the best performance in scenario I and CEE were used to ablate a large liver tumor with the diameter of 3.5 cm (scenario II) to evaluate the effectiveness of complete tumor ablation of the designed RF electrode. In scenario I, the ablation volumes of CEE, HEE, and MEE were 12.11 cm3, 33.29 cm3, and 48.75 cm3, respectively. The values of sphericity index (SI) of CEE, HEE, and MEE were 0.457, 0.957, and 0.976, respectively. The best performance was achieved by using MEE. In scenario II, the ablation volumes of MEE and CEE were 71.59 cm3 and 19.53 cm3, respectively. Also, a rounder ablation volume was achieved by using MEE compared to CEE (SI: 0.978 versus 0.596). The study concluded that: (1) compared with CEE, both MEE and HEE get larger and rounder ablation volumes due to the larger electrode–tissue interface and rounder shape of hook deployment; (2) MEE has the best performance in getting a larger and rounder ablation volume; and (3) computer simulation result shows that MEE is also able to ablate a large liver tumor (i.e., 3.5 cm in diameter) completely, which has at least 0.785 cm safety margin.
Introduction
Radiofrequency ablation (RFA) as a thermal ablative method is used to treat focal primary and secondary malignancies in organs including the liver, lungs, kidneys, bones, and adrenal glands [1–3]. It has been shown that the clinical effectiveness of RFA drops off sharply for tumor diameters above 3 cm [4]. As a result, the clinical application of RFA is limited to small tumors (i.e., <3 cm in diameter) with a single electrode [5], called small tumor problem in the study. During the RFA procedure, tissue water boiling, vaporization, and carbonization usually start sequentially to occur when the tissue temperature is over 100 , which results in an immediate increase in impedance and subsequently a drop in the radiofrequency (RF) power that can be delivered [6,7].
To overcome the small tumor problem, we designed an expandable electrode to enlarge the interface between the tissue and electrodes to distribute the power to a wider zone [8]. So far, several similar expandable electrodes have been designed for the same reason, such as the umbrella electrode [8], the “Christmas tree” electrode [9], the InCircle Monarch electrode [10], the origami electrode [11], and the bicomponent conformal electrode [12]. Yet it is still a challenge to get a large (>3 cm in diameter) and round (close to a sphere) ablation zone for those expandable electrodes due to the shape and the difficulty in the deployment of hooks.
The purpose of the study was to evaluate the ablation results of our novel expandable electrode applied in two different modes compared to a commercial expandable electrode (CEE) in terms of the enlargement of the ablation volume in a spherical shape.
Materials and Methods
Design of a New RF Electrode.
A 15-gauge CEE (RITA Medical Systems, Mountain View, CA) was chosen as a reference RF electrode [13]. The CEE contains four curved hooks that can be deployed from the electrode hollow in an umbrellalike shape (3 cm in diameter with full deployment) shown as in Fig. 1(a). In addition to the four hooks, the CEE also contains a 1-cm distal part as an active tip. Considering the shape of ablation zone generated by the CEE, an electrode with more spherical shape of deployed hooks and larger tissue–electrode interface was designed in the study. As shown in Fig. 1(b), the designed electrode has three parts, such as part A: six active hooks, part B: changing shaft, and part C: insulated part. The dimensions of the designed electrode are given in Fig. 1(b). For the designed electrode, there were two working modes considered in the present study: mode 1 (active tips: parts A and B, grounding pads: the side and bottom boundaries of liver tissue (Fig. 2). For this mode, the designed electrode was termed as monopolar expandable electrode (MEE)) and mode 2 (active tip: part A, grounding pads: part B and the side and bottom boundaries of liver tissue. In this mode, the designed electrode was termed as hybrid expandable electrode (HEE)).
Finite Element Model of RFA.
In the study, a three-dimensional finite element model (FEM) was used for RFA simulation. As shown in Fig. 2, a cylinder with 100 mm in diameter and 120 mm in height was expressed as the liver tissue. The RF electrode was placed at the center of liver with an insert depth of 60 mm.
where is the current density, is the electrical field intensity, is the electrical conductivity, and is the applied voltage.
where is the constant blood perfusion of tissue (refer to Table 1).
Property | Liver | Tumor | Active tip | Insulated shaft | Blood |
---|---|---|---|---|---|
1080 | 1045 | 6450 | 70 | 1000 | |
3455 | 3760 | 840 | 1045 | 4180 | |
0.515a | 0.60a | 18 | 0.026 | 0.49 | |
0.203a | 0.50a | 1.0108 | 1.010−5 | 0.667 | |
0.016a | 0.002a | — | — | — |
Property | Liver | Tumor | Active tip | Insulated shaft | Blood |
---|---|---|---|---|---|
1080 | 1045 | 6450 | 70 | 1000 | |
3455 | 3760 | 840 | 1045 | 4180 | |
0.515a | 0.60a | 18 | 0.026 | 0.49 | |
0.203a | 0.50a | 1.0108 | 1.010−5 | 0.667 | |
0.016a | 0.002a | — | — | — |
Evaluated at 21 .
where is the frequency factor, is the activation energy barrier, is the universal gas constant (8.314 ), and is the absolute temperature. For the liver and tumor tissues, and , respectively, and and , respectively [20]. The threshold of tissue death is , which means a 63% probability of cell death (D63) [21].
Initial and Boundary Conditions.
The initial temperature of the whole tissue domain was set at 37 (inner body temperature). A zero-voltage electric potential was applied at the grounding pads. The boundaries of the active tip were set to a variable power source determined by a proportional–integrative–derivative controller (, , and ). The up boundary of liver tissue was set as electrically isolative. The other electrical and thermal boundary conditions of the FEM were considered as continuous.
Mesh Convergence Test.
The software of comsol multiphysics 5.2a (COMSOL, Inc., Burlington, MA) was used to simulate and predict the ablation zone. The CEE was used to perform the mesh convergence test. In the convergence test, RFA was applied at the constant applied voltage of 30 V and the ablation time of 200 s in scenario I. The mesh convergence criterion was considered as the difference in the maximum temperature between two contiguous meshes smaller than 0.1%. The free tetrahedral mesh was used to discretize the physical domains of the model, which was divided as two parts: part I (active tips, insulated shaft, and the tumor sphere) and part II (the remaining parts of the model). The mesh of part II was fixed at the default “normal” in comsol; the mesh of part I was set as the default “finer” with the minimum element size of 0.5 mm and further refined until convergence criterion got satisfied. Then, the previous mesh with fewer elements was used to reduce the computational load. The meshing results of elements are given in Table 2.
Validation of the FEM Model.
The results of in vivo or ex vivo experiments of RFA using the 15-gauge CEE from Ref. [13] were chosen as the references. The same situations with the work of de Baere et al. [13] were considered in the validation, such as the dimensions of RF electrode, boundary conditions, and anatomical environments. However, it is worth mentioning that the human physical properties of liver tissue were used in the study, unlike the porcine and bovine studies of de Baere et al. Both in the work of Hall et al. [22] and Trujillo and Berjano [15] concluded that properties affecting the ablation volume in computer simulation significantly are the electrical conductivity and perfusion ratio of liver tissue. By examining the liver electrical conductivities of human, pigs, and cows (0.467–0.512, 0.465–0.528, and 0.488 , respectively [22–26]) and the perfusion ratios of the human and porcine liver (0.009–0.019 and 0.018 , respectively [22,27]), we believe that the work of de Baere et al. is reliable to be used as the reference.
where and are the diameters perpendicular to the axis of the electrode, and is the diameter along it. In the literature, de Baere et al. assumed that due to the technical difficulty of measurement. For the ablation size of the FEM, the ablation volume was calculated by using Eq. (6) (), while , , and were measured by the software along the directions mentioned previously. As shown in Table 3, there was good agreement between ablation size found for FEM and the experiments. It is also worth mentioning that (in vivo and ex vivo) in the FEM is slightly larger than its peer in the experiment. The reasons for this might be the measurement error in the experiment and the blood perfusion error in simulation. Nonetheless, it is reasonable to accept the proposed FEM of RFA as a test-bed for ablation volume estimation.
FEM | Experiment | |||
---|---|---|---|---|
Ex vivo | In vivo | Ex vivo | In vivo | |
(cm) | 4.51 | 3.70 | 3.94 ± 0.38 | 2.86 ± 0.4 |
(cm) | 3.94 | 2.55 | 3.94 ± 0.38 | 2.86 ± 0.4 |
(cm) | 3.82 | 2.75 | 3.46 ± 0.5 | 2.55 ± 0.44 |
(cm3) | 34.20 | 12.11 | 26.67 ± 9.59 | 11.16 ± 3.65 |
FEM | Experiment | |||
---|---|---|---|---|
Ex vivo | In vivo | Ex vivo | In vivo | |
(cm) | 4.51 | 3.70 | 3.94 ± 0.38 | 2.86 ± 0.4 |
(cm) | 3.94 | 2.55 | 3.94 ± 0.38 | 2.86 ± 0.4 |
(cm) | 3.82 | 2.75 | 3.46 ± 0.5 | 2.55 ± 0.44 |
(cm3) | 34.20 | 12.11 | 26.67 ± 9.59 | 11.16 ± 3.65 |
Analysis of Sphericity.
where is the maximum value among , , and .
Study Design.
A proportional–integrative–derivative controller was implemented in the comsol software to control the temperature at a hook tip (Fig. 1). The goal was to achieve a temperature of 95 for 15 min [13]. As mentioned before, two scenarios were considered in the present study. Scenario I: three RF electrodes (i.e., CEE, MEE, and HEE) were compared in terms of size and the sphericity of the ablation volume using one-compartment RFA model. Scenario II: the electrode that achieved the largest ablation volume and the highest SI was compared to the CEE in the treatment of a liver tumor with a diameter of 3.5 cm. Successful ablation was defined using the clinical standard of complete ablation of the ablation of the liver tumor and ablation of a 0.5–1 cm margin of normal tissue surrounding the tumor.
Results
Results of Scenario I.
Figures 3(a)–3(c) show the potential distribution of CEE, HEE, and MEE at the end of treatment (15 min), respectively. The temperature distribution of CEE, HEE, and MEE is shown in Figs. 3(d)–3(f), respectively. It is worth mentioning that the high-temperature areas (>90 ) are different among the three electrodes, as shown in Figs. 3(d)–3(f). For CEE and MEE, the maximum temperature (i.e., 95.5 and 99.6 , respectively) was at the area around the hooks. For HEE, however, the maximum temperature (i.e., 106.7 ) was at the area around part B of the electrode. Figure 4 shows the results of ablation zone generated by the three RF electrodes (i.e., CEE, HEE, and MEE). From Fig. 4, we can find that the ablation zones generated by HEE and MEE were much better than that generated by CEE in both the ablation volume and shape. The ablation volumes of CEE, HEE, and MEE were 12.11 cm3, 33.29 cm3, and 48.75 cm3, respectively. Compared with CEE, the ablation volumes of HEE and MEE increased by 174.9% and 302.6%, respectively. The same conclusion was also able to be found in the shape of ablation volume. As shown in Table 4, the ablation volumes generated by HEE and MEE were very close to spheres (SI = 0.957 and 0.976, respectively). However, the ablation volume generated by CEE was in an umbrella shape with SI = 0.457. MEE was found to be the best performing electrode with the ablation volume of 48.75 cm3 and SI = 0.976. Thus, MEE was chosen to be used in scenario II.
Results of Scenario II.
In scenario II, a liver tumor with the diameter of 3.5 cm was taken into account. CEE and MEE were used to treat the liver tumor in a clinical way (15 min and temperature-controlled RFA at 95 ). Table 5 shows the results of ablation zones and sphericity indices generated by CEE and MEE. Compared with CEE, the sizes of ablation zone generated by MEE were increased by 59.4%, 30.7%, and 61.8% in the diameters of , , and , respectively. The ablation volume was increased by 266.6% from 19.53 cm3 to 71.59 cm3. More spherical ablation zone was achieved by using MEE (SI = 0.978) compared with CEE (SI = 0.596). Furthermore, at least 0.785 cm ablative margin surrounding the liver tumor can be obtained in all the diameters of , , and , as shown in Table 5 and Fig. 5.
Discussion
In scenario I, both HEE and MEE can achieve larger and rounder ablation volumes compared to CEE. We attribute this to the larger electrode–tissue interfaces and the rounder shape of hooks. Furthermore, a larger and rounder ablation volume can be achieved by MEE compared to HEE (ablation volume: 48.75 cm3 versus 33.29 cm3 and SI: 0.976 versus 0.957). This phenomenon may be attributed to the following two reasons. The first reason is the role of part B of the designed electrode. In the case of HEE, the alternating current not only flows to the tissue bottom and side boundaries but also flows to part B of the designed electrode, as shown in Fig. 1(b). A smaller potential difference at the ablation boundaries can be found in HEE compared to MEE, shown in Figs. 3(b) and 3(c). Such a difference is due to a reduction in power delivered to the tissue areas outside the hooks, which results in a smaller ablation diameter perpendicular to the electrode compared with MEE. However, there is no significant difference in the ablation diameter along the electrode. Therefore, a larger and rounder ablation volume is achieved in MEE. The second reason is the position of the maximum temperature. For HEE, the maximum temperature is found in the area surrounding part B, which is different from MEE where the maximum temperature is around the hooks. This phenomenon is because of the tremendously Joule heat generated by RF power around part B due to the larger potential difference. Therefore, a greater amount of RF power is delivered to the tissue areas inner the hooks in HEE than MEE.
It is noticed that the maximum temperature (106.7 ) occurred in the area surrounding part B in the case of HEE, while the temperature at one of the hook tips was controlled at 95 . This phenomenon was believed as the reason of dramatically decrease in potential in the area surrounding part B that was set as grounding pad (Fig. 3(b)). Thus, the tissue dehydration might be the case in the same area, which can lead to a certain degree of error in the calculation of ablation volume. Due to the target temperature (95 ) and the ablation time used in the study, only a very limited area of tissue surrounding part B was dehydrated (Fig. 3(e)). So the error in ablation volume should be negligible in the study. Otherwise, it is necessary to consider the tissue phase change in the FEM as it affects the accuracy of ablation volume.
In scenario II, due to the tumor characteristics, both HEE and MEE generate larger ablation volumes than that in scenario I. For CEE, only a fan-shaped ablation zone is generated along the section view of cross-electrode-axis (Fig. 5(a)), which cannot completely cover the tumor area. Because there are only four hooks on this commercial electrode, the size of ablation zone between the two adjacent hooks () is much smaller than the size between the two diagonal hooks () (Fig. 5(c)). And the semicircle hooks of CEE cannot achieve a sphere ablation volume (Fig. 5(a)). As shown in Figs. 5(b) and 5(d), MEE generates the ablation volume with at least 0.785 cm safety margin, which supports the assertion that our novel monopolar electrode has the potential to ablate tumors larger than 3 cm completely.
A limitation of this preliminary study is the lack of experimental results. Future studies are warranted to the fabrication of the design RF electrode and its application using tissue-mimicking phantom followed by ex vivo/in vivo porcine liver tissue. The difficulty in the deployment of hook can be imaged due to its spherical and long shape. Improvements on the design of hook or the analysis of forces between hook and biological tissue during its deployment must be investigated in the future work.
Conclusion
A novel expandable electrode was designed and assessed using finite element model in the study, and the following three conclusions can be drawn:
- (1)
Compared with CEE, both MEE and HEE can achieve larger and rounder ablation volumes due to the larger electrode–tissue interface and the rounder shape of hooks.
- (2)
Compared with HEE, MEE yielded a larger and rounder ablation volume.
- (3)
MEE was able to completely ablate a large liver tumor (i.e., 3.5 cm in diameter) with at least 0.785 cm safety margin.
Funding Data
National Sciences and Engineering Research Council of Canada (NSERC) through a discovery grant to Wenjun Zhang.
National Natural Science of China (Grant No. 51375166).
Nomenclature
- =
frequency factor
- =
specific heat of liver
- =
specific heat of blood
- CEE =
commercial expandable electrode
- =
longest diameters perpendicular to the axis of the electrode
- =
shortest diameters perpendicular to the axis of the electrode
- =
diameter along the electrode
- =
electrical field intensity
- FEM =
finite element model
- HEE =
hybrid expandable electrode
- =
current density
- =
thermal conductivity of liver tissue
- =
thermal conductivity measured at 21 °C
- =
coefficient for the derivative term
- =
coefficient for the integrative term
- =
coefficient for the proportional term
- MEE =
monopolar expandable electrode
- =
energy generated by the RF electrical current
- =
energy generated due to metabolic processes
- =
universal gas constant
- RF =
radiofrequency
- RFA =
radiofrequency ablation
- SI =
sphere index
- t =
ablation time
- =
temperature of liver
- =
absolute temperature which is a function of ablation time
- =
temperature of the blood entering the tissue
- =
applied voltage
- =
ablation volume
- =
Cartesian coordinate system
- =
activation energy barrier
- =
liver density
- =
blood density
- =
electrical conductivity of liver tissue
- =
electrical conductivity measured at 21 °C
- =
degree of tissue death
- =
blood perfusion
- =
constant blood perfusion of tissue
- =
analyzed spatial domains