Calculation of arterial wall temperature in atherosclerotic arteries: effect of pulsatile flow, arterial geometry, and plaque structure
- Obdulia Ley1Email author and
- Taehong Kim1
https://doi.org/10.1186/1475-925X-6-8
© Ley and Kim; licensee BioMed Central Ltd. 2007
Received: 29 September 2006
Accepted: 01 March 2007
Published: 01 March 2007
Abstract
Background
This paper presents calculations of the temperature distribution in an atherosclerotic plaque experiencing an inflammatory process; it analyzes the presence of hot spots in the plaque region and their relationship to blood flow, arterial geometry, and inflammatory cell distribution. Determination of the plaque temperature has become an important topic because plaques showing a temperature inhomogeneity have a higher likelihood of rupture. As a result, monitoring plaque temperature and knowing the factors affecting it can help in the prevention of sudden rupture.
Methods
The transient temperature profile in inflamed atherosclerotic plaques is calculated by solving an energy equation and the Navier-Stokes equations in 2D idealized arterial models of a bending artery and an arterial bifurcation. For obtaining the numerical solution, the commercial package COMSOL 3.2 was used. The calculations correspond to a parametric study where arterial type and size, as well as plaque geometry and composition, are varied. These calculations are used to analyze the contribution of different factors affecting arterial wall temperature measurements. The main factors considered are the metabolic heat production of inflammatory cells, atherosclerotic plaque length l p , inflammatory cell layer length l mp , and inflammatory cell layer thickness d mp .
Results
The calculations indicate that the best location to perform the temperature measurement is at the back region of the plaque (0.5 ≤ l/l p ≤ 0.7). The location of the maximum temperature, or hot spot, at the plaque surface can move during the cardiac cycle depending on the arterial geometry and is a direct result of the blood flow pattern. For the bending artery, the hot spot moves 0.6 millimeters along the longitudinal direction; for the arterial bifurcation, the hot spot is concentrated at a single location due to the flow recirculation observed at both ends of the plaque. Focusing on the thermal history of different points selected at the plaque surface, it is seen that during the cardiac cycle the temperature at a point located at l/l p = 0.7 can change between 0.5 and 0.1 degrees Celsius for the bending artery, while no significant variation is observed in the arterial bifurcation. Calculations performed for different values of inflammatory cell layer thickness d mp indicate the same behavior reported experimentally; that corresponds to an increase in the maximum temperature observed, which for the bending artery ranges from 0.6 to 2.0 degrees Celsius, for d mp = 25 and 100 micrometers, respectively.
Conclusion
The results indicate that direct temperature measurements should be taken (1) as close as possible to the plaque/lumen surface, as the calculations show a significant drop in temperature within 120 micrometers from the plaque surface; (2) in the presence of blood flow, temperature measurement should be performed in the downstream edge of the plaque, as it shows higher temperature independently of the arterial geometry; and (3) it is necessary to perform measurements at a sampling rate that is higher than the cardiac cycle; the measurement should be extended through several cardiac cycles, as variations of up to 0.7 degrees Celsius were observed at l/l p = 0.7 for the bending artery.
Keywords
Background
Atherosclerosis refers to the thickening of the arterial wall; the lesion starts by adhesion and filtration of lipoproteins (LDL) and monocytes through the endothelium. The lesion grows as the monocytes are transformed into macrophages (inflammatory cells) and engulf LDL molecules to form foam cells. These factors, together with the migration of smooth muscle cells and the formation of cholesterol crystals, increase the arterial wall thickness which reduces the luminal area where the blood flow takes place.
A detailed description of the evolution, composition, and classification of atherosclerotic plaque can be found in [1, 2]. As plaque grows and changes its geometry and composition, the thin fibrous cap atheroma is disrupted due to alterations in local shear stress. The fibrous cap then breaks, releasing the plaque contents into the arterial lumen. When underlying tissue is exposed, a thrombus is formed which obstructs blood flow and produces ischemic syndromes such as stroke, angina pectoris, and myocardial infarctions (MI).
Clinical and postmortem anatomical studies reveal that atherosclerotic lesions in humans develop preferentially along the inner wall of curved segments and the outer walls of bifurcations of relatively large arteries [3, 4]. This is because the blood flow is disturbed by the occurrence of flow separation and formation of complex secondary and recirculation flows [5, 6]. These sites correspond to locations where the fluid shear stress on the vessel wall is significantly lower in magnitude than the normal physiological value [8–10]. This pattern of plaque localization is observed regardless of diet and ethnicity [4].
Cardiovascular disease (CVD) is the third major cause of death in the developed world. In the United States, about 1.1 million patients annually suffer from acute myocardial infarction (AMI) [11] and nearly two thirds of the acute coronary syndromes pass undetected [12]. As a result, early detection and follow-up of vulnerable plaque is a major challenge that cardiologists face today because it is crucial in the prevention of acute cardiac events [13].
Strong interest in understanding atherosclerosis in the scientific community is indicated by the number of publications in the area, up from about 6,000 in the 1980's to 18,000 in 2000–2005. Modeling studies underlying the physical factors of cardiovascular disease and atherosclerosis are classified in the following main areas: 1) the analysis of blood flow behavior to correlate between plaque location and wall shear stress [7, 9, 10]; 2) the calculation of circumferential stress distribution in atherosclerotic vessels to analyze the hypothesized direct correlation between maximal stress and plaque rupture [14, 15]; 3) the analysis of the transport and accumulation of macromolecules, such as LDL and albumin, on the arterial wall to explain deposit formation in terms of local wall shear stress variations [8, 16]; and 4) the study of the endothelial transport of inflammatory cells or transmigration, using confocal laser microscopy [17–19] and chemotaxis [20].
The fact that atherosclerosis is recognized as an inflammatory disease, along with the correlation observed between plaque temperature and rupture likelihood, is awakening interest in the determination of the arterial wall temperature (AWT) and the characterization of plaque vulnerability through local temperature [21]. It is reported [22] that plaques have several regions with surface temperatures that vary reproducibly by 0.2–0.3°C, and 37% have substantially warmer regions with temperature variations between 0.4 and 2.2°C. Such hot spots appear in plaques due to the increased transport and activation of macrophages or inflammatory cells [1]. These penetrate the arterial wall, form a thin layer, and (due to their high metabolic activity) produce the observed thermal inhomogeneity or hot spot.
This paper demonstrates calculations of inflamed atherosclerotic plaque temperature profiles. These calculations reveal information about unstable plaques and are used to correlate plaque geometry, activity, and evolutionary stage to temperature profiles that can be measured using novel catheters [13]. These calculations are used to study the factors affecting the arterial wall temperature measurements taken using such catheters [12, 22–24]. The information that can be extracted from the model and calculations presented here is: the best location at which to perform the measurement, an adequate measurement sampling rate in relation to cardiac cycle, and the cooling effect of blood flow during measurement. Information related to temperature gradients calculated with the model proposed herein can serve clinical research groups studying AWT [23, 24] because there are divergent opinions about the cooling effect of blood flow and other parameters which affect direct temperature measurement.
The numerical determination of arterial wall temperature distribution is necessary in order to: 1) observe the effect of vessel geometry; 2) perform parametric studies rapidly; 3) reduce animal experimentation; and 4) improve measurement strategies involving AWT contact methods. In this paper, the arterial wall temperature distribution is calculated by solving the Navier-Stokes equations and the energy equation in the arterial lumen and arterial wall. In the plaque, a heat generation term (due to the metabolic heat released by the macrophages or inflammatory cells) is considered.
Vessel geometries considered and corresponding geometrical parameters. The dimensions presented are shown in millimeters.
Recently, a similar study dealing with the estimation of plaque temperature has been published [34]. The present work, however, considers detailed plaque structures based on anatomical descriptions of atherosclerotic deposits [1]. Consideration of anatomically correct plaques is necessary in order to observe the combined effect of the convective cooling associated with blood flow as well as the metabolic production in the macrophage layer.
Methods
Blood flow in arteries
In this study, blood is treated as a homogeneous and incompressible Newtonian fluid. These assumptions are justified for the large arteries (where the blood components are small compared with the vessel diameter) and are also based on studies showing that the use of different viscosity models has only small effects on resulting blood flow [33, 35]. Fluid-Structure interaction (FSI) is also neglected as recent studies indicate that neglecting the elastic nature of the arterial wall incurs an error of approximately 10% in the velocity profiles [36–38].
The transient calculations of arterial wall temperature and blood flow are carried out under the different physiological flows for the two arteries selected (bending artery and arterial bifurcation), but the vessel geometry and waveform effects of blood flows are not tested independently. The pulsatile blood flow is governed by the Navier-Stokes equations and the continuity equation for an incompressible fluid:
∇·v = 0, (2)
Pulsatile flow waveform used at the inlet of (a) the bending artery [35] and (b) the arterial bifurcation [41] [42].
Temperature prediction at the arterial wall
Given the relationship between plaque vulnerability and local temperature, this paper shows calculations of temperature fields in large arteries containing small active plaques. To determine the temperature of the arterial wall and inflamed plaque, the atherosclerotic artery is modeled along the longitudinal direction. In this system, convection due to the luminal blood flow and heat conduction through the vessel walls are present, and the energy equation for the lumen and arterial wall regions takes the following forms
Thermophysical parameters of blood, arterial wall, plaque tissue and macrophage layer
Blood | Arterial wall | Plaque | Macrophage layer | |
---|---|---|---|---|
k i (W/m°C) | 0.549 | 0.476 | 0.484 | 0.484 |
ρ i (kg/m3) | 1050 | 1075 | 920 | 920 |
Cp i (J/kg°C) | 4390 | 3490 | 4080 | 4080 |
μ i (Pa·s) | 0.0033 |
For the thermal model, it is assumed that different layers in the vessel wall (intima, media, and adventitia) have the same thermal properties. The boundary conditions used to solve the equations (3) and (4) correspond to: (1) fixed temperature at the external vessel wall to account for the blood perfusion at the vasa vasorum (T = T a ); (2) fixed blood and tissue temperature at the vessel inlet (T = T a ), and (3) no temperature gradient at the vessel outlet ( = 0). T a is a constant that represents the arterial or core temperature and is assigned a value of 37.5°C. Finally, continuity of heat flux and temperature at the lumen/plaque interface, the plaque/vessel wall interface, and the plaque/macrophage layer interface is assumed.
The metabolic heat released by the inflamed plaque is a direct function of plaque composition and the developmental stage of the lesion. Macrophages are involved in all evolutionary stages of the lesion, but their activation and subsequent metabolic heat production varies in each stage. Macrophages embedded in atherosclerotic plaques are highly active as they are in the process of engulfing lipid molecules [2]. However, there are no reports of cellular metabolic heat q cell for active macrophages in atherosclerotic plaques. Other highly active cells are hepatocytes, for which microcalorimetry studies report a basal heat production of 300 pW per cell [46]. Microcalorimetry is commonly used to follow cell development and activation [47, 48], and it has shown that the metabolic heat production of cells is dependent on their concentration in the sample. For macrophages, cell volume as well as oxygen consumption increase proportionally [49, 50]. In this study, m was chosen to produce a maximum temperature change comparable with the reported measurement [21, 51]. For the calculations, three different values for the heat generation of the macrophage layer are used and they correspond to m = 0.05,0.1 and 0.2 W/mm3, these values are higher than the reported values of macrophages located in the lung of rabbits; higher values are assumed to consider the continuous involvement of the macrophages in engulfing lipid substances inside plaque. The value of m for the macrophage layer was approximated considering the cellular volume and metabolic heat production of a single cell by the following relationship
where V cell is the volume of a single cell and q cell is the heat produced by a single cell. It is assumed that macrophages have a volume of 30 to 60 μm3 and a cellular heat production q cell similar to that of hepatocytes.
Plaque composition, size and distribution
Plaque geometry and dimensions. d w is the arterial wall thickness, d p is the plaque thickness, d mp is the macrophage rich layer thickness, l f is the thickness of the fibrous cap. l p and l mp represent the extension or length of the plaque and the macrophage layer in the longitudinal direction, respectively. The location x i = x/l p of 4 observation points is presented. Point A corresponds to the center of the plaque (x A = 0.5), point B presents a point at x B = 0.65; point C is defined at x C = 0.7; and point D is the ending point of the plaque (x D = 1).
Dimensions of plaque and macrophage (MC) layer used for bending artery and arterial bifurcation (unit: μm)
Vessel Type | Plaque length l p | Wall thickness d w | MC layer length l mp | Plaque thickness d p | Fibrous cap thickness l f | MC layer thickness d mp |
---|---|---|---|---|---|---|
Bending | 3,000 | 270 | 1,500 | 540 | 25 | 50 |
50 | 25 | |||||
50 | 50 | |||||
50 | 100 | |||||
100 | 50 | |||||
Bifurcation | 3,000 | 335 | 1,500 | 670 | 25 | 50 |
50 | 25 | |||||
50 | 50 | |||||
50 | 100 | |||||
100 | 50 |
Calculations
Computations and mesh generation
Navier-Stokes equations and an energy equation for the two-dimensional (2D) transient case is solved using the multi-physics software COMSOL 3.2 (COMSOL Inc., SWEDEN). The software provides a complete FEA package with a unique equation-based modeling approach that allows for fully coupled multi-physics modeling in 2D and 3D. The FEM discretization of the time-dependent PDE problem produces a system of ordinary differential equations (ODE) or differential algebraic equations (DAE); COMSOL uses a version of the DAE solver DASPK [54] and a Newton-type iterative method.
The grid used is a non-structured mesh that is finer near the plaque and macrophage layer regions. The average number of triangular elements used in the calculation is 32,000. The number of elements used for the mesh is increased systematically to check for mesh independence; the mesh refinement is stopped when the maximum relative difference for the calculated parameter P for the different meshes satisfied the following condition
particularly, the maximum difference in the calculated temperature is 0.003% for the bending artery and 0.1% for the arterial bifurcation.
In addition to a grid independence analysis, a convergence verification for the initial condition is performed; this is needed to verify that the initial temperature used in the calculations corresponds to a stable temperature distribution. The calculations using the steady state temperature distribution are calculated using a velocity equal to the average velocity during one cardiac cycle (Fig. 2). Calculations are extended for several cycles using the waveform functions of Figure 2 as the inlet velocity condition; the number of cycles is selected to ensure that the temperature at several points along the plaque/lumen interface do not vary considerably between two consecutive cycles. For the arteries studied here, 6 cycles are necessary; the difference in the temperature distribution between the fifth and sixth cycles is 1.08% for the bending artery and 0.96% for the arterial bifurcation.
Temperature calculations
In this study, emphasis is given to the characterization of the temperature inhomogeneity considering the following parameters: 1) the arterial geometry and plaque location (Figure 1), 2) the extension of the macrophage layer (d mp ), and 3) the heat generation produced by the macrophage layer ( m ). This paper considers variations of d mp and m during a cardiac cycle because steady state calculations indicate that these parameters have more influence over the plaque temperature [55].
The plots presented herein indicate the temperature along the plaque/lumen surface and its variation with time. In these plots, the horizontal axis coincides with the base of the plaque, and the temperature at the plaque/lumen interface is plotted with respect to its horizontal location. In order to compare the temperature distribution at the plaque/lumen interface of the different arteries studied, the horizontal position is divided by the plaque length (l p ). The other type of plots presented here correspond to the time variation of the temperature in a selected point at the plaque/lumen surface.
As an initial condition of the transient calculation, the steady state solution is calculated with the uniform velocities of 0.08 m/s and 0.096 m/s for the bending artery and the arterial bifurcation, respectively. These values correspond to the average velocity during the cardiac cycle in each arterial geometry. The calculations are then performed for 6 cycles to attain a stationary state as discussed in the previous section.
Transient temperature distribution at the plaque/lumen interface for the bending artery at different times, produced by variations in the local heat generation in the plaque m (W/mm3). The times presented correspond to t = 4.0, 4.4 and 4.8 s. The inlet velocity profile for this arterial geometry is given in Figure 2a.
Transient temperature distribution at the plaque/lumen interface for the arterial bifurcation at different times. Temperature distribution is produced by variations in the local heat generation in the plaque m (W/mm3). The times indicated correspond to t = 4.0,4.4,4.6 and 4.8 s. The inlet velocity profile for this arterial geometry is given in Figure 2b.
For the bending artery, the representative times are chosen as: a) t = 4.0 and 5.0 s, which represent the beginning and the end of the cardiac cycle; b) t = 4.4 s, where the peak of a small forward flow occurs; and d) t = 4.8 s or a time approximately halfway through the rapid forward flow. For the arterial bifurcation, on the other hand, the representative times selected correspond to t = 4.0,4.4,4.6, and 4.8 seconds. As seen in Figure 2b, the velocity waveform for the bifurcation does not present significant variation after t = 4.5 s; however, t = 4.4 s is chosen instead of t = 4.2 s because the pulse propagation of the extreme peak at t = 4.2 s reaches the plaque region at t = 4.4 s due to the distance between the vessel inlet and the plaque location. As observed in the steady-state case [55], as the blood velocity increases, the temperature at the plaque/lumen interface decreases and vice-versa.
Results and discussion
Figure 4, shows a temporal temperature change of between 0.2°C and 1.1°C for q m = 0.05 W/mm3 and q m = 0.2 W/mm3 respectively. The maximum temperature peaks vary proportionally to the inlet blood flow, i.e. ΔT max is reduced at t = 4.4s and increases after t = 4.8 s. The maximum spatial temperature change during the cardiac cycle is located between 0.5 <l/l p < 0.7 (the solid line represents the location of ΔT max at t = 4.0 s and the dashed line indicates ΔT max at t = 4.4 s). It is observed that spatial ΔT max occurs closer to the center of the plaque when the blood velocity is reduced; it moves to the downstream edge of the plaque as blood flow increases in the last part of the cardiac cycle.
Flow circulation observed around an inflamed plaque located in an arterial bifurcation. Results shown correspond to calculations at times t = 4.0, 4.4, 4.6 and 4.8 s. The inlet velocity profile for this arterial geometry is given in Figure 2b.
Plots of ΔT(t) - ΔT(t = 0) along the plaque/lumen interface for the arterial bifurcation at different times t, for various metabolic heat production values corresponding to m = 0.2, 0.1 and 0.05 W/mm3 and inlet velocity profile of Figure 2b. The times indicated correspond to the times t = 4.4,4.6 and 4.8 s.
As in Figure 5, the vertical solid line indicates the location of spatial ΔT max at the beginning of the cycle (t = 4.0). It is seen that in the first part of the cycle (t < 4.6 s) the temperature at the front edge of the plaque drops due to convective cooling produced by an increase in blood flow velocity; as the blood velocity drops after t = 4.6 s, the temperature at the front edge of the plaque/lumen interface increases. The variation in the temperature distribution at the front edge of the plaque results from the flow circulation seen in Figure 6.
From Figures 4 and 5, three thermal regions are observed at the plaque/lumen interface and they are: (1) a front region where temperature is not affected by the presence of the macrophage layer and which changes according to the blood flow; (2) a middle region where temperature is directly affected by the macrophage layer dimensions and local convective heat transfer; and (3) a rear region where the convective heat transfer due to the blood flow is dominant over the presence of the macrophage layer. To study the regional variations of the plaque/lumen temperature, four points (A, B, C and D) are defined on the plaque lumen surface. Point A corresponds to the center of the plaque, point B presents the ending point of the middle region, point C is the beginning point of the rear region, and point D is the ending point of the plaque. The relative location of these points for the two arterial geometries studied is indicated in Figure 3.
Variation over time of plaque/lumen temperature measured at four different points; the arterial geometry corresponds to the bending artery. The inlet velocity profile for this arterial geometry is given in Figure 2a, the points where the temperature is recorded correspond to points A through D shown in Figure 3.
Variation over time of plaque/lumen temperature measured at four different points; the arterial geometry corresponds to the arterial bifurcation. The inlet velocity profile for this arterial geometry is given in Figure 2b, the points where the temperature is recorded correspond to points A through D shown in Figure 3.
Variation over time of plaque/lumen temperature measured at four different points along the interface for the bending artery. This figure indicates temperature changes produced by variations in the macrophage thickness d mp . The inlet velocity profile for this arterial geometry is given in Figure 2a, the points where the temperature is recorded correspond to points A through D, located as indicated in Figure 3.
Temperature contours around an inflamed plaque located in a bending artery. Results shown correspond to two different times during the cardiac cycle: t = 4.4, and 4.6 s. Calculations were obtained using the inlet velocity profile of Figure 2a. In this figure, the macrophage layer dimensions are m = 0.1 W/mm3, l mp = 1,500 μm, d mp = 100 μm, d p = 540 μm, and l f = 50 μm.
Finally, from Figures 8, 9 and 11, it is seen that for a given arterial geometry, some locations at the plaque/lumen interface show significant temperature variations during the cardiac cycle. This observation is important to determine the sampling rate of the arterial wall temperature (AWT) measurements performed around the plaque region. In particular, it is seen that plaque temperature measurements should be performed in at least three locations; this is because temperature distributions are characterized by three thermal regions with different effects of convective and conductive heat transfer. Measurements should be taken during one or more cardiac cycles as temperature change is affected by the inlet pulsatile flow.
Conclusion
The presence of the macrophage layer inside the plaque produces the appearance of a hot spot at the plaque/lumen interface. This work presents a mathematical model that can be used as a tool to improve direct arterial wall temperature measurements and the detection of vulnerable plaque. The location of the hot spot varies during the cardiac cycle and it is observed to be located behind the apex of the plaque. The exact location of the hot spot depends on the arterial geometry, the flow velocity, the presence of flow separation and recirculation, and the extension of the macrophage layer. The parameters that have the most influence over the plaque temperature are the macrophage layer thickness d mp and the metabolic heat accounting for the macrophage population and the inflammatory stage. The magnitude of m is assigned the values 0.05,0.1 and 0.2 W/mm3 which produces a maximum temperature change comparable with the reported measurements [21, 51]. The calculations reproduce the linear relationships observed experimentally between the macrophage layer thickness d mp and the maximum plaque temperature ΔT max .
It is observed that the best location to measure plaque temperature in the presence of blood flow is between the middle and the far edge of the plaque, as the highest spatial temperature variation occurs within this region; particularly, the maximum temperature location occurs between between 0.5 <l/l p < 0.7 for the bending artery, and is steadily located at l/l p = 0.7 for the arterial bifurcation model. It is seen that the maximum spatial temperature location travels within 0.6 mm along the longitudinal direction for the bending artery during one cardiac cycle. Compared with the spatial resolution of 0.5 mm of thermography in specially designed catheters [12, 22, 24], the path of the maximum spatial temperature location can not be ignored while directly measuring AWT during the cardiac cycles. Considering the temporal temperature change calculated for the bending artery, direct measurement of arterial wall temperature should be taken very close to the plaque/lumen surface; the blood temperature variations of 0.1°C are observed within 120 μm from the surface of the plaque center during a cardiac cycle. It is seen that the transient temperature variation at the plaque/lumen interface is different depending on the arterial geometry and the velocity waveform; finally, in order to have representative local temperature measurements, the temperature should be recorded in the same location over at least one full cardiac cycle.
The model proposed herein can be used to analyze the cooling effect of blood over the plaque surface and the effect of plaque geometry, size, and location in the artery. However, to improve the predictive nature of the model, it is necessary to consider four significant improvements: (1) realistic geometry, (2) the presence of turbulence, (3) non-newtonian characteristics of blood, and (4) Fluid-Structure interaction.
The current work employs idealized geometries derived from anatomical descriptions. This approximation is performed because realistic geometries are exceedingly complex due to the irregular shape and heterogeneous structure of atherosclerotic plaque. The geometries in this study assume a constant diameter, a smooth and round luminal surface, and a smooth boundary between the plaque and macrophage layers. Future studies should contemplate three dimensional and realistic geometries, to estimate the predictive nature of the model against direct measurements, as well as the need for realistic geometry to analyze the factors affecting arterial wall temperature.
The presence of turbulence can induce a large impact for temperature distributions which may be detected near plaques; it is expected that this effect is larger in plaques located by the arterial bifurcation or segments where large curvatures are detected. This study shows that the effect of turbulence depends on the arterial geometry and the velocity waveform. Turbulence profoundly influences the mechanisms of heat and momentum transfer in a thermal boundary layer, so future studies should introduce different turbulence models [56].
Multi-phase systems need to be adopted to minimize the potential errors which arise from neglecting the variation of blood viscosity. The theory indicates heat transfer in blood-solid flow with components (i.e. erythrocytes, plasma, etc.) is influenced by the momentum transport process. The plaque regions with low Reynolds numbers could show a denser packing of blood particles by diffusion-dominated flow. Such dense blood particles appreciably affect the convective heat transfer between both blood flow and particle components. However, blood flow (assumed as Newtonian fluid) in large arteries can be a dilute flow which keeps a sufficiently large distance between individual particles in order to avoid significant particle-particle interactions [57]. Finally, this study neglects the interaction between the fluids and structures which relate to the deformation of the arterial wall as well as the stress distribution over the plaque. Because of the hot spot formation due to the inflammation experienced by the plaque, consideration of the fluid-structure interaction and the calculation of thermal stresses produced by the hot plaque seems to be the most important characteristics missing from the proposed model.
In future studies, this model will be used to estimate the errors associated to direct measurements of AWT using the basket catheters described in [23, 24]; these catheters are formed by an array of thermocouples and are used to measure the temperature alterations in the arterial wall due to plaque activity. When catheters are introduced, the flow is disturbed, reduced, or completely occluded. This in turn should affect the temperature at the plaque surface. The present model used to analyze the alterations in AWT during the blood flow interruption or reduction may further elucidate the convective cooling effect of blood flow during measurement using catheters.
Declarations
Authors’ Affiliations
References
- Virmani R, Kolodgie FD, Burke AP, Farb A, Schwartz SM: Lessons from sudden coronary death: A comprehensive morphological classification scheme for atherosclerotic lesions. Arterioscler Thromb Vasc Biol 2000, 20: 1262–1275.View ArticleGoogle Scholar
- Jessup W, Baoutina A, Kritharides L: Macrophages in cardiovascular disease. In Macrophage. Edited by: Burke B, Lewis CE. Oxford University Press; 2002:490–522.Google Scholar
- Kerber CW, Hecht ST, Knox K, Buxton RB, Meltzer HS: Flow dynamics in fatal aneurysm of the basilar artery. American Journal of Neuroradiology 1996, 17: 1417–1421.Google Scholar
- Malek AM, Alper SL, Izumo S: Hemodynamic shear stress and its role in atherosclerosis. JAMA 1999, 282: 2035–2042. 10.1001/jama.282.21.2035View ArticleGoogle Scholar
- Motomiya M, Karino T: Flow patterns in the human carotid artery bifurcation. Stroke 1984, 15: 50–56.View ArticleGoogle Scholar
- Wada S, Karino T: Theoretical prediction of low-density lipoproteins concentration at the luminal surface of an artery with a multiple bend. Annals of Biomedical Engineering 2002, 30: 778–791. 10.1114/1.1495868View ArticleGoogle Scholar
- Bharadvaj BK, Maron RF, Giddens DP: Steady flow in a model of the human carotid bifurcation, Part I – flow visualization. Journal of Biomechanics 1982, 15: 349–362. 10.1016/0021-9290(82)90057-4View ArticleGoogle Scholar
- Darbeau MZ, Lutz RJ, Collins WE: Simulated lipoprotein transport in the wall of branched arteries. ASAIO Journal 2000, 46: 669–678. 10.1097/00002480-200011000-00006View ArticleGoogle Scholar
- Gnasso A, Irace C, Carallo C, De Franceschi MS, Motti C, Mattioli PL, Pujia A: In vivo association between low wall shear stress and plaque in subjects with asymmetrical carotid atherosclerosis. Stroke 1997, 28: 993–998.View ArticleGoogle Scholar
- Zarnis CK, Giddens DP, Bharadvaj BK, Sottiurai VS, Mabon RF, Glagov S: Carotid bifurcation atherosclerosis: Quantitative correlation of plaque localization with flow velocity profiles and wall shear stress. Circulation Research 1983, 53: 502–514.View ArticleGoogle Scholar
- 2002 Heart and Stroke Statistical Update. Dallas, TX: American Heart Association; 2002.Google Scholar
- Naghavi M, Madjid M, Gul K, Siadaty MS, Litovsky S, Willerson JT, Casscells SW: Thermography basket catheter: In vivo measurement of the temperature of atherosclerotic plaques for detection of vulnerable plaques. Catheterization in Cardiovascular Interventions 2003, 59: 52–59. 10.1002/ccd.10486View ArticleGoogle Scholar
- Diamantopoulos L: Arterial wall thermography. Journal of Interventional Cardiology 2003, 16: 261–266. 10.1034/j.1600-0854.2003.8041.xView ArticleGoogle Scholar
- Baldewsing RA, Mastik F, Schaar JA, Serruys PW, van der Steen AFW: Young's modulus reconstruction of vulnerable atherosclerotic plaque components using deformable curves. Ultrasound In Medicine & Biology 2006, 32: 201–210. 10.1016/j.ultrasmedbio.2005.11.016View ArticleGoogle Scholar
- Tang D, Yang C, Zheng J, Woodard PK, Saffitz JE, Sicard GA, Pilgram TK, Yuan C: Quantifying effects of plaque structure and material properties on stress distributions in human atherosclerotic plaques using 3D FSI models. Journal of Biomechanical Engineering 2005, 127: 1185–1194. 10.1115/1.2073668View ArticleGoogle Scholar
- Stangeby DK, Ethier CR: Coupled computational analysis of arterial LDL transport. Effects of hypertension. Computer Methods in Biomechanics and Biomedical Engineering 2002, 5: 233–241. 10.1080/10255840290010733View ArticleGoogle Scholar
- Hashimoto K, Kataoka N, Nakamura E, Asahara H, Ogasawara Y, Tsujioka K, Kajiya F: Direct observation and quantitative analysis of spatiotemporal dynamics of individual living monocytes during transendothelial migration. Atherosclerosis 2004, 177: 19–27. 10.1016/j.atherosclerosis.2004.06.016View ArticleGoogle Scholar
- Andjelkovic AV, Zochowski MR, Morgan F, Pachter JS: Qualitative and quantitative analysis of monocyte transendothelial migration by confocal microscopy and three-dimensional image reconstruction. Vitro Cell Dev Biol 2001, 37: 111–120. Publisher Full Text 10.1290/1071-2690(2001)037<0111:QAQAOM>2.0.CO;2View ArticleGoogle Scholar
- Hsiai TK, Cho SK, Reddy S, Hama S, Navab M, Demer LL, Honda HM, Ho CM: Pulsatile flow regulates monocyte adhesion to oxidized lipid--induced endothelial cells. Arterioscler Thromb Vasc Biol 2001, 21: 1770–1776.View ArticleGoogle Scholar
- Hsiai TK, Cho SK, Honda HM, Hama S, Navab M, Demer LL, Ho CM: Endothelial cell dynamics under pulsating flows: Significance of high versus low shear stress slew rates. Annals of Biomedical Engineering 2002, 30: 646–656. 10.1114/1.1484222View ArticleGoogle Scholar
- Madjid M, Naghavi M, Malik BA, Litovsky S, Willerson JT, Casscells W: Thermal detection of vulnerable plaque. American Journal of Cardiology 2002, (Suppl 90):36L-39L. 10.1016/S0002-9149(02)02962-4Google Scholar
- Stefanadis C, Toutouzas K, Tsiamis E, Vavuranakis M, Kallikazaros I, Toutouzas P: Thermography of human arterial system by means of new thermography catheters. Catheterization and Cardiovascular Interventions 2001, 54: 51–58. 10.1002/ccd.1237View ArticleGoogle Scholar
- Diamantopoulos L, Liu X, Scheerder ID, Krams R, Li S, Cleemput JV, Desmet W, Serruys PW: The effect of reduced blood-flow on the coronary wall temperature: Are significant lesions suitable for intravascular thermography? European Heart Journal 2003, 24: 1788–1795. 10.1016/S0195-668X(03)00440-8View ArticleGoogle Scholar
- Stefanadis C, Toutouzas K, Tsiamis E, Mitropoulos I, Tsioufis C, Kallikazaros I, Pitsavos C, Toutouzas P: Thermal heterogeneity in stable human coronary atherosclerotic plaques is underestimated in vivo: The "cooling effect" of blood flow. Journal of the American College of Cardiology 2003, 41: 403–408. 10.1016/S0735-1097(02)02817-6View ArticleGoogle Scholar
- Rappitsch G, Perktold K: Computer simulation of convective diffusion processes in large arteries. Journal of Biomechanics 1996, 29: 207–215. 10.1016/0021-9290(95)00045-3View ArticleGoogle Scholar
- Stangeby DK, Ethier CR: Computational analysis of coupled blood-wall arterial LDL transport. Journal of Biomechanical Engineering 2002, 124: 1–8. 10.1115/1.1427041View ArticleGoogle Scholar
- Lutostansky EM, Karner G, Rappitsch G, Ku DN, Perktold K: Analysis of hemodynamic fluid phase mass transport in a separated flow region. Journal of Biomechanical Engineering 2003, 125: 189–196. 10.1115/1.1543547View ArticleGoogle Scholar
- Huang H, Virmani R, Younis H, Burke AP, Kamm RD, Lee RT: The impact of calcification on the biomechanical stability of atherosclerotic plaques. Circulation 2001, 103: 1051–1056.View ArticleGoogle Scholar
- Rioufol G, Gilard M, Finet G, Ginon I, Boschat J, Andre-Fouet X: Evolution of spontaneous atherosclerotic plaque rupture with medical therapy: long-term follow-up with intravascular ultrasound. Circulation 2004, 110: 2875–2880. 10.1161/01.CIR.0000146337.05073.22View ArticleGoogle Scholar
- Bank AJ, Versluis A, Dodge SM, Douglas WH: Atherosclerotic plaque rupture: a fatigue process? Med Hypotheses 2000, 55: 480–484. 10.1054/mehy.2000.1096View ArticleGoogle Scholar
- Asakura T, Karino T: Flow patterns and spatial distribution of atherosclerotic lesions in human coronary arteries. Circ Res 1990, 66: 1045–1066.View ArticleGoogle Scholar
- Perktold K, Resch M, Peter RO: Three-dimensional numerical analysis of pulsatile flow and wall shear stress in the carotid artery bifurcation. Journal of Biomechanics 1991, 24: 409–420. 10.1016/0021-9290(91)90029-MView ArticleGoogle Scholar
- Filipovic N, Kojic M: Computer Simulations of blood flow with mass transport through the carotid artery bifurcation. Theoret Appl Mech 2004, 31: 1–33.View ArticleMATHGoogle Scholar
- Ten Have AG, Gijsen FJH, Wentzel JJ, Slager CJ, van der Steen AFW: Temperature distribution in atherosclerotic coronary arteries: influence of plaque geometry and flow (a numerical study). Physics in Medicine and Biology 2004, 49: 4447–4462. 10.1088/0031-9155/49/19/001View ArticleGoogle Scholar
- Johnston BM, Johnston PR, Corney S, Kilpatrick D: Newtonian blood flow in human right coronary arteries: transient simulations. Journal of Biomechanics 2006, 39: 1116–1128. 10.1016/j.jbiomech.2005.01.034View ArticleGoogle Scholar
- Zeng D, Ding Z, Friedman MH, Ethier CR: Effects of cardiac motion on right coronary artery hemodynamics. Annals of Biomedical Engineering 2003, 31: 420–429. 10.1114/1.1560631View ArticleGoogle Scholar
- Liu Y, Lai Y, Nagaraj A, Kane B, Hamilton A, Greene R, McPherson DD, Chandran KB: Pulsatile flow simulation in arterial vascular segments with intravascular ultrasound images. Medical Engineering and Physics 2001, 23: 583–595. 10.1016/S1350-4533(01)00088-1View ArticleGoogle Scholar
- Zamir M: The Physics of Pulsatile Flow. New York: AIP Press; 2000.View ArticleMATHGoogle Scholar
- Nichols WW, O'Rourke MF: McDonald's Blood Flow in Arteries: Theoretical, Experimental and Clinical Principles. Oxford University Press; 1997.Google Scholar
- Perktold K, Rappitsch G: Computer simulation of local blood flow and vessel mechanics in a compliant carotid artery bifurcation model. Journal of Biomechanics 1995, 28: 845–856. 10.1016/0021-9290(95)95273-8View ArticleGoogle Scholar
- Chen J, Lu XY: Numerical investigation of the non-newtonian pulsatile blood flow in a bifurcation model with a non-planar branch. Journal of Biomechanics 2006, 39: 818–832. 10.1016/j.jbiomech.2005.02.003View ArticleGoogle Scholar
- Gijsen FJH, Allanic E, van de Vosse FN, Janssen JD: The influence of the non-newtonian properties of blood on the flow in large arteries: unsteady flow in a 90° curved tube. Journal of Biomechanics 1999, 32: 705–713. 10.1016/S0021-9290(99)00014-7View ArticleGoogle Scholar
- Duck FA: Physical Properties of Tissue: A Comprehensive Reference Book. San Diego: Academic Press; 1990.Google Scholar
- Welch AJ, van Gemert MJC: Optical Thermal Response of Laser-Irradiated Tissue. New York: Plenum; 1995.View ArticleGoogle Scholar
- Valvano JW, Chitsabesan B: Thermal conductivity and diffusivity of arterial wall and atherosclerotic plaque. Lasers in the Life Sciences 1987, 1: 219–229.Google Scholar
- Thoren SA, Monti M, Holma B: Heat conduction microcalorimetry of overall metabolism in rabbit alveolar macrophages on monolayers and in suspensions. Biochimica et Biophysica Acta 1990, 1033: 305–310.View ArticleGoogle Scholar
- Gebreselassie D, Schoen A: Energy metabolism of non-transformed and benzpyrene-transformed 3T3 cells: A microcalorimetric study. Biochimica et Biophysica Acta 1993, 1143: 211–214. 10.1016/0005-2728(93)90145-6View ArticleGoogle Scholar
- Okuda A, Kimura G: Factors affecting heat production of rat 3Y1 flbroblasts in flow microcalorimetry. Cell Struct Funct 1988, 13: 97–104.View ArticleGoogle Scholar
- Bidani A, Brown SE, Heming TA: pH i regulation in alveolar macrophages: Relative roles of Na + - H + antiport and H + -ATPase. American Journal of Physiology-Lung Cellular and Molecular Physiology 1994, 266: L681-L688.Google Scholar
- Bidani A, Brown SE: ATP dependent pHi recovery in macrophages: Evidence for a plasma membrane H + -ATPase. American Journal of Physiology-Cell Physiology 1990, 259: C586-C598.Google Scholar
- Casscells W, Hathorn X, David M, Krabach T, Vaughn WK, McAllister HA, Bearman G, Willerson JT: Thermal detection of cellular infiltrates in living atherosclerotic plaques: Possible implications for plaque rupture and thrombosis. Lancet 1996, 347: 1447–1449. 10.1016/S0140-6736(96)91684-0View ArticleGoogle Scholar
- Verheye S, De Meyer GRY, Van Lengenhove G, Knaapen MWM, Kockx MM: In vivo temperature heterogeneity of atherosclerotic plaques is determined by plaque composition. Circulation 2002, 105: 1596–1601. 10.1161/01.CIR.0000012527.94843.BFView ArticleGoogle Scholar
- Bhatia V, Bhatia R, Dhindsa S, Dhindsa M: Imaging of the vulnerable plaque: New modalities. Southern Medical Journal 2003, 96: 1142–1147. 10.1097/01.SMJ.0000089063.76530.04View ArticleGoogle Scholar
- Brown PN, Hindmarsh AC, Petzold LR: Using Krylov methods in the solution of large-scale differential-algebraic systems. SIAM J Sci Comput 1994, 15: 1467–1488. 10.1137/0915088View ArticleMATHGoogle Scholar
- Ley O, Kim T: Numerical prediction of atherosclerotic plaque temperature as a function of plaque size and composition. Proceedings of IMECE: 2005, ASME 2005.Google Scholar
- Kays WM, Crawford ME: Convective Heat and Mass Transfer. New York: McGraw-Hill; 1993.MATHGoogle Scholar
- Naterer GF: Heat Transfer in Single and Multiphase Systems. Boca Raton, Florida: CRC Press; 2003.Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.