24 January 2000 Blood vessel architectural features and their effects on thermal phenomena
Author Affiliations +
Proceedings Volume 10297, Matching the Energy Source to the Clinical Need: A Critical Review; 1029706 (2000) https://doi.org/10.1117/12.375226
Event: Digital Optical Computing, 1990, Los Angeles, United States
Blood vessels are anatomically and mechanically complex tissue structures. The particular architecture of arteries, veins and capillaries affects their response in therapeutic heating applications. This work summarizes the signal anatomic features of blood vessels for engineers and scientists, introduces the primary underlying physical principles which describe their function and response to external heating for clinicians, and provides a few representative examples of heating applications. It is by no means complete — such a work would be at least two orders of magnitude longer — but, we hope, provides an adequate survey of the pertinent literature so that the interested can locate more substantial treatments of the subject.



Blood vessels are hollow round tubes functioning as conduits for blood flow. Blood, a tissue composed of fluid (plasma) and the red (erythrocytes) and white (inflammatory) cells, is a major player in transport and transfer of heat, gases, nutrients and metabolic wastes in living systems. The circulatory system is designed for heat and mass transfer.

The basic structural anatomy of blood vessel walls is constant but the relative amounts of each component vary considerably among the different vessels to satisfy different design criteria and constraints. The blood circulation begins at the heart that pumps blood into the systemic vascular system via the aorta (about 1 cm in diameter) and the pulmonary vascular system via the pulmonary arteries (about 6 to 7 mm in diameter). The arteries become smaller in diameter and their function changes from a major conduit to smaller delivery pipelines feeding blood to the various organs and tissues, the medium sized arteries. Once these feeder arteries enter into the organs they branch into smaller and smaller arteries to form very small vessels, the arterioles. The arterioles branch further and, ultimately, form the smallest blood vessels, the capillaries (5-10 μm in diameter) through which gases, nutrients and waste materials are exchanged. The blood flow drains from the capillaries into venules that eventually empty into larger and larger veins to connect with the great veins, the superior and inferior vena cavae and the pulmonary veins. Veins in the lower extremities have valve leaflets approximately at 10 cm intervals to prevent back flow under hydrostatic loading and a tunica media as thick as arteries for the same reason. Venous flow in the extremities is driven by muscular contractions. Venous flow in the head and trunk is driven by expansion of the thorax during inspiration. The large veins (1-3 cm in diameter) empty into the heart providing continuity in the flow system.

Histologic studies of vessels from different animals reveal significant species-dependent differences in the spatial distribution of the vessel constituents (vascular architecture) of analogous vessels; and significant architectural differences among different blood vessels in the same species. In other words, the physical and anatomic properties of one blood vessel cannot be applied to all other blood vessels, though there are global similarities.


Tissue Constituents

The basic wall structure of all blood vessels, except capillaries, is three-layered consisting of 1) the intima, the inner layer that lines the vascular lumen, 2) the media, the muscular/connective tissue middle layer and 3) the adventitia, the connective tissue outer layer. The intima is formed of the endothelium, the lining epithelium that forms the luminal surface (about 3 to 5 μm thick), and the underlying myointimal tissue that supports the endothelium. The myointimal layer measures less than 5 μm thick in most normal blood vessels and is formed of myointimal cells intermixed with a felt-work of small collagen fibers. The intima of arteries and veins is essentially the same; but the intima of capillaries consists of endothelium only.

The media is the most important wall component that controls blood flow and blood vessel conformation. To accomplish these functions, the media consists of variable amounts of contractile smooth muscle cells supported by a fibrillar collagenous scaffolding sandwiched between elastin fibers or plates. The arterial media can vary from 1-2 mm (aorta) to 10-15 μm (precapillary arterioles) thick and contains variable ratios of muscle cells, elastin and collagen. The adventitia is the outer most mural layer and is formed of two sublayers. The inner sublayer is a dense, tightly-packed feltwork of collagen and elastin fibers. This layer is tightly adherent to the media and is difficult to dissect away from the vessel. The outer layer is a looser, network of thick collagen fibers that gradually blends with the surrounding tissues. This outer layer can be removed partially by blunt dissection. The adventitia is thicker in arteries than veins and, in general, is thicker relative to the media in larger vessels than in smaller ones. Most of the great arteries also have a network of smaller blood vessels in the adventitia, called the vasa vasorum, whose purpose is to perfuse the vessel tissues1, 2, 3, 4, 5.



Endothelial cells are squamous epithelium specialized for lining the blood vessel. The endothelium is one cell layer thick (1-3 μm) and forms an essentially continuous membrane. All vessels have an endothelial layer, and in capillaries the endothelium comprises the entire vessel wall. Endothelial cells communicate with each other through gap junctions to reform and remodel the vessel walls in response to wall shear stresses; key cells act as shear stress gradient sensors. The cells assume shapes appropriate to local gradients in wall shear, and have particularly interesting morphologies around bifurcations; especially in small arteriole branches where the cells range from 80 to 100 μm long and from quite flat to nearly cuboidal6, as depicted in Figure 1.

Figure 1.

Sketch of representative endothelial cell morphology as typically seen around a branch in an arteriole6.



Smooth Muscle

Vascular smooth muscle cells are the contractile elements that control blood flow and blood pressure gradients by contraction (vasoconstriction) or relaxation (vasodilatation) in response to a plethora of physiologic stimuli and feed-back signals through the nervous system, hormonal and other biochemical signal pathways. Therefore, smooth muscle is the more prominent vessel wall constituent in those arteries and sometimes the veins that control perfusion of a region, organ or tissue. These medium sized arteries are called “muscular” arteries in distinction to the elastin-rich larger “elastic” arteries, such as the aorta, carotid and femoral arteries. Local tissue perfusion pressures are generally governed by the vasodilatation and vasoconstriction of the arterioles; smooth muscle is generally wrapped circumferentially around the outside of an arteriole at strategic locations rather than located diffusely within the vessel wall, as in large vessels. Vascular smooth muscle maintains the mural “tone” of arteries and veins to provide constant perfusion pressure gradients. The vascular smooth muscles are harnessed by collagen scaffolding, elastin plates and fibers to transmit the work of each cell to the vascular wall. The muscle cells are typically oriented circumferentially or spirally relative to the axis of the vessel.



Collagen includes a group of fibrillar proteins, or web-like proteins, that form the major constituent of connective tissue. Examples are: 1) the tendons and fibrous sheaths (fascia) of muscles, 2) skin dermis and 3) the meninges, tissue sheets that cover the brain and spinal cord. The collagens of non-diseased blood vessels form a delicate, fibrillar scaffolding that supports the vascular smooth muscles and interacts with the elastin in the vascular media. The adventitial collagens form thick fibers arranged in a felt-work that encompasses the adventitial elastin then extends to blend with the surrounding tissues. The collagens and elastin are produced by myointimal cells, fibroblasts and myofibroblasts, cells that reside in the vessel wall.

Collagen chemistry is much studied — to date at least fourteen types of collagen molecules have been described. A very effective collection of information on the collagens is the multiple volume series edited by Nimmi, the first volume of which describes collagen biochemistry7.



Elastin is a fibrillar protein that forms dense fibers and flat plates in the media and adventitia of arteries and veins. Elastin’s elastic properties allow arteries and their components to respond to the relentless pounding of pulsatile blood flow in the large arteries and to store the potential energy required to maintain blood flow during diastole. Also, elastin interacts with the more rigid, strong medial and adventitial collagen to maintain the structural integrity of the blood vessels. An almost continuous, distinct elastin membrane, the internal elastic lamina, separates the intima from the media in all systemic and pulmonary arteries. A distinct internal elastic lamina is not present in capillaries and veins. Most elastin is present in the media and the adventitia of both arteries and veins. Arterial medial elastin plates and fibers tend to be parallel to the long axes of the vessel and form more or less regular layers. The medial elastin fibers in veins have a looser, less layered distribution than those of arteries. The elastin fibers of the adventitia form a feltwork layer just adjacent to the media in both arteries and veins. The density, volume fraction and distribution of elastin fibers varies considerably among different blood vessels in the same animal and for the same vessel in different species, and also from point to point in the same vessel. For example, the media of the human aorta and its major branches contain more elastin plates and fibers than the medium-sized arteries entering the organs such as the kidney, heart and gut. On the other hand, the adventitia of medium-sized veins and arteries of pigs contain more elastin than the same vessels in mongrel dogs: the same human vessels have less adventitial elastin than either animal.


Small Vessels

Small diameter capillaries (5 to 10 μm diameter) are the final sites of mass transfer in the body. The three basic types of capillaries are: 1) continuous endothelium with micropinocytosis, 2) fenestrated (i.e. perforated) endothelium — some fenestrae are open while others are closed by diaphragms5, and 3) discontinuous endothelium and/or sinusoids. Non-fenestrated continuous capillaries have endothelial cells held tightly together by gap junctions and desmosomes. Fenestrated capillary endothelial cells have 60 to 100 nm pores, and the variety with diaphragms occurs in tissues requiring rapid fluid exchange (gastro-intestinal mucosa, kidney peritubular plexus cells etc.)5. Larger open fenestrations are typical of capillaries in the carotid and aortic bodies5. The normal pressure difference between the arterial and venous ends of the capillary is around 10 to 20 torr (mmHg). Junctional complexes between the endothelial cells occur at intervals5 in continuous endothelium and permit passive diffusion of fluids, osmotic proteins and nutrients into the extracellular space of the surrounding tissue near the arterial end and collection of fluids, osmotic proteins and waste products near the venous end.

Venules collect blood from capillaries and are of larger diameter but otherwise nearly indistinguishable. Their walls are composed similarly to capillaries consisting of endothelium essentially without a muscular shroud.

Arterioles have endothelial cells surrounded by basement membrane and a spiraling monolayer of smooth muscle cells at varying intervals and locations5. One smooth muscle cell may spiral several times around an arteriole. Contraction narrows the arteriole diameter to control flow. The adventitia of arterioles contains fibroblasts, collagen fibers and sparsely distributed muscle cells. A few delicate longitudinal elastin fibers are distributed just inside the endothelium and its basement membrane (the tunica intima of an arteriole).


Image Orientation and Histologic Stain Characteristics:

  • 1) Weigert-vanGieson elastin stain colors elastin fibers black, collagen red and smooth muscle a light yellow-tan. In the black and white images these correspond to black, dark gray and light gray, respectively.

  • 2) Mallory’s trichrome stain colors elastin pink, collagen bright blue and smooth muscle pale yellow. In the black and white images these correspond to very light gray, dark gray and light gray, respectively.

  • 3) hematoxylin and eosin stains most cellular components a light salmon color and the cell nuclei purple. In the black and white images these correspond to very light gray and dark gray, respectively.

All higher power images have the intimal surface on top and the adventitia on bottom. Collagen native and form birefringence yields bright areas on an otherwise dark field under transmission polarizing microscopy (TPM). Birefringence images in this paper are usually taken from H&E stained sections.


Middle Sized Vessels

Medium sized and smaller arteries and veins tend to be paired or arranged in groups called “vascular bundles”, frequently grouped with a major nerve, as well. The arteries have smaller outer and luminal diameters and thicker, more rigid walls than those of their companion veins. Frequently, veins are more easily recognized than arteries by their larger size and the red blood visible through their thin walls. The artery and nerve in the bundle often have similar size and appearance at first glance.

Figure 2.

The renal vein in the pig at original magnification 40X, (a) hematoxylin and eosin stain (H&E) reveals smooth muscle nuclei (darker spots) distributed throughout the vessel wall and the transmission polarizing microscopy view (TPM), (b) shows the relatively uniform collagen distribution (the bright field).


Figure 3.

The renal artery in the pig at original magnification 6.25X, (a) Elastin stain reveals smooth muscle packed media nearly devoid of elastin with elastin-rich adventitia as a dense feltwork. b) the transmission polarizing microscopy view (TPM) shows collagen loosely distributed in the outer adventitia only.


Figure 4.

The splenic artery in the pig at original magnification 40X. a) hematoxylin and eosin stain (H&E) reveals smooth muscle and myointimal cell nuclei (darker spots) distributed throughout the vessel wall, b) transmission polarizing microscopy (TPM) shows a more uniform collagen distribution (bright field) than the renal artery.


Figure 5.

The splenic vein in the pig at an original magnification of 6.25X: a) H&E stain; b) TPM delineates the adventitial collagen.



Large Vessels

The aorta has substantial structure devoted to storage of blood pulse energy in elastin bands and plates and ultimate strength in a surrounding collagen harness. In figure 7 the aorta in the pig may be seen to be a little thinner in cross-section than the canine aorta of figure 6 (note magnification 3IX compared to 25X for the canine aorta). Aorta biomechanics is an intricate problem combining both fluid mechanics forces and pressure stress forces which are not in phase.

Figure 6.

Canine aorta at an original magnification of 25X. a) Hematoxylin and eosin (H&E) shows elastin plates as medium gray with sparse smooth muscle nuclei (dark spots) at this location, b) Transmission polarizing microscopy (TPM) reveals sparse collagen in the media with a denser feltwork in the adventitia.


Figure 7.

Porcine aorta sectarian at an original magnification of 3IX. (a) Hematoxylin and eosin (H&E) stain reveals smooth muscle nuclei (dark spots), (b) Transmission polarizing microscopy (TPM) of the same section shows a loose collagen feltwork in the adventitia. (c) Weigert-vanGieson elastin stain of an adjacent section illuminates the thin elastin plate arrangement.


Figure 8 illustrates the typical architecture of the canine femoral artery and vein and the pig femoral artery for comparison. The femoral arterial media is separated from the intima by a prominent lamina propria. The media consists mainly of smooth muscle cells supported by a collagen framework and spare, thin elastin plates in the dog and thick well organized plates in the pig. In the dog, the femoral arterial adventitia is composed of a thick elastin feltwork that is almost as thick as the media. The femoral vein has a thinner wall with no internal elastic lamina and ribbons of elastin interspersed with dense collagen fibers and smooth muscle cells. A TPM image (not shown) reveals very little collagen in the media and a dense feltwork of collagen in the adventitia for the pig femoral artery section.

Figure 8.

Canine Femoral Artery (a) and Femoral Vein (b). For comparison, the femoral artery in the pig at 40X original magnification. [Elastin Stain. Bar (a) = 99 μm and (b) = 62 μm]. Note the difference in elastin distribution in the two species: the media has only a small amount of elastin in the dog and is a bit more elastin-rich and plate-like in the pig. The feltwork in the adventitia is more dense in the pig specimen, (c).


Two sections of the canine carotid artery are shown in Figure 9. The specimen shown in Figure 9a was collected near the origin of the carotid artery at the aorta and the specimen in 9b is from the middle of the common carotid. The proximal carotid artery has a thick media composed of prominent elastin plates connected by a collagenous scaffolding in which smooth muscle cells are found. The internal elastic lamina (luminal surface) is thin and mingles with the medial elastin. The adventitia (not shown) consists of a dense feltwork of collagen. On the other hand, the intimal surface of the middle carotid artery is separated from the media by a continuous (except for very small pores) elastic membrane, the internal elastic lamina. The media is composed of sheets of smooth muscle cells separated by irregular elastin plates and ribbons and by collagen. The proximal adventitia is a dense feltwork of elastin fibers that gradually intermixes with the collagen of the more peripheral adventitia.

Figure 9.

Proximal Carotid Artery (a) and Middle Carotid Artery (b) in the dog. For comparison, the subclavian artery in the pig at 40X original magnification (c) has elastin plates in the media and dense but more disorganized elastin sheets in the adventitia. The elastin is black, the collagen, medium gray, and vascular smooth muscle, light gray. [Elastin Stains. Bar = 62 μm in both (a) and (b)]


Figure 10.

Carotid Artery in the pig elastin stain, (a) and TPM view, (b) at an original; magnification of 8X. (c) is the same section at 40X original magnification. The media has elastin plates and dense but more disorganized elastin sheets in the adventitia.


Figure 11.

Iliac Artery in the pig elastin stain, (a), at an original magnification 6.25X. b) the same section at an original magnification of 40X. (c) TPM view of the same section, also at 40X. The media has elastin plates and dense and relatively organized elastin sheets in the adventitia.


The internal jugular vein in the pig (figures 12 and 13) illustrates the complimentary distribution of constituents. The H&E section,(a) reveals muscle cell nuclei, the Transmission polarizing microscopy (TPM) section, (b), shows the collagen distribution and the elastin stain, (c), identifies the elastin distribution. Smooth muscle occurs in the dark gaps in the birefringence image (b). Elastin is primarily distributed throughout the adventitia and in specific plates in the media. In the dog the internal jugular vein is also elastinized, but with a sparse distribution of thin plates, some of which may just be seen in the H&E stained section of figure 14.

Figure 12.

The internal jugular vein in the pig at original magnification 8X. (a) Elastin section shows elastin bands in the media and outer adventitia. (b) Transmission polarizing microscopy shows collagen distribution (bright field) primarily in the adventitia.


Figure 13.

The internal jugular vein in the pig at original magnification 40X. (a) H&E section, (b) Transmission polarizing microscopy shows collagen distribution (bright field), (c) Weigert vanGieson elastin stain gives a better picture of the elastin bands at this magnification.


Figure 14.

The internal jugular vein in the dog at original magnification 40X under H&E staining, a) The section is slightly thinner than in the pig. b) TPM view of the same section (also at 40X) shows a surprisingly uniform distribution of collagen in sharp contrast to the pig internal jugular vein.


Figure 15.

The femoral vein in the dog at original magnification 25X H&E stain, a) is slightly thinner than in the pig. (b) TPM view of the same section (also at 25X) shows a surprisingly uniform distribution of collagen very similar to the internal jugular vein of figure 14.


Table 1.

Dimensional data for great vessels in the pig, average and (standard deviation).

VesselMedia Thickness (mm)Adventitia Thickness (mm)
Aorta0.471 (0.035)0.229 (0.068)
Subclavian A. (thick media)0.430 (0.032)0.045 (0.02)
Subclavian A. (thin media)0.257 (0.016)0.190 (0.037)
Carotid A. (elastic region)0.402 (0.015)0.047 (0.01)
Carotid A. (muse, region)0.0461 (0.022)0.146 (0.01)
Iliac A.0.258 (0.04)0.185 (0.015)
Femoral A.0.274 (0.06)0.139 (0.025)
Mesenteric A.0.186 (0.017)0.095 (0.028)
Splenic A.0.255 (0.018)0.105 (0.02)
Renal A.0.36 (0.05)0.18 (0.05)
Int. Jugular V.0.095 (0.075)0.145 (0.03)

Table 2.

Relative elastin content for great arteries in the pig, average and (standard deviation).

VesselMedia Elastin (%)Adventitia Elastin (%)
Aorta17.6 (2.0)9.4 (2.0)
Subclavian A. (thick media)52.3 (2.7)26 to 54
Subclavian A. (thin media)44.1 (3.2)37.1 (10.8)
Carotid A. (elastic region)56.4 (1.2)21.8 (1.7)
Carotid A. (muse, region)16.4 (1.2)67.2 (2.2)
Femoral A.8.4 (2.1)80.0 (5.1)
Femoral A.22.6 (3.4)83.8 (3.2)
Mesenteric A.11.5 (1.5)61.2 (7.7)
Renal A.7.5 (0.1)47.3 (8.7)

Table 3.

Relative elastin content for great vessels in the dog, average and (standard deviation); ranges are given when several samples from different animals were analyzed.

VesselMedia Elastin (%)Adventitia Elastin (%)
Aorta (distal)21.4 (2.1)16.4 (5.4)
Aorta (proximal)19.9 (5.5)23.0(11.1)
Carotid A.9.2 to 11.4 (2.2 to 3.2)12.4 to 15.3 (2.9 to 3.0)
Femoral A.1.0 to 3.0 (0.4 to 1.7)37.4 to 47.8 (2.7 to 4.3)
Femoral V.*1.2 to 2.8 (0.4 to 1.2)13.4 to 16.0 (1.9 to 2.9)
Common Iliac A.1.6 to 2.2 (0.4 to 0.5)30.8 to 36.0 (3.2 to 4.2)
Renal A.2.0 to 3.0 (0.0)29.6 to 32.0 (3.3 to 8.3)



Interplay among the fluid mechanical and vessel structural forces is a dominant feature of vessel design, and thus is a determining phenomenon for the vessel architecture. The delicate and subtle differences in vascular architecture exist to satisfy required performance criteria at different locations in the circulation and are no accident. Even though the basic mechanical phenomena play little role in therapeutic heating per se, their effect on the architecture certainly is of singular importance, and so it is worthwhile to review them to establish perspective. Beyond that, the thermodynamic and heat transfer characteristics are important in all therapeutic heating applications, and a review is included here, as well.


Fluid Mechanics

Blood flow has been a subject of considerable investigation from the early measurements of Harvey [in the 1600's], to flow studies by Hagen8 and Poiseuille9, to a large body of recent work summarized effectively in10,11,15.

The continuity equation (conservation of mass) can be used to estimate average velocities from flowrates:


where: Q = the volume flow rate of blood (m3/s), u = the speed (m/s) and A is the cross-sectional area of the vessel (m2). A flow of 4 1pm (66.7 cm3/s) in the descending aorta, about 1 cm in diameter (area = 0.785 cm2), has an average speed of 84.9 cm/s. The typical mass flow rate for normal cardiac output in the descending aorta, as in the above example, would be about 70 g/s.

The velocity is not uniform over the cross-section of the vessel: it is maximum at the center and zero at the wall, the so-called “zero slip” condition, due to affinity between blood plasma and the intimal surface. Thus, a velocity gradient always exists at any surface including vessels. A shear stress (force) is exerted on the wall given by:


where: τ is the shear stress (N/m2), μ is the viscosity (N-s/m2) and x is the direction normal to the vessel surface (m). The surface velocity gradient, du/dx at x = 0, is the shear rate The shear force vector acts along the direction of the vessel wall. The zero slip condition is not provable by molecular theory, but there is no experiment which indicates that the zero slip assumption is invalid11.

Blood is nearly a Newtonian fluid in large and medium sized vessels at ordinary shear rates. However, the non-Newtonian and multiphase nature of blood, where the semi-solid cells are suspended in plasma, are often important determinants of its behavior at high shear, especially in vessels smaller than 0.1 mm in diameter. In those vessels, and at locations of high shear rates in larger vessels, the viscosity of blood depends on the shear, and thus it is non-Newtonian. Additionally, the pulsatile nature of blood flow adds significant complication to its description by standard fluid mechanics means.

The characteristics of fluid flow fields change dramatically when viscous forces (equation 2) are more important than momentum forces, flow due to a pressure gradient. The primary parameter used to characterize the flow state is the non-dimensional Reynolds number, the ratio of momentum forces to viscous forces16:


where: ρ = the density, u = the average velocity and d = vessel diameter. Matching the units in this relation requires some care. Reynolds numbers in excess of the “critical” value, about 2300 in cylindrical tubes, indicate that turbulent flow is likely. Low Reynolds numbers indicate that laminar flow dominates. In laminar flow the differential control volumes of fluid move in stream tubes parallel to the vessel walls and the velocity profile is parabolic across the vessel diameter and satisfies the “zero slip” condition at the wall. A thorough treatment of boundary layer theory may be found in the classical text by Schlicting17.

The time-average Reynolds numbers encountered in vessels are all below the turbulent transition range, about 1000 to 3000 in cylindrical tubes. In human aorta the average Reynolds number is in the range of 150 to 30012, but during a pulse the peak Reynolds number may exceed 500. As an example, if we take the simplified estimates from the above calculation for the aorta, u = 84.9 cm/s, d = 1 cm with an estimated viscosity for blood at low shear of 3 mN-s/m2 (estimated from Figure 10.8 of Caro et al11), the Reynolds number would be about 386. The apparent viscosity of whole blood is very sensitive to shear rate, see Caro et al11 and so while the result is near the expected range there is a considerable band of uncertainty. In the right descending coronary artery Re is the neighborhood of 200 to 30013 and in small arterioles it is around 0.016.

Steady flow in small vessels is well described by the Hagen-Poiseuille result when the Reynolds number is around 1 or less:


where: Δp = the pressure difference between the inlet and outlet of the vessel (N/m2), L is the vessel length (m) and d its diameter (m) and Q the flowrate (m3/s).

In large vessels the multiphase nature of blood is significant near the wall where the shear rate is high. S.L. Soo14 has an excellent discussion of multiphase flow fields. Red blood cells near the wall spin and tumble11 in the velocity gradient (equation 2) and this leaves an RBC deficiency near the wall. Small vessels attached to the larger vessel receive blood with a significantly lower hematocrit than that of the free stream blood.

Flow in capillaries is decidedly different in nature. In a capillary vessel the red blood cell completely fills the vessel diameter as the cells squeeze through it. The flow has been described as “slug flow” and the RBCs are subjected to considerable mechanical deformation as they traverse the capillary.

The pulsatile nature of blood flow is important in all vessels. An unsteady flow, such as vessel flow, is generally described by the Womersley number, though it is probably not sinusoidal. The dimensionless Womersley number, α, indicates the expected deviation from a Poiseuille (i.e. parabolic) velocity profile for a sinusoidally varying flow of angular frequency ω (radians/s):


The flow is “quasi-steady” for values of α less than 1.0, the velocity profile is parabolic during the pulse and the magnitude of the flow rate is determined by the instantaneous pressure gradient (equation 4). At higher α the velocity profile is distorted since inertial forces dominate as the flow is accelerated and decelerated11.

In cases where the blood may be treated approximately as a Newtonian fluid the Navier-Stokes equations may be used to calculate the flow field — this nearly invariably requires a numerical solution except in the simplest of geometries. The Navier-Stokes equations include both momentum and viscous forces on a fluid volume and are more difficult to manage as a result. A finite control volume of fluid is subjected to both principal stresses (force per unit area, N/m2) normal to the face, σjj, due to pressure gradients, and shear stresses parallel to each face, τij due to viscous forces (see figure 16).

Figure 16:

Fluid cube subjected to normal momentum stresses, σ, and shear stresses, τ, as described by the Navier-Stokes equations.


There is one equation for each vector axis — Cartesian coordinates are used in this example and gravitational (body) forces have not been included:


or, put into vector form all three equations reduce to:


where: u is the speed in the x-direction, v in the y-direction, and w in the z-direction, u is the vector velocity field, and this formulation is for an incompressible fluid (blood).


Structural Mechanics

One of the most comprehensive and instructive recent treatments of arterial mechanics is that by Humphrey19. There are others, see for example20, 21, 22, 24 which treat a variety of aspects of this intricate problem. These few paragraphs will only serve to very briefly summarize a few of the important points. The reader is commended to more serious literature on this subject.

Strain (i.e. stretch), ɛm (m/m), in materials is the elongation or deformation in the direction of an applied stress (i.e. force per unit area), ɛm (N/m2), please refer to figure 17. In metals the strain seldom exceeds about 0.002, or 0.2%19. The modulus of elasticity, Em (Young’s modulus), relates principal stress and strain:


Figure 17.

Applied force per unit area (stress, σm = F/A) and resulting elongation (strain, ɛm = ΔL/L).


The shear modulus, Gm, relates torsional strain, γm, and applied shear stress, τm:


In the general sense, a solid control volume is subject to at least as many stresses, both principal and shear, as the fluid control volume in figure 16.

An elastic material has constant values for Em and Gm under all conditions and exhibits no hysteresis when deformed. No material is elastic under all conditions — inelastic failures happen in all materials — however, vessels are much more complex than most structural materials in that they are visco-elastic. A visco-elastic material has combined viscous (i.e. flow, creep and stress relaxation while the load is applied) and also elastic properties. Creep is permanent deformation over time in response to the load while stress relaxation is a reduction in stress (applied force) over time at constant stretch. In visco-elastic vessel walls the strains are much larger in magnitude for smaller applied stresses than in structural materials. Substantial force components off-axis with respect to the orientation of tissue features, like collagen fibers and elastin plates, result in complex nonlinear strain tensors (see, for example, Humphrey19).

It is perhaps most helpful to reduce the complete tensor-based discussion to a few simple observations at the expense of mathematical rigor. As a simple model, vessels are approximately “thin walled pressure vessels” (figure 18). A closed vessel (i.e. closed at both ends) of radius, r, and thickness, t, with internal pressure, P will have two principal stress components, a circumferential stress, σC (N/m2) and a longitudinal stress, σL (N/m2) (in the absence of flow) given by:


Figure 18.

A blood vessel as a thin wall pressure vessel at static pressure P (no flow).


Notice that the static longitudinal stress is always smaller than the circumferential stress, and their relationship is fixed at a ratio of 2:1. We might reasonably expect that the vessel energy absorbing and structural components responsible for handling static stress would be spirally wound at a pitch angle near 27 degrees since tan(27°) = 0.5. Of course, vessels in situ are not closed, in the static sense, and so the static analysis doesn't predict where the longitudinal stress comes from. There will be longitudinal stresses applied, even in the very low flow case, due to momentum reaction forces in the curve of the aorta, downstream pressure drops, and (mostly) to wall shear stress when blood is flowing. So, the winding angle will probably be quite different from 27°.

The origin of most of the shear stresses on the vessel wall in situ is the fluid shear described in equation (2), though other effects also contribute. It had been generally assumed that the wall shear stress (WSS), μ ∂u/∂x, was the origin of most of the adverse effects leading to vascular compromise, but new results suggest that it may in fact be the wall shear stress gradient (WSSG), μ∇2u, which is actually responsible for the adverse effects induced by flow irregularities23.

Arteries may be considered, mechanically, as a mixture-composite since they are a solid-fluid mixture and a composite of elastin-collagen-smooth muscle19. Most larger arteries are sized so that the mean wall stress is relatively constant, even in hypertension19. The aorta and carotid arteries are built up of layers of musculo-elastic fascicles, as first suggested by Clark and Glagov25. These may be seen in each of the several arteries shown in figures 2 to 13 in the previous discussion. A fascicle consists of smooth muscle cells surrounded by elastin plates and collagen bundles between the elastin sheets or plates. The number of fascicle layers increases with vessel diameter: approximately a 5 μm increase in wall thickness for each 1 mm increase in artery diameter19. However, as was seen in the above discussion, the distribution of collagen, elastin and smooth muscle varies significantly from vessel to vessel and between species.

Vessels are formed under load. Physiologic loading includes a significant stretch of elastin sheets even at zero pressure. If an artery is transected or resected, its length shortens typically by 30 to 50%. A resected artery still has substantial residual stress. If a ring of resected artery is cut longitudinally, it will open until the net residual stress is zero; this measurement is referred to as the “opening angle”, studied in some detail by Fung and co-workers26, 27, 28, 29 (and several other papers). They concluded that the opening angle is relatively insensitive to the location of single radial cuts and that a single cut relieves almost all of the residual stress — this latter observation was concluded from the observation that the opening angles of single and multiple radial cuts were found to be the same. However, Vossoughi et al30 reported that if one separates an arterial ring into inner and outer portions (media and adventitia), the single and multiple cut opening angles are different for the two segments. And in very recent model studies by Van Dyke and Hoger31, it was shown that Fung’s observation of identical opening angles between single and several radial cuts does not imply the absence of residual stress: the models predict that even if there is residual stress the opening angles are about the same. In any event, resecting an artery results in isovolumic shortening of the elastin (and other constituents) and, perhaps, a change in the local volume fractions, which may influence the results of opening angle measurements.

An obviously significant experiment is to measure arterial diameter as a function of internal static pressure; the rate of change in circumference is directly related to the elastic modulus in that direction. Cox32 collected data from various canine arteries at their in vivo length (a critical point in view of the tendency of arteries to shorten after resection). The resulting pressure-diameter data show a nonlinear response (as would be expected); substantial differences among the vessels, also expected; and a substantial contribution of smooth muscle to the result. Other studies by Cox 33, 34, 35 are similarly illuminating.

The wall stresses in the aorta and other vessels arise from two separate phenomena. The circumferential stress is due primarily to the internal blood pressure and the longitudinal stress primarily from fluid shear as the blood flows. The maximum circumferential stress occurs at peak pressure and the maximum longitudinal shear stress at peak flow — i.e. at two different times in the cardiac cycle. Considering the cardiac beat frequency as one cycle (where the period of the pressure wave corresponds to 360°), Qiu and Tarbell12. have determined the relative phase angles of the two stresses. Defining the phase angle between pressure distention of the vessel, D, and flow induced shear, τ as ϕ(D – τ):


where Q refers to the time of maximum flow. The first term on the right hand side is the phase angle of the vessel flow “impedance”, Z (complex):


assuming that the peak pressure occurs simultaneously with the peak distention of the vessel wall. They found in detailed model studies that at high unsteadiness (i.e. at high Womersley number, a) the typical phase angles were: ϕ(P - Q) = -45° = ϕ(τ – Q), which means that ϕ(D – τ) = -90° in large arteries, and ϕ(D - τ) = - 45° in large veins.


Heat Transfer and Thermodynamics

Thermodynamics is the study of different states of matter — the net result of energy exchange and material transformation — with no real consideration of the rate at which exchange and transformation take place. Heat transfer is primarily concerned with the rate at which thermal energy is exchanged. The language and mathematical formalism are quite different in the two fields and often practitioners of the two professions do not have much to say to each other. In thermal surgery we must be concerned with both points of view.


Heat transfer

Heat transfer is one of the primary functions of the circulatory system; however it turns out that as a practical matter most of it occurs in smaller vessels. There are many texts on the fundamentals of heat transfer. Perhaps a one of the more complete treatments is by DeWitt and Incropera36, though both Rosenhow and Choi37 and Krieth and Bohn38, among many others, are also quite informative and useful. The best treatment of heat transfer in the vascular system is the two volume set edited by Shitzer and Eberhart.40

Within tissues conduction heat transfer (including convection) dominates. Radiation heat can make a significant contribution at surfaces, however. It is, in fact, the signal itself when you apply thermographic imaging to estimate surface temperature. Convection is important at vessel walls and, possibly, on exposed surfaces. Kays and Crawford39 is probably the most comprehensive treatment of convection heat transfer.

Conduction heat transfer through a uniform solid slab of material of thickness Δx with constant temperatures on both faces is:


where: q = heat flux (W/m2), k = thermal conductivity (W/m-°C), ΔT = the temperature difference (°C) across the slab, Δx = slab thickness (m) and the negative sign confirms that heat flux is from high temperature to low temperature. The conduction heat flux is actually a vector field quantity, q, in a general three-dimensional sense, and can be expressed in terms of the local temperature gradient, ∇T:


Convection heat flux, qconv (W/m2), occurs at the surfaces of blood vessels and is typically formulated as a one-dimensional flux:


where: h = the convection coefficient (W/m2-°C), Tw = the wall temperature and T is the “free stream” temperature of the fluid flowing above the surface. A flat plate in room air without any forced air flow has an effective value for h between 5 and 25 W/m2-°C.36 When the air (or other fluid) is moving various functional relations are used to adjust the convection coefficient, h, to include its velocity.

For completeness, surface infrared radiation (at body temperatures primarily between 3 and 14 μm wavelength) can be formulated in a similar fashion:


where: qrad = radiation heat flux to the environment (W/m2), ɛsurf = surface emissivity, σsb = the Stefan-Boltzmann constant (5.67 × 108 W/m2-K4) and θ = the absolute temperature (K). This relation describes only the simplest model of radiation heat to a black body environment at uniform temperature with the simplest of shape and view factors. DeWitt and Incropera36 treat the intricacies of the radiation heat problem in a readable and understandable fashion; but perhaps the most comprehensive treatment is by Siegel and Howell41.

While equation (14) is useful for a surface exposed to the atmosphere, it is not useful within the circulation, except in the largest of vessels. In a vessel the effect of the boundary layer is of prime importance — there is a thermal boundary layer similar, in most senses, to the fluid boundary layer except that it is based on the temperature gradient near the wall, dT/dx, rather than the velocity gradient of equation (2). In laminar flow, the boundary layer fills the entire vessel cross-section and the description of the overall heat transfer is more complex since heat transfer and local fluid mechanics are strongly coupled. Early treatments of this problem are Chato’s analysis for single vessels and for an artery-vein pair in close proximity42 and the study of tumor circulation effects by Chen and Holmes at about the same time43.

Most vessel-tissue and vessel-vessel heat transfer occurs in “thermally significant” vessels40. Thermally significant vessels have lengths similar to their equilibration length. The length of vessel required for equilibration depends on the velocity of the blood, u, the diameter of the vessel, d, the overall heat transfer coefficient, U (essentially, the effective value of h in equation 14) and the density and specific heat of the blood, ρb and cb respectively:


where: Le is in m, ρb is in kg/m3, cb is in J/kg-°C, u in m/s and U in W/m2-°C. Chen and Holmes estimate that the overall heat transfer coefficient may be found from:


where: kt is the thermal conductivity of the tissue and, Λ is, let’s face it, a fudge factor. A useful estimate for Λ is approximately 3 with an uncertainty of a factor of 2 (i.e. Λ is probably somewhere between 1.5 and 6 for most vessels and tissues).

Thermally significant vessels should have diameters in the neighborhood of 300 to 350 μm or so, and lengths on the order of 9 or 10 mm; larger than arterioles and venules but smaller than typical middle sized vessels. Larger vessels (both arteries and veins) have high flow velocities and the residence time of blood control volumes is too short to exchange significant heat across the diameter. Smaller vessels (arterioles, venules and capillaries) have longer residence times, but the blood comes quickly into thermal equilibrium with the surrounding tissues. Blood flow (perfusion) in these small vessels may actually be modeled as an increase in effective thermal conductivity in the direction of blood flow as a result, the so-called k-effective, keff, model. The properties of the several vascular compartments have been determined44 and are summarized in Table 4

Table 4.

Properties of vascular compartments adapted from44. However, please note that the dimensions are relative. For example, capillaries can be longer than 200 μm.

Vessel% of Vascular VolumeRadius (mm)Le (mm)1/Le
Aorta3.305.01.9 × 1050.002
Large artery6.591.540000.05
Arterial branch5.490.53000.3
Terminal branch0.550.3800.1
Not Identified1.00.17591
Terminal Vein3.300.751000.1
Venous branch29.71.23000.3
Large vein24.23.050000.04
Vena cava5.496.251.9 × 1050.002



The thermodynamic state of a material is completely described by specifying any two of the six state variables: temperature, T, pressure, P, density, ρ (or its inverse, the specific volume, v = 1/p), enthalpy, H (J), internal thermal energy, U (J), and entropy, S (J/°C). The state variables are inter-related by Maxwell’s relations. So, for example, the enthalpy is the sum of internal thermal energy and the P-v product (i.e. the “flow work”):


where: m is the mass of material and c is the specific heat. We heat blood vessels and other tissues to achieve some therapeutic change of state: promotion of tissue perfusion, analgesia, desiccation, irreversible alteration of tissue proteins, deactivation of tissue enzymes, and apoptosis are examples. The energy source raises the temperature in a fashion described by the first law of thermodynamics — (energy in - energy out) is stored in the internal thermal energy within the volume:


where ΔT is the temperature change. Energy flows across the boundaries of the control volume by heat transfer, work, mass transfer and phase transformation. We usually express these various terms in the point form, rather than in the integral form, since we are concerned with highly localized phenomena. For example, a first law continuity relation considering only net conduction heat transfer in the integral (volume) form is:


where the left hand side is the rate of increase in the internal thermal energy, the closed surface, Σ, defines the control volume of integration, and dS is the vector surface normal pointing outwards from the surface (figure 19a).

Figure 19.

Energy balance calculation for a finite control volume of tissue defined by the closed surface Σ. (a) Considering conduction heat transfer only, as in equation 20a. (b) The tissue as an open system with blood flow in, Qin (m3/s), and out, Qout (m3/s) contributing flow work and fluid enthalpy.


This relation may be reduced to the point form for numerical modeling purposes by applying the divergence theorem:


In tissues the first law must be augmented to specifically include tissue metabolic heat and the effects of blood perfusion (figure 19b). Blood flow energy terms can be significant even if the vessel is not one of the “thermally significant” vessels:


where the b subscript refers to blood thermal properties and T is the tissue temperature.

The original formulation of the first law for tissues, which included the flow enthalpy, was the so-called “bioheat equation” by Pennes45 in 1948:


where: qgen has been added for this paper and is the volumetric power deposition rate from the heating source(W/m3), qm is the metabolic heat (W/m3), w is the tissue perfusion in (kgblood/kgtissue-s), Ta the arterial temperature, T the temperature of the infinitesimal volume of surrounding tissue and, again, the b subscript refers to blood properties. Pennes results were derived by assuming that the arterial blood comes to thermal equilibrium while in the tissue volume — this is a supportable assumption in arterioles, capillaries and venules; and is useful in analogous heat transfer situations46, 47. The result has proven quite controversial in that gross simplifying assumptions were employed which are not at all compatible with the details of circulatory heat transfer as discussed above. In a seminal series of papers Weinbaum formulated an elegant alternative approach to handle the general case48, 49, but it is difficult to implement unless significant details of the precise distribution of blood perfusion and vessel dimensions are well-defined in the tissue. These issues are of substantial significance when large volumes of tissue are to be heated for analgesia or, especially, in tumor treatment by hyperthermia therapy where elevation to a precise steady state therapeutic temperature range for extended time periods is essential. In tumor heating one must precisely balance applied power deposition with local perfusion and heat conduction to achieve steady state.

For thermal surgery purposes we are applying significant external heating to irreversibly alter tissues after a significant temperature rise. Putting aside the cases of therapeutic heating at low temperatures for analgesia and the like, phase change and power generation by the surgical heat source are the most important terms and dominate metabolic heat and tissue perfusion heat transfer by several orders of magnitude. So:


where: qfg is the phase change heat (expressed as W/m3), and “surface losses” represent radiation loss and convective heat transfer into the blood and/or the atmosphere.


Thermal damage processes

The effect of elevated temperatures is to obtain an irreversible change in the vessel tissues, collectively described here as thermal damage. Kinetic models of thermal damage processes based on an Arrhenius formulation have been used for many years to describe different forms of thermal damage50, 51, 52. In standard thermal damage process descriptions the endpoint damage, Ω, is expressed as the time integral of the damage accumulation rate:


where: Ω is the damage parameter (dimensionless), C the concentration of undamaged tissue constituent, A the “frequency factor” (s-1) is a measure of the collision frequency between molecules, E is an energy barrier which must be overcome by the native state molecules in order to denature (J/mole), R the gas constant (8.3143 J/mole-°K), T the absolute temperature (K), and t is the time (s). Assumptions in this model are that: 1) a single damage process is under consideration (thus one pair of coefficients suffices), and 2) the damage process follows approximately first order kinetics. A discussion of the origins and physical significance of the various components of the damage model may be found in section 4.4 of the Handbook of Thermal Engineering written by Diller, Pearce and Valvano53.

The damage process coefficients, A and E, must be derived from experimental data and are specific to the particular damage endpoint under study. The physical significance of the damage parameter, Ω, is that it is the logarithm of the ratio of original concentration of undamaged material, C(0), to the remaining undamaged material at the end of heating, C(τ). A threshold experiment is one for which Ω = 1; A is then obtained from the intercept and E from the slope on an Arrhenius plot of ln(τ) vs. 1/T for the threshold experiments. A threshold experiment exhibits significant thermal damage since C(τ) = 36.8% — or 63.2% of the vessel wall material has been damaged. Working estimates of representative damage coefficients are contained in Table 5. Additionally, Kang et al62 have measured heat-induced changes in coronary artery wall, but have not extracted A and E values from the data.

Table 5.

Damage model parameters from various sources.

Tissue/ReferenceA(s-l)E (J/mole)Conditions
Takata et al540 T≤ 316 K
 4.3×10644.2×105316 < T ≤ 323 K
 9.3×101046.7×105T > 323 K
Welch & Polhamus553.1×10996.28×105 
Weaver & Stoll582.2×101247.83×105317 < T < 323 K
   T ≥ 323 K
Yang et al593.8×10573.85×105White
Pearce et al611.606×10453.06×105Rat skin


Water vaporization

Water is by far the most thermodynamically active tissue constituent, and for temperatures above about 90 °C water behavior dominates other observed thermal phenomena. Expansion of the water, whether or not phase change takes place, accounts for most of the internal stresses applied to tissue. Above the saturation temperature (100 °C at standard pressure), the vaporization process becomes a volumetric energy sink in which phase change can generate high local pressures and tissue stresses. At extreme heating rates in tissue and/or in free water, such as are realized by extremely short pulsed laser activations, the water may become metastable where the vaporization is rate-limited with rate coefficients (similar to A and E above) determined by the supersaturation ratio, P/Psat, and the surface free energy of nucleated vapor bubbles in the liquid phase. In the rapid ablation case the energy powering the debris plume probably comes from rapid water phase changes.


At lower temperatures water effects are still important since all tissues exposed to air will lose water vapor to the atmosphere. Below the saturation temperature at local atmospheric pressure — i.e. below 100 °C at sea level—water evaporates from exposed surfaces due to the difference between the vapor pressure (the saturation pressure) at the surface temperature and the partial pressure of water vapor in the atmosphere. The vapor pressure at temperature may be found from thermodynamic tables or estimated from a polynomial fit (original data from Keenan and Keyes63):


where: Psat is the saturation pressure (in kPa) and T is the temperature (°C). The relative humidity is the ratio of the actual air water vapor pressure to the saturation pressure of water at the air temperature. Below 100 °C (the saturation temperature) at atmospheric pressure, 100.1 kPa, water vaporization is a diffusion-limited surface loss phenomenon and depends primarily on surface characteristics such as local humidity and temperature-dependent mass diffusion coefficients. A useful means to estimate surface water evaporation is Snelling’s formula64, which was derived from experimental measurement of evaporation rates from solar ponds:


where: ζ is the evaporation rate in meters of water depth/s, u is the free stream velocity of air 7.5 m above the surface, Psat is the saturation vapor pressure of the water at temperature and Patm is the vapor pressure of water in the atmosphere (i.e. the relative humidity times the water saturation pressure at the air temperature). Surface water vaporization typically dominates surface radiation and convection losses and bulk tissue loss terms above about 60 °C65.

Boiling processes

In order to include the effects of boiling in the energy balance, the left hand side of equation (23) must be modified. Neglecting surface losses and allowing changes in pressure and density of the tissue in the enthalpy on the left hand side we obtain:


where the specific enthalpy, h (J/kg), includes internal thermal energy and flow work and qfg is the energy of phase change due to boiling (or evaporation) in (W/m3). This model may be recast in terms of the rate of local temperature rise (as in equation 23) to facilitate numerical implementation if two thermodynamic parameters can be evaluated, ∂P/∂T and ∂ρ/∂T:


where: ∂m/∂t is the mass rate of phase change per unit volume of tissue (kg/s-m3) and hfg is the phase change enthalpy (J/kg). The saturation density curve, a fifth order polynomial fit from 20 to 374 °C is approximately (original data from Keenan and Keyes63):


where: ρsat is in (kg/m3) and is T the temperature (°C). The required term ∂ρ/∂T can readily be calculated from this relation. Similarly, from the same source, the phase change enthalpy as a function of temperature is (original data from Keenan and Keyes63):


where: hfg is in (J/kg) and T is the temperature (°C).

In free liquid phase water which is not constrained by rigid or visco-elastic tissue structures, ∂P/∂T is negligibly small except for extremely short laser activations or electrosurgical arcs. The density change term cannot be neglected in free water, however. This additional term (when phase change is included) accounts for the kinetic energy imparted to the plume when water vaporization is rapid and significant. It also contains the origin of acoustic waves in tissue whether or not phase change is significant. For no vaporization, an estimate of ∂ρ/∂T can be obtained from the saturation density vs. temperature curve under the valid assumption that the liquid phase water density is primarily determined by temperature rather than local pressure.

A model based on the energy balance of equation (28), should use a small time step so that the contributions of local heat transfer may be calculated and summed with the volume power term. When the vaporization process is not rate-limited, excess energy above that required to maintain a constant temperature may be imagined to go into the vaporization process at the local temperature and pressure. The effect of vaporization on the local pressure may be calculated if an acceptable model for tissue mechanical response is available. Note that at 120 °C the saturation pressure (equation 25) is already up to 2 atmospheres (200.2 kPa) and climbs precipitously above that. When external energy, qgen, is added faster than evolved water vapor can diffuse through tissue to maintain liquid saturation (100 C) and atmospheric pressure (100.1 kPa), the tissue temperature readily exceeds 100 °C, which increases the local tissue pressure (thus the saturation temperature) and lowers the phase change enthalpy (equation 29). In a numerical model of slow heating one may usually assume that vaporization occurs at local saturation conditions (equilibrium thermodynamics applies) with rate determined by excess heating energy (that is, energy in excess of that consumed by conduction heat transfer and other losses) and thermodynamic properties.

The mechanical constraints on intracellular water are significant. Extracellular expanded water and water vapor needs only to “percolate” through spaces in the gross structure while intracellular water vapor expands cell volume putting stress on the cell membrane. The applied stresses, both intracellular and extracellular, create a thin wall pressure vessel problem. In many cases the local pressure increases cause rupture in the tissue structure, histologically observable as the “popcorn” effect — communicating vacuoles of thin remnants of cell structure in large void spaces. Elevated tissue pressures also explain the readily observed tearing phenomena in highly layered tissues. These steam vacuoles are sometimes mis-identified as adipose tissue in histologic sections owing to the similarities in appearance.

Rate-limited vaporization: nucleate boiling

The problem of extremely rapid boiling has been well studied66, 67. Under ordinary conditions there is a plentiful supply of nucleating centers for boiling to initiate and thus the equilibrium boiling approach is sufficient for most purposes. In the case of extremely short activations (lasers or RF power) in free water, or in the case of a laser source for which the primary chromophore is water, the boiling may happen in the absence of a nucleating site for the bubbles. This is described as nucleate boiling in which the external energy is used to create a bubble in uniform liquid phase water. To form a critically sized bubble a sufficient number of liquid phase molecules must jump to vapor phase. A bubble of critical radius, rcrit, has a 50% probability of growing or of collapsing back to liquid phase. This radius is small enough that the energy stored in surface tension is quite large:


where: σfg (N/m) is the surface tension at the liquid temperature and the liquid pressure, P, is less than Psat at the liquid temperature because the liquid phase is metastable (i.e. super-saturated). Two energy barriers must be overcome, similar to the one surmounted in the Arrhenius damage model of equation (24) to form a critically sized bubble. A liquid phase water molecule must jump the kinetic barrier to vapor phase (hfg) and also overcome the surface tension of the bubble-liquid interface. There are thus two exponential energy barriers in nucleation theory:


where: Jnucl. = the rate of formation of critically sized bubbles (number/s-m3), Z1 = number of molecules/unit volume = ρ NA/(0.018) for water and NA = 6.023×1023, M is the number of moles per unit mass, m is the mass per molecule = 18/NA grams for water, and b is the super-saturation ratio: b = (PSat - P)/Psat. Only a few degrees of superheat are necessary to have extremely high rates of critical bubble formation. Once a nucleating center is formed it grows to larger size depending on the relationship between the internal bubble pressure and the saturation pressure of the liquid phase at its temperature.



Tissues in general, and vessels in particular, are both semiconductors and lossy dielectric materials. At microwave frequencies (between about 500 MHz and 10 GHz) the dominant power deposition mechanism is by losses during dipole rotation of water molecules in the microwave electric field. The ISM (Industrial, Scientific and Medical) frequencies of practical use in this range are 915 MHz and 2.45 GHz. At radio frequencies (6.78, 13.56 and 27.12 MHz) and below, the dominant heating mechanism is by “Joule” or resistive heating. The frequency range typical of electrosurgery is about 500 kHz to 5 MHz, and all such applications are likewise governed by resistive heating. Heating by direct absorption in a magnetic field is not important in tissues since, for all practical purposes, the tissues have the same permeability as free space. Yet, tissues heat well in a radio frequency magnetic field. This is because the radio frequency magnetic field induces an electric field (by Faraday’s Law of Induction) and the tissues are actually heated by the dominant electric field loss mechanism.


The Volume Power Generation Term

From inspection of the Poynting Power Theorem (see for example Paris and Hurd68, or Roussy and Pearce69), the divergence of the power density in an electric field expressed in the frequency domain is:


where: S = power density (W/m2) in the electromagnetic wave, E = electric field vector (V/m), J = current density (A/m2), ω = angular frequency (r/s), D = electric flux density (C/m2), σ = electric conductivity (S/m) and ɛ = electric permittivity (F/m). The power density has both real and imaginary parts, of which the real part represents heating and the imaginary part energy storage (for wave propagation). The volume generation heating term, qgen, is:


where: ɛr” is the “loss factor”, the imaginary part of the relative permittivity, and ɛ0 is the permittivity of free space, 8.854 × 10-12 (F/m). The first term represents Joule (resistive) heating and the second term lossy dielectric heating. The power generation term for equation (23) may be calculated when the respective properties and local electric field are known.

At lower radio frequencies typical of electrosurgery the electric field may be estimated by applying the quasi-static assumption. The maximum dimensions of any practical problem are much shorter than the wavelength in tissues; thus, the Laplace equation may be solved subject to appropriate boundary conditions to estimate the quasi-static potential field, V(x, y, z):


and the electric field estimated from the gradient of the potential:


At microwave frequencies it is usually necessary to solve the complete set of Maxwell’s equations to determine the electric field in tissues. Methods for estimating electric field distributions in microwave fields have been extensively studied; a description and critical review of the various techniques may be found in Roussy and Pearce69.

As a practical matter, numerical models must be used to study heating in vessels because of their complex structure and response. It is important to keep in mind that numerical models are almost never able to accurately predict specific experimental results owing to uncertainty in physical properties and randomness in tissue structure. However, they do permit study of the governing processes on fine spatial and temporal scales which cannot be obtained in experiments. When numerical models are sufficiently well behaved they can be used to study transient micro-structural thermodynamic events, and to compare dominant processes and effects. They can be used to provide substantial insight into the clinical problem.


Electrical Permittivity and Conductivity

At microwave frequencies the electrical property of interest is the loss factor, which is dominated by the tissue water. No published values of loss factor measurements specific to vessels could be located. Some data for other tissue types are summarized in Table 6. Useful working estimates of the loss factor may be obtained from simple volume fraction mixture formula:


where: m refers to the mixture, ϑ is the volume fraction of material 1 dispersed uniformly in material 2 and ɛ* is the complex permittivity, ɛ* = ɛ' - jɛ”. Vessels are typically about 55% by volume water and the loss factor for water over these frequencies is ɛr” = 16 to 18. Losses at microwave frequencies in the non water parts of the vessel are relatively small because the protein molecules are too large to rotate in the electric fields. So, other losses may be ignored in calculations with little error.

Table 6.

Relative electric permittivity71 and conductivity72 of representative tissues.

TissueConductivity (S/m)Permittivity, ɛ'Loss Factor, ɛ″
Fat0.05 to 0.093.9 to 7.20.67 to 1.4
Muscle0.1 to 1.045 to 4813 to 13.9
Bone0.006 to 0.054.2 to 5.8a0.7 to 1.3a


Measured in bone marrow.

At radio frequencies the heating is dominated by resistive heating for which the electrical conductivity is the parameter of interest. The electrical conductivities of representative vessels and model tissues were estimated from tetrapolar impedance measurements in an approximately uniform electric field and are described by Pearce et al70 with properties of elastin and collagen determined by measurements on ligamentum nuchae and dura mater, respectively. Details of the measurement methods may be found in the reference70. The results are summarized in Table 7.

Table 7

Longitudinal and transverse electrical conductivity at room temperature. Values are averages of five samples and (standard deviation)70.

Electrical Conductivity (S/m) (S.D.)
Tissue500 kHz2 MHz
Canine Femoral Artery (trans)0.75 (0.03)0.85 (0.05)
Canine Femoral Artery (long)0.64 (0.12)0.70 (0.11)
Canine Femoral Vein (trans)0.57 (0.10)0.60 (0.11)
Canine Femoral Vein (long)1.04 (0.18)1.05 (0.17)
Bovine Dura Mater (trans)0.26 (0.05)0.27 (0.05)
Bovine Dura Mater (long)0.41 (0.17)0.43 (0.17)
Bovine Ligamentum N. (trans)0.67 (0.02)0.68 (0.02)
Bovine Ligamentum N. (long)0.99 (0.14)1.00 (0.14)

There is no significant difference between electrical conductivities at the representative electrosurgical frequencies of 500 kHz and 2 MHz, as was expected. As reflected by the standard deviations, the uncertainty in the transverse measurements is uniformly smaller than for the longitudinal measurements, as expected from preliminary estimates. Unexpectedly, there is no significant difference between the longitudinal and transverse conductivities of canine femoral artery. The other tissues studied exhibited significant anisotropy, with the longitudinal conductivity typically about twice that of the transverse conductivity.

The electrical conductivity of a semiconductor is the product of its charge per unit volume and mobility of the charge carriers:


where: p is the number density of charge carriers, number/cm3, q is the charge of an electron (C), and μp is the mobility (cm2/V-s). The mobility of charge carriers in solution increases with increasing temperature with the result that the conductivity increases about 1.5%/°C temperature rise [reference71 and others], so that:


where σ0 is the conductivity at temperature T0. Failure to account for even the difference between room temperature and body temperature can result in significant measurement errors (about 20%).



Another important tissue heating technology is laser heating. In contrast to electrical properties, the optical properties of many tissues vary over several orders of magnitude throughout the visible and infrared wavelengths. This may be a little surprising since the index of refraction is, in fact, the square root of the electric permittivity relative to water and since the previous discussion on microwave properties is water-dominated. However, at optical frequencies the energy modes which interact with the irradiation are vibrational rather than rotational, and many more of the cellular constituents participate in the absorption process. Interacting cellular constituents are “chromophores” and they absorb light wavelengths to widely varying degrees. Also, in contrast with radio and microwave frequencies, scattering is often a dominant process at optical frequencies. The relative importance of scattering and absorption is extremely sensitive to wavelength and accurate calculations of light propagation in tissue are difficult.


The Volume Power Generation Term

A few underlying concepts are important in describing laser-tissue interactions. The radiant intensity of a light beam is the power in the beam per unit solid angle (W/sr) at the angle at which it is measured — the beam irradiance (W/m2) at the tissue surface is sometimes referred to as the beam intensity, but this is incorrect. There are 4π ster-radians (sr) in a sphere — a ster-radian is a spherical version of the planar radian where there are 2π radians in a circle. The radiance at a point in tissue, L, is the radiant power per unit area per unit solid angle in the direction of the light beam (W/m2-sr). The light energy which heats the tissue may come from any angle, so the “fluence rate”, Φ (W/m2) in the tissue is the sum of all the radiance at the point, (x,y,z), in the tissue:


where: Ω denotes the solid angle over which the integral is calculated.

The volume power generation term is obtained from the local fluence rate, Φ (W/m2), and the absorption coefficient, μa (m-1) is:


An interesting phenomenon of importance in vessels is the effect of light scattering. Calculating the local fluence rate can be quite complicated when scattering is important. Perhaps the most comprehensive treatment of the physical aspects laser-tissue interaction is that of Welch and vanGemert73. Wolbarsht74 is another useful reference which contains, in addition, a thorough description of the clinical applications of lasers. There are many other treatments of the topic which can be obtained from this very meeting. The most common and effective approach to fluence rate estimation is numerical modeling since the calculations are intricate and tissues are inhomogeneous and anisotropic.

The most common application of optical measurements in clinical medicine is the use of light absorbance to determine blood chemistry. The solutions are very dilute and have been treated with various enzymes etc. to affect the optical properties in a beneficial way. The fluence rate exiting the cuvette is well described by the familiar Beer’s law if scattering and specular reflection are negligible:


where: propagation is in the z-direction, fluence rate is measured for some finite thickness z, and Φ(0) is the immediate sub-surface fluence rate. It is tempting to include the scattering effects of tissues by applying an effective “extinction” coefficient, μeff, in Beer’s law:


where: the effective attenuation coefficient, μeff2 = 3μas + μs(1-g), comes from a solution to the diffusion approximation, μs describes attenuation due to scattering, and the “anisotropy factor”, g, is commonly used to describe the expected cosine of the refraction angle73. If the light is completely scattered in the forward (z) direction g = 1, while if the light is isotropically scattered uniformly in all directions g = 0. The formulation of equation (43) can be effective if the directional nature of the scattered light is included. Scattering is primarily a refraction event and the angle of refraction depends on the angle of incidence relative to the orientation of the scattering center. In standard integrating sphere measurements it is not possible to directly measure μs. The measurements yield the “reduced scattering” coefficient, μs', one of the terms in μeff:


When scattering is important the total reflectance from a surface consists of specular reflectance from the surface (angle of incidence equals angle of reflection) plus diffuse reflectance — light radiance which propagates into the tissue, undergoes several refraction events and emerges from the surface in random directions. Scattering also extends the lateral propagation of the incident light, and may result in effective radii up to several times that of the incident light beam.

At a wavelength for which scattering is a dominant propagation effect, the fluence rate just below the surface may be as much as two to three times the fluence rate at the surface if scattering were not important for the same surface irradiation. We will see the importance of this effect in section 7.2, where a specimen of aorta was irradiated with an Argon laser at 514 nm, a highly scattered wavelength. Even larger factors can be obtained when hollow organs, such as the bladder or an air- or gas-filled blood vessel, are irradiated due the “integrating sphere” effect.


Optical Properties of Bovine Aorta

At both infrared and ultraviolet wavelengths light propagation is dominated by absorption effects. At most of the ultraviolet wavelengths the primary absorbing chromophores are apparently cellular organelles, perhaps the mitochondria. Certainly, absorption by cellular water is the dominant mechanism for wavelengths of interest longer than about 1.5 μm. The water absorption spectrum has a peak at λ = 1.9 μm (μa = 10 mm-1) and another peak just under 3 μm, at about λ = 2.9 μm (where μa = 1000 mm-1)73. Absorption at longer wavelengths is almost invariably determined by water absorption; and so the water content of the tissue dominates the absorption coefficient.

Table 8.

Optical properties of bovine aorta at selected wavelengths. Data for wavelengths shorter than 2 μm are interpolated from figures in Chan et al75. These data were collected over pressures from 0 to 200 kPa and the range of values is given at each wavelength. Both optical properties increase with increasing pressure. Longer wavelength absorption values are estimated from laboratory experience and tissue water content.

WavelengthLaser TypeAbsorption μa (mm-1)Reduced Scattering μs' (mm-1)
400 nm 0.3 to 0.63.5 to 9.5
515 nmCu Vapor0.1 to 0.23.5 to 7.5
514 nmArgon0.1 to 0.23.5 to 7.0
632 nmHe Ne0.1 to 0.23.0 to 5.5
800 nmDiode0.08 to 0.152.5 to 4.2
1.06 μmNd:YAG0.1 to 0.152.2 to 4.0
1.45 μmAbsorption peak0.8 to 2.21.5 to 3.5
2.01 μmTm:YAG5.62 
2.09 μmHo:YAG2.82 
2.94 μmEr:YAG960 
10.6 μmC0267 



The thermal properties required to compute the energy balance (equation 23) are the thermal conductivity, the specific heat, and the density. Tissues are in large part made up of water, and the thermal properties differ little from those of water for most tissues. Perhaps the most complete collection of thermal property data is that of Bowman, Cravalho and Woods76. Of all the data collected in that paper only a few of the measurements are germane to blood vessels. The most comprehensive collection of blood vessel data is by Valvano and co-workers77,78,79.

The density and specific heat always appear as a multiplicative pair. Given k and α the ρc product may quickly be found. The usual format for thermal property data collection is in terms of the thermal conductivity, k, and the thermal diffusivity, α:


Several empirical relationships are of use in estimating thermal properties in the absence of data. Spells80 determined in biological solutions that:


where: k is in (mW/cm-°C), and mw* = the mass fraction of water. Cooper and Trezek81 found an analogous relationship for tissues:


where: k is again in (mW/cm-°C), ρ is the density (g/cm3), mp* is the mass fraction of protein, and mf* is the mass fraction of fat and the coefficients have the required units. They report similar correlations for specific heat (J/g-°C) and density (g/cm3):


Table 9.

Thermal properties of vessels and related tissues near body temperature. Numbers in parentheses are standard deviations.

TissueRef.k (mW/cm-°C)α × 103 (cm2/s)
Normal aorta774.76 (0.41)1.27 (0.07)
Fatty plaque774.84 (0.44)1.28 (0.05)
Fibrous plaque774.85 (0.22)1.29 (0.03)
Calcified plaque775.02 (0.59)1.32 (0.07)
Dog Fern. A.773.91 
Dog Car. A.774.48 
Bovine Aorta, Axial784.91 (0.59) 
Bovine Aorta, Radial783.91 (0.13) 
Muscle762.8 to 5.3 
Fat761.3 to 3.7 
Bone (in vivo)768 to 30 
Plasma795.70 (0.10) 
Whole Blood794.92 (0.09) 
Water765.91 to 6.66 

Correlations for elevated temperatures are also reported by Valvano et al77. For canine femoral artery:


and for canine carotid artery:


where: k is in (mW/cm-°C), and α is in (cm2/s).



Since the early 1930's electrosurgical current at radio frequencies (RF) between about 500 kHz and 5 MHz has been used to cut fuse, seal or coagulate resected and transected blood vessels. Sealing and fusion of vessels by electrosurgical current are strongly influenced by the inhomogeneous architecture of the tissue constituents, particularly in the large arteries. Inhomogeneities in electrical properties of the constituents, smooth muscle, collagen and elastin, lead to sharp spatial gradients in volumetric power deposition which result in uneven heating. The mechanical properties of the various tissue constituents are also of considerable importance. Collagen alteration is apparently an important prerequisite to successful tissue fusion84, 85, 86, 87, 88. Electron microscopic (EM) studies suggest that the fusion process in blood vessels may be dominated by random re-entwinement of thermally dissociated adventitial collagen fibrils (Type I) during the end stage heating and early cooling phases87.


Bipolar Electrode Fields

The model space consists of the first quadrant (I) of the symmetrical system composed of two bipolar flat plate electrodes clamped transversely around a vessel of arbitrary thickness and composition, as shown in Figure 20. The electrodes are of arbitrary size and potential and an insulated thermocouple (t/c) of arbitrary size can be included. Electrical boundary conditions are: 1) zero flux at x = 0, 2) V = 0 at y = 0, 3) constant potential under the electrode and zero flux outside the electrode at y = maximum, and 4) either a mixed boundary condition (simulating leakage flux) or a zero flux boundary at x = maximum. The model space used in the model was 81H X 21V nodes capable of 9 tissue layers. The tissue layers may have arbitrary thickness and electrical conductivity and can be represented by a selectable number of node spaces. Interfaces between tissues are constrained to lie along node planes. Calculations of electrode current and power are on a per unit length basis along the Z-axis.

Figure 20.

Idealized Cartesian 2-D model space of the first quadrant of bipolar electrodes (cross-hatched) clamped symmetrically around a vessel. Thermocouple (t/c) may be included.


The electrical conductivities used in the calculations were: 1) for collagen, σc = 0.26 S/m (transverse), 2) for elastin, σe = 0.67 S/m (transverse), 3) for smooth muscle no published values are available — an estimate was formed using published values for cardiac muscle (values range from 0.11 to 0.47 S/m72) the geometric mean of σm = 0.23 S/m was used for these calculations. Transverse conductivities were chosen because the primary electric field direction was transverse to the fibers and plates. Cardiac muscle was thought to be a more effective electrical analog of smooth muscle than skeletal muscle since it is essentially devoid of the cascade of collagenous connective tissue shrouds which electrically and mechanically isolate the skeletal muscle fibers. Both models have an “intima” 10 μm thick (actual dimensions in vessels are closer to 1 to 3 μm) modeled by muscle properties and an “internal elastic lamina” 20 μm thick (close to the actual dimension) for which elastin properties were used. For both models the electrode was 0.5 mm thick and 4mm across (assigned nodes 1 to 41) with no specific thermocouple site. The major difference between the models is the location of the elastin plates: in the media for the carotid artery and in the adventitia for the femoral artery.

Figure 21.

Vessel model dimensions in the 81 X 21 model space for a) canine carotid artery and b) canine femoral artery. Tissue constituents are collagen (c), elastin (e) and smooth muscle (m). Both models have an intima 10 μm thick (bottom) and an internal elastic lamina 20 μm thick just above it (dark shade). Specific dimensions are in the text.



Canine carotid artery model.

Specific dimensions in the carotid artery model are: 1) 0.38 mm thick overall with 2) an adventitia (collagen mat) 0.17 mm thick, assigned 4 node spaces, 2) four pairs of elastin/muscle layers, each elastin layer 10 μm thick and each muscle layer 50 μm thick, assigned 2 node spaces for each layer. The potential field is shown in figure 22a while the volumetric power density (W/mm3) is in figure 22b. The maximum of the power density was to 27 W/mm3, and the gray scale covers 0 to 10 W/mm3. The electrical model converged to a maximum electrode current change of 0.001% in 10 iterations after 3460 iterations. For a single electrode potential of 50 volts (100Vrms bipolar) the current per unit length was 84.4 mA/mm and resistance per unit length was 592 Ω/mm. For 7 mm of vessel in contact with the electrode, this would correspond to a total current of 591 mA and inter-electrode resistance of 85 Ω, both reasonable values in the light of experimental results.

Figure 22.

Geometrically scaled carotid artery E-field model results, (a) Potential field (0 black to 50 V white, interpolated), (b) power density (0 black to 10 W/mm3, white, no interpolation used).


The voltage gradient is very strong at the edge of the electrode, as expected, with a correspondingly high power density there. The high edge power density is moderated by strong heat transfer effects in the thermal model, figure 23. The variation in electrical properties of the elastin and muscle in the media is plainly visible in the power density field. The heating time was 4 s. (current adjusted to 42.2 mA/mm) with a 1 s cooling period, Figures 12 and 13, and required 9 hours to calculate (for the same time step). Tissue thermal properties were assumed similar to those of water and the vessel is assumed thermally linear homogeneous and isotropic. Water vaporization (i.e. boiling) was assumed to occur at rates low enough that no significant tissue pressure increases resulted. Consequently, all vaporization processes occurred at or below 100 °C; all qgen in excess of that required to supply local heat transfer processes was assumed to apply to local equilibrium boiling at 100 °C. Tissue was allowed to increase in temperature above 100 °C only if fully desiccated.

Figure 23.

Geometrically scaled results for the carotid artery model at 4s. a) Temperature field (0 °C black to 100 °C white, interpolated), b) remaining tissue water distribution (range 0%, black, to 55%, white, by volume).


The thermal conductivity turns out to be a particularly strong influence on the results in the vertical direction since the vessels are very thin and the thermal mass of the electrode is significant. Thermal conduction in the horizontal direction is negligible. Tissue density is very close to that of water and was assumed constant.


Canine femoral artery model.

Specific dimensions for the femoral model artery are: 1) 0.46 mm thick overall with 2) an adventitia composed of four pairs of elastin/collagen layers, each elastin layer 10 μm thick and each collagen layer 50 μm thick (except for the upper-most collagen layer which was 80 mm thick), assigned 2 node spaces for each layer, and 3) the media (muscle) 0.22 mm thick, assigned 4 node spaces. The electrical model converged to a maximum electrode current change of 0.001% in 10 iterations after 1660 iterations. For a single electrode potential of 50 volts (100 v bipolar) the current per unit length was 62.1 mA/mm and resistance per unit length was 805 Ω/mm, higher than the carotid artery model. For 7 mm of vessel in contact with the electrode, this would correspond to a total current of 435 mA and inter-electrode resistance of 115 Ω, both reasonable values in the light of experimental results, and following the trend observed in the experiments.

The potential field is shown in Figure 24a while the volumetric power density (W/mm3) is in Figure 24b. The maximum of the power density was to 15.6 W/mm3, and the gray scale covers 0 to 10 W/mm3. Both the potential and power density fields have a markedly different appearance than for the carotid artery model.

Figure 24.

Geometrically scaled femoral artery E-field model results, a) Potential field (0 black to 50 V, white interpolated) and b) power density (0 black to 10 white W/mm3, no interpolation used).


Again, the voltage gradient is very strong at the edge of the electrode with a correspondingly high power density located there. The variation in electrical properties of the elastin and collagen in the adventitia is not visible in the power density field, however, unlike the carotid artery results.

Figure 25.

Geometrically scaled results for the femoral artery model, (a) Temperature field (0 to 100 °C, interpolated), (b) remaining tissue water distribution (range 0%, black, to 55%, white, by volume).



Fusion seals in larger arteries

Sealing processes in larger arteries are complicated by the vessel mechanics. Heat damaged collagen and elastin shrink in length and swell in thickness. As in laser anastomosis, tissue fusion attains a primary bond by heating collagenous structures to temperatures sufficient to unravel the fibrous array. The fibers re-entwine at lower temperatures during the cooling phase. Re-entwinement is random and either of low probability or of fundamentally differing bond type since the fused tissue is of low tensile strength compared to undamaged tissue in situ. The temperature history must not create brittle collagen at the interface by excessive desiccation in order for the fusion process to be successful.

Figure 26.

Radio frequency current sealing in carotid artery in the dog at 20X original magnification. Large arrows indicate edges of bipolar electrodes, small arrows indicate steam vacuoles, a) successful seal, H&E stain, b) failed seal, elastin stain.



Vaporization in the Aorta

The effect of water vaporization in highly compartmentalized vessels like the aorta is to create a field of steam vacuoles. Figure 27 is the result of an approximately uniform electric field with current density high enough to cause severe boiling in a post mortem sample of human aorta. The pockets created by the elastin plates of the media confine evolved water vapor and result in the steam vacuoles in the figure. The volume power generation term was larger in the elastin rich media than in the surrounding muscular media owing to differences in the effective electrical conductivity.

Figure 27.

Human aorta damaged by uniform current density sufficient to reach boiling temperatures. Intimal surface is at the top. (a) At original magnification 100X the vacuoles are confined to the elastin-rich media where the plates have confined the vapor resulting in severe vacuolization at high boiling rates, (b) At original magnification of 400X the location of the vacuoles between elastin plates may be seen. [Mallory’s trichrome stain]



Mesenteric Artery-Vein Coagulation

Sealing of a large tissue bundle, such as the mesenteric artery-vein pair shown in figure 28, is not always successful. In the figure, the tissue was damaged to coagulation levels at the edges of the bipolar electrodes but not across the cross-section, so no seal was obtained. High electrode pressure and temperatures often in excess of 100 °C are required to establish a reliable seal.

Figure 28.

Mesenteric artery-vein pair in the pig after heating with RF current from bipolar electrodes. Current concentrations at the electrode edges have damaged the collagen connective tissue (darkened masses) but no seal was accomplished in any of the vessels. Weigert-vanGieson elastin stain.





Laser Vascular Anastomosis.

Electron microscopic (EM) studies suggest that the fusion process in blood vessels may be dominated by random re-entwinement of thermally dissociated adventitial collagen fibrils (Type I) during the end stage heating and early cooling phases91. At the light microscopic level, this bonding process is reflected by the formation of an amorphous coagulum of thermally coagulated adventitial collagen at the anastomotic site.

The temperature history must not create brittle collagen at the interface by excessive desiccation92 in order for the fusion process to be successful. Feedback control of tissue fusion by measurements of tissue surface apparent temperature using, for example, thermal infra-red (IR) sensors may require differing temperature history algorithms to obtain comparable welds depending on the dimensional and structural characteristics of the tissue and on the laser wavelength used. Tissue water is the dominant thermodynamically active material and, for the CO2, Th:YAG, Er:YAG and Ho:YAG lasers, the dominant chromophore, as well.

Histologically, the anastomotic site in arteries collected within 5 minutes of the completion of the anastomosis is marked by the formation of a densely staining, glassy (hyalinized) coagulum of thermally altered adventitial collagen which forms the blood tight bond (figure 29a). However, the full extent of lethal cellular thermal damage cannot be determined until about 3 days later when the heat-killed smooth muscle cells of the medial portion of the arterial wall have undergone necrosis (Figure 29b).

Figure 29:

Longitudinally oriented light microscopic sections of rat femoral artery anastomoses obtained with a CO2 laser beam.90 a) Five minutes after anastomosis. The bond is formed by a glassy coagulum (small arrows) of adventitial collagen that is not carbonized but stains darkly with the eosin dye thus obscuring the fused cut ends (large arrow) of the vessels. The average radius of collagen hyalinization in this section is 0.22 mm. (hematoxylin and eosin -- original magnification 160X). b) Three days after anastomosis. Thermally induced smooth muscle necrosis is detectable at three days (small arrows). The radius of total smooth muscle necrosis (die back) in this specimen is 0.58 mm. The severed ends of the artery are embedded in fibrin, which usually fills the anastomotic site at this time (large arrows), (hematoxylin and eosin stain -- original magnification 80X). This figure is from reference90, used with permission.


Figure 30.

Longitudinally oriented light microscopic sections of dog brachial artery laser anastomosis harvested one week after the bonding procedure90. The boundaries of smooth muscle necrosis (small arrows) are distinct. A fold in the tissue section partially obscures the anastomotic site (large arrow). (Hematoxylin and eosin stain -- original magnification 50X). This figure is from reference90, used with permission.



Argon Laser Heating in Aorta

The effects of high local tissue pressures have been observed in many studies when, for example, highly scattered argon laser irradiation is used to ablate vessel wall. The scattering events result in high fluence rates just below the tissue surface and corresponding thermally induced stresses. An example of water vapor dominated tearing is shown in figure 31. A post mortem sample of human aorta was irradiated with an Argon laser at 514 nm. An ablation crater was created in the surface by explosive vaporization, as well as a sub-surface dissection in a region of high local fluence rate. The scattering resulted in higher fluence rates below the surface and high pressure from confined steam evolution created the dissection approximately along an elastin plate boundary. The laser beam power was 2.5W., beam radius 1mm, duration 2.5s and the maximum surface center temperature recorded by thermographic imaging at video rates was 175 °C immediately prior to explosive vaporization.

Figure 31:

Argon laser irradiation of aorta at 40X original magnification, H&E stain. The largest arrow indicates the approximate center of the 2 mm diameter laser beam. The smallest arrows indicate the approximate edge of the ablation crater (some intimal tissue is missing on the right side of the crater) and the middle sized arrows indicate a sub-surface dissection created by steam generation from high local fluence rates due to scattering processes.



The authors would like to express their appreciation to the Texas Higher Education Coordinating Board (grant number ARP-138), to Valleylab Inc., to the Whitaker Foundation, to the National Science Foundation (grant number BCS-8451123) and to the Air Force Office of Scientific Research through MURI from DDR&E (F49620-98-1-0480) for partial support of this work.



Bloom and D.W. Fawcett, Textbook of Histology, 12th Ed. Chapman and Hall, New York, 1993.Google Scholar


A. P. Spence, Basic Human Anatomy 2nd Ed., Benjamin Cummings, Menlo Park CA, 1986.Google Scholar


L.C. Junqueira, J. Carneiro, A. Contopoulos, Basic Histology, Lance, Los Altos, CA, 1975.Google Scholar


W.D. Gardner and W. A. Osburn, Structure of the Human Body, W.B. Saunders, Philadelphia, 1973.Google Scholar


P.S. Amenta, Elias-Pauly’s Histology and Human Microanatomy, John Wiley & Sons, New York, 1987.Google Scholar


M.D.S. Frame, “Shear rate gradient in arteriolar bifurcations: theoretical and experimental”, Proc. IEEE 21st Ann. Conf. on Engr. in Med. & Biol. Soc, paper no. 2.3.7-2, (proc. on CD-ROM) 1999. https://doi.org/10.1109/IEMBS.1999.802257Google Scholar


M.E. Nimmi (Ed); Collagen; Volume I, Biochemistry, CRC Press, Boca Raton, FL, 1988.Google Scholar


G. Hagen, “Uber die Bewegung des Wassers in engen zylindrischen Rohren”, Pogg. Ann. 46, 423–442, 1839. As described in17.Google Scholar


J. Poiseuille, “Récherches expérimentelles sur le mouvement des liquides dans les tubes de trés petits diamétres”, Compte Rendus, 11, 961–967 and 1041 – 1048 (1840). As described in17.Google Scholar


A.C. Burton, Physiology and Biophysics of the Circulation,Google Scholar


C.G. Caro, T.J. Pedley, R.C. Schroter and W.A. Seed The Mechanics of the Circulation, Oxford University Press, New York, 1978.Google Scholar


Y. Qiu and J.M. Tarbell, “Phase angle between circumferential strain and wall shear stresses affects endothelial cell biochemical production”, Proc. IEEE 21st Ann. Conf. on Engr. in Med. & Biol. Soc., paper no. 2.3.7-3, (proc. on CD-ROM), 1999.Google Scholar


D.R. Karolyi, S. Jin, L. Knoepp, C.M. Brophy, D.P. Giddens, “The effect of aortic coarctation on shear stresses and flow fields through the rabbit aortic bifurcation”, Proc. IEEE 21st Ann. Conf. on Engr. in Med. & Biol. Soc., paper no. 2.3.7-6, (proc. on CD-ROM) 1999. https://doi.org/10.1109/IEMBS.1999.802261Google Scholar


S.L. Soo, Fluid Dynamics of Multiphase Systems, Ginn-Blaisdell, Walthham, MA, 1967.Google Scholar


A. Noordergraaf, Circulatory System Dynamics, Academic Press, New York, 1978.Google Scholar


O. Reynolds, “An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels”, Phil. Trans. Royal Soc., London, 174, 935–982, 1883. As described in17. https://doi.org/10.1098/rstl.1883.0029Google Scholar


H. Schlicting, Boundary Layer Theory, 6th Ed., McGraw-Hill, New York, 1968.Google Scholar


M.L. Rhagavan D.A Vorp, M.P. Federle, M.S. Makaroun, M.W. Webster, “Mechanical stresses on three-dimensionally reconstructed models of abdominal aortic aneurysm”, Proc. IEEE 21st Ann. Conf. on Engr. in Med. & Biol. Soc., paper no. 2.3.7-4, (proc. on CD-ROM) 1999. https://doi.org/10.1109/IEMBS.1999.802259Google Scholar


J. Humphey, “Mechanics of the arterial wall: review and directions”, Critical Reviews in Biomedical Engineering (Begell House Publishing), 23, 1–162 (nos. 1 and 2, complete), 1995.Google Scholar


T.R. Canfield and P.B. Dobrin, “Static elastic properties of blood vessels”, in Handbook of Bioengineering, R. Salak and S. Chien Eds., McGraw-Hill, New York, 1987.Google Scholar


F.H. Silver, D.L. Christiansen and C.M. Buntin, “Mechanical properties of the aorta: a review” CRC Crit. Rev. Biomed. Engr. 17, 323–358, 1989.Google Scholar


Y.C. Fung, S.Q. Liu and J.B. Zhou, “Remodeling of the constitutive relation while a blood vessel remodels itself, Trans. ASME J. Biomech. Engr., 115, 453–459, 1993. https://doi.org/10.1115/1.2895523Google Scholar


J.G. Myers, J.A. Moore, M. Ojha, K.W. Johnston, C.R. Ethier, “Waveform dependence of wall shear stresses in a human right coronary artery”, Proc. IEEE 21st Ann. Conf. on Engr. in Med. & Biol. Soc., paper no. 2.3.7-1, (proc. on CD-ROM) 1999. https://doi.org/10.1109/IEMBS.1999.802256Google Scholar


R.T. Lee and R.D. Kamm, “Vascular mechanics for the cardiologist”, Am. J. Cardiol., 23, 1289–1295, 1994. https://doi.org/10.1016/0735-1097(94)90369-7Google Scholar


J.M. Clark and S. Glagov, “Transmural organization of the arterial media” Arteriosclerosis, 5, 19–34, 1985. https://doi.org/10.1161/01.ATV.5.1.19Google Scholar


S.Q. Liu and Y.C. Fung “Zero stress states of arteries”, Trans. ASME J. Biomech. Engr., 110, 82–84, 1988. https://doi.org/10.1115/1.3108410Google Scholar


Y.C. Fung and S.Q. Liu, “Change of residual strains in arteries due to hypertrophy caused by aortic constriction”, Circ. Res. 65, 1340–1349, 1989. https://doi.org/10.1161/01.RES.65.5.1340Google Scholar


S.Q. Liu and Y.C. Fung “Relationship between hypertension, hypertrophy and opening angle of zero-stress state of arteries following aortic constriction”, Trans. ASME J. Biomech. Engr., Ill, 325–335, 1989. https://doi.org/10.1115/1.3168386Google Scholar


H.C. Han and Y.C. Fung “Species dependence of the zero-stress state of the aorta: pig versus rat”, Trans. ASME J. Biomech. Engr., 113, 446–451, 1991. https://doi.org/10.1115/1.2895425Google Scholar


J. Vossoughi, Z. Hedjazi and F.S. Bonis, “Intimal residual stress and strain in large arteries”, In 1993 Advances in Bioengineering, New York, 1993.Google Scholar


T.J. Van Dyke and A. Hoger, “Growth-induced residual stresses and opening angles in soft biological tissues”, Proc. 36th Ann. Meeting, Soc. of Engr. Sciences, paper no. MA7-2, 1999.Google Scholar


R.H. Cox, “Comparison of carotid artery mechanics in the rat, rabbit and dog”, Am. J. Physiol., 234, H280–288, 1978.Google Scholar


R.H. Cox, “Regional variation of series elasticity in canine arterial smooth muscle”, Am. J. Physiol., 234, H542–551, 1978. https://doi.org/10.1152/ajpcell.1978.234.5.C146Google Scholar


R.H. Cox, “Passive mechanics and connective tissue composition of canine arteries”, Am. J. Physiol., 234, H533–541, 1978.Google Scholar


R.H. Cox, “Differences in mechanics of arterial smooth muscle from hind limb arteries”, Am. J. Physiol., 234, H649–656, 1978.Google Scholar


F.P. Incropera and D.P. DeWitt, Fundamentals of Heat and Mass Transfer, 4th Ed., John Wiley & Sons, New York, 1996.Google Scholar


W.M. Rosenhow and H.Y. Choi, Heat, Mass and Momentum Transfer, Prentice-Hall, Englewood Cliffs, NJ, 1961.Google Scholar


F. Kreith and M.S. Bohn, Principles of Heat Transfer, West Publ., St. Paul, MN, 1993.Google Scholar


W.M. Kays and M.E. Crawford, Convective Heat and Mass Transfer, McGraw-Hill, New York, 1980.Google Scholar


A. Shitzer and R.C. Eberhart, Heat Transfer in Medicine and Biology, Plenum Press, New York, 1985. https://doi.org/10.1007/978-1-4684-8285-0Google Scholar


R. Siegel and J.R. Howell, Thermal Radiation Heat Transfer, 3rd Ed., Taylor & Francis, Washington, DC, 1992.Google Scholar


J.C. Chato, “Heat transfer to blood vessels”, Trans, of the ASME, J. of Biomech. Engr., 102, 110–118, 1980. https://doi.org/10.1115/1.3138205Google Scholar


M.M. Chen and K.R. Holmes, “Microvasculature contributions in tissue heat transfer, In Thermal Characteristics of Tumors: Applications in Detection and Treatment, R.K. Jain and P.M. Gullino Eds. Ann. N.Y. Acad, of Sci., 335, 137–150, 1980.Google Scholar


M.M. Chen, “The tissue energy balance equation”, Chapter 7 in Heat Transfer in Medicine and Biology, Shitzer and Eberhart, Eds., Plenum Press, New York, 1985.Google Scholar


H.H. Pennes, “Analysis of tissue and arterial blood temperatures in the resting human forearm”, J. Appl. Physiol. 1, 93–122, 1948. https://doi.org/10.1152/jappl.1948.1.2.93Google Scholar


E.H. Wissler, “A mathematical model of the human thermal system”, Bull. Math, and Biophys., 26, 147–166, 1964. https://doi.org/10.1007/BF02476835Google Scholar


J.A. Stolwijk, “Mathematical model of thermoregulation”, In Physiological and Behavioral Temperature Regulation, J.D. Hardy et al Eds., pp703–721, Thomas, Springfield IL, 1971.Google Scholar


W.J. Song, S. Weinbaum and L.M. Jiji, “A theoretical model for peripheral tissue heat transfer using the bioheat equation of Weinbaum and Jiji”, Trans, of ASME, J. of Biomech. Engr., 109, 72–78, 1987. https://doi.org/10.1115/1.3138646Google Scholar


C.K. Charny, S. Weinbaum and R.L. Levin, “An evaluation of the Weinbaum-Jiji bioheat equatin for normal and hyperthermic conditions”, Trans, of ASME, J. of Biomech. Engr., 112(1), 80–87, 1990. https://doi.org/10.1115/1.2891130Google Scholar


A.R. Moritz and F.C. Henriques, “Studies in thermal injury: II the relative importance of time and surface temperature in the causation of cutaneous burns” Am. J. Pathol. 23(5), 695–720, 1947.Google Scholar


J.A. Pearce, S. Thomsen, H. Vijverberg and T.J. McMurray, “Kinetics for birefringence changes in thermally coagulated rat skin collagen”, Proc. SPIE 1876, 180–186, 1993.Google Scholar


S. Thomsen, J.A. Pearce and W. F. Cheong, “Loss of birefringence as a low temperature marker of thermal damage in tissues”, IEEE Transactions on Biomedical Engineering BME-36 (12), 1174–1179, 1989. https://doi.org/10.1109/10.42111Google Scholar


K.R. Diller, J.W. Valvano and J.A. Pearce, Bioheat Transfer, Section 4 of Chapter 4 “Applications” in Handbook of Thermal Engineering, F. Kreith, Ed., CRC Press (in press).Google Scholar


A.N. Takata et al, “Thermal model of laser induced eye damage”, Final Report USAF School of Aerospace Med., Brooks AFB, TX, Contract F41609-74-C-0005, IIT Research Institute, Chicago, IL, 1974.Google Scholar


A.J. Welch, G.D. Polhamus, “Measurement and prediction of thermal injury in the retina of the Rhesus monkey”, IEEE Trans. Biomed. Engr. BME-31, 633–644, 1984. https://doi.org/10.1109/TBME.1984.325313Google Scholar


R. Birngruber, F. Hillenkamp, V-P. Gabel, “Theoretical investigations of laser thermal retinal injury”, Health Phys., 48, 781–796, 1985. https://doi.org/10.1097/00004032-198506000-00006Google Scholar


F.C. Henriques, “Studies of thermal injury V: The predictability and significance of thermally induced rate processes leading to irreversible injury”, Arch. Pathol. 43, 489–502, 1947.Google Scholar


J.A. Weaver and A.M. Stoll, “Mathematical model of skin exposed to thermal radiation”, Aerospace Med., 40 (1), 24, 1969.Google Scholar


Y. Yang, A.J. Welch and H.G. Rylander III, “Rate process parameters of albumen”, Lasers in Surg. & Med., 11, 188–190, 1991. https://doi.org/10.1002/(ISSN)1096-9101Google Scholar


R. Agah, “Quantitative characterization of arterial tissue thermal damage”, MSE thesis, The University of Texas at Austin, 1988.Google Scholar


J.A. Pearce, S. Thomsen, H. Vijverberg and T. J. McMurray, “Kinetics for birefringence changes in thermally coagulated rat skin collagen”, Proc. SPIE, 1876, 180–186, 1993.Google Scholar


T. Kang, J. Resar, and J.D. Humphrey, “Heat-induced changes in the mechanical behavior of passive coronary arteries”, Trans. ASME J. Biomech. Engr., 117, 86–93, 1995. https://doi.org/10.1115/1.2792274Google Scholar


J.H Keenan and F.G Keyes, Thermodynamic Properties of Steam 1st Ed., John Wiley, New York, 1948.Google Scholar


W. Brutasaert, Evaporation Into the Atmosphere: Theory, History and Applications, D. Reidel, Dordrecht Netherlands, 1982. https://doi.org/10.1007/978-94-017-1497-6Google Scholar


J.H. Torres, M. Motamedi, J.A. Pearce and A.J. Welch, “Experimental evaluation of mathematical models for predicting the thermal response to laser irradiation”, Appl. Optics, 32(4), 597–606, 1993. https://doi.org/10.1364/AO.32.000597Google Scholar


V.P. Skripov, Metastable Liquids, (translated from Russian by R. Kondor), John Wiley and Sons, New York, 1974.Google Scholar


E. Gerum, J. Straub and U. Grigull, “Examination of the metastable state of liquids”, Chapter 19 in Heat Transfer in Boiling, E. Hahne and U. Grigull Eds., Academic Press, New York, 1977.Google Scholar


D. T. Paris and F.K. Hurd Basic Electromagnetic Theory McGraw-Hill, New York, 1969.Google Scholar


G. Roussy and J.A. Pearce Foundations and Industrial Applications of Microwaves and Radio Frequency Fields: Physical and Chemical Processes; John Wiley & Sons, Chichester, UK, 1995.Google Scholar


J.A. Pearce and S. Thomsen, “The effect of vessel architecture on fusion by radio frequency current, Proc. SPIE, 3249, 217–228, 1998.Google Scholar


H.P. Schwann, “Biophysics of Diathermy”, Chapter 3 of Therapeutic Heat and Cold, S. Licht Ed.; Waverly Press, Baltimore MD, 1965.Google Scholar


L.A. Geddes and L.E. Baker, “The specific resistance of biological material. A compendium of data for the biomedical engineer and physiologist”, Med. & Biol. Engr., 5, 271–293, 1967. https://doi.org/10.1007/BF02474537Google Scholar


A.J. Welch and M.J.C. vanGemert, (Eds.), Optical-Thermal Response of Laser-Irradiated Tissue, Plenum, New York, 1995. https://doi.org/10.1007/978-1-4757-6092-7Google Scholar


M.L. Wolbarsht, Laser Applications in Medicine and Biology, Plenum, New York, 1974. https://doi.org/10.1007/978-1-4615-7323-4Google Scholar


E.K. Chan, B. Sorg, D. Protsenko, M. O'Neil, M. Motamedi, and A.J. Welch, “Effects of compression on soft tissue optical properties”, IEEE J. Sel. Topics in Quant. Mech., 2(4), 943–950, 1996. https://doi.org/10.1109/2944.577320Google Scholar


H.F. Bowman, E.G. Cravalho and M. Woods, “Theory, measurement and application of thermal properties of biological materials”, Ann. Rev. of Biophys & Bioengr. 4, pp43–80, 1975. https://doi.org/10.1146/annurev.bb.04.060175.000355Google Scholar


J.W. Valvano and C. Chitsabesan, “Thermal conductivity and diffusivity of arterial wall and atherosclerotic plaque”, Lasers in the Life Sciences, 1, 219–229, 1987.Google Scholar


V.K. Dalvi, “Thermal conductivity measurement of arterial walls”, M.S. Thesis, Biomedical Enginering Program, The University of Texas at Austin, 1993.Google Scholar


T.A. Balasubramaniam and H.F. Bowman, “Thermal conductivity and thermal diffusivity of biomaterials. A simultaneous measurement technique”, J. Biomech. Engr., 99, 148–154, 1977. https://doi.org/10.1115/1.3426282Google Scholar


K.E. Spells, “The thermal conductivities of some biological fluids”, Phys. Med. & Biol., 5, 139–153, 1960. https://doi.org/10.1088/0031-9155/5/2/304Google Scholar


T.E. Cooper and G.J. Trezek, “Correlation of thermal properties of some human tissues with water content”, Aerosp. Med., 42, 24–27, 1971.Google Scholar


T.A. Magnusen and J.A. Pearce, “Temperature controlled electrosurgical vessel sealing”, Proc. IEEE-EMBS 16th Annual Meeting (CD-ROM), 1994. https://doi.org/10.1109/IEMBS.1994.415219Google Scholar


J.S. Kennedy J.A. Pearce and S. Thomsen, “Mechanisms of electrosurgical fusion for large vessel hemostasis, Proc. 7th Intl. Mtng, Soc. Min. Invas. Therapy (in J. Min. Invas. Therap., 4 Suppl 1, p19), 1995.Google Scholar


G.M. Lemole, R.R. Anderson and S. DeCoste; “Preliminary evaluation of collagen as a component in the thermally-induced ‘weld’” Proc. SPIE 1422:116–122, 1991.Google Scholar


R.A. White, G.E. Kopchok and G. White, “Vascular tissue bonding using laser surgery” Lasers in Cardiovascular Medicine and Surgery: Fundamentals and Techniques. G. Abela ed. pp361–372, Kluwer Academic Press. Boston, 1990.Google Scholar


G.E. Kopchock, R.A. White, H.J. White et al; “CO2 and argon laser vascular welding: acute histologic an thermodynamic comparison”, Lasers in Surgery and Medicine 8(6) 584–588, 1988. https://doi.org/10.1002/(ISSN)1096-9101Google Scholar


R. Schober, F. Ulrich, T. Sander, H. Durselen and S. Hessel; “Laser-induced alteration of collagen substructure allows microsurgical tissue welding” Science 232:1411–1421, 1986. https://doi.org/10.1126/science.3715454Google Scholar


M.R. Treat, C.O. Mehmet, and L.S. Bass, “New technologies and future applications of surgical lasers” Surgical Clinics of North America 72(3), 705–42, 1992. https://doi.org/10.1016/S0039-6109(16)45742-7Google Scholar


C.R. Neblett, J.R. Morris and S. Thomsen; “Laser-assisted microsurgical anastomosis” Neurosurgery 19, 914–934, 1986. https://doi.org/10.1227/00006123-198612000-00006Google Scholar


J.A. Pearce and S. Thomsen, “Kinetic models of tissue fusion processes”, Proc. SPIE 1643:251–260, 1992.Google Scholar


R. Schrober, F. Ulrich, T. Sander, H. Durselen, and S. Hessel, “Laser-induced alteration of collagen substructure allows microsurgical tissue welding”, Science 232, 1411–1421, 1986.Google Scholar


M. Quigley, H. Bailes, H. Kwan, L. Cerullo, J.T. Brown; “Aneurysm formation after low power CO2 laser assisted vascular anastomosis” Neurosurgery 18, 292–298, 1986. https://doi.org/10.1227/00006123-198603000-00006Google Scholar

© (2000) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
John A. Pearce, John A. Pearce, Sharon L. Thomsen, Sharon L. Thomsen, } "Blood vessel architectural features and their effects on thermal phenomena", Proc. SPIE 10297, Matching the Energy Source to the Clinical Need: A Critical Review, 1029706 (24 January 2000); doi: 10.1117/12.375226; https://doi.org/10.1117/12.375226

Back to Top