0
Research Papers

Design of a Novel Electrode of Radiofrequency Ablation for Large Tumors: A Finite Element Study OPEN ACCESS

[+] Author and Article Information
Zheng Fang

Tumor Ablation Group,
CISR Center,
East China University of
Science and Technology,
Shanghai 200237, China
e-mail: jasonfang2014@hotmail.com

Bing Zhang

Mem. ASME
Division of Biomedical Engineering,
University of Saskatchewan,
Saskatoon, SK S7N 5A9, Canada
e-mail: bing.zhang84@usask.ca

Michael Moser

Department of Surgery,
University of Saskatchewan,
Saskatoon, SK S7N 0W8, Canada
e-mail: mam3045@mail.usask.ca

Edwin Zhang

Division of Vascular and
Interventional Radiology,
Department of Medical Imaging,
University of Toronto,
Toronto, ON M5T 1W7, Canada
e-mail: edwinmzhang@gmail.com

Wenjun Zhang

Fellow ASME
Department of Mechanical Engineering,
University of Saskatchewan,
Saskatoon, SK S7N 5A9, Canada
e-mail: chris.zhang@usask.ca

1Corresponding author.

Manuscript received July 12, 2017; final manuscript received September 29, 2017; published online October 19, 2017. Assoc. Editor: Linxia Gu.

ASME J of Medical Diagnostics 1(1), 011001 (Oct 19, 2017) (6 pages) Paper No: JESMDT-17-2017; doi: 10.1115/1.4038129 History: Received July 12, 2017; Revised September 29, 2017

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.

FIGURES IN THIS ARTICLE
<>

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.

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.

The governing equation of the RFA procedure follows the Penne’s bioheat transfer equation [1] Display Formula

(1)ρcT(x,t)t=·kT(x,t)+ρbcbωb(Tb-T(x,t))+Qhs(x,t)+Qm(x,t)x

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]: Display Formula

(2)Qhs=J·E=σU2

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]: Display Formula

(3)k(T)=k0+0.0013(Tx,t-T0)
Display Formula
(4)σ(T)=σ01+0.02(Tx,t-T0)

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] Display Formula

(5)ωbt=ωb0forT50°C0forT>50°C

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

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

(6)Ωt=A0te-ΔERT(τ)dτ

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.

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: Display Formula

(7)V=π(d1·d2·d3)/6

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.

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] Display Formula

(8)SI=V[π·maxd1,d2,d33]/6

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

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.

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.

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.

  • National Sciences and Engineering Research Council of Canada (NSERC) through a discovery grant to Wenjun Zhang.

  • National Natural Science of China (Grant No. 51375166).

  • 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

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. [CrossRef] [PubMed]
Ahmed, M. , and Goldberg, S. N. , 2012, “ Principles of Radiofrequency Ablation,” Interventional Oncology, Springer, New York, pp. 23–37. [CrossRef]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef]
Goldberg, S. N. , 2001, “ Radiofrequency Tumor Ablation: Principles and Techniques,” Eur. J. Ultrasound, 13(2), pp. 129–147. [CrossRef] [PubMed]
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. [CrossRef]
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. [CrossRef] [PubMed]
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. [CrossRef]
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. [CrossRef]
Lau, L. , and Han, Y.-L. , 2016, “ Exploring a Novel Heating Probe Design for Tumor Ablation,” ASME J. Med. Devices, 10(3), p. 030930. [CrossRef]
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. [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef]
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.
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
Chang, I. A. , and Nguyen, U. D. , 2004, “ Thermal Modeling of Lesion Growth With Radiofrequency Ablation Devices,” Biomed. Eng. Online, 3(1), p. 27. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef]
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. [CrossRef]
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. [CrossRef]
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. [CrossRef] [PubMed]
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/
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. [CrossRef] [PubMed]
Copyright © 2018 by ASME
View article in PDF format.

References

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. [CrossRef] [PubMed]
Ahmed, M. , and Goldberg, S. N. , 2012, “ Principles of Radiofrequency Ablation,” Interventional Oncology, Springer, New York, pp. 23–37. [CrossRef]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef]
Goldberg, S. N. , 2001, “ Radiofrequency Tumor Ablation: Principles and Techniques,” Eur. J. Ultrasound, 13(2), pp. 129–147. [CrossRef] [PubMed]
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. [CrossRef]
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. [CrossRef] [PubMed]
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. [CrossRef]
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. [CrossRef]
Lau, L. , and Han, Y.-L. , 2016, “ Exploring a Novel Heating Probe Design for Tumor Ablation,” ASME J. Med. Devices, 10(3), p. 030930. [CrossRef]
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. [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef]
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.
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
Chang, I. A. , and Nguyen, U. D. , 2004, “ Thermal Modeling of Lesion Growth With Radiofrequency Ablation Devices,” Biomed. Eng. Online, 3(1), p. 27. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef] [PubMed]
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. [CrossRef]
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. [CrossRef]
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. [CrossRef]
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. [CrossRef] [PubMed]
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/
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. [CrossRef] [PubMed]

Figures

Grahic Jump Location
Fig. 1

Dimensions (in mm) and structures of two RF electrodes: (a) CEE and (b) the designed electrode

Grahic Jump Location
Fig. 2

Three-dimensional geometric model (in mm) of liver tissue with the liver tumor and the designed RF electrode

Grahic Jump Location
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

Grahic Jump Location
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)

Grahic Jump Location
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)

Tables

Table Grahic Jump Location
Table 1 Physical properties of the materials used in the model
Table Footer NoteaEvaluated at 21 °C.
Table Grahic Jump Location
Table 2 Number of elements for each scenario
Table Grahic Jump Location
Table 3 Comparisons of ablation size from the FEM and ex vivo/in vivo experiments
Table Grahic Jump Location
Table 4 Ablation volumes and SIs of three electrodes
Table Grahic Jump Location
Table 5 Ablation results of CEE and MEE for a liver tumor with 3.5 cm in diameter

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In