A Three-Dimensional Finite Element Model to Study Dynamic Changes of Tissue during Laser Interstitial Thermotherapy

Transcription

1 American Journal of Biomedical Engineering 2015, 5(3): DOI: /j.ajbe A Three-Dimensional Finite Element Model to Study Dynamic Changes of Tissue during Laser Interstitial Thermotherapy Mohammad Ehsan Salavati *, Mohammad Hossein Miran Baygi Faculty of Electrical and Computer Engineering, Tarbiat Modares University, Tehran, Iran Abstract Laser Interstitial Thermotherapy (LITT) is a surgical method for treatment of cancer in biological tissues such as liver, brain, prostate and neck. Mathematical modeling of the processes in laser-tissue interactions during laser irradiation is essential for optimal treatment planning. One of the parameters that affect the temperature and damage increase in the tissue during LITT is dynamic changes of tissue properties. In this study, a three-dimensional model was developed to model the effect of dynamic changes of tissue properties. Monte Carlo method was used for modeling the light transport in tissue. Finite Element method was used for solving Pennes and Arrhenius equations to obtain distribution of temperature and thermal damage in tissue respectively. Three conditions were considered in this study. The results show that with constant tissue properties the highest temperature and the largest damage volume occur. The temperature increase and damage volume are larger when considering dynamic changes of all properties rather than blood perfusion rate only. Keywords Laser Interstitial Thermotherapy, Dynamic changes, Finite Element Method, Monte Carlo Method, Laser-tissue interactions 1. Introduction Treatment of cancer has become one of the top priorities in health care centers worldwide and a large amount of resources are invested to improve or find new therapeutic methods for the treatment of cancer. One of the promising and alternative methods of therapy is laser-induced interstitial thermotherapy (LITT). LITT is a minimally invasive treatment to cause hyperthermia and coagulation in tumours in solid organs such as liver, brain, prostate, neck, etc. [1]. LITT treatment is usually performed by inserting a laser applicator into the tumour as shown in Figure 1. The tumour is destroyed by the immediate or delayed hyperthemic and coagulative effects, due to photon absorption and heat transfer in the tissue [2, 3]. In order to achieve a safe LITT procedure, a pre-knowledge of the final damage zone is essential in order to prevent injury to adjacent healthy tissues and sensitive structures. Several monitoring methods have been developed to achieve safe LITT procedures such as using thermocouples, fibre-optic temperature probes, color-coded Doppler * Corresponding author: me.salavati@gmail.com (Mohammad Ehsan Salavati) Published online at Copyright 2015 Scientific & Academic Publishing. All Rights Reserved sonography, X-ray imaging and magnetic resonance imaging (MRI). But all these methods have problems with accuracy or patient convenience. The best method suggested till now is MRI. But this method is expensive and errors occur with patient movements. Figure 1. General Diagram of LITT treatment Therefore, according to problems of temperature monitoring methods mentioned above, numerical modelling of the laser-tissue interaction process is still an important, economical and effective way to evaluate a wide range of parameters for a desired outcome without extensive in vivo trials. Nd:YAG lasers (1064 nm) or different types of diode lasers ( nm) are usually used as light sources for LITT treatments because of their deeper penetration depths in tissues in the NIR-range [4].

2 American Journal of Biomedical Engineering 2015, 5(3): Conventional modelling has normally included the following three steps: (1) calculation of the laser energy distribution using optical models; (2) calculation of the temperature increase using bio-heat transfer equations; (3) calculation of the extent of the thermal damage [5-12]. In order to have an accurate simulation which is closer to what happens in reality, one should consider all parameters that affect the final results. One of these effective parameters is dynamic changes of tissue. The increase in tissue temperature changes tissue properties. Thermal, optical and blood perfusion rate properties of tissue are changed when the tissue is heated [5-8]. The main purpose of this article is to study the dynamic changes of Liver properties during LITT. Similar studies have been done with 2D models and with numerical methods such as Finite Difference method [5-8]. This study adds more information about dynamic changes of tissue properties using a 3D model with Monte Carlo for light transport in tissue and Finite Element method (FEM) for calculating temperature and damage increase in tissue which are more accurate and a step closer to reality. No study has done before using Monte Carlo and FEM simultaneously for LITT before. 2. Methodology For simplicity, in this study we have used a three-layer cubic model. The tumour is the layer between two healthy tissue layers. The laser applicator is placed horizontally in the middle of the tumour layer along x-axis. The geometry of this model is shown in Figure 2. Each layer has the length and width of 4.5 cm and the height of 1.5 cm. So the model is a cube of 4.5*4.5*4.5 cm 3. Figure 2. Geometry of the 3D model 2.1. Monte Carlo Simulation of Photon Propagation in Tissues In biological tissues, light is both absorbed and scattered. Monte Carlo simulation of photon propagation offers a flexible and accurate approach toward photon transport in turbid materials, such as biological tissues; and it can easily deal with two or three dimensional problems with complex geometry [13]. Propagation steps in Monte Carlo are either fixed or variable [14-16]. Once the photon is launched, it may be scattered, absorbed, propagated undisturbed, internally reflected or transmitted out of tissue. The photon is repeatedly moved until it escapes from the tissue or it is totally absorbed. If the photon escapes from the tissue, the refection or transmission of the photon is recorded. If the photon is absorbed, the position of the absorption is recorded. Once this has been done a new photon is launched at the origin. This process is repeated until the desired numbers of photons have been propagated. The recorded reflection, transmission, and absorption profiles will approach true values. The absorption can be decided in two ways: The usual method that after each propagation step, the photon's weight is reduced by the probability of absorption or the implicit capture technique that provides some absorption information at each photon step, rather than just at times when the photon is completely absorbed [14, 15]. In this article we have used the variable step size Monte Carlo with implicit capture technique to reach desirable accuracy will less photons. At the end of the Monte Carlo simulation, an absorption matrix is obtained which is the raw matrix of total photon weight in each grid element. This matrix with the below formulation then is turned into the absorption photon probability density [14-18]: A(x,y,z) Bxyz (,, ) dx. dy. dz. N (1) Where B is the raw matrix of total photon weight, dx, dy and dz are volume elements in Cartesian coordinates and N is the number of photons used for the simulation. The obtained matrix from (1) is then multiplied to the laser power to give the laser volume density in tissue The Heat Distribution Equation in Tissue The laser light absorption results in internal energy increase in the local tissue. The standard Pennes bio-heat equation is still the most commonly used model to describe the thermal behavior of living tissue [9-12, 19-22]: T A pc.. =.[ k T] + pb. ωb. cb.[ Ta T] + S t where T is the temperature [ C], ρ is the density of tissue [kg.cm 3 ], c is the specific heat of tissue [J.kg 1.C 1 ], k is the thermal conductivity of tissue [W.cm 1.C 1 ], t is the time [s], T art is the temperature of arterial blood [C ], ω b is the perfusion rate [s -1 ], c b is the specific heat of blood [J.kg 1.C 1 ] and S is the deposited light power [W.cm 3 ] that comes from Monte Carlo The Damage Function According to the Arrhenius formulation, the degree of (2)

3 88 Mohammad Ehsan Salavati et al.: A Three-Dimensional Finite Element Model to Study Dynamic Changes of Tissue during Laser Interstitial Thermotherapy tissue damage can be quantified by the damage index Ω [20-27]: cr (,0) Ω ( r, τ ) = ln( ) cr (, τ ) τ = 0 Ea Aexp( ) d RT Where Ω is the dimensionless damage index, E a is the activation energy [J.mol 1 ], A is the frequency factor [s 1 ], R is the universal gas constant [J.mol 1.K 1 ] and τ is the total heating time [26, 27]. The activation energy E a and the frequency factor A are derived from thermodynamic variables. They can be related to the denaturation process of proteins and other cellular constituents. The above equation indicates that the damage index Ω is the logarithm of the ratio of the initial concentration of undamaged tissue to the concentration of healthy molecules after damage accumulated for the time interval t = 0 to t = τ. Therefore, Ω = 1 corresponds to the reduction in concentration of native molecules to 37% and it is usually used for the threshold of irreversible damage of the tissue The Dependences of the Tissue Properties on Tissue Temperature and Damage Although no complete temperature and damage dependent thermal properties have been reported, tissue thermal properties can be approximated according to the water content and the dynamic changes of the thermal properties of water in the range of C [9-12]: (3) ρ( T) = 1000( κ ρ w) (4) kt ( ) = 0.419( κ w) (5) ct ( ) = 4190( κ w) (6) Where w is the water mass percentage in the tissue which was considered 70% in this paper and T [ C] is the tissue temperature. k ρ, k c and k k are temperature-dependent factors for the density, specific heat and thermal conductivity of water: κ ρ ( T 20 ο C) = (7) c k κc T C 4 = ( 20 ο ) (8) 3 κk = ( T 20 ο C) (9) The dynamic changes of the blood perfusion rate with temperature and damage can be formulated as [9, 10]: ω ( T, Ω ) = ω f f (10) b b0 T u Where ω b0 is the blood perfusion rate and is ω b0 = [s -1 ] is what we used for liver in this study. f u is the undamaged fraction of the tissue is expressed as f = exp( Ω). f T is a dimensionless function that accounts u for vessel dilation at slightly elevated temperatures which can be approximated as [9, 10]: f T ο ( T 42 C) = 4 3. Model Implementation ο ο 37 C T < 42 C ο T 42 C (11) Nd:YAG lasers (1064 nm) or different types of diode lasers ( nm) are usually used as light sources for LITT treatments because of their deeper penetration depths in tissues in the NIR-range [2]. In this study we used 850 nm diode lasers. The values used for this wavelength are listed in Table 1. We have used properties of healthy and cancerous liver for this study [9, 10]. The whole process is done in 3 steps: 1) first the distribution of light in tissue is modeled by Monte Carlo. In each Monte Carlo simulation we have used photons. 2) The heat generation and transport is calculated using Pennes bio-heat equation. 3) The thermal damage is calculated using Arrhenius equation. Both Pennes bio-heat and Arrhenius equations are modeled using FEM with Comsol Multiphysics 3.5 [28] which utilizes FEM to solve the related governing equations. For irreversible damage, the damage threshold is Ω=1. In LITT, for coagulating the tissue and prevent bleeding laser powers are between 1 to 5W. The contact applicators for this purpose usually have a diameter up to 400µm. In our simulations we used the laser applicator is 10mm long with diameter of 200µm and the laser power of 5W. Table 1. Optical and Thermal Properties of Liver for 850nm Diode Laser [9, 10] Wavelength Tissue µ a[1/cm] µ s[1/cm] n C[J/kg.k] ρ[kg/m 3 ] k[w/m.k] T[ C] 850nm Healthy nm cancerous

4 American Journal of Biomedical Engineering 2015, 5(3): Results At first, the dynamic changes of blood perfusion rate is not considered and dynamic changes are considered only for heat conductivity, density and specific heat and the results are compared with constant thermal properties. Then the dynamically changing blood perfusion rate results are considered, when other properties are constant. At the end, the results for all temperature-dependent thermal properties are considered and these three conditions are compared with each other Dynamic Changes of Temperature-Dependent Tissue Properties Using equations (4) to (9) in the 3-layer cubic model, the final temperature after 60 seconds with constant and temperature-dependent properties are C and C respectively. The final damage volumes are also [cm 3 ] and [cm 3 ] respectively. In the Figures 3 and 4, the temperature increase and damage volume increase during the laser irradiation with constant and changing tissue properties are shown. As can be seen from the results and the Figure, considering the dependence of specific heat, thermal conductivity and tissue density, cause temperature and damage volume to diminish. The cause of this matter is shown in figure 5, where the changes of specific heat, thermal conductivity and tissue density with temperature increase in tissue, is plotted. Accordingly, the specific heat increased 0.36%, tissue density rose about 0.53% and thermal conductivity increased 8.2%. These results agree well with other studies [9, 10 and 12]. Figure 3. Temperature increase during laser irradiation at the center of the tissue with constant and changing thermal properties Figure 4. Damage volume increase during laser irradiation at the center of the tissue with constant and changing thermal properties Figure 5. Dynamic changes of tissue properties during laser irradiation. a) specific heat. b) density. c) heat conductivity

5 90 Mohammad Ehsan Salavati et al.: A Three-Dimensional Finite Element Model to Study Dynamic Changes of Tissue during Laser Interstitial Thermotherapy 4.2. Dynamic Changes of Blood Perfusion Rate Blood perfusion rate in tissue has a large impact on heat transfer and damage. Blood perfusion rate with its convective cooling effect of tissue during laser irradiation is an important property in LITT. It is one of the properties of tissue that changes with temperature increase. For modeling blood perfusion rate changes during LITT, we used equations (10) and (11). In this part due to high computational calculations, the simulation time is set to 40 seconds. Figure 6. Dynamic changes of tissue properties during laser irradiation. a) specific heat. b) density. c) heat conductivity The results show that final temperature and damage volume with constant properties and with temperature and damage dependent blood perfusion rate are C [cm 3 ] and C [cm 3 ] respectively. The temperature increase during laser irradiation in the center of the tissue is shown in figure 6. Due to expansion of blood vessels and increment of blood perfusion rate, cooling properties of blood in the tissue will have a greater impact. Therefore the final temperature is lower comparing with constant tissue properties. After the 30 th second, an increase in the slope of the temperature is seen. This is when the blood vessels are destroyed and coagulated due to high temperature in the tissue. The temperature expansion in line of applicator position (x axis) and tissue depth (z axis) in figure 7, also shows that with dynamic changes of blood perfusion rate, the cooling effect of blood increases and the depth of heated tissue is decreased Dynamic Changes of All Temperature-Dependent Tissue Properties Including Blood Perfusion Rate After considering all thermal and damage dependent properties such as specific heat, heat conductivity, density and blood perfusion rate, the final temperature and damage volume are C and 7.96[cm 3 ] respectively. In figure 8 the temperature increase in center point of the tissue after 40 seconds of laser irradiation is shown with 3 conditions of constant properties, dynamic changes of blood perfusion rate only and dynamic changes of all temperature and damage dependent properties. The temperature expansion in line of applicator position (x axis) and tissue depth (z axis) for these 3 conditions are also shown in figure 9. Figure 7. Temperature expansion of tissue with constant thermal properties and with changing blood perfusion rate. a) In line of applicator position (x axis). b) Tissue depth (z axis)

6 American Journal of Biomedical Engineering 2015, 5(3): Figure 8. Temperature increase during laser irradiation at the center of the tissue with constant thermal properties, dynamic changes of only blood perfusion rate and dynamic changes of all temperature and damages properties Figure 9. Temperature expansion with constant thermal properties, dynamic changes of only blood perfusion rate and dynamic changes of all temperature and damages properties.. a) In line of applicator position (x axis). b) Tissue depth (z axis) According to these results and considering three conditions studied in this article, with constant tissue properties the highest temperature and the largest damage volume occurs. The temperature increase and damage volume is higher when considering dynamic changes of all properties rather than only blood perfusion rate. The increase in the slope of temperature due to coagulation of blood vessels after the 30 th second is higher when considering dynamic changes of all temperature and damage dependent properties. In reality due to coagulation of tissue and dilation of blood vessels, these conditions happen. Hence, dynamic changes of temperature and damage dependent properties should be considered to have a more reliable result and closer to reality. 5. Conclusions and Discussion Laser Interstitial Thermotherapy (LITT) is a minimally invasive surgical method to treat cancerous tumours in biological tissues. Using lasers for surgery has many advantages over conventional surgery methods such as increased accuracy, less bleeding and less damage to the tissue. LITT is an effective treatment but has limitations. In order to achieve a safe LITT procedure, a pre-knowledge of the final damage zone is essential in order to prevent injury to adjacent healthy tissues and sensitive structures. One of the parameters that affect the temperature increase in biological tissues is dynamic changes of tissue properties during laser irradiation. In this study, a 3D model was

7 92 Mohammad Ehsan Salavati et al.: A Three-Dimensional Finite Element Model to Study Dynamic Changes of Tissue during Laser Interstitial Thermotherapy developed to model the effect of dynamic changes of tissue properties and study the final results considering these changes. For that matter, Monte Carlo method with variable step size and implicit capture technique was used to simulate photon transport in tissue. The Pennes bio-heat equation was used for heat generation in tissue. Arrhenius damage function was used to predict the damage volume in tissue. These equations were solved with Finite Element Method simultaneously using Comsol Multiphysics 3.5. Three conditions were considered for studying dynamic changes during LITT: 1) with constant properties, 2) dynamic changes of blood perfusion rate only, 3) dynamic changes of specific heat capacity, heat conductivity, density and blood perfusion rate. The results showed that with constant tissue properties the highest temperature and the largest damage volume occur. The temperature increase and damage volume is higher when considering dynamic changes of all properties rather than blood perfusion rate only. Considering all dynamic changes in tissue properties, this model provides a better understanding of the complex process of laser-tissue interactions during LITT. Hence the physician can have a better planning for LITT treatment. However, these computer simulations should be verified with in vivo experiments in controlled experimental laboratories and with measured optical and thermal properties for any specific patient. REFERENCES [1] Enevoldsen M. S., O. Skovgaard, A Combined Mathematical-Physical Model of Laser-induced Thermotherapy (LITT), European Conferences on Biomedical Optics, Therapeutic Laser Applications and Laser-Tissue Interactions (ECBO), Germany, Munich [2] Niemz M.H., Laser-Tissue Interaction: Fundamentals and Applications, Springer-verlag, Berlin, Heidlberg [3] Muller G. J., A. Roggan, Laser-Induced Interstitial Thermotherapy, The International Society for Optical Engineering, Bellingham, Washington, [4] Oth K., D. Russ, Laser Coagulation Zone Induced with the ND:YAG Laser in the Liver, Laser in Medical Science 1997,Vol.12, pp [5] Vogl T. J., R. Straub, S. Zangos, M. G. Mack, K. Eichler, M. Bottger, Laser-induced thermotherapy of malignant liver tumours: general principals, equipments, procedures, side effects, complications and results, European Journal of Ultrasound, Vol.13, pp , [6] Vogl T. J., R. Straub, S. Zangos, M. G. Mack, K. Eichler, Mohamed Nabil, Effect of laser-induced thermotherapy of liver metastases, Expert Rev, Anticancer 6, pp , [7] Vogl T. J., R. Straub, S. Zangos, M. G. Mack, MR-guided laser-induced thermotherapy (LITT) of liver tumours: experimental and clinical data, International Journal of Hyperthermia, Vol.20, No 7, pp , [8] Germer C. T., C. M. Isbert, D. Albrecht, J. P. Ritz, A. Roggan, K. J. Wolf, Laser-induced thermotherapy for the treatment of liver metastasis, Surgical Endoscopy, Springer-Verlag, Vol.12, pp , New York, [9] Iizuka M. N., A. Vitkin, M. C. Kolios, The effect of dynamic optical properties during interstitial laser photocoagulation, Phys. Med. Biol. 45, , [10] Jiang S. H., X. Zhang, Effects of dynamic changes of tissue properties during laser-induced interstitial thermotherapy (LITT), Laser in Medical Science, 19: , [11] Jiang S. H., X. Zhang, Dynamic modeling of photothermal interactions for laser-induced interstitial thermotherapy: parameter sensitivity analysis, Lasers in Medical Science, 20: ,2005. [12] Kim B. M., S.L. Jacques, S. Rastegar, S. Thomsen, M. Motamedi, Nonlinear Finite Element analysis of the role of dynamic changes in blood perfusion and optical properties in laser coagulation of tissue, IEEE Journal of Selected Topics in Quantum Electronics, Vol.2, No.4, December [13] Zhou J., Y. Zhang, J.K. Chen, A dual reciprocity boundary element method for photothermal interactions in laser-induced thermotherapy, international Journal of Heat and Mass Transfer, vol. 53,pp , [14] Jacques S. L., L. Wang, Monte Carlo Modeling of light transport in multi-layered tissue in standard C, Oregon Medical Laser Center, [15] Prahl S., Light transport in tissue, PhD Thesis in the university of Texas at Austin, [16] Jacques S. L., Modeling tissue optics using Monte Carlo modeling: a tutorial, Proc. Of SPIE, Vol.6854, [17] Kolinko V. G., F. F. M. Del Mul, J. Greve, A. V. Priezzhev, On refraction in Monte Carlo simulations of light transport through biological tissue, Med. Biol. Eng. Comp, Vol. 35, pp , [18] Wang L., S. L. Jacques, MCML-Monte Carlo modeling of light transport in multi-layered tissues, Computer Methods and Programs in Biomedicine, Vol.47, pp , [19] Arkin H., L. X. Xu, K. R. Holmes, Recent developments in modeling heat transfer in blood perfused tissues, IEEE Transaction on Biomedical Engineering, Vol. 41, No. 2, February [20] Chang I., Finite element analysis of hepatic radiofrequency ablation probes using temperature-dependent electrical conductivity, Biomedical Engineering Online, content 2/1/12, [21] Glenn N. T., S. Rastegar, S. L. Jacques, Finite Element Analysis of temperature controlled coagulation in laser irradiated tissue, IEEE transaction on biomedical engineering, Vol.43, No.1, January [22] Marqa M., P. Colin, P. Nevoux, S. R. Mordon, N. Betrouni, Focal laser ablation of prostate cancer: numerical simulation of temperature and damage distribution, Biomedical Engineering Online, content 10/1/45, 2011.

8 American Journal of Biomedical Engineering 2015, 5(3): [23] Olsrud J., R. Wirestam, B. R. R. Persson, K. G. Tranberg, Simplified treatment planning for interstitial laser thermotherapy by disregarding light transport: a numerical study, Lasers in Surgery and Medicine, Vol.25, pp , [24] Cvetkovic M., D. Poljak, A. Peratta, Thermal modeling of the human eye exposed to laser radiation, IEEE Transaction on Biomedical Engineering, Vol. 63, No. 6, February [25] Mohammed Y., J. F. Verhey, A finite element method model to simulate laser interstitial thermo therapy in anatomical inhomogeneous regions, Biomedical Engineering Online, content 4/1/2, [26] Mohammed Y., J. F. Verhey, A. Ludwig, K. Giese, Anatomical of a practical model for light and heat distribution using laser-induced thermotherapy near to a large vessel, Phys. Med. Biol., Vol. 48, pp , [27] Ma N., X. Gao, X. X. Zhang, Two-laser simulation model of laser-induced interstitial thermo-therapy, Laser in Medical Science, Vol.18, pp , [28] COMSOLAB Software, Sweden,