Modelling ultrasound-induced mild hyperthermia of hyperplasia in vascular grafts
© Brinton et al; licensee BioMed Central Ltd. 2011
Received: 9 February 2011
Accepted: 3 November 2011
Published: 3 November 2011
Expanded polytetrafluoroethylene (ePTFE) vascular grafts frequently develop occlusive neointimal hyperplasia as a result of myofibroblast over-growth, leading to graft failure. ePTFE exhibits higher ultrasound attenuation than native soft tissues. We modelled the selective absorption of ultrasound by ePTFE, and explored the feasibility of preventing hyperplasia in ePTFE grafts by ultrasound heating. Specifically, we simulated the temperature profiles of implanted grafts and nearby soft tissues and blood under ultrasound exposure. The goal was to determine whether ultrasound exposure of an ePTFE graft can generate temperatures sufficient to prevent cell growth on the graft without damaging nearby soft tissues and blood.
Ultrasound beams from two transducers (1.5 and 3.2 MHz) were simulated in two graft/tissue models, with and without an intra-graft cellular layer mimicking hyperplasia, using the finite-difference time-domain (FDTD) method. The resulting power deposition patterns were used as a heat source for the Pennes bioheat equation in a COMSOL® Multiphysics heat transfer model. 50°C is known to cause cell death and therefore the transducer powers were adjusted to produce a 13°C temperature rise from 37°C in the ePTFE.
Simulations showed that both the frequency of the transducers and the presence of hyperplasia significantly affect the power deposition patterns and subsequent temperature profiles on the grafts and nearby tissues. While neither transducer significantly raised the temperature of the blood, the 1.5-MHz transducer was less focused and heated larger volumes of the graft and nearby soft tissues than the 3.2-MHz transducer. The presence of hyperplasia had little effect on the blood's temperature, but further increased the temperature of the graft and nearby soft tissues in response to either transducer. Skin cooling and blood flow play a significant role in preventing overheating of the native tissues.
Modelling shows that ultrasound can selectively heat ePTFE grafts and produce temperatures that cause cell death on the graft. The temperature increase in blood is negligible and that in the adjacent soft tissues may be minimized by skin cooling and using appropriate transducers. Therefore, ultrasound heating may have the potential to reduce neointimal hyperplasia and failure of ePTFE vascular grafts.
There are two primary types of permanent hemodialysis accesses, native arteriovenous (AV) fistulae and expanded-polytetrafluoroethylene (ePTFE) AV grafts. When compared to the fistula, AV grafts provide vascular access with higher blood flow rates, do not require a maturation period, and can be used sooner after surgical placement. Grafts are a particularly useful form of vascular access in patients whose vasculature is insufficient to support a functional native AV fistula. The latter is important because up to 60% of AV fistulae in the United States fail to mature and become functional. Despite these advantages, only 23.3% of the permanent hemodialysis vascular accesses in the United States are AV grafts , largely due to their low one- and two-year primary patency rates, being 50% and 25%, respectively [2, 3].
The low primary patency rates of ePTFE grafts are attributed to their propensity to develop excessive cell growth in the graft lumen, preferentially near the graft-venous anastomosis site . This growth, termed neointimal hyperplasia (NH), causes lumen stenosis that predisposes to thrombosis and occlusion. If the occurrence of NH could be prevented or reduced, the synthetic graft should be an excellent form of dialysis vascular access for the large population of dialysis patients whose vasculature cannot support a functional AV fistula.
Strategies to prevent or reduce NH in synthetic AV grafts include the use of drug-eluting stents [5, 6], endothelialization of the graft lumen [7–9], antiproliferative drugs delivered using hydrogel systems to prevent cell overgrowth [5, 6, 10–12], and catheter-based gamma radiation . The stent and catheter treatments require invasive procedures, and none of the strategies has proven effective in humans.
Ultrasound has been used for a variety of therapeutic applications including muscle diathermy and bone healing , drug delivery  and mild hyperthermia [16–18]. Ultrasound exposure has been shown to reduce smooth muscle cell proliferation in vitro . Intravascular ultrasound has also been shown to reduce NH in implanted stents by 35% . Ultrasound is comprised of acoustic or vibrational waves that carry mechanical energy. The amount of acoustic energy that is deposited in a material and converted into heat is approximately proportional to the material's viscosity . Since ePTFE exhibits 5 to 10 times more ultrasound attenuation (absorption and scattering) than native soft tissues, transcutaneous delivery of ultrasound may provide a noninvasive method to selectively heat an implanted ePTFE graft, thereby preventing cell growth and NH on the graft surface [22, 23].
In this study, we modelled the selective absorption of ultrasound in ePTFE and the corresponding thermal profiles to determine whether ultrasound exposure can elevate the graft to temperatures capable of inflicting NH cell death without significantly damaging the nearby tissues and blood. Specifically, we used mathematical and computational modelling to simulate the ultrasound power absorbed and the effectual temperature profiles in implanted grafts and nearby tissues and blood under ultrasound generated from two different transducers.
To obtain a 3D pattern of the ultrasound power density deposited in the graft/tissue models, we applied the FDTD numerical technique to equations that describe the propagation of an acoustic pressure wave through a medium. The 3D models have a form composed of individual regularly spaced cuboidal voxels, each with possibly different densities, speeds of sound and attenuation properties, and with planar front faces. Acoustic waves can be described using Newton's force equation and the conservation of mass. Using finite differences these equations can be solved for the acoustic pressure and the particle velocity in the x, y and z directions. The method of finite differences involves replacing classical derivatives with discrete differences in time and on a 3D grid space. For example, the derivative of the x-direction particle velocity with respect to x becomes the difference of the particle velocity at locations x and x + Δx divided by Δx, the FDTD grid resolution in the x direction (see details in Appendix A). Using initial pressure values, the particle velocities in the x, y and z directions for the entire grid and the first time step are calculated. Then, using the newly calculated particle velocities, new pressure values for the same time step are calculated. Continuing in this manner the program, written in MATLAB, works by iterating back and forth between space and time steps solving for the particle velocities and pressure as the acoustic wave propagates through the 3D gridded models of the inhomogeneous medium.
The accuracy of the FDTD method depends on the chosen temporal and spatial step size. While smaller step sizes more closely approximate the original differential equations and tissue geometry, small step sizes result in large models and long calculation times. FDTD simulations were completed on a computer using an Intel® Core™ i7 processor at 2.67 GHz and with 11.9 GB RAM. To avoid aliasing, the model grid spacings Δx, Δy and Δz were set to 0.1 mm for the 1.5-MHz models and to 0.05 mm for the 3.2-MHz models, roughly 1/10th of the acoustic wavelength in water for the respective frequencies. The overall dimensions of the 1.5-MHz models were 2 × 2 × 4 cm, producing a model with 200 × 200 × 400 elements. Because the 0.05-mm step size in the 3.2-MHz models would produce eight times as many elements as the 1.5-MHz models (exceeding the memory limits of the computer), the overall dimensions of the 3.2-MHz model were reduced to 1.5 × 1.5 × 3 cm, resulting in 300 × 300 × 600 elements. Some clipping of the incident power occurred since the areas of the front planes of the models (1.5 × 1.5 cm and 2 × 2 cm) were slightly smaller than the extent of the beams from the transducer. This clipping caused 2.5% and 7.5% reductions in the power from the 1.5- and 3.2-MHz transducers, respectively.
Physiological models for acoustic simulation
Material properties used for acoustic FDTD and COMSOL thermal modelling.
Density ( ρ ) [kg/m 3 ]
Specific heat capacity ( Cp ) [J/kg·K]
Thermal conductivity ( k ) [W/m·K]
Perfusion ( ω ) [s -1 ]
Speed of sound ( c ) [m/s]
Attenuation (α) at 1.5/3.2 MHz [Np/m]
where α is material pressure attenuation coefficient, ρ is material mass density, c is speed of sound, p is the rms pressure and f is the frequency of the wave. During the duration of the ultrasound exposure, the magnitude of the applied acoustic power and pressure was held constant. After calculating Q from the 3D FDTD pressure solution, it was used as a heat source in the Pennes bioheat equation.
Heat transfer using COMSOL Multiphysics
The general 3D heat transfer application in COMSOL Multiphysics, a commercial finite-element method (FEM) modelling software designed to solve heat transfer, acoustic, fluid flow and other engineering problems described by differential equations, was used. The simplified graft/tissue models used for the FDTD calculations were also used in COMSOL. The thermal properties assumed for the blood, ePTFE, muscle and fat regions are listed in Table 2.
To simulate skin cooling, the front fat boundary was held constant at 20°C. Except where the blood entered the graft/tissue model (where the entering blood temperature was held constant at body temperature of 37°C), all other exterior boundary conditions were considered insulated, allowing no heat flux across exterior boundaries. Insulated boundaries provide an accurate assumption if the simulated volume is large enough that the side and rear boundaries are much farther from the ultrasound heated region than the nominal diffusion length.
Physiologically relevant blood flow volumes and velocities were included in the modelling, and the COMSOL Multiphysics software has the capability to model flow within the heat transfer application. In implanted grafts, blood flow near the anastomosis is often turbulent, and may contribute to neointimal proliferation . However, for simplicity, convection in the graft lumen was modelled using laminar blood flow. Volumetric flow rates based on clinical data were used to determine the peak and mean velocities for the laminar velocity distribution. For the graft without NH, given a luminal cross-sectional area of 28.3 mm2, an average flow rate of 412 ml/min was assumed , resulting in an average flow velocity of 0.243 m/s. Since grafts exhibiting flow rates less than 100 ml/min are considered failing, for the 1.2-mm thick NH model this flow rate was used to calculate an average flow velocity of 0.164 m/s. From the average flow velocities the laminar flow profiles were defined and implemented in the COMSOL models.
where ω is the blood perfusion rate of the tissue [s-1] and κ is a term to account for the incomplete thermal transfer to perfused blood. In our simulations κ was assumed to be zero (also assumed by Pennes and most subsequent investigators ) to model complete thermal equilibration between the blood and tissue.
Thermal damage to healthy native tissues, a major concern during hyperthermia, can be reduced with active skin surface cooling [35, 36]. The pre-exposure skin cooling was simulated in COMSOL by solving a steady-state model prior to ultrasound exposure. In this model, the skin boundary was held constant at 20°C and the blood entering the graft and perfusing the tissue was held constant at 37°C. Using the skin-cooling model, the temperatures resulting from 10-second or 30-second transient applications with the ultrasound power deposition values from the FDTD simulations were calculated. To explore the effects of intermittent ultrasound exposure, the temperatures during the cool-down period following the ultrasound exposure were simulated.
Simulated pressure and power depositions for the 1.5- and 3.2-MHz transducers
Simulated temperature profiles
1.5- and 3.2-MHz transducer modelling
Wave frequency is an important parameter in therapeutic ultrasound. Ultrasound attenuation increases with frequency; consequently, high-frequency ultrasound is used for imaging and therapy near or at the body surface, while low-frequency ultrasound is appropriate for deeper imaging and therapy [37, 38]. Frequency also has an effect on the dimensions of the focal zone. Higher frequency waves focus to smaller zones . The increase in frequency from 1.5 MHz to 3.2 MHz, a factor of 2.1, causes the focal length and width to decrease by this same factor for a given transducer geometry. In addition, the difference in the ratio of transducer focal length to diameter for the 1.5- and 3.2-MHz transducers, 1.8 and 1.4, respectively, further reduces the length of the focal zone by a factor of 1.3. These differences are evident in Figure 2. The focal zone for the 3.2-MHz transducer was small enough that power was deposited primarily in the rear graft wall (Figure 3b).
The amount of ultrasound power absorbed by a material depends on, among other things, the material's acoustic attenuation (Equation 1). The attenuation of ePTFE is almost 10 times greater than the tissues surrounding the ePTFE graft (Table 2). Consequently, the ePTFE material absorbed much more power than native tissue from the ultrasound beam. This strategic concentration of power created the high temperatures within the graft and the mild heating in the surrounding tissues (Figures 4 and 5).
The higher frequency and the smaller focal zone (resulting in higher beam intensity) of the 3.2-MHz transducer caused a higher, more concentrated power deposition pattern. Although the 3.2-MHz transducer operated at 63% of the power of the 1.5-MHz transducer it deposited a maximum Q approximately three times more than the 1.5-MHz transducer (Figures 3a and 3b). Interestingly, the 3.2-MHz transducer generated similar temperatures while depositing a higher maximum Q with less overall power than the 1.5-MHz transducer. The roughly three-fold Q of the 3.2-MHz transducer induced a faster but more localized increase in temperature than the 1.5-MHz transducer. As a result, the more localized heating produced larger spatial temperature gradients. These large gradients gave rise to greater heat transfer away from the focal region of the 3.2-MHz transducer, limiting graft and tissue temperatures to be similar to those in the 1.5-MHz transducer model after 30 seconds of exposure.
Understanding cellular response to specific temperatures provides relevance for the thermal modelling. While cellular response to thermal exposure varies with cell species and origin , both apoptosis (organized cell death) and necrosis (ruptured cellular membrane) can occur in response to thermal exposure. In general, temperatures below 43°C can cause a heat-shock-protein response but only induce apoptosis mildly [40, 41]. Samali et al.  and Harmon et al.  exposed suspended cells for one hour and for 30 minutes, respectively, and concluded the turning point from apoptosis to necrosis occurred at 45°C. Apoptosis, rather than necrosis, is the preferred cell death pathway for our target clinical application. Necrotic death causes inflammation and the attraction of macrophages that may exacerbate NH formation. Hehrlein et al. reported that vascular damage from high-temperature exposure led to increased intimal growth in canines . In addition, temperatures exceeding 47°C cause epidermal and dermal injury . Taken together, the target graft temperature should be between 45-47°C to induce apoptotic death while limiting necrotic death.
Differences in thermal exposure between the 1.5- and 3.2-MHz transducer simulations may influence the extent of cell death on implanted ePTFE grafts. For example, the 3.2-MHz transducer heated the graft to a higher temperature after 5 seconds (Figures 5b and 5d). Consequently, using the 3.2-MHz transducer the luminal graft surface may reach effective elevated temperatures more quickly, thereby heating hyperplasia at effective temperatures for a greater portion of the 30 second ultrasound pulse. The 3.2-MHz transducer modelling also showed less native tissue heating above 43°C. As a result, the 3.2-MHz transducer's more confined spatial heating may facilitate thermal exposure at the anastomosis site with less damage to adjacent tissues, which is a major advantage when compared to the 1.5-MHz transducer.
Our modelling results suggest that the ultrasound should be delivered in an intermittent manner in our target clinical application. Non-continuous ultrasound therapy helps to reduce unwanted heating of healthy tissues by allowing cooling between the ultrasound pulses. Pulsed ultrasound exposure could also help ensure that hyperplasia receives only mild heat, thereby avoiding tissue necrosis. To avoid residual temperature buildup during pulsed ultrasound exposure, a 30-second exposure should be followed by a 30-second rest, for a duty cycle of 0.5. Shorter rest times would not allow sufficient time for the tissue to return to 37°C before the next ultrasound pulse exposed the ePTFE graft. Consequently, a treatment equivalent to 30 minutes at 45°C would take twice as long, lasting up to an hour.
During the ultrasound exposure the low thermal conductivity of the tissue slowed thermal diffusion into the tissue. Thus, after 30 seconds of exposure, the rise in temperature was confined to and in close proximity of the ePTFE material (Figure 5). Furthermore, the heat that did transfer into the tissue also dissipated slowly during the cooling cycle, evident by the regions in the fat and muscle which cooled to the pre-exposure temperatures last (Figure 6). Consequently, in our target clinical application the cooling time between consecutive ultrasound pulses, which limited by the thermal conductivity of the tissue, is necessary in order to avoid thermal buildup.
Due to convection, the blood-ePTFE boundary (Figures 5a and 5c) and the blood-NH boundary (Figures 5b and 5d) did not reach temperatures capable of inducing apoptotic cell death, suggesting ultrasound heating may not prevent NH formation in NH-free grafts and may not be able to treat pre-existent NH entirely, respectively. Attempts to heat more of the NH region using longer ultrasound exposures could cause thermal damage to tissue adjacent to the graft where the heat dissipated more slowly. However, temperatures capable of inducing apoptosis were generated near the NH-ePTFE boundary (Figures 5b and 5d). As a result, repeated ultrasound exposure may reduce NH growth and delay the development of significant stenosis. In the presence of NH, a graft can still remain functional as long as the NH size is not so big that it impedes blood flow. Therefore, even though the proposed approach may not completely prevent the formation of NH or remove a pre-existing NH entirely, repeated ultrasound exposure may restrict the NH size, thereby lengthening the patency rate of the graft.
Results of our modelling showed selective ePTFE graft heating could be achieved by either the 1.5-MHz or the 3.2-MHz transducer. Both transducers can lead to temperatures capable of inducing apoptosis near the NH-ePTFE boundary. The temperature increase in blood was negligible using either transducer, while the temperature increase in the adjacent soft tissue was smaller when the 3.2-MHz transducer was used. Tissue temperature variations due to the presence or absence of NH suggest the need of patient-specific ultrasound exposure protocols. Overall, intermittent ultrasound heating may have the potential for clinical application to reduce neointimal hyperplasia and failure of ePTFE vascular grafts.
Detailed FDTD method
ū = u x ẋ + u y ẏ + u z ż is the particle velocity,
u x , u y , and u z are the respective particle velocities in the ẋ, ẏ, and ż directions,
η is the dynamic coefficient of shear viscosity of the material,
η' is the dynamic coefficient of bulk viscosity of the material,
ρ o is the average mass density of the material, and
Κ is the compressibility of the material.
The FDTD algorithm works iteratively between (A3) and (A4), as well as between the respective equations for the u y and u z velocity solutions not shown here, to calculate the pressure and velocities at each time step throughout the entire model grid space.
AKC received support from the National Institutes of Health (RO1HL67646) and the Veterans Affairs Merit Review Program. DAC received support from the University of Utah Research Foundation. YES and MRB received support from the American Heart Association (0970307N and 0765113Y).
- 2009 Annual Report of the Dialysis Outcomes and Practice Patterns Study: Hemodialysis Data 1999-2008. Arbor Research Collaborative for Health, Ann Arbor, MI
- Pisoni R, Young E, Dykstra D, Greenwood R, Hecking E, Gillespie B, Wolfe R, Goodkin D, Held P: Vascular access use in Europe and United States: Results from the DOPPS. Kidney Int. 2002, 61: 305-316. 10.1046/j.1523-1755.2002.00117.x.View ArticlePubMed
- Pisoni R, Young E, Mapes D, Keen M, Port F: Vascular access use and outcomes in the U.S., Europe, and Japan: results from the Dialysis Outcomes and Practice Patterns Study. Nephrol News Issues. 2003, 17: 38-43,47.PubMed
- Li L, Terry CM, Shiu YT, Cheung AK: Neointimal hyperplasia associated with synthetic hemodialysis grafts. Kidney Int. 2008, 74: 1247-1261. 10.1038/ki.2008.318.PubMed CentralView ArticlePubMed
- Bordenave L, Fernandez P, Remy-Zolghadri M, Villars S, Daculsi R, Midy D: In vitro endothelialized ePTFE prostheses: clinical update 20 years after the first realization. Clin Hemorheol Microcirc. 2005, 33: 227-234.PubMed
- Pawlowski K, Rittgers S, Schmidt S, Bowlin G: Endothelial cell seeding of polymeric vascular grafts. Front Biosci. 2004, 9: 1412-1421. 10.2741/1302.View ArticlePubMed
- Rotmans J, Heyligers J, Verhagen H, Velema E, Nagtegaal M, de Kleijn D, de Groot F, Stroes E, Pasterkamp G: In vivo cell seeding with anti-CD34 antibodies successfully accelerates endothelialization but stimulates intimal hyperplasia in porcine arteriovenous expanded polytetrafluoroethylene grafts. Circulation. 2005, 112: 12-18. 10.1161/CIRCULATIONAHA.104.504407.View ArticlePubMed
- Colombo A, Drzewiecki J, Banning A, Grube E, Hauptmann K, Silber S, Dudek D, Fort S, Schiele F, Zmudka K, Guagliumi G, Russell M: Randomized study to assess the effectiveness of slow- and moderate-release polymer-based paclitaxel-eluting stents for coronary artery lesions. Circulation. 2003, 108: 788-794. 10.1161/01.CIR.0000086926.62288.A6.View ArticlePubMed
- Regar E, Serruys P, Bode C, Holubarsch C, Guermonprez J, Wijns W, Bartorelli A, Constantini C, Degertekin M, Tanabe K, Disco C, Wuelfert E, Morice M: Angiographic findings of the multicenter Randomized Study With the Sirolimus-Eluting Bx Velocity Balloon-Expandable Stent (RAVEL): sirolimus-eluting stents inhibit restenosis irrespective of the vessel size. Circulation. 2002, 106: 1949-1956. 10.1161/01.CIR.0000034045.36219.12.View ArticlePubMed
- Kopecek J: Smart and genetically engineered biomaterials and drug delivery systems. Eur J Pharm Sci. 2003, 20: 1-16. 10.1016/S0928-0987(03)00164-7.View ArticlePubMed
- Masaki T, Rathi R, Zentner G, Leypoldt JK, Mohammad SF, Burns GL, Li L, Zhuplatov S, Chirananthavat T, Kim SJ, Kern S, Holman J, Kim SW, Cheung AK: Inhibition of neointimal hyperplasia in vascular grafts by sustained perivascular delivery of paclitaxel. Kidney Int. 2004, 66: 2061-2069. 10.1111/j.1523-1755.2004.00985.x.View ArticlePubMed
- Kohler TR, Toleikis PM, Gravett DM, Avelar RL: Inhibition of neointimal hyperplasia in a sheep model of dialysis access failure with the bioabsorbable Vascular Wrap paclitaxel-eluting mesh. J Vasc Surg. 2007, 45: 1029-1037. 10.1016/j.jvs.2007.01.057. discussion:1037-1028View ArticlePubMed
- Sun S, Beitler J, Ohki T, Calderon T, Schechner R, Yaparpalvi R, Berman J, Tellis V, Greenstein S: Inhibitory effect of brachytherapy on intimal hyperplasia in arteriovenous fistula. J Surg Res. 2003, 115: 200-208. 10.1016/S0022-4804(03)00253-1.View ArticlePubMed
- Hadjiargyrou M, McLeod K, Ryaby JP, Rubin C: Enhancement of fracture healing by low intensity ultrasound. Clin Orthop Relat Res. 1998, 355: S216-229.View ArticlePubMed
- Gao Z, Kennedy AM, Christensen DA, Rapoport NY: Drug-loaded nano/microbubbles for combining ultrasonography and targeted chemotherapy. Ultrasonics. 2008, 48: 260-270. 10.1016/j.ultras.2007.11.002.PubMed CentralView ArticlePubMed
- Hynynen K, Davis KL: Small cylindrical ultrasound sources for induction of hyperthermia via body cavities or interstitial implants. Int Journal Hyp. 1993, 9: 263-274. 10.3109/02656739309022539.View Article
- Hurwitz MD, Kaplan ID, Svensson GK, Hansen MS, Hynynen K: Feasibility and patient tolerance of a novel transrectal ultrasound hyperthermia system for treatment of prostate cancer. Int Journal Hyp. 2001, 17: 31-37. 2001View Article
- Hynynen K, Deyoung D: Temperature elevation at muscle-bone interface during scanned, focused ultrasound hyperthermia. Int Journal Hyp. 1988, 4: 267-279. 10.3109/02656738809051103.View Article
- Lawrie A, Brisken A, Francis A, Tayler D, Chamberlain J, Crossman D, Cumberland D, Newman C: Ultrasound enhances reporter gene expression after transfection of vascular cells in vitro. Circulation. 1999, 99: 2617-2620.View ArticlePubMed
- Fitzgerald PJ, Takagi A, Moore MP, Hayase M, Kolodgie FD, Corl D, Nassi M, Virmani R, Yock PG: Intravascular sonotherapy decreases neointimal hyperplasia after stent implantation in swine. Circulation. 2001, 103: 1828-1831.View ArticlePubMed
- Christensen D: Ultrasonic Bioinstrumentation. 1988, New York: John Wiley & Sons
- Kuhn PK, Christensen DA: Influence of Temperature Probe Sheathing Materials During Ultrasonic Heating. IEEE Trans Biomed Eng. 1986, BME-33: 536-538.View Article
- Vickers D: Hyperthermia for treatment of hyperplasia in a hemodialysis access graft. MS thesis. 2008, University of Utah, Department of Bioengineering
- Kunz K, Luebbers R: The finite difference time domain method for electromagnetics. 1993, Boca Raton: CRC Press
- Iijima F, Tsuchiya T, Nobuyki E: Analysis of characteristics of underwater sound propagation in the ocean by finite difference time domain method. Jpn J Appl Phys. 2000, 39: 3200-3204. 10.1143/JJAP.39.3200.View Article
- Hallaj I, Cleveland R: FDTD simulation of finite-amplitude pressure and temperature fields for biomedical ultrasound. J Acoust Soc Am. 1999, 105: L8-L12.
- Vasseur J, Deymier P, Chenni B, Djafari-Rouhani B, Dobrzynski L, Prevost D: Experimental and theoretical evidence for the existence of absolute acoustic band gaps in two-dimensional solid phononic crystals. Phys Rev Letters. 2001, 86: 3012-3015. 10.1103/PhysRevLett.86.3012.View Article
- Pennes HH: Analysis of tissue and arterial blood temperatures in the resting human forearm. (1948). J Appl Physiol. 1998, 85: 5-34.PubMed
- Brinton M: Ultrasound-induced thermal therapy of hyperplasia in hemodialysis access grafts. MS thesis. 2010, University of Utah, Department of Electrical and Computer Engineering
- Rayleigh JWS: The Theory of Sound. 1945, NewYork: Dover, 2: 107-
- Duck F: Physical Properties of Tissues: A Comprehensive Reference Book. 1990, San Diego: Academic Press
- Incropera F, DeWitt D: Introduction to Heat Transfer. 1996, New York: John Wiley & Sons, 3
- Paszkowiak J, Dardik A: Arterial wall shear stress: observations from the bench to the bedside. Vasc Endovascular Surg. 2003, 37: 47-57. 10.1177/153857440303700107.View ArticlePubMed
- Wissler EH: Pennes' 1948 Paper Revisited. J Appl Physiol. 1998, 85: 35-41.PubMed
- Overgaard J, Gonzalez D, Hulshof M, Arcangeli G, Dahl O, Mella O, Bentzen S: Hyperthermia as adjuvant to radiation therapy of recurrent or metastatic malignant melanoma. A multicentre randomized trial by the European Society for Hyperthermic Oncology. Int Journal Hyp. 1996, 12: 3-20. 10.3109/02656739609023685.View Article
- Corry P, Spanos W, Tilchen E, Barlogie B, Barkley H, Armour E: Combined ultrasound and radiation therapy treatment of human superficial tumors. Radiology. 1982, 145: 165-169.View ArticlePubMed
- Daftari I, Barash D, Lin S, O'Brien J: Use of high-frequency ultrasound imaging to improve delineation of anterior uveal melanoma for proton irradiation. Phys Med Biol. 2001, 46: 579-590. 10.1088/0031-9155/46/2/322.View ArticlePubMed
- Chen L, Dyson M, Rymer J, Bolton P, Young S: The use of high-frequency diagnostic ultrasound to investigate the effect of hormone replacement therapy on skin thickness. Skin Res Tech. 2001, 7: 95-97. 10.1034/j.1600-0846.2001.70207.x.View Article
- Ketis N, Hoover R, Karnovsky M: Effects of hyperthermia on cell survival and patterns of protein synthesis in endothelial cells from different origins. Cancer Res. 1988, 48: 2101-2106.PubMed
- Hikim APS, Lue Y, Yamamoto CM, Vera Y, Rodriguez S, Yen P, Soeng K, Wang C, Swerdloff R: Key apoptotic pathways for heat-induced programmed germ cell death in the testis. Endocrinology. 2003, 144: 3167-3175. 10.1210/en.2003-0175.View ArticlePubMed
- Cohen JJ, Duke RC, Fadok VA, Sellins KS: Apotosis and programmed cell death in immunity. Annu Rev Immunol. 1992, 10: 267-93. 10.1146/annurev.iy.10.040192.001411.View ArticlePubMed
- Samali A, Holmberg C, Sistonen L, Orrenius S: Thermotolerance and cell death are distinct cellular responses to stress: dependence on heat shock proteins. FEBS Lett. 1999, 461: 306-10. 10.1016/S0014-5793(99)01486-6.View ArticlePubMed
- Harmon B, Corder A, Collins R, Gobé G, Allen J, Allan D, Kerr J: Cell death induced in a murine mastocytoma by 42-47°C heating in vitro: evidence that the form of death changes from apoptosis to necrosis above a critical heat load. Int J Radiat Biol. 1990, 58: 845-858. 10.1080/09553009014552221.View ArticlePubMed
- Hehrlein C, Chuang CH, Tuntelder JR, Tatsis GP, Littmann L, Svenson RH: Effects of vascular runoff on myointimal hyperplasia after mechanical balloon or thermal laser arterial injury in dogs. Circulation. 1991, 84: 884-890.View ArticlePubMed
- DiBernardo BE, Reyes J, Chen B: Evaluation of tissue thermal effects from 1064/1320-nm laser-assisted lipolysis and its clinical implications. J Cosmet Laser Ther. 2009, 11: 62-69. 10.1080/14764170902792181.View ArticlePubMed
- Morse PM, Ingard KU: Theoretical Acoustics. 1968, Princeton: Princeton University Press, 278-279.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.