Estimating the Electrical Conductivity of Human Tissue in Radiofrequency Hyperthermia Therapy

Abstract

The use of mathematical models to study complex systems such as physical and biological phenomena allows understanding their behavior, specifically regarding variables and parameters that are difficult to obtain. Additionally, studying optimization techniques has made it possible to approximate the characteristics of these systems by correlating numerical simulations and experimentation. Radiofrequency hyperthermia therapy for cancer treatment is currently under consideration for future medical applications. However, some of its properties are complex to measure, which could prevent their control. This is the case of electrical conductivity, which depends on the induction frequency and the tissue characteristics. In this paper, radiofrequency hyperthermia therapy was simulated via the finite element method. Then, an estimation of the electrical conductivity involved in the treatment was performed using the particle swarm optimization method. The execution time and the difference between the estimated parameter and the exact value were evaluated and compared with those obtained using the Levenberg-Marquardt method. The results indicate a significant agreement between the estimated and exact values in three different cases. The Levenberg-Marquardt method has a difference of 0,1942% and a performance time of 22 minutes, whereas the particle swarm optimization method has a difference of 0,0967% and a performance time of 327 minutes. The latter performs better in terms of parameter value estimation, whereas the former has better computational times. These techniques may help medical doctors to prescribe treatment protocols and may open the possibility of devising control strategies for hyperthermia therapy as a cancer treatment.

El uso de modelos matemáticos para el estudio de sistemas complejos como los fenómenos físicos y biológicos permite comprender su comportamiento, específicamente con respecto a variables y parámetros difíciles de obtener. Adicionalmente, el estudio de técnicas de optimización ha permitido aproximar las características de estos sistemas por medio de la correlación de simulaciones numéricas y la experimentación. La terapia de hipertermia por radiofrecuencia para el tratamiento del cáncer está actualmente en consideración para su futura aplicación médica. Sin embargo, algunas de sus propiedades son difíciles de medir, lo cual impediría su control. Este es el caso de la conductividad eléctrica, que depende de la frecuencia de inducción y de las características del tejido. En este artículo se simuló la terapia de hipertermia por radiofrecuencia mediante el método de elementos finitos. Luego se realizó una estimación de la conductividad eléctrica en el tratamiento mediante el método de optimización por enjambres de partículas. Se evaluaron el tiempo de ejecución y la diferencia del valor estimado con respecto al valor exacto, y se compararon sus valores estimados con los obtenidos mediante el método de Levenberg-Marquardt. Los resultados indican una concordancia significativa entre los valores estimados y los exactos en tres casos diferentes. El método de Levenberg-Marquardt tiene una diferencia de 0,1942% y un tiempo de ejecución de 22 minutos, mientras que el método de optimización de enjambres de partículas tiene una diferencia de 0,0967% y un tiempo de ejecución de 327 minutos. Este último tiene un mejor rendimiento en términos de estimación del valor de los parámetros, mientras que el otro tiene un mejor tiempo de ejecución computacional. Estas técnicas podrían ayudar a los médicos a prescribir protocolos de tratamiento y abrir la posibilidad de diseñar estrategias de control para la terapia de hipertermia como tratamiento para el cáncer.