Novel conformal technique to reduce staircasing artifacts at material boundaries for FDTD modeling of the bioheat equation

Validation and practical use of Plan2Heat hyperthermia treatment planning for capacitive heating

BACKGROUND

Capacitive devices are used for hyperthermia delivery, initially mainly in Asia, but nowadays also increasingly in Europe. Treatment planning can be very useful to determine the most effective patient-specific treatment set-up. This paper provides a validation of GPU-based simulations using Plan2Heat for capacitive hyperthermia devices.

METHODS

Validation was first performed by comparing simulations with an analytical solution for a spherical object placed inside a uniform electric field. Resolution was 5, 2.5 or 1 mm. Next, simulations for homogeneous and inhomogeneous phantom setups were performed for Thermotron RF8 and Celsius TCS capacitive heating devices at 2.5 mm resolution. Also different combinations of electrode sizes were evaluated. Normalized SAR profiles were compared to phantom measurements from the literature. Possible clinical use of treatment planning was demonstrated for an anal cancer patient, evaluating different treatment set-ups in prone and supine position.

RESULTS

Numerical and analytical solutions showed excellent agreement. At the center of the sphere, the error was 5.1%, 2.9% and 0.2% for a resolution of 5, 2.5 and 1 mm, respectively. Comparison of measurements and simulations for both Thermotron RF8 and Celsius TCS showed very good agreement within 5% for all phantom set-ups. Simulations were capable of accurately predicting the penetration depth; a very relevant parameter for clinical application. The patient case illustrated that planning can provide insight by comparing effectiveness of different treatment strategies.

CONCLUSION

Plan2Heat can rapidly and accurately predict heating patterns generated by capacitive devices. Thus, Plan2Heat is suitable for patient-specific treatment planning for capacitive hyperthermia.

Assessment of absorbed power density and temperature rise for nonplanar body model under electromagnetic exposure above 6 GHz

The averaged absorbed power density (APD) and temperature rise in body models with nonplanar surfaces were computed for electromagnetic exposure above 6 GHz. Different calculation schemes for the averaged APD were investigated. Additionally, a novel compensation method for correcting the heat convection rate on the air/skin interface in voxel human models was proposed and validated. The compensation method can be easily incorporated into bioheat calculations and does not require information regarding the normal direction of the boundary voxels, in contrast to a previously proposed method. The APD and temperature rise were evaluated using models of a two-dimensional cylinder and a three-dimensional partial forearm. The heating factor, which was defined as the ratio of the temperature rise to the APD, was calculated using different APD averaging schemes. Our computational results revealed different frequency and curvature dependences. For body models with curvature radii of >30 mm and at frequencies of >20 GHz, the differences in the heating factors among the APD schemes were small.

Real-Time Three-Dimensional Microwave Monitoring of Interstitial Thermal Therapy

We report a method for real-time three-dimensional monitoring of thermal therapy through the use of noncontact microwave imaging. This method is predicated on using microwaves to image changes in the dielectric properties of tissue with changing temperature. Instead of the precomputed linear Born approximation that was used in prior work to speed up the frame-to-frame inversions, here we use the nonlinear distorted Born iterative method (DBIM) to solve the electric volume integral equation (VIE) to image the temperature change. This is made possible by using a recently developed graphic processing unit accelerated conformal finite difference time domain method to solve the forward problem and update the electric field in the monitored region in each DBIM iteration. Compared to our previous work, this approach provides a far superior approximation of the electric field within the VIE, and thus yields a more accurate reconstruction of tissue temperature change. The proposed method is validated using a realistic numerical model of interstitial thermal therapy for a deep-seated brain lesion. With the new DBIM, we reduced the average estimation error of the mean temperature within the region of interest from 2.5 to 1.0 for the noise-free case, and from 2.9 to 1.7 for the 2% background noise case.

Desktop exposure system and dosimetry for small scale in vivo radiofrequency exposure experiments

This paper describes a new approach to the risk assessment of exposure from wireless network devices, including an exposure setup and dosimetric assessment for in vivo studies. A novel desktop reverberation chamber has been developed for well-controlled exposure of mice for up to 24 h per day to address the biological impact of human exposure scenarios by wireless networks. The carrier frequency of 2.45 GHz corresponds to one of the major bands used in data communication networks and is modulated by various modulation schemes, including Global System for Mobile Communications (GSM), Universal Mobile Telecommunications System (UMTS), Radio Frequency Identification (RFID), and wireless local area network, etc. The system has been designed to enable exposures of whole-body averaged specific absorption rate (SAR) of up to 15 W/kg for six mice of an average weight of 25 g or of up to 320 V/m incident time-averaged fields under loaded conditions without distortion of the signal. The dosimetry for whole-body SAR and organ-averaged SAR of the exposed mice, with analysis of uncertainty and variation analysis, is assessed. The experimental dosimetry based on temperature measurement agrees well with the numerical dosimetry, with a very good SAR uniformity of 0.4 dB in the chamber. Furthermore, a thermal analysis and measurements were performed to provide better understanding of the temperature load and distribution in the mice during exposure.

An accelerated photo-magnetic imaging reconstruction algorithm based on an analytical forward solution and a fast Jacobian assembly method

We previously introduced photo-magnetic imaging (PMI), an imaging technique that illuminates the medium under investigation with near-infrared light and measures the induced temperature increase using magnetic resonance thermometry (MRT). Using a multiphysics solver combining photon migration and heat diffusion, PMI models the spatiotemporal distribution of temperature variation and recovers high resolution optical absorption images using these temperature maps. In this paper, we present a new fast non-iterative reconstruction algorithm for PMI. This new algorithm uses analytic methods during the resolution of the forward problem and the assembly of the sensitivity matrix. We validate our new analytic-based algorithm with the first generation finite element method (FEM) based reconstruction algorithm previously developed by our team. The validation is performed using, first synthetic data and afterwards, real MRT measured temperature maps. Our new method accelerates the reconstruction process 30-fold when compared to a single iteration of the FEM-based algorithm.

Graphics Processing Unit-Based Bioheat Simulation to Facilitate Rapid Decision Making Associated with Cryosurgery Training

This study focuses on the implementation of an efficient numerical technique for cryosurgery simulations on a graphics processing unit as an alternative means to accelerate runtime. This study is part of an ongoing effort to develop computerized training tools for cryosurgery, with prostate cryosurgery as a developmental model. The ability to perform rapid simulations of various test cases is critical to facilitate sound decision making associated with medical training. Consistent with clinical practice, the training tool aims at correlating the frozen region contour and the corresponding temperature field with the target region shape. The current study focuses on the feasibility of graphics processing unit-based computation using C++ accelerated massive parallelism, as one possible implementation. Benchmark results on a variety of computation platforms display between 3-fold acceleration (laptop) and 13-fold acceleration (gaming computer) of cryosurgery simulation, in comparison with the more common implementation on a multicore central processing unit. While the general concept of graphics processing unit-based simulations is not new, its application to phase-change problems, combined with the unique requirements for cryosurgery optimization, represents the core contribution of the current study.

Simulation techniques in hyperthermia treatment planning

Abstract Clinical trials have shown that hyperthermia (HT), i.e. an increase of tissue temperature to 39-44 °C, significantly enhance radiotherapy and chemotherapy effectiveness [1]. Driven by the developments in computational techniques and computing power, personalised hyperthermia treatment planning (HTP) has matured and has become a powerful tool for optimising treatment quality. Electromagnetic, ultrasound, and thermal simulations using realistic clinical set-ups are now being performed to achieve patient-specific treatment optimisation. In addition, extensive studies aimed to properly implement novel HT tools and techniques, and to assess the quality of HT, are becoming more common. In this paper, we review the simulation tools and techniques developed for clinical hyperthermia, and evaluate their current status on the path from 'model' to 'clinic'. In addition, we illustrate the major techniques employed for validation and optimisation. HTP has become an essential tool for improvement, control, and assessment of HT treatment quality. As such, it plays a pivotal role in the quest to establish HT as an efficacious addition to multi-modality treatment of cancer.

Thermal tissue damage model analyzed for different whole-body SAR and scan durations for standard MR body coils

PURPOSE

This article investigates the safety of radiofrequency induced local thermal hotspots within a 1.5T body coil by assessing the transient local peak temperatures as a function of exposure level and local thermoregulation in four anatomical human models in different Z-positions.

METHODS

To quantize the effective thermal stress of the tissues, the thermal dose model cumulative equivalent minutes at 43°C was employed, allowing the prediction of thermal tissue damage risk and the identification of potentially hazardous MR scan-scenarios. The numerical results were validated by B1 (+) - and skin temperature measurements.

RESULTS

At continuous 4 W/kg whole-body exposure, peak tissue temperatures of up to 42.8°C were computed for the thermoregulated model (60°C in nonregulated case). When applying cumulative equivalent minutes at 43°C damage thresholds of 15 min (muscle, skin, fat, and bone) and 2 min (other), possible tissue damage cannot be excluded after 25 min for the thermoregulated model (4 min in nonregulated).

CONCLUSION

The results are found to be consistent with the history of safe use in MR scanning, but not with current safety guidelines. For future safety concepts, we suggest to use thermal dose models instead of temperatures or SAR. Special safety concerns for patients with impaired thermoregulation (e.g., the elderly, diabetics) should be addressed.

Dominant factors affecting temperature rise in simulations of human thermoregulation during RF exposure

Numerical models of the human thermoregulatory system can be used together with realistic voxel models of the human anatomy to simulate the body temperature increases caused by the power absorption from radio-frequency electromagnetic fields. In this paper, the Pennes bioheat equation with a thermoregulatory model is used for calculating local peak temperatures as well as the body-core-temperature elevation in a realistic human body model for grounded plane-wave exposures at frequencies 39, 800 and 2400 MHz. The electromagnetic power loss is solved by the finite-difference time-domain (FDTD) method, and the discretized bioheat equation is solved by the geometric multigrid method. Human thermoregulatory models contain numerous thermophysiological and computational parameters--some of which may be subject to considerable uncertainty--that affect the simulated core and local temperature elevations. The goal of this paper is to find how greatly the computed temperature is influenced by changes in various modelling parameters, such as the skin blood flow rate, models for vasodilation and sweating, and clothing and air movement. The results show that the peak temperature rises are most strongly affected by the modelling of tissue blood flow and its temperature dependence, and mostly unaffected by the central control mechanism for vasodilation and sweating. Almost the opposite is true for the body-core-temperature rise, which is however typically greatly lower than the peak temperature rise. It also seems that ignoring the thermoregulation and the blood temperature increase is a good approximation when the local 10 g averaged specific absorption rate is smaller than 10 W kg(-1).

Children and adults exposed to electromagnetic fields at the ICNIRP reference levels: theoretical assessment of the induced peak temperature increase

To avoid potentially adverse health effects of electromagnetic fields (EMF), the International Commission on Non-Ionizing Radiation Protection (ICNIRP) has defined EMF reference levels. Restrictions on induced whole-body-averaged specific absorption rate (SAR(wb)) are provided to keep the whole-body temperature increase (T(body, incr)) under 1 °C during 30 min. Additional restrictions on the peak 10 g spatial-averaged SAR (SAR(10g)) are provided to prevent excessive localized tissue heating. The objective of this study is to assess the localized peak temperature increase (T(incr, max)) in children upon exposure at the reference levels. Finite-difference time-domain modeling was used to calculate T(incr, max) in six children and two adults exposed to orthogonal plane-wave configurations. We performed a sensitivity study and Monte Carlo analysis to assess the uncertainty of the results. Considering the uncertainties in the model parameters, we found that a peak temperature increase as high as 1 °C can occur for worst-case scenarios at the ICNIRP reference levels. Since the guidelines are deduced from temperature increase, we used T(incr, max) as being a better metric to prevent excessive localized tissue heating instead of localized peak SAR. However, we note that the exposure time should also be considered in future guidelines. Hence, we advise defining limits on T(incr, max) for specified durations of exposure.

Fast FFT-based bioheat transfer equation computation

This paper describes a modeling method of the tissue temperature evolution over time in hyper or hypothermia. The tissue temperature evolution over time is classically described by Pennes' bioheat transfer equation which is generally solved by a finite difference method. In this paper we will present a method where the bioheat transfer equation can be algebraically solved after a Fourier transformation over the space coordinates. As an example, we implemented this method for the simulation of a percutaneous high intensity ultrasound hepatocellular carcinoma curative treatment and compared it with the finite difference method and experimental data.

The Virtual Family--development of surface-based anatomical models of two adults and two children for dosimetric simulations

The objective of this study was to develop anatomically correct whole body human models of an adult male (34 years old), an adult female (26 years old) and two children (an 11-year-old girl and a six-year-old boy) for the optimized evaluation of electromagnetic exposure. These four models are referred to as the Virtual Family. They are based on high resolution magnetic resonance (MR) images of healthy volunteers. More than 80 different tissue types were distinguished during the segmentation. To improve the accuracy and the effectiveness of the segmentation, a novel semi-automated tool was used to analyze and segment the data. All tissues and organs were reconstructed as three-dimensional (3D) unstructured triangulated surface objects, yielding high precision images of individual features of the body. This greatly enhances the meshing flexibility and the accuracy with respect to thin tissue layers and small organs in comparison with the traditional voxel-based representation of anatomical models. Conformal computational techniques were also applied. The techniques and tools developed in this study can be used to more effectively develop future models and further improve the accuracy of the models for various applications. For research purposes, the four models are provided for free to the scientific community.