Considerations for thermal injury analysis for RF ablation devices.

BACKGROUND: The estimation of lesion size is an integral part of treatment planning for the clinical applications of radiofrequency ablation. However, to date, studies have not directly evaluated the impact of different computational estimation techniques for predicting lesion size. In this study, we focus on three common methods used for predicting tissue injury: (1) iso-temperature contours, (2) Cumulative equivalent minutes, (3) Arrhenius based thermal injury. METHODS: We created a geometric model of a multi-tyne ablation electrode and simulated thermal and tissue injury profiles that result from three calculation methods after 15 minutes exposure to a constant RF voltage source. A hybrid finite element technique was used to calculate temperature and tissue injury. Time-temperature curves were used in the assessment of iso-temperature thresholds and the method of cumulative equivalent minutes. An Arrhenius-based formulation was used to calculate sequential and recursive thermal injury to tissues. RESULTS: The data demonstrate that while iso-temperature and cumulative equivalent minute contours are similar in shape, these two methodologies grossly over-estimate the amount of tissue injury when compared to recursive thermal injury calculations, which have previously been shown to correlate closely with in vitro pathologic lesion volume measurement. In addition, Arrhenius calculations that do not use a recursive algorithm result in a significant underestimation of lesion volume. The data also demonstrate that lesion width and depth are inadequate means of characterizing treatment volume for multi-tine ablation devices. CONCLUSIONS: Recursive thermal injury remains the most physiologically relevant means of computationally estimating lesion size for hepatic tumor applications. Iso-thermal and cumulative equivalent minute approaches may produce significant errors in the estimation of lesion size.

Thermal modeling of lesion growth with radiofrequency ablation devices.

BACKGROUND: Temperature is a frequently used parameter to describe the predicted size of lesions computed by computational models. In many cases, however, temperature correlates poorly with lesion size. Although many studies have been conducted to characterize the relationship between time-temperature exposure of tissue heating to cell damage, to date these relationships have not been employed in a finite element model. METHODS: We present an axisymmetric two-dimensional finite element model that calculates cell damage in tissues and compare lesion sizes using common tissue damage and iso-temperature contour definitions. The model accounts for both temperature-dependent changes in the electrical conductivity of tissue as well as tissue damage-dependent changes in local tissue perfusion. The data is validated using excised porcine liver tissues. RESULTS: The data demonstrate the size of thermal lesions is grossly overestimated when calculated using traditional temperature isocontours of 42 degrees C and 47 degrees C. The computational model results predicted lesion dimensions that were within 5% of the experimental measurements. CONCLUSION: When modeling radiofrequency ablation problems, temperature isotherms may not be representative of actual tissue damage patterns.

Numerical evaluation of heating of the human head due to magnetic resonance imaging.

In this paper, we present a numerical model for evaluating tissue heating during magnetic resonance imaging (MRI). Our method, which included a detailed anatomical model of a human head, calculated both the electromagnetic power deposition and the associated temperature elevations during an MRI head examination. Numerical studies were conducted using a realistic birdcage coil excited at frequencies ranging from 63 to 500 MHz. The model was validated both experimentally and analytically. The experimental validation was performed at the MR test facility located at the Food and Drug Administration's Center for Devices and Radiological Health.

An alternative method for determination of low-frequency specific absorption rate patterns in homogeneous phantoms.

We propose a new application of voltage gradient measurements to determine specific absorption rate (SAR) at low frequencies where quasi-static electromagnetic conditions apply. This method, which we call the voltage gradient method, relies on direct measurement of the voltage field rather than measurement of the electric field or thermal transients. The voltage gradient method is fast and can be implemented with voltmeters of moderate cost. We tested the voltage gradient method using normal saline, in a phantom with simple geometry, and a sine wave voltage source at 5, 10, 20 and 50 kHz. Compared to the SAR measured thermally in the same phantom, the voltage gradient method produced almost identical curves when normalized. When the results of the voltage gradient method were scaled to the same power level used for the thermal SAR, the agreement was compatible with typical thermal SAR accuracy.