Numerical Modeling of Compression-Controlled Low-level Laser Probe for Increasing Photon Density in Soft Tissue
Numerical Modeling of Compression-Controlled Low-level Laser Probe for Increasing Photon Density in Soft Tissue
Journal of the Optical Society of Korea. 2011. Dec, 15(4): 321-328
Copyright ©2011, Optical Society of Korea
  • Received : July 07, 2011
  • Accepted : October 10, 2011
  • Published : December 25, 2011
Export by style
Cited by
About the Authors
Kiwoon Kwon
Taeyoon Son
Changmin Yeo
Byungjo Jung
Various methods have been investigated to increase photon density in soft tissue, an important factor in low-level laser therapy. Previously we developed a compression-controlled low-level laser probe (CCLLP) utilizing mechanical negative compression, and experimentally verified its efficacy. In this study, we used Bezier curves to numerically simulate the skin deformation and photon density variation generated by the CCLLP. In addition, we numerically modeled changes in optical coefficients due to skin deformation using a linearization technique with appropriate parameterization. The simulated results were consistent with both human in vivo and porcine ex vivo experimental results, confirming the efficacy of the CCLLP.
Ultraviolet (UV), visible, and near infrared (NIR) light have been extensively used in therapeutic and diagnostic modalities [1 - 3] . However, the depth of light penetration is still limited to a few centimeters because of the light scattering property of tissue, resulting in low photon densities.
Various methods such as increasing the power and diameter of the incident beam [3 - 5] , application of optical clearing agents [6 , 7] , and skin compression [8 - 10] have been suggested to increase the photon density in tissues. In a previous study, we developed a compression-controlled low-level laser probe (CCLLP) to increase the photon density in tissues by applying negative compression [11] . Figure 1 represents the CCLLP, which consists of a vacuum hole to supply negative compression, a hollow brace for guiding and supporting the optical fiber, and an acrylic probe body to maintain the negative compression applied to the tissue. We experimentally validated the efficacy of the CCLLP by measuring the photon density in ex vivo porcine skin samples [11] .
In another previous study, we mathematically modeled and simulated the effects of tissue compression [10] . We simulated tissue compression in two ways: (1) non-uniform skin deformation with constant optical coefficients and (2) uniform skin deformation with variable optical coefficients depending on changes in skin thickness. However, neither of these approaches accurately simulated the mechanism of the CCLLP.
Although several numerical simulations for non-flat boundaries or layers have been reported [12 , 13] , the geometry
Lager Image
The compression-controlled low-level laser probe (CCLLP) constructed with radii of (a) 1, (b) 1.5, and (c) 2 cm.
of the outer boundaries and inner layers similar to the skin affected by the CCLLP has not been described. The skin region affected by the CCLLP swells outward except for the area under the optical fiber probe. In the previous uniform compression study [10] , the photon density did not increase in the swollen areas but increased in the region of skin under the optical fiber probe. In this study, which is a fundamental generalization of the simulation reported in [10] , we numerically simulated such competing effects between the swollen and compressed regions in the CCLLP.
In three subjects, we experimentally measured the change in skin thickness due to negative compression and used the values for the numerical simulations. We represented the deformation of the skin surface associated with the CCLLP using three quadratic Bezier curves. We also numerically modeled the changes in optical coefficients with changes in skin thickness using a linearization method. We obtained simulation results as functions of the modeling parameters, radius of the CCLLP, and negative compression, and compared these with results from an ex vivo experiment using a scaling factor.
The photon density in tissues was modeled using the radiative transfer equation or its diffusion approximation [14] . The diffusion approximation in the frequency domain was represented by the following elliptic partial differential equation in the region of interest Ω :
Lager Image
where r0 is some point on the boundary of Ω . For the notations used in the equations, see the Nomenclature section before the Reference section at the end of this paper.
To solve (1) using the finite element method (FEM), the geometrical domain Ω was assumed to be bounded, and appropriate boundary conditions were assigned to the boundary Ω . Let us assume that Ω = [- w , w ] × [-d, 0]. That is to say, the region of interest is a rectangular region in the two-dimensional xz plane modeling a few skin layers with depth d , deposited in the direction of the z-axis. Then Ω is composed of four straight lines:
  • ∂topΩ =[-w, w]×{0}: the top of the skin contacting the air
  • ∂bottomΩ =[-w, w]×{-d}: the bottom of the skin region of interest
  • ∂leftΩ ={-w}×{-d,0}: the left of the skin region of interest
  • ∂rightΩ ={-w}×{d,0}: the right of the skin region of interest
Note that only top Ω contacts air, while the other boundaries contact skin. Thus, we imposed the Neumann boundary condition on left Ω , right Ω , and bottom Ω and the Robin boundary condition on top Ω to reflect differences in refractive index mismatching [14] . Therefore, we obtained the following boundary value problem for modeling light propagation in skin:
Lager Image
Let us divide Ω into small triangles or rectangles, which are called finite elements, and { ui } i =1,… Nb be the corresponding finite element bases spanning piecewise-linear or -bilinear function spaces. The finite element formulation for (1) with the boundary condition (2) using the finite element bases is then given as follows:
Lager Image
The finite element solution for (3) was found using COMSOL 3.4 (COMSOL AB, Stockholm, Sweden). In the simulation, the photon density, or light fluence rate, was measured at the bottom ∂ bottom Ω of the domain of interest and at a distance 7.5 cm away. However, we could neglect the 7.5 cm distance in the air for the photon density measurement since light is mostly absorbed and scattered in the tissue rather than in the air.
The optical properties of human tissues have been reported in several studies [3 , 15] . Based on these studies, a five-layer skin model was used for a 633 nm light [3] .
- 2.1. Skin Deformation by Negative Compression
For simplicity, we modeled the domain of interest Ω in two dimensions by neglecting the y-axis. When there is no compression, Ω is given as follows:
Lager Image
where Ω i and di are the domain of interest and the thickness of the i -th skin layer, respectively, dx is the dimension of the skin layer along the x -axis, and t 5 = d . The thicknesses of layers di ( i = 1,2,3,4,5) are presented in Table 1 with their optical coefficients. The overall five-layer structure is depicted in Figure 2 .
We then formulated the layers Ω i ( i = 1,2,3,4,5) with negative compression. Let 2 K and 2 O be the length of the region contacting the entire acrylic probe body of the CCLLP and
Five-layer skin model for 633 nm light used to numerically model the photon density profile in tissue
Lager Image
Five-layer skin model for 633 nm light used to numerically model the photon density profile in tissue
Lager Image
The comparison of the five-layer skin structure before (black line) and after deformation (red line on the surfaces of the first and second layers).
the length of the region contacting the guide for the optical fiber, respectively. Let C be the midpoint between K and O . If negative compression is applied throughout the regions [ -K , -O ] and [ O , K ] on the skin surface, these regions become swollen by the maximum height difference shown in Figure 3 for the applied negative compression. On the other hand, in the region [- O , O ], the skin surface should be retracted between the points (- O , 0) and ( O , 0) because of the swollen surfaces surrounding it. For more detailed evaluations, elastic or visco-elastic analysis could be used [16 , 17] ; however, in this study, we used three quadratic Bezier curves for the regions [- K , - O ], [- O , O ], and [ O , K ]. The data shown in Figure 4 were used as the maxima of the two Bezier curves for [- K , - O ] and [ O , K ], and the Bezier curve for [- O , O ] was constructed to join the three curves without any discontinuity of tangential lines. Let li be the upper surface of i -th layer ( i = 1,2,3,4,5), which is a straight line before the application of negative compression. For this purpose, control points were chosen as follows:
Lager Image
Lager Image
In vivo skin deformation due to negative compression at (a) 0, (b) 5, (c) 10, (d) 15, and (e) 20 kPa applied to the skin surface of a subject using the CCLLP with a radius of 1 cm.
Lager Image
The changes in skin thickness as a function of negative compression and CCLLP radius (1, 1.5, and 2 cm).
where C was chosen as K + O/2 to cause the two Bezier curves for [- K , - O ] and [ O , K ] to be symmetric; S = 2 O/ ( K - O ) was chosen so that the three Bezier curves joined without any discontinuity of tangent lines, especially at - O and O ; and si(i = 1,2,3,4,5) is the maximum skin surface change for i -th layer. We chose w = 3.5 cm, K = 2 cm, and O = 0.2 cm, based on the sizes of the skin sample and CCLLP.
Using the definition of a quadratic Bezier curve, we obtained the following detailed representations for li, i = 1,2,3,4,5:
Lager Image
Therefore, the thicknesess of the i layers are as follows:
Lager Image
Assuming the thickness changes for the third, fourth, and fifth layers are relatively small and the change for second layer is half of the change in the first layer, we let
s 3 = s 4 = s 5 = 0 and
Lager Image
- 2.2. Change in Optical Coefficients with Negative Compression
In the CCLLP, negative compression is applied to the skin surrounding the central area [–O,–K ]×{0}⋃[ K,O ]×{0} in contact with the probe, and positive compression is applied at the center [? K,K ]×{0} due to skin tension. The change in optical coefficients with uniform and partial compression of the skin was previously studied and reported [10] . Let
Lager Image
μ a , μ s '
be the optical coefficients for the i-th layer ( i = 1,2) when the thickness of the layer is di and
Lager Image
respectively. Using an analysis similar to that in [10] , we assumed that the relationship between the thickness and the absorption coefficient is:
Lager Image
The coefficient α is the optical modeling coefficient between 0 and 1, which related the reduction of skin penetration depth to the increase of optical coefficients. The more fundamental meaning of the coefficient is discussed in the “Discussion” section.
Equation (8) is linearly approximated as follows:
Lager Image
Similarly, we obtained the following relationship between the thickness and the reduced scattering coefficient:
Lager Image
Given that the diffusion coefficient is defined as a constant multiplied by the reciprocal of the sum of the absorption and reduced scattering coefficients, using (9) and (10) we obtained:
Lager Image
In (9), (10), and (11), we used first-order Taylor series approximations.
Using (7) and (11), we obtained an explicit representation of the diffusion coefficient as follows:
Lager Image
We chose the value for s1 based on the data shown in Figure 4 , and si(i = 1,2,3,4,5) were chosen such that
Lager Image
decreased with the ratio of the depth from the surface to the entire depth of the skin. Thus, (12) becomes:
Lager Image
and (9) becomes:
Lager Image
Figure 1 shows CCLLPs with radii of 1, 1.5, and 2 cm. In Figure 3 , Skin deformation by negative compression was investigated in the left arms of three subjects. Figure 4 shows the change in skin thickness averaged across the three subjects during negative compression as a function of the radius of the CCLLP. As described in the section II.1, the change in skin thickness was modeled using three Bezier curves with continuous tangential lines.
For the FEM analysis of the CCLLP, a finite element mesh was applied to the modeled geometry of the skin. We used meshes shown in Figure 5 for the skin deformation geometry to obtain the photon density distribution by FEM. The mesh is finer in the upper layers and coarser in the lower layers because skin deformation is expected to be less pronounced in the lower layers.
The optical coefficients for the five layers used in modeling the skin structure are shown in Table 1 [3] . The resulting photon density distributions at a depth of 0.2 cm below the skin surface are presented in Figure 6 as functions of negative compression, the thickness-coefficient
Lager Image
Mesh grids used for finite element analysis of (a) the entire skin sample and (b) the region around the optical fiber probe with a negative compression of 5 kPa and a CCLLP radius of 2 cm.
relation parameter α which is represented in (8), and the radius of the CCLLP.
In our previous study [11] , we investigated the efficacy of CCLLP by measuring the photon density distribution in ex vivo porcine samples. In Figure 7 , the photon density ratio, which is defined as the ratio of the photon density distribution between a given negative compression and no compression, was compared to the experimental data for a CCLLP with a radius of 1 cm. To compare the results of this study with the ex vivo data, the maximum change in skin thickness in vivo for human subjects was halved in the numerical simulation. The experimental and simulation results for various values of α (0.5,0.6,0.7, and 0.8) was shown in Figure 7 .
The photon density in the tissue due to negative compression was influenced by two successive and competing physical phenomena: (1) the change in skin thickness due to compression and (2) the change in optical coefficients due to the change in skin thickness. In previous studies, the correlation between the change in skin thickness and negative compression was mathematically modeled using elastic equations [18 , 19] . In this study, the correlation was experimentally obtained from three subjects, and the data ( Figure 4 ) was used to mathematically simulate the mechanism of the CCLLP. The changes in skin thickness for all three subjects were approximately proportional to the logarithm of the negative compression, which corresponds well with the results of previous studies [18 , 19] .
Skin surface deformation by the CCLLP was experi
Lager Image
Photon density distribution at the base of the skin sample as a function of the (a) negative compression (α = 0.1, radius = 1 cm), (b) parameter α (radius = 1 cm, negative compression = 5 kPa), and (c) CCLLP radius (α = 0.1, negative compression = 5 kPa).
Lager Image
The comparison of experimental photon density ratios [11] and modeling results with respect to the various values of the optical coefficient modeling parameter α . The photon density ratio is defined as the ratio of the photon density distribution with negative compression (0, 5, 10, 15, or 20 kPa) to that without compression (0 kPa)
mentally measured as a function of negative compression ( Figure 3 ). The region of skin contacting the optical fiber probe was compressed for two reasons: (1) this region of the skin was pressed down and thus could not swell up, and (2) the skin in the contacted region was pushed aside to the surrounding region, which swells outward. In other words, the skin inside the CCLLP, except for that in contact with the optical fiber probe, swells outward, whereas the contacted region is compressed inward. Thus, the skin surface can be modeled by: (1) two straight lines representing the two flat areas of skin outside the CCLLP, (2) two concave Bezier curves representing the raised skin region inside the CCLLP not contacting the optical fiber probe, and (3) one convex Bezier curve representing the compressed skin region contacting the optical fiber probe. The two concave Bezier curves and the convex Bezier curve were constructed to satisfy the C1 assumption at the two points where the curves meet, representing the boundary of the optical fiber probe. The C1 assumption states that the tangential lines representing the skin surface, as well as the skin surface itself, are continuously connected across the optical fiber probe inside the CCLLP. In other words, the skin surface deformation due to negative compression should be smooth and not have any discontinuities along the skin surface. The resulting skin surface deformation was shown in Figure 5 .
The optical coefficients are known to be inversely proportional to skin penetration depth [8 , 9] . However, the relationship has been experimentally investigated with only a limited number of samples and is not well established for human or animal skin samples. Therefore, the relationship was mathematically modeled [10] when uniform compression was applied to the skin surface, corresponding to α = 0.5 in (8). Analysis of the changes in the optical coefficients with skin thickness was generalized from a uniform geometry [18] to a more complex geometry ( Figure 5 ). Detailed changes in the optical coefficients with the generalized parameter α were computed using (13) and (14).
The skin thickness at the center decreased because of negative compression, resulting in increased photon density. At the same time, by the result of [8 , 9 , 10] , we assumed that the optical coefficients (including the absorption coefficient (8),(9), the reduced scattering coefficient (10), and the diffusion coefficient (11)) per unit skin penetration depth increased through the optical modeling coefficient α as the skin thickness decreased, and thus interfered with the increase in photon density. Previous experimental [8 , 9] and simulation [10] results have shown that the increase in photon density due to the reduction in skin thickness is greater than the decrease due to changes in optical coefficients. However, this study modeled a more complicated phenomenon than the uniform compression case. In the skin deformation model, the skin regions inside the CCLLP and those contacting the optical fiber probe are raised and compressed, respectively. Therefore, the effects of skin deformation on photon density in these two regions are in opposite directions. A previous study using ex vivo porcine samples [11] presented the efficacy of the CCLLP in enhancing photon density. These results were verified by numerical simulation results ( Figure 6 ) which used data on in vivo changes in human skin thickness. From these results, we conclude that the decrease in skin thickness directly under the optical fiber probe is more important than the increase in skin thickness in the surrounding areas. In addition, the negative compression, value for α , and radius of the CCLLP in the present results were all larger than in the previous study.
The peak intensities of photons in the ex vivo porcine samples [11] were compared with the simulated results. Changes in skin thickness were not obtained in the previous study because of the instability of the physical state of the ex vivo porcine skin samples. Therefore, skin deformation was measured in the arms of three human subjects, and changes in the skin thickness were halved to enable comparisons with the ex vivo porcine skin samples.
The CCLLP was developed to increase the photon density in tissues by applying negative compression. Skin deformation due to the CCLLP was numerically modeled using three Bezier curves. The finite element simulation used experimental data on changes in skin thickness from three human subjects, skin deformation modeling, and changes in optical coefficients as a function of the optical modeling coefficient α . The photon density increased when the negative compression and radius of the CCLLP increased. The simulated results at various values of α were similar to the ex vivo experimental results. Although we obtained favorable result using the optical modeling coefficient, we have not experimental verification about the optimal choice for the coefficient, which will be our next study subject.
Φ : The photon density profile [W/cm 2 ]
q : The light source function [W/cm 3 ]
R : The Robin function, which is the photon density profile with respect to the Dirac delta source [W/cm 2 ]
μ a : The absorption coefficient [cm -1 ]
μ s : The scattering coefficient [cm -1 ]
g : The anisotropy factor
μ s ' : The reduced scattering coefficient [cm -1 ] [μ s ' = (1-g) μ s]
κ : The diffusion coefficient [cm] [κ = (3(μ a μ s '))-1]
n : The refractive index
a : The reflection coefficient on the surface
c : The speed of light in a vacuum [cm/sec] (c = 2.99×10 10 cm/sec)
ω : The angular modulation frequency of the source intensity [sec -1 ]
Ω : The region of interest (skin structure)
∂ Ω : The boundary of Ω (∂ Ω = ∂ top Ω ∂ left Ω ∂ right Ω + ∂ bottom Ω )
top Ω : The top boundary of Ω
bottom Ω : The bottom boundary of Ω
left Ω : The left boundary of Ω
right Ω : The right boundary of Ω
N b : The number of finite element basis functions
This study was supported by a grant from the Korea Healthcare Technology Research and Development Project, Ministry for Health, Welfare and Family Affairs, Republic of Korea (B090033), and the Dongguk University Research Fund of 2010, Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Education, Science and Technology (2010-0004047).
Jobsis F. F. (1977) “Noninvasive, infrared monitoring of cerebral and myocardial oxygen sufficiency and circulatory parameters,” Science 198 1264 - 1267
Gratton G. , Maier J. S. , Fabiani M. , Mantulin W. W. , Gratton E. (1994) “Feasibility of intracranial near-infrared optical scanning,” Psychophysiology 31 211 - 215
Tuchin V. V. 2000 Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis SPIE Press WA, USA,
Grängesberg, Sweden J. , Hode L. 2002 Laser Therapy-clinical Practice and Scientific Background Prima Books Grängesberg, Sweden
Carroll L. , Humphreys T. R. (2006) “LASER-tissue interactions,” Clin. Dermatol. 24 2 - 7
Tuchin V. V. , Xu X. , Wang R. K. (2002) “Dynamic optical coherence tomography in studies of optical clearing, sedimentation, and aggregation of immersed blood,” Appl. Opt. 41 258 - 271
Kang H. , Son T. , Yoon J. , Kwon K. , Nelson J. S. , Jung B. (2008) “Evaluation of laser beam profile in soft tissue due to compression, glycerol, and micro-needling,” Lasers Surg. Med. 40 570 - 575
Chan E. K. , Sorg B. , Protsenko D. , O'Neil M. , Motamedi M. , Welch A. J. (1996) “Effects of compression on soft tissue optical properties,” IEEE J. Select. Topics Quantum Electron. 2 943 - 943
Shangguan H. , Prahl S. A. , Jacques S. L. , Casperson L. W. , Gregory K. W. (1998) Pressure effects on soft tissues monitored by changes in tissue optical properties,” Proc. SPIE 3254 366 - 371
Kwon K. , Son T. , Lee K. J. , Jung B. (2009) “Enhancement of light propagation depth in skin: cross-validation of mathematical modeling methods,” Lasers Med. Sci. 24 605 - 615
Yeo C. , Park J. , Son T. , Lee Y. , Jung B. (2009) “A pressure applied low-level laser probe to enhance laser photon density in soft tissue,” J. Biomed. Eng. Res. 30 18 - 22
Melinger I. V. , Matcher S. J (2001) “Modelling the sampling volume for skin blood oxygenation measurements,” Med. Biol. Eng. Comput. 39 44 - 50
Okada E. , Firbank M. , Schweiger M. , Arridge S. R. , Cope M. , Delpy D. T. (1997) “Theorecitcal and experimental investigation of near-infrared light propagation in a model of the adult head,” Appl. Opt. 36 21 - 31
Arridge S. R. (1999) “Optical tomography in medical imaging,” Inverse. Probl. 15 R41 - 93
Van Gemert M. J. C. , Jacques S. L. , Sterenborg H. J. C. M. , Star W. M. (1989) “Skin optics,” IEEE T. Bio-Med. Eng. 36 1146 - 1154
Childer C. C , McCoy C. W. , Nigg H. N. , Stansly P. A. , Roger M. E. Florida Pest Management Guide: Rust Mites, Spider Mites, and Other Phytophagous Mites (University of Florida Cooperative Extension Service Institute of Food and Agricultural Sciences Gainesville, FL, USA
Franco M. P. , Mulder M. , Gilman R. H. , Smits H. L. (2007) “Human brucellosis,” Lancent Infect. Dis. 7 775 - 786
Childers M. A. , Franco W. , Nelson J. S. , Aguilar G. (2007) “Laser surgery of port wine stains using local vacuum prezssure: changes in skin morphology and optical properties (Part I),” Lasers Surg. Med. 39 108 - 117
Franco W. , Childers M. , Nelson J. S. , Aguilar G. (2007) “Laser surgery of port wine stains using local vacuum [corrected] pressure: changes in calculated energy deposition (Part II),” Lasers Surg. Med. 39 118 - 127