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 °C 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 [13]. 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 °C, 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)).

Fig. 1
Dimensions (in mm) and structures of two RF electrodes: (a) CEE and (b) the designed electrode
Fig. 1
Dimensions (in mm) and structures of two RF electrodes: (a) CEE and (b) the designed electrode
Close modal
Fig. 2
Three-dimensional geometric model (in mm) of liver tissue with the liver tumor and the designed RF electrode
Fig. 2
Three-dimensional geometric model (in mm) of liver tissue with the liver tumor and the designed RF electrode
Close modal

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.

The governing equation of the RFA procedure follows the Penne’s bioheat transfer equation [1]
(1)
where ρ(kg m-3) is the density, c (J kg-1 K-1) is the specific heat, T(x,t)(°C) is the temperature, k (W m-1K-1) is the thermal conductivity, ρb(kg m-3) is the blood density, cb(J kg-1 K-1) is the specific heat of the blood, ωb(s-1) is the blood perfusion, Tb is the temperature of the blood entering the tissue (assumed to be 37 °C), x={x,y,z} in the Cartesian coordinate system, denotes the analyzed spatial domains, Qm(x,t)(W m-3) is the energy generated due to metabolic processes which has been neglected since metabolic heat source is magnitude less than the spatial heat, and Qhs(x,t)(W m-3) generated by the RF electrical current as follows [14]:
(2)

where J (Am-2) is the current density, E (Vm-1) is the electrical field intensity, σ(Sm-1) is the electrical conductivity, and U(V) is the applied voltage.

The values of the physical properties in the FEM adopted from the literature [9,15,16] are tabulated in Table 1. The present study assumed the thermal conductivity and electrical conductivity of both liver and tumor tissues at the temperature below 100 °C as follows [17]:
(3)
(4)
where k0 and σ0 are the thermal conductivity and electrical conductivity, respectively, measured at the temperature at 21 °C (T0). Besides, when the temperature is above 100 °C, the thermal conductivity was taken as a constant. The electrical conductivity decreases by 2 orders of the magnitude from 100 °C to 105 °C, and a constant value at the temperature above 105 °C [15]. A piecewise function was used to describe the situation of blood perfusion in liver and tumor tissues [18]
(5)

where ωb0 is the constant blood perfusion of tissue (refer to Table 1).

Table 1

Physical properties of the materials used in the model

PropertyLiverTumorActive tipInsulated shaftBlood
ρ(kgm-3)108010456450701000
c(Jkg-1K-1)3455376084010454180
k(Wm-1K-1)0.515a0.60a180.0260.49
σ(Sm-1)0.203a0.50a1.0×1081.0×10−50.667
ωb(s-1)0.016a0.002a
PropertyLiverTumorActive tipInsulated shaftBlood
ρ(kgm-3)108010456450701000
c(Jkg-1K-1)3455376084010454180
k(Wm-1K-1)0.515a0.60a180.0260.49
σ(Sm-1)0.203a0.50a1.0×1081.0×10−50.667
ωb(s-1)0.016a0.002a
a

Evaluated at 21 °C.

In the present study, the induced thermal damage was evaluated by the Arrhenius model for the estimation of ablation zones [19]
(6)

where A(s-1) is the frequency factor, ΔE(Jmol-1) is the activation energy barrier, R (Jmol-1K-1) is the universal gas constant (8.314 Jmol-1K-1), and Tτ (K) is the absolute temperature. For the liver and tumor tissues, A=7.390×1039 s-1 and 3.247×1043s-1, respectively, and ΔE=2.577×105Jmol-1 and 2.814×105J mol-1, respectively [20]. The threshold of tissue death is Ωt=1, 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 °C (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 (KP=32, KI=1, and KD=0.001). 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.

Table 2

Number of elements for each scenario

ScenarioModelNumber of elements
One-compartment modelCEE231,483
MEE426,934
HEE426,934
Two-compartment modelCEE231,483
MEE426,934
ScenarioModelNumber of elements
One-compartment modelCEE231,483
MEE426,934
HEE426,934
Two-compartment modelCEE231,483
MEE426,934

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 S m-1, respectively [2226]) and the perfusion ratios of the human and porcine liver (0.009–0.019 and 0.018 s-1, respectively [22,27]), we believe that the work of de Baere et al. is reliable to be used as the reference.

The ablation volume was calculated approximately by using the following equation:
(7)

where d1 and d2 are the diameters perpendicular to the axis of the electrode, and d3 is the diameter along it. In the literature, de Baere et al. assumed that d1=d2 due to the technical difficulty of measurement. For the ablation size of the FEM, the ablation volume was calculated by using Eq. (6) (Ωt1), while d1, d2, and d3 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 d1 (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.

Table 3

Comparisons of ablation size from the FEM and ex vivo/in vivo experiments

FEMExperiment
Ex vivoIn vivoEx vivoIn vivo
d1 (cm)4.513.703.94 ± 0.382.86 ± 0.4
d2 (cm)3.942.553.94 ± 0.382.86 ± 0.4
d3 (cm)3.822.753.46 ± 0.52.55 ± 0.44
V (cm3)34.2012.1126.67 ± 9.5911.16 ± 3.65
FEMExperiment
Ex vivoIn vivoEx vivoIn vivo
d1 (cm)4.513.703.94 ± 0.382.86 ± 0.4
d2 (cm)3.942.553.94 ± 0.382.86 ± 0.4
d3 (cm)3.822.753.46 ± 0.52.55 ± 0.44
V (cm3)34.2012.1126.67 ± 9.5911.16 ± 3.65

Analysis of Sphericity.

In the present study, the sphericity of the ablation volume was also investigated. A near-spherical shape is desirable because most tumors are more or less round [28]. Clinically, it is important to know the shape of the proposed ablation zone so that the electrodes can be placed in such a way as to optimize the ablation zone relative to the tumor. A sphericity index (SI) was used to assess the shape of the ablation volumes for the three RF electrodes [29]
(8)

where maxd1,d2,d3 is the maximum value among d1, d2, and d3.

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 °C 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 °C) 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 °C and 99.6 °C, respectively) was at the area around the hooks. For HEE, however, the maximum temperature (i.e., 106.7 °C) 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.

Fig. 3
Potential distributions (in V) of (a) CEE, (b) HEE, and (c) MEE and temperature distributions (in °C) of (d) CEE, (e) HEE, and (f) MEE
Fig. 3
Potential distributions (in V) of (a) CEE, (b) HEE, and (c) MEE and temperature distributions (in °C) of (d) CEE, (e) HEE, and (f) MEE
Close modal
Fig. 4
Liver tissue death rates of the three RF electrodes: section view crossing the axis of (a) CEE, (b) HEE, and (c) MEE, respectively, and section view perpendicular to the axis of (d) CEE, (e) HEE, and (f) MEE, respectively, at the dashed line (the black circle means the liver tumor)
Fig. 4
Liver tissue death rates of the three RF electrodes: section view crossing the axis of (a) CEE, (b) HEE, and (c) MEE, respectively, and section view perpendicular to the axis of (d) CEE, (e) HEE, and (f) MEE, respectively, at the dashed line (the black circle means the liver tumor)
Close modal
Table 4

Ablation volumes and SIs of three electrodes

CEEHEEMEE
V (cm3)12.1133.2948.75
SI0.4570.9570.976
CEEHEEMEE
V (cm3)12.1133.2948.75
SI0.4570.9570.976

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 °C). 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 d1, d2, and d3, 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 d1, d2, and d3, as shown in Table 5 and Fig. 5.

Fig. 5
Tumor and liver tissue death rates in scenario II: section view crossing the axis of (a) CEE and (b) MEE, respectively, and section view perpendicular to the axis of (c) CEE and (d) MEE, respectively, at the dashed line (the black circle means the liver tumor)
Fig. 5
Tumor and liver tissue death rates in scenario II: section view crossing the axis of (a) CEE and (b) MEE, respectively, and section view perpendicular to the axis of (c) CEE and (d) MEE, respectively, at the dashed line (the black circle means the liver tumor)
Close modal
Table 5

Ablation results of CEE and MEE for a liver tumor with 3.5 cm in diameter

Electroded1 (cm)d2 (cm)d3 (cm)V (cm3)SIComplete ablation
CEE3.973.183.1419.530.596No
MEE5.195.075.0871.590.978Yes
Electroded1 (cm)d2 (cm)d3 (cm)V (cm3)SIComplete ablation
CEE3.973.183.1419.530.596No
MEE5.195.075.0871.590.978Yes

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 °C) 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 °C. 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 °C) 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 (d2) is much smaller than the size between the two diagonal hooks (d1) (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. (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. (2)

    Compared with HEE, MEE yielded a larger and rounder ablation volume.

  3. (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

A =

frequency factor

c =

specific heat of liver

cb =

specific heat of blood

CEE =

commercial expandable electrode

d1 =

longest diameters perpendicular to the axis of the electrode

d2 =

shortest diameters perpendicular to the axis of the electrode

d3 =

diameter along the electrode

E =

electrical field intensity

FEM =

finite element model

HEE =

hybrid expandable electrode

J =

current density

k =

thermal conductivity of liver tissue

k0 =

thermal conductivity measured at 21 °C

KD =

coefficient for the derivative term

KI =

coefficient for the integrative term

KP =

coefficient for the proportional term

MEE =

monopolar expandable electrode

Qhs =

energy generated by the RF electrical current

Qm =

energy generated due to metabolic processes

R =

universal gas constant

RF =

radiofrequency

RFA =

radiofrequency ablation

SI =

sphere index

t =

ablation time

T =

temperature of liver

Tτ =

absolute temperature which is a function of ablation time

Tb =

temperature of the blood entering the tissue

U =

applied voltage

V =

ablation volume

x =

Cartesian coordinate system

ΔE =

activation energy barrier

ρ =

liver density

ρb =

blood density

σ =

electrical conductivity of liver tissue

σ0 =

electrical conductivity measured at 21 °C

Ω  =

degree of tissue death

ωb =

blood perfusion

ωb0 =

constant blood perfusion of tissue

=

analyzed spatial domains

References

1.
Zhang
,
B.
,
Moser
,
M. A.
,
Zhang
,
E. M.
,
Luo
,
Y.
,
Zhang
,
H.
, and
Zhang
,
W.
,
2014
, “
Study of the Relationship Between the Target Tissue Necrosis Volume and the Target Tissue Size in Liver Tumours Using Two-Compartment Finite Element RFA Modelling
,”
Int. J. Hyperthermia
,
30
(
8
), pp.
593
602
.
2.
Ahmed
,
M.
, and
Goldberg
,
S. N.
,
2012
, “
Principles of Radiofrequency Ablation
,”
Interventional Oncology
, Springer, New York, pp.
23
37
.
3.
Ni
,
Y.
,
Mulier
,
S.
,
Miao
,
Y.
,
Michel
,
L.
, and
Marchal
,
G.
,
2005
, “
A Review of the General Aspects of Radiofrequency Ablation
,”
Abdom. Imaging
,
30
(
4
), pp.
381
400
.
4.
Zhang
,
B.
,
Moser
,
M. A.
,
Zhang
,
E. M.
,
Luo
,
Y.
,
Liu
,
C.
, and
Zhang
,
W.
,
2016
, “
A Review of Radiofrequency Ablation: Large Target Tissue Necrosis and Mathematical Modelling
,”
Phys. Med.
,
32
(
8
), pp.
961
971
.
5.
Romero-Méndez
,
R.
, and
Berjano
,
E.
,
2017
, “
An Analytical Solution for Radiofrequency Ablation With a Cooled Cylindrical Electrode
,”
Math. Probl. Eng.
,
2017
, p.
9021616
.
6.
Goldberg
,
S. N.
,
2001
, “
Radiofrequency Tumor Ablation: Principles and Techniques
,”
Eur. J. Ultrasound
,
13
(
2
), pp.
129
147
.
7.
Zhang
,
B.
,
Moser
,
M. A.
,
Luo
,
Y.
,
Zhang
,
E. M.
, and
Zhang
,
W.
,
2014
, “
Evaluation of the Current Radiofrequency Ablation Systems Using Axiomatic Design Theory
,”
Proc. Inst. Mech. Eng., Part H
,
228
(
4
), pp.
397
408
.
8.
Mulier
,
S.
,
Miao
,
Y.
,
Mulier
,
P.
,
Dupas
,
B.
,
Pereira
,
P.
,
De Baere
,
T.
,
Lencioni
,
R.
,
Leveillee
,
R.
,
Marchal
,
G.
,
Michel
,
L.
, and
Ni
,
Y.
,
2006
, “
Electrodes and Multiple Electrode Systems for Radio Frequency Ablation: A Proposal for Updated Terminology
,”
Adv. Exp. Med. Biol.
,
574
, pp.
57
73
.
9.
Singh
,
S.
, and
Repaka
,
R.
,
2017
, “
Temperature-Controlled Radiofrequency Ablation of Different Tissues Using Two-Compartment Models
,”
Int. J. Hyperthermia
,
33
(
2
), pp.
122
134
.
10.
Huo
,
Y. R.
,
Pillai
,
K.
,
Akhter
,
J.
, and
Morris
,
D. L.
,
2016
, “‘
Edgeboost’: A Novel Technique to Extend the Ablation Zone Lateral to a Two-Probe Bipolar Radiofrequency Device
,”
Cardiovasc. Interventional Radiol.
,
39
(
1
), pp.
97
105
.
11.
Lau
,
L.
, and
Han
,
Y.-L.
,
2016
, “
Exploring a Novel Heating Probe Design for Tumor Ablation
,”
ASME J. Med. Devices
,
10
(
3
), p.
030930
.
12.
Wang
,
Z.
,
Luo
,
H.
,
Coleman
,
S.
, and
Cuschieri
,
A.
,
2016
, “
Bicomponent Conformal Electrode for Radiofrequency Sequential Ablation and Circumferential Separation of Large Tumors in Solid Organs: Development and In Vivo Evaluation
,”
IEEE Trans. Biomed. Eng.
,
64
(
3
), pp.
699
705
.
13.
de Baere
,
T.
,
Denys
,
A.
,
Wood
,
B. J.
,
Lassau
,
N.
,
Kardache
,
M.
,
Vilgrain
,
V.
,
Menu
,
Y.
, and
Roche
,
A.
,
2001
, “
Radiofrequency Liver Ablation: Experimental Comparative Study of Water-Cooled Versus Expandable Systems
,”
AJR Am. J.·Roentgenol.
,
176
(
1
), pp.
187
192
.
14.
Jamil
,
M.
, and
Ng
,
E. Y.
,
2015
, “
Quantification of the Effect of Electrical and Thermal Parameters on Radiofrequency Ablation for Concentric Tumour Model of Different Sizes
,”
J. Therm. Biol.
,
51
, pp.
23
32
.
15.
Trujillo
,
M.
, and
Berjano
,
E.
,
2013
, “
Review of the Mathematical Functions Used to Model the Temperature Dependence of Electrical and Thermal Conductivities of Biological Tissue in Radiofrequency Ablation
,”
Int. J. Hyperthermia
,
29
(
6
), pp.
590
597
.
16.
Zhang
,
B.
,
Moser
,
M. A.
,
Zhang
,
E. M.
,
Luo
,
Y.
, and
Zhang
,
W.
,
2015
, “
Numerical Analysis of the Relationship Between the Area of Target Tissue Necrosis and the Size of Target Tissue in Liver Tumours With Pulsed Radiofrequency Ablation
,”
Int. J. Hyperthermia
,
31
(
7
), pp.
715
725
.
17.
Dodde
,
R. E.
,
Miller
,
S. F.
,
Geiger
,
J. D.
, and
Shih
,
A. J.
,
2008
, “
Thermal-Electric Finite Element Analysis and Experimental Validation of Bipolar Electrosurgical Cautery
,”
ASME J. Manuf. Sci. Eng.
,
130
(
2
), p.
021015
.
18.
Duan
,
B.
, and
Chui
,
C.-K.
,
2016
,
Multiscale Modeling of Liver Bio-Impedance and Frequency Control for Radiofrequency Ablation
,”
IEEE Region 10 Conference (
TENCON
), Singapore, Nov. 22–25, pp.
1532
1535
.
19.
Henriques
,
F. C.
, Jr.
,
1947
, “
Studies of Thermal Injury: Automatic Recording Caloric Applicator and Skin-Tissue and Skin-Surface Thermocouples
,”
Rev. Sci. Instrum.
,
18
(
9
), pp.
673
680
.
20.
Reddy
,
G.
,
Dreher
,
M. R.
,
Rossmann
,
C.
,
Wood
,
B. J.
, and
Haemmerich
,
D.
,
2013
, “
Cytotoxicity of Hepatocellular Carcinoma Cells to Hyperthermic and Ablative Temperature Exposures: In Vitro Studies and Mathematical Modelling
,”
Int. J. Hyperthermia
,
29
(
4
), pp.
318
323
.
21.
Chang
,
I. A.
, and
Nguyen
,
U. D.
,
2004
, “
Thermal Modeling of Lesion Growth With Radiofrequency Ablation Devices
,”
Biomed. Eng. Online
,
3
(
1
), p.
27
.
22.
Hall
,
S. K.
,
Ooi
,
E. H.
, and
Payne
,
S. J.
,
2015
, “
Cell Death, Perfusion and Electrical Parameters Are Critical in Models of Hepatic Radiofrequency Ablation
,”
Int. J. Hyperthermia
,
31
(
5
), pp.
538
550
.
23.
Wang
,
H.
,
He
,
Y.
,
Yan
,
Q.
,
You
,
F.
,
Fu
,
F.
,
Dong
,
X.
,
Shi
,
X.
, and
Yang
,
M.
,
2015
, “
Correlation Between the Dielectric Properties and Biological Activities of Human Ex Vivo Hepatic Tissue
,”
Phys. Med. Biol.
,
60
(
6
), pp.
2603
2617
.
24.
Macchi
,
E. G.
,
Gallati
,
M.
,
Braschi
,
G.
, and
Persi
,
E.
,
2014
, “
Dielectric Properties of RF Heated Ex Vivo Porcine Liver Tissue at 480 kHz: Measurements and Simulations
,”
J. Phys. D: Appl. Phys.
,
47
(
48
), p.
485401
.
25.
Gabriel
,
C.
,
Peyman
,
A.
, and
Grant
,
E.
,
2009
, “
Electrical Conductivity of Tissue at Frequencies Below 1 MHz
,”
Phys. Med. Biol.
,
54
(
16
), pp.
4863–4878
.
26.
Surowiec
,
A.
,
Stuchly
,
S.
, and
Swarup
,
A.
,
1985
, “
Radiofrequency Dielectric Properties of Animal Tissues as a Function of Time following Death
,”
Phys. Med. Biol.
,
30
(
10
), pp.
1131–1141
.
27.
Schutt
,
D. J.
, and
Haemmerich
,
D.
,
2008
, “
Effects of Variation in Perfusion Rates and of Perfusion Models in Computational Models of Radio Frequency Tumor Ablation
,”
Med. Phys.
,
35
(
8
), pp.
3462
3470
.
28.
Masuda
,
Y.
,
Foruzan
,
A. H.
,
Tateyama
,
T.
, and
Chen
,
Y. W.
,
2010
,
Automatic Liver Tumor Detection Using EM/MPM Algorithm and Shape Information
,”
Second IEEE International Conference on Software Engineering and Data Mining
(
SEDM
), Chengdu, China, June 23–25, pp.
692
695
.http://ieeexplore.ieee.org/abstract/document/5542834/
29.
Karampatzakis
,
A.
,
Kuhn
,
S.
,
Tsanidis
,
G.
,
Neufeld
,
E.
,
Samaras
,
T.
, and
Kuster
,
N.
,
2013
, “
Antenna Design and Tissue Parameters Considerations for an Improved Modelling of Microwave Ablation in the Liver
,”
Phys. Med. Biol.
,
58
(
10
), pp.
3191
3206
.