Influence of catheter design on lumen wall temperature distribution in intracoronary thermography

Intracoronary thermography is a currently used vulnerable plaque detection method. We studied how catheter design and catheter location influence the temperature readings, and thus its capacity to detect vulnerable plaques. Finite element calculations were performed on geometries representing the coronary artery, the vulnerable plaque and the catheter. Catheter material, diameter and location with respect to the plaque were varied. Both flow and no-flow situations were studied. Maximal lumen wall temperature difference without a catheter (DeltaT=0.12 degrees C, flow=75 cm(3) min(-1)) was considered the reference. Presence of a 1.0mm nitinol catheter right under the plaque increased DeltaT to 0.14 degrees C, whereas a 1.0 mm polyurethane catheter increased DeltaT to 0.51 degrees C. The location at which a thermosensitive element should be placed for most optimal temperature readings during a pullback was shown to lie at the catheter edge for the nitinol catheter and at 1.1mm from the catheter edge for the polyurethane catheter. Temperature readings decreased to background temperature when the catheter was in close proximity but not overlapping the plaque. DeltaT decreased approximately by 70% when a gap of 0.2 mm existed between the catheter and the lumen wall. Occlusion of blood flow increased DeltaT values in all cases, but most pronounced for nitinol catheters. A polyurethane catheter increased the temperature readings, since its heat conductivity is lower than that of blood, which makes it a very good choice for heat source detection. Catheter design can contribute to enhanced temperature readings and thus can enable more optimal vulnerable plaque detection.     

Calculation of arterial wall temperature in atherosclerotic arteries: effect of pulsatile flow,arterial geometry,and plaque structure

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.     

Analysis of oxygen transport from pulsatile, viscous blood flow to diseased coronary arteries of man

Oxygen transport to multiple “non-obstructive” plaque regions in main coronary arteries of man was examined by numerically solving the oxygen transport equation for convective and diffusive processes in the lumen for actual variations of blood flow rate and the velocity field during the cardiac cycle. Oxygen transport to the wall varied significantly along the arterial section, was strongly dependent upon the various flow regions that occurred, and varied considerably during the cardiac cycle. A drastic reduction in oxygen transport to the arterial wall occurred at the incipient separation location on the back side of a plaque where it is believed that the lumen side resistance to oxygen transport is at least an order of magnitude greater than the inner avascular wall resistance, and therefore the availability of oxygen for cellular respiration is essentially boundary layer controlled. In vivo measurements with oxygen microelectrodes in animals are needed to learn more about variations of oxygen transport in plaque regions, in particular on the back side of plaques where hypoxia may occur.     

Computer model of fetal-maternal heat exchange in sheep

We constructed and used a mathematical model of maternal-fetal heat exchange in the sheep to explore the effects of changes in certain parameters on steady-state fetal temperatures and to determine whether the fetus in the model has any potential to control its own temperature. The model took into account both fetal and placental heat production and exchange of heat in the placenta, across the fetal skin, via amniotic fluid, and through the uterine wall. The maternal ewe was assumed to be a constant temperature heat sink. Changes in placental or fetal heat production were calculated to change the ratio of heat exiting across the placenta or fetal skin significantly but to have little effect on fetal core temperature, e.g., a rise of only 0.8 degrees C was predicted after a twofold increase in fetal heat production. Fetal placental blood flow was calculated to affect fetal temperature the most of any flow, a reduction to zero causing fetal temperature to rise 5.0 degrees C. Changes in heat conductances between fetal skin and amniotic fluid, or between amniotic fluid and uterine wall, had minimal effect on fetal temperature. From the model calculations here and because heat exchange within the sheep placenta has previously been calculated to be extremely efficient, we conclude that the fetal sheep has little ability to control its temperature by changes in heat dissipated through extraplacental pathways. Thus the model predicts an effective heat clamp that closely links fetal to maternal temperature.     

Theoretical study on flow-dependent concentration polarization of low density lipoproteins at the luminal surface of a straight artery

It is suspected that physical and fluid mechanical factors play important roles in the localization of atherosclerotic lesions and intimal hyperplasia in man by affecting the transport of cholesterol in flowing blood to arterial walls. Hence, we have studied theoretically the effects of various physical and fluid mechanical factors such as wall shear rate, diffusivity of low density lipoproteins (LDL), and filtration velocity of water at the vessel wall on surface concentration of LDL at an arterial wall by means of a computer simulation of convective and diffusive transport of LDL in flowing blood to the wall of a straight artery under conditions of a steady flow. It was found that under normal physiologic conditions prevailing in the human arterial system, due to the presence of a filtration flow of water at the vessel wall, flow-dependent concentration polarization (accumulation or depletion) of LDL occurs at a blood/endothelium boundary. The surface concentration of LDL at an arterial wall takes higher values than that in the bulk flow in that vessel, and it is affected by three major factors, that is, wall shear rate, gamma w, filtration velocity of water at the vessel wall, Vw, and the distance from the entrance of the artery, L. It increases with increasing Vw and L, and decreasing gamma w hence the flow rate. Thus, under certain circumstances, the surface concentration of LDL could rise locally to a value which is several times higher than that in the bulk flow, or drop locally to a value even lower than a critical concentration for the maintenance of normal functions and survival of cells forming the vessel wall. These results suggest the possibility that all the vascular phenomena such as the localization of atherosclerotic lesions and intimal hyperplasia, formation of cerebral aneurysms, and adaptive changes of lumen diameter and wall structure of arteries and veins to certain changes in hemodynamic conditions in the circulation are governed by this flow-dependent concentration polarization of LDL which carry cholesterol.     

Multiphysics simulation of blood flow and LDL transport in a porohyperelastic arterial wall model

Atherosclerosis localizes at a bend andor bifurcation of an artery, and low density lipoproteins (LDL) accumulate in the intima. Hemodynamic factors are known to affect this localization and LDL accumulation, but the details of the process remain unknown. It is thought that the LDL concentration will be affected by the filtration flow, and that the velocity of this flow will be affected by deformation of the arterial wall. Thus, a coupled model of a blood flow and a deformable arterial wall with filtration flow would be invaluable for simulation of the flow field and concentration field in sequence. However, this type of highly coupled interaction analysis has not yet been attempted. Therefore, we performed a coupled analysis of an artery with multiple bends in sequence. First, based on the theory of porous media, we modeled a deformable arterial wall using a porohyperelastic model (PHEM) that was able to express both the filtration flow and the viscoelastic behavior of the living tissue, and simulated a blood flow field in the arterial lumen, a filtration flow field and a displacement field in the arterial wall using a fluid-structure interaction (FSI) program code by the finite element method (FEM). Next, based on the obtained results, we further simulated LDL transport using a mass transfer analysis code by the FEM. We analyzed the PHEM in comparison with a rigid model. For the blood flow, stagnation was observed downward of the bends. The direction of the filtration flow was only from the lumen to the wall for the rigid model, while filtration flows from both the wall to the lumen and the lumen to the wall were observed for the PHEM. The LDL concentration was high at the lumenwall interface for both the PHEM and rigid model, and reached its maximum value at the stagnation area. For the PHEM, the maximum LDL concentration in the wall in the radial direction was observed at the position of 3% wall thickness from the lumenwall interface, while for the rigid model, it was observed just at the lumenwall interface. In addition, the peak LDL accumulation area of the PHEM moved about according to the pulsatile flow. These results demonstrate that the blood flow, arterial wall deformation, and filtration flow all affect the LDL concentration, and that LDL accumulation is due to stagnation and the presence of filtration flow. Thus, FSI analysis is indispensable.     

Computational analysis of coupled blood-wall arterial LDL transport

The transport of macromolecules, such as low density lipoproteins (LDLs), across the artery wall and their accumulation in the wall is a key step in atherogenesis. Our objective was to model fluid flow within both the lumen and wall of a constricted, axisymmetric tube simulating a stenosed artery, and to then use this flow pattern to study LDL mass transport from the blood to the artery wall. Coupled analysis of lumenal blood flow and transmural fluid flow was achieved through the solution of Brinkman's model, which is an extension of the Navier-Stokes equations for porous media. This coupled approach offers advantages over traditional analyses of this problem, which have used possibly unrealistic boundary conditions at the blood-wall interface; instead, we prescribe a more natural pressure boundary condition at the adventitial vasa vasorum, and allow variations in wall permeability due to the occurrence of plaque. Numerical complications due to the convection dominated mass transport process (low LDL diffusivity) are handled by the streamline upwind/Petrov-Galerkin (SUPG) finite element method. This new fluid-plus-porous-wall method was implemented for conditions typical of LDL transport in a stenosed artery with a 75 percent area reduction (Peclet number=2 x 10(8)). The results show an elevated LDL concentration at the downstream side of the stenosis. For the higher Darcian wall permeability thought to occur in regions containing atheromatous lesions, this leads to an increased transendothelial LDL flux downstream of the stenosis. Increased transmural filtration in such regions, when coupled with a concentration-dependent endothelial permeability to LDL, could be an important contributor to LDL infiltration into the arterial wall. Experimental work is needed to confirm these results.     

Role of the immune response in atherosclerosis and acute coronary syndromes

Cardiovascular diseases represent one of the most important causes of death in the world. The underlying pathogenetic process is atherosclerosis which leads to the progressive reduction of the arterial lumen and therefore to the ischemia of the perfused organs. Atherogenesis results from the interaction between the biology of the arterial wall and the various stress stimuli present in the circulating blood. The first steps of atherogenesis occur very early, already during the fetal life. Those arterial segments that are subjected to the initiating causes (including hemodynamic stress) show altered endothelial permeability and allow the infiltration of macromolecules, like lipoproteins, in the subintimal space. At the same time, the smooth muscle cells that are subjected to the same local factors produce proteoglycans able to bind lipoproteins and to promote their oxidation. Oxidized lipoproteins induce the expression of chemokines and adhesion molecules on the luminal surface of the endothelium, which then allow the local recruitment of monocytes-macrophages and T lymphocytes. This is a local inflammatory process that, in theory, should contribute to reestablish the homeostasis of the vascular wall by promoting the elimination of injured tissue and its repair. Unfortunately, for unknown reasons, the immuno-inflammatory reaction persists and autoamplifies, the various components of the immune response finally contributing to the pathogenesis of atherosclerosis as well as of atherosclerotic complications.     

Radial heat and water transport across the airway wall

The thermal profiles in the airways of healthy human volunteers and patients with asthma differ after cessation of hyperpnea. The asthmatic patients rewarm their airways more rapidly. To identify thermal properties and processes that could account for the difference between these populations, we developed a model describing the radial transport of heat and water across the trachea. A distinctive feature of the model is a variable parameter describing blood supply to the mucosal and submucosal layers. Simulations performed with the model are initiated by a breath-hold maneuver and are propagative in time. Blood perfusion rates in the airway wall, the thickness of the layer of airway surface liquid, and the mucosa-submucosa thickness, all thought to be more pronounced in asthmatic patients, were varied by changing model parameters and initial conditions. Increasing the thickness of the liquid layer by more than an order of magnitude had little effect on the temperature or water content in the airway lumen. Doubling the blood flow to the mucosa-submucosa resulted in a slight increase in airway temperature. When this effect was coupled, however, with an increase in the thickness of the mucosa-submucosa layer, the increase in temperature was more pronounced. Because the bronchial circulation is the major source of heat to the airway, these results indicate that differences in airway wall thickness coupled with differences in the magnitude or responsiveness of the bronchial microcirculation could account for the differences in intra-airway temperature between the two populations.     

Influence of arterial wall compliance on the pressure drop across coronary artery stenoses under hyperemic flow condition

Hemodynamic endpoints such as flow and pressure drop are often measured during angioplasty procedures to determine the functional severity of a coronary artery stenosis. There is a lack of knowledge regarding the influence of compliance of the arterial wall-stenosis on the pressure drop under hyperemic flows across coronary lesions. This study evaluates the influence in flow and pressure drop caused by variation in arterial-stenosis compliance for a wide range of stenosis severities. The flow and pressure drop were evaluated for three different severities of stenosis and tested for limiting scenarios of compliant models. The Mooney-Rivlin model defined the non-linear material properties of the arterial wall and the plaque regions. The non-Newtonian Carreau model was used to model the blood flow viscosity. The fluid (blood)-structure (arterial wall) interaction equations were solved numerically using the finite element method. Irrespective of the stenosis severity, the compliant models produced a lower pressure drop than the rigid artery due to compliance of the plaque region. A wide variation in the pressure drop was observed between different compliant models for significant (90% area occlusion) stenosis with 41.0, 32.1, and 29.8 mmHg for the rigid artery, compliant artery with calcified plaque, and compliant artery with smooth muscle cell proliferation, respectively. When compared with the rigid artery for significant stenosis the pressure drop decreased by 27.7% and 37.6% for the calcified plaque and for the smooth muscle cell proliferation case, respectively. These significant variations in pressure drop for the higher stenosis may lead to misinterpretation and misdiagnosis of the stenosis severity.     

Computational study of pulsatile blood flow in prototype vessel geometries of coronary segments

The spatial and temporal distributions of wall shear stress (WSS) in prototype vessel geometries of coronary segments are investigated via numerical simulation, and the potential association with vascular disease and specifically atherosclerosis and plaque rupture is discussed. In particular, simulation results of WSS spatio-temporal distributions are presented for pulsatile, non-Newtonian blood flow conditions for: (a) curved pipes with different curvatures, and (b) bifurcating pipes with different branching angles and flow division. The effects of non-Newtonian flow on WSS (compared to Newtonian flow) are found to be small at Reynolds numbers representative of blood flow in coronary arteries. Specific preferential sites of average low WSS (and likely atherogenesis) were found at the outer regions of the bifurcating branches just after the bifurcation, and at the outer-entry and inner-exit flow regions of the curved vessel segment. The drop in WSS was more dramatic at the bifurcating vessel sites (less than 5% of the pre-bifurcation value). These sites were also near rapid gradients of WSS changes in space and time – a fact that increases the risk of rupture of plaque likely to develop at these sites. The time variation of the WSS spatial distributions was very rapid around the start and end of the systolic phase of the cardiac cycle, when strong fluctuations of intravascular pressure were also observed. These rapid and strong changes of WSS and pressure coincide temporally with the greatest flexion and mechanical stresses induced in the vessel wall by myocardial motion (ventricular contraction). The combination of these factors may increase the risk of plaque rupture and thrombus formation at these sites.     

Modeling the large-scale geometry of human coronary arteries

Two principles suffice to model the large-scale geometry of normal human coronary arterial networks. The first principle states that artery diameters are set to minimize the power required to distribute blood through the network. The second principle states that arterial tree geometries are set to globally minimize the lumen volume. Given only the coordinates of an arterial tree's source and "leaves", the model predicts the nature of the network connecting the source to the leaves. Measurements were made of the actual geometries of arterial trees from postmortem healthy human coronary arteriograms. The tree geometries predicted by the model look qualitatively similar to the actual tree geometries and have volumes that are within a few percent of those of the actual tree geometries. Human coronary arteries are therefore within a few percent of perfect global volume optimality. A possible mechanism for this near-perfect global volume optimality is suggested. Also, the model performs best under the assumption that the flow is not entirely steady and laminar.     

An analytical model of temperature regulation in human head

An analytical model of human brain temperature regulation is proposed. The model describes the distribution of brain temperature as a function of internal and external parameters, such as temperature of the incoming arterial blood, blood flow, oxygen consumption rate, ambient temperature, and heat exchange with the environment. It is shown that substantial changes in human brain temperature can be accomplished only through changes in the temperature of the incoming arterial blood or substantial suppression of blood flow. Other parameters can lead only to temperature changes near the brain surface.     

Is arterial wall-strain stiffening an additional process responsible for atherosclerosis in coronary bifurcations?: an in vivo study based on dynamic CT and MRI

Coronary bifurcations represent specific regions of the arterial tree that are susceptible to atherosclerotic lesions. While the effects of vessel compliance, curvature, pulsatile blood flow, and cardiac motion on coronary endothelial shear stress have been widely explored, the effects of myocardial contraction on arterial wall stress/strain (WS/S) and vessel stiffness distributions remain unclear. Local increase of vessel stiffness resulting from wall-strain stiffening phenomenon (a local process due to the nonlinear mechanical properties of the arterial wall) may be critical in the development of atherosclerotic lesions. Therefore, the aim of this study was to quantify WS/S and stiffness in coronary bifurcations and to investigate correlations with plaque sites. Anatomic coronary geometry and cardiac motion were generated based on both computed tomography and MRI examinations of eight patients with minimal coronary disease. Computational structural analyses using the finite element method were subsequently performed, and spatial luminal arterial wall stretch (LW(Stretch)) and stiffness (LW(Stiff)) distributions in the left main coronary bifurcations were calculated. Our results show that all plaque sites were concomitantly subject to high LW(Stretch) and high LW(Stiff), with mean amplitudes of 34.7 ± 1.6% and 442.4 ± 113.0 kPa, respectively. The mean LW(Stiff) amplitude was found slightly greater at the plaque sites on the left main coronary artery (mean value: 482.2 ± 88.1 kPa) compared with those computed on the left anterior descending and left circumflex coronary arteries (416.3 ± 61.5 and 428.7 ± 181.8 kPa, respectively). These findings suggest that local wall stiffness plays a role in the initiation of atherosclerotic lesions.     

Thermal behavior of biological tissues—A general analysis

A general analysis is presented for the thermal behavior of a biological tissue. Energy transport by the circulatory system is assumed to be represented by a modified Fick's law. General boundary conditions are assumed for the two-dimensional model and solutions are obtained for rectangular, cylindrical, and spherical geometries. The effects of blood perfusion rate, metabolic rate, arterial temperature and heat exchange with the environment are considered. Results indicate a region of almost constant temperature in the deeper layers of the tissue and reaffirm the important role which blood flow plays in maintaining homeostasis.     

The influence of vascular anatomy on carotid artery stenting: A parametric study for damage assessment

Carotid artery stenting is emerging as an alternative technique to surgery for the treatment of symptomatic severe carotid stenosis. Clinical and experimental evidence demonstrates that both plaque morphology and biomechanical changes due to the device implantation can be possible causes of an unsuccessful treatment. In order to gain further insights of the endovascular intervention, a virtual environment based on structural finite element simulations was built to emulate the stenting procedure on generalized atherosclerotic carotid geometries which included a damage model to quantify the injury of the vessel. Five possible lesion scenarios were simulated by changing both material properties and vascular geometrical features to cover both presumed vulnerable and stable plaques. The results were analyzed with respect to lumen gain and wall stresses which are potentially related to the failure of the procedure according to previous studies. Our findings show that an elliptic lumen shape and a thinner fibrous cap with an underlying lipid pool result in higher stenosis reduction, while large calcifications and fibrotic tissue are more prone to recoil. The shielding effect of a thicker fibrous cap helps to reduce local compressive stresses in the soft plaque. The presence of a soft plaque reduces the damage in the healthy vascular structures. Contrarily, the presence of hard plaque promotes less damage volume in the fibrous cap and reduces stress peaks in this region, but they seem to increase stresses in the media-intima layer. Finally the reliability of the achieved results was put into clinical perspective.     

Effect of continuous arterial blood flow in patients with rotary cardiac assist device on the washout of a stenosis wake in the carotid bifurcation: a computer simulation study

In recipients of rotary blood pumps for cardiac assist, the pulsatility of arterial flow is considerably diminished. This influences the shear stress patterns and streamlines in the arterial bed, with potential influence on washout and plaque growth. These effects may be aggravated in the recirculation area of stenoses, and therefore, exclude patients with atherosclerosis from the therapy with these devices. A numerical study was performed for the human carotid artery bifurcation with the assumption of a massive stenosis (75% reduction of cross-section area) in the carotid bulb. Four different flow time patterns (no support to full pump support) were applied. Flow patterns and particle residence time within the recirculation region were calculated, once within the relevant volume behind the stenosis and and once within a small region directly at the posterior heel of the stenosis. The flow patterns showed a considerable radial vorticity behind the stenosis. Mean particle residence time in the whole recirculation region was 15% less for high pump support (nearly continuous flow) compared to the natural flow pattern (0.19s compared to 0.22s), and nearly identical for the small heel region (0.28 to 0.27s). The flow simulation demonstrates, that even in the case of a pre-existing stenosis, the local effects of continuous flow on particle residence times are rather minimal (as was shown previously for intact arterial geometries). Therefore, from the point of macroscopic flow field analysis, continuous flow should not enhance the thromboembolic risk in ventricular assist device recipients.