Ahmet Yasin ÇEİNTKAYA, Şaban Selim ŞEKER

Study of temperature distribution in light tissue interaction using the FEM

In this paper, the numerical modeling of interaction between light sources and living biological tissue is studied. One of the most important results of physical interactions of light with biological tissue are thermal changes, which are a consequence of absorption of light. The severity of the effects is dependent upon several factors, including exposure duration, wavelength of the beam, beam energy, tissue type, and the area that is exposed to the beam. The temperature response of tissue irradiation is governed by the extensively used Pennes bioheat equation in a 3-dimensional space. The Pennes equation is solved initially for a single-layer human skin model, thanks to which temperature distribution can be predicted using different laser light radiation burns at steady state. Furthermore, a multilayer model is used to investigate the effects of the energy of the beam and exposure duration. The numerical solutions are obtained by the finite element method (FEM), in which the geometry studied is divided into a finite element mesh. This method gives a better description of the geometry with a smaller number of nodes when compared with the finite difference time domain (FDTD) method. Using the FDTD method, which may require spatial step ranges to achieve convergence, makes the computational simulation time increase dramatically. On the contrary, less memory space and disk space with shorter run times make the FEM more advantageous than other numerical solution techniques and that is why this method has been chosen for this study. Temperature contours and penetration depths are plotted in 3-dimensional spaces. In this work it is demonstrated that by adjusting the exposure duration, wavelength of the beam, energy of the beam, and the area and type of tissue, temperature effects of light on tissue can be altered.

PDF

___

[1] Fanjul-Velez F, Romanov OG, Arce-Diego JL. Efficient 3D numerical approach for temperature prediction in laser irradiated biological tissues. Comput Biol Med 2009; 39: 810-817.
[2] Pennes HH. Analysis of tissue and arterial blood temperatures in the resting human forearm. J Appl Physiol 1948; 1: 93-122.
[3] Yoon G, Welch AJ, Motamedi M, Gemert MV. Development and application of three-dimensional light distribution model for laser irradiated tissue. IEEE J Quantum Elect 1987; 23: 1721-1733.
[4] Vo-Dinh T. Biomedical Photonics Handbook. Boca Raton, FL, USA: CRC Press, 2003.
[5] Wang L, Jacques SL, Zheng L. MCML-Monte Carlo modeling of light transport in multi-layered tissues. Comput Meth Prog Bio 1995; 47: 131-146.
[6] Kolios MC, Worthington AE, Sherar MD, Hunt JW. Experimental evaluation of two simple thermal models using transient temperature analysis. Phys Med Biol 1998; 43: 3325-3340.
[7] Verhey JF, Mohammed Y, Ludwig A, Giese K. Implementation of a practical model for light and heat distribution using laser-induced thermotherapy near to a large vessel. Phys Med Biol 2003; 48: 3595-3610.
[8] Zhou J, Zhang Y, Chen JK. An axisymmetric dual-phase-lag bioheat model for laser heating of living tissues. Int J Therm Sci 2009; 48: 1477-1485.
[9] Ng EY, Tan HM, Ooi EH. Prediction and parametric analysis of thermal profiles within heated human skin using the boundary element method. Philos T Roy Soc A 2010; 368: 655-678.
[10] Citkaya AY, Seker SS. Modeling and simulation of temperature distribution in laser-tissue interaction. In: PIERS Progress In Electromagnetics Research Symposium; 2023 March 2011; Marrakesh, Morocco. Cambridge, MA, USA: PIERS. pp. 844-847.
[11] Cheong WF, Prahl SA, Welch AJ. A review of the optical properties of biological tissues. IEEE J Quantum Elect 1990; 26: 2166-2185.
[12] Bashkatov AN, Genina EA, Kochubey VI, Tuchin VV. Optical properties of human skin, subcutaneous and mucous tissues in the wavelength range from 400 to 2000 nm. J Phys D Appl Phys 2005; 38: 2543-2555.
[13] Thompson MS. Photodynamic therapy utilizing interstitial light delivery combined with spectroscopic methods. PhD, Lund University, Sweden, 2004.
[14] Klose AD, Larsen EW. Light transport in biological tissue based on the simplified spherical harmonics equations. J Comput Phys 2006; 220: 441-470.