Patient-Specific Artery Shrinkage and 3D Zero-Stress State in Multi-Component 3D FSI Models for Carotid Atherosclerotic Plaques Based on <i>In Vivo</i> MRI Data

Image-based computational models for atherosclerotic plaques have been developed to perform mechanical analysis to quantify critical flow and stress/strain conditions related to plaque rupture which often leads directly to heart attack or stroke. An important modeling issue is how to determine zero stress state from in vivo plaque geometries. This paper presents a method to quantify human carotid artery axial and inner circumferential shrinkages by using patient-specific ex vivo and in vivo MRI images. A shrink-stretch process based on patient-specific in vivo plaque morphology and shrinkage data was introduced to shrink the in vivo geometry first to find the zero-stress state (opening angle was ignored to reduce the complexity), and then stretch and pressurize to recover the in vivo plaque geometry with computed initial stress, strain, flow pressure and velocity conditions. Effects of the shrink-stretch process on plaque stress/strain distributions were demonstrated based on patient-specific data using 3D models with fluid-structure interactions (FSI). The average artery axial and inner circumferential shrinkages were 25% and 7.9%, respectively, based on a data set obtained from 10 patients. Maximum values of maximum principal stress and strain increased 349.8% and 249% respectively with 33% axial stretch. Influence of inner circumferential shrinkage (7.9%) was not very noticeable under 33% axial stretch, but became more noticeable under smaller axial stretch. Our results indicated that accurate knowledge of artery shrinkages and the shrink-stretch process will considerably improve the accuracy of computational predictions made based on results from those in vivo MRI-based FSI models.     

Fluid-structure interaction modeling of aneurysmal arteries under steady-state and pulsatile blood flow: a stability analysis

Tortuous aneurysmal arteries are often associated with a higher risk of rupture but the mechanism remains unclear. The goal of this study was to analyze the buckling and post-buckling behaviors of aneurysmal arteries under pulsatile flow. To accomplish this goal, we analyzed the buckling behavior of model carotid and abdominal aorta with aneurysms by utilizing fluid-structure interaction (FSI) method with realistic waveforms boundary conditions. FSI simulations were done under steady-state and pulsatile flow for normal (1.5) and reduced (1.3) axial stretch ratios to investigate the influence of aneurysm, pulsatile lumen pressure and axial tension on stability. Our results indicated that aneurysmal artery buckled at the critical buckling pressure and its deflection nonlinearly increased with increasing lumen pressure. Buckling elevates the peak stress (up to 118%). The maximum aneurysm wall stress at pulsatile FSI flow was (29%) higher than under static pressure at the peak lumen pressure of 130 mmHg. Buckling results show an increase in lumen shear stress at the inner side of the maximum deflection. Vortex flow was dramatically enlarged with increasing lumen pressure and artery diameter. Aneurysmal arteries are more susceptible than normal arteries to mechanical instability which causes high stresses in the aneurysm wall that could lead to aneurysm rupture.     

Carotid geometry effects on blood flow and on risk for vascular disease

It has been widely observed that atherosclerotic diseases occur at sites with complex hemodynamics, such as artery bifurcations, junctions, and regions of high curvature. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present work, 3D pulsatile blood flow through a model of the carotid artery bifurcation was simulated using a finite volume numerical method. The goal was to quantify the risk of atherogenesis associated with different carotid artery geometries. A risk scale based on the average WSS on the sinus wall of the internal carotid artery was proposed-a scale that can be used to quantify the effect of the carotid geometry on the relative risk for developing vascular disease. It was found that the bifurcation angle and the out-of-plane angle of the internal carotid artery affect the formation of low stress regions on the carotid walls. The main conclusions are: (a) larger internal carotid artery angles (theta(IC)) generally increase the frequency and the area of blood recirculation and lower the WSS on the sinus wall, hence increasing the risk of plaque build-up; (b) off-plane angles were found to lower the WSS on the sinus for geometries with theta(IC)25 degrees . Larger off-plane angles generally increase the danger of plague build-up; (c) for theta(IC) < 25 degrees , the off-plane angle does not have an obvious effect on the hemodynamic WSS; (d) symmetric bifurcations were found to increase the WSS on the sinus wall and ease the risk of vascular disease.     

A theoretical and non-destructive experimental approach for direct inclusion of measured collagen orientation and recruitment into mechanical models of the artery wall

Gradual collagen recruitment has been hypothesized as the underlying mechanism for the mechanical stiffening with increasing stress in arteries. In this work, we investigated this hypothesis in eight rabbit carotid arteries by directly measuring the distribution of collagen recruitment stretch under increasing circumferential loading using a custom uniaxial (UA) extension device combined with a multi-photon microscope (MPM). This approach allowed simultaneous mechanical testing and imaging of collagen fibers without traditional destructive fixation methods. Fiber recruitment was quantified from 3D rendered MPM images, and fiber orientation was measured in projected stacks of images. Collagen recruitment was observed to initiate at a finite strain, corresponding to a sharp increase in the measured mechanical stiffness, confirming the previous hypothesis and motivating the development of a new constitutive model to capture this response. Previous constitutive equations for the arterial wall have modeled the collagen contribution with either abrupt recruitment at zero strain, abrupt recruitment at finite strain or as gradual recruitment beginning at infinitesimal strain. Based on our experimental data, a new combined constitutive model was presented in which fiber recruitment begins at a finite strain with activation stretch represented by a probability distribution function. By directly including this recruitment data, the collagen contribution was modeled using a simple Neo-Hookean equation. As a result, only two phenomenological material constants were required from the fit to the stress stretch data. Three other models for the arterial wall were then compared with these results. The approach taken here was successful in combining stress-strain analysis with simultaneous microstructural imaging of collagen recruitment and orientation, providing a new approach by which underlying fiber architecture may be quantified and included in constitutive equations.     

Aortic valve leaflet wall shear stress characterization revisited: impact of coronary flow

Computational characterizations of aortic valve hemodynamics have typically discarded the effects of coronary flow. The objective of this study was to complement our previous fluid–structure interaction aortic valve model with a physiologic coronary circulation model to quantify the impact of coronary flow on aortic sinus hemodynamics and leaflet wall shear stress (WSS). Coronary flow suppressed vortex development in the two coronary sinuses and altered WSS magnitude and directionality on the three leaflets, with the most substantial differences occurring in the belly and tip regions.     

Baroreceptor influence on a spinal cardiovascular reflex

A stretch of the walls of the thoracic aorta, performed in vagotomized cats without obstructing aortic flow, induces increases in heart rate, myocardial contractility, and arterial pressure. These reflex responses are still present after high spinal section. Cats under chloralose-urethane anesthesia were vagotomized and one carotid sinus was isolated and perfused with arterial blood at constant flow. The contralateral carotid sinus nerve and both aortic nerves were sectioned. A stretch of the walls of the thoracic aorta between the 7th and 10th intercostal arteries induced a reflex increase in mean arterial pressure 29 +/- 2 mmHg (mean +/- SE). Stepwise increases of carotid sinus pressure (CSP) or electrical stimulation of the carotid sinus nerve induced stepwise decreases of this reflex response. At maximal baroreceptor stimulation (CSP 212 +/- 9 mmHg) the reflex response to aortic stretch was reduced by 42%. These experiments show that this spinal cardiovascular reflex is at least partially under the inhibitory control of the baroreceptor input.     

Opposition of rapid baroreceptor resetting by prostanoids in rabbits

Arterial baroreceptors reset rapidly within minutes during acute hypertension; baroreceptor pressure threshold (Pth) is increased and the pressure-baroreceptor activity relation is shifted to the right. The purpose of the present study was to determine if prostacyclin (PGI2) or other prostanoids, released during acute hypertension modulate the magnitude of baroreceptor resetting. Baroreceptor activity was recorded from the vascularly-isolated carotid sinus during distension of the sinus with slow pressure ramp in rabbits anesthetized with chloralose. Pressure-activity curves were generated after holding carotid sinus pressure for 10-15 min from 30 to 100 mmHg. In control, the elevation of holding pressure increased Pth from 44+/- to 65+/-5 mmHg (p < 0.05, n = 12). In the presence of PGI2 (20 microM), Pth averaged 43+/-4 and 45+/-3 mmHg (n = 12) after holding pressure at 30 and 100 mmHg, respectively. In the control group before exposing the carotid sinus to indomethacin, an elevation of holding pressure increased Pth from 49+/-2 to 71+/-3 mmHg (p < 0.05, n = 12). After inhibition of the endogenous formation of prostanoids with indomethacin (20 microM), Pth increased by a significantly greater extent from 61+/-2 to 90+/-3 mmHg (p < 0.05, n = 12) with the increase in holding pressure. The slope of the pressure-activity curve (baroreceptor gain) was not influenced by the change in holding pressure. It was increased significantly by PGI2, while decreased by indomethacin. Neither the change in holding pressure nor PGI2 affected the circumferential wall strain of carotid sinus over a wide range of pressure alteration. The results suggest that PGI2 or other prostanoids released during acute hypertension sensitizes baroreceptors and provides a negative feedback mechanism that opposes and limits the magnitude of rapid baroreceptor resetting.     

The impact of imperfect frame deployment and rotational orientation on stress within the prosthetic leaflets during transcatheter aortic valve implantation

TAVI devices are manufactured with cylindrical frames. However, the frames are rarely cylindrical post-deployment since deformation due to localised under expansion can be induced by calcified material on the native valve leaflets exerting irregular forces upon the frame. Consequently, the leaflets within a deformed TAVI device may undergo elevated stress during operation, which may lead to premature device failure.Using computational analysis a complete TAVI device model was simulated undergoing deployment into an aortic root model derived from CT data for a patient with severe calcific aortic stenosis, followed by a pressure simulated cardiac cycle. The complete analysis was performed eight times, each with the device at a different rotational orientation relative to the native valve, with an increment spacing of 15°.The TAVI device frames consistently featured significant distortions associated with bulky calcified material at the base of the non-coronary sinus. It was found that the average von Mises stress in the prosthetic valves was only increased in one of the cases relative to an idealised device. However, the maximum von Mises stress in the prosthetic valves was elevated in the majority of the cases.Furthermore, it was found that there were preferable orientations to deploy the prosthetic device, in this case, when the prosthetic leaflets were aligned with the native leaflets. As device orientation deviated from this orientation, the stresses in the valve increased because the distance between the prosthetic commissures decreased. This potentially could represent a sufficient increase in stress to induce variation in device lifespan.     

Laser Doppler anemometer measurements of pulsatile flow in a model carotid bifurcation

Hemodynamics at the human carotid bifurcation is important to the understanding of atherosclerotic plaque initiation and progression as well as to the diagnosis of clinically important disease. Laser Doppler anemometry was performed in a large scale model of an average human carotid. Pulsatile waveforms and physiologic flow divisions were incorporated. Disturbance levels and shear stresses were computed from ensemble averages of the velocity waveform measurements. Flow in the common carotid was laminar and symmetric. Flow patterns in the sinus, however, were complex and varied considerably during the cycle. Strong helical patterns and outer wall flow separation waxed and waned during each systole. The changing flow patterns resulted in an oscillatory shear stress at the outer wall ranging from -13 to 9 dyn cm-2 during systole with a time-averaged mean of only -0.5 dyn cm-2. This contrasts markedly with an inner wall shear stress range of 17-50, (mean 26) dyn cm-2. The region of transient separation was confined to the carotid sinus outer wall with no reverse velocities detected in the distal internal carotid. Notable disturbance velocities were also time-dependent, occurring only during the deceleration phase of systole and the beginning of diastole. The present pulsatile flow studies have aided in identifying hemodynamic conditions which correlate with early intimal thickening and predict the physiologic level of flow disturbances in the bulb of undiseased internal carotid arteries.     

Surrounding tissues affect the passive mechanics of the vessel wall: theory and experiment

The stress and strain in the vessel wall are important determinants of vascular physiology and pathophysiology. Vessels are constrained radially by the surrounding tissue. The hypothesis in this work is that the surrounding tissue takes up a considerable portion of the intravascular pressure and significantly reduces the wall strain and stress. Ten swine of either sex were used to test this hypothesis. An impedance catheter was inserted into the carotid or femoral artery, and after mechanical preconditioning pressure-cross-sectional area relations were obtained with the surrounding tissue intact and dissected away (untethered), respectively. The radial constraint of the surrounding tissue was quantified as an effective perivascular pressure on the outer surface of the vessel, which was estimated as 50% or more of the intravascular pressure. For carotid arteries at pressure of 100 mmHg, the circumferential wall stretch ratio in the intact state was approximately 20% lower than in the untethered state and the average circumferential stress was reduced by approximately 70%. For femoral arteries, the reductions were approximately 15% and 70%, respectively. These experimental data support the proposed hypothesis and suggest that in vitro and in vivo measurements of the mechanical properties of vessels must be interpreted with consideration of the constraint of the surrounding tissue.     

Multi-patient study for coronary vulnerable plaque model comparisons: 2D/3D and fluid–structure interaction simulations

Several image-based computational models have been used to perform mechanical analysis for atherosclerotic plaque progression and vulnerability investigations. However, differences of computational predictions from those models have not been quantified at multi-patient level. In vivo intravascular ultrasound (IVUS) coronary plaque data were acquired from seven patients. Seven 2D/3D models with/without circumferential shrink, cyclic bending and fluid–structure interactions (FSI) were constructed for the seven patients to perform model comparisons and quantify impact of 2D simplification, circumferential shrink, FSI and cyclic bending plaque wall stress/strain (PWS/PWSn) and flow shear stress (FSS) calculations. PWS/PWSn and FSS averages from seven patients (388 slices for 2D and 3D thin-layer models) were used for comparison. Compared to 2D models with shrink process, 2D models without shrink process overestimated PWS by 17.26%. PWS change at location with greatest curvature change from 3D FSI models with/without cyclic bending varied from 15.07% to 49.52% for the seven patients (average = 30.13%). Mean Max-FSS, Min-FSS and Ave-FSS from the flow-only models under maximum pressure condition were 4.02%, 11.29% and 5.45% higher than those from full FSI models with cycle bending, respectively. Mean PWS and PWSn differences between FSI and structure-only models were only 4.38% and 1.78%. Model differences had noticeable patient variations. FSI and flow-only model differences were greater for minimum FSS predictions, notable since low FSS is known to be related to plaque progression. Structure-only models could provide PWS/PWSn calculations as good approximations to FSI models for simplicity and time savings in calculation.

     

Unloading arterial baroreceptors causes neurogenic hypertension

We developed a new model to examine the role of arterial baroreceptors in the long-term control of mean arterial pressure (MAP) in dogs. Baroreceptors in the aortic arch and one carotid sinus were denervated, and catheters were implanted in the descending aorta and common carotid arteries. MAP and carotid sinus pressure (CSP) averaged 104 +/- 2 and 102 +/- 2 mmHg (means +/- 1 SE), respectively, during a 5-day control period. Baroreceptor unloading was induced by ligation of the common carotid artery proximal to the innervated sinus (n = 6 dogs). MAP and CSP averaged 127 +/- 7 and 100 +/- 3 mmHg, respectively, during the 7-day period of baroreceptor unloading. MAP was significantly elevated (P < 0.01) compared to control, but CSP was unchanged. Heart rate and plasma renin activity increased significantly in response to baroreceptor unloading. Removal of the ligature to restore normal flow through the carotid resulted in normalization of all variables. Ligation of the carotid below a denervated sinus (n = 4) caused a significant decrease in CSP but no systemic hypertension. These results indicate that chronic unloading of carotid baroreceptors can produce neurogenic hypertension and provide strong evidence that arterial baroreceptors are involved in the long-term control of blood pressure.     

Fluid–structure interaction analysis of turbulent pulsatile flow within a layered aortic wall as related to aortic dissection

Turbulent pulsatile flow and wall mechanics were studied numerically in an axisymmetric three-layered wall model of a descending aorta. The transport equations were solved using the finite element formulation based on the Galerkin method of weighted residuals. A fully-coupled fluid–structure interaction (FSI) analysis was utilized in this investigation. We calculated Von Mises wall stress, streamlines and fluid pressure contours. The findings of this study show that peak wall stress and maximum shear stress are highest in the media layer. The difference in the elastic properties of contiguous layers of the wall of the aorta probably determines the occurrence of dissection in the media layer. Moreover, the presence of aortic intramural hematoma is found to have a significant effect on the peak wall stress acting on the inner layer.     

Anomalous hemodynamic effects of a self-expanding intracranial stent: Comparing in-vitro and ex-vivo models using ultra-high resolution MicroCT based CFD

Previous research on the effects of intracranial stents on arterial hemodynamics has involved computational hemodynamics (CHD) simulations applied to artificially generated stent models. In this study, accurate geometric reconstructions of in-vitro (PTFE tube) and ex-vivo (canine artery) deployed stents based on ultra-high resolution MicroCT imaging were used. The primary goal was to compare the hemodynamic effects of deployment in these two different models and to identify flow perturbations due to deployment anomalies such as stent malapposition and strut prolapse, important adverse mechanics occurring in clinical practice, but not considered in studies using idealized stent models.Ultra-high resolution MicroCT data provided detailed visualization of deployment characteristics allowing for accurate in-stent flow simulation. For stent cells that are regularly and symmetrically deployed, the near wall flow velocities and wall shear stresses were similar to previously published results derived from idealized models. In-stent hemodynamics were significantly altered by misaligned or malapposed stent cells, important effects not realistically captured in previous models. This research shows the feasibility and value of an ex-vivo stent model for MicroCT based CHD studies. It validates previous in-vitro studies and further contributes to the understanding of in-stent hemodynamics associated with adverse mechanics of self-expanding intracranial stents.     

Reflex mechanisms for changes in renal nerve activity during positive end-expiratory pressure

This study was designed to investigate the interaction between carotid sinus baroreceptors and cardiopulmonary receptors in the reflex control of renal nerve activity (RNA) during positive end-expiratory pressure (PEEP) in anesthetized dogs. PEEP at two different levels (10 and 20 cmH2O) was applied to the following groups: animals with neuraxis intact (I group, n = 12); vagal and aortic nerve denervated animals with carotid sinus nerves intact (V group, n = 6); carotid sinus denervated animals with vagal and aortic nerves intact (SD group, n = 6); and carotid sinus denervated animals also having severed vagal and aortic nerves (SAV group, n = 12). Mean blood pressure (MBP), central venous pressure, and mean airway pressure were also simultaneously measured. In the I group, no significant alterations in RNA occurred during PEEP at both levels, even when MBP fell significantly. Although the drop in MBP in the SD group was similar to that in the I group, RNA decreased significantly 10 s after intervention at both PEEP levels, followed by a recovery of RNA toward the control level. In contrast, a significant increase in RNA, which continued until the end of PEEP, appeared in the V group immediately after each intervention. In the SAV group, RNA responses to PEEP, which were observed in the other groups, were abolished. These results provide evidence that during PEEP, renal nerve activity is modified by an interaction between carotid sinus baroreceptors and cardiopulmonary receptors; excitatory effects occur via carotid sinus nerves and inhibitory effects occur via vagal afferents.     

Impedance-based outflow boundary conditions for human carotid haemodynamics

In this study, we develop structured tree outflow boundary conditions for modelling the human carotid haemodynamics. The model geometry was reconstructed through computerised tomography scan. Unsteady-state computational fluid dynamic analyses were performed under different conditions using a commercial software package ADINA R&D, Inc., (Watertown, MA, USA) in order to assess the impact of the boundary conditions on the flow variables. In particular, the results showed that the peripheral vessels massively impact the pressure while the flow is relatively unaffected. As an example of application of these outflow conditions, an unsteady fluid-structure interaction (FSI) simulation was carried out and the dependence of the wall shear stress (WSS) on the arterial wall compliance in the carotid bifurcation was studied. In particular, a comparison between FSI and rigid-wall models was conducted. Results showed that the WSS distributions were substantially affected by the diameter variation of the arterial wall. In particular, even similar WSS distributions were found for both cases, and differences in the computed WSS values were also found.     

Cyclic Bending Contributes to High Stress in a Human Coronary Atherosclerotic Plaque and Rupture Risk: In Vitro Experimental Modeling and Ex Vivo MRI-Based Computational Modeling Approach

Many acute cardiovascular syndromes such as heart attack and stroke are caused by atherosclerotic plaque ruptures which often happen without warning. MRI-based models with fluid-structure interactions (FSI) have been introduced to perform flow and stress/strain analysis for atherosclerotic plaques and identify possible mechanical and morphological indices for accurate plaque vulnerability assessment. In this paper, cyclic bending was added to 3D FSI coronary plaque models for more accurate mechanical predictions. Curvature variation was prescribed using the data of a human left anterior descending (LAD) coronary artery. Five computational models were constructed based on ex vivo MRI human coronary plaque data to assess the effects of cyclic bending, pulsating pressure, plaque structure, and axial stretch on plaque stress/strain distributions. In vitro experiments using a hydrogel stenosis model with cyclical bending were performed to observe effect of cyclical bending on flow conditions. Our results indicate that cyclical bending may cause more than 100% or even up to more than 1000% increase in maximum principal stress values at locations where the plaque is bent most. Stress increase is higher when bending is coupled with axial stretch, non-smooth plaque structure, or resonant pressure conditions (zero phase angle shift). Effects of cyclic bending on flow behaviors are more modest (21.6% decrease in maximum velocity, 10.8% decrease in flow rate, maximum flow shear stress changes were < 5%). Computational FSI models including cyclic bending, plaque components and structure, axial stretch, accurate in vivo measurements of pressure, curvature, and material properties should lead to significant improvement on stress-based plaque mechanical analysis and more accurate coronary plaque vulnerability assessment.     

Diameter-dependent axial prestretch of porcine coronary arteries and veins

The pressure-diameter relation (PDR) and the wall strain of coronary blood vessels have important implications for coronary blood flow and arthrosclerosis, respectively. Previous studies have shown that these mechanical quantities are significantly affected by the axial stretch of the vessels. The objective of this study was to measure the physiological axial stretch in the coronary vasculature; i.e., from left anterior descending (LAD) artery tree to coronary sinus vein and to determine its effect on the PDR and hence wall stiffness. Silicone elastomer was perfused through the LAD artery and coronary sinus trees to cast the vessels at the physiologic pressure. The results show that the physiological axial stretch exists for orders 4 to 11 (> 24 μm in diameter) arteries and orders -4 to -12 (>38 μm in diameter) veins but vanishes for the smaller vessels. Statistically, the axial stretch is higher for larger vessels and is higher for arteries than veins. The axial stretch λ(z) shows a linear variation with the order number (n) as: λ(z) = 0.062n + 0.75 (R(2) = 0.99) for artery and λ(z) = -0.029n + 0.89 (R(2) = 0.99) for vein. The mechanical analysis shows that the axial stretch significantly affects the PDR of the larger vessels. The circumferential stretch/strain was found to be significantly higher for the epicardial arteries (orders 9-11), which are free of myocardium constraint, than the intramyocardial arteries (orders 4-8). These findings have fundamental implications for coronary blood vessel mechanics.     

Wall stress of the cervical carotid artery in patients with carotid dissection: a case-control study

Spontaneous internal carotid artery (ICA) dissection (sICAD) results from an intimal tear located around the distal carotid sinus. The mechanisms causing the tear are unknown. This case-control study tested the hypotheses that head movements increase the wall stress in the cervical ICA and that the stress increase is greater in patients with sICAD than in controls. Five patients with unilateral, recanalized, left sICAD and five matched controls were investigated before and after maximal head rotation to the left and neck hyperextension after 45° head rotation to the left. The anatomy of the extracranial carotid arteries was assessed by magnetic resonance imaging and used to create finite element models of the right ICA. Wall stress increased after head movements. Increases above the 80th and 90th percentile were located at the intimal side of the artery wall from 7.4 mm below to 10 mm above the cranial edge of the carotid sinus, i.e., at the same location as histologically confirmed tears in patients with sICAD. Wall stress increase did not differ between patients and controls. The present findings suggest that wall stress increases at the intimal side of the artery wall surrounding the distal edge of the carotid bulb after head movements may be important for the development of carotid dissection. The lack of wall stress difference between the two groups indicates that the carotid arteries of patients with carotid dissection have either distinct functional or anatomical properties or endured unusually heavy wall stresses to initiate dissection.     

The ratio of systemic and coronary fractions of left ventricular output in reflex pressor reactions

In acute experiments the participation of coronary and systemic fractions was studied during suppression of the carotid sinus baroreceptors by the occlusion of the carotid arteries and stimulation of the tibial nerve afferent fibers. In most tests systemic fraction was reduced in carotid arteries occlusion and increased in tibial nerve stimulation. The coronary fraction was always increased. The cardiac output (the sum of systemic and coronary fractions) was steady with blood pressure increase by 20 to 70%. The role of coronary fraction in the mechanism of homeometric regulation of the heart is discussed.