# Mathematical modeling of endovenous laser treatment (ELT)

- Serge R Mordon
^{1}Email author, - Benjamin Wassmer
^{2}and - Jaouad Zemmouri
^{2}

**5**:26

**DOI: **10.1186/1475-925X-5-26

© Mordon et al; licensee BioMed Central Ltd. 2006

**Received: **14 February 2006

**Accepted: **25 April 2006

**Published: **25 April 2006

## Abstract

### Background and objectives

Endovenous laser treatment (ELT) has been recently proposed as an alternative in the treatment of reflux of the Great Saphenous Vein (GSV) and Small Saphenous Vein (SSV). Successful ELT depends on the selection of optimal parameters required to achieve an optimal vein damage while avoiding side effects. Mathematical modeling of ELT could provide a better understanding of the ELT process and could determine the optimal dosage as a function of vein diameter.

### Study design/materials and methods

The model is based on calculations describing the light distribution using the diffusion approximation of the transport theory, the temperature rise using the bioheat equation and the laser-induced injury using the Arrhenius damage model. The geometry to simulate ELT was based on a 2D model consisting of a cylindrically symmetric blood vessel including a vessel wall and surrounded by an infinite homogenous tissue. The mathematical model was implemented using the Macsyma-Pdease2D software (Macsyma Inc., Arlington, MA, USA). Damage to the vein wall for CW and single shot energy was calculated for 3 and 5 mm vein diameters. In pulsed mode, the pullback distance (3, 5 and 7 mm) was considered. For CW mode simulation, the pullback speed (1, 2, 3 mm/s) was the variable. The total dose was expressed as joules per centimeter in order to perform comparison to results already reported in clinical studies.

### Results

In pulsed mode, for a 3 mm vein diameter, irrespective of the pullback distance (2, 5 or 7 mm), a minimum fluence of 15 J/cm is required to obtain a permanent damage of the intima. For a 5 mm vein diameter, 50 J/cm (15W-2s) is required. In continuous mode, for a 3 mm and 5 mm vein diameter, respectively 65 J/cm and 100 J/cm are required to obtain a permanent damage of the vessel wall. Finally, the use of different wavelengths (810 nm or 980 nm) played only a minor influence on these results.

### Discussion and conclusion

The parameters determined by mathematical modeling are in agreement with those used in clinical practice. They confirm that thermal damage of the inner vein wall (tunica intima) is required to achieve the tissue alterations necessary in order to lead the vein to permanent occlusion. However, in order to obtain a high rate of success without adverse events, the knowledge of the vein diameter after tumescent anesthesia is recommended in order to use the optimal energy. As clearly demonstrated by our calculations, both pulsed and continuous mode operations of the laser can be efficient. An interesting observation in our model is that less amount of energy is required in pulsed mode than in continuous mode. Damaging the vein sequentially along its entire length may lead to permanent occlusion. However, the pulsed mode requires a very precise positioning of the fiber after each pullback and the duration of the treatment is much longer. For these reasons, continuous irradiation seems to be preferred by most clinicians. This model should serve as a useful tool to simulate and better understand the mechanism of action of the ELT

## Introduction

Lower-extremity venous insufficiency is a common medical condition afflicting 25% of women and 15% of men in the United States and in Europe. Great saphenous vein (GSV) reflux is the most common underlying cause of significant varicose veins. Traditional treatment of GSV reflux has been surgical removal of the GSV. Although surgical ligation and stripping of the GSV has been the most durable treatment, it is associated with significant peri-operative morbidity. Less-invasive surgical treatments including high ligation of the GSV at the saphenofemoral junction (SFJ) have been attempted in the hope that gravitational reflux would be controlled while the vein is preserved for possible use as a bypass graft. Unfortunately, ligation of the GSV alone usually results in recurrent varicose veins. Even when high ligation has been combined with phlebectomy of varicose tributaries or retrograde sclerotherapy, recurrence has been the rule. Therefore, when it is determined that GSV reflux is the principal underlying problem, treatment should involve eliminating this source of reflux with ablation of any associated incompetent venous segments [1].

In an attempt to reduce morbidity and improve recovery time, several minimally invasive techniques have been developed as alternatives to surgery in the last few years. Endovenous laser treatment (ELT) is one of the most promising of these new techniques [2–4]. In 1999, Boné first reported the delivery of endoluminal laser energy [5].

Numerous studies have since demonstrated that this technique is both safe and efficacious. Several wavelengths have been proposed, respectively 810, 940, 980, 1064 and 1320 nm [6–10] with 810, 940 and 980 the most commonly used. At these wavelengths, power is usually set between 10 and 15 W. The energy is administered endovenously, either in a pulsed fashion (pulse duration: 1 to 3 s with fiber pull back in 3 to 5 mm increments every 2 seconds) or continuously with a constant pullback of the laser fiber (pullback velocity ranging from 1 to 3 mm). At these parameters, doses applied range from 20 J/cm to 140 J/cm [11, 12]. These doses induce an heating of the vein wall which is necessary to cause collagen contraction and destruction of endothelium. This stimulates vein wall thickening leading to luminal contraction, venous thrombosis and vein fibrosis [13]. Since tumescent anesthesia is always delivered, patients feel no pain during endovenous laser ablation at the suggested or commonly used laser parameters. The pain that patients feel occurs 5–8 days following the procedure and is related to the inflammation resulting from a successful endovenous ablation (i.e. wall thickening). It is not related to the presence or degree of ecchymosis nor is it the result of non-target laser damage to perivenous tissue. However, if greater doses of energy are delivered, the treatment is becoming painful. Excessive and nonspecific thermal damage has led to 67% of patients complaining of pain along the treated vein for 1 week in one study using the 940-nm intravascular laser [14].

Consequently, successful ELT must achieve permanent damage of the vessel wall, with both the endothelium (tunica intima) and also the outer layer (tunica externa) requiring sufficient heating. Conversely, to avoid side effects "collateral" damage of the perivenous tissue must be avoided. In this context, reports concerning temperatures reached inside vessel lumen during ELT are controversial. In the model studied by Zimmet and Min, peak temperatures (maximum: 50°C) of perivenous tissues generated during ELT (810 nm, 12 W, 1.5s-pulse duration) were unlikely to cause permanent damage to the perivenous tissues. Using another experimental model, Weiss showed that 100% of the laser (810 nm, 12 W, 1s-pulse duration) treated veins showed perforations.

Mathematical modeling of ELT may provide a better understanding of the ELT process and help determine optimal dosage as a function of vein diameter. Although many other models address the laser treatments of blood vessels (Port Wine Stains for example) or Interstitial laser therapy, the mathematical modeling of ELT is yet to be proposed [15–19].

The aim of this paper is to present a mathematical model using dynamic tissue changes based upon the Arrhenius damage model. Numerical simulations are compared to data previously reported in the literature. Optimal parameters emerging from these calculations can be taken in account for improving ELT clinical application.

## Materials and methods

### Mathematical modeling

An opto-thermal model of ELT consists of calculations of light distribution, temperature rise and the extent of thermal damage. The following sections describe the manner in which each stage has been implemented in our calculations.

### Geometrical description of the model

The geometry used to simulate ELT was based on a 2D model consisting of a cylindrical blood vessel (radius: R) including a vessel wall (Thickness: T) and surrounded by infinite and homogenous tissue. Calculations were performed for different vein diameters (3, 5 mm) and at different distances from the center of the vein: in the tunica intima, in the tunica externa and at 1, 2 and 3 mm from the tunica externa.

### Light distribution in tissue

The light emitted from the fiber inserted in the vein was modeled as an isotropically radiating point source. As previously proposed by Lizuka et al, spatial distribution has been considered to be dominated by scattering processes [20]. The light irradiance rate (W.mm^{-2}) of an isotropic point source emitting P_{laser} (W) within an infinite homogeneous medium can be expressed as

where: *P*
_{
Laser
}(W) : power of the light source

*μ*
_{
eff
}(mm^{-1}) : effective attenuation coefficient

*r* (mm) : radial distance from the source

*D* (mm): optical diffusion distance

*μ*
_{
eff
}is determined by the following equation

where *μ*
_{
a
}(mm^{-1}) : absorption coefficient in tissue

*μ*'_{
s
}(mm^{-1}) : reduced scattering coefficient:
= *μ*
_{
s
}(1 - *g*)

*μ*
_{
s
}(mm^{-1}) : scattering coefficient

*g* : anisotropy factor incorporating the effects of directionally dependent scattering.

*D (mm)* is determined by the following equation:

*r* is defined by the following equation

Where: x (mm): transverse dimension

z (mm): longitudinal dimension

The absorbed power density (W.mm^{-3}) is expressed as follows (Welch 1984):

*P*
_{
abs
}= *μ*
_{
a
}φ(r) (4)

The first laser pulse is always applied at coordinates (0,0) in figure 1. When using several pulses, the relative position of the fiber inside the vein is given by:

*z*' = *z* - *z*
_{
inc
}

Where: z_{inc} (mm) is the absolute position of the fiber inside the vein. This position is calculated for each pulse by taking into account the pull-back distance.

When simulating a continuous irradiation performed with a progressive pull-back of the fiber, the relative position of each irradiation is obtained by taking into account the pull-back speed:

*z*
_{
inc
}= *t* × *v*

where : v (mm.s^{-1}) is the pull back speed.

### Calculation of temperature rise

Absorption of light in tissue causes a local elevation in temperature. Tissue heat transfer due to the deposited light is described by the bioheat transfer equation as proposed by Zhang et al [21].

Where

*T* (*r*,*t*): temperature (°K)

ρ : density of tissue (g mm^{-3})

*C* : specific heat of tissue (J. g^{-1}.°K^{-1})

*C*
_{
p
}= *C*·ρ: *heat capacity* (J.mm^{-3}.°K^{-1})

*k* = thermal conductivity of tissue (W. mm^{-1}.°K^{-1})

*r* = radial distance (mm)

*t* = time (s)

listing of physical parameters used for numerical simulation

Blood | Vessel wall | Perivenous tissue | References | |
---|---|---|---|---|

| 0.20 (810 nm) 0.28 (980 nm) | 0.1 | 0.030 | [32,33] |

| 0.70 (810 nm) 0.6 (980 nm) | 2.0 | 1.0 | [32,33] |

| 0,86 (980 nm) | 0.79 | 0.30 | [32,34] |

| 3.82 | 3.78 | 3.78 | [33,34] |

| 1.05.10 | 1.05.10 | 1.05.10 | [33,34] |

| 5.6.10 | 5.6.10 | 5.6.10 | [33,34] |

E | 4.48.10 | 4.30.10 | 4.30.10 | [34,35] |

A (s | 7.6.10 | 5.6.10 | 5.6.10 | [34,35] |

### Phase transition of blood

In our analysis, we considered that the maximum temperature of blood would not exceed 100°C for all laser parameters. Water being the main constituent of blood, the latent heat of water was included in the calculations. It takes 4.18 J to raise the temperature of one gram of water by one degree C. When considering a 50 mm^{3} cylinder (diameter: 4 mm, length: 4 mm), the energy required to reach 100°C is 13 J. The latent heat of vaporization of water is 2260 kJ/kg at 100°C. The relationship between pressure and temperature is given by the phase transition diagram of water. At normal pressure: 1.10^{5} Pa and body temperature 37°C, water is in a liquid state. When the temperature increases, two phenomena can be observed: 1) pressure increases and the water stays in its liquid state, 2) pressure remains constant and water reaches its gaseous state, (steam or vapor) [22] To reach the gaseous state of this 50 mm^{3} cylinder, (steam or vapor) 119 J will be required. It is therefore unlikely that the blood will boil during ELT.

If we consider the vessel to be a closed but deformable container, water evaporation will not occur, and consequently, water will return to its initial liquid state as soon as the temperature drops back to below 100°C. In conclusion, the assumption is that blood temperature stays around 100°C during laser irradiation is both valid and confirmed by previous studies [23].

### Damage function

Thermal damage in cells and tissue can be described mathematically by a first-order thermal-chemical rate equation, in which temperature history determines damage. Damage is considered to be a unimolecular process, where native molecules transform into a denatured/coagulated state through an activated state leading to cell death. Damage is quantified using a single parameter Ω, which ranges on the entire positive real axis and is calculated from the Arrhenius law [24]. Damage Ω is dimensionless, exponentially dependent on temperature, and linearly dependant on time of exposure.

where *A* (s^{-1}) is the frequency factor,

*E*
_{
a
}(J. mole^{-1}) is the activation energy,

*R* (J. mole^{-1}.°K^{-1}) is the universal gas constant,

The activation energy *E*
_{
a
}and the frequency factor *A* are derived from thermodynamic variables. They describe the denaturation process of proteins and other cellular constituents. *A* ranges from 10^{40}s^{-1} to 10^{105}s^{-1}, and *E*
_{
a
}from 10^{5}J/mole to 10^{6}J/mole [25]. Values used for calculations are reported in Table 1.

*T* (°K) is the temperature.

Ω can be determined by the following equation

where: *C*
_{0} is the concentration of the undamaged molecules at the beginning

*C(t)* is the concentration of the undamaged molecules at time *τ*.

The equation indicates that the measure of damage describes the probability of tissue being destroyed. It is the logarithm of the ratio of the initial concentration of undamaged tissue to the concentration once damage has accumulated, for the time interval t = 0 to t = *τ*. Therefore, Ω = 1 corresponds to an irreversible damage of 100% of the affected cells.

### Numerical implementation

The mathematical model was implemented using the Macsyma-Pdease2D software (Macsyma Inc., Arlington, MA, USA). This Finite Element CAD Software needs to specify the Partial Differential Equations, variables, geometry, and boundary conditions; PDEase2D creates both numerical output tables and plots. PDEase2D generates and refines the element grid, adaptively selects time steps, and iterates until it attains convergence in nonlinear problems. And because PDEase2D does automatic error analysis, you do not need to make several runs with different meshes to verify convergence. You can choose whether or not to override the automatic defaults

For numerical simulations, parameters commonly used for ELT were used. A 600 μm laser fiber was considered. Two different vessel diameters were evaluated: 3 mm and 5 mm. Laser power (10 & 15 W), pulse duration (1, 2, 3 s), delay (or off-phase) between pulses (2, 3 and 4 s) and pullback distances (3, 5 and 7 mm) were the variables. In this case, thermal energy is applied along the length of the vein by withdrawing the laser fiber in 3, 5 or 7 mm increments.

For CW mode simulation, laser power (10, 15 W) and pullback speed (1,2 and 3 mm/s) were the variables.

An irregular 10 × 60 finite element grid was used. The time steps were 0.1s. The tolerances used to converge the solution were 10^{-3}. The vein wall thickness was considered to be 0.4 mm. The initial temperature was set at 37°C. The listing of physical parameters used for numerical simulation is reported in Table 1.

The product of pullback rate and power yields the total dose of energy delivered to a vein during treatment. Several studies have suggested that this parameter is a major determinant of treatment outcome [26–28]. In order to compare the results obtained through mathematical modeling, the energy per centimeter (J/cm) was calculated.

## Results

### Light distribution

### Wavelength

### Pulsed mode

### CW mode

## Discussion

To date, mathematical modeling of ELT has never been proposed. This task was performed to assist in providing a better understanding of the ELT process and possibly to determine the optimal dosage as a function of vein diameter. Our model remains a mathematical model, implying that errors may appear owing to the considerations and simplifications required to realize it. Generally, such errors appear because of inaccuracy of the optical, thermal, and damage properties that are critical points in the model's set of equations. In fact, these properties play a key role in the accuracy of the results achieved. Many methods have been presented to calculate these properties but still we see differences in the values presented by the different groups, which reflect the difficulty of measuring these properties. The problem is increased by the reliance of the properties on different variables (temperature, damage) over time. This makes the deviation neither linear nor regular [24].

- 1)
One of the main problems remains the knowledge of the vein diameters during ELT treatment. Since the vein diameter is considerably reduced after tumescent anesthesia, it should be systematically measured. A recent study by Desmyttere et al have demonstrated that after tumescent anesthesia, the vein diameter was usually reduced down to 5 mm or less [27].

- 2)
In contrast to the mode of action of VNUS closure (radiofrequency) where a significant shrinkage of the vessel wall is observed, Proebstle has clearly demonstrated that, when performing ELT, permanent occlusion, reported at 3 months or later, can be obtained by thermal damage of the tunica intima inner vein wall only [29]. This observation is confirmed by the histological study performed by Corcos et al. They showed that when permanent occlusion was observed, the endothelium and intima were always damaged and that success was independent of the vessel wall thickness [30].

The results of mathematical modeling for 810 and 980 nm shows 980 nm leads to a slightly greater damage of the tunica intima when compared to 810 nm. This is owed to its better absorption by blood at 980 nm. However, because of the inaccuracy on vein diameter, one can consider that the choice of the wavelength between 810 and 980 nm has no influence on the results. This is confirmed by the literature where 810, 940, 980 nm were used for ELT with similar parameters [7, 9, 10]. Proebstle et al performed an in vitro study to evaluate the role of intravascular blood for the effective transfer of thermal damage to the vein wall through absorption of laser energy with 810, 940 and 980 nm. Similar results were obtained with 810 nm, 940 nm and 980 nm [29].

When using the continuous mode (810 nm, 14 W), Timperman, has treated 100 veins with an average energy of 95 J/cm (range, 57–145 J/cm; SD: 16 J/cm). Follow-up and success at 1 week was 100%, 96% at 3 month follow-up and finally 95% at 9 month follow-up [31]. Using 58 J/cm, in another series, Timperman reported only a 76% complete vein ablation rate. Similarly, Theivacumar has confirmed that, of all parameters, energy per cm was the main determinant of successful LSV ablation by ELT [26]. The parameters used by Timperman are very similar to those determined by our simulations. In continuous mode, for a 3 mm vein diameter, 50 J/cm are required to damage the vessel wall. For a 5 mm diameter, 100 J/cm are required.

In pulsed mode, it is often difficult to obtain details concerning parameters used for ELT. If power and pulse duration are usually reported, information concerning speed of fiber withdrawal is usually missing. Results reported in the literature are highly variable and tend to prove that this parameter is not well controlled by the operator. In a study performed in 476 limbs (810 nm, 12 W, 1s and a pull back distance varying from 2 mm to 2.8 mm giving respectively 61 J/cm and 43 J/cm), Theivacumar et al observed that occlusion rates were significantly greater at higher energy levels [26]. Using 50 J/cm, in a series of 56 limbs in 41 patients, Mozes et al reported complete resolution at 3-month follow up [11]. Finally, Proebstle et al on a series of 77 patients, performed ELT with a median energy delivery of 23.4 J/cm (range of 11.8 to 35.5) Using these parameters, at 3 months post ELT, 10% of GSVs were found open by color Doppler examination [12]. These results are similar to our calculations in that it was determined that 15 J/cm is required to obtain a permanent damage of a 3 mm vein diameter and 50 J/cm for a 5 mm vein diameter.

An interesting observation is that less energy is required in pulsed mode than in continuous mode because the vein is not heated along its entire length. However, as illustrated in figures 9 and 10, damaging the vein sequentially along its entire length may lead to permanent occlusion.

At last, if steam formation during ELT has already been reported, the interpretation given by the authors is inappropriate [6]. Steam bubbles originating from boiling blood cannot be the pathophysiological mechanism of action of ELT. The steam produced by absorption of laser energy by the blood is a tiny fraction of the energy necessary to damage the vein wall and cannot be the primary mechanism of injury to the vein with endovenous laser. The carbonization and tract within the vein walls seen by histology following endovenous laser can only be the result of direct contact between the laser fiber tip and the vein wall. Venous caliber reduction is maybe due to collagen shrinking by simple heating, but media contraction is only obtained by selective 810–980 nm irradiation.

Consequently, the parameters determined through calculation in our model, especially the ones concerning the different outcome in relation with the different diameters, seem to confirm and explain the observations emerged by the clinical practice

## Conclusion

The parameters determined by mathematical modeling are in agreement with those used in clinical practice. They confirm that thermal damage of the inner vein wall (tunica intima) is required to achieve the tissue alterations necessary in order to lead the vein to permanent occlusion. However, in order to obtain a high rate of success without adverse events, the knowledge of the vein diameter after tumescent anesthesia is recommended in order to use the optimal energy. As clearly demonstrated by our calculations, both pulsed and continuous mode operations of the laser can be efficient. An interesting observation in our model is that less amount of energy is required in pulsed mode than in continuous mode. Damaging the vein sequentially along its entire length may lead to permanent occlusion. However, the pulsed mode requires a very precise positioning of the fiber after each pullback and the duration of the treatment is much longer. For these reasons, continuous irradiation seems to be preferred by most clinicians.

This model should serve as a useful tool to simulate and better understand the mechanism of action of the ELT

## Declarations

### Acknowledgements

The authors wish to thank Pascal Servell for careful reading of the manuscript.

## Authors’ Affiliations

## References

- Min RJ, Khilnani N, Zimmet SE:
**Endovenous laser treatment of saphenous vein reflux: long-term results.***J Vasc Interv Radiol*2003,**14**(8):991–996.View ArticleGoogle Scholar - Beale RJ, Mavor AI, Gough MJ:
**Minimally invasive treatment for varicose veins: a review of endovenous laser treatment and radiofrequency ablation.***Int J Low Extrem Wounds*2004,**3**(4):188–197. 10.1177/1534734604272245View ArticleGoogle Scholar - Morrison N:
**Saphenous ablation: what are the choices, laser or RF energy.***Semin Vasc Surg*2005,**18**(1):15–18.View ArticleGoogle Scholar - Mordon S:
**L'endolaser: une alternative à la phlébectomie ambulatoire.***Réalités Thérapeutiques en Dermato-Vénérologie*2001,**109:**26–31.Google Scholar - Boné Salat C:
**Tratamiento endoluminal de las varices con laser de diodo: estudio preliminar.***Rev Patol Vasc*1999,**5:**35–46.Google Scholar - Proebstle TM, Lehr HA, Kargl A, Espinola-Klein C, Rother W, Bethge S,
*et al*.:**Endovenous treatment of the greater saphenous vein with a 940-nm diode laser: thrombotic occlusion after endoluminal thermal damage by laser- generated steam bubbles.***J Vasc Surg*2002,**35**(4):729–736. 10.1067/mva.2002.121132View ArticleGoogle Scholar - Chang CJ, Chua JJ:
**Endovenous laser photocoagulation (EVLP) for varicose veins.***Lasers Surg Med*2002,**31**(4):257–262. 10.1002/lsm.10103View ArticleGoogle Scholar - Goldman MP, Mauricio M, Rao J:
**Intravascular 1320-nm laser closure of the great saphenous vein: a 6- to 12-month follow-up study.***Dermatol Surg*2004,**30**(11):1380–1385. 10.1111/j.1524-4725.2004.30431.xGoogle Scholar - Oh CK, Jung DS, Jang HS, Kwon KS:
**Endovenous laser surgery of the incompetent greater saphenous vein with a 980-nm diode laser.***Dermatol Surg*2003,**29**(11):1135–1140. 10.1046/j.1524-4725.2003.29353.xGoogle Scholar - Min RJ, Zimmet SE, Isaacs MN, Forrestal MD:
**Endovenous laser treatment of the incompetent greater saphenous vein.***J Vasc Interv Radiol*2001,**12**(10):1167–1171.View ArticleGoogle Scholar - Mozes G, Kalra M, Carmo M, Swenson L, Gloviczki P:
**Extension of saphenous thrombus into the femoral vein: a potential complication of new endovenous ablation techniques.***J Vasc Surg*2005,**41**(1):130–135. 10.1016/j.jvs.2004.10.045View ArticleGoogle Scholar - Proebstle TM, Krummenauer F, Gul D, Knop J:
**Nonocclusion and early reopening of the great saphenous vein after endovenous laser treatment is fluence dependent.***Dermatol Surg*2004,**30**(2):174–178. 10.1111/j.1524-4725.2004.30051.xGoogle Scholar - Weiss RA:
**Endovenous techniques for elimination of saphenous reflux: a valuable treatment modality.***Dermatol Surg*2001,**27**(10):902–905. 10.1046/j.1524-4725.2001.01999-2.xGoogle Scholar - Proebstle TM, Gul D, Kargl A, Knop J:
**Endovenous laser treatment of the lesser saphenous vein with a 940-nm diode laser: early results.***Dermatol Surg*2003,**29**(4):357–361. 10.1046/j.1524-4725.2003.29085.xGoogle Scholar - Sturesson C, Andersson-Engels S:
**Mathematical modelling of dynamic cooling and pre-heating, used to increase the depth of selective damage to blood vessels in laser treatment of port wine stains.***Phys Med Biol*1996,**41**(3):413–428. 10.1088/0031-9155/41/3/006View ArticleGoogle Scholar - Babilas P, Shafirstein G, Baumler W, Baier J, Landthaler M, Szeimies RM,
*et al*.:**Selective photothermolysis of blood vessels following flashlamp-pumped pulsed dye laser irradiation: in vivo results and mathematical modelling are in agreement.***J Invest Dermatol*2005,**125**(2):343–352.Google Scholar - Anvari B, Rastegar S, Motamedi M:
**Modeling of intraluminal heating of biological tissue: implications for treatment of benign prostatic hyperplasia.***IEEE Trans Biomed Eng*1994,**41**(9):854–864. 10.1109/10.312093View ArticleGoogle Scholar - Beacco CM, Mordon SR, Brunetaud JM:
**Development and experimental in vivo validation of mathematical modeling of laser coagulation.***Lasers Surg Med*1994,**14**(4):362–373.View ArticleGoogle Scholar - Glenn TN, Rastegar S, Jacques SL:
**Finite element analysis of temperature controlled coagulation in laser irradiated tissue.***IEEE Trans Biomed Eng*1996,**43**(1):79–87. 10.1109/10.477703View ArticleGoogle Scholar - Iizuka MN, Vitkin IA, Kolios MC, Sherar MD:
**The effects of dynamic optical properties during interstitial laser photocoagulation.***Phys Med Biol*2000,**45**(5):1335–1357. 10.1088/0031-9155/45/5/317View ArticleGoogle Scholar - Zhang R, Verkruysse W, Aguilar G, Nelson JS:
**Comparison of diffusion approximation and Monte Carlo based finite element models for simulating thermal responses to laser irradiation in discrete vessels.***Phys Med Biol*2005,**50**(17):4075–4086. 10.1088/0031-9155/50/17/011View ArticleGoogle Scholar - Chang IA, Nguyen UD:
**Thermal modeling of lesion growth with radiofrequency ablation devices.***Biomed Eng Online*2004,**3**(1):27. 10.1186/1475-925X-3-27View ArticleGoogle Scholar - Shafirstein G, Baumler W, Lapidoth M, Ferguson S, North PE, Waner M:
**A new mathematical approach to the diffusion approximation theory for selective photothermolysis modeling and its implication in laser treatment of port-wine stains.***Lasers Surg Med*2004,**34**(4):335–347. 10.1002/lsm.20028View ArticleGoogle Scholar - Mohammed Y, Verhey JF:
**A finite element method model to simulate laser interstitial thermo therapy in anatomical inhomogeneous regions.***Biomed Eng Online*2005,**4**(1):2. 10.1186/1475-925X-4-2View ArticleGoogle Scholar - Welch A, Van Gemert M:
*Optical-Thermal Response of Laser-Irradiated Tissue.*New York: Plenum Press; 1995.View ArticleGoogle Scholar - Theivacumar N, Beale RJ, Mavor AI, Gough MJ:
**Factors influencing the effectiveness of Endovenous Laser Treatment (EVLT) for varicose veins due to saphenofemoral (SF) and long saphenous (LSV) reflux.**Ireland**TVSoGBa, editor.**In*The Vascular Society Yearbook*. London, UK; 2005:40.Google Scholar - Desmyttere J, Grard C, Mordon S:
**A 2 years follow-up study of endovenous 980 nm laser treatment of the great saphenous vein: role of the blood content in the GSV.***Medical Laser Application*2005,**20:**283–289. 10.1016/j.mla.2005.08.003View ArticleGoogle Scholar - Mundy L, Merlin TL, Fitridge RA, Hiller JE:
**Systematic review of endovenous laser treatment for varicose veins.***Br J Surg*2005,**92**(10):1189–1194. 10.1002/bjs.5142View ArticleGoogle Scholar - Proebstle TM, Sandhofer M, Kargl A, Gul D, Rother W, Knop J,
*et al*.:**Thermal damage of the inner vein wall during endovenous laser treatment: key role of energy absorption by intravascular blood.***Dermatol Surg*2002,**28**(7):596–600. 10.1046/j.1524-4725.2002.01309.xGoogle Scholar - Corcos L, Dini S, De Anna D, Marangoni O, Ferlaino E, Procacci T,
*et al*.:**The immediate effects of endovenous diode 808-nm laser in the greater saphenous vein: morphologic study and clinical implications.***J Vasc Surg*2005,**41**(6):1018–24. 10.1016/j.jvs.2005.03.002View ArticleGoogle Scholar - Timperman PE:
**Prospective evaluation of higher energy great saphenous vein endovenous laser treatment.***J Vasc Interv Radiol*2005,**16**(6):791–794.View ArticleGoogle Scholar - Spinelli L, Torricelli A, Pifferi A, Taroni P, Danesini GM, Cubeddu R:
**Bulk optical properties and tissue components in the female breast from multiwavelength time-resolved optical mammography.***J Biomed Opt*2004,**9**(6):1137–1142. 10.1117/1.1803546View ArticleGoogle Scholar - Cheong W, Welch A:
**Analysis of optical and thermal behavior of laser assisted balloon angioplasty.***IEEE Trans Biomed Eng*1999,**36**(12):1233–1243. 10.1109/10.42118View ArticleGoogle Scholar - Pfefer TJ, Choi B, Vargas G, McNally KM, Welch AJ:
**Pulsed laser-induced thermal damage in whole blood.***J Biomech Eng*2000,**122**(2):196–202. 10.1115/1.429642View ArticleGoogle Scholar - Pearce J, Thomsen S:
*Rate process analysis of thermal damage. Thermal reponse of laser irradiated tissue.*Edited by: Welch A. New York, USA: Plenum Press; 1995.Google Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.