Antenna design and it is interactions with the liver tissues were simulated with Computer Simulation Technology (CST) software. CST [16] is a powerful simulation platform for electromagnetic problems. It provides multiple simulation modules, two of them were used in this work; the CST MICROWAVE STUDO and the CST MPHYSICS STUDO. The first one was used to calculate the antenna parameters and the specific absorption rate (SAR) in the liver and tumor, while the second module was used to calculate the temperature elevation due to the absorption of electromagnetic energy.
A. Antenna Design
The antenna was designed from a 50Ω coaxial cable, with copper inner and outer conductors separated by Polytetrafluoroethylene (commonly known by its trade name, Teflon) as a dielectric material. The antenna is a 3.4 mm in diameter with Teflon to act as a catheter. The antenna also covered by Teflon tape for easy insertion and removal after hyperthermia. Antenna geometry parameters, slot spacing, and floating sleeve length were chosen based on the effective wavelength in human liver tissue at 2.45 GHz, which calculated using the Eq. 1 [8]:
$${{\lambda }}_{\text{e}\text{f}\text{f}}= \frac{\text{c}}{\text{f}\sqrt{{{\epsilon }}_{\text{r}}} }$$
1
Where c is the speed of light in free space measured in (m/s), f is the operating frequency of the microwave generator (2.45 GHz), and \({{\epsilon }}_{\text{r}}\) = 43.03 is the relative permittivity of human liver tissue at the operating frequency; this yielded the effective wavelength of approximately 18.667 mm. The equation only provides a very crude approximation for the design. Generally, slot spacing, floating sleeve length corresponds to 0.25\({{\lambda }}_{\text{e}\text{f}\text{f}}\), and \({{\lambda }}_{\text{e}\text{f}\text{f}}\) respectively, which are chosen and then optimized to achieve localized power deposition near the distal tip of the antenna. Figure 4 shows the schematic diagram of the designed floating sleeve antenna (dimensions in mm). The diameters of each layer are illustrated in Table 1.
Table.1. Dimensions of the designed antenna
Dimension
|
Value (mm)
|
Diameter of the inner conductor
|
0.5
|
Outer Diameter of the dielectric material
|
1.6260
|
Outer diameter of outer conductor
|
1.7886
|
Inner diameter of the sleeve
|
2.5
|
Outer diameter of the sleeve
|
3.2
|
Diameter of the inner tip
|
2.6
|
Diameter of the outer tip
|
3.4
|
In the antenna design, the pattern of a linear dipole with an overall length less than the half-wavelength (L < λ/2) is insensitive to the frequency [17]. Thus, the length is chosen to be greater than the half-wavelength. The equation that used to calculate the wavelength is:
$${\lambda }=\frac{\text{c}}{\text{f}}$$
2
Where c is the speed of light in free space and f is the operating frequency of the microwave generator (2.45 GHz), this result in half-wavelength of 61.18mm.
A. The Original Liver Model
A realistic 3D liver model was downloaded from 3D-IRCADb-01 database [18], which composed of 3D Computed Tomography (CT) scans of 10 women and 10 men with hepatic tumors. In this study, a man liver with blood vessels and a tumor (2\(\times\) 3cm) was used. This model was used to generate 6 different models as will be explained later.
B. Assigning Material Properties
The electrical properties (dielectric constant and conductivity) and the thermal properties for different liver tissues were assigned according to the literature [19] [8, 9, 12, 20–27], and the CST material library. All dielectric properties were chosen at the operating frequency (2.45 GHz).
The assigning of metabolic properties was one of the main challenges in this work, since up to our knowledge; there are no any reported values for the metabolic rate of hepatic tumors. The literature cites the changing of tumors metabolic activities in comparison to the normal cells, and these changes support the acquisition and maintenance of malignant properties [28]. Gorbach et al. used the infrared camera to obtain thermal images of brain tumors, the obtained temperature distribution showed that glial origins of brain tumors have 0.5 to 2 \(\text{℃}\) temperature difference compared to surrounding brain tissues which is due to high metabolic processes [29]. In addition, a study of Mital et al. discussed that local temperatures of the skin over a tumor were significantly higher (about 2–3 degrees) than normal skin temperatures. This was due to convection effects associated with increased metabolism around the tumor [21].
In compare to breast and brain, liver has a higher metabolic rate [30]. Based on this, the basal metabolic rate in tumor was chosen to be 10 times of the surrounding tissues (120000 W/m3), with 3400 W/\({\text{m}}^{3}\)/K as perfusion rate. Table 2 presents the assigned material properties.
Table 2
Material properties used Bio-heat equation (at 2.45 GHz)
Property
|
Liver
|
Blood
|
Tumor
|
Relative permittivity (Ɛr)
|
43.03
|
58.30
|
48.16
|
Electrical conductivity (σ)
|
1.69
|
2.54
|
2.09
|
Density (\(\rho\))
|
1060
|
1000
|
1045
|
Specific heat capacity (cb)
|
3600
|
4180
|
4200
|
Thermal conductivity (k)
|
0.512
|
0.49
|
0.60
|
Blood flow coefficient
|
68000
|
1e + 006
|
3400
|
Basal metabolic rate BMR
|
12000
|
-
|
120000
|
Increasing of biological tissues temperature depends on the spatial distribution of the electromagnetic fields, the governing thermodynamics, and the thermal constitutive parameters of the biological system. The thermal energy calculations in this study were based on Pennes’ suggestion, expressed in the following equation:
$${\rho }\text{c}\frac{\partial \text{T}}{\partial \text{t}}=\text{k}\frac{{\partial }^{2}\text{T}}{\partial {\text{x}}^{2}}+\text{k}\frac{{\partial }^{2}\text{T}}{\partial {\text{y}}^{2}}+ \text{k}\frac{{\partial }^{2}\text{T}}{\partial {\text{z}}^{2}}+ {\text{w}}_{\text{b}}{\text{c}}_{\text{b}}\left({\text{T}}_{\text{a}}-\text{T}\right)+{\text{Q}}_{\text{m}}+{\text{Q}}_{\text{r}}(\text{x},\text{y},\text{z},\text{t})$$
3
Where \(\text{T}=\text{T}(\text{x},\text{y},\text{z},\text{t})\) is the temperature elevation (\(\text{℃}\)), ρ the physical density of the tissue (\(\text{k}\text{g}/{\text{m}}^{3}\)), c the specific heat of the tissue (\(\text{J}/\text{k}\text{g}/\text{℃}\)), k the tissue thermal conductivity (\(\text{W}/\text{m}\text{℃}\)), \({\text{w}}_{\text{b}}\) the blood volumetric perfusion rate (\({\text{k}\text{g}/\text{m}}^{3}/\text{s}\)), \({\text{c}}_{\text{b}}\) the specific heat of blood (\(\text{J}/\text{k}\text{g}/\text{℃}\)), and \({{\text{T}}_{\text{a}}=\text{T}}_{\text{a}}(\text{x},\text{y},\text{z},\text{t})\) the average temperature elevation of the arteries (\(\text{℃}\)). \({\text{Q}}_{\text{m}}\) is the mechanism for modeling physiological heat generation (\(\text{W}/{\text{m}}^{3}\)) and \({\text{Q}}_{\text{r}}\) the regional heat delivered by the source (\(\text{W}/{\text{m}}^{3}\)). The term\({\text{w}}_{\text{b}}{\text{c}}_{\text{b}}\left({\text{T}}_{\text{a}}-\text{T}\right)\), which is the perfusion heat loss (\(\text{W}/{\text{m}}^{3}\)), is always considered in the case of tissues which have a high degree of perfusion, such as a liver. In general, \({\text{w}}_{\text{b}}\) is assumed to be uniform throughout the tissue. However, its value may increase with heating time because of vasodilation and capillary recruitment [2].
C. The Whole System (Antenna And Liver) Simulation
The simulation started in CST MICROWAVE STUDO, where the antenna was designed. Then, the liver model was imported and materials properties were assigned. The transient solver of CST MICROWAVE STUDO was used to calculate the antenna parameters and the Specific Absorption Rate (SAR). The thermal solver in CST MPHYSICS STUDIO was used to solve the Pennes’ bioheat equation taking into account the effects of living tissue such as metabolic mechanisms and blood flow influence.
The initial temperature of 25\(\text{℃}\) and 37\(\text{℃}\) for the surrounding ambient and tissues respectively was set. While, open boundary conditions were set for both solvers. The results were taken with 10 minutes simulation time.
Figure 5 shows the whole system with different models. The designed antenna was tested firstly in a liver model without a tumor or blood vessels (Model A), then in a model with a tumor of about 2\(\times\)3cm and without any blood vessels (Model B). After that, the antenna is tested in a complete model with the same tumor and blood vessels (Model C). To test the effect of both size and location of the tumor, the location of the original tumor was changed (Model D) and the tumor was shrieked to (1.5\(\times\)1.5cm) (Model E). Finally, a model with a spherical tumor of 1.5 cm diameter (Model F) was tested. An Electromagnetic-Thermal coupling platform was initiated and in each time, SAR and temperature distributions were recorded.