Quality assurance guidelines for interstitial hyperthermia.

Quality assurance (QA) guidelines are essential to provide uniform execution of clinical hyperthermia treatments and trials. This document outlines the clinical and technical consequences of the specific properties of interstitial heat delivery and specifies recommendations for hyperthermia administration with interstitial techniques. Interstitial hyperthermia aims at tumor temperatures in the 40-44 °C range as an adjunct to radiation or chemotherapy. The clinical part of this document imparts specific clinical experience of interstitial heat delivery to various tumor sites as well as recommended interstitial hyperthermia workflow and procedures. The second part describes technical requirements for quality assurance of current interstitial heating equipment including electromagnetic (radiative and capacitive) and ultrasound heating techniques. Detailed instructions are provided on characterization and documentation of the performance of interstitial hyperthermia applicators to achieve reproducible hyperthermia treatments of uniform high quality. Output power and consequent temperature rise are the key parameters for characterization of applicator performance in these QA guidelines. These characteristics determine the specific maximum tumor size and depth that can be heated adequately. The guidelines were developed by the ESHO Technical Committee with participation of senior STM members and members of the Atzelsberg Circle.

Comparison of MR-thermography and planning calculations in phantoms.

A systematic comparison of three-dimensional MR (magnetic resonance) thermography and planning calculations in phantoms for the hyperthermia (HT) SIGMA-Eye applicator. We performed 2 x 6 experiments in a homogeneous cylindrical and a heterogeneous elliptical phantom by adjusting 82 different patterns with different phase control inside an MR tomograph (Siemens Magnetom Symphony, 1.5 Tesla). For MR thermography, we employed the proton resonance frequency shift method with a drift correction based on silicon tubes. For the planning calculations, we used the finite-difference time-domain (FDTD) method and, in addition, modeled the antennas and the transforming network. We generated regions according to a segmentation of bones and tissue, and used an interpolation technique with a subgrid of 0.5 cm size at the interfaces. A Gauss-Newton solver has been developed to adapt phases and amplitudes. A qualitative agreement between the planning program and measurements was obtained, including a correct prediction of hot spot locations. The final deviation between planning and measurement is in the range of 2-3 W/kg, i.e., below 10%. Additional HT phase and amplitude adaptation, as well as position correction of the phantom in the SIGMA-Eye, further improve the results. HT phase corrections in the range of 30-40 degrees and HT amplitude corrections of +/- 20-30% are required for the best agreement. The deviation /MR-FDTD/, and the HT phase/amplitude corrections depend on the type of phantom, certain channel groups, pattern steering, and the positioning error. Appropriate agreement between three-dimensional specific absorption rate distributions measured by MR-thermography and planning calculations is achieved, if the correct position and adapted feed point parameters are considered. As long as feed-point parameters are uncertain (i.e., cannot be directly measured during therapy), a prospective planning will remain difficult. However, we can use the information of MR thermography to better predict the patterns in the future even without the knowledge of feed-point parameters.

Clinical evaluation and verification of the hyperthermia treatment planning system hyperplan.

PURPOSE: A prototype of the hyperthermia treatment planning system (HTPS) HyperPlan for the SIGMA-60 applicator (BSD Medical Corp., Salt Lake City, Utah, USA) has been evaluated with respect to clinical practicability and correctness. MATERIALS AND METHODS: HyperPlan modules extract tissue boundaries from computed tomography (CT) images to generate regular and tetrahedral grids as patient models, to calculate electric field (E-field) distributions, and to visualize three-dimensional data sets. The finite difference time-domain (FDTD) method is applied to calculate the specific absorption rate (SAR) inside the patient. Temperature distributions are calculated by a finite-element code and can be optimized. HyperPlan was tested on 6 patients with pelvic tumors. For verification, measured SAR values were compared with calculated SAR values. Furthermore, intracorporeal E-field scans were performed and compared with calculated profiles. RESULTS: The HTPS can be applied under clinical conditions. Measured absolute SAR (in W/kg), as well as relative E-field scans, correlated well with calculated values (+/-20%) using the contour-based FDTD method. Values calculated by applying the FDTD method directly on the voxel (CT) grid, were less well correlated with measured data. CONCLUSION: The HyperPlan system proved to be clinically feasible, and the results were quantitatively and qualitatively verified for the contour-based FDTD method.

Scanning E-field sensor device for online measurements in annular phased-array systems.

PURPOSE: A measurement device for noninvasive and simultaneous control of antennas during regional radiofrequency (rf) hyperthermia and, subsequently, the estimation of the power distribution in the interior of patients are essential preconditions for further technological progress. Aiming at this, the feasibility of an electro-optical electric field sensor was investigated during clinical rf hyperthermia. MATERIAL AND METHODS: The electro-optical electric field (E-field) sensor is based on lithiumniobate crystals and the Mach-Zehnder interferometer structure, and was tested in an earlier phantom study. For this study, a mechanical scanning device was developed allowing the registration of the E-field during clinical application. Data were recorded along a curve in the water bolus of the SIGMA 60 applicator of the annular phased-array system BSD-2000 (BSD Medical Corp., Salt Lake City, UT) close to the base points of the flat biconical dipole antennas. The results were compared with modeling calculations using the finite-difference time-domain (FDTD) method. For the latter, different antenna models were assumed. For systematic registration of the E-field curves in amplitude and phase, we employed an elliptical lamp phantom with fat-equivalent ring (filled with saline solution) and an elliptical polyacrylamide phantom with acrylic glass wall. Further measurements were carried out during the treatment of 5 patients with 20 hyperthermia treatments. RESULTS: Data of both phantom and patient measurements can be satisfactorily described by the FDTD method, if the antenna model is refined by taking into account the conical form of the dipoles and the special dielectric environment of the feeding point. Phase deviations can be entered ex posteriori for correction in the calculation algorithm. A comparison of amplifier power measurement (forward and backward power) and bolus E-field scans near the antenna base points demonstrates that E-field measurements between antennas and patient are a necessity for the appropriate characterization of antenna radiation properties. These measurements are sensitive to variations of the lossy medium in position and shape, and can be correctly predicted with current models. However, the differences between different patients are moderate and unspecific in both calculations and measurements, with fluctuations at maximum of 30 degrees in phases and 40% in amplitudes. CONCLUSIONS: The measurement method presented here turned out to be a practical tool for online registration of E-fields in phases and amplitudes along arbitrary curves in a water bolus or phantom. It can be utilized to evaluate antenna design and modeling calculations and leads, thus, to a better understanding of complicated multiantenna systems. In clinical routine, it can be employed as input for patient-specific hyperthermia planning and, finally, for the realization of online control with subsequent optimization of the power distribution in the patient.

Electric field distributions in a phased-array applicator with 12 channels: measurements and numerical simulations.

In this paper we examine the SIGMA-Eye hyperthermia applicator (BSD Medical Corp., Salt Lake City, Utah 84119) with respect to the control of electric field distributions. This applicator is equipped with 12 pairs of antennas fed by 12 amplifiers, allowing the individual adjustment of phase and power for each of them. Measurements were conducted using phantoms with well-defined electrical properties. Specific electro-optical sensors, capable of measuring both electric field amplitudes and phases, have been developed, and a system for data acquisition and analysis has been set up. In its initial state the applicator appeared not to be satisfactorily matched at 100 MHz for the phantom used, with return losses up to 20% in power. By tuner readjustments we achieved values below 5%. For various settings of the amplifiers' control parameters we measured field distributions, both in the phantom and in the surrounding water bolus. The experimental results were compared with numerical simulations based on finite difference and finite element methods. Measured and calculated electric fields exhibit deviations of 10% on average, allowing, in principle, a satisfactory prediction of fields by numerical simulations or as well by on-line measurements at selected locations of the applicator at antenna proximity. However, to obtain this satisfactory agreement a modification of the control parameters in the calculations (phases and amplitudes in the feed points of the antennas) was necessary. The origin of these problems is mainly attributed to cross-talk phenomena and other characteristics of the transforming network, which need to be scrutinized further for a full understanding.

Antenna arrays in the SIGMA-eye applicator: interactions and transforming networks.

OBJECTIVES: In multiantenna applicators such as the SIGMA-60 or SIGMA-Eye, which consist of 4 or 12 pairs of antennas shunt to 4 or 12 amplifiers ("antenna couplets"), phases and amplitudes in the feed points of these antennas under certain conditions can significantly differ from the values selected at the multichannel amplifier (forward parameters), mainly due to coupling. In the SIGMA-Eye, this interaction is particularly affected by the transforming networks between the generators and the feed points, thus hampering the control of the feed point parameters. In this work, we perform measurements at existing applicators, present a formalism to describe the facts numerically, and investigate modifications of the transforming networks to improve the performance. METHODS AND MATERIALS: We prepared an experimental setup for the SIGMA-Eye applicator that is fed by forward waves of a 12-channel amplifier system. In this setup, we made the water bolus, the interior of the tissue-equivalent phantom, and the entire transforming network accessible for measuring probes. Then, we constructed various alternative transforming networks such as Pawsey loops, LC matching networks, and power dividers and compared them with the original matching network of the SIGMA-Eye applicator. In particular, we utilized a high-resistive probe to determine the disturbances and influences caused by some channels with respect to some selected feed points of the SIGMA-Eye dipoles. RESULTS: In the original SIGMA-Eye applicator, the influences of coupling channels on the phases and voltages in the feed point of a particular antenna are largest for adjacent longitudinal channels. Here, the +/- 10 degrees phase shift and +/- 30% voltage change were observed if the reference channel (i.e., the disturbed channel) and disturbing channel are equally powered. The changes eminently increased to -30 degrees to + 100 degrees phase shift and -80% to +50% voltage change if the reference channel is fed with much lower power (four to eight-fold) than the disturbing channel. The disturbance from distant channels is less but still significant, reaching shifts of -10 degrees to +50 degrees and -50% to +20%, respectively. Using Pawsey loops instead of the original ferrite rings in the SIGMA-Eye network, the efficacy of the baluns was improved by a more than a factor of 4. Using an LC matching network, dependencies on frequency and external arrangements can be reduced significantly. Applying a power divider circuit, the coupling between antennas combined to one channel is considerably diminished (down to <-25 dB). CONCLUSION: Coupling between resonators (pairs of antennas including the matching network) reduces the control of the SIGMA-Eye applicator, i.e., it causes deviations between the selection of forward parameters at the amplifier and the total actual parameters in the feed points of the antennas. Modified transformation networks can improve the control, in particular by reducing sheath currents and asymmetries. There is a linear but variable relationship between selected (amplifiers) and actually given (feed points) parameters. This linear mapping (described by a matrix) and its characteristics need further investigation.

Visualization and registration of three-dimensional E-field distributions in annual-phased-array applicators.

A testing system is presented allowing registration, digitization, and evaluation of three-dimensional power distributions rendered by annular-phased-array applicators in homogeneous liquid media. The system is based on a lamp phantom originally developed to visualize power distributions. Now the brightness distribution is registered via a charge-coupled device camera and transferred to a PC-based evaluation system outside the shielding room. An appropriate mechanical coupling of camera and sensor matrix probing the phantom was built in order to keep optical image conditions constant under movement. For visualization and evaluation commercially customized software was employed. The evaluation of the system shows the linearity between sensor signal and power density magnitude to be sufficient for evaluation and graphical representation of three-dimensional data sets. In a first practical application the testing system was employed to evaluate dependencies of power distributions as a function of frequency and phase settings on temperatures and, subsequently, the relevance of those results for clinical hyperthermia in a SIGMA-60 applicator (BSD-2000 system). Now, the system is ready to evaluate more complex multiantenna array applicators like the SIGMA-Eye applicator. The measuring system is particularly suitable for a fast comparison of APA applicators applied for a homogeneous medium. Implications for heterogeneous structures (like in patients) are then possible via modeling calculations.

3-D computation of E fields by the volume-surface integral equation (VSIE) method in comparison with the finite-integration theory (FIT) method.

An algorithm has been developed for calculation of 3-dimensional E fields by the volume-surface integral equation (VSIE) method. Integration over surface elements is performed by elementary analytical formulas, assuming a linear interpolation of surface charges. Grid points at electrical interfaces are split off, well considering the E field behavior at these contours, specifically at sharp bends and multimedia junctions. Averaging procedures are utilized in order to avoid undefined or infinite values at critical points. The VSIE is solved by iteration using GMRES ("general minimum residuum") solver on a SUN workstation SPARC-IPX or Cray XMP, whereby convergence speed decreases considerably as the heterogeneity of the problem increases. Computation time (e.g., 20 min on a supercomputer for approximately 30,000 cells) needs to be reduced by further code development. Results for 3-D test cases (plane wave illuminating a layered cylinder) generally agree well with the finite-integration-theory (FIT) method if high E field gradients occur perpendicular to electrical boundaries. The VSIE method predicts slightly higher E fields only in critical regions. On the other hand, the FIT method at present is more efficient with respect to computation time for large domains with high cell numbers (> 100,000 cells).

[Part-body hyperthermia with a radiofrequency multiantenna applicator under online control in a 1.5 T MR-tomograph]. / Teilkörperhyperthermie mit einem Radiofrequenz-Multiantennen-Applikator unter online Kontrolle in einem 1.5 T MR-Tomographen.

Objective of this study is the integration of a multiantenna applicator for part-body hyperthermia (BSD 2000/3D) in a 1.5 T MR-tomograph (Siemens Magnetom Symphony) in order to perform noninvasive MR monitoring in real time to increase safety and effectiveness of heat treatments. The positioning unit is mechanically coupled to the MR gantry from the back side and the body coil is utilised for imaging. For that purpose, the hyperthermia antenna system (100 MHz, 1.500 W) and the MR receiver (63.9 MHs) have to be decoupled in terms of high frequency (filter) and electromagnetically (emc). The processing of MR data sets is performed in a hyperthermia planning system. A simultaneous operation of radiofrequency hyperthermia and MR system is possible at clinically relevant power levels. MR imaging is used for tumor-diagnostics (standard spin echo sequences), for hyperthermia planning (T1-weighted gradient echo sequences in equal- and opposed-phase techniques), and for temperature measurements according to the proton resonance frequency method (PRF method, phase evaluation registration using a gradient echo sequence with long echo time). In 33 patients with advanced pelvic and abdominal tumors we performed 150 heat sessions under MR monitoring. For 70% of these patients a visualisation of temperature sensitive data during treatment was possible. The evaluated difference images represent a superposition of real temperature -increase and a (temperature-induced) perfusion elevation. The -hybrid approach renders development of part body hyperthermia possible as an MR-controlled intervention in radiology.

Image-guided volumetric modulated arc therapy for breast cancer: a feasibility study and plan comparison with three-dimensional conformal and intensity-modulated radiotherapy.

OBJECTIVE: To test the feasibility of volumetric modulated arc therapy (VMAT) in breast cancer and to compare it with three-dimensional conformal radiotherapy (3D-CRT) as conventional tangential field radiotheraphy (conTFRT). METHODS: 12 patients (Stage I, 8: 6 left breast cancer and 2 right breast cancer; Stage II, 4: 2 on each side). Three plans were calculated for each case after breast-conserving surgery. Breast was treated with 50 Gy in four patients with supraclavicular lymph node inclusion, and in eight patients without the node inclusion. Multiple indices and dose parameters were measured. RESULTS: V95% was not achieved by any modality. Heterogeneity index: 0.16 (VMAT), 0.13 [intensity-modulated radiotherapy (IMRT)] and 0.14 (conTFRT). Conformity index: 1.06 (VMAT), 1.15 (IMRT) and 1.69 (conTFRT). For both indices, IMRT was more effective than VMAT (p=0.009, p=0.002). Dmean and V20 for ipsilateral lung were lower for IMRT than VMAT (p=0.0001, p=0.003). Dmean, V2 and V5 of contralateral lung were lower for IMRT than VMAT (p>0.0001, p=0.005). Mean dose and V5 to the heart were lower for IMRT than for VMAT (p=0.015, p=0.002). CONCLUSION: The hypothesis of equivalence of VMAT to IMRT was not confirmed for planning target volume parameter or dose distribution to organs at risk. VMAT was inferior to IMRT and 3D-CRT with regard to dose distribution to organs at risk, especially at the low dose level. ADVANCES IN KNOWLEDGE: New technology VMAT is not superior to IMRT or conventional radiotherapy in breast cancer in any aspect.

Methods and potentials of magnetic resonance imaging for monitoring radiofrequency hyperthermia in a hybrid system.

INTRODUCTION: Non-invasive thermometry (NIT) is a valuable and probably indispensable tool for further development of radiofrequency (RF) hyperthermia. A hybridization of an MRI scanner with a hyperthermia system is necessary for a real-time NIT. The selection of the best thermographic method is difficult, because many parameters and attributes have to be considered. METHODS: In the hybrid system (Siemens Symphony/BSD-2000-3D) the standard methods for NIT were tested such as T1, diffusion (ADC: apparent diffusion coefficient) and proton-resonance-frequency shift (PFS) method. A series of three-dimensional datasets was acquired with different gradient-echo sequences, diffusion-weighted EPI spin-echo sequences and calculated MR-temperatures in the software platform AMIRA-HyperPlan. In particular for the PFS-method, corrective methods were developed and tested with respect to drift and other disturbances. Experiments were performed in phantoms and the results compared with direct temperature measurements. Then the procedures were transferred to clinical applications in patients with larger tumours of the lower extremity or the pelvis. RESULTS: Heating experiments and MR-thermography in a homogeneous cylindrical phantom give an excellent survey over the potentials of the methods. Under clinical conditions all these methods have difficulties due to motion, physiological changes, inhomogeneous composition and susceptibility variations in human tissues. The PFS-method is most stable in patients yielding reasonable MR temperature distributions and time curves for pelvic and lower extremity tumours over realistic treatment times of 60-90 min. Pooled data exist for rectal tumour recurrencies and soft tissue sarcomas. The fat tissue can be used for drift correction in these patients. T1 and diffusion-dependent methods appear less suitable for these patients. The standard methods have different sensitivities with respect to the various error sources. The advantages and pitfalls of every method are discussed with respect to the literature and illustrated by the phantom and patient measurements. CONCLUSIONS: MR-controlled RF hyperthermia in a hybrid system is well established in phantoms and already feasible for patients in the pelvic and lower extremity region. Under optimal conditions the temperature accuracy might be in the range of 0.5 degrees C. However a variety of developments, especially sequences and post-processing modules, are still required for the clinical routine.

Strategies for optimized application of annular-phased-array systems in clinical hyperthermia.

A theoretical framework is presented for optimized heating of deep-seated tumours by phase and amplitude steering. The optimization problem for a specific tumour and perfusion case results in a functional dependency between power-level and maximum obtainable therapeutic efficiency. Different optimization criteria and strategies are outlined, which cause an increase of power or thermal dose in the tumour. Three tumour models (central pelvic tumour, eccentric abdominal tumour with or without necrosis) are analysed in detail. The simulation studies predict that appreciable parts of these tumours (50-100%) can be heated efficiently (42.5-43 degrees C) within the range of available and clinically tolerated power levels (1-5 kW/m), if tumour perfusion is less than 20-25 ml/100 g min. Some improvements are obtained by increasing the number of independent channels (from four to eight) and by the application of time-dependent (complementary) power-deposition patterns.

Simulation studies promote technological development of radiofrequency phased array hyperthermia.

A treatment planning program package for radiofrequency hyperthermia has been developed. It consists of software modules for processing three-dimensional computerized tomography (CT) data sets, manual segmentation, generation of tetrahedral grids, numerical calculation and optimisation of three-dimensional E field distributions using a volume surface integral equation algorithm as well as temperature distributions using an adaptive multilevel finite-elements code, and graphical tools for simultaneous representation of CT data and simulation results. Heat treatments are limited by hot spots in healthy tissues caused by E field maxima at electrical interfaces (bone/muscle). In order to reduce or avoid hot spots suitable objective functions are derived from power deposition patterns and temperature distributions, and are utilised to optimise antenna parameters (phases, amplitudes). The simulation and optimisation tools have been applied to estimate the improvements that could be reached by upgrades of the clinically used SIGMA-60 applicator (consisting of a single ring of four antenna pairs). The investigated upgrades are increased number of antennas and channels (triple-ring of 3 x 8 antennas and variation of antenna inclination. Significant improvement of index temperatures (1-2 degrees C) is achieved by upgrading the single ring to a triple ring with free phase selection for every antenna or antenna pair. Antenna amplitudes and inclinations proved as less important parameters.

Development and testing of SAR-visualizing phantoms for quality control in RF hyperthermia.

A new prototype of an elliptical standard phantom with fat-equivalent walls and a lamp matrix for SAR (specific absorption rate) visualization has been developed. This paper outlines the manufacture of solid components based upon either polyester resin or epoxy resin, as well as the adjustment of their electrical conditions (epsilon r, sigma) by admixtures of carbon and/or aluminium powder. Visualizing sensors (LED = light-emitting diodes, miniature lamps) are evaluated with respect to their transformation of electric field strength into light. Standard SAR patterns of the hyperthermia system BSD-2000 have been semiquantitatively assessed by the visualizing technique (power stepping method) and quantitatively by E field sensor scans. Extracted iso-SAR distributions are in good agreement with E field sensor scans performed with a lamp sensor coupled to a fibre or using a dipole probe with high resistive leads. The requirement for periodic quality control of SAR patterns of RF (radio frequency) hyperthermia systems is demonstrated. Comparisons between techniques are given, specifically with respect to the LED phantom of Schneider and van Dijk.

[Determinant factors and disturbances in controlling power distribution patterns by the hyperthermia-ring system BSD-2000. 2. Measuring techniques and analysis]. / Einflussfaktoren und Störeffekte bei der Steuerung von Leistungsverteilungen mit dem Hyperthermie-Ringsystem BSD-2000. Teil 2: Messtechnische analyse]

Clinical observables and phantom measurements (part 1) have suggested that the control of power deposition patterns can still be improved for the hyperthermia system BSD-2000. This is addressed to system-specific phase errors as well as inadequacies of phase selection (target point method), which might be corrected by modifications of the manufacturer. Furthermore, frequency-dependent physical effects (coupling, mode excitation) are existing, which might cause distortions and asymmetries of current distribution on antennas and consequently deteriorate the power deposition pattern (e.g. focussing capability). The application of a network analyzer system is described in order to determine electrical material constants, phase errors, coupling coefficients, reflection coefficients and current distributions on antennas. The analysis of the measurement datas suggests that ring-applicator has a variable frequency-optimum (supposed around 80 ... 95 MHz) characterized by minimal coupling and asymmetries.

[The influencing factors and interfering effects in the control of the power distributions with the BSD-20000 hyperthermia ring system. 1. The clinical observables and phantom measurements]. / Einflussfaktoren und Störeffekte bei der Steuerung von Leistungsverteilungen mit dem Hyperthermie-Ringsystem BSD-2000. Teil 1: Klinische Observablen und Phantommessungen.

A new generation of annular-phased-array systems BSD-2000 has been clinically applied in a pilot study. Therapeutic intratumoral temperatures greater than 42 degrees C were obtained in 15/15 sessions with six patients. However, the control of power deposition pattern has to be improved in order to increase the therapeutic gain and to guarantee an efficient therapy. A retrospective analysis of clinical phenomena has been performed by phantom set-ups because the power deposition pattern cannot be determined during therapy. Phantom measurement techniques are outlined, specifically phantom materials and visualization of power distributions. The problem of focus balance and frequency choice is illustrated by self-developed phantoms (liquid crystal sheets, light-emitting-diode-arrays). Especially, the limitation of modeling calculations is demonstrated.

Electromagnetic phased arrays for regional hyperthermia: optimal frequency and antenna arrangement.

This paper investigates the effects of the three-dimensional arrangement of antennae and frequency on temperature distributions that can be achieved in regional hyperthermia using an electromagnetic phased array. It compares the results of power-based and temperature-based optimization. Thus, one is able to explain the discrepancies between previous studies favouring more antenna rings on the one hand and more antennae per ring on the other hand. The sensitivity of the results is analysed with respect to changes in amplitudes and phases, as well as patient position. This analysis can be used for different purposes. First, it provides additional criteria for selecting the optimal frequency. Secondly, it can be used for specifying the required phase and amplitude accuracy for a real phased array system. Furthermore, it may serve as a basis for technological developments in order to reduce both types of sensitivities described above.

Influence of patient models and numerical methods on predicted power deposition patterns.

BACKGROUND: A hyperthermia planning system has been developed for generating patient and applicator models as well as calculating and visualizing E-field and temperature distributions. Significant dependencies on models and algorithms have been found. METHODS: Computerized tomography (CT) data sets are first transformed into so called 'labelled CT-volume'-data sets of equal resolution, which are used for segmentation. The first type of patient model obtained subsequently is based on regions with specified electrical properties representing tissues or organs (so called 'region-based model'). The second patient model renders a direct transformation of Hounsfield Units (HU) to electrical constants (so called 'HU-based model'). The FDTD-method (finite difference time domain) is then applied on a cubic lattice employing either an auxiliary 'sub-cubic lattice' (for HU-based segmentation) or a tetrahedron grid (for region-based segmentation) to assign the electrical properties, both representing the anatomy of the patient. E-field distributions are corrected by a post-processing procedure with respect to the geometry of interfaces defined by the tetrahedron grid. For comparison, the VSIE method (volume surface integral equation) is performed on the same tetrahedron grid. The applicator model assumes eight half-wavelength dipole antennas fed with constant voltages with water as background medium. RESULTS: For both numerical methods (FDTD, VSIE) the resulting antenna input impedances as well as the current distributions along the antennas were quite similar and almost insensitive to the particular geometry model (region-based, HU-based). In contrast to that, the power deposition patterns in the interior of the patient depended strongly on those models. Major differences can be related to different labels of the tissue type bone in the HU-based model in comparison to the definition via regions. Conversely, comparable results were obtained using the VSIE method and the FDTD method on the region-based patient model with a posteriori correction at the tetrahedron grid points. SAR (specific absorption rate) elevations up to a factor of 10 were predicted when employing region-based models. Those peaks might correspond to specific toxicity of electromagnetic radiation clinically known as hot spot phenomena or musculo-skeletal syndromes. Conversely, HU-based models generated quite homogeneous power deposition patterns with fluctuations of at most factor 2. CONCLUSION: The methods employing region-based geometry models such as the VSIE method and FDTD method in conjunction with a posteriori correction at tissue interfaces result in comparable E-field distributions for regional hyperthermia. Due to its shorter calculation time, the FDTD method is currently used in the clinic. Predictions derived from HU-based models without prior corrections of tissue specifications are not always supported by clinical experience.