Antenna design for microwave hepatic ablation using an axisymmetric electromagnetic model
© Bertram et al; licensee BioMed Central Ltd. 2006
Received: 21 October 2005
Accepted: 27 February 2006
Published: 27 February 2006
An axisymmetric finite element method (FEM) model was employed to demonstrate important techniques used in the design of antennas for hepatic microwave ablation (MWA). To effectively treat deep-seated hepatic tumors, these antennas should produce a highly localized specific absorption rate (SAR) pattern and be efficient radiators at approved generator frequencies.
Methods and results
As an example, a double slot choked antenna for hepatic MWA was designed and implemented using FEMLAB™ 3.0.
This paper emphasizes the importance of factors that can affect simulation accuracy, which include boundary conditions, the dielectric properties of liver tissue, and mesh resolution.
Liver cancer is a significant worldwide public health issue. The disease has a mortality rate of 100% at 5 years in untreated cases  and results in the deaths of more than one million people each year worldwide [2–5]. Although liver cancer can be treated successfully by surgical resection of the malignant tissue, approximately 90% of patients with the disease are ineligible for the procedure due to factors such as insufficient hepatic reserve and the close proximity of tumors to blood vessels [1, 6]. One promising alternative for these patients is hepatic microwave ablation (MWA), an experimental procedure in which an antenna is inserted percutaneously or during surgery  to induce cell necrosis through the heating of deep-seated tumors. Unlike other alternatives to resection such as RF ablation and cryoablation, MWA systems are able to produce large lesions in the presence of blood perfusion and are not restricted by tissue charring [5, 7]. MWA also has favorable one, two, and three year survival rates of 96%, 83%, and 73% .
The many perceived advantages of microwave ablation have driven researchers to develop innovative antennas to effectively treat deep-seated, nonresectable hepatic tumors. These designs have focused largely on thin, coaxial-based interstitial antennas [9–11], which are minimally invasive and capable of delivering a large amount of electromagnetic power. These antennas can usually be classified as one of three types (dipole, slot, or monopole) based on their physical features and radiative properties.
To assist in antenna design for MWA, many researchers have employed the use of mathematical models rooted in computational electromagnetics (CEM), a discipline that employs numerical methods to describe propagation of electromagnetic waves. These methods center around the formulation of discrete solutions to the fundamental electromagnetic equations collectively referred to as Maxwell's equations. Three principal techniques exist within CEM. The first of these, the finite-difference time-domain (FDTD) method, is based on the Yee algorithm  and uses finite difference approximations of the time and space derivatives of Maxwell's curl equations to create a discrete three-dimensional representation of the electric and magnetic fields. This method has been widely used to numerically evaluate the electromagnetic radiation patterns of antennas in tissue [13–17], although long computation times are generally required. Another commonly used CEM technique is method of moments (MoM) [18, 19], in which approximate numerical solutions to integral equations are formulated in the frequency domain to determine an unknown current distribution for an antenna. This distribution can then be subsequently extended to yield the antenna's radiation pattern.
An alternative to the two techniques above is the finite element method (FEM), which has been extensively used in simulations of cardiac and hepatic radio-frequency (RF) ablation [20, 21]. FEM models can provide users with quick, accurate solutions to multiple systems of differential equations and as such, are well suited to heat transfer problems like ablation. In this study, we have adapted an existing axisymmetric electromagnetic FEM model that was implemented using FEMLAB™ 3.0  to demonstrate proper techniques for the design of an antenna for hepatic MWA.
The FEM model used in this study was adapted from a coaxial slot antenna general model, developed by COMSOL for microwave cancer therapy . In this model, the electric and magnetic fields associated with the time-varying transverse electromagnetic (TEM) wave generated by the microwave source propagating in a coaxial cable in the z-direction was expressed in 2D axially symmetric cylindrical coordinates as
where r o and r 1are the outer and inner radii of the coaxial cable (m), P in is the input power (W), ε rd is the relative permittivity of the dielectric, is the wave impedance in the dielectric of the coaxial cable (Ω), is the intrinsic impedance (Ω), ε 0 = 8.854 × 10-12 is the permittivity of free space (F/m), μ 0 = 4π × 10-7 is the permeability of free space (H/m), ω = 2π f is the angular frequency (rad/s), f is the frequency (Hz), k = 2π/λ is the propagation constant (m-1), and λ is the wavelength (m). For interstitial coaxial-based antennas during MWA, the magnetic field is purely azimuthal. The electric field is in the radial direction only inside the coaxial cable and in both radial and the axial direction inside the tissue. This allows for the antenna to be modeled using an axisymmetric transverse magnetic (TM) wave formulation, in which the source was modeled as a low reflection boundary
with excitation magnetic field H φ 0= C/Zr. Such axisymmetric FEM models are highly desirable as they dramatically reduce computation time.
Careful examinations of both the antenna's specific absorption rate (SAR) pattern and frequency-dependent reflection coefficient in tissue are essential for the optimization of antennas for hepatic MWA. SAR represents the electromagnetic power deposited per unit mass in tissue (W/kg) and can be defined mathematically as
where σ is tissue conductivity (S/m) and ρ is tissue density (kg/m3) . For the treatment of deep-seated hepatic tumors, the SAR pattern of an interstitial antenna should be highly localized near the distal tip of the antenna. Antenna efficiency can be quantified using the frequency-dependent reflection coefficient, which can be expressed logarithmically as
where P r indicates reflected power (W). The frequency where the reflection coefficient is minimum is commonly referred to as the resonant frequency and should be approximately the same as the operating frequency of the generator used. Antennas operating with high reflection coefficients (especially at higher power levels) can cause overheating of the feedline possibly leading to damage to the coaxial line or due to the thin outer conductor damage to the tissue .
Materials and methods
Several antenna designs have been shown to be effective for the treatment of deep-seated hepatic tumors. Recently, Saito et al  presented a design for a double slot antenna that is capable of achieving a higher degree of SAR localization than a standard single slot antenna and has been used in initial clinical trials. Another popular design, the cap-choke antenna [29–31], was originally designed for cardiac MWA applications but has been also found through independent investigation to be an effective applicator for hepatic MWA. In this design, the choke acts as a balun (which isolates the antenna and feedline to provide a balanced output) to prevent current backflow along the axial length of the antenna, resulting in high SAR localization near the distal tip of the antenna [32–34].
where c is the speed of light in free space (m/s), f is the operating frequency of the microwave generator (2.45 GHz), and ε r = 44.4 is the relative permittivity of bovine liver tissue at the operating frequency; this yielded 18.4 mm for the effective wavelength. However, because the catheter and its thickness also affect the optimal geometry and performance of the antenna, equation 7 only provides a very crude approximation for the design. Fig. 2 shows that slot spacing, choke offset, and choke length correspond to 0.25λ eff, 0.5λ eff, and λ eff respectively, which were chosen to achieve localized power deposition near the distal tip of the antenna. Tip length and catheter thickness were adjusted to achieve resonance.
Frequency-independent material parameters.
Coaxial cable dielectric
To determine the frequency response of the antenna, the model was simulated in FEMLAB™ at multiple discrete frequencies between 0.8 and 10 GHz using 50 MHz increments and the dielectric properties for bovine liver tissue shown in Fig. 1. At each frequency, a boundary integration of power outflow along the source was performed using FEMLAB™ to determine the net power delivered to the antenna. Subsequently, (6) was used to calculate the reflected power P r and frequency dependent reflection coefficient Γ(f).
FEMLAB™ mesh parameters.
Element growth rate
Maximum element size scaling factor
Mesh curvature cut off
Mesh curvature factor
To validate the performance of the FEMLAB™ antenna model, a double slot choked antenna was constructed from UT-85 semirigid coaxial cable and gold-plated SMA connectors obtained from Microstock, Inc. This antenna was physically immersed in saline without a catheter and connected to an Agilent E8364A performance network analyzer (PNA) to measure the frequency-dependent reflection coefficient of the antenna. The dielectric properties of saline were measured using the same dielectric spectroscopy procedure mentioned above for liver tissue. These properties were subsequently used in a separate FEMLAB™ model to simulate the reflection coefficient of the antenna in saline without a catheter and were compared against measured results. A catheter was not used for the antenna in both measurement and FEMLAB model in order to avoid the possible effects by the catheter. Antenna performance is quite sensitive to the catheter thickness, but for practical reasons, the thickness could not be controlled exactly during the course of antenna construction. The comparison with the antenna without a catheter would still well serve the purpose to validate the computer model to experimental measurements.
Fig. 3(b) shows that the double slot choked antenna is able to provide a degree of power localization similar to those produced by other recently designed antennas. Volumes for 50% and 90% energy dissipation are 7.55 cm3 and 99.3 cm3 respectively. This initial analysis also shows that it may be capable of providing more localized power deposition than the earlier double slot antenna and standard cap-choke antenna, and should therefore be considered as a viable design for future research. Fig. 4 shows that the antenna is resonant at 2.45 GHz with a low reflection coefficient of approximately -22 dB. Reducing catheter thickness by 0.1 mm slightly improved localization but decreased the resonant frequency of the antenna by approximately 400 MHz. Such effects were expected because the catheter is an important part of the antenna structure. However more studies are apparently necessary to learn the exact effects of the catheter and the antenna geometry with catheter should be optimized together for the desired resonant frequency and SAR localization.
Several factors should be considered when designing and simulating antennas for hepatic MWA. Due to phase considerations, both antenna dimensions and the thickness of the catheter should be chosen and optimized based on the effective wavelength in tissue. As this quantity is highly dependent on the relative permittivity of the medium, it is also essential that electromagnetic models of antennas for hepatic MWA employ accurate dielectric properties of liver tissue. Great care should be taken when measuring these properties or using previously published data, as the dielectric properties of liver tissue are highly dispersive and change with tissue type , water content , and temperature .
In addition, tradeoffs that are not well understood exist between performance metrics such as SAR and the reflection coefficient. From our experience, it is usually more desirable to achieve a localized power deposition than a highly resonant antenna since external impedance matching networks can be used to improve antenna efficiency . However, our additional experiments have shown that due to thermal conduction, a localized SAR pattern will only result in controlled tissue heating when treatment duration is minimized. Not only through the tissue, heat can also be effectively conducted by the antenna metal body and result in uncontrolled tissue heating along the antenna body for long treatment durations even with localized SAR pattern of the antenna.
Compared with other CEM software packages such as XFDTD™ and FEKO™, FEMLAB™ provides a user-friendly graphical interface for quick and reliable antenna design using axisymmetric electromagnetic models that can be integrated with MATLAB™ for increased functionality. Also, as the software was designed for multiphysics applications, it should also be ideal for future thermal simulations of hepatic MWA that are based on the Pennes bio-heat equation  and are dependent on both the electromagnetic and thermal properties of tissue. Initial work has also shown that electromagnetic models implemented using FEMLAB™ can be used in genetic algorithms  to better optimize antennas.
An antenna for hepatic MWA was quickly and accurately simulated using an axisymmetric electromagnetic model implemented in FEMLAB™ 3.0. This model is ideal for analyzing the SAR patterns and reflection coefficient of antennas used in hepatic MWA, two metrics which are commonly used for the evaluation of such antennas. Using this model, a double slot choked antenna with performance comparable to current antennas was designed and evaluated.
The authors would like to thank Mariya Lazebnik and Prof. Susan Hagness at the University of Wisconsin-Madison for their assistance with dielectric spectroscopy measurements and validation of this antenna. This work was supported by National Institute of Health (NIH) grant DK58839.
- McGahan JP, Dodd GD: Radiofrequency ablation of the liver: Current status. Am J Roentgenol 2001, 176: 3–16.View ArticleGoogle Scholar
- Chinn SB, Lee FT, Kennedy GD, Chinn C, Johnson CD, Winter TC, Warner TF, Mahvi DM: Effect of vascular occlusion on radiofrequency ablation of the liver: Results in a porcine model. Am J Roentgenol 2001, 176: 789–795.View ArticleGoogle Scholar
- Haemmerich D, Staelin ST, Tungjitkusolmun S, Lee FT, Mahvi DM, Webster JG: Hepatic bipolar radio-frequency ablation between separated multiprong electrodes. IEEE Trans Biomed Eng 2001, 48: 1145–1152. 10.1109/10.951517View ArticleGoogle Scholar
- Lau WY, Leung TWT, Yu SCH, Ho SKW: Percutaneous local ablative therapy for hepatocellular carcinoma - A review and look into the future. Ann Surg 2003, 237: 171–179. 10.1097/00000658-200302000-00005Google Scholar
- Wright AS, Mahvi DM, Haemmerich DG, Lee FTJ: Minimally invasive approaches in management of hepatic tumors. Surg Technol Int 2003, 11: 144–153.Google Scholar
- Cha C, Lee FT: Rationale for the combination of cryoablation with surgical resection of hepatic tumors. Journal of Gastrointestinal Surgery 2001, 5: 206–213. 10.1016/S1091-255X(01)80034-2View ArticleGoogle Scholar
- Wright AS, Lee FT, Mahvi DM: Hepatic microwave ablation with multiple antennae results in synergistically larger zones of coagulation necrosis. Ann Surg Oncol 2003, 10: 275–283. 10.1245/ASO.2003.03.045View ArticleGoogle Scholar
- Lu MD, Chen JW, Xie XY, Liu L, Huang XQ, Liang LJ, Huang JF: Hepatocellular carcinoma: US-guided percutaneous microwave coagulation therapy. Radiology 2001, 221: 167–172.View ArticleGoogle Scholar
- Schaller G, Erb J, Engelbrecht R: Field simulation of dipole antennas for interstitial microwave hyperthermia. IEEE Trans Microwave Theory Tech 1996, 44: 887–895. 10.1109/22.506448View ArticleGoogle Scholar
- Longo I, Gentili GB, Cerretelli M, Tosoratti N: A coaxial antenna with miniaturized choke for minimally invasive interstitial heating. IEEE Trans Biomed Eng 2003, 50: 82–88. 10.1109/TBME.2002.807320View ArticleGoogle Scholar
- Brace CL, van der Weide DW, Lee FT, Laeseke PF, Sampson LA: Analysis and experimental validation of a triaxial antenna for microwave tumor ablation. Microw Symp Dig 2004 IEEE MTT-S Int 2004, 3: 1437–1440.View ArticleGoogle Scholar
- Yee KS: Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media. IEEE Trans Ant Prop 1966, AP-14: 302–307.Google Scholar
- Sullivan D: Three-dimensional computer simulation in deep regional hyperthermia using the finite-difference time-domain method. IEEE Trans Microw Theory Tech 1990, 38: 204–211. 10.1109/22.46431View ArticleGoogle Scholar
- Chen JY: Numerical simulation of annular-phased arrays of dipoles for hyperthermia of deep-seated tumors. IEEE Trans Biomed Eng 1992, 39: 209–216. 10.1109/10.108122View ArticleGoogle Scholar
- Cherry PC, Iskander MF: Calculations of heating patterns of an array of microwave interstitial antennas. IEEE Trans Biomed Eng 1993, 40: 771–779. 10.1109/10.238461View ArticleGoogle Scholar
- Gentili GB, Leoncini M, Trembly BS, Schweitzer SE: FDTD electromagnetic and thermal analysis of interstitial hyperthermic applicators. IEEE Trans Biomed Eng 1995, 42: 973–980. 10.1109/10.464371View ArticleGoogle Scholar
- Camart JC, Despretz D, Chive M, Pribetich J: Modeling of various kinds of applicators used for microwave hyperthermia based on the FDTD method. IEEE Trans Microw Theory Tech 1996, 44: 1811–1818. 10.1109/22.539939View ArticleGoogle Scholar
- Cerri G, Deleo R, Primiani VM: Thermal end-fire interstitial applicator for microwave hyperthermia. Ieee T Microw Theory 1993, 41: 1135–1142. 10.1109/22.238538View ArticleGoogle Scholar
- Sabariego RV, Landsea L, Obelleiro F: Synthesis of an array antenna for hyperthermia applications. IEEE Trans Magn 2000, 36: 1696–1699. 10.1109/20.877769View 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
- Haemmerich D, Chachati L, Wright AS, Mahvi DM, Lee FT, Webster JG: Hepatic radiofrequency ablation with internally cooled probes: Effect of coolant temperature on lesion size. IEEE Trans Biomed Eng 2003, 50: 493–500. 10.1109/TBME.2003.809488View ArticleGoogle Scholar
- Electromagnetics Module Model Library In FEMLAB 30. Stockholm, Sweden:, COMSOL AB; 2004:32.
- Clibbon KL, Mccowen A: Efficient computation of SAR distributions from interstitial microwave antenna arrays. IEEE Trans Microw Theory Tech 1994, 42: 595–600. 10.1109/22.285065View ArticleGoogle Scholar
- Nevels RD, Arndt GD, Raffoul GW, Carl JR, Pacifico A: Microwave catheter design. IEEE Trans Biomed Eng 1998, 45: 885–890. 10.1109/10.686796View ArticleGoogle Scholar
- Popovic D, McCartney L, Beasley C, Lazebnik M, Okoniewski M, Hagness SC, Booske J: Precision open-ended coaxial probes for in-vivo and ex-vivo dielectric spectroscopy of biological tissues at microwave frequencies. IEEE Trans Microw Theory Tech 2004.Google Scholar
- Anderson J, Sibbald C, Stuchly S: Dielectric measurements using a rational function model. IEEE Trans Microwave Theory Tech 1994, 42: 199–204. 10.1109/22.275247View ArticleGoogle Scholar
- Stuchly SS, Sibbald CL, Anderson JM: A new aperture admittance model for open-ended waveguides. IEEE Trans Microw Theory Tech 1994, 42: 192–198. 10.1109/22.275246View ArticleGoogle Scholar
- Saito K, Yoshimura H, Ito K, Aoyagi Y, Horita H: Clinical trials of interstitial microwave hyperthermia by use of coaxial-slot antenna with two slots. IEEE Trans Microw Theory Tech 2004, 52: 1987–1991. 10.1109/TMTT.2004.832005View ArticleGoogle Scholar
- Lin JC, Wang YJ: The cap-choke catheter antenna for microwave ablation treatment. IEEE Trans Biomed Eng 1996, 43: 657–660. 10.1109/10.495286View ArticleGoogle Scholar
- Pisa S, Cavagnaro M, Bernardi P, Lin JC: A 915-MHz antenna for microwave thermal ablation treatment: physical design, computer modeling and experimental measurement. IEEE Trans Biomed Eng 2001, 48: 599–601. 10.1109/10.918599View ArticleGoogle Scholar
- Bernardi P, Cavagnaro M, Lin JC, Pisa S, Piuzzi E: Distribution of SAR and temperature elevation induced in a phantom by a microwave cardiac ablation catheter. IEEE Trans Microw Theory Tech 2004, 52: 1978–1986. 10.1109/TMTT.2004.832031View ArticleGoogle Scholar
- Prasad S, King RWP: Experimental study of embedded insulated antennas. IEEE Trans Ant Prop 1982, 30: 1013–1017. 10.1109/TAP.1982.1142896View ArticleGoogle Scholar
- Chiu HM, Mohan AS, Weily AR, Guy DJ, Ross DL: Analysis of a novel expanded tip wire (ETW) antenna for microwave ablation of cardiac arrhythmias. IEEE Trans Biomed Eng 2003, 50: 890–899. 10.1109/TBME.2003.813541View ArticleGoogle Scholar
- Icheln C, Krogerus J, Vainikainen P: Use of balun chokes in small-antenna radiation measurements. IEEE Trans Instrum Meas 2004, 53: 498–506. 10.1109/TIM.2004.823299View ArticleGoogle Scholar
- Stauffer PR, Rossetto F, Prakash M, Neuman DG, Lee T: Phantom and animal tissues for modelling the electrical properties of human liver. Int J Hyperthermia 2003, 19: 89–101. 10.1080/0265673021000017064View ArticleGoogle Scholar
- Foster KR, Schwan HP: Dielectric-properties of tissues and biological-materials - a critical-review. Crit Rev Biomed Eng 1989, 17: 25–104.Google Scholar
- Chin L, Sherar M: Changes in dielectric properties of ex vivo bovine liver at 915 MHz during heating. Phys Med Biol 2001, 46: 197–211. 10.1088/0031-9155/46/1/314View ArticleGoogle Scholar
- Jones KM, Mechling JA, Strohbehn JW, Trembly BS: Theoretical and experimental SAR distributions for interstitial dipole antenna-arrays used in hyperthermia. Ieee T Microw Theory 1989, 37: 1200–1209. 10.1109/22.31079View ArticleGoogle Scholar
- Pennes HH: Analysis of tissue and arterial blood temperatures in the resting human forearm. J Appl Phys 1948, 1: 93–122.Google Scholar
- Johnson JM, Rhamat-Samii Y: Genetic algorithms in engineering electromagnetics. IEEE Ant Prop Mag 1997., 39: Google Scholar
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.