ap-4-10.dvi Acta Polytechnica Vol. 50 No. 4/2010 Simulation of Hyperthermic Treatment Using the Matrix of Stripline Applicators B. Vrbová, L. Víšek Abstract This paper describes the design of a microwave stripline applicator for hyperthermic treatment, and the design of an anatomically based biological model, which is a necessary part of hyperthermia treatment planning for measuring the distribution of SAR. In this paper we compare the SAR distribution in a cylindrical homogeneous agar phantom (which has similar characteristics to biological tissue) and in an anatomically based biological model of the femur (which has been developed from a computer tomography scan) using a matrix of two applicators of the same type. Keywords: microwave thermotherapy, SAR, cancer treatment, anatomically based models, hyperthermia applicators. 1 Introduction Thermotherapy is a method based on differences in the behavior of healthy tissue and tumor tissue under enhanced temperatures. One of the most important methods ofmicrowavethermotherapy ishyperthermia. This method works in a temperature interval between 41 and 45◦C inwhich cancer cells are destroyed,while healthy cells are able to survive up to 45◦C. If the temperature is increased above 45◦C, coagulation of healthy cells occurs. The self protectivemechanism of tumors more than 2 cm in diameter already fails at a temperature of 41◦C. The blood flow in tumor cells decreaseswith increasing temperature, andso the tem- perature in tumor cells increases even more rapidly. Finally, the tumor tissue is destroyed. The duration of a single hyperthermic treatment should not exceed 50 minutes. The level of the hyperthermic dose de- pends on temperature and time [1]. Ahyperthermic set consists of ahighpowergenera- tor, thermal sensorsand, especially, theapplicator. An electromagnetic wave is generated and the applicator delivers it to the biological tissue, where it is absorbed, because the biological tissue is a lossy dielectric ma- terial. Applicators are mainly designed for a working frequency of 434 MHz. This is one of the ISM (Indus- trial, Scientific, Medical) dedicated frequencies [2]. The distribution of the electromagnetic field in the biological tissue can be calculated using a simple, ho- mogeneous model of tissue, which consists of only one type of tissue with defined electric parameters. How- ever, the real tissue is much more complicated. The electromagnetic field diffuses through several types of tissue, so simulationswith anatomicalmodels are used for more realistic and more accurate results. These models can also be used for hyperthermia treatment planning. Treatment planning is a very im- portant procedure in hyperthermia. Before the pa- tient is treated by hyperthermia, several calculations are performed to find out the distribution of the ab- sorbed power in this particular case. The best type of applicator for the treatment is selected according to the location of the tumor. In cases where the tu- mor is located near to vital organs or other critical areas, calculations with anatomicalmodels of patients are necessary. 2 Hyperthermia applicator The basic parameters of a TEM wave depend mainly on the dielectric parameters of the media through which it propagates (complex permittivity ε) and on the working frequency, but not on the cross section dimensions or on the type of transmission line [3]. The applicator discussed hereworks at a frequency of 434MHz, and ismade fromahighly conductivema- terial (copper). The sides of the applicator are made fromacrylic glass. TheTEMwave is transferredalong a section of the microwave stripline transmission line with cross-sectiondimensionsof 50×30mmand length 80mm. Thehorn section of the applicator is 80mm in length. The horn aperture is 120×80mm (Fig. 1) [4]. Fig. 1: Model of a stripline applicator 106 Acta Polytechnica Vol. 50 No. 4/2010 The length of the applicator is equal to twice the wave length, depending on the relative permittivity of the biological tissue (εr = 78). To avoid reflection of the wave back to the generator, the applicator is filled with a suitable dielectric material. In our case, the applicator is filled with distilled water, because 70 % of biological tissue consists of water. A water bolus is inserted between the applicator and the tissue to improve the transmission of the elec- tromagnetic energy into the tissue. This improves the thermal profile in the biological tissue. Impedance matching is very important for pre- venting reflections of electromagnetic energy back to the generator. Reflected energy could destroy the high power output generator. The basic condition for impedance matching is that VSWR has to be lower than 2. The impedance matching of the applicator was simulated using the FDTD simulator (e.g. SEM- CAD X EM Field simulator from SPEAG, Schmid & Partner Engineering AG, Switzerland [6]) in order to find the position and the length of the exciting probe. The reflection coefficient at working frequency s11 is equal to −32.3 dB, which represents a value of V SW R = 1.05. The length of the exciting probe is equal to 27 mm and the distance from the short- circuited plane is 15 mm. 3 Anatomically based biological model Ananatomically basedbiologicalmodel is essential for numerical dosimetry. The numerical model is usually developed from CT scans (Fig. 2). Fig. 2: Computed tomography scan Dosimetry is used in the design and application of hyperthermia applicators. Dosimetry quantifies the magnitude and distribution of the absorbed electro- magnetic energy within biological objects exposed to electromagnetic fields. InRF, the dosimetric quantity, referred to as the specific absorption rate (SAR), is defined as the rate at which energy is absorbed per unit mass. SAR is determined not only by the inci- dent electromagnetic waves but also by the electrical and geometric characteristics of the irradiated subject and nearby objects. It is related to the internal elec- tric field strength and also to the electric conductivity and the density of the tissues. It is therefore a suit- able dosimetric parameter, even when a mechanism is determined to be “athermal”. SAR distributions are usually determined from calculations on human mod- els. Rapid progress with computers has enabled high- level numericaldosimetry tobeperformedwith the aid of high-resolution biological models. The finite differ- ence time domain (FDTD) method is currently the mostwidely used numerical RFdosimetrymethod [5]. In order to develop a model for numerical dosime- try, the original gray-scale data must be interpreted into tissue types. This process is known as segmenta- tion. Segmentation of medical images involves parti- tioning the data into contiguous regions representing individual anatomical objects. This is a prerequisite for further investigations in many computer-assisted medical applications, e.g. individual therapy plan- ning and evaluation, diagnosis, simulation and image- guided surgery. Segmentation is adifficult taskbecause it is inmost cases very hard to separate the object from the image background. This is due to the characteristics of the imaging process as well as the grey-value mappings of the objects themselves. The most common medi- cal imageacquisitionmodalities includeCT(computer tomography) and MRI (magnetic resonance imaging) images. MRI or CT provides gray-scale image data in the form of many transverse slices, at a designated spacing, from the head to the feet of the biological body. The resolution in each slice is of the order of several millimeters. The gray-scale images are first rescaled to produce appropriate voxels. Each voxel in the images then is identified rigorously as belonging to one type of tissue. This is donebyassigningto eachvoxela red-green-blue code that identifies thediscrete tissue type of that par- ticularvoxel. Thisprocess canbeperformedwithcom- mercial software [7]. All identified transverse images are then combined to obtain a three-dimensional nu- merical model. Fine adjustment is generally required to connect each slice smoothly in the three orthogonal planes (axial, sagittal and coronal). DICOM CT scans for a femur model were ob- tained from Bulovka University Hospital. The femur model has resolution 2 mm, meaning a voxel size of 2 × 2 × 2 mm. Each voxel was assigned to one of 3 different tissue types, i.e. muscle, fat and bone (Fig. 3). 107 Acta Polytechnica Vol. 50 No. 4/2010 Fig. 3: Anatomical model of a femur 4 Matrix composition of two applicators When the tumours cover a large area of the human body, e.g. on the arm, leg, back or stomach area, we can use applicators of bigger dimensions, or we can create a matrix of several smaller applicators. In our case, we chose to use a matrix. We chose the matrix composition of two applicators of the same type. In the first case we put two applicators on a cylindrical agar phantom next to each other (Fig. 4). Fig. 4: Matrix of two applicators on a cylindrical agar phantom This homogeneous agar phantom represents a fe- mur from a human body, and its dielectric properties are shown in the following table (Tab. 1). The radius of the cylindrical agar phantom is 8 cm. In the second case, we placed a matrix of two ap- plicators on an anatomically-based biological model (Fig. 5), the dielectric parameters of which are listed in Table 1. Table 1: Dielectric properties at a frequencyof 434MHz [8] Name Conductivity [S/m] Relative permittivity Agar 0.80 54.00 Cortical Bone 0.09 13.07 Muscle 0.80 56.86 Fat 0.04 5.56 Fig. 5: Matrix composition of two applicators on an anatomically-based biological model 5 Results The results of the simulations of SAR distribution are shown in the following pictures. The maximum SAR distribution is situated in the middle of the two appli- cators in the homogeneous agar phantom (Fig. 6). Fig. 6: SAR distribution in the agar phantom in the lon- gitudinal cutting plane A comparison of the SAR distribution in the agar phantom and in the anatomic model shows that the SARdistributions acquire theirmaximumat the same locations, but the shape of the SARdistribution in the anatomical phantom is influenced by the fat and bone (Fig. 7). 108 Acta Polytechnica Vol. 50 No. 4/2010 Fig. 7: SAR distribution in the anatomical model in the longitudinal cutting plane Fathasa lowerpermittivityvalue thanmuscle, and therefore a small part of the energy is absorbed in the fat andmost of the energygoes into another layer, into the muscle. Muscle behaves as a lossy environment, so the energy is absorbed there. Bones also affect the shape of the SAR distribution, as can be seen in the transversal section of the femur (Fig. 8). Fig. 8: SAR distribution in the transversal cutting plane of the anatomic model Bones, like fat, have a lowpermittivity value (they contain a small amount of water), but the SAR value is almost equal to zero in this place, because bone be- haves as a lossless environment andmost of the energy goes through the bones. Fig. 9 shows the shape of the SAR distribution in the transversal layer of the homogeneous agar phan- tom. Simulations of SAR distribution shows that mi- crowave stripline applicators can be used to treat tu- mors locatedunder the surface of tissue. Simulation of the SAR distribution in an agar phantom is used for testing applicators. However, in cases of hyperther- mic treatment planning where the tumor is located near to vital organs or other critical areas, calcula- tions must be made with an anatomical model of the patient. Fig. 9: SAR distribution in the transversal cutting plane of an agar phantom 6 Conclusions Aset of severalapplicators canbeused inclinical prac- tice to treat tumors covering a large area of the human body. In future the treatment of tumors can be im- proved by setting the amplitude and the phase of the applicators. Acknowledgement This researchhasbeen supportedby theGrantAgency of theCzechRepublic project: “Non-standardapplica- tion of physical fields – analogy,modelling, verification and simulation” (102/08/H081). References [1] Falk, H. M., Issels, R. D.: Hyperthermia in Oncology. International Journal of Hyperthermia, Vol. 17, No. 1, 1–18, 2001. [2] Vrba, J.: Medical Applications of Microwaves. Prague (Czech Republic), 2007. ISBN 978-80-01-02705-9. [3] Vrba, J.: Introduction to Microwave Technology. Prague (Czech Republic), 2007. ISBN 978-80-01-03670-9. [4] Vrbová, B.: Diploma thesis. Microwave stripline applicator for local thermotherapy. Prague, 2009. [5] Fujiwara, O., Wang, J.: Electromagnetics in Biology, Springer Japan, chapter on Radiofre- quency Dosimetry and Exposure Systems, 2006, p. 223–225. [6] Schmid & Partner Engineering AG, [Online], URL: http://www.semcad.com/, 2009. [7] 3D-DOCTOR, 3D Imaging, Modelling, Rendering and Measurement Software, [Online], URL: http://www.ablesw.com/3d-doctor, 2009. 109 Acta Polytechnica Vol. 50 No. 4/2010 [8] Gabriel, C.: Compilation of the Dielectric Proper- ties ofBodyTissuesatRFandMicrowaveFrequen- cies, Brooks Air Force Technical Report, AL/OE- TR-1996-0037. About the authors Barbora VRBOVÁ was born in Nové Zámky, Slo- vakia, onMarch28, 1984. She receivedherMScdegree in biomedical engineering from the Czech Technical University in Prague in 2009. She deals with radio- metric methods for verifying the biological effects of EM fields. Lukáš VÍŠEKwas born inVysokéMýto in 1982 and receivedhisMScdegree fromtheCzechTechnicalUni- versity in Prague in February 2006. He is currently a postgraduate student at the Department of Electro- magnetic Field at the Faculty of Electrical Engineer- ing. His present work is on developing an exposure system for unrestrained small animalswhichwill serve for research on the non-thermal effects of electromag- netic fields and hyperthermic applicators. Barbora Vrbová Lukáš Víšek E-mail: vrbovbar@fel.cvut.cz, visekluk@fel.cvut.cz Dept. of Electromagnetic Field Faculty of Electrical Engineering Czech Technical University Technická 2, 166 27 Praha, Czech Republic 110