Next Article in Journal
Combined Effect of Spirulina Platensis and Punica Granatum Peel Extacts: Phytochemical Content and Antiphytophatogenic Activity
Next Article in Special Issue
Measurement of Structural Loads Using a Novel MEMS Extrinsic Fabry–Perot Strain Sensor
Previous Article in Journal
Fuzzy Rules to Help Predict Rains and Temperatures in a Brazilian Capital State Based on Data Collected from Satellites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Study on Pressure Damage Based on Clinical Purpura during the Laser Irradiation of Port Wine Stains with Real Complex Vessels

1
State Key Laboratory of Multiphase Flow in Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China
2
State-Province Joint Engineering Lab of Fluid Transmission System Technology, Zhejiang Sci-Tech University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Submission received: 18 November 2019 / Revised: 7 December 2019 / Accepted: 9 December 2019 / Published: 13 December 2019
(This article belongs to the Special Issue Biomedical Optics and Engineering)

Abstract

:
Port wine stains (PWSs) are congenital dermal vascular lesions composed of a hyperdilated vasculature. Purpura represented by local hemorrhage from water vaporization in blood during laser therapy of PWS is typically considered a clinical feedback, but with a low cure rate. In this study, light propagation and heat deposition in skin and PWSs is simulated by a tetrahedron-based Monte Carlo method fitted to curved bio-tissues. A curvature-corrected pressure damage model was established to accurately evaluate the relationship between purpura-bleeding area (rate) and laser therapy strategy for real complex vessels. Results showed that the standard deviation of Gaussian curvature of the vessel wall has negative relation with the fluence threshold of vessel rupture, but has positive relation with the effective laser fluence of vessel damage. This finding indicated the probable reason for the poor treatment of PWS, that is, considering purpura formation as a treatment end point (TEP) only leads to partial removal of vascular lesions. Instead, appropriate purpura area ratio with marked effects or rehabilitation should be adopted as TEP. The quantitative correlation between the fluence of a pulsed dye laser and the characteristics of vascular lesions can provide personalized and precise guidance for clinical treatments.

1. Introduction

Port wine stains (PWSs) are congenital vascular birthmarks characterized by overabundant dermal capillaries or venules, and these marks occur in approximately 0.3% of newborns [1]. The effective treatment strategy of PWSs focuses on the laser irradiation of the affected portions of the skin [2,3]. In accordance with selective photothermolysis (SP) theory [4], laser irradiation can irreversibly induce thermal injury to abnormal blood vessels while minimally damaging the surrounding skin tissue.
A pulsed dye laser (PDL) with a wavelength of 585 or 595 nm is considered the golden standard for the treatment of PWS lesions [2,5]. In comparison with a low absorption efficiency and an unapparent response of the regions exposed to lasers at the near-infrared band, the strong absorption of PDL energy by hemoglobin can easily result in water vaporization in blood and local hemorrhage, which is clinically manifested as purpura [2,6,7], a real-time feedback of SP. According to more than 100 experiments, a laser-irradiated region remains perfused in the absence of purpura of the irradiated vessels. In clinical practice, lesions with purpura formation undergo discoloration 1 month after treatment, whereas lesions in a nonpurpura region remain dark. Therefore, purpura is typically considered a qualitative clinical treatment end point (TEP) and prognostic indicator for lesion removal [8].
The extent of purpura formation has not been regarded as an accurate indicator of treatment outcome. A too low or too high therapeutic dose can lead to vascular recanalization or scarring, and clinical therapy gradually becomes limited in improving the cure rate (only approximately 50%) of PWS nevus [9]. Purpura is characterized by blood spilling out of the skin and the mucosa after vascular rupture. A model involving “purpura-bleeding area (rate)–realistic vessel geometry–laser therapy strategy” should be established to discover the mechanism of vessel damage represented by purpura and improve the curative effect [10]. Nevertheless, minimal attention has been paid to this issue.
Blood coagulation and vessel rupture are regarded as thermal damages and modeled by the damage integral Ω, which increases exponentially with tissue temperature [11,12,13,14]. After the temperature field is solved with the bioheat conduction equation, Ω can be obtained by calculating the integral related to temperature, and blood coagulation can be represented as high-Ω regions [15]. As a subsequent effect of blood coagulation, complete vessel constriction has also been modeled by Ω, although its reasonable damage threshold should be clarified [8]. However, this single-parameter damage model considering thermal blood coagulation is inapplicable to vessel rupture, which is marked by purpura and has been improperly regarded as a type of thermal damage. In animal experiments utilizing a PDL, vessel wall rupture is accompanied with large bubble formation from the change in the water phase in blood [8,16]. Vessel rupture is caused by excessive pressure on a blood vessel wall after vaporization, not by heat from blood cells to the blood vessel wall. Pressure produced by gasification probably damages the blood vessel wall before heat transfers from blood cells to the blood vessel wall. Therefore, vessel rupture should be modeled by using a given pressure threshold and not by heat damage Ω.
Zhang et al. [14] introduced a heat transfer equation with a pressure-dependent vaporization factor to predict temperature distribution, including water phase change in laser surgery in accordance with SP theory. However, the luminal pressure of the blood vessel in this model is only used to correct the temperature field and Ω rather than be regarded as an independent parameter to evaluate the vessel rupture potential, which cannot be predicted effectively in this manner. Our group [7] proposed a rupture potential index based on the wall pressure (RPIP) to estimate the vessel rupture. A double-parameter criterion comprising RPIP and temperature integral (Ω) is also suggested to investigate the competition between the pressure damage (vessel rupture) and thermal damage (vasoconstriction) of the irradiated vessels. Our results show that the damage characteristics of the PDL and the infrared laser are fundamentally distinguished. In the PDL, vessels constantly reach the pressure threshold prior to the heat damage threshold, and the vessel undergoes rupture instead of complete constriction. With infrared laser, uniform energy is deposited in the vessel, thereby inducing much higher probability of constriction than hemorrhage. The superiority of PDL is revealed by verifying the rupture as a more severe damage type than the constriction because of the more drastic pressure variation at the vessel wall.
In clinical settings, the complex structure of the PWS vasculature results in an unpredictable rupture position and bleeding volume, leading to erroneous diagnosis by physicians. A numerical study on real blood vessels, such as aorta and celiac artery, has shown that the stress on the vessel wall often reaches the extremum at the local saddle point [17]. This finding indicates that the bending degree of the vessel surface significantly affects the pressure distribution of the vascular wall. Bucking and tortuosity occur frequently in the patient-specific PWS vessel structure, thereby changing the blood hydrodynamics and seriously affecting the rupture [18]. Available skin models usually adopt ideal cylinders as pathological vessels, but they do not represent real clinical lesions. Heat deposition in individual vessels by light absorption in the blood can be significantly reduced by competition from a branched or surrounding vasculature, which is also not considered by a simplified geometric model [19]. Consequently, an anatomically realistic geometric model of the vasculature should be established to accurately simulate the force on vessel walls, and the curvature of these walls should be investigated to quantitatively characterize the relationship between vascular morphology and pressure damage.
Before a theoretical study is conducted, a real vascular lesion should be geometrically reconstructed. Computer-reconstructed biopsy shows the detailed structure of vascular lesions, but it causes unnecessary injury to patients’ skin [20]. Noninvasive skin imaging technologies, such as optical coherence tomography and laser speckle imaging, are reliable in distinguishing normal or diseased tissues on a large scale but have a low spatial resolution in 3D imaging of thin tissue layers, such as skin tissues [21]. As a type of in vivo imaging, a dorsal skin chamber (DSC) model [16,22,23,24,25,26] is effective in investigating the thermal response of blood vessels. The geometric feature of blood vessels in the DSC model is nearly parallel, thereby allowing 3D reconstruction from 2D blood vessel images to reveal the mechanism of real vascular injuries, especially vascular rupture and hemorrhage. Mu et al. [27] performed a preliminary 3D reconstruction of real blood vessels from the image processing of a DSC model through microscopic observation. However, a vessel rupture based on vessel wall curvature is not considered, leading to an unclear relationship between the bleeding area of blood vessels (purpura) and the parameters of lasers.
In the present study, a curvature-corrected pressure damage (CCPD) model was proposed to accurately predict the rupture of a real radiant vessel through 3D vascular reconstruction and establish the relationship between purpura-bleeding area (rate) and laser therapy strategy. The model was validated by comparing its results with the experimental results of blood vessels in the DSC model with different curvatures. The effects of vessel curvature, burial depth, and epidermal melanin content (i.e., human skin color) on the parameters of the 595 nm PDL treatment for PWSs were systematically studied to establish a quantitative relationship of the parameters of PDL therapy with the burial depth and curvature of vascular lesions and the melanin content of the epidermis. This study provided a guide on the clinical treatment of PWS.

2. Materials and Methods

2.1. Problem Description and Basic Assumptions

Human skin consists of two main layers, i.e., the superficial epidermal layer and the deep dermal layer. PWSs result from the malformation of a significant amount of blood vessels buried within the dermis. In the treatment of vascular lesions, laser energy is mainly absorbed by hemoglobin in blood. However, epidermal melanin is another chromophore for light absorption, and its content influences the risk of epidermal injury. According to SP theory [4], laser irradiation should maximally induce an irreversible injury to abnormal blood vessels while minimizing the risk of injury in the surrounding dermis and epidermis. The major assumptions in this work are as follows:
(1)
The epidermis of a rat was removed in the DSC model for microscopic observation [12,13]. For convenience, the epidermal layer was excluded from the simulation to validate the theoretical model by the in vivo results from the DSC model. Afterward, the whole skin model comprising the epidermis and the dermis was adopted to study the effects of epidermal melanin content and vessel depth on laser fluence. The surfaces of the epidermis and the dermis were considered planar.
(2)
The tetrahedron-based Monte Carlo (TMC) method was employed to simulate light propagation and energy deposition because a tetrahedron mesh is body fitted for a complex vessel shape [28]. The optical properties of different tissues are assumed to be constant and uniform within a given component. The following parameters were disregarded during the short laser pulse duration: changes in the density, heat capacity, and thermal conductivity of skin tissues; metabolic heat generation; and biochemical reactions.
(3)
The epidermis and the dermis were treated as solid phases because of the low water content [14], and only the energy equation was involved. Thermal damage represented by temperature was used to characterize the safety of the epidermis and the dermis.
(4)
The purpura area was assumed to be equal to the bleeding area in the DSC model and directly proportional to the volume fraction of blood vessels in the lesion under the same treatment parameters in accordance with the approximate relationship between the clinical treatment and the animal experiment. The purpura image in clinical settings was analyzed to provide the appropriate area ratio of purpura for setting the laser parameters in simulation.

2.2. Governing Equations

The complete two-layer skin model contained the epidermis, the dermis, and the blood vessels (Figure 1). A 1.5 mm × 1.5 mm × 1.0 mm computational domain of the skin with a reconstructed blood vessel buried at a depth of 0.5 mm was used. The thicknesses of the epidermis and the dermis were set respectively as He = 60 μm and Hd = 940 μm based on statistical average values [15]. As mentioned in Section 2.1, in validating the damage model, the reconstructed blood vessels in accordance with the DSC model were embedded in the skin without the epidermis.
In the present work, light propagation and energy deposition in the skin tissue were simulated by using our in-house TMC code [28] through which the body-fitted material grid of the tetrahedra was generated by assigning a value representing a tissue type. No interface tetrahedral cells existed, thereby avoiding the error of photon reflection/refraction in the voxel-based Monte Carlo method [29]. Photon propagations in TMC were expressed as probability distributions that describe the step size of photon movement between sites of photon–tissue interaction. Once a photon has reached an interaction site, a fraction of the photon weight was absorbed in the local tetrahedron. The “random walk” of photon could be simulated repeatedly until a sufficient quantity of photons was absorbed in tetrahedral cells to establish reliable trends of energy deposition. In this work, the laser beam was a flat top profile with a diameter of 1.0 mm. 108 photons were launched at z = 0 μm with their x and y determined by the random number restricted in the range of the beam. The directional cosines (μx, μy, μz) were initialized to (0, 0, 1) and the photon weight is initialized to 1.
Energy deposition in different cells determines the heat source and thus strongly influences mass and heat flow in tissues. The following numerical formulation was implemented in the finite-volume-based CFD package FLUENT 17.0 (Ansys Inc., Canonsburg, U.S.A.) to simulate blood flow and heat transfer in vessels. The tetrahedral grid used in CFD computation was the same as that in TMC, and the energy deposition computed by TMC was treated as heat source to the governing equations. Considering that blood is an incompressible Newtonian fluid, we can write a single set of governing equations, including continuity, momentum, and energy equations in the portion of the vessel, as follows:
ρ t + ( ρ u ) = 0 ,
( ρ u ) t + ( ρ u u ) = p + [ μ ( u + u T ) ] + ρ g + F ,
( ρ c p T ) t + ( ρ c p u   T ) = ( k   T ) + S L + S E ,
where F is the surface tension, SL is the energy source computed by TMC through laser pulse, and SE is the energy source term related to evaporation or condensation.
By contrast, the epidermis and the dermis were treated as a solid phase, and only heat transfer caused by laser heating was involved:
( ρ c p T ) t + ( ρ c p u   T ) = ( k   T ) + S L .
The liquid–vapor mass transfer caused by laser heating (evaporation and condensation) is governed by the vapor transport equation [7], and the inbuilt VOF module in FLUENT 17.0 was adopted to capture the interface when bubbles formed in the blood.

2.3. Boundary Conditions

In laser dermatologic surgery, cryogen spray cooling (CSC) is applied to avoid unwanted damage from epidermal burning due to the melanin absorption of light energy. In this study, the refrigerant R134a was sprayed to the epidermis before laser pulses were administered to protect the epidermis. The spray duration was tc = 100 ms, and the spray distance was 30 mm [30]. After CSC was conducted (t > tc), the natural convection of air was considered the boundary condition. The standard convection condition was used to describe the cooling effect of CSC on the skin surface (z = 0):
k e T z | z = 0 = { h c ( T T c ) ,   t t c h ( r * , τ ) ( T T a ) ,   t > t c ,
where ke is the thermal conductivity of the epidermis, Ta = 25 °C is the ambient temperature, and ha = 10 W·m−2·K−1 is the convective heat transfer coefficient of air. Tc of −48.5 °C [30] is the temperature of the refrigerant liquid film. The dimensionless surface heat transfer coefficient h*(r*, τ) based on a large number of experimental data was expressed as follows [30]:
h ( r , τ ) = { h o ( τ ) 0 r 0.4 5 ( 1 r ) 3 h o ( τ )   0.4 < r 1.0 τ 1.0 h ( r , τ ) = { h o ( τ ) 0 r 0.2 ( τ + 1 ) h o ( τ ) + ( 5 r τ 1 ) [ 0.09 ( τ 1 ) h o ( τ ) ] ( 4 τ ) 0.2 ( τ + 1 ) < r 1.0     1.0 < τ < 4.0 h ( r , τ ) = h o ( τ ) 0 r 1.0 τ 4.0 ,
where nondimensional time τ = t/tmax, and nondimensional space r* = r/rspray; the expression of the dimensionless convective heat transfer coefficient at the center of refrigerant h0*(τ) is presented as follows [30]:
h o ( τ ) = { τ τ 1.0 1.0 0.35 ( τ 1 ) 1.0 < τ 3.0 0.3 0.02 ( τ 3 ) 3.0 < τ 8.0 0.2 0.0125 ( τ 8 ) τ > 8.0 .
Adiabatic conditions were adopted for the side and bottom boundaries. The initial skin temperature was assumed to be 35 °C, and the inlet velocity of the blood was set as 1 × 10−3 m/s [31].

2.4. CCPD Model

The morphological characteristics of a real curved blood vessel can be represented by Gaussian curvature. The two eigenvalues k1 and k2 of Weingarten transformation at any point on a surface are called the principal curvatures of surface S at that point [32]. Gaussian curvature KG is defined as the product of the principle curvatures [32] as follows:
K G = k 1 k 2 .
In this work, the vessel wall was discretized by a triangular mesh. Meyer et al. [32] obtained the discrete operator of Gaussian curvature by using the geometric parameters of the one-ring neighborhood mesh:
K G ( x i ) = ( 2 π j = 1 n θ j ) / A v o r ,
where j is the interior angle of the jth triangle at the vertex xi, n is the number of triangles around the vertex, and Avor is the area of Voronoi region [33,34] in the one-ring neighborhood of the vertex (Figure 2):
A v o r = 1 8 j N 1 ( i ) ( cot α j + cot β j ) x i x j ,
where N1(i) is the point set of all vertices (except i itself) in the one-ring neighborhood of vertex i, and αj and βj are the two angles opposite to the edge in the two triangles sharing the edge (xi, xj) as depicted in Figure 2.
The computation error of Gaussian curvature could be controlled to less than 1% by using this algorithm [32].
After the Gaussian curvature was precisely calculated, virtual pressure PGa caused by Gaussian curvature was introduced to consider the effect of the anatomical geometry of blood vessels on rupture, and this parameter was added to the wall pressure P:
P = P + P G a
where P’ is the KG-corrected pressure on the vessel wall. Statistical analysis [35] showed that KG on the vessel wall follows a negative linear relationship with the additional wall stress caused by curvature (Pearson’s r = −0.437). The pressure on the vessel wall at the moment of laser irradiation was considered proportional to the stress according to our previous work [7]. Therefore, PGa under laser irradiation can be assumed to be linear with the Gaussian curvature on the vessel wall:
P G a = A K G + B
where A and B are constants.
In this study, an ideal vessel had a cylindrical surface with a Gaussian curvature of zero. Thus, the intercept B should be zero when the cylindrical surface was considered a standard without additional pressure. The Gaussian curvatures of convex and concave surfaces were positive and negative, respectively, and the pressure capacity of the concave surface was less than that of the convex surface. Therefore, the expression of the virtual pressure caused by the curvature could be written as:
P G a = P / P w , m a x K G
where Pw,max is the maximum wall pressure of blood vessels. The model based on P’ is named as the CCPD model, which differs from the pressure damage (PD) model based on P [7]. The pressure threshold to vessel damage was chosen as 2500 Pa through analogy analysis on pressure thresholds of rupture on different vessels in our previous work [7]. The rupture area of the blood vessel (Sbr) could be determined by the area with P’ higher than 2500 Pa. The validation of the CCPD model through rupture area under a given threshold of rupture pressure is discussed in the next section.
The region of purpura or bleeding had a recognizable contour that could be captured by Image Pro Plus. Koster et al. [36] and Tanghetti et al. [37] conducted an immediate histopathological analysis of patients with PWS after laser treatment and found that Sbr follows a linear relationship with the area of purpura (Spu). Therefore, the functional relationship between Sbr and Spu could be obtained by linearly fitting between the computation and the animal experiment. The linear fit according to our damage model is also discussed in the next section.

3. Validation of CCPD Model

3.1. 3D Reconstruction of Real Blood Vessels

Considering the approximate parallelism of the blood vessel with the observation plane in the DSC model, the blood vessel axis was regarded on a horizontal plane. The dorsal blood vessels of male Sprague–Dawley rats weighing 100–120 g were used as a target to be irradiated. To make the observation chamber, the dorsal skin was lifted and fixed between two symmetrical titanium frames containing a central annular window (Ø = 15 mm). The entire epidermal, dermal, and subcutaneous connective tissues of one side of the skin in the central annular window were carefully removed to expose the vascular structure (Figure 3). The fascia on the exposed site was carefully removed to improve the imaging. The rats were placed on a specially engineered microscope translator stage containing three pins that interlock with holes in the implanted optical chamber so as to minimize motional disturbances by the breathing of the anesthetized rats during imaging. Proper targeted blood vessel was oriented by mounting the xyz-translator stage. Then, the skin images were captured with a CCD camera (DP71, Olympus Corp., Tokyo, Japan) connected to the microscope (10 × objective, NA = 0.30, BX51, Olympus Corp.) at a resolution of 1360 × 1024. The images were simply processed and displayed by cellSens standard software. All of the procedures involving the animal experiments were approved by the Institutional Animal Care and Use Committees of Xi’an Jiaotong University.
The procedure of 3D reconstruction of real blood vessels is shown in Figure 4. After a snapshot of the DSC model was taken (Figure 4a), the binary image of the main blood vessels was transformed (Figure 4b), and the skeleton and contour lines of the blood vessels were extracted from the single 2D image (Figure 4c). After a smooth and loft operation (Figure 4d), the 3D blood vessel model was generated on the basis of the skeleton line and the contour [27] (Figure 4e). Only one image instead of multislices was used for reconstruction; thus, complex operations, such as identifying and assembling blood vessel regions between slices, were avoided, and the reconstruction steps were greatly simplified. The generated 3D geometric model was imported into GAMBIT and discretized by tetrahedral meshes.

3.2. Grid Independence Test

According to the characteristic of blood vessel geometry in the DSC model[5], the reconstructed blood vessels buried at a depth of 0.5 mm with an average diameter of 100 ± 10 μm were embedded in the skin model to verify the rupture damage mode and validate the CCPD model (Figure 5). As mentioned in Section 2.2, only the dermal layer and the blood vessel in it were considered and compared with those in the in vivo experiment. The skin depth was 1 mm, and the cross section of the skin was 1.5 mm × 1.5 mm–2 mm × 2 mm, depending on the size of each blood vessel. A 595 nm PDL with a pulse width of 6 ms and a spot diameter of 1 mm was used in the in vivo experiment and the numerical simulation. The optical and thermal properties of dermis and blood corresponding to the wavelength of 595 nm are presented in Table 1.
Six different tetrahedron grids with average interval sizes of 20, 15, 12.5, 10, 7.5, and 5 μm in the vessel region were tested for the reconstructed blood vessel (Figure 5) to check the grid independence of the present algorithm. Figure 6 shows the local Gaussian curvature at the vessel bifurcation. Figure 7 depicts the Gaussian curvature variation with different grid intervals of points A and B on the vessel wall (Figure 6). The average deviation in GC between the results with cellular sizes of 7.5 and 5 μm at the two points was less than 1.1%. The grid independence test showed that 7.5 μm was the appropriate size for temperature and pressure computation. Therefore, the grid with a cell size of 7.5 μm was adopted in the rest of our simulations.

3.3. Validation of the CCPD Model

Both the skin model and the damage model are considered as important factors in the computation. In this section, three computational models are made by combination from either of the two blood vessel geometries and the two damage models. For simplicity, the three models were referred to as models A, B, and C, which are a PD model with an ideal vessel (a straight cylinder), a PD model with a real vessel, and a CCPD model with a real vessel. With P = 2500 Pa as the pressure damage threshold [7], the location and area of the vessel rupture in the three models were computed and compared with those of the experiments.
Based on our computation, the standard deviation of the Gaussian curvature (StdGC) of a single vessel ranged from 10 to 80. In this study, the blood vessels with a diameter of 100 ± 10 μm were divided into three groups in terms of their StdGC: small (0 ≤ StdGC < 25), medium (25 ≤ StdGC < 50), and large (StdGC ≥ 50). For each group, one typical vessel was chosen to analyze pressure distribution. Figure 8 shows the experimental snapshots and the reconstructed blood vessels. For the convenience of model comparison, the position of the ideal cylindrical vessel was adjusted to coincide with the reconstructed vessel (Figure 9).
Figure 10, Figure 11 and Figure 12 show the wall pressure distribution of the three models. Model A has the largest rupture area, whereas model B has the smallest rupture area because the tortuosity of the real vessel disturbs the blood flow, which restricts the growth of large bubbles or breaks them into small bubbles. Therefore, the contact area between the region of high pressure (gas phase) and the blood vessel wall reduces, leading to a small rupture area. The rupture area of the real vessel in model C is between those of models A and B, and the rupture site agrees with the large Gaussian curvature spot on the wall surface, such as the saddle point. Figure 8c,f,i illustrate that the high-pressure position and area in model C (the present CCPD model) are more consistent with those of rupture and bleeding in the animal experiments than those in the two other models. This finding confirms the validity and feasibility of the new model in predicting vessel hemorrhage.
Figure 13 shows the relationship between Sbr obtained by the three models and Spu measured by the experiment. For a certain laser fluence, models A and B respectively show the largest and smallest rupture areas, whereas model C has moderate-sized rupture area.
As shown in Figure 13, the CCPD model (model C) predicts a more significant linear relationship between Sbr and Spu, which is in accordance with the histopathological analysis of patients with PWS [36,37]. After comprehensive consideration with our analysis in Figure 11, Figure 12 and Figure 13, the CCPD model is the best predictive model. The linear relationship between Sbr and Spu by model C fitted in Figure 13 is expressed as:
S p u = 1.01 S b r + 0.0279 .

4. Results and Discussion

In this section, the effects of different vessel curvatures, burial depths, and epidermal melanin contents (i.e., human skin color) on the parameters of the PDL treatment of PWSs were systematically studied to determine the reasonable area ratio of the vascular purpura obtained from clinical observations. The depths of blood vessels in our simulation were set as 0.25, 0.375, and 0.5 mm corresponding to the superficially, mid-, and deeply buried blood vessels, respectively. Vessel models were selected as those with small, medium, and large StdGC, as shown in Figure 8. Laser pulse duration and laser spot diameter were fixed at tp = 6 ms and D = 1 mm, respectively.
The melanin volume fractions (fm) in the epidermis were set as 1%, 3%, 5%, 7%, and 9% corresponding to people with white, light yellow, dark yellow, brown, and black skin, respectively [38]. Assuming that melanin is evenly distributed in the epidermis, we present the optical and thermal parameters of the epidermis containing melanin in Table 2. The effect of melanin content on the absorption coefficient of the epidermis was introduced in the form of volume fraction. The same optical and thermal properties of the vessel and the dermis in Table 1 were adopted.

4.1. Energy Threshold of Vessel Hemorrhage

The area of vessel rupture is closely related to purpura caused by vessel hemorrhage. For the convenience of statistical analysis of the degree of purpura in clinical patients during laser therapy, the area ratio of purpura (bleeding) Srpu is introduced as:
S r p u = S p u N S s p o t ,
where Spu is the total area of the purpura after laser treatment, N is the total number of laser pulses, and Sspot is the area of laser spot. For the DSC model in which the epidermis was removed, Srpu shows the area ratio of bleeding.
Fifty patients were treated with a 595 nm PDL in Laser Treatment Center, Department of Dermatology, Second Affiliated Hospital of Xi’an Jiaotong University to establish the relationship between Srpu predicted through numerical simulation and the treatment criterion based on clinical Srpu. The study included 34 female and 16 male PWS patients, with phototypes II-IV on the Fiztpatrick scale and an average age of 28. Spu was measured immediately after each patient was treated (Figure 14a). Sspot is the circular area with a spot diameter of 1 mm, and N is related to the morphological characteristics of vascular lesions.
The reasonable range of Srpu after treatment could be obtained in terms of the degree of disease recovery at the return visit after 1 month of treatment (Figure 14b). If the laser fluence was too high, the bleeding area was also too large, and side effects such as erythema, scar, or pigmentation appeared after treatment. If the re-examination showed that the scar remained visible, a patient was included in the group with side effects. If no side effects were found during the return visit, patients were divided into four groups based on the clearance rate. Therefore, a total of five groups were set as follows: (a) no effect (clearance rate < 30%); (b) weak effect (30% ≤ clearance rate < 60%); (c) marked effect (60% ≤ clearance rate < 90%); (d) rehabilitation (clearance rate ≥ 90%); and (e) side effects (erythema, scar, or pigmentation).
Only a suitable laser fluence could achieve an acceptable therapeutic effect. The medical community in dermatology has reached a consensus on the effectiveness in treating PWS nevus, that is, the effective cure rate is equal to the sum of the rate of marked effect and rehabilitation. Figure 15 shows the mean and confidence interval of Srpu of the five groups. Mean Srpu increases almost uniformly as the laser fluence increases, and the curative effect changes gradually from no effect to side effects. Ideal Srpu should be higher than Srpu of marked effects (c) and lower than Srpu of side effects (d), that is, higher than 13.8% and lower than 25.1%, as shown in Figure 15. The volume fraction of lesions in patients with moderate PWS nevus is approximately 20%, which is twice the volume fraction of the DSC vessels; hence, the appropriate range of Srpu in computation based on the DSC model should be between 6.9% and 12.6%. Laser fluence leading to 6.9% ≤ Srpu ≤ 12.6% in the DSC model experiment is defined as the effective laser fluence Eeff, i.e., the fluence threshold of appropriate vessel hemorrhage.

4.2. Effect of Laser Transport in Skin with PWS Vessel by TMC

The effects of 595 nm PDL transport in the skin with the curved PWS vessel by TMC were studied on the two layered skin model with a domain of 1.5 mm × 1.5 mm × 1.0 mm. The thicknesses of the epidermis and the dermis were set respectively as He = 60 μm and Hd = 940 μm and a reconstructed blood vessel buried in the skin at a depth of 0.5 mm based on statistical average values [15].
Figure 16 shows laser energy deposition in skin with blood vessel of model C (shown in Section 3.3) and melanin content of 3% in epidermis (fm = 3%). Photon deposition distribution in the cross section at y = 0.75 mm is shown in Figure 16a. Due to high absorbance of the blood vessel and epidermis, energy deposition concentrates in the two domains. The axis of the real blood vessel is circuitous and arbitrarily positioned, therefore its cross-section shows an ellipse not directly below the laser beam. It is apparent that the top right part of the vessel, which faces the photon initial direction, collects more energy than the bottom left part.
Figure 16b presents the corresponding lateral energy deposition profile. The line graph includes four parts: epidermis, dermis, blood vessel and dermis (Figure 16b). Laser energy was mostly absorbed by the blood vessel, then the epidermis and least the dermis. In each of the four parts, the energy deposition decreases with the depth except for the bottom of the blood vessel. The slight increase at the bottom of the vessel is caused by incident and deposition of the scattered light from the deep layer of the dermis. For people with the light-yellow skin (fm = 3%), the energy absorbance of the epidermis is over 300 J/cm2. In this condition, whether spray cooling is necessary depends on the results of heat transfer calculation in the next Section 4.3.

4.3. Effect of Melanin Content in the Epidermis on Laser Fluence Selection

The effects of melanin content in the epidermis were investigated from 1% to 9% on a real vessel under 595 nm PDL irradiation. Figure 17 shows Srpu of a vessel buried in the median layer (0.375 mm) at fm of 1% or 5%. The laser fluence of the vessel rupture increases with the melanin content. For example, in the blood vessel of stdGC = 41.3, the bleeding threshold increases from 4 J/cm2 to 6 J/cm2 when fm increases from 1% to 5% because the light absorption capacity of the epidermis increases with melanin content (Table 2). Subsequently, the relevant absorption ratio of the blood vessel decreases. Therefore, a high laser fluence is required to rupture the targeted blood vessel. With an increase in laser fluence, Srpu of the vessel with different values of stdGC continues to increase and gradually achieves a significant effect (6.9% < Srup < 12.6%). However, when the laser fluence continues to increase, Srup rapidly exceeds the upper limit (12.6%), and the skin experiences a side effect.
For fm = 1%, three vessels with different values of stdGC simultaneously reach the rupture threshold at E = 4 J/cm2. Among them, the vessels with a large stdGC have a large Gaussian curvature at the saddle point of the vessel surface. As a consequence, rupture and bleeding occur easily. When the laser fluence exceeds 5 J/cm2, the rupture area of the vessel with a small stdGC increases faster than that with a large stdGC because the lumen of the vessel with a small stdGC is relatively straight. As a result, the disturbance to the blood flow and the generated bubbles is weak. The bubbles grow and come in contact with the vessel wall without restriction, leading to a large high-pressure area on the vessel wall and rupture. At fm = 1%, when the laser fluence reaches 9 J/cm2, the gap in the rupture area of the three vessels is sharply narrowed because high energy greatly reduces the curvature effect. Similar results are obtained at fm = 5%.
In laser dermatologic surgery, CSC with R134a is an efficient technique to avoid epidermal damage among white people because of light energy absorption by melanin. However, even with CSC protection, people with dark pigmentation may suffer from epidermal damage. Previous studies showed that thermal damage occurs when the central temperature of the epidermis Tec is higher than 80 °C at the end of laser irradiation [40]. Therefore, the PDL fluence threshold of epidermal injury Ee,R134a is defined as the laser fluence leading to Tec = 80 °C with R134a cooling. Ee,R134a increases with an increase in the cooling effect of cryogen. Eeff should be higher than Ee,R134a to ensure the safety of the epidermis in the treatment of vascular lesions.
Figure 18 shows Eeff at different melanin contents and a vessel depth of 0.375 mm. Eeff varies by 1–2 J/cm2 with the degree of curvature regardless of fm. With the present CSC technology, for white (fm = 1%) and light-yellow (fm = 3%) people, the epidermis remains undamaged when Eeff is used. For dark-yellow people (fm = 5%), Eeff of a large-stdGC vessel exceeds Ee,R134a. Therefore, Ee,R134a should be improved by strengthening the cooling capacity of cryogen. For brown (fm = 7%) and black (fm = 9%) people, Eeff is higher than Ee,R134a. The near-infrared laser instead of PDL should be used to treat darkly pigmented patients because of its weak absorption of melanin.

4.4. Effect of Vessel Depth on the Selection of Laser Fluence

Deep PWS vessels are usually resistant to laser surgery. In this study, the effect of vessel depth on the PDL treatment of PWS nevus was investigated. Figure 19 shows Srpu of a vessel buried in the superficial layer (0.25 mm) or the deep layer (0.5 mm) when fm is 3%. A shallow buried blood vessel leads to rupture and bleeding at low fluence. When H decreases from 0.5 mm to 0.25 mm, the bleeding threshold decreases significantly from 8 J/cm2 to 2 J/cm2 because the deeper the buried vessel is, the thicker the dermis will be for light propagation to reach the vessel. In this regard, more photon energy is absorbed and scattered by the dermis, so the energy required for vessel rupture consequently increases.
At the same burial depth, the rupture thresholds and Eeff of the vessels with different stdGC vary. For the vessel depth of H = 0.25 mm and H = 0.5 mm, the rupture fluence of the vessels with a large stdGC (64.5) is lower than that with a medium or small stdGC, whereas Eeff is 1–2 J/cm2 higher. As the laser fluence increases, the vessels with a large stdGC initially rupture, subsequently experience purpura, and finally reach Eeff. The occurrence of purpura does not mean that the vessel has been properly treated. The optimal treatment parameters can only be determined by the appropriate Srpu. Simply taking the occurrence of purpura as the TEP leads to an inadequate treatment and a high recurrence rate of vessels with a large stdGC. This observation reveals one of the reasons why clinically inexperienced physicians fail to treat many patients with PWS. Therefore, in the treatment of vessels with a large stdGC, the therapeutic effect should be determined on the basis of appropriate Srpu (higher than Srpu of marked effect and lower than Srpu of side effects), and the required laser fluence should be higher than that of vessels with a small stdGC.
Figures 21 and 22 depict Eeff of laser at different melanin contents and vessel depths of 0.25 and 0.5 mm. Eeff differs by 1–2 J/cm2 with different values of stdGC. For white people (fm = 1%), superficial vessels (H = 0.25 mm) can be treated reasonably without CSC, whereas deep vessels (H = 0.5 mm) need CSC.
As shown in Figure 20, when the melanin content is equal to 7%, Ee,R134a in the superficial vessel (H = 0.25 mm) is lower than Eeff. Therefore, the cooling effect of the epidermis should be improved. For median (H = 0.375 mm) and deep (H = 0.5 mm) vessels, insufficient epidermal cooling occurs when the melanin contents are 5% and 3%, respectively. Therefore, for darkly pigmented skin, R134a cannot meet the requirement for treating median and deep blood vessels. Zhou et al. [30] found that an R134a spray can be substituted with low-boiling-point cryogen R404A (with a standard boiling point of −46.1 °C) to improve the surface cooling effect. Therefore, R404A can be used to cool the epidermis in the treatment of PWS lesions in darkly pigmented people and protect the epidermis from damage while removing vascular lesions.
The depth of blood vessels influences the selection range of Eeff. For a large-stdGC vessel (64.5), when the vessel depth is 0.25 or 0.375 mm, Eeff is 2 J/cm2 (shown as the height of the red rectangle in Figure 19 and Figure 21), whereas it reduces mostly to 1 J/cm2 when the vessel depth is 0.5 mm (shown as the height of the red rectangle in Figure 21). The difference in Ee,R134a is important because its absolute value is not high, especially when vessels are buried at a depth of 0.25 mm (5–9 J/cm2).
The simulated data in Figure 18, Figure 20 and Figure 21 are presented in Figure 22.
The relationship of effective fluence with melanin, vessel depth, and curvature can be obtained by least square fitting as follows:
E e f f = 15.536 f m 0.197 H 1.048 s t d G C 0.0861 ,
where vascular depth H is measured in millimeters. The error between the simulation results of Eeff and the predicted values of the fitting formula is less than ±16%, as shown in Figure 22. This correlation can predict the laser fluence threshold of vessel hemorrhage in the treatment of vascular lesions by a 595 nm PDL. The application conditions of Equation (16) are 1% ≤ fm ≤ 9%, 0.25 mm ≤ H ≤ 0.5 mm, and 20.5 ≤ stdGC ≤ 64.5.
In the simulation of light transport in turbid medium like skin tissues, the establishment of scattering model is important to our TMC method. For the dynamic scattering medium, special consideration should be done to build an appropriate scattering model. In the animal experiments by which our CCPD model was validated, motional disturbances by the breathing of the anesthetized rats was mostly eliminated by the fixation of translator stage and the optical chamber. But the effect of blood flow to light scattering was neglected in this work. Recently, Yang et al. [41] established an angular spectrum model to trace the field propagation during the entire optical phase conjugation process in the presence of dynamic scattering media, which can be used in the establishment of scattering model in our subsequent research. In addition, optical scattering prevents light from being focused to the target vessel through thick dermal layer at depths greater than ~1 mm [42]. Although the vessel is buried at a depth of ~0.5 mm in our simulation, the scattering effect of the skin needs to be carefully considered. In our future work, the digital optical phase conjugation based wavefront shaping technique [42] will be introduced to improve the scattering model for simulation of light propagation in skin with deep buried blood vessel.

5. Conclusions

In this study, the curvature corrected pressure damage (CCPD) model was proposed to examine the relationship between bleeding area (purpura area) and laser fluence in the treatment of real PWS vessels and accurately evaluate the location and area of a vessel rupture. The main conclusions were as follows:
(1)
Vascular lesions with a high stdGC are prone to hemorrhage and purpura formation but have difficulty in achieving Srpu to damage vascular lesions. This observation reveals one possible reason for the poor treatment of PWS in clinic, that is, purpura formation is simply regarded as the TEP, which leads to insufficient treatment and a high recurrence rate of vascular lesions with a large stdGC. An appropriate Srpu should be adopted as the TEP that is higher than Srpu of marked effects and lower than Srpu of side effects. With the ability to provide appropriate Srpu, Eeff is recommended to precisely treat vascular lesions.
(2)
Our simulation results reveal the quantitative correlations of Eeff of PDL with epidermal melanin content, vessel curvature, and depth for individualized and precise guidance on clinical treatments. To PWS patients with different melanin levels, different treatment strategies are suggested. PDL should be used carefully in treating light-yellow people (fm = 3%). Cryogen R404A with a good cooling capacity is suggested to protect the epidermis with fm = 5%. For brown (fm = 7%) and black (fm = 9%) people, the implementation of Eeff of blood vessels in PDL severely harms the epidermis; thus, near-infrared lasers are recommended because they minimally absorb melanin.

Author Contributions

Conceptualization, H.J., B.C. and D.L.; Methodology, H.J.; Software, H.J.; Validation, H.J., B.C.; Formal analysis, D.L.; Investigation, H.J.; Data curation, H.J.; Writing—original draft preparation, H.J.; writing—review and editing, B.C.; supervision, B.C.; Project administration, B.C.; Funding acquisition, B.C.

Funding

This work was supported by the National Natural Science Foundation of China (51727811), Zhejiang Provincial Natural Science Foundation of China (LZ15E090002) and Science Foundation of Zhejiang Sci-Tech University (ZSTU) (19022101-Y, 11130231711926).

Acknowledgments

Thanks for the technical support of the Second Affiliated Hospital of Xi’an Jiaotong University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aguilar, G.; Choi, B.; Broekgaarden, M.; Yang, O.; Yang, B.; Ghasri, P.; Chen, J.K.; Bezemer, R.; Nelson, J.S.; van Drooge, A.M.; et al. An overview of three promising mechanical, optical, and biochemical engineering approaches to improve selective photothermolysis of refractory port wine stains. Ann. Biomed. Eng. 2012, 40, 486–506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Passeron, T.; Maza, A.; Fontas, E.; Toubel, G.; Vabres, P.; Livideanu, C.; Mazer, J.-M.; Rossi, B.; Boukari, F.; Harmelin, Y.; et al. Treatment of port wine stains with pulsed dye laser and topical timolol: A multicenter randomized controlled trial. Br. J. Dermatol. 2014, 170, 1350–1353. [Google Scholar] [CrossRef] [PubMed]
  3. Smith, T.K.; Choi, B.; Ramirez-San-Juan, J.C.; Nelson, J.S.; Osann, K.; Kelly, K.M. Microvascular blood flow dynamics associated with photodynamic therapy, pulsed dye laser irradiation and combined regimens. Lasers Surg. Med. 2006, 38, 532–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Anderson, R.R.; Parrish, J.A. Selective photothermolysis: Precise microsurgery by selective absorption of pulsed radiation. Science 1983, 220, 524–527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Li, D.; Chen, B.; Wu, W.J.; Wang, G.X.; He, Y.L.; Ying, Z.X. Experimental study on the vascular thermal response to visible laser pulses. Lasers Med. Sci. 2015, 30, 135–145. [Google Scholar] [CrossRef] [PubMed]
  6. Kelly, K.M.; Choi, B.S.; Motosue, A.; Jung, B.; Khan, M.H.; Jc, R.S.J.; Nelson, J.S. Description and analysis of treatments for port-wine stain birthmarks. Arch. Facial Plast. Surg. 2005, 7, 287–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Jia, H.; Chen, B.; Li, D. Criteria of pressure and thermal damage during laser irradiation of port wine stains: Which is dominant to vascular lesions? Int. J. Heat Mass Transf. 2019, 132, 848–860. [Google Scholar] [CrossRef]
  8. Li, D.; Farshidi, D.; Wang, G.X.; He, Y.L.; Kelly, K.M.; Wu, W.J.; Chen, B.; Ying, Z.X. A comparison of microvascular responses to visible and near-infrared lasers. Lasers Surg. Med. 2014, 46, 479–487. [Google Scholar] [CrossRef] [Green Version]
  9. Goldman, M.P. Cutaneous and Cosmetic Laser Surgery; Elsevier Pte Ltd.: Singapore, 2006. [Google Scholar]
  10. Goldberg, D.J. Laser and Lights; Saunders/Elsevier: Amsterdam, The Netherlands, 2009. [Google Scholar]
  11. Chen, B.; Thomsen, S.L.; Thomas, R.J.; Oliver, J.; Welch, A.J. Histological and modeling study of skin thermal injury to 2.0 μm laser irradiation. Lasers Surg. Med. 2008, 40, 358–370. [Google Scholar] [CrossRef]
  12. Jia, W.; Choi, B.; Franco, W.; Lotfi, J.; Majaron, B.; Aguilar, G.; Nelson, J.S. Treatment of cutaneous vascular lesions using multiple-intermittent cryogen spurts and two-wavelength laser pulses: Numerical and animal studies. Lasers Surg. Med. 2007, 39, 494–503. [Google Scholar] [CrossRef] [Green Version]
  13. Welch, A.J.; Gemert, M.J.C.V. Optical-Thermal Response of Laser-Irradiated Tissue; Springer: Dordrecht, The Netherlands, 2011; pp. 579–589. [Google Scholar]
  14. Zhang, J.Z.; Zhang, X.X.; Audette, M. A photothermal model of selective photothermolysis with dynamically changing vaporization temperature. Lasers Med. Sci. 2011, 26, 633. [Google Scholar] [CrossRef] [PubMed]
  15. Li, D.; Wang, G.X.; He, Y.L.; Kelly, K.M.; Wu, W.J.; Wang, Y.X.; Ying, Z.X. A two-temperature model for selective photothermolysis laser treatment of port wine stains. Appl. Therm. Eng. 2013, 59, 41–51. [Google Scholar] [CrossRef] [PubMed]
  16. Suthamjariya, K.; Farinelli, W.A.; Koh, W.; Anderson, R.R. Mechanisms of microvascular response to laser pulses. J. Investig. Dermatol. 2004, 122, 518–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Shen, Y.; Zhang, A.; Xu, L.X. A study on mechanical damage of tumor microvasculature induced by alternate cooling and heating. J. Therm. Sci. Eng. Appl. 2009, 1, 031002. [Google Scholar] [CrossRef]
  18. Han, H.-C.; Chesnutt, J.K.W.; Garcia, J.R.; Liu, Q.; Wen, Q. Artery buckling: New phenotypes, models, and applications. Ann. Biomed. Eng. 2013, 41, 1399–1410. [Google Scholar] [CrossRef]
  19. Smithies, D.J.; Butler, P.H. Modelling the distribution of laser light in port-wine stains with the Monte Carlo method. Phys. Med. Biol. 1995, 40, 701. [Google Scholar] [CrossRef]
  20. Pfefer, T.J.; Barton, J.K.; Smithies, D.J. Modeling Laser Treatment of Port Wine Stains With a Computer-Reconstructed Biopsy. Lasers Surg. Med. 1999, 24, 151–166. [Google Scholar] [CrossRef]
  21. Daly, S.M.; Leahy, M.J. ‘Go with the flow’: A review of methods and advancements in blood flow imaging. J. Biophotonics 2013, 6, 217–255. [Google Scholar] [CrossRef]
  22. Barton, J.K.; Frangineas, G.; Pummer, H.; Black, J.F. Cooperative phenomena in two-pulse, two-color laser photocoagulation of cutaneous blood vessels. Photochem. Photobiol. 2001, 73, 642–650. [Google Scholar] [CrossRef]
  23. Babilas, P.; Shafirstein, G.; Bäumler, W.; Baier, J.; Landthaler, M.; Szeimies, R.M.; Abels, C. Selective Photothermolysis of Blood Vessels Following Flashlamp-Pumped Pulsed Dye Laser Irradiation: In Vivo Results and Mathematical Modelling Are in Agreement. J. Investig. Dermatol. 2005, 20, 343–352. [Google Scholar] [CrossRef]
  24. Li, D.; Li, R.; Jia, H.; Chen, B.; Wu, W.; Ying, Z. Experimental and numerical investigation on the transient vascular thermal response to multi-pulse Nd:YAG laser. Lasers Surg. Med. 2017, 49, 852–865. [Google Scholar] [CrossRef] [PubMed]
  25. Moy, A.J.; White, S.M.; Indrawan, E.S.; Lotfi, J.; Nudelman, M.J.; Costantini, S.J.; Agarwal, N.; Jia, W.; Kelly, K.M.; Sorg, B.S.; et al. Wide-field functional imaging of blood flow and hemoglobin oxygen saturation in the rodent dorsal window chamber. Microvasc. Res. 2011, 82, 199–209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Palmer, G.M.; Fontanella, A.N.; Shan, S.; Hanna, G.; Zhang, G.; Fraser, C.L.; Dewhirst, M.W. In vivo optical molecular imaging and analysis in mice using dorsal window chamber models applied to hypoxia, vasculature and fluorescent reporters. Nat. Protoc. 2011, 6, 1355–1366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Mu, L.Z.; Shao, H.W.; He, Y.; Oda, T.; Jia, X.M. Construction of anatomically accurate finite element models of the human hand and a rat kidney. J. Mech. Med. Biol. 2012, 11, 1141–1164. [Google Scholar] [CrossRef]
  28. Jia, H.; Chen, B.; Li, D.; Zhang, Y. Boundary discretization in the numerical simulation of light propagation in skin tissue: Problem and strategy. J. Biomed. Opt. 2015, 20, 25007. [Google Scholar] [CrossRef] [PubMed]
  29. Pfefer, T.J.; Barton, J.K.; Chan, E.K.; Ducros, M.G. A three-dimensional modular adaptable grid numerical model for light propagation during laser irradiation of skin tissue. IEEE J. Sel. Top. Quantum Electron. 1996, 2, 934–942. [Google Scholar] [CrossRef] [Green Version]
  30. Li, D.; Chen, B.; Wu, W.J.; Wang, G.X.; He, Y.L. Multi-scale modeling of tissue freezing during cryogen spray cooling with R134a, R407c and R404a. Appl. Therm. Eng. 2014, 73, 1489–1500. [Google Scholar] [CrossRef]
  31. Xiang, H.; Chen, B.; Wu, W.J.; Zhang, Y.; Jia, H. An integral MPS model of blood coagulation by laser irradiation: Application to the optimization of multi-pulse Nd:YAG laser treatment of port-wine stains. Int. J. Heat Mass Transf. 2017, 114, 1220–1233. [Google Scholar] [CrossRef]
  32. Meyer, M.; Desbrun, M.; Schröder, P.; Barr, A.H. Discrete Differential-Geometry Operators for Triangulated 2-Manifolds. Vis. Math. 2003, 35–57. [Google Scholar]
  33. Bowyer, A. Computing Dirichlet tessellations. Comput. J. 1981, 24, 162–166. [Google Scholar] [CrossRef] [Green Version]
  34. Watson, D.F. Computing the n-dimensional Delaunay tessellation with application to Voronoi polytopes. Comput. J. 1981, 24, 167–172. [Google Scholar] [CrossRef] [Green Version]
  35. Ruiz, d.G.S.; Cazón, A.; Antón, R.; Finol, E.A. A Methodology for Verifying Abdominal Aortic Aneurysm Wall Stress. J. Biomech. Eng. 2017, 139, 011006. [Google Scholar] [CrossRef] [PubMed]
  36. Koster, P.H.; van der Horst, C.M.; van der Wal, A.C. Histologic evaluation of skin damage after overlapping and nonoverlapping flashlamp pumped pulsed dye laser pulses: A study on normal human skin as a model for port wine stains. Lasers Surg. Med. 2001, 28, 176–181. [Google Scholar] [CrossRef] [PubMed]
  37. Tanghetti, E.; Sherr, E.A.; Sierra, R.; Mirkov, M. The effects of pulse dye laser double-pass treatment intervals on depth of vessel coagulation. Lasers Surg. Med. 2006, 38, 16–21. [Google Scholar] [CrossRef]
  38. Milanic, M.; Majaron, B. Three-dimensional Monte Carlo model of pulsed-laser treatment of cutaneous vascular lesions. J. Biomed. Opt. 2011, 16, 128002. [Google Scholar] [CrossRef]
  39. Li, D.; Chen, B.; Wang, G.X.; Wu, W.; He, Y.L. A New Theoretical Model of Selectively Photothermolysis to Aid Laser Treatment of Poor Responding Port Wine Stain Blood Vessels. Begel House Inc 2014, 6993–7006. [Google Scholar]
  40. Li, D.; Chen, B.; Wu, W.J.; Ying, Z.X. Influence of laser wavelength on the thermal responses of port wine stain lesions in light, moderate and heavy pigmented skin. Appl. Therm. Eng. 2017, 117, 193–203. [Google Scholar] [CrossRef]
  41. Yang, J.; Shen, Y.; Liu, Y.; Hemphill, A.S.; Wang, L.V. Focusing light through scattering media by polarization modulation based generalized digital optical phase conjugation. Appl. Phys. Lett. 2017, 111, 201108. [Google Scholar] [CrossRef]
  42. Yang, J.; Li, J.; He, S.; Wang, L.V. Angular-spectrum modeling of focusing light inside scattering media by optical phase conjugation. Optica 2019, 6, 250–256. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Computational domain of the two-layer skin model with a reconstructed blood vessel.
Figure 1. Computational domain of the two-layer skin model with a reconstructed blood vessel.
Applsci 09 05478 g001
Figure 2. Finite-volume region on a one-ring triangulated surface of xi by using Voronoi cells. For the point inside each adjacent triangle, the surface area with the circumcenters is recognized as the local Voronoi cell (region surrounded by the blue lines).
Figure 2. Finite-volume region on a one-ring triangulated surface of xi by using Voronoi cells. For the point inside each adjacent triangle, the surface area with the circumcenters is recognized as the local Voronoi cell (region surrounded by the blue lines).
Applsci 09 05478 g002
Figure 3. Dorsal skin chamber model.
Figure 3. Dorsal skin chamber model.
Applsci 09 05478 g003
Figure 4. Procedure of the 3D reconstruction of real blood vessels. (a) Taking a photo of the DSC model; (b) extracting the main blood vessels and transforming them to a binary image; (c) extracting the contour and skeleton lines; (d) smoothing the line and outputting it in a vector form; (e) conducting 3D angiogenesis.
Figure 4. Procedure of the 3D reconstruction of real blood vessels. (a) Taking a photo of the DSC model; (b) extracting the main blood vessels and transforming them to a binary image; (c) extracting the contour and skeleton lines; (d) smoothing the line and outputting it in a vector form; (e) conducting 3D angiogenesis.
Applsci 09 05478 g004
Figure 5. Geometric skin model with a 3D reconstructed blood vessel. (a) Vessel image of the dorsal skin chamber model (the main vessel of reconstruction is indicated by the yellow arrow). (b) Skin model with a 3D reconstructed vessel.
Figure 5. Geometric skin model with a 3D reconstructed blood vessel. (a) Vessel image of the dorsal skin chamber model (the main vessel of reconstruction is indicated by the yellow arrow). (b) Skin model with a 3D reconstructed vessel.
Applsci 09 05478 g005
Figure 6. Local Gaussian curvature obtained from the mesh with different intervals.
Figure 6. Local Gaussian curvature obtained from the mesh with different intervals.
Applsci 09 05478 g006
Figure 7. Influence of the mesh on the calculation of the Gaussian curvature on the bifurcated vessels.
Figure 7. Influence of the mesh on the calculation of the Gaussian curvature on the bifurcated vessels.
Applsci 09 05478 g007
Figure 8. Responses of blood vessels irradiated with a 595 nm laser with a pulse width of 6 ms (yellow arrow refers to the irradiated reconstructed vessel) (a) Before irradiation (stdGC = 20.5); (b) reconstructed model (stdGC = 20.5); (c) after irradiation (stdGC = 20.5); (d) before irradiation (stdGC = 41.3); (e) reconstructed model (stdGC = 41.3); (f) after irradiation (stdGC = 41.3); (g) before irradiation (stdGC = 64.5); (h) reconstructed model (stdGC = 64.5); and (i) after irradiation (stdGC = 64.5).
Figure 8. Responses of blood vessels irradiated with a 595 nm laser with a pulse width of 6 ms (yellow arrow refers to the irradiated reconstructed vessel) (a) Before irradiation (stdGC = 20.5); (b) reconstructed model (stdGC = 20.5); (c) after irradiation (stdGC = 20.5); (d) before irradiation (stdGC = 41.3); (e) reconstructed model (stdGC = 41.3); (f) after irradiation (stdGC = 41.3); (g) before irradiation (stdGC = 64.5); (h) reconstructed model (stdGC = 64.5); and (i) after irradiation (stdGC = 64.5).
Applsci 09 05478 g008
Figure 9. Reconstructed vessel models with different curvatures and corresponding cylindrical vessel models (the dark red vessel is the 3D reconstructed vessel, and the light red vessel is the corresponding cylindrical vessel). The position of the ideal vessel was adjusted to coincide with the reconstructed vessel.
Figure 9. Reconstructed vessel models with different curvatures and corresponding cylindrical vessel models (the dark red vessel is the 3D reconstructed vessel, and the light red vessel is the corresponding cylindrical vessel). The position of the ideal vessel was adjusted to coincide with the reconstructed vessel.
Applsci 09 05478 g009
Figure 10. Pressure distribution on the wall of small curvature (stdGC = 20.5 mm−2) vessels after laser irradiation (E = 11 J/cm2, tp = 6 ms).
Figure 10. Pressure distribution on the wall of small curvature (stdGC = 20.5 mm−2) vessels after laser irradiation (E = 11 J/cm2, tp = 6 ms).
Applsci 09 05478 g010
Figure 11. Pressure distribution on the wall of the median curvature (stdGC = 41.3 mm−2) vessels after laser irradiation (E = 12 J/cm2, tp = 6 ms).
Figure 11. Pressure distribution on the wall of the median curvature (stdGC = 41.3 mm−2) vessels after laser irradiation (E = 12 J/cm2, tp = 6 ms).
Applsci 09 05478 g011
Figure 12. Pressure distribution on the wall of the median curvature (stdGC = 64.5 mm−2) vessels after laser irradiation (E = 13 J/cm2, tp = 6 ms).
Figure 12. Pressure distribution on the wall of the median curvature (stdGC = 64.5 mm−2) vessels after laser irradiation (E = 13 J/cm2, tp = 6 ms).
Applsci 09 05478 g012
Figure 13. Relationship between Sbr predicted by three damage models and Spu measured by the experiments under different incident energy densities.
Figure 13. Relationship between Sbr predicted by three damage models and Spu measured by the experiments under different incident energy densities.
Applsci 09 05478 g013
Figure 14. Purpura and curative effect of patients with port wine stains and treated with a 595 nm pulsed dye laser (a) Purpura occurs immediately after irradiation. (b) Return visit after 1 month of treatment (marked effect).
Figure 14. Purpura and curative effect of patients with port wine stains and treated with a 595 nm pulsed dye laser (a) Purpura occurs immediately after irradiation. (b) Return visit after 1 month of treatment (marked effect).
Applsci 09 05478 g014
Figure 15. Relationship between the therapeutic effect of PWS and Srpu (Data source: Laser Treatment Center, Department of dermatology, Second Affiliated Hospital of Xi’an Jiaotong University).
Figure 15. Relationship between the therapeutic effect of PWS and Srpu (Data source: Laser Treatment Center, Department of dermatology, Second Affiliated Hospital of Xi’an Jiaotong University).
Applsci 09 05478 g015
Figure 16. Energy deposition in skin with blood vessel model C (fm = 3%, 595 nm PDL, 12 J/cm2). (a) Cross-sectional maps of laser energy deposited at y = 0.75 mm. (b) Line graph of energy deposition versus depth at the vessel center (shown as white dashed line in Figure 16a) for simulations of skin with single blood vessel.
Figure 16. Energy deposition in skin with blood vessel model C (fm = 3%, 595 nm PDL, 12 J/cm2). (a) Cross-sectional maps of laser energy deposited at y = 0.75 mm. (b) Line graph of energy deposition versus depth at the vessel center (shown as white dashed line in Figure 16a) for simulations of skin with single blood vessel.
Applsci 09 05478 g016
Figure 17. Purpura rate of blood vessels at a depth of 0.375 mm (the laser fluence of vessel rupture is marked by the black circles): (a) fm = 1% (b) fm = 5%.
Figure 17. Purpura rate of blood vessels at a depth of 0.375 mm (the laser fluence of vessel rupture is marked by the black circles): (a) fm = 1% (b) fm = 5%.
Applsci 09 05478 g017
Figure 18. Eeff of laser at different melanin contents with the vessel buried at a depth of 0.375 mm (marked by solid rectangles).
Figure 18. Eeff of laser at different melanin contents with the vessel buried at a depth of 0.375 mm (marked by solid rectangles).
Applsci 09 05478 g018
Figure 19. Purpura rate of blood vessels with melanin content fm = 3%. (a) Vessel burial depth H = 0.25 mm (b) Vessel burial depth H = 0.5 mm
Figure 19. Purpura rate of blood vessels with melanin content fm = 3%. (a) Vessel burial depth H = 0.25 mm (b) Vessel burial depth H = 0.5 mm
Applsci 09 05478 g019
Figure 20. Eeff of laser at different melanin contents in vessels buried at a depth of 0.25 mm (marked by solid rectangles).
Figure 20. Eeff of laser at different melanin contents in vessels buried at a depth of 0.25 mm (marked by solid rectangles).
Applsci 09 05478 g020
Figure 21. Eeff of laser at different melanin contents in vessels buried at a depth of 0.5 mm (marked by solid rectangles).
Figure 21. Eeff of laser at different melanin contents in vessels buried at a depth of 0.5 mm (marked by solid rectangles).
Applsci 09 05478 g021
Figure 22. Simulation results of Eeff with different parameters of vessel and epidermal melanin.
Figure 22. Simulation results of Eeff with different parameters of vessel and epidermal melanin.
Applsci 09 05478 g022
Table 1. Optical and hydrodynamic properties of dermis and blood [7].
Table 1. Optical and hydrodynamic properties of dermis and blood [7].
PropertyDermisBlood
Optical propertiesAbsorption coefficient, µa/cm−12049.3
Scattering coefficient, µs/cm−1460466
Anisotropy index, g0.80.995
Refractive index, n1.371.33
Thermal propertiesDensity, ρ/kg·m−310901060
Thermal conductivity, k/kW·m−1·K−10.410.55
Specific heat, c/J·kg−1·K−135003600
Table 2. Optical and thermal properties of the epidermis irradiated by a 595 nm pulsed dye laser [39].
Table 2. Optical and thermal properties of the epidermis irradiated by a 595 nm pulsed dye laser [39].
PropertyEpidermis (No Melanin)MelaninEpidermis with fm Melanin
Absorption coefficient, µa/cm−10.371402.63μa,m × fm + μa,e × (1 − fm)
Scattering coefficient, µs/cm−1470------470
Anisotropy index, g0.79------0.79
Refractive index, n1.37------1.37
Density, ρ/kg·m−31120------1120
Thermal conductivity, k/kW·m−1·K−10.34------0.34
Specific heat, c/J·kg−1·K−13200------3200

Share and Cite

MDPI and ACS Style

Jia, H.; Chen, B.; Li, D. Theoretical Study on Pressure Damage Based on Clinical Purpura during the Laser Irradiation of Port Wine Stains with Real Complex Vessels. Appl. Sci. 2019, 9, 5478. https://0-doi-org.brum.beds.ac.uk/10.3390/app9245478

AMA Style

Jia H, Chen B, Li D. Theoretical Study on Pressure Damage Based on Clinical Purpura during the Laser Irradiation of Port Wine Stains with Real Complex Vessels. Applied Sciences. 2019; 9(24):5478. https://0-doi-org.brum.beds.ac.uk/10.3390/app9245478

Chicago/Turabian Style

Jia, Hao, Bin Chen, and Dong Li. 2019. "Theoretical Study on Pressure Damage Based on Clinical Purpura during the Laser Irradiation of Port Wine Stains with Real Complex Vessels" Applied Sciences 9, no. 24: 5478. https://0-doi-org.brum.beds.ac.uk/10.3390/app9245478

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop